with the successful launch in 1999 of the _ chandra x - ray observatory _
@xcite came the opportunity to explore the x - ray universe with greater angular resolution and sensitivity in the  energy range than ever before .
one of the primary goals of the mission is to perform surveys of the extragalactic sky at these higher energies .
the longest _ chandra _ pointings thus far are the _ chandra _ deep field - north  @xcite and -south  @xcite , with 2  ms and 1  ms exposure times , respectively .
these fields have explored faint sources , with many objects reaching below the luminosity range of active galactic nuclei ( agn ) to include  normal " galaxies .
though the deep fields are pushing to lower fluxes than previously attainable , they cover only a small area ( @xmath4 deg@xmath1 each ) , and do not have good statistics for the flux range in which the dominant contribution to the  x - ray background arises . to get such statistics , a larger area of the sky
must be covered by combining multiple _
_ pointings .
the _ chandra _ archive provides the opportunity to explore multiple , mid - depth ( @xmath5 ks ) exposures . in we presented the serendipitous extragalactic x - ray source identification ( sexsi ) program design , our x - ray analysis techniques and a catalog of over one thousand  sources , as well as some initial x - ray results and comparisons with previous work . with data from 27 archival _
fields , sexsi is designed to focus on a large - area ( @xmath6 deg@xmath1 ) sample of x - ray sources in the intermediate flux range ( @xmath7 ) .
our goal is to determine spectroscopically the redshift distributions and nature of hard ( ) x - ray source populations from which the bulk of the x - ray background arises .
the campaign of single - color photometry presented here attempts a minimum depth for every field of @xmath8 , consistent with identifying counterparts for which classifiable spectra can be obtained with a 10-m telescope in a @xmath9 hour integration .
this strategy maximizes the number of hard sources identified spectroscopically in our fields for our fixed observing campaign .
other programs are also aiming to understand the x - ray source population at medium depth with degree - scale _ chandra _ and _ xmm - newton _ surveys .
examples of such surveys include hellas2xmm ( e.g. , * ? ? ?
* ; * ? ? ?
* ) and the _ chandra _ multiwavelength project ( champ ; e.g. , * ? ? ?
@xcite presents a comprehensive list of the many ongoing x - ray background surveys .
the sexsi program is unique among the _ chandra _ serendipitous surveys in concentrating on obtaining a uniform spectroscopic survey restricted to the hard source population .
the champ , in contrast , is a comprehensive program of imaging and spectroscopy to follow - up both soft and hard sources in a large number of publicly - accessible extragalactic _
the sexsi spectroscopic catalog will therefore contain a larger fraction of obscured agn than champ .
eight of the fields are common to both surveys , and for one of these champ has published some spectroscopic followup @xcite . even for the common fields , sexsi plans deeper spectroscopy and should identify a greater number of hard x - ray sources . between the large number of distinct fields and the different focus for ground - based followup , the sexsi survey provides an essential complement to current surveys of the sources responsible for the x - ray background .    in this paper
, we present optical @xmath3-band imaging of the sexsi _ chandra _ fields , as well as the methodology used to derive x - ray source identifications and a catalog of @xmath3-band counterpart magnitudes .
in addition , we explore relationships between the x - ray flux , hardness ratio , and @xmath3-band magnitude . a third paper in this series ( eckart et al .
, in preparation ; paper iii ) will include our optical spectroscopic results , with redshifts and classifications for @xmath10 of these objects , and a discussion of the luminosity distribution , redshift distribution , and composition of the sample .
we have imaged the sexsi fields using the palomar 60-inch ( p60 ) and 200-inch ( p200 ) telescopes , the mdm 2.4-meter ( mdm-2.4 ) and 1.3-meter ( mdm-1.3 ) telescopes , and the keck  i telescope .
much of the early optical imaging for this program was done with small field - of - view , single - ccd array cameras .
for example , the p60 ccd13 camera has a 12.6@xmath1112.6  field of view .
this requires several pointings to map a typical @xmath12  _ chandra _ field , making multi - color imaging very time consuming . since spring 2000 , both the p200 and mdm-2.4 have had large - format ( @xmath13 ) ccd cameras available , allowing us to cover nearly all of a _
field in a single pointing .
the p200 large format camera ( lfc ; * ? ? ? * ) uses the sloan digital sky survey ( sdss ) filter set ( @xmath14 , @xmath15 , @xmath16 , and @xmath17 ; * ? ? ?
* ) , while , for our mdm-2.4 8k camera @xcite imaging , we used johnson - cousins filters ( @xmath18 , @xmath19 , @xmath3 , and @xmath20 ) .
in each filter band , our goal has been to obtain images with a limiting vega magnitude of @xmath21 , the practical limit for obtaining counterpart spectra at the keck telescope in reasonable integration times .
because the large - format cameras were not available for the first 1.5 years of the sexsi project , some fields are imaged only in @xmath3 band ( due to observing time constraints ) , while the fields observed later in the program generally have multicolor data .
for the purposes of this paper , we use only the @xmath3-band data .
table  [ tbl : cameras ] summarizes the characteristics of the telescopes and imaging instruments used in this program .
we reduced the optical images using standard techniques , relying upon iraf tasks for the preliminary steps .
median - combined bias frames , taken the same night as the science data , were subtracted from the images . for flatfielding , we used either dome flats or sky flats generated by combining all the images from a night , including twilight flats , with min / max rejection to remove stars .
we used standard techniques to mask out cosmic rays and bad pixels prior to combining data .    for astrometric calibration of most of the data
, we used the dophot software package to locate all nonsaturated , point - like objects in individual exposures . comparing with the usno - a2.0 catalog @xcite
, we then used the object - matching program starmatch ( by doug reynolds ) to both align the images before stacking and to obtain an astrometric solution for the final , combined image .
the standard deviations of the astrometric solutions are generally @xmath22 in each axis .
we reduced images using the same methodology for both the lfc data from the p200 telescope and the 8k data from the mdm-2.4 telescope , although slightly different suites of software were required .
we relied upon the mscred package within iraf , designed for reduction of large - format , multi - array optical data from the kitt peak / cerra tololo mosaic cameras , following the informative reduction pipeline outlined by b.  jannuzi . for the lfc images this necessitated first changing the data into the requisite format using the lfcred package created by m.
huntmph / lfcred . ] .
we relied on a variety of methods to provide photometric calibration of the data , depending upon the conditions under which data for a sexsi field were taken . throughout ,
unless otherwise noted , we refer our magnitudes to the vega - based system .
if conditions were photometric , we determined magnitude zeropoints from @xcite fields observed the same night , at an airmass close to that of the field in question . in cases where a non - photometric image overlaps a photometric image from another night , we determined the zeropoint magnitude by matching the photometry of overlapping regions .    in cases where we obtained no photometric images of a field , we relied on shallow sky surveys to provide magnitude zeropoints .
where available , we used @xmath15 photometry from the sdss first data release @xcite , converting from ab magnitudes to vega magnitudes with @xmath23 alternatively , we used the usno - b1.0 catalog @xcite to establish the zeropoint .
typical photometric calibration errors for this technique were 0.2 mag .
for the purposes of this paper , we make the assumption that the sdss @xmath15 and johnson - cousins @xmath3 filters are similar . convolving the respective filter transmission curves with the @xcite composite sdss quasar spectrum
redshifted from @xmath24 to @xmath25 , we find that this assumption provides a @xmath26 mag systematic uncertainty in our photometry .
performing the same exercise on the @xcite composite galaxy templates , we find systematic offsets @xmath27 mag for @xmath24 , with larger systematic offsets at higher redshifts .
for example , the s0 composite template at @xmath28 has a 0.39  mag difference between the @xmath15 and @xmath3 vega magnitudes due to the 4000    break redshifting out of the @xmath15 filter .
though imprecise , our assumption that the @xmath15 and @xmath3 filters are similar is adequate to establish the faint source limit of each image , and to plan spectroscopic observing runs .      for extraction of sources from our combined , calibrated optical images
, we used the sextractor code of @xcite to generate catalogs for both the full optical images and for the smaller , central portions of each image .
because they avoid the ccd edges where the increased noise leads to spurious sources , these latter catalogs are used to derive photometric depths .
quoted magnitudes refer to the sextractor magbest , which usually reports photometry within elliptical apertures @xcite , but reverts to a corrected isophotal magnitude if a source is strongly contaminated by nearby sources .
to determine limiting magnitudes for each image , we made histograms of the @xmath3-band source density in 0.5 magnitude bins .
figure  [ fig : magcutoff ] shows an example histogram from sexsi field q2345 . for each sexsi field image
, we compare the measured number counts @xmath29 to that derived by @xcite from deep subaru imaging of the _ hubble _ deep field - north : @xmath30 where @xmath29 has units of number degree@xmath31 0.5 mag@xmath32 . for @xmath3-band magnitudes in the range @xmath33 , @xcite finds @xmath34 and @xmath35 . converting our @xmath3-band photometry into the ab system with equation  1
, we plot our measured number counts for each image with the published , deep - field fit ( see figure  [ fig : magcutoff ] ) . in order to derive a crude depth for our images , we compare our number counts in each 0.5  mag bin to equation  2 ( at bin center ) .
limiting magnitudes are defined by the inter - bin flux where our number counts go from greater than 80% to less than 80% of the @xcite value : this provides an @xmath36 completeness depth for each image . by the nature of our algorithm , then , these limits are good to @xmath37 mag .    table  [ tbl : images ] summarizes the @xmath3-band optical imaging we have obtained for the sexsi fields .
the average x - ray to optical offsets , @xmath38 and @xmath39 , are discussed in   [ sec : matching ] below .
we derived @xmath3-band galactic extinction for the central position of each x - ray image using the ned galactic extinction calculator . checking the variance of the ned extinction values across the _ chandra _ field of view , we find that the typical error on the extinction is @xmath40 mag .
three of our fields , however , have significantly higher extinctions than the other fields , and thus have an uncertainty of @xmath41 mag in their extinction values for individual sources .
these fields are flagged in table  [ tbl : images ] .
approximately half of the sexsi fields were observed with a single , uniform depth image ; the other half rely on imaging from multiple cameras and telescopes over several observing runs .
of the fifteen fields with a single optical image , six reach @xmath3-band depths @xmath42 , and ten reach @xmath3-band depths @xmath43 .
the shallowest field , cl  0442 + 0202 , only reaches @xmath44 , while the deepest field , cl  0848 + 4454 , reaches @xmath45 . of the twelve fields with imaging from multiple pointings ,
eight have sections reaching depths of at least @xmath46 .
a histogram of imaging depth for each source is shown in figure  [ fig : depth_hist ] .
we match x - ray and optical sources in several steps .
the procedure is slightly different depending on whether there is one optical image or multiple optical images per _ chandra _ field .    for fields with a single optical image
, we first identify the closest optical source within 4  of each _ chandra _ source position . excluding sources displaced by more than 1.5@xmath47 from the mean offset , we use these identifications to calculate a first estimate of the average offset between the _ chandra _ and usno astrometry . for the following iterations ,
the search radius used to identify optical counterparts depends upon the off - axis angle ( oaa ) of the x - ray source : we use the larger of 1.5  and psf/3 , where psf is the full - width , half - maximum of the _ chandra _ point spread function at each oaa .
again , we exclude sources displaced by more than 1.5@xmath47 from the mean as we iteratively correct the x - ray source positions .
we continue iterating until the corrections to the offsets are less than the corresponding standard deviation .
two to three iterations are typically required , with the biggest shift applied after the initial matching . as shown in table  2 ,
the typical standard deviation on this astrometric correction is @xmath48 0.3 , comparable to the accuracy of the optical astrometry .    for fields with multiple optical images
we modify this algorithm slightly .
some images cover an area with few ( or one ) x - ray sources , which , if we were to follow the procedure described above , could lead to incorrect offsets and matches . instead , for a given _
field , we first take each optical image and find the closest match within a 4  radius of all _ chandra _ source positions that fall within the image .
we record these astrometric differences . for x - ray sources falling on multiple optical images
, we use the data from the image with the deepest limiting magnitude . as before , we use this list of multi - image astrometric differences to calculate an average _ chandra _ to usno offset , again eliminating sources greater than 1.5@xmath47 from the mean offset . this astrometric correction
is then applied to all _ chandra _ positions for the sexsi field considered . a second pass at optical identifications
is then made with the oaa - dependent matching radius , using the larger of 1.5  and psf/3 .
of the 998 sexsi sources with optical coverage , 655 used a 1.5  search radius , while 343 used the psf/3 search radius .
we expect few false matches due to an optical source randomly overlapping the matching search area .
the number of false matches per field depends upon both the number of x - ray sources detected in the _ chandra _ field , and the depth of the optical image .
we tested several fields from table  [ tbl : images ] for false matches by finding the total area covered by the 1.5  radius search circles and multiplying that area by the optical source density .
we find that for medium - depth x - ray and optical images , the number of false sources detected per field is generally less than one .
for example , we predict 0.6 false matches for q2345 , 0.7 false matches for cl 0442 + 0202 , and 1.3 false matches for 1156 + 295 , where these numbers are for the 1.5  search radius sources only . for the 34% of the sources with larger
oaas and thus larger search areas we expect a larger fraction of false matches , though still fewer than 3 per field for most fields . for example , we predict 1.3 false matches from the @xmath49  match radius sources in q2345 , 1.4 false matches in cl 0442 + 0202 , and 2.4 false matches in 1156 + 295 . for several of the fields , including the three example fields just mentioned , we shifted all x - ray source positions by 1 , producing a `` fake '' x - ray catalog , and ran the matching routine again with the average x - ray - to - optical offset forced to zero .
this process , repeated with several different shifts of the x - ray source positions , gives a check on the number of false detections predicted above .
for all of the fields tested the number of false detections are consistent with the values we calculate .
our spectroscopic results are also consistent with this matching scheme and will be discussed further in paper iii .
once matching is complete , we produce first - draft optical postage stamp cutouts for our entire hard x - ray catalog ( see figure  [ fig : postage ] ) .
each stamp is centered on the astrometrically - corrected x - ray source position , and the size of the _ chandra _ psf - dependent search radius is shown by the ( centered ) solid circle . a dashed circle , located at the un - corrected _ chandra _ position , illustrates the need for this offset correction .
arrows point to identified optical counterparts , while the absence of an arrow signals a non - detection . if the magnitude found by sextractor is fainter than the limiting optical image , we annotate the photometry with an asterisk . as in paper
cl " flag is added if a source is potentially within 1  mpc of a target cluster center .
next we visually inspect postage stamps for each x - ray source , flagging sources with potentially inaccurate photometry as required .
these flags , detailed in
[ sec : catalog ] , identify sources which either ( 1 ) have multiple optical identifications within the psf - dependent search area , or ( 2 ) are , or are near , a saturated source in our optical image .
a final matching iteration is then done , excluding the newly flagged sources .
we produce new optical catalogs and postage stamp cutouts , omitting stamps for the handful of saturated sources where counterpart magnitudes are drawn from the literature and the 65 sexsi sources that either lack optical coverage in our data or have an unknown limiting magnitude due to nearby bright - source contamination ( optical flag = 6 , see   [ sec : catalog ] ) .
an example of six postage stamp cutouts is included as figure  [ fig : postage ] ; the entire catalog of postage stamps is provided in the online version of the manuscript .    for the 262 sources with @xmath50 , 160 ( 61% ) have identified counterparts , while for the 434 sources with @xmath51 , 291 ( 67% ) have identified counterparts , and for the 167 sources with @xmath52 , 124 ( 74% ) have identified counterparts .
our total sample of 947 sources with unambiguous photometry ( no contamination ) identifies 603 counterparts ( 64% ) .
in table  [ tbl : catalog ] we present the catalog of 1034 hard - band sexsi sources  the table is published in its entirety in the electronic version of the _
astrophysical journal_. columns  1  7 present x - ray source information for easy reference , while the optical photometric data are presented in columns  8  15 .
complete x - ray source information is presented in paper  i , table  4 .
the x - ray source positions in table  [ tbl : catalog ] are corrected for mean optical to x - ray offsets .
note that since the source names ( column 1 ) , identical to the source names in paper  i , are derived from the hard - band x - ray images , the refined positions of columns  2  3 will not exactly match those of column  1 ( though mean offsets are typically less than 1 ) .
column 4 lists the off - axis angle ( oaa , _
i.e. _ , the angular distance , in arcmin , of the source position from the telescope aim point ) .
the  flux ( in units of @xmath53 ) and detection snr are shown in columns  5  6 , while column  7 gives the hardness ratio , @xmath54 , where @xmath55 and @xmath56 are the counts cm@xmath31 in the  and  bands , respectively . here ,
as distinct from paper i , we record the hardness ratio derived from the net soft x - ray counts recorded at the hard - band source position when there was not a significant soft - band source detected ( in paper i these cases are reported as @xmath58 ) .
in addition , for a subset of these cases , when the soft - band counts recorded at the hard - band source position were less than twice the soft - band background counts , the @xmath59 is considered a lower - limit , flagged as such in the catalog , and set to @xmath60 , where @xmath61 soft - band background counts .
column  8 contains an optical flag code essential to interpretation of the optical data : 0 = no optical coverage , 1 = a solid optical i d , 2 = no optical counterpart ( the magnitude listed is then a lower limit ) , 3 = saturated in the sexsi optical image ( @xmath3-band magnitude taken from the guide star catalog ii  @xcite or a secondary source in vizier database ) , 4 = a solid i d but @xmath3-band magnitude affected by nearby bright source , 5 = more than one optical source in x - ray error circle ( the source with the smallest positional offset is recorded in the table ) , and 6 = lower limit ( no optical counterpart ) but area contaminated by nearby bright source so limiting magnitude is unknown .
column  9 is the @xmath3-band magnitude of the optical counterpart , with its error shown in column  10 .
column  11 gives the limiting magnitude for the image from which the optical counterpart , or a lower limit thereto , was derived .
column  12 lists @xmath62 in arcsec , where @xmath63 is the astrometrically - corrected x - ray source ra ( from column 2 ) and @xmath64 is the optical counterpart ra ; column 13 is the corresponding value for declination .
column 14 , @xmath65 ( @xmath66 ) , is the x - ray to optical position difference , also in arcsec .
following @xcite and @xcite , the logarithmic x - ray - to - optical flux ratio ( column  15 ) is given by the relation @xmath67 derived from the kron - cousins @xmath3-band filter transmission function .
column 16 presents the sextractor class_star parameter @xcite for basic star - galaxy separation , with values ranging from 0.0 for significantly extended sources to 1.0 for sources consistent with perfectly stellar psfs .
the stellarity value of the 8 bright sources with optflag = 3 ( from the literature ) is set to 1.00 .
we present the stellarity data with the caveat that the values should only be used for broad separation , for example , of sources near 0 verses 1 , but not for detailed quantitative analysis .
note that for all of the optical counterpart data columns , special attention must be paid to the optical flag . for example , if the code is a 2 ( optical counterpart is not detected ) , column  9 describes a lower limit , columns  10 , 12  14 , and 16 have no data , and the x - ray - to - optical flux ratio in column 15 is a lower limit .
we now discuss the results of the optical identifications of the hard x - ray source counterparts and how the counterpart properties relate to the x - ray properties . throughout this section ,
only sources with a solid optical counterpart or upper limit thereto are plotted and analyzed ; we exclude sources with @xmath3-band magnitudes or limits contaminated by a nearby bright source and sources with more than one optical source within the identification radius .
this provides a sample of 947 hard x - ray sources with unambiguous counterpart photometry , of which 603 are identifications and 344 have limits to the counterpart optical magnitude .
figure  [ fig : r_hflux_scatter ] presents the @xmath3-band magnitudes , or limits thereto , plotted as a function of  x - ray flux .
dashed lines show constant x - ray - to - optical flux ratios .
the bulk of the sexsi sources have @xmath68 and are fainter than @xmath69 .
shallow , wide - area x - ray surveys , such as the _ asca _ hard x - ray survey reported by @xcite , find that the majority of extremely bright x - ray sources ( @xmath70 ) are agn and have @xmath71 .
many of the sexsi sources are likely fainter and/or more distant analogs .    as has been found in other surveys ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , sexsi detects a population of sources overluminous in the x - ray for their optical magnitudes .
this population is found over the entire x - ray flux range sampled .
two explanations have been commonly discussed to explain large values of @xmath72 : ( i ) extremely high - redshift agns that might have bright x - ray fluxes but faint @xmath3-band fluxes due to absorption from the lyman transitions of hydrogen along our line of sight , or ( ii ) heavy obscuration by material in the host galaxy . the latter interpretation for the majority of sources
is supported by near - ir studies , which find a large fraction to be very red early - type galaxies at @xmath73 @xcite .
in addition , a number of _ type  ii quasars _ , which would have similar properties , have been identified in the deepest _ chandra _ surveys ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) .    by virtue of the large area surveyed ,
sexsi contains numerous sources with @xmath74 at moderate x - ray fluxes , @xmath75 .
many of these sources have @xmath76 , well within the reach of 8- to 10-meter class telescopes for spectroscopic follow - up .
these are brighter examples of the new x - ray overluminous population .
deep survey versions of these sources are typically extremely faint  as an example , @xcite report on seven `` extreme x - ray / optical '' ( exo ) sources in the 2  msec _ chandra _ deep field - north ( cdf - n ) which , despite having extremely robust detections in the _ chandra _ data ( @xmath77 counts ) , remain undetected ( @xmath78 ) in the great observatories origins deep survey imaging of the cdf - n @xcite .
@xcite  report on the spectra of 13 sources with @xmath79 from the hellas2xmm survey .
they find 8 narrow - lined sources , with @xmath80 erg s@xmath32 .
our survey contains 109 such sources with @xmath79 of which 64 have an optical limiting magnitude and @xmath81 have spectral data .
these high @xmath82 sources provide a useful catalog for future infrared surveys .
we also find a number of sources with very bright ( @xmath83 ) optical counterparts .
the bulk of these sources ( 7 out of 8) are identified as stars in the literature .
figure  [ fig : hfluxhist ] presents a histogram of the number of sources in sexsi and in the cdf - n  @xcite as a function of hard x - ray flux and split by @xmath3-band magnitude at @xmath84 .
this magnitude was chosen to separate our sources approximately in half . the small number of sources from images with @xmath85 were not included in this plot or in figures  [ fig : hrhist ] ,  [ fig : meanhr ] , and  [ fig : fxfohr ] .
the majority of the sexsi sources have @xmath86 , a flux level that lies between the _ asca _ and _ bepposax _ sensitivity limits and the _ chandra _ deep survey capability .
this important flux range corresponds to the regime in which the  relation changes slope and from which the bulk of the  x - ray background arises @xcite . compared to the deep surveys , sexsi has nearly an order of magnitude more sources at this x - ray depth , approximately half of which have @xmath87 , making them easy spectroscopic targets for 8- to 10-meter class telescopes .
figure  [ fig : hrhist ] shows the hardness ratio histogram of the sexsi sources , again split at @xmath84 .
the optically - brighter sources are peaked at a low hardness ratio , while the fainter sources have a much harder , and broader , distribution , showing that the optically - fainter portion of the sample has a higher fraction of sources with a flat x - ray spectral slope .
the large peak of the @xmath88 sources is near @xmath89 , which corresponds to a power - law photon index ( @xmath90 ) between 1.8 and 1.9 , a typical value for unobscured broad - line agn .
this suggests that the majority of the optically - brighter sources are broad - lined agn .
spectroscopic followup ( paper iii ) will be able to specifically address this hypothesis .
the harder and broader distribution of the @xmath91 sources indicates that obscuration at the source is likely involved .
@xcite  shows a similar trend in the 1 ms cdf - n data .
figure  [ fig : meanhr ] , which displays the mean hardness ratio as a function of hard - band flux separately for bright and faint counterpart magnitudes , again emphasizes that distinct populations are contributing to the overall source counts .
the mean values for this figure are calculated using source @xmath59s regardless of if they are flagged as lower limits .
setting all such flagged sources to the upper extreme , @xmath92 , shifts the points only slightly , if at all . for reference , the @xmath59 of a @xmath93 spectrum with varying intrinsic absorbing column density and redshift
is presented in the adjacent panel .
this shows that modest - redshift , heavily - obscured sources produce the highest @xmath59s . for sources fainter than @xmath84 ,
the mean hardness ratio is essentially constant at @xmath94 over two orders of magnitude in x - ray flux . for those sources with brighter magnitudes ,
however , there is a significant decline in mean hardness ratio , from @xmath95 at @xmath96 to @xmath97 at @xmath98 . in a simple - minded two - component model with median hardness ratios of @xmath99 and @xmath100 for the two source classes ,
the faint sources are dominated by the hard population with perhaps an admixture of 10% soft sources at all flux levels , while for the brighter sources , the ratio of hard to soft sources changes from 20:80 at @xmath96 to 80:20 at @xmath98 .
this is qualitatively consistent with the  relations derived in paper i , where we showed that hard sources follow a constant power - law over three orders of magnitude in hard - band flux , while softer sources show a distinct break at @xmath101 . in this picture , the soft sources are the standard , largely unobscured agn at all redshifts , while the harder component represents the lower mean - redshift , lower - luminosity populations revealed in the deep surveys ( e.g. , * ? ? ?
* ; * ? ? ?
figure  [ fig : fxfohr ] shows the hardness ratio as a function of x - ray to optical flux ratio for the sexsi sample , again with the @xmath59 of a @xmath93 spectrum with varying intrinsic absorbing column density as a function of redshift presented in the adjacent panel .
the main plot shows a trend toward larger hardness ratios as @xmath82 increases .
this is in general consistent with the notion that the increasing @xmath82 results from larger absorbing columns which attenuate the @xmath3-band flux much more severely than the hard x - ray flux ( e.g. , * ? ? ?
it is interesting , however , that the most extreme sources with @xmath79 are not all hard , but exhibit a wide range of @xmath59s , implying that obscuration may not be the sole explanation for the dim optical counterparts .
@xcite  find that for non - broad - lined sources ( e.g. , sources with optical spectra showing narrow agn emission lines or early - type galaxy absorption lines ) there is a linear correlation between @xmath72 and @xmath102 .
this relationship implies that some of the high @xmath82 sources with low @xmath59s are high luminosity , high redshift , narrow - lined agn .
we also explored the optical source densities in the vicinity of the x - ray sources .
for each optical counterpart , we counted the number of sources in the surrounding region and compared the source density to the overall field density to calculate overdensities .
we employed 20 , 30 , 40 , and 60  radius circles around each optical counterpart and counted sources with @xmath3-magnitude of @xmath103  and @xmath104  for each of the sources .
we found no significant overdensities around the x - ray sources .
we present @xmath3-band imaging of  @xmath105% of the  x - ray sources in the sexsi survey .
we describe our optical data reduction and the x - ray to optical source matching algorithms employed , and present a catalog of @xmath3-band counterpart photometry and astrometry .
while the power of the sexsi sample is fully realized only with the addition of our collection of optical counterpart spectra , the photometric identification of the sources is an important step that provides clues to the composition of the source population .
we find that by splitting the sources at @xmath84 and analyzing each group s x - ray properties , we begin to see the emergence of what appear to be physically distinct populations .
for example we find that the hardness ratio distribution of the optically brighter sources is sharply peaked near @xmath106 , typical of unobscured agn , while the fainter , @xmath91 sources have a much broader and harder spectral distribution .
the findings of this paper suggest that the analysis of optical spectral data from counterparts fainter than @xmath84 ( attainable with 10-meter - class telescopes ) , will be an essential part of exploring the x - ray background composition .
this research has made use of the nasa / ipac extragalactic database ( ned ) which is operated by the jet propulsion laboratory , california institute of technology , under contract with the national aeronautics and space administration ( nasa ) .
some of the data presented herein were obtained at the w.m .
keck observatory , which is operated as a scientific partnership among the california institute of technology , the university of california and nasa .
the observatory was made possible by the generous financial support of the w.m .
keck foundation .
the authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of mauna kea has always had within the indigenous hawaiian community .
we are most fortunate to have the opportunity to conduct observations from this mountain .
this work has made use of the usnofs image and catalog archive operated by the united states naval observatory , flagstaff station .
the guide star catalog was produced at the space telescope science institute under u.s .
government grant with data based on photographic data obtained using the oschin schmidt telescope on palomar mountain and the uk schmidt telescope .
this research has made use of nasa s astrophysics data system .
this work has been supported by nasa nag5 - 6035 ( djh ) , as well as by a small _ chandra _ archival grant .
the work of ds was carried out at the jet propulsion laboratory , california institute of technology , under a contract with nasa .
llccl 8k & mdm 2.4 m & 0.18 & @xmath107 & @xcite + ccd13 & palomar 60  & 0.37 & @xmath108 & + cosmic & palomar 200  & 0.29 & @xmath109 & @xcite + echelle & mdm 1.3 m & 0.50 & @xmath110 & + echelle & mdm 2.4 m & 0.28 & @xmath111 & + lfc & palomar 200  & 0.18 & @xmath112 & @xcite + lris & keck i & 0.22 & @xmath113 & @xcite + templeton & mdm 1.3 m & 0.50 & @xmath114 & + [ tbl : cameras ]    lcccclllcc ngc 891 & 51 & -0.9 @xmath115 0.3 &  0.6 @xmath115 0.2 & 0.17 & 02 22 40 & + 42 26 16 & 8k & 24.3 & 1.6 + awm 7 & 48 &  0.5 @xmath115 0.3 &  0.2 @xmath115 0.1 & 0.30 & 02 54 45 & + 41 40 10 & 8k & 24.2 & 1.7 + xrf 011130 & 30 &  0.4 @xmath115 0.2 & -0.5 @xmath115 0.2 & 0.26 & 03 05 28 & + 03 48 52 & lfc & 23.2 & 1.2 + ngc 1569 & 97 & -0.3 @xmath115 0.2 &  0.6 @xmath115 0.3 & 2.00 & 04 29 08 & + 64 45 42 & echelle ( 1.3 ) & 21.5 & 1.7 + & & & & & 04 29 42 & + 64 42 35 & lris & 21.5 & 0.9 + & & & & & 04 29 53 & + 64 39 47 & ccd13 & 21.0 & 1.4 + & & & & & 04 30 36 & + 64 51 06 & lris & 23.0 & 1.2 + & & & & & 04 31 15 & + 64 51 06 & lris & 22.5 & 0.9 + 3c 123 & 47 & -0.1 @xmath115 0.1 &  0.0 @xmath115 0.0 & 2.61 & 04 36 52 & + 29 40 08 & 8k & 21.9 & 1.4 + cl 0442 + 0202 & 44 & -0.2 @xmath115 0.2 & -1.2 @xmath115 0.3 & 0.41 & 04 42 17 & + 02 03 35 & lfc & 21.1 & 1.2 + cl 0848 + 4454 & 186 & -0.3 @xmath115 0.1 & -1.2 @xmath115 0.1 & 0.08 & 04 48 49 & + 44 54 09 & lfc & 24.4 & 1.1 + rx j0910 & 171 & -0.3 @xmath115 0.1 & -2.3 @xmath115 0.2 & 0.05 & 09 10 33 & + 54 19 38 & 8k & 24.0 & 1.6 + 1156 + 295 & 49 & -0.5 @xmath115 0.2 & -1.1 @xmath115 0.3 & 0.05 & 11 59 32 & + 29 16 19 & 8k & 24.5 & 1.5 + ngc 4244 & 49 & -0.2 @xmath115 0.2 &  0.0 @xmath115 0.1 & 0.06 & 12 16 46 & + 37 51 46 & ccd13 & 22.9 & 1.8 + & & & & & 12 17 34 & + 37 46 54 & echelle ( 1.3 ) & 23.4 & 1.8 + & & & & & 12 17 54 & + 37 55 15 & ccd13 & 23.4 & 1.5 + ngc 4631 & 59 &  0.4 @xmath115 0.3 & -0.3 @xmath115 0.2 & 0.05 & 12 41 29 & + 32 36 46 & ccd13 & 23.0 & 1.6 + & & & & & 12 41 39 & + 32 30 44 & ccd13 & 23.0 & 2.5 + & & & & & 12 42 03 & + 32 40 53 & echelle ( 2.4 ) & 23.0 & 1.2 + & & & & & 12 42 21 & + 32 37 09 & ccd13 & 23.0 & 1.7 + & & & & & 12 42 59 & + 32 32 49 & templeton & 23.0 & 1.4 + hcg 62 & 49 &  0.1 @xmath115 0.1 & -0.1 @xmath115 0.1 & 0.13 & 12 53 04 & @xmath11609 14 20 & echelle ( 2.4 ) & 22.4 & 1.2 + & & & & & 12 53 06 & @xmath11609 05 48 & echelle ( 2.4 ) & 23.4 & 1.0 + & & & & & 12 53 04 & @xmath11609 14 20 & ccd13 & 22.4 & 1.6 + rx j1317 & 111 & -1.6 @xmath115 0.2 & -1.6 @xmath115 0.2 & 0.03 & 13 17 13 & + 29 11 32 & lfc & 23.5 & 1.5 + bd 1338 & 38 & -1.3 @xmath115 0.2 & -1.1 @xmath115 0.1 & 0.04 & 13 37 44 & + 29 25 24 & ccd13 & 22.5 & 1.3 + & & & & & 13 37 53 & + 29 29 44 & ccd13 & 22.5 & 1.7 + & & & & & 13 38 15 & + 29 36 51 & ccd13 & 22.5 & 2.0 + & & & & & 13 38 26 & + 29 24 53 & ccd13 & 22.5 & 1.5 + rx j1350 & 58 &  0.7 @xmath115 0.3 & -2.5 @xmath115 0.5 & 0.03 & 13 50 39 & + 60 04 11 & lfc & 22.5 & 0.9 + 3c 295 & 23 &  0.2 @xmath115 0.2 & -0.3 @xmath115 0.3 & 0.05 & 14 11 01 & + 52 21 18 & ccd13 & 22.5 & 2.2 + & & & & & 14 11 24 & + 52 13 36 & ccd13 & 21.5 & 2.8 + & & & & & 14 11 46 & + 52 05 42 & ccd13 & 22.0 & 2.3 + grb 010222 & 18 & -0.6 @xmath115 0.2 &  0.1 @xmath115 0.1 & 0.06 & 14 51 59 & + 43 08 54 & ccd13 & 21.9 & 2.1 + & & & & & 14 52 13 & + 43 01 06 & ccd13 & 22.9 & 1.4 + & & & & & 14 52 43 & + 42 55 12 & ccd13 & 22.4 & 1.6 + & & & & & 14 52 59 & + 43 06 07 & ccd13 & 22.9 & 1.3 + qso 1508 & 89 & -0.6 @xmath115 0.2 & -0.2 @xmath115 0.1 & 0.03 & 15 09 50 & + 57 04 16 & lfc & 24.0 & 1.2 + mkw 3s & 57 &
-0.8 @xmath115 0.1 &  0.1 @xmath115 0.1 & 0.09 & 15 21 30 & + 07 47 43 & ccd13 & 23.4 & 1.2 + & & & & & 15 21 37 & + 07 39 08 & ccd13 & 22.9 & 1.1 + & & & & & 15 22 12 & + 07 48 36 & ccd13 & 22.9 & 1.3 + & & & & & 15 22 13 & + 07 47 21 & ccd13 & 22.9 & 1.6 + ms 1621 & 30 & -1.2 @xmath115 0.1 &  0.3 @xmath115 0.2 & 0.09 & 16 23 06 & + 26 36 33 & ccd13 & 22.4 & 2.3 + & & & & & 16 23 13 & + 26 30 22 & cosmic & 24.4 & 1.4 + & & & & & 16 23 26 & + 26 38 28 & cosmic & 23.9 & 1.6 + & & & & & 16 23 42 & + 26 43 49 & ccd13 & 22.9 & 2.3 + & & & & & 16 24 03 & + 26 35 32 & cosmic & 23.9 & 1.2 + grb 000926 & 32 & -0.7 @xmath115 0.2 & -0.6 @xmath115 0.3 & 0.08 & 17 03 19 & + 51 47 58 & 8k & 22.9 & 1.2 + rx j1716 & 52 & -0.5 @xmath115 0.1 &  1.2 @xmath115 0.1 & 0.09 & 17 15 36 & + 67 19 45 & ccd13 & 22.4 & 1.7 + & & & & & 17 15 55 & + 67 05 59 & cosmic & 24.4 & 1.3 + & & & & & 17 16 39 & + 67 13 27 & cosmic & 23.9 & 1.4 + & & & & & 17 17 11 & + 67 01 40 & cosmic & 24.4 & 1.2 + & & & & & 17 17 42 & + 67 12 36 & ccd13 & 22.9 & 1.3 + & & & & & 17 17 56 & + 67 09 08 & cosmic & 23.9 & 1.6 + & & & & & 17 18 04 & + 67 12 01 & ccd13 & 22.4 & 2.2 + ngc 6543 & 46 & -0.5 @xmath115 0.2 & -0.2 @xmath115 0.1 & 0.12 & 17 57 44 & + 66 44 40 & ccd13 & 22.9 & 1.3 + & & & & & 17 58 31 & + 66 34 02 & ccd13 & 22.9 & 1.4 + & & & & & 18 00 25 & + 66 30 52 & ccd13 & 22.9 & 1.2 + xrf 011030 & 47 & -0.2 @xmath115 0.2 & -0.4 @xmath115 0.2 & 1.09 & 20 44 02 & + 77 21 06 & lfc & 21.4 & 2.2 + ms 2053 & 44 & -1.1 @xmath115 0.1 &  0.4 @xmath115 0.1 & 0.22 & 20 55 48 & @xmath11604 31 10 & ccd13 & 22.8 & 1.7 + & & & & & 20 55 59 & @xmath11604 35 47 & cosmic & 23.8 & 1.4 + & & & & & 20 56 09 & @xmath11604 44 01 & cosmic & 23.8 & 1.4 + & & & & & 20 56 30 & @xmath11604 38 17 & lris & 23.8 & 1.0 + & & & & & 20 56 33 & @xmath11604 30 53 & ccd13 & 22.3 & 2.1 + & & & & & 20 56 42 & @xmath11604 41 28 & cosmic & 23.8 & 1.2 + rx j2247 & 49 & -0.8 @xmath115 0.1 &  1.0 @xmath115 0.1 & 0.15 & 22 47 21 & + 03 39 25 & 8k & 23.4 & 1.7 + q2345 & 74 &  0.0 @xmath115 0.1 & -0.3 @xmath115 0.2 & 0.07 & 23 48 25 & + 01 01 12 & 8k & 23.4 & 1.5 + [ tbl : images ]    lccrrrrclccrrcrrc j022142.6 + 422654 & 02 21 42.75 & + 42 26 53.5 & 9.49 & 23.10 & 2.83 & -0.33 & 1 & 18.73 & 0.16 & 24.3 & 1.0 & 1.8 & 2.1 & -0.64 & 0.66 + j022143.6 + 421631 & 02 21 43.72 & + 42 16 31.0 & 8.33 & 74.80 & 8.81 & -0.28 & 6 &   &  &  &  &  &  &  &
 + j022151.6 + 422319 & 02 21 51.76 & + 42 23 18.7 & 6.17 & 6.55 & 2.06 & -0.16 & 2 & 24.3 &  & 24.3 &  &  &  & 1.04 &  + j022204.9 + 422338 & 02 22 05.08 & + 42 23 37.7 & 4.24 & 8.73 & 2.74 & 0.84 & 1 & 20.37 & 0.16 & 24.3 & -0.2 & 0.5 & 0.5 & -0.41 & 0.02 + j022205.1 + 422213 & 02 22 05.21 & + 42 22 12.7 & 3.45 & 6.82 & 2.38 & 0.07 & 2 & 24.3 &  & 24.3 &  &  &  & 1.07 &
 + j022207.1 + 422918 & 02 22 07.19 & + 42 29 18.2 & 8.93 & 18.90 & 3.84 & -0.33 & 1 & 21.56 & 0.16 & 24.3 & 2.3 & -0.8 & 2.4 & 0.40 & 0.96 + j022210.0 + 422956 & 02 22 10.08 & + 42 29 55.7 & 9.38 & 17.30 & 3.52 & -0.62 & 1 & 18.88 & 0.16 & 24.3 & 0.2 & -0.6 & 0.6 & -0.71 & 0.03 + j022210.8 + 422016 & 02 22 10.93 & + 42 20 16.1 & 2.17 & 3.63 & 1.53 & -0.46 & 1 & 21.57 & 0.17 & 24.3 & 0.6 & -0.4 & 0.8 & -0.31 & 0.95 + j022211.7 + 421910 & 02 22 11.79 & + 42 19 10.1 & 2.55 & 6.87 & 2.36 & 0.14 & 6 &   &  &  &  &  &  &  &  + j022215.0 + 422341 & 02 22 15.12 & + 42 23 41.0 & 3.15 & 15.40 & 3.88 & -0.41 & 1 & 22.52 & 0.18 & 24.3 & 0.1 & 0.0 & 0.1 & 0.70 & 0.64 + j022215.1 + 422045 & 02 22 15.19 & + 42 20 44.5 & 1.32 & 42.90 & 7.49 & -0.42 & 1 & 17.33 & 0.16 & 24.3 & -0.3 & 0.1 & 0.3 & -0.94 & 0.03 + j022215.5 + 421842 & 02 22 15.63 & + 42 18 41.7 & 2.46 & 6.79 & 2.34 & 0.28 & 1 & 20.33 & 0.16 & 24.3 & -0.4 & 0.1 & 0.4 & -0.54 & 0.03 + j022219.3 + 422052 & 02 22 19.40 & + 42 20 51.6 & 0.54 & 6.38 & 2.31 & 0.37 & 1 & 23.82 & 0.22 & 24.3 & -0.9 & 0.6 & 1.1 & 0.83 & 0.38 + j022224.3 + 422138 & 02 22 24.45 & + 42 21 38.4 & 0.91 & 64.60 & 8.92 & -0.44 & 1 & 17.41 & 0.16 & 24.3 & 0.0 & 0.0 & 0.0 & -0.73 & 0.98 + j022225.2 + 422451 & 02 22 25.33 & + 42 24 50.9 & 4.07 & 128.00 & 12.99 & -0.34 & 1 & 19.31 & 0.16 & 24.3 & 0.1 & 0.2 & 0.2 & 0.33 & 0.98 + j022226.5 + 422154 & 02 22 26.63 & + 42 21 54.4 & 1.35 & 6.48 & 2.32 & -0.35 & 2 & 24.3 &  & 24.3 &  &  &  & 1.04 &
 + j022232.5 + 423015 & 02 22 32.61 & + 42 30 14.6 & 9.61 & 53.60 & 6.50 & 0.12 & 1 & 20.88 & 0.16 & 24.3 & 0.4 & 0.4 & 0.5 & 0.58 & 0.77 + j022236.3 + 421730 & 02 22 36.45 & + 42 17 30.2 & 4.23 & 12.00 & 3.30 & 0.01 & 1 & 22.38 & 0.18 & 24.3 & 0.3 & -0.9 & 1.0 & 0.53 & 0.42 + j022236.7 + 422858 & 02 22 36.88 & + 42 28 57.7 & 8.57 & 12.90 & 3.00 & -0.68 & 1 & 22.05 & 0.17 & 24.3 & -1.5 & 2.1 & 2.6 & 0.43 & 0.81 + j022259.1 + 422434 & 02 22 59.18 & + 42 24 33.6 & 7.77 & 32.50 & 5.43 & 0.49 & 2 & 24.3 &  & 24.3 &  &  &  & 1.74 &  + j022334.0 + 422212 & 02 23 34.13 & + 42 22 11.6 & 13.34 & 37.50 & 4.47 & 0.15 & 2 & 24.3 &  & 24.3 &  &  &  & 1.81 &  + [ tbl : catalog ]


the study of the electronic structure of materials usually begins with noninteracting electrons due to the vast number of particles involved .
many - body methods then provide a hierarchy of corrections to account for the coulomb interaction to various order .
many - body approaches such as gw@xcite and the bethe - salpeter equation ( bse)@xcite are frequently and successfully employed in the calculation of the electronic structure and excitations of materials . though accurate and physically sound , these many - body methods can become cumbersome and impractical for large systems due to the steep scaling of the numerical cost versus the system size .    alternatively ,
density - functional theory ( dft ) and time - dependent density - functional theory ( tddft)@xcite are popular methods for calculating electronic ground states and excitations , respectively , and are widely used in chemistry , physics , materials science , and other areas .
density - functional methods solve the many - body problem by constructing a noninteracting system which reproduces the electronic density of the interacting , physical system .
the favorable balance between accuracy and efficiency makes the resulting dft and tddft schemes unrivaled for large but finite system sizes.@xcite    considerable effort has been spent to replicate this success of tddft for periodic solids.@xcite generally speaking , tddft works very well for simple metallic systems , where the excitation spectrum is dominated by collective plasmon modes .
the reason is that common local and semilocal exchange - correlation ( xc ) functionals are based on the homogeneous electron liquid as reference system , which is an ideal starting point to describe electrons in metals .
the situation is more complicated in insulators and semiconductors .
the first problem that comes to mind is that of the band gap , which is typically strongly underestimated by most popular xc functionals of dft . in principle
, tddft provides a mechanism to obtain the correct band gap,@xcite but this puts very strong demands on the xc kernel @xmath0 ( it has to simulate a discontinuity , which requires a strong frequency dependence ) .
the second difficulty are excitonic effects .
it is a well - known fact that standard local and semilocal xc functionals do not produce any excitonic binding;@xcite again , the proper choice of @xmath0 is crucial .
there are many examples in the literature of successful tddft calculations of excitonic effects , using exact exchange , @xcite an effective xc kernel engineered from the bse,@xcite a meta - gga kernel,@xcite and a recent ` bootstrap ' xc kernel.@xcite these kernels all have in common that they have a long spatial range ; however , it has also been shown that certain excitonic features can be equally well reproduced by simple short - range kernels .
@xcite this calls for further explanation .    due to the complexity of real solids ,
the question of the general requirements for excitonic binding in tddft has been difficult to analyze . as a first step towards a simplified tddft approach for excitons ,
a two - band model was recently developed , which was used to test the performance of simple xc kernels for calculating excitonic binding energies in several iii - v and ii - vi semiconductors.@xcite in this paper we will push this reductionist approach further and propose a minimal tddft model for excitons .
our model is one - dimensional ( 1d ) and uses two simple kronig - penney - type bands as input .
we show that the minimal model reproduces and reveals many aspects of excitons .
the model is accessible and relatively easy to implement , and it can be used to identify important aspects of the xc functional for excitonic effects .
it clearly shows that excitons are collective excitations of the many - body system : the appropriate phase - coherent mixing of the single - particle excitations is accomplished via a coupling matrix featuring @xmath0 .
the properties of this coupling matrix are analyzed and compared with its bse counterpart .    in textbooks
, excitons are usually introduced through the two - body wannier equation , which describes an electron and a hole which interact via a screened coulomb potential .
while this arguably constitutes the simplest model for excitons , it is based on several drastic assumptions which are not fulfilled in general.@xcite we discuss how and under what circumstances a wannier - like equation emerges from our minimal tddft model .
the paper is structured as follows .
we give an introduction to wannier excitons in sec .
[ sect : bkgnd ] and to tddft for solids in sec .
[ sect : tddft ] , followed by a description of the minimal model in sec .
[ sect : model ] .
section [ sect : results ] then presents results for the minimal exciton model comparing tddft with the bse , and discusses various implications . in sec .
[ sect : wannier ] we show how the wannier equation emerges from tddft .
conclusions are given in in sec .
[ sect : app : details ] .
details of the bse method are provided in appendix [ sect : app : bse ] .
atomic units ( @xmath1 ) are used throughout unless otherwise stated , and we will only consider spin - unpolarized systems .
the electronic structure of crystalline solids is described by the bloch theory , where the electrons move in a periodic effective single - particle potential which reflects the crystal symmetry . as a result ,
the electronic states form energy bands . in insulators and semiconductors , electronic excitations take place between the occupied ( valence ) and unoccupied ( conduction ) bands .
these interband transitions can be described within a simple independent - particle approach based on fermi s golden rule ; one thus obtains a reasonable qualitative account of the optical properties in these materials.@xcite        however , experiments reveal that there are important modifications to this picture , as illustrated schematically in fig . [
fig : bkgnd : illustrationexcitons ] . above the band gap , the spectrum appears strongly enhanced , and below the band gap one may find discrete absorption peaks known as bound excitons .
the origin of these modifications are coulomb interactions : the simple picture of independent single - particle excitations is replaced by a more complex scenario where these excitations are dynamically coupled .
excitonic effects are ubiquitous in nature , and occur in 3d , 2d and 1d systems alike.@xcite the details of the excitonic modifications to the noninteracting spectrum depend on the dimensionality of the system.@xcite        the standard textbook explanation of excitons is based on the simple picture of an electron - hole pair held together by coulomb interactions , see fig .
[ fig : bkgnd : wannierexciton ] .
one thus arrives at a two - body problem similar to the positronium atom , with a center - of - mass momentum @xmath2 and relative motion described by the wannier equation:@xcite
@xmath3\psi_\nu({{\bf r}})=e_\nu\psi_\nu({{\bf r } } ) , \label{eqn : bkgnd : wannier}\ ] ] where @xmath4 is the reduced mass , defined as @xmath5 .
@xmath6 and @xmath7 are the effective masses of the conduction band electrons and the valence band holes .
@xmath8 is the coulomb interaction between the electron and the hole , divided by the static dielectric constant of the system . @xmath9 and @xmath10 are the excitonic wave function and binding energy , respectively .    in 1d systems ,
the coulomb interaction is ill - defined and requires , in general , some parametrized form;@xcite we will use the following soft - coulomb interaction : @xmath11 where @xmath12 and @xmath13 are parameters . in the following
, we will use @xmath14 throughout .
the wannier equation ( [ eqn : bkgnd : wannier ] ) has the form of a hydrogenic schrdinger equation , so it possesses a rydberg series ( even for 1d cases with soft - coulomb interaction ) with infinitely many eigenvalues below the band gap , and a continuum above the gap.@xcite the excitons in the wannier model for 3d and 2d cases both enhance the optical absorption near the band gap , while the presence of excitons in 1d systems suppresses the optical absorption just above the band gap .
sham and rice@xcite showed how the wannier equation can be derived from first principles starting from the bse , under the assumption that the bohr radius of the exciton is much bigger than the lattice constant .
the resulting wannier picture of excitons appears clear and intuitive , but this simplicity is somewhat deceptive . in reality , excitons are a dynamical many - body phenomenon and require a subtle coordination and cooperation of many single - particle transitions between two bands . in the following , we will develop a model based on tddft which will illustrate the true physical nature of excitons , but which will remain sufficiently transparent to allow a simple interpretation of the collective many - body effects that are responsible for the excitonic binding .
tddft@xcite is an in principle exact approach for electron dynamics , based on the uniqueness of the mapping between the time - dependent electronic density and the external potential.@xcite the key equation of tddft is the time - dependent kohn - sham ( tdks ) equation : @xmath15\phi_i({{\bf r}},t ) , \label{eqn : bkgnd : tddft}\end{aligned}\ ] ] where @xmath16 are the tdks orbitals of the noninteracting kohn - sham system which reproduces the density of the real interacting system . @xmath17 and @xmath18 are the external potential of the physical system and the hartree potential , respectively , and the xc potential @xmath19 is the only piece that needs to be approximated in practice .    the excitation spectrum of a system can be calculated via time propagation of eq .
( [ eqn : bkgnd : tddft ] ) following a suitably chosen initial perturbation .
alternatively , one can obtain excitation energies and optical spectra directly from linear - response tddft , using the so - called casida equation:@xcite @xmath20 equation ( [ eqn : bkgnd : casida ] ) is a generalized eigenvalue equation .
one obtains the optical transition frequencies from the eigenvalues @xmath21 , and the corresponding eigenvectors tell us how the kohn - sham single - particle transitions are mixed to form the transitions of the interacting system .
the matrices @xmath22 and @xmath23 in eq .
( [ eqn : bkgnd : casida ] ) are defined as @xmath24 where @xmath25,@xmath26,@xmath27,@xmath28 are labels for ground - state kohn - sham orbitals , and the @xmath29 s are the associated kohn - sham orbital energies . @xmath30 in eq .
( [ eqn : bkgnd : ab ] ) is defined as @xmath31 where the @xmath32 s are the ground - state kohn - sham orbitals , and @xmath33 is the hartree - exchange - correlation ( hxc ) kernel , defined as a fourier transform of @xmath34 the xc kernel @xmath35 has to be approximated in practice .
@xmath36 and @xmath37 make up the eigenvector in eq .
( [ eqn : bkgnd : casida ] ) , and they describe excitations and de - excitations , respectively . a commonly used approximation to eq .
( [ eqn : bkgnd : casida ] ) , known as tamm - dancoff approximation ( tda),@xcite is to set @xmath38 , so that the casida equation reduces to @xmath39\rho_{mn}(\omega)=\omega\rho_{ij}(\omega ) .
\label{eqn : bkgnd : tddftworking}\ ] ] this decouples excitations and de - excitations , and the computational cost is reduced .
there are situations , for instance for molecular excitations of open - shell systems,@xcite in which the tda is preferred in practice over the full casida equation ( [ eqn : bkgnd : casida ] ) .
we find that the tda can also be advantageous for excitons ( see sec .
[ sect : results ] ) .    by considering only a single kohn - sham transition in eq .
( [ eqn : bkgnd : casida ] ) one arrives at the small - matrix approximation ( sma):@xcite @xmath40 where @xmath41 is the kohn - sham transition frequency to be corrected .
one can further simplify this by making the tda , which yields the single - pole approximation ( spa ) : @xmath42 the sma and spa are valid when the considered excitation is far away from other transitions in the system . though they are not usually accurate enough for real calculations , their simplicity makes them very useful for theoretical analysis and development .    in periodic solids ,
the kohn - sham orbitals are labeled with the band index @xmath25 and the wavevector @xmath43 and have the bloch form : @xmath44 it is in principle possible to adapt the casida equation ( [ eqn : bkgnd : casida ] ) for the case of orbitals of the form ( [ eqn : bkgnd : blochform]),@xcite and use this to calculate the excitation spectrum .
however , to obtain the optical spectrum in solids it is more convenient to calculate the macroscopic dielectric function:@xcite @xmath45 where the @xmath46 s are reciprocal lattice vectors , and @xmath47 is the linear response function .
@xmath48 can be expressed as@xcite @xmath49 where @xmath50 labels the solutions of eq .
( [ eqn : bkgnd : tddftworking ] ) for an extended system .    beyond linear response ,
few real - time tddft calculations exist for periodic solids.@xcite instead of directly solving the tdks equation ( [ eqn : bkgnd : tddft ] ) , the tdks orbitals can be expanded in terms of ground - state kohn - sham bloch functions as @xmath51 the time - dependent density matrix is defined as @xmath52 the equation of motion of the density matrix is then @xmath53 , \label{eqn : bkgnd : tddmeom}\ ] ] with the tdks hamiltonian matrix @xmath54 defined by @xmath55 where @xmath56 is the volume of the unit cell , and @xmath57 is the tdks hamiltonian of eq . ( [ eqn : bkgnd : tddft ] ) .
this density - matrix approach has been used to derive the tddft version of the semiconductor bloch equations.@xcite in this formalism we only consider vertical transitions , where the bloch wavevector @xmath2 does not change during the dynamics .
nonvertical excitations are not considered , since they involve indirect ( i.e. , phonon - assisted ) processes which we ignore here .
solids are formally described by the many - body schrdinger equation . since exact solutions are not possible , it is instructive to resort to model systems to eliminate undesired details of the many - body system and provide clear illustrations of the specific features one is interested in .
the wannier equation for excitons presented in sec . [ sect : bkgnd ] is such a model .    while intuitive , the wannier model assumes the electron - hole interaction as given , so the excitonic effects are already built in by default
however , this does not explain under what conditions one expects to see the formation of excitons , and the many - body nature of excitonic effects remains hidden .
we therefore propose a minimal model for excitons which lowers the abstraction level of the wannier model , and where excitonic effects show up without any _ ad hoc _ assumptions .    for excitations near the band gap ,
a reasonable approximation is to use a two - band model , i.e. , only to consider the highest valence band @xmath58 and the lowest conduction band @xmath59 .
this means that we only need to consider those elements of the time - dependent density matrix @xmath60 [ eq .
( [ tddm ] ) ] for which @xmath61 and @xmath62 ; the latter index will be dropped in the following .    for the case
when a small perturbative electric field is applied to the system , it is sufficient , to lowest order in the perturbation , to consider only the time evolution of the off - diagonal part of the density matrix,@xcite @xmath63 .
one then obtains from eq .
( [ eqn : bkgnd : tddmeom ] ) @xmath64 where @xmath65 , and @xmath66 . here
, @xmath67 denotes the matrix elements of @xmath68 , defined similarly to eq .
( [ hmat ] ) .
there is no external perturbation in eq .
( [ eqn : bkgnd : rhocvintime ] ) ; we assume free propagation of the system .
our reasoning is that individual excitations can be viewed as the eigenmodes of the system and do not depend on the specific form of the perturbative field .
fourier transformation of eq .
( [ eqn : bkgnd : rhocvintime ] ) gives @xmath69 where the factor 2 accounts for the spin , and @xmath70    equation ( [ eqn : bkgnd : rhocv ] ) is equivalent to the sma for finite systems . while the sma only refers to the transition between one individual occupied and one individual unoccupied orbital , eq . ( [ eqn : bkgnd : rhocv ] ) considers the transitions between the valence and the conduction band as a whole . ignoring the coupling between excitations and de - excitations by setting @xmath71 ( i.e. , making the tda ) , one arrives at the solid - state analog of the spa : @xmath72\rho_{{\mathbf{k}}'}^{cv}(\omega ) = \omega\rho_{\mathbf{k}}^{cv}(\omega ) .
\label{eqn : bkgnd : rhocvspa}\ ] ] this is the central equation which we will use to describe excitonic effects .
it requires as input the kohn - sham bloch functions for the valence and the conduction band of an insulator or semiconductor . to keep things as simple as possible
, we will consider the kronig - penney ( kp ) model@xcite rather than a real material .
the kp model is a 1d noninteracting system with a periodic potential of square wells . within the unit cell
@xmath73 $ ] , the potential is @xmath74        a typical example for the band structure of the kp model is plotted in fig .
[ fig : bkgnd : kppot ] . despite its simple appearance ,
the kp model is very versatile : by varying the values of the lattice constant @xmath75 , barrier width @xmath76 , and barrier height @xmath77 , a wide range of band gaps and band curvatures can be achieved .
a square - well potential of finite depth does not support an infinite number of bound states as the coulomb potential does ; but this , in fact , closely reflects the reality of the effective potential felt by the valence electrons in a solid , which is relatively shallow due to the screening of the bare nuclear charges by the core electrons . in practice , many solid - state calculations account for this screening by using muffin - tin potentials or pseudopotentials .
@xcite the kp model can be viewed as an elementary version of this approach .    in the following , we always choose the first two bands to be fully occupied and the higher bands to be empty .
we then make the two - band approximation for band 2 and 3 , since the shapes of these bands resemble the highest valence band and lowest conduction band in direct - gap materials such as gaas .
the bands in the kp model are sufficiently well separated , so the two - band approximation is justified .    to establish a connection with tddft
, we assume the solution to the noninteracting kp model as our ground - state kohn - sham system .
in other words , the potential in eq .
( [ eqn : bkgnd : kppot ] ) represents the _ exact _ kohn - sham potential @xmath78 which corresponds to a physical system whose external potential @xmath79 is uniquely determined thanks to the hohenberg - kohn theorem of dft . for our purpose , it is not necessary to know what this external potential looks like .
the kohn - sham bloch functions can then be determined in an elementary fashion.@xcite    for comparison , we will also carry out bse calculations in our model system ( see appendix [ sect : app : bse ] for technical details ) .
bse calculations are typically based on ground - state quasiparticle states obtained from the gw method.@xcite this is because the single - particle gap in gw is usually closer to experiment than the approximate kohn - sham gap .
however , in our case this distinction is not important because we use the given kp band structure as input for both bse and tddft .    in summary , our minimal tddft model for excitons
consists of the following two ingredients :    \(1 ) a two - band model for the vertical transitions between the highest valence band and the lowest conduction band , see eq .
( [ eqn : bkgnd : rhocvspa ] ) ;    \(2 ) the band structure from a 1d kp model .
of course , the model is not complete without a choice for the xc kernel @xmath35 .
this will be discussed below .
the exact xc kernel @xmath80 is unknown and must be approximated ; we restrict ourselves to adiabatic kernels that have no frequency dependence .
the adiabatic local - density approximation ( alda ) , as well as all semilocal , gradient - corrected xc kernels , are known to be unable to describe excitonic effects.@xcite the exact xc kernel has a long - range decay of @xmath81 , which is absent in all ( semi)local xc kernels derived from the uniform electron gas .
this long - range part is thought to be essential for excitons.@xcite    the long - range behavior of @xmath35 depends on the dimensionality , and in our 1d model system we define the following long - ranged ( or ` soft - coulomb ' ) xc kernel : @xmath82 we also consider an extremely short - ranged contact xc kernel : @xmath83 these model xc kernels depend on the constants @xmath84 and @xmath85 , which we will treat as fitting parameters in the following .
the idea is to tune the parameters in the model xc kernels so that bound excitons are produced , and to align the lowest bound exciton in the tddft spectrum with the lowest bound exciton in the bse spectrum .
results for the imaginary part of the dielectric function are presented in figs .
[ fig : results : tddftoneexciton ] and [ fig : results : tddfttwoexcitons ] .
we find that both the long - ranged @xmath86 and the short - ranged @xmath87 produce bound excitons , and thus , strictly speaking , the long - range behavior of the xc kernel is not really required for excitonic effects .
the bse results in figs . [ fig : results : tddftoneexciton ] and [ fig : results : tddfttwoexcitons ] show several identifiable bound excitons , in agreement with the rydberg series predicted by the wannier model ; the number of visible bound excitons somewhat depends on the numerical resolution in momentum space .    for the kp model parameters of fig .
[ fig : results : tddftoneexciton ] , we find that the adiabatic tddft can only bind a single excitonic state .
for other kp parameters ( specifically those in which the lowest conduction band is above the barrier ) , tddft produces two excitons , see fig .
[ fig : results : tddfttwoexcitons ] , which agree well with the lowest two excitons in bse .
there are additional , higher - lying bound excitons in bse which are very faint and difficult to resolve numerically .
for all the kp systems we tested , we never found more than two bound excitons with tddft .
this indicates the limitations of the adiabatic xc kernels used here .     +    as mentioned in sec .
[ sect : model ] , the wannier model does not clearly demonstrate that excitons are collective excitations .
since the wannier model assumes a single electron - hole pair picture , one can not immediately see that excitonic excitations are composed of a coherent superposition of many single - particle excitations . in our minimal model
, we solve the eigenvalue equation ( [ eqn : bkgnd : rhocvspa ] ) , and the eigenvectors @xmath88 ( which depends on @xmath21 parametrically ) describe how the transitions between noninteracting orbitals form the transitions in the interacting system .
@xmath89 is the percentage of a noninteracting transition in the transition of the interacting system .
two typical cases are plotted in figs .
[ fig : results : excitoniceigenvect ] and [ fig : results : nonexcitoniceigenvect ] .
figure [ fig : results : excitoniceigenvect ] clearly shows that excitons are collective excitations which are formed by mixing a wide distribution of single - particle transitions . as expected , the lowest exciton eigenfunction is nodeless and the second excitonic eigenfunction has a single node . with purely parabolic bands ,
the results from the wannier model would be recovered , as we will show below .
in contrast , the transitions in the continuum shown in fig .
[ fig : results : nonexcitoniceigenvect ] have a strong single - particle character ( the two peaks arise from the @xmath90 degeneracy in the kp model ) .    equation ( [ eqn : bkgnd : rhocvspa ] ) is equivalent to the spa for finite systems , which ignores the coupling between excitations and de - excitations ( tda ) .
we also investigated what happens when we do not make the tda , i.e. , when we work instead of eq .
( [ eqn : bkgnd : rhocvspa ] ) with the full equation for the two - band model , eq . ( [ eqn : bkgnd : rhocv ] ) .
as long as we describe relatively weakly bound excitons that are not too far below the band gap , we find that the difference between the two methods is very minor .
however , we also discovered that , under rare circumstances , eq .
( [ eqn : bkgnd : rhocv ] ) can lead to tddft excitonic binding energies that are purely imaginary . in our minimal model ,
such instabilities arise when the interaction strength @xmath12 in eq .
( [ eqn : results : fxcsc ] ) and ( [ eqn : results : fxccont ] ) increases so that the excitonic binding energy becomes greater than the band gap .
this situation is comparable to the well - known triplet instability in tddft , for which the tda generally leads to an overall better behavior;@xcite for excitons binding energies in our minimal model , we draw similar conclusions .
the bse scheme is commonly implemented within an adiabatic scheme ( see appendix [ sect : app : bse ] ) ; as we have seen , it produces a series of bound excitons .
since we assumed that the kp model is the kohn - sham ground state in tddft and the gw quasiparticle ground state in bse , the difference between tddft and bse becomes easily comparable , since the central equation to be solved have the same form , eq .
( [ eqn : bkgnd : rhocvspa ] ) .
the @xmath91 coupling matrices for tddft and bse are @xmath92 where @xmath93 is the hartree kernel ( the 1d soft - coulomb interaction ) , and @xmath94 is the screened interaction .
aside from the change from @xmath35 to @xmath94 , the most prominent difference between bse and tddft is the order of the indices for @xmath94 in eq .
( [ eqn : results : fbse ] ) . since the noninteracting ground - state wave functions have the bloch form ( [ eqn : bkgnd : blochform ] ) , we can see from eq .
( [ eqn : bkgnd : fijmn ] ) that @xmath95 has a strong @xmath96 dependence ; in fig .
[ fig : results : bsecoupling ] this shows up as a dominance along the diagonal .
by contrast , this @xmath96 dependence is clearly absent in @xmath97 , as demonstrated in fig .
[ fig : results : tddftcoupling ] .
@xmath98 and @xmath99 with only vertical transitions can be expressed in momentum space as @xmath100 @xmath101 the xc matrix ( [ eqn : results : fxctrans ] ) only depends on the long - range ( @xmath102 ) behavior of its momentum space representation @xmath103 , while the @xmath94 matrix ( [ eqn : results : wtrans ] ) also depends on other @xmath104 values in its momentum space representation @xmath105 .
it is impossible to find an adiabatic @xmath35 that reproduces the bse coupling matrix as in fig .
[ fig : results : bsecoupling ] , since @xmath105 has an extra degree of freedom over @xmath106 .
one can only hope to reproduce a portion of the bse coupling matrix with adiabatic tddft ( as pointed out in ref . ) , or make the xc kernel frequency dependent so that the information from the @xmath104-dependence in @xmath105 is mapped into the frequency dependence in @xmath107 .    considering the nature of the objects involved in this mapping , a highly nontrivial frequency dependence in @xmath35
is required to reproduce a series of bound excitons .
for example , one can easily construct an @xmath35 which reproduces a given series of bound excitons by using a different contact kernel in the region @xmath108 surrounding each exciton at frequency @xmath109 : @xmath110\theta[(\omega_i+\eta_i^+)-\omega].\ ] ] here , the @xmath111 s are parameters which are adjusted so that a tddft calculation with the frequency - independent kernel @xmath112 would produce @xmath109 as the lowest excitonic binding energy .
such an @xmath35 is of course completely ad hoc , but the fact that the excitonic series can be reproduced in this way demonstrates that the inclusion of the frequency dependence would greatly improve the flexibility of the tddft scheme .    on the other hand , within the adiabatic approximation
the characteristics of the @xmath113 coupling matrix are important for excitonic effects . to emphasize this ,
we now show that in very special cases the number of discrete excitonic eigenvalues can be derived .
consider a real matrix @xmath114 , where @xmath115 is a diagonal real matrix with @xmath116 , and @xmath117 has the symmetry @xmath118 . within second order perturbation theory ,
there is at most one discrete eigenvalue of @xmath114 in the limit where @xmath119 become continuous.@xcite though this case does not correspond to the matrices that would occur in real calculations , it indicates the close relationship between the properties of the coupling matrix and excitonic effects .
it is also possible to derive properties of the discrete eigenvalues if @xmath120 and @xmath121 are completely decoupled in the xc kernel .
owing to the symmetry @xmath122 implied by eq .
( [ eqn : bkgnd : fijmn ] ) , such separable kernels can only have the form @xmath123 for an excitation below the band gap with frequency @xmath21 , we can show@xcite that it must satisfy @xmath124 where the sum is carried out over the first brillouin zone ( fbz ) .
equation ( [ eqn : results : eqnsepkernel ] ) shows that eq .
( [ eqn : results : sepkernel ] ) must have the negative sign in order to have bound excitons . the left - hand side of eq .
( [ eqn : results : eqnsepkernel ] ) is monotonically increasing with @xmath21 , so for separable kernels of the form of eq .
( [ eqn : results : sepkernel ] ) , there is only one bound excitonic solution .    as shown in fig .
[ fig : results : tddftcoupling ] , tddft coupling matrices lack the strong dependence of @xmath125 as in bse coupling matrices . expanding the tddft coupling matrices into a power series of separable matrices and truncating at the first order would be a reasonable approximation , explaining why tddft produces fewer bound excitons ( if any at all ) than many - body methods such as bse .
the contact xc kernel and the soft - coulomb xc kernel in eqs .
( [ eqn : results : fxcsc ] ) and ( [ eqn : results : fxccont ] ) have the following simple form in momentum space : @xmath126 where @xmath127 , @xmath128 and @xmath128 are reciprocal lattice vectors , and @xmath129 is the modified bessel function of the second kind .
it is customary to refer to the matrix elements where @xmath130 as ` head ' , @xmath131 or @xmath132 as ` wings ' , and @xmath133 as ` body ' .
the 3d coulomb potential has the form @xmath134 in momentum space .
however , in 1d systems there is no real coulomb interaction which behaves as @xmath135 for @xmath136 , and one has to use the soft - coulomb interaction instead .
though there are many flavors of the soft - coulomb interaction , they all have the same @xmath137 behavior for @xmath136 .
however , the linear response function @xmath138 always behaves as @xmath135 for @xmath136 and does not depend on the dimensionality .
this renders quantities like the macroscopic dielectric function ( [ eqn : bkgnd : macrodielec ] ) ill - defined for strictly 1d systems .
furthermore , the bootstrap xc kernel@xcite and other xc kernels that depend on the cancellation of the 3d coulomb @xmath135 singularity will not work as designed in strictly 1d and 2d systems .
therefore @xmath139 shown in fig .
[ fig : results : tddftoneexciton ] and [ fig : results : tddfttwoexcitons ] are calculated at a small but finite @xmath140 .
the coupling matrix @xmath141 can be written in momentum space as @xmath142 \nonumber\\ & & \times{\left\langle c , k \left| e^{igx } \right| v , k \right\rangle}{\left\langle v , k ' \left| e^{-ig'x } \right| c , k ' \right\rangle}. \label{eqn : app : detail : fvcvc}\end{aligned}\ ] ]    for any xc kernel that behaves as @xmath135 for @xmath136 , one can further simplify the calculation by ignoring the so - called local field effects,@xcite i.e. instead of summing over @xmath128 and @xmath143 in eq .
( [ eqn : app : detail : fvcvc ] ) , only the head is considered . in 3d systems ,
a prominent example is the long - range kernel @xmath144 , which is obtained as an effective xc kernel with only head matrix elements from inverting the bse of contact excitons.@xcite    on the other hand , any xc kernel that diverges more slowly than @xmath135 for @xmath136 changes the spectrum only through the local field effects , i.e. all @xmath128 and @xmath143 must be summed in eq .
( [ eqn : app : detail : fvcvc ] ) . in other words ,
effective xc kernels with only the head are not feasible in strictly 1d systems due to the asymptotic behavior of the soft - coulomb potential discussed above . for 1d systems with @xmath131
, we have @xmath145 and @xmath35 s with 1d long - range behavior such as the soft - coulomb kernel also behave as @xmath146 .
considering eq .
( [ eqn : app : detail : fvcvc ] ) , these asymptotic properties imply that the head and wing contributions to @xmath91 always vanish in strictly 1d systems for physically meaningful xc kernels . due to these dimensionality restrictions ,
the xc kernel changes the strictly 1d system only through the local field effects .    in 3d the head contribution to the coupling matrix
@xmath147 is much more important than the local field effects , which is the reason that long - range kernels ( with nonzero head ) work much better than local xc kernels ( with vanishing head ) such as alda . in our strictly 1d model system , the head contribution is zero even for the bse , and thus the long - range kernel does not outperform local kernels such as the contact kernel .
these peculiarities only occur when one considers strictly 1d and 2d systems . in a more realistic picture ,
one encounters quasi-2d systems@xcite and quasi-1d systems ( such as quantum wires with finite radius or nanotubes@xcite ) , in which the movement of electrons is confined in certain directions such that the transverse motion can be averaged in comparison with the longitudinal motion .
though these systems show low - dimensional characteristics in various properties due to confinement , in the limit of @xmath136 they eventually differ from strictly low - dimensional systems .
our minimal model and the wannier exciton picture can be connected by considering the fourier transform of eq .
( [ eqn : bkgnd : rhocvspa ] ) .
we define an effective two - body potential @xmath148 via the fourier transform of @xmath149 : @xmath150 where @xmath151 is a direct lattice vector .
the fourier transform of the density matrix is @xmath152 since wannier excitons extends over many lattice constants , we approximate @xmath151 as a continuous variable @xmath153 . assuming the effective mass approximation , eq .
( [ eqn : bkgnd : rhocvspa ] ) becomes @xmath154 where @xmath155 is the excitonic binding energy , and the integration is carried out over all space .
we call eq .
( [ eqn : bkgnd : tddftwannier ] ) the tddft wannier equation , since it has the same form as eq .
( [ eqn : bkgnd : wannier ] ) . with proper choice of the approximated xc kernel ,
the nonlocal effective electron - hole interaction potential @xmath148 supports bound excitonic states .    since the bse and the tddft are formally similar within the minimal model , eq .
( [ eqn : bkgnd : rhocvspa ] ) can also be applied to the bse results .
[ fig : results : veh ] shows the effective interaction potential @xmath148 for tddft and bse .
\(a ) + ( b )    the tddft wannier equation provides an intuitive way of describing the effective nonlocal electron - hole interaction , and of explaining why adiabatic tddft usually has fewer excitons than bse and the wannier model . however , in most cases the tddft wannier equation is not suitable for quantitative use due to the approximations involved .
the approximation where we take the lattice vector @xmath151 as a continuous variable assumes that the exciton radius is much larger than the lattice constant ; this works fine in most cases we tested .
but the effective mass approximation where @xmath156 is approximated by @xmath157 is only good for transitions near the band gap , thus requiring these transitions of the noninteracting system to dominate the exciton , which is equivalent to the exciton extending over many lattice constants .
one obtains the @xmath158 term in eq .
( [ eqn : bkgnd : tddftwannier ] ) from @xmath157 in the limit where the lattice constant @xmath159 , and this approximation is not valid for most systems .
although @xmath148 is a nonlocal potential , in most cases we find that the @xmath148 s for both tddft and bse are dominated by the diagonal part , so the exciton problem is in analogy to one - body systems .
[ fig : results : vehdiag ] shows the diagonal part of @xmath148 , which can be taken as the effective one - body potential .
the wannier model in 1d has the soft - coulomb potential , which supports an infinite number of bound excitons ( the soft - coulomb interaction is fitted so that the binding energy of the first exciton matches that of the bse ) .
we find in general that the diagonal parts of @xmath148 for both bse and tddft are much more shallow than the soft - coulomb potential , and the tddft one is more narrow than the bse one .
thus , bse and tddft are not able to produce a complete excitonic rydberg series , and tddft in general produces fewer bound excitons than the bse .
the tddft wannier equation is not suitable for quantitative use for most of our model systems , despite the success of the wannier model in describing real semiconductors.@xcite since the approximations involved in eq .
( [ eqn : bkgnd : tddftwannier ] ) require that the exciton radius is large compared to the lattice constant , this suggests that this discrepancy is due to the special nature of 1d systems : namely , for similar effective masses the exciton radius in 1d is much smaller than in 3d and 2d.@xcite
the purpose of this paper was to construct a transparent and accessible minimalist model system that produces excitonic effects in a non-_ad hoc _ fashion using tddft .
the model , as presented here , is not intended to be a testing ground for xc kernels .
thus despite the dimensionality restrictions for the strictly 1d system , our results in sec .
[ sect : results ] carry over to 3d systems .    with our minimal model ,
we show that adiabatic tddft is capable of producing bound excitons through the local field effect even when the xc kernel is local in space , provided the strength of the kernel is strong enough .
this statement is still true in 3d ; however , due to the non - vanishing head contribution of the exact @xmath35 , we expect that the deviation of the effective interaction strength of a local @xmath35 from the real , nonlocal @xmath35 becomes larger than the in our strictly 1d model . in this sense
the long - range kernel , though very favorable , is not a necessary condition for excitonic effects .
we show the connection between tddft and the wannier model for excitons by deriving the tddft wannier equation , which describes a real - space system featuring a nonlocal effective electron - hole interaction .
such a connection intuitively demonstrates how adiabatic tddft generally produces fewer bound excitons than bse , and does not have a complete rydberg series .
the eigenvectors of the excitonic excitations in the minimal model clearly show their collective nature , which is not obvious from the wannier model alone .
excitonic instabilities may show up in tddft with approximate xc kernels , and this suggests that the tda tends to be more reliable for excitons than the formally exact method .    the frequency dependence of the exact xc kernel , @xmath160 , is usually ignored . despite the fact that adiabatic xc kernels have met with considerable recent success in producing optical spectra of insulators and semiconductors ( see the discussion in the introduction ) , they are incapable of producing excitonic rydberg series .
our model system gives an explanation for why this is the case .
this failure of the adiabatic approximation for @xmath35 is quite different from that which is responsible for the inability of adiabatic tddft to produce double excitations in finite systems or certain classes of charge - transfer excitations.@xcite this calls for continuing efforts in the search for nonadiabatic xc kernels for excitons .
this work was supported by nsf grant no . dmr-1005651 .
electrons and holes near the fermi surface are well described in the quasiparticle picture .
the quasiparticle green s function @xmath128 is related to that of the noninteracting system , @xmath161 , through the use of the self - energy @xmath162 : @xmath163 where the arguments denote sets of space and time variables .
a widely used approximation for the self - energy is the gw approximation : @xcite @xmath164 where @xmath94 is the screened interaction , @xmath165 @xmath166 is the bare coulomb interaction , and the inverse dielectric function @xmath167 is obtained as @xmath168 the linear response function @xmath138 can be calculated by its lehmann representation : @xmath169 where @xmath170 are quasiparticle states , @xmath171 are quasiparticle energies , and @xmath172 are occupation numbers . in practice ,
evaluating @xmath138 through eq .
( [ eqn : app : bse : chi ] ) can be quite time - consuming , and @xmath138 is often calculated with the plasmon - pole model.@xcite however , our minimal model is simple enough to allow us to use eq .
( [ eqn : app : bse : chi ] ) directly .
the gw quasiparticle green s function obtained from eq .
( [ eqn : app : bse : gwsigma ] ) misses important dynamical many - body effects , such as the electron - hole ( excitonic ) interaction .
the two - particle green s function includes these effects .
the bse@xcite describes the relation between the four - point polarization function @xmath173 of an interacting system and the corresponding object of the quasiparticle system : @xmath174 in which @xmath175 is @xmath176 and assuming the gw approximation , the kernel @xmath129 is @xmath177 here , @xmath178 denotes the coulomb interaction with the long - range part removed.@xcite in practice the bse is often solved in the transition space , which is spanned by single - particle excitations .
a four - point function such as @xmath179 then becomes @xmath180 where the @xmath32 s can be any complete basis set .
( [ eqn : app : bse : l ] ) in the transition space becomes @xmath181^{-1}_{(ij)(mn)}(f_m - f_n),\ ] ] where the excitonic hamiltonian matrix is @xmath182 we make the adiabatic approximation for @xmath183 and arrive at the following eigenvalue problem : @xmath184 and @xmath185 can be expressed in terms of these eigenvectors by @xmath186 only the transitions between the valence and conduction bands contribute . within our two - band model , the excitonic hamiltonian has the following block matrix form : @xmath187 ignoring the off - diagonal part in eq .
( [ eqn : app : bse : hexc ] ) is equivalent to the tda .    as shown in sec .
[ sect : results ] , it is possible that instabilities show up in the full bse results when the underlying ground - state calculation is not exact .
such instabilities in the minimal model are an artifact originating from the assumption that the solution of the kp model constitutes the ground state of the many - body system .
however , this is not a matter of great concern in practice .    in principle
, the transition space spans all possible combinations of valence and conduction orbitals , including nonvertical transitions connecting different bloch wavevectors .
the kernel @xmath188 of the bse in momentum space , eq .
( [ bsekernel ] ) , has the following ingredients : @xmath189 and @xmath190 only the excitations with the same momentum transfer @xmath191 are coupled due to the @xmath192 functions in eq .
( [ eqn : app : bse : v ] ) and ( [ eqn : app : bse : w ] ) , so we only need to include vertical transitions in the calculations for optical properties .