a common problem in statistical physics is that the probability distribution functions ( pdfs ) are always too complicated for practical calculations and we need to replace them by more tractable approximations .
a possible solution is to identify a family of trial distributions @xmath0 , where x are parameters characterized systems and select the member of the family that is closest  to the exact  distribution @xmath1 .
the problem , of course , is that it is not clear what one means by ` closest ' .
one could minimize    @xmath2 ^{2},\ ] ]    but why this particular functional and not another ? and also , why limit oneself to an approximation by a single member of the trial family ? why not consider a linear combination of the trial distributions , some kind of average over the trial family ?
but then , how should we choose the optimal weight assigned to each @xmath3 ? we propose to tackle these questions using the method of maximum relative entropy ( which we abbreviate as me ) and information geometry @xcite .
the me method , which is developed in @xcite-@xcite , has historical roots in the earlier method of maximum entropy that was pioneered by e. t. jaynes and is commonly known as maxent @xcite .
the me method is designed for updating probabilities from arbitrary priors for information in the form of arbitrary constraints and it includes bayes rule and the older maxent as special cases @xcite , @xcite .    the purpose of this paper is to develop a me based method to generate optimal approximations ( brief accounts of some of the results discussed below have previously been presented in @xcite and @xcite ) .
the general formalism , which is the main result of this paper , is developed in section [ sec2 ] . in section
[ sec2 - 1 ] we justify the use of relative entropy as the unique and natural criterion to select the preferred approximation , which is labeled by some parameters .
the optimal approximation is obtained in section [ sec2 - 3 ] by marginalizing over the variational parameters .
a suitably weighed average over the whole family of trial distributions with the optimal weight provides an optimal approximation .
in the second part of the paper we demonstrate the proposed formalism by applying it to simple classical fluids , a well studied field in the past barkerhenderson76-@xcite . to approximate the behavior of simple fluids we chose trial distributions that describe hard spheres @xcite-@xcite ( section  [ sec3 ] ) .
the me formalism is first used ( section  [ sec4 ] ) to select the preferred value of the hard - sphere diameter .
this is equivalent to applying the bogoliubov variational principle and reproduces the results obtained by mansoori et al .
@xcite whose variational principle was justified by a very different argument .
an advantage of the variational or the me methods over the perturbative approaches such as barker and henderson ( bh ) @xcite and of weeks , chandler and anderson ( wca ) @xcite is that there is no need for ad hoc criteria dictating how to separate the intermolecular potential into a strong short range repulsion and a weak long range attraction . on the other hand ,
a disadvantage of the standard variational approach is that it fails to take the softness of the repulsive core into account . at high temperatures
this leads to results that are inferior to the perturbative approaches .    in the standard variational
a single preferred value of the hard - sphere diameter is selected .
but , as discussed in @xcite and caticha00 , in the me method non - preferred values are not completely ruled out .
this allows us ( in section  [ sec5 ] ) to marginalize over hard sphere diameter to obtain an optimal hard - sphere approximation with a suitable weighting .
that leads to significant improvements over the standard variational method .    in section  [ sec6 ]
we test our method by comparing its predictions for a lennard - jones model for argon with molecular dynamics simulation data ( verlet68 , @xcite ) .
we find that the me predictions for thermodynamic variables and for the radial distribution function are considerable improvements over the standard bogoliubov variational result , and are comparable to the perturbative results barkerhenderson76 @xcite .
( for a recent discussion of some of the strengths and limitations of the perturbative approach see @xcite . ) . despite the shortcomings of perturbation theory
it remains very popular because it provides quantitative insights at a much lower computational cost than dynamical simulations . for recent applications to the glass transition and other more complex systems
see @xcite - @xcite .
although in this work the me is not applied to such complex problems one may fully expect that the information theory based me method will yield insights not only the thermodynamic behavior of complex systems but also about the approximation methods needed to analyze them . finally , our conclusions and some remarks on further improvements are presented in section  [ sec7 ] .
consider a system with microstates labeled by @xmath4 ( for example , the location in phase space or perhaps the values of spin variables ) .
let the probability that the microstate lies within a particular range @xmath5 be given by the intractable canonical distribution @xmath6where @xmath7 is partition function , in which @xmath8 is the hamiltonian and free energy @xmath9 is defined through this partition function .
the goal is to generate an approximation @xmath10 that is optimal in the sense that is `` closest '' to the `` exact '' distribution @xmath11 .
it includes two steps .
the first is to select a preferred trial distribution from a family of trials @xmath12 , each member of the trial family being labeled by one or more parameters @xmath13 .
more generally , one could define the trial family in a non - parametric way by specifying various constraints .
the second is to marginalize over parameters @xmath13 to obtain an optimal distribution .      *
relative entropy as the selection criterion . * the selection of the preferred approximation is achieved by ranking the distributions @xmath10 according to increasing _ preference _ , a real number @xmath14 $ ] which we call the `` entropy '' of @xmath15 .
the numbers @xmath16 $ ] are such that if @xmath17 is preferred over @xmath18 , then @xmath19>\mathcal{d}[p_{2}]$ ] .
thus , by design , the `` preferred '' approximation @xmath15 is that which maximizes the entropy @xmath16 $ ] .
next we determine the functional form of @xmath16 $ ] .
this is the general rule that provides the criterion for preference ; in our case it defines what we mean by the `` closest '' or `` preferred '' approximation .
the basic strategy @xcite is one of induction : ( 1 ) if a general rule exists , then it must apply to special cases ; ( 2 ) if in a certain special case we know which is the best approximation , then this knowledge can be used to constrain the form of @xmath16 $ ] ; and finally , ( 3 ) if enough special cases are known , then @xmath20 $ ] will be completely determined .
the known special cases are called the axioms  of me and they reflect the conviction that whatever information was originally codified into the exact @xmath11 is important and should be preserved
. the selected trial distribution should coincide with the exact one as closely as possible and one should only tolerate those minimal changes that are demanded by the information that defines the family of trials .
three axioms and their consequences are listed below . detailed proofs and more extensive comments are given in caticha03 and @xcite .
* axiom 1 : locality*. _ local information has local effects .
_ if the constraints that define the trial family do not refer to a certain domain @xmath21 of the variable @xmath4 , then the conditional probabilities @xmath22 need not be revised , @xmath23 .
the consequence of the axiom is that non - overlapping domains of @xmath4 contribute additively to the entropy : @xmath20=\int dq\,f(p(q))$ ] where @xmath9 is some unknown function .
* axiom 2 : coordinate invariance . * _ the ranking should not depend on the system of coordinates . _ the coordinates that label the points @xmath4 are arbitrary ; they carry no information .
the consequence of this axiom is that @xmath20=\int dq\,p(q)f(p(q)/m(q))$ ] involves coordinate invariants such as @xmath24 and @xmath25 , where the function @xmath26 is a density , and both functions @xmath27 and @xmath28 are , at this point , unknown .
next we make a second use of axiom 1 ( locality ) .
when there are no constraints at all and the family of trials includes the exact @xmath11 the selected trial should coincide with @xmath11 ; that is , the best approximation to @xmath11 is @xmath11 itself .
the consequence is that up to normalization the previously unknown density @xmath26 is the exact distribution @xmath11 .
* axiom 3 :  consistency for independent subsystems*. _ when a system is composed of subsystems that are independent it should not matter whether the approximation procedure treats them separately or jointly .
_ specifically , if @xmath29 , and the exact distributions for the subsystems , @xmath30 and @xmath31 , are respectively approximated by @xmath32 and @xmath33 , then the exact distribution for the whole system @xmath34 should be approximated by @xmath35 .
this axiom restricts the function @xmath28 to be a logarithm .
the overall consequence of these axioms is that the trial approximations @xmath36 should be ranked relative to the exact @xmath11 according to their ( relative ) entropy , @xmath37=-\int dq\,p(q)\log \frac{p(q)}{p(q)}.   \label{s[p|p]}\]]the derivation has singled out the relative entropy @xmath38 $ ] as _ the unique functional to be used for the purpose of selecting a preferred approximation_. other functionals , may be useful for other purposes , but they are not a generalization from the simple cases described in the axioms above
.    * remark . *
suppose a member of a family of trial canonical distributions @xmath12 with hamiltonian @xmath39 that are conditional probability distributions and depend on parameters @xmath40 , are given by @xmath41where @xmath42 , in which free energy @xmath43 is also defined .
the preferred trial is then selected by maximizing @xmath44 $ ] . substituting eq .
( [ p(q ) ] ) and eq .
( [ ptheta(q ) ] ) into eq . ( [ s[p|p ] ] )
gives , @xmath45=\beta \left ( \langle h_{\theta } -h\rangle _ { \theta } -f_{\theta } + f\right ) \
, ,   \label{s[ptheta|p]}\]]where @xmath46 refers to averages computed with the distribution @xmath47 .
the inequality @xmath44\leq 0 $ ] , can then be written as @xmath48thus , maximizing @xmath44 $ ] is equivalent to minimizing the quantity @xmath49 .
this form of the variational principle and its use to generate approximations is well known .
it is usually associated with the name of bogoliubov @xcite and it is the main technique to generate mean field approximations for discrete systems of spins on a lattice .      the extent to which the preferred @xmath13 is preferred over other values ( @xcite , @xcite )
is expressed by the probability of @xmath13 , @xmath50 .
the original me problem of assigning a probability to @xmath4 is now broadened into assigning probabilities to both @xmath4 and @xmath13 .  in this section
we use me again to find the preferred joint distribution @xmath51 .
note that this is the kind of problem where the bayesian interpretation of probabilities is essential . within a frequentist interpretation it makes no sense to talk about @xmath50 or about @xmath52 because @xmath13 is not a random variable ; the value of @xmath13 is unknown but it is not random .    to proceed we must ask a question . what is the prior distribution , that is , what do we know about @xmath4 and @xmath13 before the trial family is specified ?
the joint prior @xmath53 can be expressed according to the product rule as , @xmath54 , where @xmath55 is conditional probability of observing system in state @xmath4 given parameter @xmath56 .
our goal is to determine the preferred @xmath57 that is closest to the prior @xmath53 that reflects our initial knowledge about @xmath4 and ignorance about the @xmath13s .
initially we know nothing about @xmath13 , not even how it is related to @xmath4 .
the prior that represents this state of knowledge is a product@xmath58indeed , when @xmath53 is a product no correlations between @xmath13 and @xmath4 are introduced which means that information about @xmath4 tells us nothing about @xmath13 and vice versa .
the first factor in @xmath53 reflects our prior knowledge about @xmath4 : the distribution for @xmath4 is known to be the exact @xmath11 .
the second factor reflects our complete ignorance about @xmath13 : we choose @xmath59 to be as uniform as possible .
our method applies whether @xmath13 is a discrete or a continuous variable .
when @xmath13 is a continuous variable .
then the uniform distribution @xmath60 is such that makes equal volumes in @xmath13 space equally likely . to define these volumes we apply method of information geometry @xcite and note that distances in @xmath61-space are uniquely defined because the @xmath13s are labels on probability distributions .
cases where @xmath13 is a discrete variable are simpler .
the relevant entropies involve sums over @xmath62 rather than integrals and the natural uniform distribution is @xmath63 . in what follows we concentrate on the more challenging continuous case .
the unique distance between @xmath13 and @xmath64 is given by the fisher - rao metric @xcite , @xmath65 , where @xmath66accordingly , the volume of a small region @xmath67 is @xmath68 , where @xmath69 is the determinant of @xmath70 .
up to an irrelevant normalization , the distribution @xmath59 that is uniform in @xmath13 is given by @xmath71 .
the preferred approximation @xmath57 to the joint distribution @xmath72 is then obtained maximizing the entropy@xmath73=-\int dq\,d\theta \,p(\theta ) p_{\theta } ( q)\,\log   \frac{p(\theta ) p_{\theta } ( q)}{\gamma ^{1/2}(\theta ) p(q ) } , \label{sigma[p]}\]]by varying @xmath50 subject to @xmath74 .
the final result for the probability that @xmath13 lies within the small volume @xmath75 is @xmath76}\gamma ^{1/2}(\theta ) d\theta ,   \label{p(theta)}\]]where @xmath77 $ ] is given in eq .
( [ s[ptheta|p ] ] ) and @xmath78 is a normalization constant .
note also that the density @xmath79 $ ] is a scalar function and the presence of the jacobian factor @xmath80 makes eq .
( [ p(theta ) ] ) manifestly invariant under changes of the coordinates @xmath13 .
( [ p(theta ) ] ) expresses the degree to which values of @xmath13 away from the preferred value are ruled out ; it tells us that the preferred value of @xmath13 is that which maximizes the probability density @xmath79 $ ] .    finally , now that we have determined the preferred joint distribution @xmath81 we can marginalize @xmath13 and use the average @xmath82as the best approximation we can construct out of the given trial family .
this approximation is expected to be better than any individual @xmath47 for the same reason that the mean is expected to be a better estimator than the mode  it minimizes the variance .
this concludes the first part of our paper . to summarize : our main results consist in the justification of the relative entropy eq .
( [ s[p|p ] ] ) as the uniquely natural functional to select the preferred approximations and the derivation of a quantitative measure of the degree to which the various trials are preferred , eq . ( [ p(theta ) ] ) .
the final result for the best approximation is eq .
( [ pbar ] ) .
next we illustrate how this me formalism is used in a specific example : simple fluids .
a simple fluid composed of @xmath83 single atom molecules is described by the hamiltonian @xmath84where @xmath85 and the many - body interactions are approximated by a pair interaction , @xmath86 where @xmath87 . the probability that the positions and momenta of the molecules lie within the phase space volume @xmath88 is given by canonical distribution @xmath89where @xmath90 . for fluids dominated by pair interactions most thermodynamic quantities of interest
can be written in terms of the one- and two - particle density distributions @xmath91where @xmath92and @xmath93the two - particle correlation function , @xmath94measures the extent to which the structure of liquids deviates from complete randomness .
if the fluid is homogeneous and isotropic @xmath95 and @xmath96 where @xmath97 is the bulk density and @xmath98 is the radial distribution function ( rdf ) .
then , the pressure is given by @xmath99where @xmath100 @xcite-@xcite .      to account for the short - distance repulsion we consider a family of trials composed by distributions that describe a gas of hard spheres of diameter @xmath101 .
for each @xmath102 the hamiltonian is @xmath103with @xmath104where @xmath105and the corresponding probability distribution is @xmath106the partition function and the free energy @xmath107 are @xmath108 @xmath109 @xmath110 two objections that can be raised for choosing @xmath111 as trials are , first , that they do not take the long - range interactions into account ; and second , that the actual short range potential is not that of hard spheres .
these are points to which we will return later .
a third objection , and this is considerably more serious , is that the exact hard - sphere rdf is not known .
however , it can be calculated within the approximation of percus and yevick ( py ) for which there exists an exact analytical solution ( @xcite , @xcite , @xcite ) which is reasonably simple and in good agreement with numerical simulations over an extended range of temperatures and densities , except perhaps at high densities .
there are several successful proposals @xcite to improve upon the py rdf but they also represent an additional level of complication .
the simpler py rdf is sufficiently accurate for our current objective  to illustrate the application and study the broad features of the me approach .
the py rdf can be written in terms of the laplace transform of @xmath112 @xcite , @xmath113 } ,   \label{g(s)}\]]where @xmath114 is a dimensionless variable @xmath115 , @xmath116 , \]]@xmath117and @xmath118 is the packing fraction , @xmath119the rdf @xmath120 is obtained from the inverse transform using residues @xcite .
the equation of state can then be computed in two alternative ways , either from the `` pressure '' equation or from the `` compressibility '' equation but , since the result above for @xmath120 is not exact , the two results do not agree .
it has been found that better agreement with simulations and with virial coefficients is obtained taking an average of the two results with weights 1/3 and 2/3 respectively .
the result is the carnahan - starling equation of state , @xcite - kalikmanov02 @xmath121the free energy , derived by integrating the equation of state , is @xmath122 , \mathrm{\ } \label{hs - free energy}\]]where @xmath123 , and the entropy is @xmath124it must be remembered that these expressions are not exact .
they are reasonable approximations for all densities up to almost crystalline densities ( about @xmath125 ) .
however , they fail to predict the face - centered - cubic phase when @xmath118 is in the range from @xmath126 up the close - packing value of @xmath127 .
as discussed in section  [ sec2 ] , the trial @xmath128 that is closest  to the exact  @xmath129 is found by maximizing the relative entropy @xmath130 $ ] , eq .
( [ s[p|p ] ] ) , with @xmath131 given by eq .
( [ phs ] ) and @xmath132 given by eq .
( [ exact dist ] ) . according to eq .
( [ f_u - general ] ) , it is equivalent to minimize @xmath133over all diameters @xmath102 , where @xmath134 is computed with @xmath111 .
thus , the variational approximation to the free energy is @xmath135where @xmath136 is the preferred diameter .    to calculate @xmath137
use @xmath138 \mathrm{~,}\]]where @xmath139 .
but @xmath140 for @xmath141 while @xmath142 for @xmath143 , therefore @xmath144with @xmath145where we have assumed that the fluid is isotropic and homogeneous , @xmath146 , and introduced the hard - sphere rdf @xmath147note that the approximation does not consist of merely replacing the exact free energy @xmath9 by a hard - sphere free energy @xmath148 which neglects the effects of long range attraction ; @xmath9 is approximated by @xmath149 which includes attraction effects through the @xmath150 term in eq .
( [ f_u ] ) .
this addresses the first of the two objections mentioned earlier : the real fluid with interactions given by @xmath151 is not being replaced by a hard - sphere fluid .
the internal energy is approximated by @xmath152 and not by @xmath153 .    to calculate @xmath150
it is convenient to write it in terms of @xmath154 , the inverse laplace transform of @xmath155 , @xmath156for example , for a lennard - jones potential , @xmath157 ,   \label{lj pot}\]]we have @xmath158 .\]]then , using equations eq .
( [ f_u ] ) and ( [ g(s ) ] ) gives @xmath159finally , it remains to minimize @xmath137 in eq .
( [ f_u ] ) to determine the preferred diameter @xmath136 .
this is done numerically in an explicit example for argon in section  [ sec6 ] .
the me method as pursued in the last section has led us to determine a preferred hard - sphere diameter .
it fails to take the softness of the repulsive core into account . as discussed in section  sec2 - 3 , our best assessment of the distribution of @xmath160
is given by the marginal over @xmath102 , @xmath161 the corresponding best approximation to the rdf is obtained using eq .
( [ numberdensity ] ) , ( [ rdf - general ] ) , and ( [ rdf - hs ] ) @xmath162 by averaging over all hard - sphere diameters we are effectively describing a soft - core potential . since @xmath163 takes into account soft - core effects while @xmath164 does not , we expect that it will lead to improved estimates for all other thermodynamic quantities .
however , we should emphasize that the distribution over the hard - sphere diameters @xmath165 is not being introduced in an ad hoc way in order to fix  the variational method .
the introduction of @xmath165 is mandated by the me method ( section  [ sec2 - 3 ] ) .
the distribution of diameters is given by eq .
( [ p(theta ) ] ) @xmath166 } } { \zeta } \gamma ^{1/2}\left ( r_{d}\right ) dr_{d}=\frac{e^{-\beta f_{u}}}{\zeta _ { u}}\gamma ^{1/2}\left ( r_{d}\right ) dr_{d } ,   \label{pd(rd)}\]]where @xmath167 = \beta \left ( f - f_{u}\right ) $ ] , the partition functions @xmath78 and @xmath168 are given by @xmath169and the natural distance @xmath170 in the space of @xmath102s is given by the fisher - rao metric , @xmath171    a convenient way to calculate the fisher - rao metric is to express it as a second derivative of the entropy eq .
( [ s[p|p ] ] ) of @xmath172 relative to @xmath173 , @xmath174 \right\vert _ { r_{d}^{\prime } = r_{d } } , \label{det g}\]]where @xmath175 = \beta \left [ \left .
{ r_{d}^{\prime } } ^{r_{d}}-\langle \left
. u_{hs}\right\vert _
{ r_{d}^{\prime } } ^{r_{d}}\rangle _ { r_{d}^{\prime } } \right ] , \]]and @xmath176 is the average over @xmath177 . as we argued above eq .
( [ f_u ] ) the expectation of the potential energy @xmath178 vanishes because the product @xmath179 vanishes for both @xmath180 and @xmath181 .
similarly , @xmath182 when @xmath183 .
however , when @xmath184 the expectation @xmath185 diverges .
the divergence is a consequence of the unphysical nature of the hard - sphere model . for more realistic continuous potentials
the distance between @xmath186 and @xmath102 is the same as the distance between @xmath187 and @xmath102 .
we can then always choose @xmath188 and define @xmath189 in eq .
( [ det g ] ) as the limit @xmath190 .
then , using eq .
( [ hs - free energy ] ) for @xmath148 , we have @xmath191    to summarize , the distribution of diameters @xmath165 is given by eq .
( [ pd(rd ) ] ) with @xmath137 given by equations ( [ f_u ] ) , ( [ hs - free energy ] ) , ( [ f_u - lt ] ) and @xmath192 given by eq . ( [ det g - exact ] ) . our best approximation to the `` exact '' @xmath193
is the @xmath194 given in eq .
( [ pbar(q ) ] ) . however , there is a problem .
since the free energy @xmath137 is an extensive quantity , @xmath195 , for large @xmath83 the distribution @xmath196 is very sharply peaked at the preferred diameter @xmath136 . when choosing a _ single _ preferred diameter for a macroscopic fluid sample we find that me confers overwhelming probability to the preferred value .
this is not surprising .
the same thing happens when we calculate the global temperature or density of a macroscopic sample and yet local fluctuations can be important .
the question then , is whether local fluctuations are relevant to the particular quantities we want to calculate .
we argue that they are .    from the very definition of @xmath197 as the probability that given an atom at a certain place
another atom will be found at a distance @xmath198 , it is clear that @xmath197 refers to purely local behavior and should be influenced by local fluctuations . to the extent that the preferred diameter @xmath136 depends on temperature and density we expect that local temperature and/or density fluctuations would also induce local diameter fluctuations .    for the purpose of calculating @xmath197 the system is effectively reduced to the small number of atoms @xmath199 in the local vicinity of the reference atom at the origin . in order to develop a systematic , fully me method for the determination of the effective number of particles @xmath199
that are locally relevant , one needs to use trial probability distributions that allow inhomogeneities in the hard sphere diameters .
yet determination of such trial probability distributions requires further investigations .
since we are interested in demonstrating me formalism for optimizing approximations in this work only , we consider a rather simple approach to estimate @xmath199 . based on our me approach to a mean field approximation for fluids @xcite ,
we have shown the rdf is given by @xmath200for a sufficiently dilute gas , the results of the mean field approximation are comparable with experimental results .
the @xmath201 is approximately close to the true , exact rdf @xmath202 for @xmath203 , the lennard - jones parameter .
a fluid of hard spheres gives @xmath204 for @xmath205 and can not reproduce the behavior of eq .
( [ estimaterule ] ) . however , once we recognize that we can use a statistical mixture , eq .
( [ pbar(q ) ] ) , we can tune the size @xmath206 of the cell and thereby change the width of @xmath207 so that the rdf @xmath208 of eq .
( [ ghs bar ] ) reproduces the known short - distance behavior of eq .
( [ estimaterule ] ) .
one of the difficulties in testing theories about fluids against experimental data is that it is not easy to see whether discrepancies are to be blamed to a faulty approximation or to a wrong intermolecular potential .
this is why theories are normally tested against molecular dynamics numerical simulations where there is control over the intermolecular potential . in this section
we compare me results against simulation results @xcite for a fluid of monoatomic molecules interacting through a lennard - jones potential , eq .
( [ lj pot ] ) .
the parameters @xmath209 and @xmath210 ( the depth of the well , @xmath211 , and the radius of the repulsive core , @xmath212 , respectively ) are chosen to model argon : @xmath213 @xmath214 and @xmath215 .
* the free energy @xmath137 .
* figure [ fig1].(a ) shows the free energy @xmath216 as a function of hard - sphere diameter @xmath102 for argon at a fixed density of @xmath217 for different temperatures .
figure [ fig1].(b ) shows @xmath218 as a function of @xmath102 for several densities at fixed @xmath219 @xmath220 . since the critical point for argon is at @xmath221 @xmath220 and @xmath222 all these curves , except that at @xmath223 @xmath220 , lie well within the liquid phase .
the increase of @xmath216 for high values of @xmath102 is due to short range repulsion between the hard spheres described by @xmath224 .
the increase for low @xmath102 is due to the lennard - jones short - range repulsion as described by @xmath225 .
the preferred @xmath102 is that which minimizes @xmath137 and depends both on temperature and density .
the preferred diameter decreases as the temperature increases because atoms with higher energy can penetrate deeper into the repulsive core .
the dependence with density is less pronounced .    * the distribution of diameters @xmath165 . * in section  [ sec5 ] we argued that the effective number of molecules that is relevant to the local structure of the fluid is not the total number of molecules in the system @xmath83 , but a smaller number , @xmath206 . in fig [ fig3].(a ) we plot the distribution of diameters @xmath226 for different temperatures , for a fixed fluid density of @xmath227 , and for an arbitrarily chosen @xmath228 . as expected the distribution shifts to higher diameters as the temperature decreases .
notice also that the distribution becomes narrower at lower temperatures in agreement with the fact that a hard - sphere approximation is better at low @xmath229 @xcite .
figure [ fig3].(b ) shows that increasing @xmath199 ( with fixed density @xmath97 ) decreases the width of @xmath165 ( solid lines ) and induces a slight shift of the whole distribution .
this is due to the dependence @xmath230 of the fisher - rao measure @xmath231 in [ det g - exact ] .
figure [ fig3].(b ) also explores the influence of @xmath232 by comparing the actual distributions @xmath165 ( solid lines ) with the distributions @xmath233 ( dotted lines ) which are obtained by setting @xmath234 in eq .
( [ pd(rd ) ] ) . the effect of @xmath235 is to shift the distribution slightly to higher @xmath102 .
* the radial distribution function .
* we are finally ready to calculate the radial distribution @xmath98 for argon .
we start by estimating the number of molecules @xmath236 that are locally relevant ; as explained earlier we choose @xmath236 so that our best approximation @xmath208 , eq .
( [ ghs bar ] ) , reproduces the known short - distance behavior @xmath201 , eq .
( [ estimaterule ] ) , for @xmath237 . we have found that the estimates for @xmath238 need not be very accurate but that they must be obtained for each value of the temperature and density . in fig [ fig5 ]
we show an example of the short - distance behavior of @xmath239 for three values of @xmath199 at @xmath219 @xmath220 and @xmath217 ; using a chi - square fit in the range from @xmath240 to @xmath241 @xmath242 the selected best value of @xmath199 is around 38000 .    in figures [ fig6].(a)-(d )
we compare three different ways to calculate the rdf .
the solid line is verlet s molecular dynamics simulation verlet68 ; it plays the role of experimental data against which we compare our theory .
the dotted line is @xmath243 for the hard - sphere fluid with preferred diameter @xmath136 .
this curve , calculated from inverse of eq .
( [ g(s ) ] ) , is also the result of the variational method and coincides with the me result for a macroscopically large @xmath244 .
the dashed line is the averaged @xmath163 of the extended me analysis
. figures .
[ fig6].(a)-(c ) were plotted at three different temperatures @xmath219 , @xmath245and @xmath246 k at the density @xmath217 .
figure [ fig6].(d ) we changed the density and the temperature to @xmath247 and @xmath248 @xmath249 .
the agreement between the me curve and verlet s data is good .
the vast improvement over the simpler variational method calculation is clear .
one might be tempted to dismiss this achievement as due to the adjustment of the parameter @xmath199 but this is not quite correct : @xmath236  has not been adjusted , it has been calculated by fitting @xmath163 to optimal mean field rdf @xmath250 for @xmath203 . indeed , despite the fact that the hard - sphere trial solutions that we employ are mere approximations , the functional form of the whole curve @xmath163 in eq .
( [ ghs bar ] ) is reproduced quite well .
however , one may note that the agreement between the me prediction and verlet s data becomes worse when the fluid density is decreased or the temperature is increased .
the reason has been spelled out by in the studies of wca @xcite they demonstrate that both the repulsive and attractive forces contribute to the fluid structure when fluid is at low and moderate densities ( @xmath251 ) . however , when the fluid density is high enough ( @xmath252 ) , the repulsive force becomes dominant . because the hard - sphere approximation does not include the attractive force , the hard - sphere rdf does not take the attractive force into account , and this error propagates into our me prediction .
the same discrepancy is also revealed in the wca theory for low density @xcite .    *
the equation of state .
* finally we use the rdf to calculate the equation of state from the pressure equation , eq .
( [ eq of state ] ) . in fig [ fig10 ]
we compare the equation of state derived from the @xmath197 obtained from verlet s simulation with calculations using the eme and variational methods and the perturbative theories of barker and henderson @xcite and of weeks , chandler and anderson @xcite , at @xmath253 @xmath249 .
the eme results constitute a clear improvement over the plain variational calculation . for low densities
all four methods agree with each other but differ from the simulation .
a better agreement in this region would probably require a better treatment of two - particle correlations at long distances . at intermediate densities the best agreement
is provided by the eme and bh results , while the wca theory seems to be the best at high densities .
also shown in fig [ fig10].(a ) are experimental data on argon @xcite .
the discrepancy between the experimental curve and the verlet simulation is very likely due to the actual potential not being precisely of the lennard - jones type .    in fig
[ fig10].(b ) we plot the eme equation of state for three different isotherms ( @xmath254 , @xmath255 and @xmath256 ) . to compare to the simulation of hansen and verlet @xcite we plot @xmath257 ( rather than @xmath258 ) as a function of density @xmath259 because this kind of plot exhibits the characteristic van der waals loop that signals the liquid - gas transition as the temperature drops .
a more exhaustive exploration lies , however , outside the scope of this paper .
the goal of this paper has been to use the eme method to generate approximations and show that this provides a generalization of the bogoliubov variational principle .
this addresses a range of applications that lie beyond the scope of the traditional maxent . to test the method we considered simple classical fluids .
when faced with the difficulty of dealing a system described by an intractable hamiltonian , the traditional approach has been to consider a similar albeit idealized system described by a simpler more tractable hamiltonian .
the approach we have followed here departs from this tradition : our goal is not to identify an approximately similar hamiltonian but rather to identify an approximately similar probability distribution .
the end result of the eme approach is a probability distribution which is a sum or an integral over distributions corresponding to different hard - sphere diameters .
while each term in the sum is of a form that can be associated to a real hard - sphere gas , the sum itself is not of the form @xmath260 , and can not be interpreted as describing any physical system .
as far as the application to simple fluids is concerned the results achieved in this paper represent progress but further improvements are possible by using better approximations to the hard - sphere fluid and by choosing a broader family of trial distributions .
an important improvement would be to use trial probability distributions that allow inhomogeneities in the hard sphere diameters . this would lead to a systematic , fully eme method for the determination of the effective number of particles @xmath199 that are locally relevant .
many perturbative approaches to fluids had been proposed , and a gradual process of selection over many years of research led to the optimized theories of bh and wca .
the variational approach was definitely less satisfactory than these bestperturbation theories . with our work
, however , the situation has changed : the eme - improved variational approach offers predictions that already are competitive with the best perturbative theories . and , of course , the potential for further improvements of the eme approach remains , at this early date , far from being exhausted .
the authors acknowledge r. scheicher and c .- w .
hong for their valuable assistance and advice with the numerical calculations .
n. n. encov : _ statistical decision rules and optimal inference , transl . math .
monographs , vol .
53 _ , am .
( providence , 1981);l .
l. campbell , proc . am .
( 1986 ) 135 ; for applications to statistics see s. amari , _ differential - geometrical methods in statistics _ , springer - verlag ( new york , 1985 ) ; s. amari and h. nagaoka , _ methods of information geometry _ , oxford ( providence , 2000 ) ; for a brief derivation see a. caticha , change , time and information geometry , in : a mohammad - djafari ( ed ) , bayesian methods and maximum entropy in science and engineering , aip conf .
proc . * 568 * , 2001 ( arxiv : math - ph/0008018 ) .
a. caticha , relative entropy and inductive inference , in : g. erickson , y. zhai ( eds ) , bayesian inference and maximum entropy methods in science and engineering , aip conf . proc . * 707 * , 2004 ( arxiv : physics/0311093 ) ; a. caticha and a. giffin , updating probabilities , in : k.h .
knuth , a. caticha , a. giffin , c. c. rodriguez , j. l. j. center ( eds . ) , bayesian inference and maximum entropy methods in science and engineering , aip conf .
proc . * 954 * , 2007 ( arxiv : physics/0608185v1 ) .
a. caticha , lectures on probability , entropy , and statistical physics , in maxen08 , so paulo , 2008 ( arxiv.org/abs/0808.0012 ) ; an application to data analysis is given in a. caticha and r. preuss , phys .
e * 70 * ( 2004 ) 0461271 .
e. t. jaynes , phys . rev .
* 106 * ( 1957 ) 620 ; e. t. jaynes phys
* 108 * ( 1957 ) 171 ; e. t. jaynes in : r d rosenkrantz ( ed ) , _
e. t. jaynes : papers on probability , statistics and statistical physics _ , dordrecht ( reidel , 1983 ) ; e. t. jaynes : _ probability theory : the logic of science _ , cambridge university press ( cambridge , 2003 ) .
tseng , a. caticha , maximum entropy approach to a mean field theory for fluids . in : c. j. williams ( ed ) , bayesian inference and maximum entropy methods in science and engineering , aip conf .
* 659 * , 2003 ( arxiv : cond - mat/0212198 ) .
tseng , a. caticha , maximum entropy approach to the theory of simple fluids , in : g. erickson and y. zhai ( eds ) , bayesian inference and maximum entropy methods in science and engineering , aip conf .
707 * , 2004 ( arxiv : cond - mat/0310746 ) .
caticha a , maximum entropy , fluctuations and priors , in : a mohammad - djafari ( ed ) , bayesian inference and maximum entropy methods in science and engineering , aip conf
. proc . * 568 * , 2001 ( arxiv : math - ph/0008017 ) .                                 as a function of hard - sphere diameter @xmath101 for argon at a density of @xmath261 for different temperatures .
the best @xmath102 is that which minimizes @xmath137 .
( b ) : @xmath137 as a function of @xmath102 for argon at @xmath219 @xmath220 for different densities.,width=480,height=336 ]     for argon for several temperatures at density @xmath261 for @xmath228 .
( b ) : @xmath262 for various @xmath238 at @xmath219 @xmath220 and @xmath263 . by setting @xmath264 ( dotted lines )
we see that the effect of the @xmath265 factor is to cause a slight shift of the distribution.,width=480,height=336 ]        ; ( b ) verlet s molecular dynamics simulation ; and ( c ) the improved eme analysis , for argon at ( a ) : density @xmath263 , temperature @xmath219 @xmath267 , and effective particle number @xmath268 .
( b ) : @xmath263 , @xmath269 @xmath267 , and @xmath270 .
( c ) : @xmath261 , @xmath271 @xmath267 , and @xmath272 .
( d ) : @xmath273 , @xmath248 @xmath267 , and @xmath274.,width=192,height=192 ]      ; ( b ) verlet s molecular dynamics simulation ; and ( c ) the improved eme analysis , for argon at ( a ) : density @xmath263 , temperature @xmath219 @xmath267 , and effective particle number @xmath268 .
( b ) : @xmath263 , @xmath269 @xmath267 , and @xmath270 .
( c ) : @xmath261 , @xmath271 @xmath267 , and @xmath272 .
( d ) : @xmath273 , @xmath248 @xmath267 , and @xmath274.,title="fig:",width=192,height=192 ] +      ; ( b ) verlet s molecular dynamics simulation ; and ( c ) the improved eme analysis , for argon at ( a ) : density @xmath263 , temperature @xmath219 @xmath267 , and effective particle number @xmath268 .
( b ) : @xmath263 , @xmath269 @xmath267 , and @xmath270 .
( c ) : @xmath261 , @xmath271 @xmath267 , and @xmath272 .
( d ) : @xmath273 , @xmath248 @xmath267 , and @xmath274.,width=192,height=192 ]      ; ( b ) verlet s molecular dynamics simulation ; and ( c ) the improved eme analysis , for argon at ( a ) : density @xmath263 , temperature @xmath219 @xmath267 , and effective particle number @xmath268 .
( b ) : @xmath263 , @xmath269 @xmath267 , and @xmath270 .
( c ) : @xmath261 , @xmath271 @xmath267 , and @xmath272 .
( d ) : @xmath273 , @xmath248 @xmath267 , and @xmath274.,width=192,height=192 ]     @xmath267 . also shown
are levelt s experimental results .
( b ) : @xmath275 versus the reduced density @xmath276 calculated using the eme method ( solid line ) and compared to the hansen - verlet simulation for three different isotherms .
the graph shows the appearance of the liquid - gas van der waals loop as the temperature drops .
, width=480,height=336 ]


more than 20 years ago , hong , ou , and mandel showed experimentally that mutually strongly quantum correlated ( entangled ) photon pairs can be emitted in the nonlinear process of parametric down - conversion @xcite at the single - photon level .
for the occurrence of entangled photon pairs , the spontaneous character of the process is important .
entanglement of two photons comprising a photon pair might occur for various physical quantities like frequencies , emission angles ( wave - vectors ) , or polarizations .
perhaps most interestingly , entangled photon pairs manifest themselves in time domain , where both photons are detected within a relatively narrow time window .
this is a direct consequence of the ` point ' character of the emission of a photon pair in the time domain .
the width of the time window characterized by an entanglement time is typically on the order of several hundreds of fs , and is experimentally observable using hong - ou - mandel interferometer .    in the time that has intervened since the original predictions , entangled photon pairs
have been used for numerous experiments demonstrating both fundamental aspects of their physical properties @xcite and their potential for applications .
these also include tests of bell inequalities @xcite , quantum teleportation @xcite , the generation of greenberger - horne - zeilinger states @xcite or quantum computation @xcite that can exploit photon pairs . quantum cryptography with photon pairs @xcite
may be considered the most important application .
we mention that applications to metrology have also been suggested @xcite .
most researchers in the field have focused their attention for the most part on bulk nonlinear crystals , pumped by intense laser beams in type i and type ii configurations .
bright sources of polarization - entangled photon pairs have been fabricated using two nonlinear type i crystals , mutually rotated by 90 degrees @xcite .
the nonlinear processes in an optical cavity have also been used to enhance the photon - pair generation rate @xcite .
in addition , periodically - poled materials may be used to increase the photon - pair generation rate in materials where phase matching is not naturally available @xcite . at present ,
even sources of photon pairs pumped by laser diodes have been developed @xcite .
on the fabrication front , techniques of structures composed of nonlinear thin layers ( of width of several tens or hundreds of nm ) have developed to the point that useful structures with pre - defined properties may be achieved rather easily @xcite .
these nonlinear photonic - band - gap structures @xcite are very promising as sources of photon pairs , as has recently been shown in @xcite . despite the small amount of nonlinear material embedded inside them
, they can generate photon pairs with relatively high efficiency thanks to the constructive interference that involves both the pump and the down - converted fields , their spatial inhomogeneity nothwithstanding .
the enhancement of the photon - pair generation rate has been predicted to be several hundreds and even thousands of times larger than photon - pair generation rates in nonlinear , bulk materials @xcite .
moreover , the spectral and spatial characteristics of the down - converted fields depend on details of the structure , a fact that might be used to control the process , at least to some extent . for example
, photon pairs with very narrow spectra may be obtained from suitable structures .
we note that along the same vein , four - wave mixing in photonic - band - gap nonlinear fibers is also promising as a modern source of photon pairs @xcite .    in this paper , we present a quantum model of photon - pair emission in a nonlinear , one - dimensional photonic - band - gap structure , based upon a perturbative solution of the schrdinger equation .
this model extends those developed for bulk nonlinear materials in @xcite .
it has been shown in @xcite that this approach is compatible with models based on methods of classical nonlinear optics ( multiple - scale spatial and temporal expansion methods ) with a specific kind of stochastic averaging concerning emitted spectra . the quantum model , however , provides a complete description of the nonlinear process .    the paper is organized as follows . in sec
ii , the model is presented in three steps .
description of a quantum field in a layered medium in subsec .
iia is followed by the description of nonlinear quantum interactions in subsec .
iib . in subsec .
iic measurable characteristics of the generated down - converted fields are determined .
the behavior of physical quantities characterizing a photon pair are discussed in sec .
iii , both for cw and pulsed pumping regimes .
iv contains our conclusions .
we consider a stack of @xmath0 nonlinear layers in which spontaneous parametric down - conversion may occur . a sketch of the system is shown in fig .
the @xmath1th layer begins at @xmath2 and ends at @xmath3 , its length is denoted as @xmath4 , @xmath5 .
linear indices of refraction of pump , signal , and idler fields in the @xmath1th layer are denoted as @xmath6 , @xmath7 , and @xmath8 , respectively .
symbols @xmath9 , @xmath10 , and @xmath11 ( @xmath12 , @xmath13 , and @xmath14 ) mean indices of refraction in front of ( beyond ) the sample .
the symbol @xmath15 is used for the nonlinear tensor of @xmath1th layer ; symbols @xmath16 , @xmath17 , and @xmath18 stand for pump- , signal- , and idler- field wave - vectors in the @xmath1th layer .
the structure is pumped by a classical strong pump field that propagates under the angle @xmath19 with respect to @xmath20 axis ( see fig .  [ fig2 ] ) .
its wave - vector @xmath21 lyes in the @xmath22 plane , whereas its electric - field amplitude @xmath23 is perpendicular to the wave - vector @xmath21 .
the pump - field positive - frequency electric - field amplitude at the input of the structure at @xmath24 is denoted as @xmath25 and can be conveniently described using a positive - frequency electric - field amplitude spectrum @xmath26 determined as follows : @xmath27 if we assume a gaussian time profile and a linear polarization of the incident electric - field amplitude @xmath28 in the direction rotated by an angle @xmath29 with respect to an incident te - wave polarization direction @xmath30 ( @xmath31 denotes an incident tm - wave polarization direction , see fig .
[ fig2 ] ) we have : @xmath32 \nonumber \\   & & \mbox { } \times
\xi_p \exp \left ( - \frac{1+ia_p}{\tau_p^2}t^2 \right )   \exp(-i\omega_p^0 t ) , \nonumber \\   & = & { \bf e}_{p_f , te}^{(+)}(z_0,t ) + { \bf   e}_{p_f , tm}^{(+)}(z_0,t ) , \label{2}\end{aligned}\ ] ] where @xmath33 is the pump - pulse amplitude , @xmath34 pulse duration , @xmath35 central frequency , and @xmath36 denotes a chirp parameter of the pulse .
the spectrum @xmath26 determined by eq .
( [ 1 ] ) is given as : @xmath37 \nonumber \\   & & \mbox { } \hspace{-2.5 cm } \times \xi_p \frac{\tau_p}{\sqrt{2(1+ia_p ) } }    \exp \left ( - \frac{\tau_p^2}{4(1+ia_p ) } ( \omega_p-\omega_p^0)^2 \right )    \nonumber \\   & = & { \bf e}_{p_f , te}^{(+)}(z_0,\omega_p ) + { \bf   e}_{p_f , tm}^{(+)}(z_0,\omega_p).\end{aligned}\ ] ] on the other hand , the following spectrum corresponds to cw pumping : @xmath38 \nonumber \\   & & \mbox { } \times \xi_p
\delta(\omega_p-\omega_p^0 ) ; \label{4}\end{aligned}\ ] ] @xmath39 means a dirac delta function .    the pump field incident on the structure is scattered at each boundary inside the structure , so as to achieve a certain profile along the @xmath20 axis , provided that the nonlinear interaction does not lead to pump - field depletion . scattering of the pump field is conveniently described using its decomposition into monochromatic waves .
the positive - frequency electric - field amplitude @xmath40 of a monochromatic component at frequency @xmath41 with polarization @xmath42 ( @xmath43 te , tm ) can be written as follows @xcite : @xmath44     \nonumber \\    & & \mbox { } \hspace{-2cm}+ \sum_{l=1}^{n } { \rm rect}_{z_{l-1},z_l}(z ) \nonumber \\    & & \mbox { } \hspace{-2 cm } \times \left [     a_{p_f,\alpha}^{(l)}(\omega_p){\bf e}_{p_f,\alpha}^{(l)}(\omega_p )     \exp(i{\bf k}_{p , z}^{(l)}(z - z_{l-1 } ) ) \right .
\nonumber \\    & & \mbox { } \hspace{-2cm}\left .
+ a_{p_b,\alpha}^{(l)}(\omega_p )     { \bf e}_{p_b,\alpha}^{(l)}(\omega_p )     \exp(-i{\bf k}_{p , z}^{(l)}(z - z_{l-1 } ) ) \right ]     \nonumber \\    & & \mbox { } \hspace{-2cm}+ { \rm rect}_{z_n,\infty}(z ) \nonumber \\    & & \mbox { } \hspace{-2cm}\times \left [     a_{p_f,\alpha}^{(n+1)}(\omega_p ) { \bf e}_{p_f,\alpha}^{(n+1)}(\omega_p )     \exp(i{\bf k}_{p , z}^{(n+1)}(z - z_n ) ) \right .
\nonumber \\    & & \mbox { } \hspace{-2cm}\left .
+    a_{p_b,\alpha}^{(n+1)}(\omega_p )     { \bf e}_{p_b,\alpha}^{(n+1)}(\omega_p )     \exp(-i{\bf k}_{p , z}^{(n+1)}(z - z_n ) ) \right ] ;     \nonumber \\    & & \alpha = te , tm ; \label{5}\end{aligned}\ ] ] the function @xmath45 equals one for @xmath46 and is zero otherwise .
polarization vectors of @xmath47-waves in @xmath1th layer are denoted as @xmath48 and @xmath49 for forward- and backward - propagating fields with respect to @xmath20 axis , respectively , and they are frequency dependent . the symbol @xmath50 denotes a @xmath20 component of the pump - field wave - vector @xmath51 in @xmath1th layer and is determined by the expression @xmath52 where the wave - vector @xmath51 propagates in the @xmath53th layer under the angle @xmath54 with respect to @xmath55 axis .
the angles @xmath54 fulfill the snell law at the boundaries , i.e. @xmath56 @xmath57 .    the symbol @xmath58 is identified as an electric - field amplitude of the pump field at frequency @xmath59 with polarization @xmath42 incident on the structure from the left - hand side , whereas the symbol @xmath60 describes an electric - field amplitude of the pump field at frequency @xmath41 with polarization @xmath42 incident from the right - hand side .
the remaining amplitudes @xmath61 and @xmath62 are determined using relations at the boundaries and free - field propagation inside the layers : @xmath63 we note that the coefficients @xmath64 and @xmath65 for @xmath66 describe the corresponding electric - field amplitudes at the beginning of @xmath53th layer .    assuming te and tm waves , the boundary transfer matrices @xmath67 and @xmath68 have the form : @xmath69 @xmath70 and @xmath71 .
the free - field propagation matrices @xmath72 can be written as : @xmath73    the positive - frequency electric - field operators @xmath74 and @xmath75 for the signal and idler fields can be decomposed into te- and tm - wave contributions @xmath76 ,
@xmath77 , @xmath78 , and @xmath79 and expressed as follows @xcite : @xmath80    \exp(-i\omega_m t )   \nonumber \\   & & \hspace{-17 mm } = \hat{\bf e}_{m , te}^{(+)}(z , t ) + \hat{\bf   e}_{m , tm}^{(+)}(z , t )    \nonumber \\   & & \hspace{-17 mm } = \frac{1}{\sqrt{2\pi } } \int_{0}^{\infty }    d\omega_m \hat{\bf e}_m^{(+)}(z,\omega_m ) \nonumber \\   & & \hspace{-17 mm } = \frac{1}{\sqrt{2\pi } } \int_{0}^{\infty }    d\omega_m \left [ \hat{\bf e}_{m , te}^{(+)}(z,\omega_m ) +    \hat{\bf e}_{m , tm}^{(+)}(z,\omega_m ) \right ] ,    \nonumber \\   & & m = s , i .
\label{11}\end{aligned}\ ] ] the permittivity of vacuum is denoted as @xmath81 ; @xmath82 means speed of light in vacuum , @xmath83 is the reduced planck constant , and @xmath84 denotes the area of the transverse profile of a beam .
the symbols @xmath85 and @xmath86 mean polarization vectors of mode @xmath87 of @xmath42 wave propagating forward and backward with respect to the @xmath20 axis .
annihilation operators of signal [ @xmath88 , @xmath89 and idler [ @xmath90 , @xmath91 photons in @xmath42 wave introduced in eq .
( [ 11 ] ) can be expressed in the same way as the pump - field amplitude @xmath92 : @xmath93    \nonumber \\   & & \hspace{-8 cm } \mbox { } + \sum_{l=1}^{n } { \rm rect}_{z_{l-1},z_l}(z ) \nonumber \\   & & \hspace{-8 cm } \mbox { } \times \left [    \hat{a}_{m_f,\alpha}^{(l)}(\omega_m ) { \bf e}_{m_f,\alpha}^{(l)}(\omega_m )    \exp(i{\bf k}_{m , z}^{(l)}(z - z_{l-1 } ) ) \right .
\nonumber \\   & & \hspace{-8 cm } \mbox { }   \left .
+ \hat{a}_{m_b,\alpha}^{(l)}(\omega_m )    { \bf e}_{m_b,\alpha}^{(l)}(\omega_m ) \exp(-i{\bf k}_{m , z}^{(l)}(z - z_{l-1 } ) ) \right ]    \nonumber \\   & & \hspace{-8cm}\mbox { } + { \rm rect}_{z_n,\infty}(z ) \nonumber \\   & & \hspace{-8 cm } \mbox { } \times \left [    \hat{a}_{m_f,\alpha}^{(n+1)}(\omega_m ) { \bf e}_{m_f,\alpha}^{(n+1)}(\omega_m )    \exp(i{\bf k}_{m , z}^{(n+1)}(z - z_n ) ) \right .
\nonumber \\   & & \hspace{-8cm}\mbox { } \left .
+ \hat{a}_{m_b,\alpha}^{(n+1)}(\omega_m )    { \bf e}_{m_b,\alpha}^{(n+1)}(\omega_m ) \exp(-i{\bf k}_{m , z}^{(n+1)}(z - z_n ) )    \right],\nonumber \\   & & \hspace{-6 cm } m = s , i ; \;\ ; \alpha = te , tm .
\label{12}\end{aligned}\ ] ] the polarization vectors @xmath94 and @xmath95 give the polarization directions of mode @xmath96 with @xmath42 wave in @xmath1th layer propagating forward and backward , respectively , whereas @xmath97 is a @xmath20 component of the wave vector @xmath98 of this mode : @xmath99 the angle @xmath100 characterizes the direction of propagation of mode @xmath87 with respect to the @xmath20 axis .
the angles @xmath101 are given by the snell law at the boundaries , i.e. @xmath102 @xmath103 , where @xmath104 stands for the angle of incidence of mode @xmath87 .
the operators @xmath105 and @xmath106 obey the following commutation relations : @xmath107     & = & \delta_{\alpha,\alpha ' } \delta_{m , m ' } \delta_{l , l ' } \delta(\omega_m-\omega'_m ) ,     \nonumber \\   { } [ \hat{a}_{m_f,\alpha}^{(l)}(\omega_m ) ,     \hat{a}_{m'_f,\alpha'}^{(l')}(\omega'_m ) ] & = & 0 ,     \nonumber \\   { } [ \hat{a}_{m_b,\alpha}^{(l)}(\omega_m ) ,     \hat{a}_{m'_b,\alpha'}^{(l')\dagger}(\omega'_m ) ]     & = & \delta_{\alpha,\alpha'}\delta_{m , m ' } \delta_{l , l ' } \delta(\omega_m-\omega'_m ) ,     \nonumber \\   { }
[ \hat{a}_{m_b,\alpha}^{(l)}(\omega_m ) ,     \hat{a}_{m'_b,\alpha'}^{(l')}(\omega'_m ) ] & = & 0 ,     \nonumber \\   { } [ \hat{a}_{m_f,\alpha}^{(l)}(\omega_m ) ,     \hat{a}_{m'_b,\alpha'}^{(l')}(\omega'_m ) ] & = & 0 ,     \nonumber \\   { } [ \hat{a}_{m_f,\alpha}^{(l)}(\omega_m ) ,     \hat{a}_{m'_b,\alpha'}^{(l')\dagger}(\omega'_m ) ] & = & 0 ;     \nonumber \\     & & \hspace{-4.5 cm } m , m'=s , i ; \;\ ; \alpha = te , tm ; \ ; \ ; l , l'=0,\ldots , n+1 .\end{aligned}\ ] ]    the photonic - band - gap structure imposes the following relations among the operators @xmath105 and @xmath106 acting in @xmath1th layer : @xmath108 transfer matrices @xmath109 at the boundaries and free - field propagation matrices @xmath110 for @xmath111 and @xmath112 are defined in the same way as those given in eqs .
( [ 9 ] ) and ( [ 10 ] ) for the pump - field amplitudes .
the hamiltonian @xmath113 describing spontaneous parametric down - conversion in a nonlinear medium of volume @xmath114 at time @xmath115 can be written as : @xmath116 , \label{17}\end{aligned}\ ] ] where @xmath117 denotes a third - order tensor of nonlinear coefficients , the symbol @xmath118 is shorthand of the tensor @xmath119 with respect to its three indices , and @xmath120 stands for a hermitian conjugated term .
the negative - frequency electric - field operators @xmath121 for @xmath122 have been introduced in eq .
( [ 17 ] ) ( @xmath123 ) . decomposing the interacting fields into te and tm waves , and using the inverse fourier transformation in eq .
( [ 17 ] ) , we arrive at : @xmath124 .
\label{18}\end{aligned}\ ] ] integrations over the variables @xmath125 and @xmath126 in eq .
( [ 17 ] ) impose conditions for @xmath125 and @xmath126 component of wave - vectors ; @xmath127 [ @xmath128 is assumed ] and @xmath129 .
the latter @xmath39-function provides the following relation between the angles @xmath130 and @xmath131 of modes @xmath132 and @xmath133 in the @xmath1th layer : @xmath134 ,    \nonumber \\   & & l=0,\ldots , n+1 . \label{19}\end{aligned}\ ] ]    solution of the schrdinger equation to first order in nonlinear perturbation together with the assumption of incident vacuum state @xmath135 in signal and idler fields provides the output state @xmath136 of signal field with @xmath137 polarization and idler field with @xmath138 polarization in the form : @xmath139    \nonumber \\   & & \hspace{-1cm}\mbox { } \times    { \rm sinc}\left[\frac{1}{2}(k_m^{(l ) } - k_n^{(l ) } - k_o^{(l ) } ) l_l \right ]
\delta(\omega_p-\omega_s-\omega_i )    \nonumber \\   & & \hspace{-1 cm } \mbox { } \times    \hat{a}_{n,\beta}^{(l)\dagger}(\omega_s ) \hat{a}_{o,\gamma}^{(l)\dagger}(\omega_i )    |{\rm vac } \rangle .
\label{20}\end{aligned}\ ] ] the wave - vectors @xmath140 introduced in eq .
( [ 20 ] ) are defined as @xmath141 and @xmath142 for @xmath143 .
the operators @xmath144 and @xmath145 for @xmath42 waves in mode @xmath87 in @xmath1th layer can be expressed in terms of the operators @xmath146 and @xmath147 .
these relations can be , e.g. , written in the form : @xmath148 \nonumber \\   & & \hspace{-3cm}\mbox { }   \times \pmatrix { 1/({\cal s}_{m,\alpha})_{11}(\omega_m ) &    - ( { \cal s}_{m,\alpha})_{12}(\omega_m ) /    ( { \cal s}_{m,\alpha})_{11}(\omega_m ) \cr 0 & 1 } \nonumber \\   & & \hspace{-3 cm } \mbox { } \times    \pmatrix{\hat{a}_{m_f,\alpha}^{(n+1)}(\omega_m ) \cr   \hat{a}_{m_b,\alpha}^{(0)}(\omega_m ) } ,    \nonumber \\   & & \hspace{-2 cm } m = s , i ; \;\ ; \alpha = te , tm ; \;\ ; l=1,\ldots , n .
\label{21}\end{aligned}\ ] ] the matrix @xmath149 used in eq .
( [ 21 ] ) describes the propagation of @xmath42 wave in field @xmath87 through the whole structure : @xmath150 , \nonumber \\    & &   m = s , i ; \;\ ; \alpha = te , tm .
\label{22}\end{aligned}\ ] ]    similarly , the pump - field amplitudes @xmath64 and @xmath151 for @xmath42 waves in @xmath1th layer can be determined from the amplitudes of the incident fields @xmath152 and @xmath153 as follows : @xmath154 \nonumber \\   & & \hspace{-3cm}\mbox { }   \times \pmatrix { 1 & 0 \cr   - ( { \cal s}_{p,\alpha})_{21}(\omega_p ) /({\cal s}_{p,\alpha})_{22}(\omega_p ) &   1/({\cal s}_{p,\alpha})_{22}(\omega_p ) } \nonumber \\   & & \hspace{-3 cm } \mbox { } \times    \pmatrix{a_{p_f,\alpha}^{(0)}(\omega_p ) \cr a_{p_b,\alpha}^{(n+1)}(\omega_p ) } ,    \nonumber \\   & & \hspace{-.5 cm } \alpha = te , tm ; \;\ ;    l=1,\ldots , n .
\label{23}\end{aligned}\ ] ] the matrix @xmath155 describes the propagation of a classical pump field with polarization @xmath42 through the whole structure , i.e. @xmath156 .
\label{24}\ ] ]    we note that the expression in eq .
( [ 20 ] ) for the output state @xmath157 including relations written in eqs .
( [ 21]-[24 ] ) can be formally recast into a compact form using the so - called left - to - right ( @xmath158 ) and right - to - left ( @xmath159 ) modes introduced in the classical electro - magnetic theory of layered structures ( for details , see @xcite ) .
fields exiting the structure at @xmath160 are described by @xmath161 functions , whereas @xmath162 functions are appropriate for fields exiting at @xmath163 . in classical theory
, this corresponds to the picture in which every nonlinear layer is a source ( emitting dipole ) of photon pairs @xcite and the fields leaving the structure are given as a sum of contributions from all layers .
we assume that the outgoing signal ( idler ) field is detected using an analyzer with the polarization that forms an angle @xmath164 ( @xmath165 ) with respect to the te - wave polarization direction ( see fig .
[ fig2 ] ) . in order to get the right wave - function
, we transform the operators of the outgoing signal ( idler ) fields into the basis with the polarization vectors given by angles @xmath166 ( @xmath165 ) and @xmath167 ( @xmath168 ) using the following formulas : @xmath169    substituting relations in eqs . ( [ 21 ] ) , ( [ 23 ] ) , and ( [ 25 ] ) into the expression @xmath170 determined using eq .
( [ 20 ] ) , we arrive at the output state @xmath171 describing a signal photon polarized along the angle @xmath166 and an idler photon polarized along the angle @xmath165 : @xmath172 . \label{26}\end{aligned}\ ] ] the function @xmath173 introduced in eq .
( [ 26 ] ) has the meaning of probability amplitude of having a signal photon in field @xmath87 polarized along the angle @xmath174 and its entangled idler twin in field @xmath175 polarized along the angle @xmath165 at the output from the photonic - band - gap structure .
we are interested only in the part @xmath176 of the output state @xmath177 given in eq .
( [ 26 ] ) that describes the generated photon - pairs . including the time evolution of the free - fields outside the photonic - band - gap structure , we can write : @xmath178 where @xmath179 the creation operators @xmath180 , @xmath181 , introduced in eq .
( [ 28 ] ) describe the linearly polarized photons outside the photonic - band - gap structure with polarization angles @xmath166 and @xmath165 , and are given as : @xmath182 we note that we do not explicitly express the dependence of quantities on the polarization angles @xmath166 and @xmath165 below .
the mean number @xmath183 of photon pairs that have a signal photon in the frequency interval @xmath184 around frequency @xmath185 and its twin idler photon in frequency interval @xmath186 around frequency @xmath187 in field @xmath188 is given by @xmath189 where the operator @xmath190 , the density of photons , is defined as @xmath191 using eq .
( [ 28 ] ) , the expression for @xmath183 in eq .
( [ 30 ] ) can be written as : @xmath192    we get the following expressions for , e.g. , the mean number @xmath193 of signal photons in frequency interval @xmath194 around frequency @xmath185 : @xmath195    the overall number @xmath196 of photon pairs emitted into field @xmath197 is determined by the expression @xmath198    the signal - field energy spectrum @xmath199 of field @xmath197 can be easily determined using the expression @xmath200    the energy spectra @xmath201 ( @xmath202 ; @xmath203 ) characterizing the outgoing down - converted fields without an inclusion of pair entanglement can be evaluated using the signal - field spectra @xmath204 determined in eq .
( [ 35 ] ) ; the idler - field spectra @xmath205 are determined analogously ; @xmath206    the properties of the down - converted fields in the time domain can be conveniently described using a two - photon amplitude @xmath207 giving the probability amplitude of detecting a signal photon at time @xmath208 and an idler photon at time @xmath209 : @xmath210 assuming the state @xmath211 given in eq .
( [ 28 ] ) , the expression in eq .
( [ 37 ] ) for the two - photon amplitude @xmath212 can be rearranged into the form : @xmath213 where the ` fourier transform ' @xmath214 of the function @xmath173 has been introduced : @xmath215    the photon flux of , e.g. , the signal photons , @xmath216 at time @xmath208 over the transverse profile of area @xmath84 in field @xmath188 is defined as @xmath217 the photon flux @xmath218 can be determined using the function @xmath173 ; @xmath219 if the spectrum of the idler field is narrow , we may use the alternative expression @xmath220    entanglement of the signal and idler photons in a pair can be detected in a hong - ou - mandel interferometer ( see fig .  [ fig3 ] ) .    in order to achieve this type of interference between the signal and idler photons , we rotate polarizations of both photons such that they are the same , and subsequently introduce a relative time delay @xmath221 between the two photons
. then two photons arrive at a 50/50  % beam - splitter whose output ports are monitored by detectors . the measured coincidence - count rate @xmath222 is given by the number of simultaneously detected photons at detectors @xmath223 and @xmath224 , placed at the output ports of the beam - splitter in a given time interval .
there occurs quantum interference between two paths leading to a coincidence count ; either a signal photon is detected by detector @xmath225 and its twin idler photon by detector @xmath226 or vice versa .
the normalized coincidence - count rate @xmath227 in this interferometer can be expressed as follows : @xmath228 where @xmath229 \label{44}\end{aligned}\ ] ] and @xmath230 the symbol @xmath231 stands for real part , and the two - photon amplitude @xmath212 has been introduced in eq .
( [ 37 ] ) . using the expressions in eqs .
( [ 44 ] ) and ( [ 45 ] ) , we arrive at : @xmath232 ,    \nonumber \\    & &   \\
r_0^{mn } & = & \frac{\hbar^2}{8\epsilon_0 ^ 2 c^2    { \cal b}^2 } \int_{0}^{\infty } d\omega_s \ , \int_{0}^{\infty } d\omega_i \ ,    \ ; \omega_s \omega_i   |\phi^{mn}(\omega_s,\omega_i)|^2 .
\nonumber \\   & & \end{aligned}\ ] ]      if cw pumping is considered , the spectrum is described using eq .
( [ 4 ] ) . as a consequence
the following expressions are obtained : @xmath233 in this case the second power of the dirac delta function @xmath234 is formal , so that the above expression should be replaced by : @xmath235 where the period of nonlinear interaction goes from @xmath236 to @xmath237 .
expressions for the physical quantities determined above have to be normalized by @xmath238 , which indicates that these quantities are related to the period of @xmath239 .
as an example , we study the properties of a nonlinear photonic - band - gap structure composed of 25 nonlinear layers of gan of thickness of 117  nm among which there are 24 linear layers of aln with a thickness of 180  nm . material characteristics of gan and aln can be found , e.g. , in @xcite .
this structure is resonant for the pump field at wavelength of 664.5  nm , designed to correspond to the first resonance peak near the band edge .
this situation favors the efficient generation of photon pairs at roughly double the pump wavelength , at angles that correspond to the first , second , and third resonance peaks , respectively , for the down - converted fields .
we assume that the gan crystallographic @xmath125 axis coincides with the @xmath20 axis of propagation ( see fig .
[ fig1 ] ) . in this configuration
there is strong nonlinear interaction between te components of the pump , signal , and idler beams .
the nonlinear coefficient @xmath240 for this geometry is chosen to be 10 pmv@xmath241 .
we first consider pumping by a cw laser tuned at first resonance peak , i.e. @xmath242  nm . a photon pair whose signal photon is emitted along the angle @xmath243 may be described by a two - photon amplitude @xmath244 with typical shapes both in time and spectral domains . in the spectral domain , the two - photon amplitude @xmath245 can be written in the form @xmath246 , where the function @xmath247 is linearly proportional to the spectral amplitude of the signal field .
this is a consequence of stationarity of the emitted down - converted fields .
the probability @xmath248 of detecting a signal photon at time @xmath249 and an idler photon at time @xmath250 is shown in fig .
[ fig4 ] for mode @xmath251 ( both photons exit the structure at @xmath160 ) .
the shape of the probability @xmath252 can be easily understood when we consider the fact that both photons in a pair are born together at the same instant of time , and then propagate independently inside the structure .
this independent zig - zag propagation increases the average time delay between the two photons , which takes into account all possible realizations of the random process of propagation of a photon that undergoes multiple bounces inside the structure .
the random zig - zag propagation causes interference between different paths .
this leads to a typical structure with local minima and maxima along the direction @xmath253 ; as shown in fig .
[ fig4]b . also , the greater the difference @xmath254 the smaller the value of the two - photon amplitude @xmath255 because the propagation with more zig - zags is less probable .
we note that typical stripes along the direction @xmath256 reflect stationarity of the process .
typical spectral properties of this structure are documented in fig .
[ fig5 ] , where the signal - field energy spectrum @xmath257 is depicted as a function of the angle @xmath258 for field @xmath251 .
the spectrum @xmath259 giving the probability of emission of a signal photon at frequency @xmath185 along the angle @xmath260 is a complex function in its variables because it is built up by interference of fields coming from the multilayer structure .
a maximum photon - pair generation rate is obtained for @xmath261  deg for degenerate frequencies ( @xmath262 ) .
photons from a generated photon pair can exit the structure also at @xmath24 ; their energy spectra @xmath263 , @xmath264 and @xmath265 are similar to that shown in fig .
[ fig5 ] . this is because efficient generation of photon pairs occurs provided both the signal and idler fields are resonant in the structure so that the generated photons exit the structure at @xmath24 or @xmath160 with comparable probabilities .
the signal - field intensity transmission @xmath266 as a function of normalized frequency @xmath267 and angle @xmath243 is shown in fig .
[ fig6 ] for comparison .    the angle - dependence of the signal - field fwhm ( full - width at half maximum ) @xmath268 of the energy spectrum in wavelengths and value of energy spectrum @xmath269 at the central frequency @xmath270
are shown in fig .
high photon - pair generation rates are observed for the frequency - degenerate case ( @xmath271 ) for values of the angle @xmath258 approximately equal to 14 , 35 , and 55 degrees . for these values of the angle @xmath258 the peaks of the signal - field intensity transmission @xmath266 cross the value @xmath185 equal to @xmath272 ( see fig .
[ fig6 ] ) .
the first and highest peaks in the signal - field energy spectra @xmath273 ( see fig .
[ fig7]b ) correspond to the case when both the signal and idler photons are tuned at the first resonance near the band edge . at this frequency ,
localization of the down - converted fields is maximum , and the nonlinear process is strongest .
the other two sets of peaks in the signal - field energy spectra @xmath273 correspond to the second and third peaks in the intensity transmission @xmath266 . at these frequencies ,
the nonlinear interaction is somewhat weaker due to smaller field intensities , i.e. , slightly worse localization of the down - converted fields .
the spectral fwhms for these frequency - degenerate emissions lye in the interval from 10  nm to 15  nm ( see fig .
[ fig7]a ) . for these values of the angle @xmath243
the best available constructive interference in the structure occurs . in a relatively narrow frequency interval the signal and idler fields
are enhanced by constructive interference . outside this frequency interval , destructive interference occurs . for the other values of the angle
@xmath258 constructive interference is less pronounced , the range of frequencies having constructive interference is broader , and spectral shapes with several peaks may occur ( see fig .
[ fig5]b ) .
we thus obtain down - converted fields with spectral fwhms reaching even 80 nm ( see fig .
[ fig7]a ) . in this case
the fwhm indicates the entire active width , which might include several peaks . however , photon - pair generation rates are low in this case .    quantum correlations ( entanglement ) between the signal and idler photons in a photon pair are visible in the profile of the normalized coincidence - count rate @xmath274 in a hong - ou - mandel interferometer ( see fig .
[ fig8 ] for mode @xmath275 ) .    for values of the angle @xmath258 with strong constructive interference ( @xmath276  deg ) this profile
is given by a relatively broad dip ( see fig .  [
fig8]b ) . on the other hand ,
oscillating curves with several local minima and maxima characterize profiles of the normalized coincidence - count rates @xmath274 for other values of the angle @xmath258 . a typical profile for this case is shown in fig .  [ fig8]b .
these profiles reflect broader spectra of the down - converted fields along these angles as well as a nonzero difference of the central frequencies @xmath270 and @xmath277 of the signal and idler fields .
period of these oscillations is proportional to @xmath278 .
position @xmath279 of the central dip , its fwhm @xmath280 , and visibility @xmath281 of the normalized coincidence - count rate @xmath227 as functions of the angle @xmath258 are shown in fig .
[ fig9 ] .
we can see in fig .
[ fig9]a that the maximum overlap of the signal and idler photon fields is reached for a nonzero mutual time delay @xmath221 for nonsymmetric modes @xmath282 and @xmath283 .
this means that the signal and idler photons exit the structure with a small mutual time delay .
this nonzero mutual time delay means that interference effects inside the structure prefer the generation of photon pairs from nonlinear layers positioned closer to one edge of the structure , or in any case approximately co - located . the visibility @xmath281 of the coincidence - count rate @xmath227 equal to one is naturally observed for symmetric modes @xmath275 and @xmath284 .
values of the visibility @xmath281 better than 0.9 can also be reached in nonsymmetric modes @xmath282 and @xmath283 provided that the fields are generated along the angles @xmath258 in the vicinity of those with a strong constructive interference . the structure of quantum correlations between the signal and idler fields in time domain is complex and prevents the visibility @xmath281 from reaching higher values along the angles @xmath258 , where there are broad spectra of the down - converted fields .
the width @xmath285 of the dip in the normalized coincidence - count rate @xmath274 determines the entanglement time , i.e. the time interval in which both photons can be detected .
the fwhm @xmath285 of the central dip for the structure ranges from 50 fs for values of the angle @xmath258 with broad spectra of the down - converted fields , to 150 - 250 fs for values of the angle @xmath243 at where strong constructive interference occurs ( see fig .  [ fig9]c ) .
we now assume that the structure is pumped by an ultrashort pulse with an unchirped gaussian profile [ see eq.([2 ] ) ] having pulse duration @xmath286 equal to 200  fs , and central ( carrier ) frequency @xmath35 at the wavelength of 664.5  nm .
the squared modulus @xmath287 fits well into the peak of resonance of the pump - field intensity transmission @xmath288 , as depicted in fig .
[ fig10 ] .    only signal and idler photons with frequencies @xmath185 and @xmath289 for which the sum @xmath290 lyes within the pump - pulse spectrum may be generated .
the allowed values of the frequencies @xmath185 and @xmath187 are indicated in the graph of the probability @xmath291 of emitting a signal photon at frequency @xmath185 and its twin at frequency @xmath187 . a typical shape of the probability latexmath:[$    @xmath251 .
we note that the shape of the probability @xmath293 in fig .
[ fig11 ] may be considered as being composed of many curves defined above the lines @xmath294 , and it may be assumed to be linearly proportional to those characterizing cw pumping at the frequency @xmath295 .
the probability @xmath252 of detecting a signal photon at time @xmath208 and its twin at time @xmath250 depicted in fig .
[ fig12 ] shows that both signal and idler photons occur in sharp time windows as a consequence of pumping by an ultrashort pulse .
the probability @xmath252 has a typical droplet shape that reflects the fact that the longer the photons from a pair propagate inside the structure , the greater the average difference of the occurrence times @xmath208 and @xmath250 of the photons .
the structure with local minima and maxima farther from the diagonal in fig .
[ fig12 ] reflects interference in the layered structure .    because the pump - pulse amplitude spectrum @xmath296 fits well into the peak of resonance of the structure , and the emitted frequencies of the down - converted fields have also to fit into the structure , the energy spectra of the signal and idler fields with pulsed pumping are very similar to those obtained for cw pumping
also , the behavior of photon pairs generated by pulsed pumping in the hong - ou - mandel interferometer is similar to that observed with cw pumping ; only the visibilities @xmath281 of the normalized coincidence - count rates @xmath274 are slightly worse .
the signal and idler fields are now emitted in the form of ultrashort , pulsed fields in multimode fock states .
a typical dependence of the signal - field photon flux @xmath297 at time @xmath208 and angle @xmath258 is shown in fig .
[ fig13 ] .
[ fig13 ] suggests that the photon flux @xmath298 typically spreads to longer times as a consequence of zig - zag movement , i.e. multiple scattering , of photons inside the structure .
time delay @xmath299 and fwhm @xmath300 of the pulsed signal field as well as photon flux @xmath301 at the center of the pulsed field as they depend on the angle @xmath260 are depicted in fig .
[ fig14 ] .
the time delay @xmath299 of the signal pulse with respect to the pump pulse ranges from 70 fs to 200 fs .
it reaches its maximum values for values of the angle @xmath258 with strong constructive interference ( @xmath276  deg ) .
fwhm @xmath300 of the signal pulse goes from 270  fs to 320  fs .
the narrower the spectrum of the signal field ( given by @xmath302 in wavelengths ) the greater the signal - field pulse duration @xmath300 ( compare figs .  [ fig7]a and [ fig14]b ) .
the largest value of the photon flux @xmath303 can be reached for the angle @xmath258 equal to 13.8  deg ( see fig .  [
fig14]c ) .
these results clearly show that the considered photonic - band - gap structure maintains a pulsed character of the nonlinear process and allows generation of pulsed down - converted fields with an ultrashort duration .
such pulsed fields are required in many experiments with photon pairs in which time synchronization of several photon pairs is necessary .
the model that we have developed is one dimensional with respect to the spatial coordinates , and does not include all aspects ( especially those related to transverse profiles of the interacting optical fields ) of the nonlinear process in a real structure .
for this reason we do not determine the absolute values of photon - pair generation rates .
instead , we judge the efficiency of the suggested structure ( stemming from constructive interference of the interacting optical fields ) with respect to an ideal reference structure which fully exploits the nonlinearity , but does not rely on interference .
this reference structure has formally all linear indices of refraction equal to one , so effects on the boundaries between layers are suppressed .
it is also assumed that the nonlinear process is fully phase matched .
the orientations of nonlinear layers , as well as polarizations of the interacting optical fields , are such that the greatest nonlinear effect occurs . a photon pair emitted from this structure
is described by the following output state @xmath304 [ compare eq .
( [ 20 ] ) ] : @xmath305 where @xmath306 [ @xmath307 denotes a creation operator of a photon in the signal [ idler ] field .
the function @xmath308 used in eq .
( [ 50 ] ) gives the maximum value from the elements of a tensor in its argument .    the relative photon - pair generation rate @xmath309 of a pair with a signal photon at frequency @xmath185 in mode @xmath197 is then given as : @xmath310 where the signal - field energy spectrum @xmath311 is given in eq .
( [ 35 ] ) , and the signal - field energy spectrum @xmath312 of the reference structure is determined in the same way using the output state @xmath304 written in eq .
( [ 50 ] ) .
the relative photon - pair generation rates @xmath313 , @xmath314 , @xmath315 , and @xmath316 for the signal - field frequency @xmath317 as a function of the angle @xmath258 are shown in fig .
[ fig15 ] assuming cw pumping .
the relative photon - pair generation rates for pulsed pumping are lower in comparison with cw pumping , because frequencies at the edges of the pump - field spectrum have lower efficiencies than those in the middle . despite this , the considered pulsed pumping
has the overall relative photon - pair generation rate @xmath318 equal to 9 , whereas @xmath319 for cw pumping for @xmath261  deg .
the structure that we consider is approximately 12 times more efficient than the reference structure , provided we consider all the emitted photon pairs ( given by the overall relative photon - pair generation rate @xmath318 in fig .  [ fig15 ] ) .
even when we restrict ourselves only to the mode @xmath251 , we still have an enhancement of the nonlinear process by a factor of three .
for comparison purposes , one layer of gan of thickness @xmath320  nm ( i.e. , containing the same amount of nonlinear material as our structure ) has the overall relative photon - pair generation rate @xmath318 equal to 0.09 for values of the angle @xmath243 less than 20  deg .
the inclusion of linear layers of aln inside the structure and the accompanying interference effects thus increase the photon - pair generation rates by two orders of magnitude .
we have developed a quantum model of spontaneous parametric down - conversion ( generating entangled photon pairs ) in one - dimensional , nonlinear , layered media . using the model we have determined a two - photon amplitude , and we have provided measurable characteristics of the down - converted fields : marginal signal and idler energy spectra , time - dependent photon fluxes of the signal and idler fields , coincidence - count interference patterns in the hong - ou - mandel interferometer , and photon - pair generation rates .
a specific structure was suggested as an efficient source of photon pairs , and interference effects of the interacting fields can enhance the photon - pair generation rates by as much as two hundred times .
the widths of the spectra of the down - converted fields and entanglement time of photons in a pair depend strongly on the angle of emission , and vary from 10  nm to 80  nm and from 50  fs to 250  fs , respectively .
pumping the structure with an ultrashort pulse of time duration of several hundreds of fs can lead to the generation of pulsed signal and idler fields extending over several hundreds of fs .
hong , z.y .
ou , and l. mandel , phys . rev . lett . * 59 * , 2044 ( 1987 ) .
l. mandel , e. wolf , _ optical coherence and quantum optics _
( cambridge univ . press .
cambridge , 1995 ) .
j. peina , z. hradil , and b. juro , _ quantum optics and fundamentals of physics _
( kluwer , dordrecht , 1994 ) .
d. bouwmeester , j .- w .
pan , k. mattle , m. eibl , h. weinfurter , and a. zeilinger , nature * 390 * , 575 ( 1997 ) .
d. bouwmeester , j .- w .
pan , m. daniell , h. weinfurter , and a. zeilinger , phys .
. lett . * 82 * , 1345 ( 1999 ) .
d. bouwmeester , a. ekert , and a. zeilinger ( eds . _ the physics of quantum information _ ( springer , berlin , 2000 ) .
d. bru , n. ltkenhaus , in _
applicable algebra in engineering , communication and computing _ vol .
10 ( springer , berlin , 2000 ) ; p. 383 .
a. migdall , physics today * 1 * , 41 ( 1999 ) .
kwiat , e. waks , a.g .
white , i. appelbaum , and p.h .
eberhard , phys .
a * 60 * , r773 ( 1999 ) .
kwiat , nature * 412 * , 866 ( 2001 ) .
y. nambu , k. usami , y. tsuda , k. matsumoto , k. nakamura , phys .
rev . a * 66 * , 033816 ( 2002 ) .
shapiro and n.c .
wong , j. opt .
b : quantum semiclass . opt .
* 2 * , l1 ( 2000 ) .
kuklewicz , m. fiorentino , g. messin , f.n.c .
wong , and j. shapiro , opt .
express * 13 * , 127 ( 2005 ) .
p. trojek , c. schmid , m. bourennane , h. weinfurter , and c. kurtsiefer , opt .
express * 12 * , 276 ( 2004 ) .
m. bertolotti , c.m .
bowden , and c. sibilia , _ nanoscale linear and nonlinear optics _ , aip vol .
560 ( aip , melville , 2001 ) .
e. yablonovitch , phys .
58 , 2059 ( 1987 ) .
s. john , phys .
58 , 2486 ( 1987 ) .
joannopoulos , r.d .
meade , and j.n .
winn , _ photonic crystals _
( princeton university press , princeton , 1995 ) .
k. sakoda , _ optical properties of photonic crystals _ ( springer , newyork , 2005 ) , 2nd ed .
vamivakas , b.e.a .
saleh , a.v .
sergienko , and m.c .
teich , phys .
a 70 , 043810 ( 2004 ) .
m. centini , j. peina jr . , l. sciscione , c. sibilia , m. scalora , m.j .
bloemer , and m. bertolotti , phys .
rev . a * 72 * , 033806 ( 2005 ) .
x. li , p.l .
voss , j.e .
sharping , and p. kumar , phys .
lett . * 94 * , 053601 ( 2005 ) .
j. fulconis , o. alibart , w.j .
wadsworth , p. s. russell , and j. g. rarity , opt .
express * 13 * , 7572 ( 2005 ) .
t.e . keller and m.h .
rubin , phys .
a 56 , 1534 ( 1997 ) .
j. peina jr , a.v .
sergienko , b.m .
jost , b.e.a .
saleh , and m.c .
teich , phys .
a * 59 * , 2359 ( 1999 ) .
g. di giuseppe , l. haiberger , f. de martini and a.v .
sergienko , phys .
a * 56 * , r21 ( 1997 ) .
grice , r. erdmann , i.a .
walmsley , and d. branning , phys .
a * 57 * , r2289 ( 1998 ) .
p. yeh , _ optical waves in layered media _ ( wiley , new york , 1988 )
. w. vogel , d.g.welsch , s. walentowicz , _ quantum optics _
( wiley - vch , weinheim , 2001 ) .
g. daguanno , n. mattiucci , m. scalora , m.j .
bloemer , and a.m. zheltikov , phys .
e * 70 * , 016612 ( 2004 ) .
sanford , a.v .
davidov , d.v .
tsvetkov , a.v .
dmitriev , s. keller , u.k .
mishra , s.p .
denbaars , s.s .
park , j.y .
han and r.j .
molnar , j. apl . phys . * 97 * , 053512 ( 2005 ) .
j. miragliotta , d.k .
winckenden , t.j .
kistenmacher , w. a. bryden , j. opt .
b * 10 * , 1447 ( 1993 ) .