Nov 1, 2012 - formulation uses the Cauchy stress and the Green-Lagrange strain without omission of higher order .... equation of motion for an Euler-Bernoulli beam were derived under the con- dition of finite deformation. ... beam problem, we obtain
Abstract. The creation of transuranium nuclides under pulsed neutron fluxes of thermonuclear explosions is investigated within the kinetic model of the astrophysical r-process, taking into account the time dependence of the external parameters and th
Aug 7, 2015 - numbers whose Galois conjugates all lie in the local fields Lp = C for each p â S, ... Lp/Qp will be a given finite normal extension for each p â S. LS will denote the field of all algebraic numbers all of whose Galois con- jugates
The investigations of the quantum oscillation effects of Shubnikovâde Haas  and of de Haasâvan Alphen , as well as a number of high frequency phenomena in strong magnetic fields [4,5] allow for a complete determination of the form of FS. Th
decision-maker can cooperate in the normal operation mode of the robotic system or in the time sharing mode with the ... application fields may be different â from manipulation to autonomous problem solving and necessary ... peculiarities of soluti
May 8, 2008 - such as Î structures can be equally suitable. We begin by preparing each of the optically active qubits, which term bro- kers following Ref. , in state .... a minimum acceptable fidelity of 1 â 10â3, thus there is a 'no-go' par
significant differences in hedging performance between short and long hedgers. Thus, tail ... John Cotter, Director of Centre for Financial Markets, Department of Banking and Finance, University College. Dublin ... Accounting and Finance, Dublin Inst
In this paper I consider learning Bayesian network structures on a finite set of discrete variables, under ... a known distribution, which is equivalent to learning from an infinite data set. Section 5 considers the more ..... edge Management (CIKM'9
Abstract: Extreme conditions in natural flows are examined, starting with a turbulent big bang. A hydro-gravitational-dynamics cosmology model is adopted. Planck-Kerr turbulence instability causes Planck-particle turbulent combustion. Inertial-vortex
Apr 11, 2016 - USA [E-Mail:[email protected]]. 2 Institute for Physical Chemistry ... our DFT results, where extensive variation of the persistence length, the chain length and the polymer concentration in the bulk solution will .... proces
May 3, 2010 - capturing, improvement of biocompatibility of drugs, etc.7,11â15 The rich variety of practical uses in our daily life ..... at least in the range of parameters investigated in this paper. information on the ..... G. de Gennes, Scaling
Oct 24, 2013 - keep making things hotter and hotter, or keep squeezing things harder and harder . These questions ... quark matter , it was recently found that photons can form a Bose-Einstein condensate . In fact, ... (or not) confinemen
Oct 20, 2014 - ever, the separate species co-location suggests the ten- dency for H and C to demix under these conditions. Fig- ure 10 shows the decomposition pathway for the molar mixture (Ï0 = 1.676g/cm3, x = 0.19). The additional xenon increases
Nov 19, 2014 - to visualize and quantify metabolic and physiological processes in the human body. The use of .... way of describing the interaction between these compartments one associates a constant capable to .... Moreover, throughout the whole pa
Dec 20, 2017 - and Computer Engineering and Beckman Institute, Univerisity of Illinois at. Urbana-Champaign, Urbana, IL, 61801 .... methods are moderate in practice, by adding small noise or pixel translations, etc. The learned ..... n4 = 40, c4 = 3;
Dec 28, 2014 - Generalities about hypersurfaces of revolution. Let M be an embed- ded hypersurface of revolution in Rn+1 obtained by rotation of a curve c(u) = (r(u),z(u)) in the plane x1xn+1 around the axis xn+1. Let N =(Ëz(u),âËr(u))/|Ëc(u)| t
and risk faced by a financial trader on a trading floor. The true risk of an accident is hard to recreate for airline pilots or truck drivers. A soldier can only experience the reality of combat on the battlefield. In each of these examples, the work
Oct 2, 2014 - density fluctuations varies less than 10% for molecules that do not form hydrogen bonds. Our findings ... A number of indirect methods have been proposed to estimate dynamic correlation volumes. [15,16,17 ..... Third-order spectrum of P
I. Blaschzyk and I. Steinwart/Classification Rates under Refined Margin Conditions. 8 that means X0 = âXX1, with Hdâ1(X0) > 0 and let X0 be (d â 1)-rectifiable. Let L be the classification loss and let for fixed n â¥ 1 and Ï â¥ 1 the bounds
Apr 18, 2005 - 2UVSOR Facility, Institute for Molecular Science, Okazaki 444-8585, Japan. 3Research Reactor Institute, Kyoto University, Osaka 590-0494, Japan. 4Department of ..... 2004A0231-NSa-np and 2004B0545-NSa-np).  J. Rossat-Mignod et al.,
remains legal in the Minkowski pattern. The results of the last section were ...... where matrices. ( ). ,. Ë. E H. S r are composed from the Mie scattering operators and complex-angle spatial rotation operators . Eqs. (39) imply that the frame
May 1, 2017 - The matrices will be bi-linear objects acting on a space M â V . In coordinate notation, we can define the matrices as follows . R : M â V ..... Radford, Introduction to the Quantum Yang-Baxter equation and Quan- tum Groups: An a
Sep 29, 2017 - constitutive relations where explored in the context of metamaterials but they frequently assume a specific geometry. [32â36]. Also ... are looking at, and the considered geometry. ...... in textbooks concerned with generalized funct
Jan 30, 2013 - price and vice versa for selling) or switching their trading to a different stock where returns are potentially higher (Fig. 1). ... ability of resources as they choose to buy or sell their stocks (37). Thus, if a trader is willing to
Intensity-dependent dispersion under conditions of electromagnetically induced transparency in coherently prepared multi-state atoms. Andrew D. Greentree1,2, Derek Richards3 , J. A. Vaccaro4 , A. V. Durrant1 , S. R. de Echaniz5 , D. M. Segal5 , and J. P. Marangos5
arXiv:quant-ph/0209067v2 19 Mar 2003
Quantum Processes Group, Department of Physics and Astronomy, The Open University, Walton Hall, Milton-Keynes, MK7 6AA, United Kingdom 2 Centre for Quantum Computer Technology, School of Physics, The University of New South Wales, Sydney, NSW 2052, Australia 3 Quantum Processes Group, Department of Applied Mathematics, The Open University, Walton Hall, Milton-Keynes, MK7 6AA, United Kingdom 4 Department of Physics and Astronomy, The University of Hertfordshire, College Lane, Hatfield AL10 9AB, United Kingdom and 5 Quantum Optics and Laser Science, Blackett Laboratory, Imperial College of Science, Technology and Medicine, Prince Consort Road, London SW7 2BW, United Kingdom (Dated: April 2, 2018) Interest in lossless nonlinearities has focussed on the the dispersive properties of Λ systems under conditions of electromagnetically induced transparency (EIT). We generalize the lambda system by introducing further degenerate states to realize a ‘Chain Λ’ atom where multiple coupling of the probe field significantly enhances the intensity dependent dispersion without compromising the EIT condition. PACS numbers: 42.50.Gy, 42.50.Hz, 32.80.-t, 42.65.-k
There has been much interest given lately to the enhancement of optical nonlinearities in Electromagnetically Induced Transparency (EIT). Most of the work has focussed on the three-state system in the Λ configuration, which has provided some dramatic examples of nonlinear optical effects. Examples include ultra-slow , stopped  and superluminal  group velocities, coherent sideband generation  etc. All these nonlinear processes depend on the creation of coherent superpositions of the ground states with accompanying loss of absorption, and such mechanisms were described in . Thorough reviews of EIT and its properties can be found in  and . Recent investigations of nonlinear optics at the few or single photon levels have identified four state systems where the probe field simultaneously couples two transitions in the N configuration. Examples of applications for such work include photon blockade  and two photon absorptive switches . The classical precursors to such experiments have also been performed . Other experiments on the N scheme have been performed by ´ Entin et al. . In order to realize larger nonlinear effects, Zubairy et al.  suggested an extension where the more usual N configuration was extended to a system with an arbitrary (even) number of states where all the states are resonantly coupled except on the final transition where detuning is present. This scheme shows enhanced nonlinearities of not only χ3 but also higher order susceptibilities. One problem with this scheme and the standard N scheme is to do with the need to balance the required nonlinearity and decoherence in the system. To enhance the nonlinearity present in the systems it is important for the detuning of the final probe field to
be minimized, however decreasing the detuning increases the amount of the final excited state which is mixed into the coherent superposition state, resulting in an increase in decoherence and optical losses. This problem is to some extent circumvented by the absorptive switch of Harris and Yamomoto  by exploiting such losses, and in photon blockade by using the cavity nonlinearity to prevent absorption of the final photon. Still the increase in decoherence proves to be a difficulty in experimental precursors to these processes and causes problems in travelling wave configurations. An alternative multistate configuration for investigating EIT enhanced nonlinearities is the tripod configuration, studied recently by Paspalakis and Knight  and earlier considered by Morris and Shore . This system has many of the advantages of the N system, but by maintaining superposition states of the three ground states also avoids the excess decoherence of the N system. Morris and Shore  also mentioned a multi-leg extension of the tripod scheme. Here we present an alternative extension to the standard Λ configuration which we term as the Chain Λ configuration, depicted in Fig. 1. We start with a Λ atom [Fig. 1(a)] with ground states |g1 i and |g2 i and excited state |e1 i. The g1 − e1 (g2 − e1 )transition is excited by a probe (coupling) field of frequency ωp (ωc ) and detuning ∆p = ωp − ωe1 g1 (∆c = ωc − ωe1 g2 ). The probe (coupling) Rabi frequency is P (C). We assume for convenience that the transition frequencies ωe1 g1 and ωe1 g2 are equal and that appropriate selection rules ensure that the probe and coupling fields only interact with their designated states (a concrete example of how this can be achieved is described below).
n X H = (j − 1) (∆p − ∆c ) |gj i hgj | ~ j=1
n−1 X j=1
P3 … Cn |g3〉
n−1 X j=1
n−1 X j=1
FIG. 1: Energy level configurations for Chain Λ atoms. (a) is the usual Λ system, (b) is the 5 state Chain Λ or M system, (c) shows the generalization to higher number of states.
The five-state Chain Λ is an M system and is illustrated in Fig. 1(b). Here the ground states are labelled |g1 i, |g2 i, |g3 i, and the excited states |e1 i, |e2 i. The probe field excites the g1 − e1 and g2 − e2 transitions simultaneously with Rabi frequencies P1 and P2 respectively (different labels are applied to take account of the different coupling strengths of the transitions), whilst the coupling field excites the g2 − e1 and g3 − e2 transitions with Rabi frequencies C1 and C2 . We note that a complementary work by Matsko et al. , which looks at Faraday rotation and Kerr nonlinearities in the Λ, N and M schemes, has been performed which confirms some of our predictions about nonlinearities in these systems. An early study of the M scheme in the context of degenerate two level systems was also performed by Morris and Shore 
The Chain Lambda atom with 2n − 1 states is shown in Fig. 1(c). The ground states are denoted |g1 i, ..., |gn i, and the excited states |e1 i, ..., |en−1 i. The probe field excites the gj − ej transition with Rabi frequency Pj and the coupling field excites the gj+1 − ej transition with Rabi frequency Cj .
In order to gain insight into the problem we consider first the Hamiltonian for the Chain Λ atom with 2n − 1 states. Using the rotating wave approximation this can
[j∆p − (j − 1) ∆c ] |ej i hej | Pj (|gj i hej | + |ej i hgj |) Cj (|gj+1 i hej | + |ej i hgj+1 |) .
We ignore the decay from the excited states in our analytical analysis in order to gain simple expressions for the dressed states of the field-atoms system and thus gain a clearer understanding of the problem. Furthermore we shall concentrate our analysis on the optical nonlinearities which are present in the vicinity of the dark state, which are relatively insensitive to decay. This is evident by direct comparisons between numerical solutions of the complete master equation with decay, and our decay-free analytic expressions. The Hamiltonian can be conveniently expressed as a tridiagonal matrix with state ordering |g1 i, |e1 i, |g2 i, ..., |gn i, 0 P1 0 P1 ∆p C1 0 C ∆ 1 cp H = 0 0 P2 ~
0 0 P2 .. . ..
. (n − 1) ∆cp Pn 0 Pn n∆cp + ∆c Cn 0 Cn n∆cp
where we have introduced ∆cp = ∆p − ∆c . Following the approach taken by Kuang et al.  and Zubairy et al. , we first calculate the eigenvectors of the Hamiltonian. These can be written as |Di i = αi,g1 |g1 i + αi,e1 |e1 i + αi,g2 |g2 i + · · · + αi,gn |gn i , where i varies from 1 to (2n − 1). It is clearly not possible to give general solutions for the αi ’s for atoms with more than 3 states, although one may simply derive numerical results. However, if we invoke the adiabatic hypothesis  and assume that the probe detuning is small (∆p ≪ Pi , Ci ), the probe field is turned on slowly, and the coupling field resonant, then we may assume that the system evolves solely into the dressed state with energy closest to 0. For convenience we denote this state E (n) where n is the number of states in the Chain Λ D0
atom. Using MAPLE , and the simplification that Pi = P, CEi = C, i = 1, 2 . . . n we have derived expressions (n) for D0 . These are presented in table I as unnormal-
ized quantities and where Ω2 = C 2 + P 2 . Note that the
3 E (3) D0 C P
α0,g1 α0,e1 ∆p ΩC2 α0,g2 α0,e2
E (5) D0
C2 P2 C 2 (Ω2 +P 2 ) ∆p P (Ω4 −C 2 P 2 ) C −P C(2Ω2 −P 2 ) ∆p Ω4 −C 2 P 2
E (7) D0
C3 −P 3 C(Ω2 −P 2 )(Ω2 +2P 2 ) ∆p Ω2 (Ω4 −C 2 P 2 ) C2 P2 2C 2 Ω2 ∆p P Ω4 −C 2 P 2 ) ( C −P C (Ω4 +2C 4 )
(2C 2 P 2 −Ω4 ) 1
TABLE I: Unnormalized coefficients of |D0 i for 3, 5 and 7 state Chain Λ atoms.
E (3) results for D0 are equivalent to those which appear in Ref.  and were also used in Ref. . Short of directly creating an artificial atomic structure in a quantum well type material, or in an optical lattice, it is important to investigate whether the required Chain Λ structure is naturally present in any materials. A simple approximate realization to the Chain Λ configuration is obtained by exciting an F = n to F ′ = n transition in an atomic vapor where, as before, 2n− 1 is the number of states in the Chain Λ. We illustrate this for the five-state Chain Λ atom in Fig. 2, although the scheme generalizes in an obvious manner. We assume that the coupling field is σ − polarized (and therefore excites transitions from mz = n to mz = n − 1) and the probe σ + polarized (exciting transitions from mz = n to mz = n + 1). Notice that without any preparation there are two systems here, first the required M system (bold lines in Fig. 2), and second an undesired W system (dashed lines). In order to select the M over the W we first apply the coupling field, this has the effect of optically pumping the population into the mz = −2 state. Next the probe beam is turned on sufficiently slowly to ensure that the system evolves adiabatically to the desired dark state. We also note that the W system does not have any dark states. So even without the adiabatic state preparation, any population in the W system will eventually be optically pumped into the desired dark state of the M system. We also note that the presence of M systems have been identified in conjunction with Λ systems at least twice before in Refs.  . It is important to notice that our method of realizing Chain Λ systems and performing experiments with them, is not significantly more complex than standard experiments on simple Λ systems, all that is necessary to achieve the enhanced nonlinearities is the appropriate choice of transition. To study intensity-dependent dispersion it is necessary to extract the susceptibility at the probe frequency as a function of small probe detuning χ=κ
n X µ2gjej j=1
ρgj ej = κ
n X µ2gjej j=1
α∗0,gj α0,ej ,
where κ = 2πN /ǫ0 ~, N being the atomic density and
FIG. 2: A possible way to realize the M system using an F = 2 to F ′ = 2 transition. The coupling field is σ − polarized and turned on before the σ + polarized probe. In this way the desired M system (bold) is favored over the W system (dashed).
µgj ej the dipole moment of the gj − ej transition. The ∗ denotes complex conjugation. In order to calculate realistic values of the dispersion which would be attainable in standard experiments with alkali atoms (for example in a magneto-optical trap or vapor cell), we have taken N = 3×1015m−3 , µgj ej = 2×10−29Cm and Γ = 5.6MHz. In order to derive simple results for the nonlinear dispersions, we shall assume the coupling constants for all transitions to be equal, i.e. µgi ei = µ and so the probe (coupling) field Rabi frequencies are the same for all transitions Pi = P (Ci = C). Using MAPLE we can then derive expressions for the intensity dependent dispersion, R, for Chain Λ atoms of varying number of states: R(3) = R(5) = R(7) = R(n)
C2 , Ω4 C 2 Ω4 + 2P 2 Ω2 − 2P 4
, 2 (Ω4 − P 2 Ω2 + P 4 ) C 2 Ω8 + 4P 4 Ω4 − 8P 6 Ω2 + 4P 8 2
Ω4 (Ω4 − 2P 2 Ω2 + 2P 4 ) ! n X 1 ∂ ρg e . = P ∂∆p i=1 i i
We note that our results for R(3) are compatible with the intensity dependent group velocities derived in Ref. . The analytically determined dispersions are plotted in Fig. 3 as a function of P/Γ for C = 0.25Γ. In Fig. 4 we present a comparison between the analytical results for the intensity dependent dispersion and that obtained by solving the density matrix equation with decay given in the Appendix. The good agreement shows that our analytic approach is justified. If one simply performs a Taylor series expansion on our results for the R’s (which are non-perturbative), it is easy to show that the Chain Λ systems exhibit nonlinearities to all orders in P/C. The linear dispersion has previously been identified as being important in EIT systems (see for example Harris et al. ), and it is clear from Fig. 3 that the linear dispersion is identical for all Chain Λ atoms. From this we may conclude that in the limit of weak probe fields, the probe field cannot couple the levels
5 level 1.6
7 level 9 level
5 0 0
P/Γ FIG. 3: Graphs showing dispersion (times linewidth squared), RΓ2 as a function of P/Γ with C/Γ = 0.25 for Chain Λ systems of 3,5,7 and 9 states.
FIG. 5: Graphs showing log10 of group velocity, vg as a function of P/Γ with C/Γ = 0.25 for Chain Λ systems of 3,5,7 and 9 states.
the value of this ratio at the maximum then for P/C < β the dispersion will be monotonically increasing with increasing P , and for P/C > β it monotonically decreases. It is clear that reciprocal results will be obtained for the corresponding group velocities. For 3, 5, 7 and 9 state atoms, the values of β are 0, 0.476, 0.698 and 0.804 (to 3 significant figures) respectively. One material property dependent on the dispersion is the group velocity, which is
18 16 14
12 10 8 6
4 2 0 0
P/Γ FIG. 4: Graph showing the comparison between intensity dependent dispersion calculated using the full master equation (solid line) and the analytical approach using dressed states (dashed line) in a five-state system.
in any fashion other than the simple Λ scheme. However as the probe intensity is increased, higher-order process begin to turn on. In order to rigorously determine the order of the nonlinearities present, one should construct an effective Hamiltonian, following methods presented by, for example, Zubairy et al.  or Klimov et al. . This has not yet been performed for the Chain Λ system and it is hoped that such investigations will shed more light on the nonlinear optical properties of these systems. An interesting feature to note in the dispersion calculations is the Rabi frequency ratio which provides the maximum dispersion. In this simple case, the position of the maximum depends only upon the ratio P/C . If β is
c 1 + ωp ∂ℜ(η) ∂∆p
√ where η = 1 + χ is the complex refractive index. It is important to realize that in systems with nonlinear dispersions as large as those for Chain Λ systems, the group velocity may be a poor parameter. This is because for realistic propagation through an optically thick medium, the intensity dependence of the medium will alter the shape of a simple gaussian pulse. There are however, other experiments sensitive to the group velocity which may be considered, for example bichromatic excitation of the probe beam to generate a beat note  or the use of a frequency modulated (rather than the more usual amplitude modulated) probe signal. In Fig. 5 we present intensity dependent group velocities for 3, 5, 7 and 9 state Chain Λ atoms corresponding to the dispersion calculations presented in Fig. 3. The analytical results given above cannot provide a complete description of the physically realizable problem, because of the different Clebsch-Gordan coupling between the states involved in the transitions. If we use the scheme suggested in Fig. 2 for the couplings and define our coupling strengths relative to the coupling in the first Λ system (i.e. g1 − e1 − g2 ) then we may write down the expressions for the other couplings in terms of these
P2 C2 P3 C3 µg2 e2 µg3 e3
5-state p 3/2 P1 p 2/3 C1 p
3/2 µg1 e1 −
stand propagation effects and especially the effect of such high nonlinearities on pulse shape. Such analysis goes beyond our simple picture and will be the focus of future work.
TABLE II: Relevant Rabi frequency ratios for 5 and 7 state Chain Λ systems 40
One of the authors (AG) would like to acknowledge useful discussions with Dr A.B. Matsko (Jet Propulsion Laboratory) and financial support from the EPSRC (UK).
7 level 20
5 level 10
P1/Γ II. FIG. 6: Graphs of intensity dependent dispersions as a function of P1 , calculated using realistic ratios of coupling strengths, i.e. with the relative Pi and Ci ’s calculated using the Clebsch-Gordan ratios in Table II.
quantities . These are summarized in table II for the 5 and 7 state Chain Λ atoms, and it is easy to generalize to higher orders. Results of numerical calculations of intensitydependent dispersions for 5- and 7-state Chain Λ systems are presented in Fig. 6. Comparing Figs. 3 and 6, shows that despite the differences in the values obtained for the dispersions, there is only minimal change to the overall shape of the intensity dependent curves. We have shown that there exist interesting nonlinear properties for Chain Λ atoms, and in particular we have focussed on the intensity-dependent dispersion as a measure for these nonlinear optical properties. The nonlinearity of these systems increases as the number of atomic states increases, whilst the EIT transparency is maintained. It therefore appears likely that such multistate systems will be useful in the search for new quantum non-linear optical materials. Our analysis has been confined to the optically thin regime. Clearly, in a full study, which would include intensity-dependent group velocities in such highly nonlinear media, it is important to under-
In this appendix we present the density matrix equations of motion for the 5-state Chain Λ atom (M scheme). The equations to be solved are:
where Γn is the n photon dephasing and Γ = Γe1 = Γe2 is the total decay rate from either excited state. A three-dimensional plot showing −ℑ(χ) as a function of probe detuning, ∆p /Γ and probe Rabi frequency, P/Γ
FIG. 7: Three-dimensional plot showing −ℑ[χ] (proportional to the probe absorption) as a function of ∆p /Γ and P/Γ for P1 = P2 = P , C1 = C2 = C = Γ/4, Γ = 5.6MHz for the five-state Chain Λ system.
for P1 = P2 = P , C1 = C2 = C = Γ/4, Γ = 5.6MHz and other parameters as above, for the five-state Chain Λ system is presented in Fig. 7.
 M. M. Kash, V. A. Sautenkov, A. S. Zibrov, L. Hollberg, G. R. Welch, M. D. Lukin, Y. Rostovtsev, E. S. Fry, and M. O. Scully, Phys. Rev. Lett. 82, 5229-5231 (1999); L. V. Hau, S. E. Harris, Z. Dutton, and C. H. Behroozi, Nature (London) 397, 594-8, (1999); J. P. Marangos, ibid. 397, 559-560 (1999).  O. Kocharovskaya, Y. Rostovtsev, and M. O. Scully, Phys. Rev. Lett. 86, 628 (2001); C. Liu, Z. Dutton, C. H. Behroozi, and L. V. Hau, Nature (London) 409, 490 (2001); D. F. Phillips, A. Fleischhauer, A. Mair, R. L. Walsworth, and M. D. Lukin, Phys. Rev. Lett. 86, 783 (2001); A. V. Turukhin, V. S. Sudarshanam, M. S. Shahriar, J. A. Musser, B. S. Ham, and P. R. Hemmer, ibid. 88, 023602, (2002).  A. M. Steinberg and R. Y. Chiao, Phys. Rev. A 49, 2071 (1994); L. J. Wang, A. Kuzmich, and A. Dogariu, Nature (London) 406, 277 (2000); A. Dogariu, A. Kuzmich, and L. J. Wang, Phys. Rev. A 63, 053806 (2001).  H. Schmidt and R. J. Ram, Appl. Phys. Lett. 76, 3173 (2000), A. F. Huss, N. Peer, R. Lammeggar, E. A. Korsunsky, and L. Windholz, Phys. Rev. A 63, 013802 (2001).  S. E. Harris, J. E. Field, and A. Imamo˘ glu, Phys. Rev. Lett. 64, 1107 (1990).  E. Arimondo, Progress in Optics (Amsterdam: Elsevier 1996) Vol. XXXV, edited by E. Wolf, p. 257; J.P. Marangos, J. Mod. Opt. 45, 471-503 (1998).  A. B. Matsko, O. Kocharovskaya, Y. Rostovtsev, G. R. Welch, A. S. Zibrov, and M. O. Scully, Adv. At. Mol. Opt. Phys. 46, 191 (2001).  A. Imamoglu, H. Schmidt, G. Woods, and M. Deutsch, Phys. Rev. Lett., 79, 1467 (1997); S. Rebi´c, S. M. Tan, A. S. Parkins, and D. F. Walls, J. Opt. B: Quantum Semiclass. Opt. 1, 490 (1999); K. M. Gheri, W. Alge, and P. Grangier, Phys. Rev. A 60, R2673 (1999); A. D. Greentree, J. A. Vaccaro, S. R. de Echaniz, A. V. Durrant and J. P. Marangos, J.Opt. B: Quantum Semiclass. Opt. 2, 252-259 (2000); S. Rebi´c, A. S. Parkins, and S. M. Tan, Phys. Rev. A 65, 043806 (2002); S. Rebi´c, A. S.
Parkins, and S. M. Tan, ibid. 65, 063804 (2002).  S. E. Harris and Y. Yamamoto, Phys. Rev. Lett. 81, 3611 (1998).  M. Yan, E. G. Rickey, and Y. Zhu, Opt. Lett. 26, 548 (2001); M. Yan, E. G. Rickey, and Y. Zhu, Phys. Rev. A 64, 013412 (2001); M. Yan, E. G. Rickey, and Y. Zhu, ibid. 64, 041801 (2001); M. Yan, E. G. Rickey, and Y. Zhu, J. Mod. Opt. 49, 675 (2002).  S. R. de Echaniz, A. D. Greentree, A. V. Durrant, D. M. Segal, J. P. Marangos, and J. A. Vaccaro, Phys. Rev. A 64, 013812 (2001).  V. M. Entin, I. I. Ryabtsev, A. E. Boguslavskii, and I. ´ M. Beterov, Pis’ma. Zh. Eksp. Teor. Fiz. 71, 257 (2000) [JETP Lett. 71, 175 (2000)].  M. S. Zubairy, A. B. Matsko, and M. O. Scully, Phys. Rev. A 65, 043804 (2002).  E. Paspalakis and P.L. Knight, J. Mod. Opt. 49, 87 (2002).  J.R. Morris and B.W. Shore, Phys. Rev. A 27, 906 (1983).  A.B. Matsko, I. Novikova, G.R. Welch, and M.S. Zubairy, Enhancement of Kerr nonlinearity via multi-photon coherence, preprint arXiv:quant-ph/0207141.  L. Kuang, G. Chen, and Y. Wu, Giant non-linearities accompanying electromagnetically induced transparency, preprint arXiv:quant-ph/0103152 v2 (2001).  MAPLE 8, Waterloo Maple Inc., 2002.  S. E. Harris, J. E. Field, and A. Kasapi, Phys. Rev. A 46, R29 (1992).  A. B. Klimov, L. L. S´ anchez-Soto, A. Navarro, and E. C. Yustas, Effective Hamiltonians in quantum optics: a systematic approach, preprint arXiv: quant-ph/0209013 v1 (2002).  There are many good routines available on the internet. We used one at http://www.ph.surrey.ac.uk/˜phs3ps/cgjava.html to calculate our results.