the near - infrared galaxy luminosity function is an important characteristic of the local galaxy population .
it is a much better tracer of evolved stars , and hence of the total stellar content of galaxies , than optical luminosity functions which can be dominated by young stellar populations and are also strongly affected by dust extinction .
hence , infrared luminosities can be much more directly related to the underlying stellar mass of galaxies and so knowledge of the present form and evolution of the infrared galaxy luminosity function places strong constraints on the history of star formation in the universe and on galaxy formation models ( cole and references therein ) .
the local k - band luminosity function has been estimated from optically selected samples by mobasher , sharples and ellis ( ) , szokoly ( ) and loveday ( ) and from k - band surveys by glazebrook ( ) , and gardner ( ) .
the existing k - band surveys are small .
the largest , by gardner et al . , covers only 4  deg@xmath12 and contains only 510 galaxies .
the recent survey of loveday covers a much larger solid angle . in this survey
the redshifts were known in advance of measuring the k - band magnitudes and this was exploited by targetting bright and faint galaxies resulting in an effective sample size much larger than the 345 galaxies actually measured . however , like all optically selected samples , it suffers from the potential problem that galaxies with extremely red infrared to optical colours could be missed . in this paper we combine the 2-micron all sky survey ( ) with the 2df galaxy redshift survey ( ) to create an infrared selected redshift survey subtending  deg@xmath12 .
currently the sky coverage of both surveys is incomplete , but already the overlap has an effective area of  deg@xmath12 . within this area
the redshift survey is complete to the magnitude limit of the catalogue and so constitutes a complete survey which is 50  times larger than the previous largest published infrared selected redshift survey .
a new catalogue of a similarly large area , also based on , has very recently been analysed by kochanek ( ) .
they adopt isophotal rather than total magnitudes and concentrate on the dependence of the luminosity function on galaxy morphology .    in section  [ sec :
selection ] we briefly describe the relevant properties of the  and  catalogues .
section  [ sec : matching ] is a detailed examination of the degree to which the matched
galaxies are a complete and representative subset of the  catalogue .
section  [ sec : mags ] examines the calibration of the  total magnitudes and section  [ sec : compl ] demonstrates that the  catalogue and the inferred luminosity functions are not affected by surface brightness selection effects . in section  [ sec : pop ] we present the method by which we compute k - corrections and evolutionary corrections and relate the observed luminosities to the underlying stellar mass .
the estimation methods and normalization of the luminosity functions are described briefly in section  [ sec : lf_est ] .
our main results are presented and discussed in section  [ sec : results ] .
these include estimates of the j and ( k - short ) luminosity functions , the @xmath9  and j@xmath9  colour distributions as a function of absolute magnitude and the distribution of spectral type .
we also estimate the stellar mass function of galaxies , which can be integrated to infer the fraction of baryons in the universe which are in the form of stars .
we conclude in section  [ sec : conc ] .
the data that we analyze are the extended source catalogue from the second incremental release of the 2-micron all sky survey ( http://pegasus.phast.umass.edu ) and the galaxy catalogue of the 2df galaxy redshift survey (  http://www.mso.anu.edu.au/2dfgrs ) . here
, we present the relevant properties of these two catalogues and investigate their selection characteristics and level of completeness .
the  is a ground - based , all - sky imaging survey in the j , h and  bands .
details of how extended sources are identified and their photometric properties measured are given by jarrett ( ) . the detection sensitivity ( 10@xmath13 ) for extended sources is quoted as 14.7 , 13.9 and 13.1 magnitudes in j , h and  respectively .
the complete survey is expected to contain 1  million galaxies of which approximately 580,000 are contained in the second incremental data release made public in march 2000 .
the  is selected in the photographic  band from the apm galaxy survey ( @xcite , , ) and subsequent extensions to it , that include a region in the northern galactic cap ( maddox in preparation ) .
the survey covers approximately  deg@xmath12 consisting of two broad declination strips .
the larger is centred on the sgp and approximately covers @xmath14@xmath15@xmath16@xmath15@xmath17 , @xmath18@xmath19@xmath20@xmath19@xmath21 ; the smaller strip is in the northern galactic cap and covers @xmath22@xmath15@xmath16@xmath15@xmath23 , @xmath24@xmath19@xmath20@xmath19@xmath25 .
in addition , there are a number of randomly located circular 2-degree fields scattered across the full extent of the low extinction regions of the southern apm galaxy survey .
there are some gaps in the  sky coverage within these boundaries due to small regions that have been excluded around bright stars and satellite trails .
the  aims to measure the redshifts of all the galaxies within these boundaries with extinction - corrected  magnitudes brighter than 19.45 .
when complete , at the end of 2001 , 250,000 galaxy redshifts will have been measured . in this paper
we use the  redshifts obtained prior to .
the overlap of the two surveys is very good .
there are some gaps in the sky coverage due to strips of the sky that were not included in the  second incremental release , but overall a substantial fraction of the  deg@xmath12 of the  is covered by .
the homogeneity and extensive sky coverage of the combined dataset make it ideal for studies of the statistical properties of the galaxy population .
= 8 truecm      here we consider whether all the  galaxies within the survey region have  counterparts and assess the extent to which the fraction of galaxies with measured redshifts represents an unbiased sub - sample .
the astrometry in both  and  is , in general , very good and it is an easy matter to match objects in the two catalogues . we choose to find the closest pairs within a search radius equal to three quarters of the semi - major axis of the j - band image ( denoted j_r_e in the  database ) . scaling the search radius in this way helps with the matching of large extended objects .
this procedure results in the identification of  counterparts for  of the  objects , when at random one would only expect to find a handful of such close pairs .
moreover , the distribution of separations shown in fig .
[ fig : seps ] peaks at 0.5  arcsec , with only 3% having separations greater than 3  arcsec .
a significant part of this tail comes from the most extended objects as is evident from the dotted histogram in fig .
[ fig : seps ] which shows objects with semi - major axes larger than 12  arcsec .
thus , we can be very confident in these identifications .
the  objects for which we have found secure  counterparts amount to % of the  extended sources that fall within the boundary of the . as discussed below , a more restrictive criterion that includes only sources fainter than j@xmath26@xmath27 that are confidently classified as galaxies by , increases the fraction with  matches to % .
the remaining % are missed for well understood reasons ( star - galaxy classification : % ; merged or close images : % ; miscellaneous : % ) , none of which ought to introduce a bias .
this is confirmed explicitly , in the middle row of fig .  [
fig : hists ] , by the close correspondence between the photometric properties of the missed % and those of the larger matched sample .
hence , in estimating luminosity functions no significant bias will be introduced by assuming the matched sample to be representative of the full population .
furthermore , the distribution shown in the bottom row of fig .
[ fig : hists ] shows that the subset of galaxies for which we have measured redshifts is a random sample of the full matched   catalogue .
this summary is the result of a thorough investigation , which we describe in the remainder of this section , into the reasons why % of the  sources are missed and whether their omission introduces a bias in the properties of the matched sample .
= 15.0 truecm    we first consider objects in the  catalogue which based on their images and colours are not confidently classified as galaxies . in the  database a high e_score or g_score indicates a high probability that the object is either not an extended source or not a galaxy .
a cc_flag@xmath280 indicates an artifact or contaminated and/or confused source .
for detailed definitions of these parameters we refer the reader to jarrett ( ) .
rejecting all objects which have either e_score@xmath151.4 , g_score@xmath151.4 or cc_flag@xmath280 removes just % of the total .
however , removing these reduces significantly the fraction of the sample that does not match with the  catalogue , from % to % .
thus , it is likely that about 30% of the  objects which have e_score@xmath151.4 , g_score@xmath151.4 or cc_flag@xmath280 are not galaxies .
the  may contain a tail of very red objects that are too faint in the -band to be included in the @xmath1919.45 sample .
[ fig : colours ] shows the distribution of -j colours for the matched objects with j@xmath1914.7 .
( here , the j - band magnitude we are using is the default magnitude denoted j_m in the database . in section  [ sec : mags ]
we will consider the issue of what magnitude definition is most appropriate for estimating the luminosity function . )
the vertical dashed line indicates the colour at which this sample starts to become incomplete due to the @xmath1919.45 magnitude limit of .
the colour distribution cuts off sharply well before this limit , suggesting that any tail of missed very red objects is extremely small . in other words
the  is sufficiently deep that even the reddest objects detected at the faintest limits of  ought to be detected in .    in the top row of fig .
[ fig : hists ] we compare the distributions of magnitude , colour and surface brightness for the matched and missed  objects . in general , the properties of the missed subset overlap well with those of the much larger matched subset .
however , we do see that the distributions for missed objects contain tails of bright and blue objects .
it is quite likely that this is due to the  extended source catalogue being contaminated by a small population of saturated or multiple stars .
the dotted histograms in the top row of fig .
[ fig : hists ] show the distributions of magnitude , colour and surface brightness for the bright subset of the missed objects with j@xmath1912 . here
we clearly see bimodal colour and surface brightness distributions .
the blue peak of the colour distribution is consistent with that expected for stars ( see jarrett ) .
excluding these bright , j@xmath1912 , objects which are clearly contaminated by stars reduces the fraction of missed  objects from % to % .
the magnitude , colour and surface brightness distributions for this remaining % are shown in the middle panel of fig .
[ fig : hists ] .
we see that the missed objects are slightly under - represented at the faintest magnitudes and also slightly bluer on average than the matched sample , while the distribution of surface brightness is almost indistinguishable for the two sets of objects .
these differences are small and so will introduce no significant bias in our luminosity function estimates . to elucidate the reasons for the remaining missed % of  objects we downloaded 100 1@xmath291  arcmin images from the stsci digitized sky survey ( dss ) centred on the positions of a random sample of the missed  objects . in each image we plotted a symbol to indicate the position of any  galaxies within the 1@xmath291  arcmin field .
we also plotted symbols to indicate the positions and classifications of all images identified in the apm scans from which the  catalogue was drawn , down to a magnitude limit of @xmath30 .
these images are classified as galaxies , stars , merged images ( galaxy+galaxy , galaxy+star or star+star ) or noise .
this set of plots allows us to perform a census of the reasons why some  objects are not present in the  survey .
= 8.5 truecm    the main cause for the absence of objects in the  is that the apm has classified these objects as stars .
these amount to % of the missed sample ( % of the full  sample ) . in some cases ,
the dss image shows clearly that these are stars and in others that they are galaxies .
however , the majority of these objects can not easily be classified from the dss images .
thus , they could be galaxies that the apm has falsely classified as stars or stars that  has falsely classified as galaxies .
the first possibility is not unexpected since the parameters used in the apm star - galaxy separation algorithm were chosen as a compromise between high completeness and low contamination such that the expected completeness is around 95% with 5% stellar contamination ( @xcite ) .
it is hard to rule out the possibility that this class of object does not include a substantial fraction of stars , but if so , their presence appears not to distort the distribution of colours shown in fig  [ fig : colours ] .
another % of the random sample ( % of the full  sample ) are classified by the apm as mergers or else consist of two close images in the dss but are classified by the apm as a single galaxy offset from the  position .
the remaining % of the random sample ( % of the full  sample ) are missed for a variety of reasons including proximity to the diffraction spikes of very bright stars and poor astrometry caused by the presence of a neighbouring unclassified image .
= 9.0 truecm    = 18.0 truecm    the  extended source database provides a large selection of different magnitude measurements . in the previous section we used the default magnitudes ( denoted j_m and k_m in the database )
these are magnitudes defined within the same circular aperture in each waveband . for galaxies brighter than @xmath26@xmath31 ,
the aperture is the circular -band isophote of 20  mag arcsec@xmath32 and for galaxies fainter than @xmath2614 it is the circular j - band isophote of 21  mag arcsec@xmath32 .
these are not the most useful definitions of magnitude for determining the galaxy luminosity function .
since we are interested in measuring the total luminosity and ultimately the total stellar mass of each galaxy , we require a magnitude definition that better represents the total flux emitted by each galaxy .
we consider kron magnitudes ( kron ) and extrapolated magnitudes .
kron magnitudes ( denoted j_m_e and k_m_e in the  database ) are measured within an aperture , the kron radius , defined as 2.5 times the intensity - weighted radius of the image .
the extrapolated magnitudes ( denoted j_m_ext and k_m_ext in the  database ) are defined by first fitting a modified exponential profile , @xmath33 $ ] , to the image from 10  arcsec to the 20  mag / arcsec@xmath12 isophotal radius , and extrapolating this from the kron radius to 4 times this radius or 80  arcsec if this is smaller ( jarrett private communication ) .
note that improvements are being made to the extended source photometry algorithms developed and employed  team and so in the final  data release the definitions of the kron and extrapolated magnitudes may be slightly modified ( jarrett private communication ) .
[ fig : mags ] compares the default , kron and extrapolated magnitudes in the j and  bands for the matched   catalogue .
the upper panel shows that while the median offset between the j - band isophotal default magnitudes and the pseudo - total kron magnitudes is small there is a large spread with some galaxies having default magnitudes more than 0.5  magnitudes fainter than the kron magnitude .
the kron magnitudes are systematically fainter than the extrapolated magnitudes by between approximately 0.1 and 0.3  magnitudes .
this offset is rather larger than expected : if the kron radius is computed using a faint isophote to define the extent of the image from which the intensity weighted radius is measured , then the kron magnitudes should be very close to total .
for an exponential light profile ( @xmath34 ) , the kron radius should capture 96% of the flux , while for an @xmath35 law ( @xmath36 ) , 90% of the flux should be enclosed . in other words ,
the kron magnitude should differ from the total magnitude by only 0.044 and 0.11  magnitudes in these two cases .
however , the choice of isophote is a compromise between depth and statistical robustness . in the case of the  second incremental release , an isophote of 21.7(20.0 )  mag arcsec@xmath32 in j ( )
was adopted ( jarrett private communication ) .
these relatively bright isophotes , particularly the -band isophote , could lead to underestimates of the kron radii and fluxes for lower surface brightness objects and plausibly accounts for much of the median offset of 0.3  magnitudes seen in fig .
[ fig : mags ] between the -band kron and extrapolated magnitudes .
this line of reasoning favours adopting the extrapolated magnitudes as the best estimate of the total magnitudes , but , on the other hand , the extrapolated magnitudes are model - dependent and have larger measurement errors .    to understand better the offset and scatter in the magnitudes
we have compared a subset of the  data with the independent k - band photometry of loveday ( 2000 ) .
the pointed observations of loveday have better resolution than the 2mass images and good signal - to - noise to a much deeper isophote .
this enables accurate , unbiased kron magnitudes to be measured .
note that the offset between the  -band and the standard k - band used by loveday is expected to be almost completely negligible ( see carpenter ) .
the left hand panels of fig .
[ fig : loveday ] compare these measurements with the corresponding  kron and extrapolated magnitudes .
the right hand panels show -band kron and extrapolated magnitudes computed by taking the  j - band kron and extrapolated magnitudes and subtracting the j@xmath9  colour measured within the default aperture .
these indirect estimates are interesting to consider as they combine the profile information from the deeper j - band image with the j@xmath9colour measured within the largest aperture in which there is good signal - to - noise .
the straight lines plotted in fig .  [ fig : loveday ] show simple least squares fits and the slope and zero - point offset of these fits are indicated on each panel along with bootstrap error estimates . also shown in the inset panels
is the distribution of residual magnitude differences about each of the fits and a gaussian fit to this distribution .
the rms of these residuals and a bootstrap error estimate is also given in each panel .
from these comparisons we first see that all the fits have slopes entirely consistent with unity , but that their zeropoints and scatters vary .
the zero - point offsets , @xmath37 , between both the  kron magnitude measurements and those of loveday confirm that the  kron magnitudes systematically underestimate the galaxy luminosities . in the case of the direct -band
magnitudes the offset is @xmath38  magnitudes . in the case of the kron magnitudes inferred from the deeper j - band image profiles , the offset is reduced to @xmath39  magnitudes . conversely the  extrapolated magnitudes are systematically brighter than the loveday kron magnitudes by @xmath40 and  @xmath41  magnitudes , where one would expect an offset of only @xmath42 to @xmath43 due to the difference in definition between ideal kron and true total magnitudes . for both estimates of the extrapolated magnitude and for the directly estimated kron magnitude
the scatter about the correlation is approximately 0.14  magnitudes and we note a slight tendency for the scatter to increase at faint magnitudes . the magnitude estimate that best correlates with the loveday measurements is the kron magnitude estimated from the  j - band kron magnitude and the default aperture j@xmath9  colour . here
the distribution of residuals has a much reduced scatter of only 0.1  magnitudes and has very few outliers .
our conclusions from the comparison of kron magnitudes is that it is preferable to adopt the -band magnitude inferred from the j - band kron or extrapolated magnitude by converting to the -band using default aperture colour , rather than to use the noisier and more biased direct -band estimates . with this definition , we find that the  kron magnitudes slightly underestimate the galaxy luminosities while the extrapolated magnitudes slightly overestimate the luminosities , particularly at faint fluxes
. we will present results for both magnitude definitions , but we note that to convert to total magnitudes we estimate that the  kron magnitudes should be brightened by @xmath44@xmath45 magnitudes and the extrapolated magnitudes dimmed by @xmath46@xmath43 magnitudes .      here
we define the magnitude limited samples which we will analyze in section  [ sec : lf_est ] and test them for possible incompleteness in both magnitude and surface brightness . for the kron and extrapolated magnitudes ,
the  catalogue has high completeness to the nominal limits of j@xmath1914.7 and @xmath1913.9 .
however , to ensure very high completeness and avoid any bias in our luminosity function estimates , we made the following more conservative cuts .
for the kron magnitudes , we limited our sample to either j@xmath1914.45 or @xmath1913.2 , and for the extrapolated magnitudes to either j@xmath1914.15 or @xmath1912.9 .
these choices are motivated by plots such as the top panel of fig .
[ fig : mags ] . here
the isophotal default magnitude limit of j@xmath19@xmath47 is responsible for the right hand edge to the distribution of data points .
one sees that this limit begins to remove objects from the distribution of kron magnitudes for j@xmath48@xmath49 .
an indication that the survey is complete to our adopted limits is given by the number counts shown in fig  [ fig : counts ] , which only begin to roll over at fainter magnitudes .
= 8.5 truecm    = 8.0 truecm    more rigorously , we have verified that the samples are complete to these limits by examining their @xmath50 distributions . here
, @xmath51 is the redshift of a galaxy in the sample , @xmath52 is the maximum redshift at which this galaxy would satisfy the sample selection criteria , and @xmath53 is the survey volume that lies at redshift less than @xmath54 . if the sample is complete and of uniform density , @xmath55 is uniformly distributed within the interval 0 to  1 . to evaluate @xmath56 we made use of the default k+e corrections described in the following section , but the results are not sensitive to reasonable variations in the assumed corrections or in the cosmology .
the solid histograms in fig .
[ fig : v_vmax ] show these distributions for each of our four magnitude limited samples .
note that the -band magnitudes are those inferred from the j - band magnitudes and aperture colours .
the dashed histograms in the lower panels show the corresponding distributions for the directly measured -band magnitudes . in all these cases
we have computed @xmath57 simply from the imposed apparent magnitude limits and have ignored any possible dependence of the catalogue completeness on surface brightness .
if the samples were incomplete the symptom one would expect to see is a deficit in the @xmath58 distributions at large @xmath58 and hence a mean @xmath59@xmath19@xmath60 .
there is no evidence for such a deficit in these distributions .
in fact each has a mean @xmath59 slightly greater than @xmath60 .
the slight gradient in the @xmath58 distribution is directly related to the galaxy number counts shown in fig  [ fig : counts ] , which are slightly steeper than expected for a homogeneous , non - evolving galaxy distribution .
a similar result has been found in the bright -band counts ( @xcite ) .
the -band result has variously been interpreted as evidence for rapid evolution , systematic errors in the magnitude calibration , or a local hole or underdensity in the galaxy distribution ( @xcite ) . here
we note that the gradient in the @xmath58 distributions ( and also in the galaxy counts ) becomes steeper both as one switches from kron to the less reliable extrapolated magnitudes and as one switches from the j - band data to the lower signal - to - noise -band data .
this gives strong support to our decision to adopt the -band magnitudes derived from the j - band kron and extrapolated magnitudes and aperture j@xmath9  colours .
it also cautions that the mean @xmath59@xmath15@xmath60 can not necessarily be taken as a sign of evolution or a local underdensity , but may instead be related to the accuracy of the magnitude measurements .
the comparison to the observations of loveday ( ) shows no evidence for systematic errors in the magnitudes , but does not constrain the possibility that the distribution of magnitude measurement errors may become broader or skewed at fainter magnitudes .
such variations would affect the @xmath58 distributions and could produce the observed behaviour .
we conclude by noting that while the shift in the mean @xmath59 is statistically significant , it is nevertheless quite small for the samples we analyze and has little effect on the resulting luminosity function estimates .
= 8.0 truecm    we now investigate explicitly the degree to which the completeness of the  catalogue depends on surface brightness by estimating @xmath57 as a function of both absolute magnitude and surface brightness .
this is an important issue : if the catalogue is missing low - surface brightness galaxies our estimates of the luminosity function will be biased .
the approach we have taken follows that developed in cross ( ) for the .
we estimate an effective central surface brightness , @xmath61 , for each observed galaxy assuming an exponential light distribution , that the kron magnitudes are total and that the kron radii are exactly five exponential scale - lengths .
this is then corrected to redshift @xmath62 using @xmath63 to account for redshift dimming and k+e corrections ( c.f section  [ sec : pop ] ) .
note that in the  catalogue , galaxies with estimated kron radii less than @xmath64  arcsec , have their kron radii set to 7 arcsec .
this will lead us to underestimate the central surface brightnesses of these galaxies , but this will only affect high surface brightness objects and will not affect whether a galaxy can or can not be seen .
the distribution in the @xmath65@xmath66 plane of our kron j - band selected sample is shown by the points in fig .
[ fig : sb ] .    cross ( ) use two different methods to estimate the value of @xmath57 associated with each position in this plane .
the first method uses the visibility theory of phillipps , davies & disney ( ) .
we model the selection characteristics of the  extended source catalogue by a set of thresholds .
the values appropriate in the j - band are a minimum isophotal diameter of @xmath67  arcsec at an isophote of @xmath68  mag arcsec@xmath32 , and an isophotal magnitude limit of j@xmath19@xmath47 at an isophote of @xmath69  mag arcsec@xmath32 ( jarrett et al .
in addition , we impose the limits in the kron magnitude of @xmath70@xmath19j@xmath19@xmath71 that define the sample we analyze .
we then calculate for each point on the @xmath65@xmath66 plane the redshift at which a such a galaxy will drop below one or other of these selection thresholds and hence compute a value of @xmath57 .
the results of this procedure are shown by the contours of constant @xmath57 plotted in fig .
[ fig : sb ] .
note that these estimates of @xmath57 are only approximate since we have made the crude assumption that all the galaxies are circular exponential disks .
in addition , the diameter and isophotal limits are only approximate and vary with observing conditions .    the second method developed by cross ( )
consists of making an empirical estimate of @xmath57 in bins in the @xmath65@xmath66 plane .
they look at the distribution of observed redshifts in a given bin and adopt the 90@xmath72 percentile of this distribution to define @xmath56 and hence @xmath57 .
it is more robust to use the 90@xmath72 percentile rather than the 100@xmath72 percentile and the effect of this choice can easily be compensated for when estimating the luminosity function ( @xcite ) .
note that in our application to the  data we do not apply corrections for incompleteness or the effects of clustering .
the result of this procedure is to confirm that for the populated bins , the @xmath57 values given by the visibility theory are a good description of the data .    = 8.5 truecm    in fig .
[ fig : sb ] we see that the distribution of galaxies in the @xmath65@xmath66 plane is well separated from the low surface brightness limit of approximately @xmath68  mag arcsec@xmath32 where the @xmath57 contours indicate that the survey has very little sensitivity .
thus , there is no evidence that low - surface brightness galaxies are missing from the  catalogue .
furthermore , in the region occupied by the observed data , the @xmath57 contours are close to vertical indicating that there is little dependence of @xmath57 on surface brightness . the way in which the @xmath58 distribution is modified by including this estimate of the surface brightness dependence is shown by the dotted histogram in the top - left panel of fig .
[ fig : v_vmax ] .
its effect is to increase the mean @xmath58 slightly , suggesting that this estimate perhaps overcorrects for the effect of surface brightness selection .
even so , the change in the estimated luminosity function is negligible as confirmed by the three estimates of the kron j - band luminosity function shown in fig .
[ fig : lf_sb ] .
these are all simple @xmath73 estimates , but with @xmath57 computed either ignoring surface brightness effects or using one of the two methods described above .
these luminosity functions differ negligibly , indicating that no bias is introduced by ignoring surface brightness selection effects .
= 8.5 truecm    = 8.5 truecm    the primary aim of this paper is to determine the present - day j and -band luminosity functions and also the stellar mass function of galaxies . since the  survey spans a range of redshift ( see fig .
[ fig : dn_dz ] ) , we must correct for both the redshifting of the filter bandpass ( k - correction ) and for the effects of galaxy evolution ( e - correction ) . in practice ,
the k and e - corrections at these wavelengths are both small and uncertainties in them have little effect on the estimated luminosity functions .
this is because these infrared bands are not dominated by young stars and also because the  survey does not probe a large range of redshift .
we have chosen to derive individual k and e - corrections for each galaxy using the stellar population synthesis models of bruzual & charlot ( ) . we have taken this approach not because such detailed modelling is necessary to derive robust luminosity functions , but because it enables us to explore the uncertainties in the derived galaxy stellar mass functions , which are , in fact , dominated by uncertainties in the properties of the stellar populations .
the latest models of bruzual & charlot ( ) provide , for a variety of different stellar initial mass functions ( imfs ) , the spectral energy distribution ( sed ) , @xmath74 , of a single population of stars formed at the same time with a single metallicity , as a function of both age , @xmath75 , and metallicity , @xmath76 . we convolve these with an assumed star formation history , @xmath77 , to compute the time - evolving sed of the model galaxy ,
@xmath78 we take account of the effect of dust extinction on the seds using the ferrara ( ) extinction model normalized so that the v - band central face - on optical depth of the milk - way is 10 .
this value corresponds to the mean optical depth of @xmath79 galaxies in the model of cole ( ) which employs the same model of dust extinction .
we assume a typical inclination angle of 60 degrees which yields a net attenuation factor of 0.53 in the v - band and 0.78 in the j - band . by varying the assumed metallicity , @xmath76 , and star formation history ,
we build up a two - dimensional grid of models .
then , for each of these models , we extract tracks of @xmath9  and j@xmath9  colours and stellar mass - to - light ratio as a function of redshift .
our standard set of tracks assumes a cosmological model with @xmath7 , @xmath8 , hubble constant @xmath80 km s@xmath81 mpc@xmath81 , and star formation histories with an exponential form , @xmath82/\tau)$ ] . here ,
@xmath83 is the age of the universe at redshift @xmath54 and the galaxy is assumed to start forming stars at @xmath84 .
for these tracks , we adopt the kennicutt imf ( @xcite ) and include the dust extinction model .
the individual tracks are labelled by a metallicity , @xmath76 , which varies from @xmath85 to @xmath86 and a star formation timescale , @xmath87 , which varies from @xmath88 gyr to @xmath89 gyr .
examples of these tracks are shown in fig .
[ fig : tracks ] , along with the observed redshifts and colours of the  galaxies .
we can see that the infrared j@xmath9  colour depends mainly on metallicity while the @xmath9  colour depends both on metallicity and star formation timescale .
thus , the use of both colours allows a unique track to be selected .
note from the bottom panel that , for all the tracks , the k+e correction is always small for the range of redshift spanned by our data .
we can gauge how robust our results are by varying the assumptions of our model . in particular , we vary the imf , the dust extinction and cosmological models , and include or exclude the evolutionary contribution to the k+e correction . also , we consider power - law star formation histories , @xmath90^{-\gamma}$ ] , as an alternative to the exponential model .
the results are discussed at beginning of section  [ sec : results ] .
the procedure for computing the individual galaxy k+e corrections is straightforward . at the measured redshift of a galaxy
, we find the model whose @xmath9  and j@xmath9  colours most closely match that of the observed galaxy .
having selected the model we then follow it to @xmath62 to predict the galaxy s present - day j and -band luminosities and also its total stellar mass .
we also use the model track to follow its k+e correction to higher redshift in order to compute @xmath56 , the maximum redshift at which this galaxy would have passed the selection criteria for inclusion into the analysis sample .
we use both the simple @xmath73 method and standard maximum likelihood methods to estimate luminosity functions .
we present schechter function fits computed using the sty method ( sandage , tammann & yahil ) and also non - parametric estimates using the stepwise maximum likelihood method ( swml ) of efstathiou , ellis & peterson ( ) .
our implementation of each of these methods is described and tested in norberg ( ) .
the advantage of the maximum likelihood methods is that they are not affected by galaxy clustering ( provided that the galaxy luminosity function is independent of galaxy density ) .
by contrast the @xmath73 method , which makes no assumption about the dependence of the luminosity function with density , is subject to biases produced by density fluctuations .
the two maximum likelihood methods determine the shape of the luminosity function , but not its overall normalization .
we have chosen to normalize the luminosity functions by matching the galaxy number counts of jarrett ( in preparation ) .
these were obtained from a 184  deg@xmath12 area selected to have low stellar density and in which all the galaxy classifications have all been confirmed by eye .
the counts are reproduced in fig .
[ fig : counts ] . by using the same 7  arcsec aperture magnitudes as jarrett ( in preparation ) and scaling the galaxy counts in our redshift survey
, we deduce that the effective area of our redshift catalogue is @xmath91  deg@xmath12 .
note that normalizing in this way by - passes the problem of whether or not some fraction of the missed objects are stars .
[ fig : counts ] also shows the kron and extrapolated magnitude j and  counts of the   redshift survey . in the lower panel , these counts are seen to be in agreement with the published k - band counts of gardner ( ) and glazebrook ( ) .
we also checked the normalization using the following independent estimate of the effective solid angle of the redshift survey .
for galaxies in the  parent catalogue brighter than @xmath92 , we computed the fraction that have both measured redshifts and match a  galaxy . for a faint @xmath93
this fraction is small as the  catalogue is much deeper than the  catalogue , but as @xmath93 is made brighter , the fraction asymptotes to the fraction of the area of the  parent catalogue covered by the joint redshift survey . by this method
we estimate that the effective area of our redshift catalogue is @xmath94  deg@xmath12 , which is in good agreement with the estimate from the counts of jarrett .
it should be noted that for neither of these estimates of the effective survey area do the quoted uncertainties take account of variations in the number counts due to large scale structure . to estimate the expected variation in the galaxy number counts within the combined -  survey due to large scale structure we constructed an ensemble of mock catalogues from the @xmath95cdm hubble volume simulation of the virgo consortium ( evrard ; evrard in preparation ; http://www.physics.lsa.umich.edu/hubble-volume ) .
mock catalogues constructed from the virgo hubble volume simulations ( baugh in preparation ) can be found at http://star-www.dur.ac.uk/@xmath96cole/mocks/hubble.html .
we simply took these catalogues and sampled them to the depth of  over a solid angle of @xmath97  deg@xmath12 .
to this magnitude limit we found an rms variation in the number of galaxies of 15% .
we took this to be a realistic estimate of the uncertainty in the  number counts and propagated this error through when computing the error on the normalization of the luminosity function .
[ cols= " < , > , > , > , > , > , > , > , > , > , > , > , > , > , > , > , > , > , > " , ]     [ tab : madau ]    it is interesting to compare our values with what is inferred by integrating the observational estimates of the mean star formation history of the universe .
[ fig : madau ] shows observational estimates for one particular choice of cosmology and imf and illustrates how the rates are sensitive to the assumed dust extinction . by fitting a smooth curve through these estimates
, we can calculate the mass of stars formed by the present day and how this depends on the imf and assumed dust extinction . the upper smooth curve shown in fig .
[ fig : madau ] is of the form @xmath98yr@xmath81mpc@xmath3 , where @xmath99 .
the data points uncorrected for dust extinction are fit with @xmath100 . as for our estimates
above , we assume that no mass goes into forming brown dwarfs and multiply the star formation rate by @xmath101@xmath9@xmath102 , where @xmath102 is the recycled fraction , so as to form an estimate of the mass locked up in stars .
values of @xmath103 estimated in this way are listed in table  [ tab : madau ] .
the values in this table are for an @xmath7 , @xmath8 cosmology , but they are insensitive to this choice .
they depend slightly on the assumed metallicity of the stellar population and would be 10% lower if half solar , rather than solar metallicity were assumed .
note that the @xmath104 values inferred from the star formation history of the universe scale differently with the assumed hubble constant than those inferred above from the ir luminosity functions . for @xmath105 our estimates from   become @xmath106 for the kennicutt imf and @xmath107 for the salpeter imf .
comparison with table  [ tab : madau ] shows that these values are consistent with those inferred from the cosmic star formation history only if the dust correction assumed in the latter is modest , e(b - v)@xmath108 .
this value is 50% smaller than the value preferred by steidel ( ) .
the new generation of very large surveys currently underway make it possible to characterize the galaxy population with unprecedented accuracy . in this paper , we have combined two such large surveys , the infrared imaging 2mass and the 2df galaxy redshift survey and  [ tab : smfswml ] are available at http://star-www.dur.ac.uk/@xmath96cole/2dfgrs-2mass . ] to obtain a complete dataset which is more than an order of magnitude larger than previous datasets used for statistical studies of the near - infrared properties of the local galaxy population .
we have used this combined catalogue to derive the most precise estimates to date of the galaxy j and -band luminosity functions and of the galaxy stellar mass function .    characterizing the near - infrared properties of galaxies
offers several advantages .
firstly , the near - infrared light is dominated by established , old stellar populations rather than by the recent star formation activity that dominates the blue light .
thus , the j and k - band luminosity functions reflect the integrated star formation history of a galaxy and , as a result , provide particularly important diagnostics of the processes of galaxy formation . for the same reason , the distribution of stellar mass in galaxies the galaxy stellar mass function can be derived from the near - infrared luminosities in a relatively straightforward way , with only a weak model dependence .
finally , corrections for dust extinction as well as k - corrections are much smaller in the near - infrared than in the optical .
due to the size of our sample , our determination of the j- and -band galaxy luminosity functions have , for the most part , smaller statistical errors than previous estimates .
furthermore , since our sample is infrared - selected , our estimates are free from any potential biases that might affect infrared luminosity functions derived from optically - selected samples .
we find that the j- and -band galaxy luminosity functions are fairly well described by schechter functions , although there is some evidence for an excess of bright galaxies relative to the best - fit schechter functional form . in general , the swml estimates are a truer representation of the luminosity functions .
our k - band estimates are in overall agreement with most previous determinations , but have smaller statistical errors .
the exception is the k - band luminosity function inferred from the near infrared sdds photometry ( @xcite ) .
the difference between the k - band luminosity function we infer from their data and our own estimate is too large to be explained by photometric differences .
the difference between the two estimates is better described by a difference in overall number density of a factor of @xmath109 .
a similar discrepancy is seen in the -band between the sdss and luminosity function estimates ( see @xcite and @xcite ) .
the suspicion is that the uncertainty in the overall normalization of the luminosity functions induced by large - scale structure within the large , but finite , survey volumes could be to blame .
however , the errors that we quote for the -luminosity functions already include an estimate of this sampling uncertainty as derived from realistic mock catalogues . a similar exercise for a catalogue with the same area and depth as that of blanton ( ) indicates that the required overdensity of a factor of @xmath109 is unlikely .
so probably there is more than one contributory factor at work and the hope is that these will be identified as the surveys progress .    using our j - band luminosity function , @xmath9  and j@xmath9  colours and simple galaxy evolutionary tracks
, we have obtained the first estimate of the galactic stellar mass function derived directly from near - infrared data .
we find that this mass function is also fairly well described by a schechter form .
an integral over the stellar mass function gives @xmath104 , the universal mass density locked up in luminous stars and stellar remnants , expressed in terms of the critical density .
@xmath104 is a key component of the overall inventory of baryons in the universe .
an accurate determination of this quantity is essential for detailed comparisons with other quantities of cosmological interest such as the total baryonic mass , @xmath110 , inferred from big bang nucleosynthesis considerations ( e.g. burles & tytler 1998 ) , the cosmic star formation rate ( e.g. steidel ) , and the cosmic evolution of the gas content of the universe ( storrie - lombardi & wolfe 2000 ) .
the statistical uncertainty in our estimate of @xmath104 is about 15% , several times smaller than the best previous determination by salucci & persic ( ) . in fact , the errors in @xmath104 are dominated by systematic uncertainties associated with the choice of stellar imf in the galaxy evolution model .
a kennicutt imf gives @xmath111 while a salpeter imf gives @xmath11 .
for @xmath105 , these values correspond to less than 11% of the baryonic mass inferred from big bang nucleosynthesis ( @xcite ) .
our values of @xmath104 today are only consistent with recent determinations of the integrated cosmic star formation if the correction for dust extinction is modest .
we thank john lucey for the initial suggestion to look at the  database .
we thank tom jarrett for advice and detailed explanations of the 2mass data reduction procedure .
we also thank ian smail for help in manipulating dss images , jon loveday for supplying his luminosity function data in electronic form and simon white for useful suggestions .
we are grateful to the virgo consortium for allowing us to use the hubble volume simulation data prior to publication.the redshift data used here were obtained with the 2-degree field facility on the 3.9 m angloaustralian telescope ( aat ) .
we thank all those involved in the smooth running and continued success of the 2df and the aat .
smc acknowledges a pparc advanced fellowship , iprn an snsf and ors studentship and csf a leverhulme research fellowship .
the 2df galaxy redshift survey was made possible through the dedicated efforts of the staff at the anglo - australian observatory , both in creating the two - degree field facility and supporting it on the telescope .


within the last decade an increasing interest in molecular electronics has developed , with the expectation of realising molecular diodes and transistors .
this is based on the progress in manipulation techniques , which now allow the controlled attachment of atomic scale structures like molecules to mesoscopic leads . with these new devices
one is able to determine the conductance properties of molecular structures . explaining and predicting the electronic behaviour of such devices is an essential step towards their design and use as nano - scale electronic circuits .    to this end a number of theoretical studies have been performed with the aim of reproducing measured @xmath0 characteristics .
these studies differ in the way they take the electronic levels of the molecules , their modification by the coupling to the leads and the change of electrostatic potential due to bias into account .
semi - empiric methods@xcite have been used , as well as first principles techniques,@xcite the latter being restricted to systems of moderate size .
the wide range of experimentally observed behaviour ( see section [ sec : experiment ] ) suggests , that not only the structure of the molecule , but also the details of the device fabrication affect the conduction properties of molecular devices .
the crucial step is the deposition of molecules onto the surface of the lead .
as this is achieved by self - assembly the amount of adsorbed molecules and their individual positions can not be exactly controlled and therefore remains unknown .
a satisfactory understanding of the interplay between geometrical alignment of the molecules and measured conductance properties has thus not yet been achieved ( for a recent review see e.g.  nitzan et al.@xcite ) .    in this paper
we study the way in which changing the geometrical alignment of the mono - layer has an influence on the conduction properties of a molecular device . in
so doing we can rule out a number of explanations which have previously been considered@xcite to explain the occurrence of additional structure in the conductance - voltage ( @xmath1 ) characteristics .
the outline is the following : first we summarise some of the recent experimental findings .
then the method we use ( based on esqc@xcite ) for calculating the conductance properties of molecular devices , is discussed .
calculations for the conductance of a sam , being build of 1,4-phenylene diisocyanide ( pdi ) and sandwiched between gold leads are then presented .
the results for three qualitatively different geometrical constitutions of the mono - layer are compared with experimental data . by this
we can rule out all but one and conclude that the only geometrical alignment , which gives rise to several peaks in the conductance curve , is a mono - layer with additional molecules perturbing the periodicity .
the devices built to study conductance properties of molecular structures differ not only in amount and chemical structure of the molecules in use but also in the way these are attached to metallic or semi - conducting leads
. single or few molecules are accessible in mechanically controllable break junctions ( mcb ) and with the scanning tunneling microscope ( stm ) .
many molecules are involved in sandwiched self - assembled mono - layer ( sam ) experiments .
the observed properties depend on the exact geometry of the device .
the conductance differs in orders of magnitude and the qualitative voltage dependency of the current ranges from simple ohmic behaviour to negative differential resistance ( ndr).@xcite    in the past reed et al.@xcite have measured the electrical conductance of a self - assembled molecular mono - layer bridging an mcb at room temperature .
molecules of 1,4-benzene dithiol ( i.e.  having two thiol end - groups , which are known to couple strongly to au - atoms ) were used and the @xmath1 characteristic was found to be symmetric with one peak in the voltage range of @xmath2v . they measured a current of the order of @xmath3na at a bias voltage of @xmath4v , which they claim is produced by transport through one single active molecule .
reichert et al.@xcite also used an mcb with molecules having two thiol end groups , but being considerably longer .
the measured current amplitude was about @xmath5na at @xmath6v , i.e.  although the molecule was more than twice as long , the current was ten times larger .    with a different setup , where a sam is sandwiched between two metallic leads ,
chen et al.@xcite have found negative differential resistance ( ndr ) , namely one peak at @xmath4v in the @xmath0 curve .
the molecule under investigation had one thiol end - group only and was attached to au - leads at both ends .
the measurements were taken at room temperature and the measured current maximum was of the order of @xmath6na .
only recently , sandwiched sam devices at @xmath7k were studied,@xcite where a benzene ring with two cyanide instead of thiol end - groups was used .
the measurements exhibited currents of the order of @xmath8na .
the @xmath1 characteristic for this molecule revealed more structure , in form of three to five peaks within a voltage range of @xmath6v .
such a behaviour was not observed with previous devices containing other molecules .
in the literature there has been presented quite a number of techniques to calculate non - equilibrium electronic transport through molecular systems attached to mesoscopic leads .
usually the landauer formalism is applied , which describes current as elastic electron transmission and therefore requires the transmission function @xmath9 . to this end
one needs a framework that allows for a description of the molecular device on the atomic level .
this involves not only the molecules themselves , but also the surface and bulk region of the leads .
quantum chemistry provides such a framework , and one can choose the level of theory according to the size of the system under consideration and the computational effort one is willing to spend .    using a quantum chemistry method , the transmission function can be obtained from an effective one - particle hamiltonian , which is an appropriate description for strong coupling of the molecules to the leads ( as in the case of covalent binding ) .
the methods differ in the generation of the one - particle hamiltonian , which might be based on semi - empirical grounds @xcite or on first - principles and self - consistent techniques .
@xcite    a different approach,@xcite taking many - particle effects explicitly into account , uses a master equation with transition rates calculated perturbatively using the golden rule .
this approach is appropriate for weak coupling .
we use the landauer formalism , as the molecules are assumed to be chemically bonded to the gold contacts ( i.e.  strong coupling ) , together with the semi - empiric extended hckel quantum - chemistry method esqc.@xcite the molecular structure is optimised@xcite beforehand .
this approach , though limited as compared to more sophisticated quantum chemistry methods , is yet justified , because we want to gain a qualitative understanding of a many molecule experiment , which can not be described by first - principle techniques , as the number of atoms involved is beyond the practical limitations of to - date computer resources .      according to the landauer formula , current along a defect region is the result of electron transmission from the source to the drain lead .
this is described by the transmission function @xmath9 and net electron transfer happens at all those energies , where there are more states occupied in the source lead than in the drain .
the difference in occupation is a result of the applied bias @xmath10 , which raises the chemical potential @xmath11 of the source lead above the one of the drain lead @xmath12 and thus changes the distribution function @xmath13 . @xmath14
the landauer formula is valid under the condition that transport is coherent across the molecule , which is plausible as the typical mean - free path of electrons within metals is of the order of @xmath5 nm , while the molecular gap between source and drain lead is only about @xmath15 nm in length .
furthermore it is assumed that the coupling of the leads to macroscopic reservoirs is reflection - less , and that the tunneling electrons equilibrate only deep inside the leads .
this ensures that the distribution function for incoming electrons in a lead can be taken to be spatially constant , even close to the molecular region .
( colour online ) partitioning of the system into three parts : the two semi - infinite leads @xmath16 ( surrounded by boxes ) and the molecular region @xmath17 between them.,scaledwidth=50.0% ]    the system is formally partitioned into three regions @xmath18 , two of them ( @xmath16 ) containing the semi - infinite leads , the third one ( @xmath19 ) being the finite region containing all molecules as well as a few surface layers of each lead ( see fig .  [
fig : system ] ) .
we use periodic boundary conditions in the directions perpendicular to the surface normal .    by a tight binding approximation
, the infinite - dimensional hamiltonian of the entire system can be composed of quantum chemical one - particle block hamiltonians of finite dimension : @xmath20 the first summation describes the isolated molecular region , by an on - site energy @xmath21 and a hopping term . the indices @xmath22 and @xmath23 run over the orbital basis set within that region .
the next two summations describe the isolated leads , labelled by @xmath24 .
layer by layer , starting with the surface layer @xmath25 , the first term accounts for intra - layer interactions , while the second one describes the interaction between layers .
( the size of each layer is chosen such that only adjacent layers have non - zero interaction .
it therefore depends on the details of the tight - binding approximation . )
finally the last term describes the coupling between the molecular region and each lead .
note that only the first layer @xmath25 contributes to that term and that there is no interaction between different leads .
these are only formal restrictions , as parts of each lead can be included into the molecular region .
the determination of the transmission function involves two steps .
first the conduction properties of the isolated leads have to be calculated .
thereby each lead will be decomposed into conducting and non - conducting in- and out - going channels .
these correspond to propagating and evanescent solutions moving in one of two possible directions respectively . in a second step
, the channels are connected to each other via the molecular region .
this is described by the scattering matrix and the transmission function is finally obtained by summing up the contribution from each channel .
the calculation can be performed either using green s function techniques@xcite or equivalently@xcite using esqc@xcite , which is a scattering - matrix approach .
we present the details of the calculation in the second scheme , as individual contributions from each channel to the transmission function can then be easily studied .
first we will restrict our attention to the semi - infinite lead hamiltonians , which do not have to be identical .
the hamiltonians of equation for one lead @xmath26 are layer independent , if one assumes periodicity , i.e.  @xmath27 and @xmath28 . using bloch s theorem
one can reduce the infinite dimensional system of equations to an @xmath29-matrix equation ( @xmath30 being the number of orbital basis functions in one layer ) @xmath31 with @xmath32 , @xmath33 , and @xmath34 is the overlap matrix between orbitals in layer @xmath35 and layer @xmath36 for cases when one does not deal with an orthonormal basis set ( otherwise @xmath37 ) .
with @xmath38 we denote the lattice spacing . defining @xmath39
one can easily see that equation is an @xmath29 quadratic eigenvalue equation .
it can be transformed into a @xmath40 linear eigenvalue equation : @xmath41    \end{aligned}\ ] ] ( where we have dropped the energy dependency for notational ease ) .
this layer - to - layer propagator @xmath42 reduces the problem of finding solutions for the entire isolated lead hamiltonian to specifying the wave function coefficients at two adjacent layers @xmath43 and @xmath44 only .
this is due to the fact that given these coefficients any other pair of coefficients @xmath45 can now be obtained via : @xmath46 all possible solutions at energy @xmath47 can be decomposed into independent channels , by solving for the eigenvalues of equation .
these eigenvalues come in pairs such that for each eigenvalue @xmath48 , there exists a corresponding eigenvalue @xmath49 satisfying the relation @xmath50 , as can be seen by transposing equation .
eigenvalues with @xmath51 , i.e.  complex @xmath52 , belong to exponentially diverging solutions ( see equation ) .
these are of course non physical , as long as the lead is infinite . in semi - infinite leads
however ( which we are dealing with ) , exponentially decaying coefficients at the boundary will contribute to the surface wave function and must not be neglected .
the contribution from a single channel to the net current can not directly be seen from equation .
it depends on the current density associated with a solution to the schrdinger equation @xmath53 and is obtained via the continuity equation .
the probability amplitude @xmath54 for a stationary solution is constant in time @xmath55 because @xmath56 and @xmath57 are hermitian . for
the probability amplitude at all layers between @xmath58 and @xmath59 one therefore has @xmath60      \langle l_2 , l_2 + 1| \gamma\rangle \\\nonumber      & -       \langle\gamma |l_1 - 1,l_1\rangle        \frac{i}{\hbar }        \left[\begin{array}{cc }          0           & -h\\          h^\dagger & 0        \end{array}\right ]      \langle l_1 - 1 , l_1| \gamma\rangle ,    \end{aligned}\ ] ] with the projectors @xmath61 .
this gives rise to the definition of the current operator @xmath62 for layer @xmath35 as @xmath63      \langle l , l+1|.\ ] ] now let both @xmath64 and @xmath65 be solutions at fixed energy @xmath47 with the eigenvalues @xmath66 and @xmath67 respectively . because the expectation value for @xmath62 is layer independent ( equation [ eq : currentoperator2 ] ) one has : @xmath68 this equation describes the connection between the current properties of a solution @xmath64 and its eigenvalue @xmath69 .
we summarise the results of a detailed analysis of this equation , which is given in appendix [ ap : current ] .
each channel @xmath70 can be assigned a current value @xmath71 , defined as @xmath72 where we have used the layer independence of @xmath62 in simply writing @xmath73 .
channels with eigenvalue modulus @xmath74 , i.e.   evanescent waves have zero current value .
they therefore do not contribute to the current .
( yet they are important at the surface , as already mentioned above . )
only channels with an eigenvalue of modulus one ( @xmath75 ) contribute to the current .
the sign of @xmath71 determines the direction of charge transport .
therefore the eigenvalues can be sorted into incoming and outgoing solutions ( @xmath48 and @xmath49 respectively ) , according to the following scheme : if @xmath76 it represents an incoming evanescent solution , if @xmath77 then it is outgoing evanescent . only if @xmath78 it belongs to a propagating solution ,
which is incoming for @xmath79 and outgoing otherwise .
we now define @xmath80 and @xmath81 as the two @xmath29 diagonal matrices composed of all incoming and outgoing eigenvalues @xmath82 .
the @xmath83-matrix @xmath84 , which diagonalises @xmath42 : @xmath85,\ ] ] has the following quadratic block form : @xmath86.\ ] ] after this transformation into the diagonal basis of the propagator , we can easily obtain all physically relevant solutions of the infinite lead by specifying the amplitudes of all propagating waves at one lattice site .      up to now
, we have considered the isolated leads only .
these are now assumed to be each coupled to the molecular defect region and thereby indirectly coupled to one another .
we are interested in stationary solutions , which consist of an incoming propagating wave in one lead , being scattered among all the accessible outgoing channels ( propagating and evanescent ones ) .
this information is contained in the scattering matrix @xmath57 @xmath87 which determines the wave amplitudes of all outgoing waves @xmath88 given the incoming ones @xmath89 .
( 0,0 )     ( 2499,1221)(439,-799 ) ( 631,254)@xmath90 ( 721,-736)@xmath91 ( 2476,-736)@xmath92 ( 2476,254)@xmath93    in appendix [ ap : smatrix ] it is shown , that @xmath57 is of the form @xmath94    it is important to notice , that the scattering matrix is always quadratic , because in each lead there are the same amount of incoming and outgoing channels .
this is opposed to the transfer matrix @xmath95 , which determines the amplitudes of in- and outgoing waves @xmath96 in one lead given the in- and outgoing waves @xmath97 of a second lead : @xmath98 this matrix is quadratic only if both leads have the same number of channels .
it then is of the form@xcite @xmath99 and the relation to the scattering matrix is@xcite @xmath100 methods calculating the scattering matrix via the transfer matrix@xcite fail , if two types of leads are used , because @xmath101 is then no longer quadratic and can not be inverted .
therefore one commonly takes source and drain lead to be identically constituted .
but even in such cases , the method becomes numerically unstable , with increasing distance between the molecular region and one lead .
this is because the matrix elements of @xmath101 and @xmath102 ( in equation ) diverge exponentially , with increasing lead separation .
taking the inverse of @xmath101 is therefore a numerically critical procedure .
both these problems are avoided by the direct calculation of the scattering matrix , which we present in appendix [ ap : smatrix ] .
this calculation is well defined without any restrictions to the number of leads and their composition .
that means that it is not necessary to restrict to identical leads .
furthermore it allows a numerically stable determination of the scattering matrix , even for weak coupling .
the transmission function is the sum over the contributions from each combination of incoming channel in the source lead to outgoing channel in the drain lead : @xmath103 .
the relation between scattering matrix @xmath57 and these transmission function elements is @xmath104 where @xmath105 is the lower left block of @xmath57 as defined in equation .
the weighting with velocity factors comes about because the scattering matrix @xmath106 does not relate current densities , but wave amplitudes .
the current densities are obtained from these wave amplitudes by multiplication with the corresponding velocity factor @xmath107 .
the factor @xmath71 in the denominator normalises the transmission function to be exactly one for perfect transmission .
if @xmath108 then @xmath109 , because incoming evanescent waves have zero amplitude at the surface .
the total current is made up of the contribution from each @xmath52-state@xcite : @xmath110 here the summation over @xmath52 has been transformed into an integral over @xmath47 with a factor of 2 accounting for spin @xmath111 .
low temperature experiments with pdi sam s sandwiched between two metallic leads show several peaks in the @xmath1-diagram.@xcite the typical voltage differences of these peaks are in the range of @xmath112 ( i.e.  there are about 5 peaks within @xmath113 and @xmath114 ) .
the commonly adopted explanation for the occurrence of such peaks is the following .
each molecular orbital that enters the energy window , which is opened by the applied voltage , enables resonant tunneling .
this increases the conductance and therefore results in a peak within the @xmath1-diagram .    typically , the energy gap between molecular orbitals is in the range of @xmath115 . in other words ,
for applied voltages up to @xmath114 there should be only a single accessible orbital per molecule , giving rise to only a single peak in the @xmath1-diagram . therefore the following question arises : are there geometrical alignments of the molecules such that the additional peaks in the @xmath1 diagram can also be explained by resonant tunneling through molecular orbitals ?      during the device fabrication , the step under least experimental control is the adsorption of the molecules onto the leads .
therefore the exact geometrical alignment of the molecular sam and , at least in the sandwich geometry , also the atomic shape of the top metallic lead , is not exactly known .
one therefore has to expect not only one specific but rather quite a variety of molecular alignments to be produced .
as one is interested in the conduction properties of the resulting device , it is important to understand the influence of each type of geometrical alignment to the transmission function .    to this end , we have investigated three such possible alignments , which will be discussed separately .    1 .   by adding metal - atoms on top of the molecular mono - layer in a random way
, there might occur metallic clusters on top of the mono - layer .
these affect the electronic configuration of the molecules individually , and might therefore have an influence to the transmission function . 2 .   in a sam experiment
, there is not just one molecule , but rather a few hundred molecules involved .
if the contribution to the transmission function was different for each molecule , then @xmath9 would change qualitatively , with a change in size of the mono - layer .
we therefore analyse how the transmission functions depends on the number of molecules involved .
3 .   if the molecular mono - layer is not strictly periodic , then there will be defects .
for example can the distance between two molecules be reduced , such that inter - molecular bonds can be build .
each of these defects will have a specific electronic structure and will therefore influence the transmission function .      in the sandwich geometry ,
first the bottom metallic lead is created .
then the molecular mono - layer is adsorbed on top of it by self - assembly .
finally the top metallic lead is build upon the molecular mono - layer .
the exact shape of neither metallic surface is known and may be anything but flat and regular .    it is likely that the surface atoms of the top metallic lead build up clusters on top of the molecular layer ( as for example in fig .
[ fig : singleandcluster ] b ) . which influence do they have on the electronic configuration of the molecule they are in contact with ? and do the clusters act as small molecules with new electronic levels ?
the influence of an au cluster on the molecular electronic structure is twofold .
first it introduces new electronic levels , and second the existing molecular electronic levels will be shifted , by an amount which depends on the strength of the coupling between cluster and molecule .
the latter effect will be observed as a shifted peak in the transmission function , only if the coupling between cluster and molecule is different to the coupling between top electrode and molecule . for clusters similar to the one shown in fig .
[ fig : singleandcluster ] ( b ) , this is however not the case .
the energetic peak positions are identical , as can be seen in fig .
[ fig : singleandcluster ] ( c ) .    furthermore , there are no additional peaks , which one might have expected because of the additional electronic levels of the cluster .
the explanation of their absence is the following : an electronic level gives rise to a peak in the transmission function only , if the corresponding orbital wave function overlaps with both the top and bottom electrode .
the overlap with the electrode the cluster is attached to ( say top electrode ) is of course large .
the overlap with the bottom electrode consists of two parts . the direct overlap and the indirect overlap via the molecule .
the direct overlap is negligible due to the large spatial separation .
the indirect overlap depends on the molecular orbital wave function .
if the energy of the cluster level does not coincide with a molecular energy level , then there is no indirect overlap . only if two levels coincide , the indirect coupling is large , but in that case , there already exists a transmission peak due to the molecule itself .
( 0.98,0.72 ) ( 0,0 )   ( a ) ( colour online ) structure for a molecule without cluster .
( b ) ( colour online ) sturcture for a molecule with a gold cluster on top .
( c ) transmission function @xmath9 for both structures .
energy scale is relative to the homo - lumo gap , such that @xmath116 corresponds to the middle of the gap .
, title="fig : " ] ( 0,0.72 )    ( 0.98,0.72 ) ( 0,0 )   ( a ) ( colour online ) structure for a molecule without cluster .
( b ) ( colour online ) sturcture for a molecule with a gold cluster on top .
( c ) transmission function @xmath9 for both structures .
energy scale is relative to the homo - lumo gap , such that @xmath116 corresponds to the middle of the gap .
, title="fig : " ] ( 0,0.72 )     +    ( 0.98,0.98 ) ( 0,0 )   ( a ) ( colour online ) structure for a molecule without cluster .
( b ) ( colour online ) sturcture for a molecule with a gold cluster on top .
( c ) transmission function @xmath9 for both structures .
energy scale is relative to the homo - lumo gap , such that @xmath116 corresponds to the middle of the gap .
, title="fig : " ] ( 0,0.98 )    therefore if transmission is already suppressed by the molecule ( at all off - resonant energies ) , it can either be further reduced by off - resonant tunneling through the cluster , or it can ( at best ) be left unchanged by resonant tunneling through the cluster . under no circumstances can transmission , once suppressed by the molecule , be afterwards increased by the cluster .
this in turn means , that metallic clusters can not give rise to additional peaks in the transmission spectrum .
( colour online ) the structure for one , two , three , and four molecules adsorbed within an au-@xmath117 super - cell .
this setup was used to test the sum rule.,title="fig : " ]   ( colour online ) the structure for one , two , three , and four molecules adsorbed within an au-@xmath117 super - cell .
this setup was used to test the sum rule.,title="fig : " ]   ( colour online ) the structure for one , two , three , and four molecules adsorbed within an au-@xmath117 super - cell .
this setup was used to test the sum rule.,title="fig : " ]   ( colour online ) the structure for one , two , three , and four molecules adsorbed within an au-@xmath117 super - cell .
this setup was used to test the sum rule.,title="fig : " ]    how do we expect the transmission function @xmath118 for @xmath22 periodically arranged molecules to look like ?
as long as the inter - molecular interactions are small ( compared to the intra - molecular ones ) the molecular levels of each molecule will not be significantly changed .
furthermore as the mono - layer consists of only one kind of molecule , all of them will have the same electronic structure .
therefore we expect each molecule to contribute the same amount to the transmission function : @xmath119 , where @xmath22 runs over all @xmath120 adsorbed molecules .
( 1.0,1.0 ) ( 0,0 )   ( a ) transmission function for one , two , three , and four pdi molecules .
( b ) when plotted against each other , the transmission functions reveal a linear relationship : @xmath121 ( bottom).,title="fig : " ] ( 0,1.0 )    ( 1,1 ) ( 0,0 )   ( a ) transmission function for one , two , three , and four pdi molecules .
( b ) when plotted against each other , the transmission functions reveal a linear relationship : @xmath121 ( bottom).,title="fig : " ] ( 0,1 )    we calculated the transmission function for @xmath120=1 to 4 molecules within an au super - cell of size @xmath117 ( the structures are shown in fig .
[ fig : onetofourmolecules ] ) .
the distance between the molecules is chosen to be a multiple of the closest au - au separation @xmath122 ( @xmath123 , with @xmath124 ) . to our knowledge ,
the parameters of the pdi - sam mono - layer have never been determined experimentally , which is why we have to assume the above values .
however stm studies@xcite and also theoretical calculations@xcite have been performed for alkanethiol mono - layers , and these parameters motivated our choice .
independent from the number of molecules present , the transmission functions have the same amount of peaks , at identical energetic positions ( see fig .  [
fig : fourmoleculespectrum ] a ) .
this result is also obtained for all larger distances of the molecules , where the inter - molecular interaction is even smaller . furthermore ,
the sum rule is indeed fulfilled , as shown in fig .
[ fig : fourmoleculespectrum ] ( b ) .
each @xmath118 is plotted against @xmath125 for @xmath126 .
the calculated transmission values ( 300 discrete values each ) clearly show a linear correlation .
the straight lines are linear fits to the data , and their slope does very well agree with the theoretically expected value of @xmath127 .
the deviation is below @xmath128 , as can be seen in table [ tab : fitparameters ] , where we summarise all the fitted values for @xmath129 .
.[tab : fitparameters ] the fitted values for @xmath130 together with their deviation from the theoretical value @xmath131 and a measure for the quality of the fit @xmath111 , where @xmath132 for @xmath133 discrete energy values @xmath134 . [ cols="<,<,^,^,^",options="header " , ]     we conclude the following : a mono - layer , where the inter - molecular distance is large enough to not let inter - molecular interactions play a significant role , has the same number of distinct electronic levels as a single molecule .
these levels are then highly degenerate .
a @xmath1-diagram will therefore have the same number of peaks .
only the net current will be increased by a factor @xmath130 compared to the single molecule case .
the mere fact , that one deals with a mono - layer instead of a single molecule does not imply that the transmission function changes qualitatively .
we have seen , that one does not observe additional peaks in the transmission function , as long as the inter - molecular influence is small . and this is the case for distances which occur in typical sam structures.@xcite we now investigate cases , where the molecular interactions are not negligible .
this occurs for example , when the periodic structure of the mono - layer is perturbed by an additional molecule , such that a molecular cluster is formed .
it is sufficient to study the transmission function of an isolated cluster only , because we have already seen that molecules in the periodic sam arrangement do not influence each other .
the sum of the transmission function for the periodic sam and the transmission function for the molecular cluster is , due to the sum rule , the total transmission function for defect and sam .
we study the influence of a shorter distance between two , three , and four molecules on the transmission spectrum and relate it to the discrete energies of the isolated molecules .
the molecules are now separated by @xmath135 which corresponds to the au - au atom spacing . the atomic structure for this calculation
is shown in fig .
[ fig : molecularinteraction ] ( a ) , the resulting transmission functions in fig .
[ fig : molecularinteraction ] ( b ) and ( c ) .
( 1.0,0.25 ) ( 0,0 )   ( a ) ( colour online ) two , three , and four molecules with a shorter inter - molecular distance .
( b ) the transmission functions for two and three molecules .
( c ) the transmission functions for three and four molecules .
in contrast to all previous cases , the peaks are shifted with respect to each other and there are also additional peaks .
these changes are due to the increase in inter - molecular interaction , which alters the electronic levels.,title="fig : " ] ( 0.33,0 )   ( a ) ( colour online ) two , three , and four molecules with a shorter inter - molecular distance .
( b ) the transmission functions for two and three molecules .
( c ) the transmission functions for three and four molecules .
in contrast to all previous cases , the peaks are shifted with respect to each other and there are also additional peaks .
these changes are due to the increase in inter - molecular interaction , which alters the electronic levels.,title="fig : " ] ( 0.66,0 )   ( a ) ( colour online ) two , three , and four molecules with a shorter inter - molecular distance .
( b ) the transmission functions for two and three molecules .
( c ) the transmission functions for three and four molecules .
in contrast to all previous cases , the peaks are shifted with respect to each other and there are also additional peaks .
these changes are due to the increase in inter - molecular interaction , which alters the electronic levels.,title="fig : " ] ( 0,0.25 )    ( 1.0,1.0 ) ( 0,0 )   ( a ) ( colour online ) two , three , and four molecules with a shorter inter - molecular distance .
( b ) the transmission functions for two and three molecules .
( c ) the transmission functions for three and four molecules .
in contrast to all previous cases , the peaks are shifted with respect to each other and there are also additional peaks .
these changes are due to the increase in inter - molecular interaction , which alters the electronic levels.,title="fig : " ] ( 0,1.0 )    ( 1.0,1.0 ) ( 0,0 )   ( a ) ( colour online ) two , three , and four molecules with a shorter inter - molecular distance .
( b ) the transmission functions for two and three molecules .
( c ) the transmission functions for three and four molecules .
in contrast to all previous cases , the peaks are shifted with respect to each other and there are also additional peaks .
these changes are due to the increase in inter - molecular interaction , which alters the electronic levels.,title="fig : " ] ( 0,1.0 )    by reducing the molecular separation from @xmath136  to @xmath137 , the transmission function qualitatively changes .
the number of peaks roughly doubles .
the new peak positions are different from the ones we have obtained in the previous calculations . and
this time , the peak positions do depend on the number of molecules involved .
this is an important point , because if there are several molecular clusters with different molecular distances , then they all give rise to peaks at different energy values .
the resulting transmission function is the sum of the individual functions and will thus contain far more peaks , than the transmission function for the non perturbed periodic layer .
we now show that the new peaks are a result of the increase in molecular interaction due to the decrease in spatial separation . for non - interacting molecules ,
the molecular energies are identical and therefore degenerate .
an interaction between molecules breaks this degeneracy and therefore new energy levels occur . by performing a diagonalisation of
the molecular hamiltonian ( without leads ) one can determine the levels of the molecular cluster .    in fig .
[ fig : moleculesandenergies ] we have again plotted the transmission function for three and four molecules , this time together with the discrete energy levels of the corresponding molecular cluster . the inset is identical to fig .
[ fig : molecularinteraction ] ( c ) , while the plot itself is a magnification , to better resolve the discrete energy levels ( which are shown as points along the transmission function ) .
each of the transmission peaks is related to at least one discrete energy value .
but not each energy value can be related to a peak in the transmission function .
why is that ?
the discrete energies can only give rise to new peaks in the transmission function , if they are not suppressed by a weak coupling to one of the leads .
all levels which are not related to any peak belong to this category .
if the position of the peak is shifted away from a corresponding energy level , then this is due to the coupling between molecules and leads .
this coupling is absent in the diagonalisation of the molecular hamiltonian , but present in the calculation of the transmission function .
magnification of the transmission functions for three and four closely spaced molecules .
additionally the discrete energy levels of the system without leads are plotted as points along the transmission function . to each peak there
belongs at least one discrete energy level .
a detailed discussion is given in the main text .
inset : transmission function ( original scale ) for three and four closely spaced molecules ( identical to fig .
[ fig : molecularinteraction ] c ) . ]
finally we show , that the additional peak structure in the transmission function for a scenario with an increased inter - molecular interaction gives rise to a number of steps in the @xmath0-curve .
figure [ fig : defectcurrent ] contains an @xmath0 calculation for a molecular structure containing all three molecular clusters shown in fig .
[ fig : molecularinteraction ] ( a ) . in this calculation
, the bias voltage @xmath138 enters as a shift of the fermi levels for source and drain lead : @xmath139 .
the molecular energy has been set to @xmath140 , where @xmath141 is the zero bias displacement of the molecular levels and @xmath142 , because of the symmetric coupling to the leads .
compared to the experiments@xcite the number of steps in the @xmath0curve is well reproduced by our calculation .
the obtained current is at least one order of magnitude larger than the experimental values.@xcite this is a phenomenon common to all theoretical methods based on the landauer formula.@xcite a satisfactory explanation for this discrepancy as well as for the broad range of experimentally observed current values has not yet been found .
@xmath0calculation for a molecular region containing all three molecular clusters shown in fig .
[ fig : molecularinteraction ] ( a ) .
there are three distinct steps within the voltage range of 1v . ]
we have shown , that the peak structure of the transmission function is robust against changes in the number of adsorbed molecules , as long as the distance between molecules is considerably large ( @xmath143 ) . and also does the exact shape of the top metallic lead not influence the qualitative structure of the transmission function . only if the distance between molecules becomes so small , that inter - molecular interactions are no longer negligible ( which is below 6  in our case ) ,
does the transmission function undergo a qualitative change .
namely an additional peak structure occurs .
how does this finding compare to the experimental data ? as we have pointed out in section [ sec : experiment ] , only in devices using molecules with two cyanide end - groups a more or less random peak structure was observed in the @xmath1 characteristic@xcite . in other devices , molecules with at least one thiol end - group are typically used .
these show significantly less peak structure .
we therefore give the following interpretation : the thiol end - group is known to bind strongly to au atoms .
it is therefore likely , that thiol - based mono - layers stably adsorb to gold leads .
resulting periodic structures are then robust against distortions .
the conductance of such structures is proportional to the corresponding single molecule conductance , i.e.  the number of molecules involved changes the absolute value of the current only , not the peak structure .
the random like peak structure in devices , made up of cyanide based molecules suggests , that there are a some molecular clusters present in the mono - layer .
these clusters might occur , because the binding of a cyanide end - group to au is considerably weaker compared to that of a thiol end - group , and weaker binding results in a less robust periodic structure .
we would like to thank xavier bouju , as well as udo beierlein for helpful discussions .
the current properties of each channel can be related to the corresponding eigenvalue .
we start from equation : @xmath144 let s first consider @xmath145 , i.e.  @xmath146 , i.e.  @xmath147 . for each channel with eigenvalue @xmath148 one then must have @xmath149 , i.e.   this channel does not itself carry any current .
this is consistent with our terminology of an evanescent wave . if , however , @xmath150 , then @xmath151 is purely imaginary , because @xmath152 is an anti - hermitian operator .
we can therefore define the velocity of a propagating wave to be @xmath153 .
now we consider the case of two different solutions @xmath154 and define @xmath155 . if their eigenvalues do not satisfy @xmath156 , then the current between these two solutions is zero @xmath157 .
so let s assume @xmath158 . because if @xmath159 then @xmath160 , a current can flow between an evanescent left going wave and an evanescent right going wave .
but if we restrict ourselves to solutions with finite amplitudes in a semi - infinite lead , then either the left or right going wave amplitude must be zero .
therefore evanescent waves do neither carry a current themselves nor do they exchange current with other channels , that is they do not at all contribute to the net current .
finally we are left with the case @xmath161 , with @xmath162 .
this is equivalent to @xmath163 , i.e.  the case of degenerate eigenvalues . therefore propagating waves to degenerate eigenvalues
do exchange current .
that in turn means , that the current of a superposition of two such waves does not necessarily equal the sum of the two individual currents , which is problematic as we want to express the total current as a sum of independent channels .
however , the propagating and evanescent waves were obtained by diagonalising the propagator @xmath42 .
this transformation is unique up to rotations in every degenerate eigenvalue subspace . because @xmath73 is anti - hermitian we can diagonalise these subspaces and the resulting diagonal elements
will be purely imaginary .
so the net current may be written as a summation over all the individual contributions of propagating channels , only if these subspace rotations are performed .
summarising we have shown that the transformation @xmath84 diagonalising the propagator @xmath42 ( i.e  @xmath164 ) can be chosen such that the transformation @xmath165 of the current operator is diagonal in the subspace of propagating waves with purely imaginary diagonal elements .
all the other diagonal entries are zero and the only non - zero non - diagonal elements belong to evanescent waves in opposite directions .
the part of the hamiltonian containing the molecular region and its coupling to the leads can be written as @xmath166 |\psi\rangle = 0.\ ] ] ( using this order for the coefficients it is straight forward to extend all formulas to the general case of more than two leads . )
the indices 1 and 2 indicate source and drain lead surface layers , while the index 0 is used for the molecular region .
@xmath167 are the coupling matrices from source / drain to the molecules .
we now transform into the basis of incoming and outgoing channels , i.e. we apply @xmath168 ,      \quad\text{with }      u^i=\left[\begin{array}{cc }          u^i _
>       & u^i _
<     \\          u^i_>\lambda^i _ > & u^i_<\lambda^i_<\\          \end{array}\right]\ ] ] from the right to equation : @xmath169,\ ] ] with @xmath170 the first and third column act on the surface layer of the incoming channels , the second and fourth act on outgoing ones , while the fifth column , acting on the molecular region , remains unchanged .
the scattering matrix expresses the outgoing channel amplitudes in terms of the incoming ones .
therefore we split the matrix of equation into two parts , one containing the outgoing columns , the other one containing the incoming ones as well as the molecular column : @xmath171\quad\text{and}\\      m_\text{in } & : =      \left [        \begin{array}{cc }          a^1 _ > & 0     \\          0      & a^2_>\\          b^1 _ > & b^2_>\\        \end{array }      \right ]
.    \end{aligned}\ ] ] the first matrix @xmath172 is a square matrix and by inverting it , we obtain the scattering matrix @xmath173                                    m.  j.  frisch , g.  w.  trucks , h.  b.  schlegel , g.  e.  scuseria , m.  a.  robb , j.  r.  cheeseman , v.  g.  zakrzewski , j.  a.  montgomery , jr . , r.  e.  stratmann , j.  c.  burant , s.  dapprich , j.  m.  millam , a.  d.  daniels , k.  n.  kudin , m.  c.  strain , o.  farkas , j.  tomasi , v.  barone , m.  cossi , r.  cammi , b.  mennucci , c.  pomelli , c.  adamo , s.  clifford , j.  ochterski , g.  a.  petersson , p.  y.  ayala , q.  cui , k.  morokuma , d.  k.  malick , a.  d.  rabuck , k.  raghavachari , j.  b.  foresman , j.  cioslowski , j.  v.  ortiz , a.  g.  baboul , b.  b.  stefanov , g.  liu , a.  liashenko , p.  piskorz , i.  komaromi , r.  gomperts , r  .l .
martin , d  .j .
fox , t.  keith , m  .a .
al - laham , c.  y.  peng , a.  nanayakkara , c.  gonzalez , m.  challacombe , p.  m  .w .
gill , b.  g.  johnson , w.  chen , m.  w.  wong , j.  l.  andres , m.  head - gordon , e.  s.  replogle and j.  a.  pople , gaussian 98 ( revision a.7 ) , gaussian , inc .
, pittsburgh pa , 1998 .