quantum communication and computation offer technological advances by performing information processing with quantum mechanics  @xcite . in order to perform these tasks over large distances , photons are the obvious choice to carry quantum information .
however , transmission over large distances is dramatically hampered by attenuation in optical fibers .
unlike classical fiber optic communication , amplifiers placed periodically along the transmission channel can not be used to overcome this loss owing to the  no - cloning " theorem  @xcite .
fortunately , quantum repeaters have been proposed to resolve this issue by using entanglement shared between adjacent nodes and joint measurements to create entanglement between the start node and the terminal node  @xcite .
a crucial element in this scheme is the ability to efficiently store and retrieve single photons using a quantum memory  @xcite . to date
, the essential ingredients of a photonic quantum memory have been demonstrated using ensembles of ultracold alkali atoms  @xcite , ensembles of ultracold atoms in cavities  @xcite , warm vapors of alkali atoms  @xcite , ultracold single atoms in cavities  @xcite , and solid - state systems composed of rare - earth dopants in crystals  @xcite . specifically , retrieval efficiencies as high as 73@xmath0 and storage times as long as 3.2  ms have been simultaneously demonstrated using ultracold atoms  @xcite .
the combination of large optical depths and long ground state hyperfine coherence make ultracold atomic ensembles an attractive platform for optical quantum memories .
while there are quantum memory schemes which use probabilistically - generated spin waves from spontaneous raman scattering in atomic ensembles , it has been shown that schemes based on fast ( rates approaching the ghz scale ) single photon or entangled photon sources can provide better performance  @xcite .
thus , one would ideally like a source of quantum light states that is on - demand and bright , meaning it can produce these photons on fast timescales ( broadband ) in a triggered fashion .
in addition , these broadband photons should be indistinguishable so they bunch perfectly on a beamsplitter  @xcite .
a promising source of such on - demand , non - classical light states are epitaxially grown , gaas - based semiconductor quantum dots ( qds )  @xcite .
the short spontaneous emission lifetime of a qd , more than an order of magnitude shorter than that of an alkali atom , and the fact that it is embedded in a robust optoelectronic material make it an attractive candidate as a quantum light source . as an individual two - level system
, the qd naturally emits one photon when it is excited and collection of this photon can be very efficient by proper design of the surrounding dielectric  @xcite .
single photon count rates can in principle approach 1 ghz , with higher rates possible by taking advantage of a purcell enhancement .
in addition , qds have been shown to emit polarization entangled photons by means of a cascaded decay  @xcite . finally , qd - generated photons can have a high degree of indistinguishability  @xcite . aside from these optical properties ,
the internal spin states of charged quantum dots have received considerable attention for quantum information processing  @xcite .
spin coherence times as along as 1  @xmath2s have been measured  @xcite and because the spin state can be entangled with the polarization of an emitted photon  @xcite , charged quantum dots are a natural candidate for a solid - state qubit that interacts strongly with light .
while qds are promising sources of single or entangled photons , they must be coupled to a high quality quantum memory in order to be a viable source for long distance quantum information processing . here , we investigate the storage and retrieval of a broadband , qd - generated single photon with a 1 ns lifetime in an ultracold , dense ensemble of @xmath1rb atoms .
taking previously measured experimental parameters , we find that the total efficiency ( @xmath3 , where @xmath4 and @xmath5 are the storage and retrieval efficiencies ) can exceed 28@xmath0 for storage and backwards retrieval of a photon with a 1  ns lifetime in a @xmath1rb ensemble with an on - resonance optical depth of 150 .
because the bandwidth of the qd photon can approach the excited - state hyperfine splitting in @xmath1rb , this result was obtained by extending the @xmath6-system theory of gorshkov _ et al _ @xcite to a four - level system . using the gradient ascent approach outlined in ref .
@xcite , we determine the maximum efficiency as well as the optimal control pulse for photon storage .
we find that the control pulse can be easily generated using a 12  mw laser diode in contrast to the broadband , raman - based scheme in ref .
@xcite where a mode - locked laser is required . finally , we consider the effects of excess dephasing and spectral wandering of the qd optical transition and show that the memory efficiency remains fairly robust .
these results show that with existing technology , qd - generated photons can be reliably interfaced with ultracold atomic ensembles , paving the way for their future use in quantum information as well as opening interesting avenues in the study of hybrid quantum systems .
this paper is arranged as follows : in sec .
[ sec:3l ] , we briefly review the physics of photon storage and retrieval in a @xmath6-system following the treatment of gorshkov _ et al _ @xcite . in sec .
[ sec:4l ] , we extend this treatment to a four - level atom , which is relevant for broadband photon storage in atomic systems with non - negligible hyperfine structure in the excited state .
section [ sec : qd ] implements the model to study storage and retrieval of a broadband single photon emitted by a quantum dot using a @xmath1rb ensemble . in sec .
[ sec : dephase ] , deleterious effects resulting from imperfect indistinguishability of the quantum dot photon are discussed . in sec .
[ sec : highod ] we investigate how the efficiencies increase for very high optical depth .
the work is concluded and summarized in sec .
[ sec : conc ] and details of the calculations and numerical implementation are discussed in the appendices a and b. appendix c discusses the role of four - wave mixing ( fwm ) in the storage process .
photon storage in a three level , @xmath6-type system has been the subject of many articles and reviews  @xcite . here
, we briefly touch on the main points and restrict the discussion to those schemes described by the theory of gorshkov _ et al _ , namely those based on electromagnetically induced transparency ( eit ) , off - resonant raman , or photon echo interactions  @xcite . in that work
, these three schemes were shown to yield similar storage efficiencies for the same optical depth of the @xmath6 medium and are in this sense equivalent
. however , there are tradeoffs in the actual physical implementation for these schemes which will be discussed later .    the basic picture of photon storage in such a system is shown in fig .
[ fig : fig1 ] .
each atom in the ensemble is composed of two ground states @xmath7 and @xmath8 that can be optically coupled to the excited state @xmath9 .
these two transitions can be individually addressed using a sufficient ground state splitting or by selection rules so as to avoid cross - coupling .
the storage procedure starts with all of the atoms initialized into @xmath7 .
then a quantum field @xmath10 which is to be stored , addresses the @xmath11 transition while a classical control field @xmath12 addresses the @xmath13 transition ( see fig .
[ fig : fig1](a ) ) .
the control pulse facilitates the transfer of the quantum field to a spatially - dependent coherence of the two ground states ( often called a  spin wave " ) @xmath14 as depicted in fig .
[ fig : fig1](b ) . the exact physical mechanism that accomplishes this task depends on the scheme . in eit storage
, the control pulse dynamically creates a transparency window in the absorption profile of the @xmath11 transition while adiabatically reducing the group velocity of the quantum field to 0  @xcite . in raman storage
, the control pulse enables the quantum field to be absorbed into @xmath8 by a two - photon raman transition  @xcite . for photon echo storage ,
the quantum field promotes an atom to the excited state @xmath9 and the control pulse transfers the excitation to @xmath8 by performing a fast @xmath15 pulse . in all cases , the quantum field transfers one atom from @xmath7 to @xmath8
, however there is no knowledge of which atom . thus the excitation is coherently distributed over the entire ensemble and it is this collective behavior that makes the process efficient .    the quantum field is now stored in this collective ground state coherence and is therefore sensitive to decoherence processes .
these processes , which can include magnetic field fluctuations or atomic motion , set the limit for the spin wave coherence time , and hence , the duration the photon can be stored .
notably , the spin wave is more robust against decoherence than many other multi - particle entangled states  @xcite
. the spin wave can be re - converted into a photon by using another control pulse ( see fig .  [
fig : fig1](c)-(d ) ) .
et al _ showed that the optimal retrieval is simply the time - reverse of storage ; so - called  backwards retrieval "  @xcite .
this process creates a photon propagating in the opposite direction to its initial propagation .
forwards retrieval is also possible , but will not be treated here .    for all of these storage schemes , the dynamics of the interaction between the two optical fields and the atoms
is described by treating the signal field on the @xmath11 transition quantum mechanically while treating the @xmath13 control field semiclassically .
the atomic level structure with relevant energy scales for photon storage in a three - level @xmath6 system is shown in fig .
[ fig : fig2 ] .
the ensemble of @xmath6 systems is composed of @xmath16 atoms distributed over a volume of length @xmath17 and cross - sectional area @xmath18 with linear density @xmath19 . in this analysis
, we ignore the motion of the atoms and thereby restrict the discussion to ultracold ensembles . as was shown in ref .
@xcite , the important parameter governing the efficiency of the storage process is the optical depth , defined as @xmath20 where @xmath21 is the single - photon coupling constant for the transition and @xmath22 is the decay of the coherence of the excited state ( for a purely radiative decay with no additional dephasing @xmath23 , the spontaneous emission rate ) . note that the @xmath24 used here and in ref .
@xcite is equal to half of the optical depth as usually defined ( @xmath25 where @xmath26 is the resonant absorption cross - section of a single atom and @xmath27 is the number density ) .
the classical control field is described by a rabi frequency envelope @xmath28 centered at frequency @xmath29 .
as also noted in app .
a , this rabi frequency is defined to be 1/2 of the usually defined rabi frequency .
the quantum field is described by a slowly varying envelope operator @xmath30 centered at frequency @xmath31 .
in addition to @xmath30 , there are two other operators required to describe the dynamics : the polarization operator @xmath32 and the spin - wave operator @xmath33 where @xmath34 are slowly - varying collective atomic operators defined in app .
a. it was shown in ref .
@xcite that in order to determine normally - ordered quantities such as efficiencies , one can neglect quantum noise operators and treat all dynamical variables as complex numbers .
the equations of motion governing @xmath35 , @xmath36 , and @xmath37 are derived by calculating the heisenberg equation of motion for all dynamical variables using the dipole and rotating wave approximations .
then , two further approximations are made related to the fact that the quantum field is weak .
first , it is assumed that almost all of the atoms remain in the ground state @xmath7 during the whole process .
second , only terms to linear order in @xmath35 are retained . under these approximations
( see ref .
@xcite and app .
a for details ) , the following equations of motion are obtained :    @xmath38    these equations of motion use a coordinate system @xmath39 , where @xmath40 is a dimensionless length parameter ( @xmath41 $ ] ) and @xmath42 is the time in a co - moving reference frame .
compared to the results of ref .
@xcite , eqs .
[ eqn:3leome]-[eqn:3leoms ] also include relative dipole moments @xmath43 .
these are defined as dipole moments relative to that of the two - level , cycling transition @xmath44 where a measurement of the optical depth would take place .
this sets the natural scale of the atom s dipole strength and allows for easy comparison between different physical implementations of levels @xmath7 , @xmath8 , and @xmath9 ( see app .
a for details ) .    for storage
, we would like to take a quantum field with initial envelope @xmath45 ( non - zero on the interval @xmath46 $ ] ) and map it into a spin wave @xmath47 using a classical pulse @xmath12 .
the boundary conditions for the dynamical variables are @xmath48 and @xmath49 .
we then want to compute the efficiency of this mapping once @xmath50 at @xmath51 .
if the envelope of the initial quantum field is normalized ( @xmath52 ) , the storage efficiency @xmath4 is given by  @xcite @xmath53    for a given input field and optical depth of the medium , one would like to determine the optimal classical pulse shape so as to maximize the storage efficiency . in ref .
@xcite it was shown that in the adiabatic limit , @xmath54 , an analytic solution for @xmath12 could be found and that the efficiency scaled as @xmath55 .
for broadband photon storage , the adiabatic limit is not necessarily met and the optimal @xmath12 must be found numerically .
_ used a gradient ascent algorithm in ref .
@xcite to numerically optimize @xmath4 and found that the results matched the analytical solution in the adiabatic limit .
in addition , it was shown that photons beyond the adiabatic regime ( @xmath56 ) could also be stored efficiently by using this optimization technique .
this is the approach we will use in the following analysis and now briefly review .
because we take a numerical approach , we now restrict the discussion to near - resonant storage schemes and exclude off - resonant raman based storage which would require a much larger computational domain so as to capture the rapidly varying detuning ( @xmath57 ) . gorshkov _
et al _ mentioned in ref .
@xcite how one can obtain the optimal control pulse for raman - based storage from the optimal control pulse found for resonant storage , but we do not pursue this here .
we do , however , note that the control pulse required for raman storage can be several orders of magnitude more intense than that required for resonant schemes  @xcite for the same optical depth @xmath24 and hence the same efficiency .
a gradient ascent algorithm simply starts with a trial solution and proceeds to the optimal solution by moving along the gradient of the quantity to be maximized . at each step in the algorithm ,
the gradient is determined and the control pulse is updated .
mathematically , this replacement rule is @xmath58 where @xmath59 is the step size parameter and @xmath60 is the quantity to be maximized . in the case of photon storage , we want to maximize the storage efficiency subject to the constraint that all of the dynamical variables fulfill the equations of motion ( eq .  [ eqn:3leome]-[eqn:3leoms ] ) .
thus , @xmath60 takes the form @xmath61 + c.c.\right\ } \nonumber \\ & + \int_0^t d\tau \int_0 ^ 1 d\tilde{z}\ , \big\ { \bar{p}_1^*\big [ \left(-\partial_{\tau}+i \delta_g-\gamma \right ) p_1+i \tilde{\mu}_{1s } \omega(\tau ) s   \nonumber \\ & + i \tilde{\mu}_{1 g } \sqrt{d \gamma } \mathcal{e } \big ] + c.c . \big\ } \nonumber \\ & + \int_0^t d\tau \int_0 ^ 1 d\tilde{z}\,\big\ { \bar{s}^*\big[-\partial_{\tau}s+ i \left(\delta_g - \delta_s \right ) s \nonumber \\ & + i \tilde{\mu}_{1s}\omega^*(\tau)p_1   \big ] + c.c .
\big\ } , \label{eqn:3lj}\end{aligned}\ ] ] where the first term is @xmath4 and lagrange multipliers @xmath62 , @xmath63 , and @xmath64 have been introduced to include the equations of motion .
the maximum storage efficiency is found when @xmath60 is stationary with respect to variations in all dynamical variables and @xmath12 .
requiring stationarity with respect to variations in @xmath35 , @xmath36 , and @xmath37 results in equations of motion and boundary conditions for the lagrange multipliers @xmath62 , @xmath63 , and @xmath64 . the equations of motion are @xmath65 with boundary conditions @xmath66 and @xmath67 . as pointed out in ref .
@xcite , these are exactly the equations of motion and boundary conditions for backwards retrieval .
thus , solving these equations will yield the retrieval efficiency , @xmath5 , and the total efficiency @xmath68 for storage followed by backwards retrieval .
explicitly , these quantities can be determined using @xmath69 where @xmath70 is the output quantum field .
the variation of @xmath60 with respect to variations in @xmath12 ( the gradient ) can be identified from eq .
[ eqn:3lj ] as @xmath71}. \label{eqn:3lgrad}\ ] ] thus , the prescription for obtaining the maximum storage efficiency and optimal control pulse for a given input quantum field and optical depth is as follows .
first , take a trial control pulse and solve the equations of motion eq .
[ eqn:3leome]-[eqn:3leoms ] with the storage boundary conditions to obtain @xmath72 and @xmath14 .
then , solve eq .
[ eqn:3leombe]-[eqn:3leombs ] using the boundary conditions for backwards retrieval ( note that these equations run backwards in time and space ) to obtain @xmath73 and @xmath74 .
now the gradient can be determined from eq .
[ eqn:3lgrad ] and the control pulse can be updated using the replacement rule .
this process is then repeated until the desired tolerance of the storage ( or total ) efficiency is obtained .     in blue ( red ) .
solid curves are guides to the eye . ]    as an example we consider the case of an idealized rb gas ( @xmath75 , @xmath76 mhz ) storing a quantum - dot photon with @xmath77 where @xmath78 is the heaviside step function and @xmath79 is the spontaneous emission lifetime ( taken to be 1  ns ) .
we take both fields to be exactly on resonance @xmath80 . as a function of @xmath24 , the storage and total efficiencies
are determined by performing gradient ascent for each value of @xmath24 .
details and limitations of the numerical implementation can be found in appendix  [ app : num ] .
the results of the numerical calculation are shown in fig .
[ fig : fig3 ] where the storage ( blue ) and total ( red ) efficiencies are plotted as a function of @xmath24 .
both curves reach an asymptotic value of @xmath8196.0@xmath0 for @xmath82 limited by the finite computational domain and the instantaneous rise of the quantum field , which is non - physical but chosen for calculational simplicity ( see app .  [
app : num ] ) .     determined using gradient ascent .
the rabi frequencies are scaled by @xmath76 mhz . ]    in addition to determining the efficiencies , it is instructive to see how the optimal control pulse @xmath12 changes as a function of @xmath24 .
figure  [ fig : fig4 ] shows the optimal control pulse for a few selected values of @xmath24 and it can be seen that the required peak rabi frequency @xmath83 increases with @xmath24 . for these parameters , we find that @xmath84 .
furthermore , @xmath12 starts to become much more heavily weighted at the start of the pulse . both of these observations can be understood as the optimal storage scheme slowly changes from a photon echo type storage towards an eit - based storage as the adiabatic limit is approached .
lastly , the widths of these control pulses are not ultrafast , i.e. they do not require a mode - locked laser but rather are achievable with more flexible techniques such as direct intensity modulation .
the combination of relatively high total efficiencies and realistic control pulses make it seem that storage of a quantum dot photon in an atomic gas like rb is certainly feasible with current technology .
importantly , the peak rabi frequency @xmath83 of the control pulse surpasses 50@xmath22 even for relatively low @xmath24 .
for such large values of @xmath83 , the more complicated level structure of the storage medium will begin to play a significant role in the dynamics .
for instance , the excited - state splittings of the commonly - used @xmath85 transitions of @xmath1rb are around 100 mhz ( @xmath86 ) . as a result ,
the three - level @xmath6 treatment is insufficient to describe broadband , on - resonance storage and more levels must be taken into account .
we now extend the treatment of sec .  [ sec:3l ] to include another excited state that can couple to both the control field and the input quantum field as shown in fig .
[ fig : fig5 ] . because the maximum rabi frequency of the control pulse can exceed the excited state level splitting @xmath87
, an atom can go from @xmath7 to @xmath8 via @xmath9 or @xmath88 . depending on the dipole matrix elements @xmath43 , these paths can interfere constructively or destructively .    following the same procedure as in sec .
[ sec:3l ] and adding another polarization @xmath89 results in the following equations of motion ( see app .
a for details ) :    @xmath90 \label{eqn:4leome}\\ \partial_{\tau}p_1 & = \left(i \delta_g-\gamma \right ) p_1+i \tilde{\mu}_{1s } \omega(\tau ) s + i \tilde{\mu}_{1 g } \sqrt{d \gamma } \mathcal{e } \label{eqn:4leomp1}\\ \partial_{\tau}p_2 & = \left(i\delta_g -i\delta_e-\gamma\right ) p_2 + i \tilde{\mu}_{2s } \omega(\tau ) s \nonumber \\ & + i \tilde{\mu}_{2 g } \sqrt{d \gamma } \mathcal{e } \label{eqn:4leomp2}\\ \partial_{\tau}s & = i \left(\delta_g - \delta_s \right ) s + i\omega^*(\tau)\left [ \tilde{\mu}_{1s}p_1 + \tilde{\mu}_{2s}p_2 \right ] \label{eqn:4leoms}\end{aligned}\ ] ]    where new relative dipole moments @xmath91 and @xmath92 have been introduced to represent the couplings to the additional excited state .
as the spin wave @xmath14 and the field @xmath93 are unchanged , the definitions of the storage , retrieval , and total efficiencies remain the same .
furthermore , the boundary conditions for storage are the same as before except @xmath94 is added .
we can also find the four - level version of @xmath60 @xmath95 + c.c.\big\ } \nonumber \\ & + \int_0^t d\tau \int_0 ^ 1 d\tilde{z}\ , \big\ { \bar{p}_1^*\big[(-\partial_{\tau}+i \delta_g-\gamma ) p_1+i \tilde{\mu}_{1s } \omega(\tau ) s   \nonumber \\ & + i \tilde{\mu}_{1 g } \sqrt{d \gamma } \mathcal{e } \big ] + c.c . \big\ } \nonumber \\ & + \int_0^t d\tau \int_0 ^ 1 d\tilde{z}\ , \big\ { \bar{p}_2^*\big [ \left(-\partial_{\tau}+i\delta_g - i\delta_e-\gamma   \right ) p_2   \nonumber \\ & + i \tilde{\mu}_{2s } \omega(\tau ) s+
i \tilde{\mu}_{2 g } \sqrt{d \gamma } \mathcal{e } \big ] + c.c . \big\ } \nonumber \\ & + \int_0^t d\tau \int_0 ^ 1 d\tilde{z}\,\big\ { \bar{s}^*\big[-\partial_{\tau}s+ i
\left(\delta_g - \delta_s \right ) s \nonumber \\ & + i\omega^*(\tau)\left ( \tilde{\mu}_{1s}p_1+\tilde{\mu}_{2s}p_2   \right ) \big ] + c.c .
\big\ } , \label{eqn:4lj}\end{aligned}\ ] ] in order to obtain the optimal control and efficiencies . by requiring
that @xmath60 is stationary with respect to variations in @xmath35 , @xmath36 , @xmath96 , and @xmath37 , we find the equations of motion @xmath97 for the lagrange multipliers and their boundary conditions ( @xmath66 and @xmath98 ) .
the gradient along @xmath12 is now modified to @xmath99 .
\label{eqn:4lgrad}\end{aligned}\ ] ]    it is clear that if @xmath100 , all of the expressions from sec .
[ sec:3l ] are recovered . in the limit where @xmath101 , then @xmath96 can be adiabatically eliminated and the results of sec .
[ sec:3l ] are again recovered .
thus , the extension to a four level system reproduces the results of the three level case in the appropriate limits .        to gain some insight into the ramifications of including a fourth level in an intermediate regime of @xmath87 , we directly compare the results of the gradient ascent algorithm applied to the three level case to that for two scenarios of the four level case . for all calculations , we again store a resonant quantum dot photon with @xmath77 using a resonant control pulse ( @xmath80 ) . for the parameters of the storage medium
we take @xmath76 mhz , @xmath102 , and @xmath103 .
the first four level scenario ( 4l+ ) we consider is @xmath104  mhz and @xmath105 while the second scenario ( 4l- ) has the same excited state energy splitting but @xmath106 .
the optimal control pulses resulting from the gradient ascent optimization are plotted in fig .
[ fig : fig6 ] for each scenario . not only does the inclusion of a fourth level dramatically change the optimal control pulse , but the sign of the relative dipole moment also has a large effect .
more importantly , the computed storage ( total ) efficiencies are 73.6@xmath0 , 77.6@xmath0 , and 43.5@xmath0 ( 63.4@xmath0 , 65.7@xmath0 , and 26.3@xmath0 ) for the 3l , 4l+ , and 4l- scenarios respectively .
if one naively uses the optimal control for 3l in the 4l+ case , the storage ( total ) efficiency is only 56.5@xmath0 ( 48.4@xmath0 ) , while for the 4l- case 20.8@xmath0 ( 4.8  @xmath0 ) is obtained . from this analysis
it is clear that additional excited states must be taken into account and a careful choosing of those levels is required for any implementation of broadband photon storage where the peak rabi frequency of the control pulse is comparable to the excited state splitting .    .
the slight asymmetry about @xmath107 is due to the finite accuracy of the numerics ( see appendix  [ app : num ] ) . ]    another question one could ask is if there is an optimal choice of the optical detuning @xmath108 in the four level system . to address this
, we again apply the gradient ascent optimization to the 4l+ configuration described above , but we vary the detuning @xmath109 from @xmath110 to @xmath111 in steps of @xmath112 .
the storage efficiency is plotted in fig .
[ fig : fig7 ] as a function of @xmath109 .
surprisingly , the efficiencies are smallest near the mid - point between @xmath9 and @xmath88 at @xmath113 , while the maxima occur near @xmath114 and @xmath115 .
nonetheless , the change in storage efficiency over the entire range is within @xmath116 , and there is not much to be gained or lost by changing the detuning at least for this set of @xmath43 .
we have also verified that varying the detuning has a small effect ( @xmath117 variation ) when using the dipole moments for the relevant @xmath1rb transitions . with these mathematical and numerical tools in hand
, we can address storage of a quantum dot photon using the actual atomic levels of @xmath1rb .
because of the long hyperfine ground - state coherence time and the ability to create cold , dense gases , @xmath1rb has proven to be a natural choice for implementation of quantum memory schemes .
indeed , both the 5@xmath118p@xmath119 and 5@xmath118p@xmath120 excited state manifolds corresponding to the @xmath121 and @xmath85 optical transitions have been used for quantum memory applications . before assigning @xmath1rb hyperfine states
@xmath122 to the various states of the four level model , the problem of frequency matching of a qd transition to the @xmath85/@xmath121 transitions at 780/795  nm must be addressed .
the most commonly studied self - assembled qds are composed of in@xmath123ga@xmath124as islands embedded in a gaas matrix , which are made to emit light in the 900 - 1000  nm band . in such qds ,
the so - called  wetting layer " ( a thin quantum well ) defines a barrier near 850  nm , below which no confined states exist in the dot .
therefore , the photons produced by these qds can not be stored in @xmath1rb , but cs transitions ( 852/895 nm ) are close to within reach and @xmath125yb@xmath126 ( 935 nm ) has a transition compatible with qds  @xcite .
nonetheless , efficient quantum frequency conversion techniques have been demonstrated using qds  @xcite such that the single photon produced by a qd could be frequency translated to another wavelength without destroying its quantum characteristics .
such techniques could bridge the frequency gap between in@xmath123ga@xmath124as qd transitions and those of @xmath1rb .
in addition , there are qds composed of gaas embedded in algaas and ingaas qds embedded in algaas which have been shown to emit in the 780  nm - 795  nm region  @xcite .
the optical properties of these qds are much less well - known compared to standard ingaas / gaas qds and detailed investigations are ongoing .
thus , either use of quantum frequency conversion techniques or proper material choices permit the study of broadband qd photon storage and retrieval in @xmath1rb .
notably , photons from a qd have been made to interact with @xmath1rb atoms in a recent experiment  @xcite . in this case
the atoms acted as a passive medium whose dispersion near the @xmath85 transition was used to reduce the group velocity of the qd photons . in this present work ,
we are interested in actively manipulating the atoms to controllably store and retrieve the qd photon .
we proceed by determining which hyperfine states to assign to the four levels of the model .
there are two hyperfine ground - state manifolds in @xmath1rb with total angular momentum @xmath127 and @xmath128 separated by 6.835 ghz .
the choice of which ground - states within these manifolds to use is determined by several factors .
firstly , cross - coupling of the control field and the quantum field can lead to unwanted processes such as four - wave mixing ( fwm ) , which is described in detail in app .
[ app : fwm ] . in order to mitigate cross - coupling
, it is advantageous to use ground - states that are widely separated in energy and couple to the excited state through perpendicular polarizations .
this becomes especially important at high optical depth and large control fields  @xcite .
one can choose ground states from the same @xmath129 manifold , but then the energy splitting is limited to what can be obtained by zeeman shifting the levels .
for the purposes of qd photon storage , we choose one ground state in @xmath127 and one in @xmath128 .
furthermore , since both states must be coupled to a common excited state , the difference between the @xmath130 values must be between -2 and 2 .
this sets the first condition for choosing the ground states .    as discussed in sec .
[ sec : intro ] , the long coherence between the hyperfine ground states of @xmath1rb make it very attractive as a quantum memory .
this long coherence is usually limited experimentally by magnetic field fluctuations that cause small zeeman shifts of the @xmath130 levels .
therefore , choosing @xmath130 levels that are common mode to magnetic field fluctuations for the ground states @xmath7 and @xmath8 is crucial to achieving long storage times .
notably , it is not possible to find two states in the same @xmath129 manifold that are common mode to magnetic field fluctuations .
however , there are two non - degenerate choices for the ground states from different @xmath129 manifolds that are common mode to magnetic field fluctuations ( to first order ) and can be coupled to the same excited state .
these levels are the clock states @xmath131 and @xmath132 and the states @xmath133 and @xmath134 ( @xmath135 and @xmath136 is a degenerate choice , but is not magnetically - trappable ) .
for each pair of ground states , there are two optical transitions that are of interest ; the @xmath85 transition to the @xmath137p@xmath120 states near 780  nm and the @xmath121 transition to the @xmath137p@xmath119 states near 795  nm .
we will investigate these transitions separately in the following sections .
the excited states of the @xmath85 transition of @xmath1rb are composed of four hyperfine manifolds @xmath138 0 , 1 , 2 , and 3 . because the two ground states must share a common excited state , only states in @xmath139 and @xmath140 are eligible .
for the pair of ground states @xmath133 and @xmath134 , there are then two possible excited states ; @xmath141 and @xmath142 as shown in blue in fig .
[ fig : fig8 ] .
rb that are eligible for implementation of a quantum memory . ]    in the following , we always work on resonance with one of the excited states ( @xmath143 ) and let the value of @xmath87 change signs .
that is to say , @xmath7 and @xmath8 couple resonantly to @xmath9 and are detuned to @xmath88 .
two choices to assign @xmath7 and two choices for @xmath9 result in four unique configurations for each set of four levels . to determine which configuration will yield the highest storage efficiency , it is instructive to consider the results of sec .
[ sec:3l ] .
because the storage efficiency depends critically on the optical depth @xmath24 , and the effective optical depth on the @xmath11 transition is @xmath144 , it would seem optimal to choose the configuration where this relative dipole moment is the largest .
a secondary concern is that @xmath145 should be large to avoid extremely intense control pulses .
the four relative dipole moments are listed in table  [ tab : d2_1_mu ] ( see app .
a for details ) , where it can be seen that the largest moment is for the @xmath146 transition .
.relative dipole moments @xmath43 for @xmath1rb @xmath85 transitions with ground states @xmath133 and @xmath134  @xcite .
[ cols="^,^,^",options="header " , ]     [ tab : d1_2_eff ]    . ]
configuration 2 yields comparable efficiencies to the best configuration for @xmath85 , with only a minor increase .
if we perform gradient ascent for this configuration but neglect the additional excited state ( configuration 3 ) , we obtain almost identical efficiencies .
comparison of the two optimal control pulses ( fig .
[ fig : fig11 ] ) shows that they are almost exactly the same , except the four level control has a small modulation at the frequency of the excited state splitting .
if we use the three level control on the four level system ( configuration 4 ) , we obtain almost exactly the same efficiencies although the control pulse is much simpler to generate .
this analysis simply implies that the large excited state splitting almost completely nullifies the effects of the additional excited state , as anticipated in sec .
[ sec:4l ] .
consequently , the efficiencies we find in this case match well with those predicted by fig .
[ fig : fig3 ] if @xmath147 is used as the effective optical depth of the storage transition .
up to this point , our analysis has assumed that the photons emitted by a qd are completely indistinguishable , which is an important property in photon - based quantum information processing  @xcite .
however , due to their dynamic solid - state environment the photons produced by excitons in quantum dots may not have fourier transform - limited spectra  @xcite .
this leads to an imperfect photon indistinguishability between subsequently emitted photons . up to now , the highest measured two
photon visibility is @xmath148  @xcite , corresponding to a linewidth that is @xmath149 times larger than the transform - limited linewidth @xmath150 . while it is possible to approach @xmath151 by careful sample selection and resonant excitation , in this section the effect of excess linewidth broadening on the storage efficiency
is calculated .
the excess linewidth of quantum dot generated photons can be caused by two distinct physical processes .
the first , sometimes referred to as spectral wandering , is a slow process that causes changes in the carrier frequency of the photons from shot to shot .
this simply means that subsequent photons have slightly different carrier frequencies due to changes in the qd environment on timescales longer than the spontaneous emission lifetime @xmath79 .
the second , referred to as pure dephasing , is a perturbation of the qd s energy levels on timescales shorter than @xmath79 .
this leads to a time - dependent phase within each shot . in a time - averaged spectral measurement
, these effects can both produce a lorentzian lineshape with a linewidth greater than the fourier limit . here
, we treat these cases separately but show that they cause exactly the same effect on storage and retrieval efficiencies in a quantum memory for a given amount of added linewidth @xmath152 .    in the case of spectral
wandering , the waveform of each photon remains @xmath77 but the carrier frequency @xmath153 is drawn from a probability distribution @xmath154 .
the time - averaged spectrum of these photons is then an integral over all realizations of @xmath153 @xmath155 where @xmath156 is the transform limited spectrum centered about @xmath153 .
if @xmath154 is a lorentzian distribution with linewidth @xmath152 , the time - average spectrum is a lorentzian with linewidth @xmath157 .
the average storage and total memory efficiencies can be calculated similarly by determining how the efficiencies depend on @xmath158 and then integrating over all realizations .
we numerically calculate the storage and total efficiencies as a function of @xmath159 using the optimal control pulse found for configuration 4 of @xmath85 storage with @xmath160 .
all other parameters are kept fixed ( @xmath161 , @xmath102 , @xmath162  ns ) and the results are plotted in fig .
[ fig : fig12 ] .
storage with ground states @xmath132 and @xmath131 as a function of @xmath159 .
all other parameters are kept fixed . ]    using @xmath163 and @xmath164 , the efficiencies as a function of @xmath152 can be obtained by @xmath165 where @xmath166 is taken to be a normalized lorentzian distribution of width @xmath152 centered at @xmath167 . performing this integration results in the plot shown in fig .
[ fig : fig13 ] .
storage with ground states @xmath132 and @xmath131 as a function of @xmath152 resulting from spectral wandering . ]
notably , spectral wandering causing a factor of 2 increase in the total linewidth ( @xmath168 ) roughly leads to a factor of 2 reduction in the storage and total efficiencies .
while the initial drop is steep , the total efficiency remains above 5@xmath0 even if the total linewidth is more than 7 times the transform limit , showing that the storage process is relatively robust to the effects of spectral wandering .
the effects of fast , pure dephasing can also be calculated numerically . in this case , the waveform remains at the carrier frequency @xmath153 but the slowly varying amplitude is modified to @xmath169 which includes a time - dependent phase @xmath170 . in the simplest model of pure dephasing
, @xmath170 is driven by a markovian langevin force @xmath171 characterized by @xmath172 and @xmath173 with diffusion constant @xmath174 .
this force will cause phase diffusion , resulting in a lorentzian time - averaged spectrum with total linewidth @xmath175 .
the effect on the storage and total efficiency can be determined by calculating the efficiencies for several trajectories of @xmath170 and averaging .
storage with ground states @xmath132 and @xmath131 as a function of @xmath176 resulting from pure dephasing ( red , open circles ) .
uncertainties in @xmath177 are given by the standard deviation of the ensemble .
for comparison , results from spectral wandering ( fig .  [ fig : fig13 ] ) are shown as blue lines . ]    using the same configuration as for the spectral wandering investigation , we performed storage and retrieval calculations for 100 phase trajectories for several values of @xmath178 .
the average storage and total efficiencies are plotted as red circles in fig .
[ fig : fig14 ] , with the uncertainties given by the standard deviation . for comparison , results from spectral
wandering ( fig .  [ fig : fig13 ] ) are shown as blue lines .
evidently , both spectral wandering and pure dephasing lead to the same reduction in efficiency for the same amount of @xmath152 , indicating that the efficiencies depend only on the time - averaged spectrum and not on the physical mechanism causing an excess linewidth .
a recent experiment by sparkes _ et al . _ has demonstrated an optical depth of 1000 in an ultracold @xmath1rb gas  @xcite using spatial and temporal dark spots .
we therefore investigate how the efficiencies increase for such high optical depth . if the optical depth is increased to @xmath24=500 in the optimization using the clock ground states ( @xmath131 and @xmath132 ) and configuration 4 of @xmath85 storage , the storage ( total ) efficiency is increased to 51  @xmath0 ( 34  @xmath0 ) . if instead configuration 2 of @xmath121 storage with clock ground states is used , the optimization yields a storage ( total ) efficiency of 82  @xmath0 ( 76  @xmath0 ) .
notably , the efficiencies are not dramatically increased for @xmath85 storage while @xmath121 storage is much more promising .
there are two factors that contribute to the worse performance of @xmath85 ; the small excited - state splitting and the large relative dipole element between the storage state @xmath8 and the unwanted excited state @xmath88 .
naively , one would expect an increase in efficiency according to fig .
[ fig : fig3 ] , but that is only true if the additional excited state does nt play a large role in the storage process
. however , as the optical depth increases so does the required control pulse amplitude and therefore the coupling to unwanted excited states . for the case of @xmath85 storage ,
the extra excited state can not be neglected and limits the achievable efficiencies even for large optical depth . on the other hand , for @xmath121 storage of a 1  ns photon at @xmath179 we find very high efficiencies , comparable to the efficiency of 80  @xmath0 reported in ref .  @xcite for much lower bandwidth photons .
in conclusion , we have calculated the efficiency with which a quantum dot generated single photon can be stored and retrieved from an optically - thick @xmath1rb ensemble .
our calculations take into account the multi - level structure of @xmath1rb by extending the standard three - level model of an atomic ensemble quantum memory . using an optical depth of 150 , the storage ( total )
efficiency can reach 46@xmath0 ( 28@xmath0 ) for a photon resulting from the 1  ns spontaneous excitonic decay in a quantum dot .
importantly , this storage can be performed using control pulses obtained from a diode laser rather than requiring pulse energies only achievable with ultrafast , mode - locked laser sources .
increasing the optical depth , for example by using bose - condensed ensembles , an optical cavity , or advanced trapping techniques can increase the efficiencies to near unity for storage on the @xmath121 transition .
in addition , we have studied the effects of spectral diffusion and pure dephasing of the quantum dot generated photons on the storage efficiency and shown that a factor of 2 increase in the time - averaged photon linewidth roughly leads to a factor of 2 reduction in the efficiency . thus , storage and retrieval of single photons from a quantum dot in an @xmath1rb ensemble is feasible with demonstrated experimental parameters even in the presence of non - ideal properties of the quantum dot .
integration of quantum dot sources with atomic ensemble quantum memories may lead to high - speed quantum networks for communication or distributed computation as well as entanglement between collective atomic degrees of freedom and the spin of an electron or hole confined in the quantum dot .
we thank a.  v. gorshkov , m. fleischhauer , and k. srinivasan for helpful discussions .
m.  t.  r. acknowledges support from a marie curie international incoming fellowship .
this work was also supported by the nccr quantum science and technology .
the equations of motion eq .
[ eqn:4leome]-[eqn:4leoms ] ( and by extension eq .  [ eqn:3leome]-[eqn:3leoms ] )
are derived by considering an ensemble of @xmath16 motionless , four - level atoms interacting with a quantum field and a classical field closely following the treatment of ref .
@xcite , whose notation we also adopt for the most part .
the hamiltonian of this system can be expressed as @xmath180 where @xmath181 and the interaction between the light fields and the atoms in the dipole approximation is @xmath182 here , the operators @xmath183 change the internal state of the @xmath184th atom from @xmath185 to @xmath186 and @xmath187 is the dipole moment of an atom for the @xmath188 transition .
the total electric field @xmath189 is composed of a classical field @xmath190 and a quantum field @xmath191 .
the @xmath192-propagating classical field with polarization orientation @xmath193 can be written as @xmath194},\ ] ] where @xmath195 is the envelope function and @xmath196 is the carrier frequency .
we have assumed that the classical pulse propagates with a group velocity of @xmath197 , which is valid if almost all of the atomic population remains in state @xmath7 .
the quantum field is taken to be a sum of modes @xmath198 centered about frequency @xmath199 with polarization @xmath200 and cross - sectional area @xmath18 @xmath201 where h.c .
denotes the hermitian conjugate . in this treatment ,
the quantum field only drives the @xmath11 and @xmath202 transitions while the classical control field drives the @xmath13 and @xmath203 transitions . as shown in app .
[ app : fwm ] , cross - coupling of the control field leads to a four - wave mixing process which can reduce the storage efficiency , but it is safe to neglect for the parameter regime considered here .
applying this assumption and making the rotating wave approximation allows eq .
[ eq : appa_hint ] to be written as the sum of the interaction with the classical field @xmath204 and the quantum field @xmath205 where @xmath206 are the rabi frequencies associated with the classical field and @xmath207 are the couplings to the quantum field . note that the rabi frequencies @xmath208 are defined differently compared to the standard definition such that @xmath209 .
in order to treat the ensemble as a continuous density distribution , we divide the ensemble into thin slices of thickness @xmath210 such that the quantum field can be taken to be constant over this range while also ensuring that the number of atoms in a slice @xmath2111 .
then , we define slowly varying operators @xmath212 using these operators , we can rewrite @xmath213 as @xmath214 + \int d\omega\ , \hbar \omega \hat{a}^{\dag}_{\omega } \hat{a}_{\omega}\end{aligned}\ ] ] and the two parts of the interaction as @xmath215 \\
\mathcal{h}_{\textrm{int , q } } & = -\hbar \int_0^l dz \ ,
n(z ) \big\{g_{1 g } \hat{\mathcal{e}}(z , t ) \hat{\sigma}_{1g}(z , t ) \nonumber \\ & + g_{2 g } \hat{\mathcal{e}}(z , t ) \hat{\sigma}_{2g}(z , t)+ \textrm{h.c . }
\big\},\end{aligned}\ ] ] where @xmath19 is the linear number density of atoms along the length of the ensemble and we have assumed the cross - sectional area of the field @xmath18 matches that of the atomic cloud .    the dynamics are determined by finding the heisenberg equations of motion for the operators .
using the commutation relations @xmath216 & = \delta(\omega-\omega ' ) \\ \left [ \hat{\sigma}_{\alpha \beta}(z , t),\hat{\sigma}_{\lambda \rho}(z',t ) \right]&= \frac{\delta(z - z')}{n(z ) } \big\{\delta_{\beta \lambda}\hat{\sigma}_{\alpha\rho}-\delta_{\alpha \rho}\hat{\sigma}_{\lambda\beta } \big\}\end{aligned}\ ] ] one obtains @xmath217 \\
\partial_t \hat{\sigma}_{g1 } & = i \delta_g \hat{\sigma}_{g1}+i \omega_{1s } \hat{\sigma}_{gs } \nonumber \\ & + i \hat{\mathcal{e } } \big [ g_{1g}(\hat{\sigma}_{gg}-\hat{\sigma}_{11})-g_{2 g } \hat{\sigma}_{21 } \big]\\ \partial_t \hat{\sigma}_{g2 } & = i ( \delta_g-\delta_e ) \hat{\sigma}_{g2}+i \omega_{2s } \hat{\sigma}_{gs } \nonumber \\ & + i \hat{\mathcal{e } } \left[g_{2g}(\hat{\sigma}_{gg}-\hat{\sigma}_{22})-g_{1 g } \hat{\sigma}_{12 } \right]\\ \partial_t \hat{\sigma}_{s1 } & = i \delta_s \hat{\sigma}_{s1}+i\omega_{1s } ( \hat{\sigma}_{ss}-\hat{\sigma}_{11 } )   \nonumber \\ & - i\omega_{2s}\hat{\sigma}_{21}+i g_{1 g } \hat{\mathcal{e } } \hat{\sigma}_{sg}\\ \partial_t \hat{\sigma}_{s2 } & = i ( \delta_s-\delta_e ) \hat{\sigma}_{s2}+i\omega_{2s } ( \hat{\sigma}_{ss}-\hat{\sigma}_{22 } )   \nonumber \\ & - i\omega_{1s}\hat{\sigma}_{12}+i g_{2 g } \hat{\mathcal{e } } \hat{\sigma}_{sg}\\ \partial_t \hat{\sigma}_{12 } & = -i \delta_e \hat{\sigma}_{12}+i(\omega_{2s}\hat{\sigma}_{1s}-\omega_{s1}\hat{\sigma}_{s2 } ) \nonumber\\ & + i(g_{2 g } \hat{\mathcal{e}}\hat{\sigma}_{1g}-g_{g1}\hat{\mathcal{e}}^{\dag}\hat{\sigma}_{g2})\\ \partial_t \hat{\sigma}_{gs } & = i(\delta_g-\delta_s)\hat{\sigma}_{gs}+i(\omega_{s1}\hat{\sigma}_{g1}+\omega_{s2}\hat{\sigma}_{g2 } ) \nonumber\\ & - i\hat{\mathcal{e}}(g_{1g}\hat{\sigma}_{1s}+g_{2g}\hat{\sigma}_{2s})\\ \partial_t \hat{\sigma}_{gg } & = i ( g_{g1}\hat{\mathcal{e}}^{\dag}\hat{\sigma}_{g1}-g_{1g}\hat{\mathcal{e}}\hat{\sigma}_{1 g } \nonumber \\ & + g_{g2}\hat{\mathcal{e}}^{\dag}\hat{\sigma}_{g2}-g_{2g}\hat{\mathcal{e}}\hat{\sigma}_{2 g } ) \\
\partial_t \hat{\sigma}_{ss } & = i ( \omega_{s1 } \hat{\sigma}_{s1}-\omega_{1s } \hat{\sigma}_{1s } \nonumber \\ & + \omega_{s2 } \hat{\sigma}_{s2}-\omega_{2s } \hat{\sigma}_{2s } ) \\ \partial_t \hat{\sigma}_{11 } & = i ( \omega_{1s}\hat{\sigma}_{1s}-\omega_{s1}\hat{\sigma}_{s1 } ) \nonumber \\ & + i(g_{1g}\hat{\mathcal{e}}\hat{\sigma}_{1g}-g_{g1}\hat{\mathcal{e}}^{\dag}\hat{\sigma}_{g1 } ) \\ \partial_t \hat{\sigma}_{22 } & = i ( \omega_{2s}\hat{\sigma}_{2s}-\omega_{s2}\hat{\sigma}_{s2 } ) \nonumber \\ & + i(g_{2g}\hat{\mathcal{e}}\hat{\sigma}_{2g}-g_{g2}\hat{\mathcal{e}}^{\dag}\hat{\sigma}_{g2})\end{aligned}\ ] ] where the time and spatial dependencies have been neglected for brevity .
these equations can be reduced considerably by making one simplifying assumption ; the quantum field is weak .
the first consequence of this assumption is that almost all atoms remain in @xmath7 for the duration of the dynamics .
secondly , we keep only terms that are linear in @xmath218 @xcite . under these assumptions ,
the equations of motion are reduced to @xmath217 \\
\partial_t \hat{\sigma}_{g1 } & = i \delta_g \hat{\sigma}_{g1}+i\omega_{1s } \hat{\sigma}_{gs}+i g_{1 g } \hat{\mathcal{e } } \\
\partial_t \hat{\sigma}_{g2 } & = i ( \delta_g-\delta_e ) \hat{\sigma}_{g2}+i\omega_{2s } \hat{\sigma}_{gs}+i g_{2 g } \hat{\mathcal{e } } \\
\partial_t \hat{\sigma}_{gs } & = i ( \delta_g-\delta_s)\hat{\sigma}_{gs}+i(\omega_{s1 } \hat{\sigma}_{g1}+\omega_{s2 } \hat{\sigma}_{g2}).\end{aligned}\ ] ] we now introduce the polarization operators @xmath32 and @xmath89 as well as the spin wave operator @xmath33 . in addition , we move to a new coordinate system @xmath39 where @xmath42 is the time in a co - moving reference frame and @xmath40 is a dimensionless length . inserting these definitions into the equations of motion
yields @xmath219 where a factor of @xmath220 has been absorbed into @xmath218 and we introduced the optical depths @xmath221 .
we have also assumed that both polarizations decay at the same rate @xmath22 and that @xmath222 is real .
finally , for notational convenience we compare all transition dipole moments @xmath223 to that of the two - level cycling transition ( @xmath224 ) . in this way
, we make the substitutions @xmath225 and @xmath226 in the equations of motion .
this enables easy comparison between different state configurations and also sets @xmath227 to where it will be measured experimentally .
it also enables @xmath208 to be connected to a light intensity through the two - level relation @xmath228 , where @xmath229 is the saturation intensity . using the value of @xmath229 for the cycling transition  @xcite
, the peak power can be related to the peak rabi frequency @xmath83 by @xmath230w_o^2 ( \omega_m/\gamma)^2 \label{eq : powerrabi}\end{aligned}\ ] ] where @xmath59 is the wavelength of the transition and @xmath231 is the @xmath232 waist of a gaussian beam .
the relevant values for the @xmath1rb cycling transition @xmath233 of the @xmath85 line have been inserted to obtain eq .
[ eq : powerrabi ] .
similarly , the pulse energy can be determined by @xmath234  wm@xmath235w_o^2 \int_0^t d\tau |\omega(\tau)/\gamma|^2 $ ] . from ref .
@xcite , we have @xmath236 where @xmath237 is the reduced dipole moment for the @xmath85 transition .
to put the @xmath121 relative dipole moments @xmath43 in units of @xmath44 , the values in tables  [ tab : d1_1_mu ] and [ tab : d1_2_mu ] have been multiplied by the factor  @xcite @xmath238    further , as discussed in sec .
[ sec:3l ] of the main text , the operators can be treated as complex numbers and their associated quantum noise can be neglected because we are interested in computing the expectation values of normally - ordered operators . making the substitutions and dropping the operator notation yields the equations of motion as found in sec .
[ sec:3l]-[sec:4l ] , @xmath90 \\ \partial_{\tau}p_1 & = \left(i\delta_g-\gamma \right ) p_1+i \tilde{\mu}_{1s } \omega(\tau ) s + i \tilde{\mu}_{1
g } \sqrt{d \gamma } \mathcal{e}\\ \partial_{\tau}p_2 & = \left(i\delta_g -i\delta_e-\gamma\right ) p_2 + i \tilde{\mu}_{2s } \omega(\tau ) s \nonumber \\ &
+ i \tilde{\mu}_{2 g } \sqrt{d \gamma } \mathcal{e } \\ \partial_{\tau}s & = i \left(\delta_g - \delta_s \right ) s + i\omega^*(\tau)\left [ \tilde{\mu}_{1s}p_1 + \tilde{\mu}_{2s}p_2 \right].\end{aligned}\ ] ]
numerical solutions of the equations of motion are obtained using the method of finite differences .
the time - space grid is composed of 9@xmath239 ( 3000 by 3000 ) points and the domain of @xmath240 is chosen such that the optimized control pulses tend toward 0 at @xmath241 .
this ensures that the drop to @xmath242 outside of the domain is smooth .
the gradient ascent algorithm is implemented using a dynamic step size @xmath59 to guarantee quick convergence .
the step size is determined using an inexact line search such that @xmath59 is initialized to a large value ( 1000@xmath22 ) at each step of the ascent and the increase of @xmath243 and its gradient are calculated at the next step .
if this step does not meet the wolfe conditions  @xcite , the step size is reduced geometrically until they are satisfied .
the optimization proceeds until @xmath243 has not been increased by more than 0.001@xmath243 compared to the average of the three previous values .
this tolerance was estimated by considering the errors resulting from the numerical integration .
errors in the reported efficiencies are @xmath2441  @xmath0 which was determined by examining the variation in efficiencies for perturbations of the time - space grid . for the parameter set
@xmath80 , @xmath102 , @xmath245  mhz , the gradient ascent optimization took approximately 20 minutes on a standard computer .
the pure dephasing process was simulated by first obtaining @xmath171 on a finite grid using ( pseudo-)random numbers drawn from a normal distribution .
then , @xmath170 was calculated using @xmath246 .
averaging over many trajectories yielded the correct diffusion behavior @xmath247 ^ 2 \rangle = d_{\phi } \tau$ ] .
the extracted @xmath4 and @xmath243 shown in fig .
[ fig : fig14 ] were obtained by averaging the results of 100 different phase trajectories for each value of @xmath174 .     for the heaviside model and for the three - level ( 3l ) model ( @xmath248  ps ) of the qd photon . ]
the slowly varying photon waveform @xmath77 is not physical due to the infinitely sharp rise of @xmath78 .
a more realistic model might include the non - zero temporal width of the excitation pulse or fast loading of the qd from another excited state .
for the latter , @xmath249 for a three - level model with instantaneous excitation of an ancillary excited state that loads the qd with rate @xmath250 .
for comparison , we have implemented this form of @xmath251 with @xmath252  ps  @xcite using configuration 4 of @xmath85 storage with @xmath253 .
we obtain roughly the same efficiencies ( @xmath254 , @xmath255 ) and peak control rabi frequency ( @xmath256 ) as found in table  [ tab : d2_2_eff ] .
both input photons and control pulses are plotted in fig .
[ fig : figb_1 ] . because the curves and efficiencies are extremely similar
, we conclude that the infinitely sharp rise of @xmath78 in the simple photon waveform does not dramatically influence the results and is therefore a sufficiently representative choice .
in the preceding analysis , cross - coupling of the control field @xmath12 to the @xmath11 transition was neglected due to the large ground - state hyperfine splitting of @xmath2576.835  ghz . because the control field can be quite strong and the optical depth quite large
, this cross - coupling can lead to detrimental effects that reduce the storage and retrieval efficiencies in practice  @xcite .
of course , a proper choice of ground states and optical polarizations can eliminate cross - coupling completely .
for example , choosing ground states whose @xmath130 values differ like @xmath258 and @xmath259 and using circularly - polarized light allows cross - coupling to be neglected . on the other hand ,
configurations such as @xmath260 and @xmath261 are coupled to the excited state by light of the same polarization . in this case , the control beam can off - resonantly drive the @xmath11 transition as shown in fig .  [
fig : figc_1 ] and coherently generate a stokes field @xmath262 in a four - wave mixing ( fwm ) type of process  @xcite .
the stokes field can interfere with the spin - wave created from storage of the quantum field @xmath35 and lead to reduced storage and retrieval efficiencies .    in ref .
@xcite , the stokes field was taken into account theoretically and was shown to match experimental results quite well .
following the same approach here , we obtain the stokes - modified equations of motion for a three - level system @xmath263 where @xmath262 is the stokes field and the off - resonant interaction @xmath264 has been adiabatically eliminated , leaving an effective coupling between @xmath262 and @xmath37 .
in addition , this interaction induces time - dependent light shifts of + @xmath265 and -@xmath266 for states @xmath9 and @xmath7 respectively . using these equations of motion with the optimized control field for on - resonance storage of a qd - generated photon with configuration 3 of @xmath121 storage with ground states @xmath131 and @xmath132 ( fig .
[ fig : fig11 ] ) , one obtains the same storage and total efficiencies as in table  [ tab : d1_2_eff ] .
this result indicates that the effect of fwm in this storage scheme is negligible .
another way to determine the relative effect of fwm is to consider the ratio of the last two terms in eq .
[ eqn:3lfwms ] .
the ratio of the fwm term to the normal @xmath267 term should roughly scale as @xmath268 , which is approximately 2@xmath269 for typical parameters considered here .
compared to typical parameters found in ref .
@xcite where the onset of fwm effects was measured , this ratio is about three orders of magnitude smaller so fwm can be safely neglected .
60ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]
+ 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty @noop _
( ,  ,  ) link:\doibase    10.1038/nature08812 [ * * ,   ( ) ] @noop * * ,   ( ) link:\doibase 10.1103/physrevlett.81.5932 [ * * ,   ( ) ] link:\doibase 10.1038/nature07127 [ * * , ( ) ] link:\doibase 10.1103/revmodphys.83.33 [ * * ,   ( ) ] link:\doibase 10.1038/35106500 [ * * ,   ( ) ] link:\doibase 10.1038/nphoton.2009.231 [ * * ,   ( ) ] link:\doibase 10.1038/nature04315 [ * * ,   ( ) ] link:\doibase 10.1038/nature06670 [ * * ,   ( ) ] link:\doibase    10.1038/nphys1153 [ * * ,   ( ) ] link:\doibase 10.1038/nphys1152 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.98.183601 [ * * ,   ( ) ] link:\doibase 10.1038/nature04327 [ * * , ( ) ] link:\doibase 10.1038/nature05136 [ * * ,   ( ) ] link:\doibase    10.1038/nphoton.2010.30 [ * * ,   ( ) ] link:\doibase 10.1038/nature09997 [ * * ,   ( ) ] link:\doibase    10.1038/nature09662 [ * * ,   ( ) ] link:\doibase    10.1038/nature09719 [ * * ,   ( ) ] link:\doibase 10.1038/nphys2324 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.76.050301 [ * * ,   ( ) ] @noop * * ,   ( ) link:\doibase 10.1038/nphoton.2007.46 [ * * ,   ( ) ] link:\doibase 10.1038/nphoton.2007.227 [ * * ,   ( ) ] link:\doibase    10.1038/nphoton.2009.287 [ * * ,   ( ) ] link:\doibase 10.1063/1.3615051 [ * * ,   ( ) ] link:\doibase    10.1038/nature04446 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.103.217402 [ * * , ( ) ] link:\doibase    10.1103/physrevlett.103.167402 [ * * ,   ( ) ] link:\doibase 10.1038/nnano.2012.262 [ * * ,   ( ) ] link:\doibase 10.1126/science.1126074 [ * * ,   ( ) ] link:\doibase    10.1038/nature06472 [ * * ,   ( ) ] link:\doibase 10.1038/nphys2078 [ * * ,   ( ) ] link:\doibase 10.1126/science.1173684 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.75.081302 [ * * ,   ( ) ] link:\doibase 10.1038/nature11577 [ * * ,   ( ) ] link:\doibase 10.1038/nature11573 [ * * ,   ( ) ] @noop * * ,   ( ) link:\doibase 10.1103/physreva.76.033805 [ * * ,   ( ) ]
link:\doibase 10.1103/physreva.77.043806 [ * * ,   ( ) ] link:\doibase 10.1103/revmodphys.75.457 [ * * , ( ) ] link:\doibase    10.1103/physreva.75.011401 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.76.033804 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.80.062330 [ * * ,   ( ) ] link:\doibase    10.1038/nphoton.2010.221 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.107.083602 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.109.147405 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.66.153316 [ * * ,   ( ) ] @noop * * , ( ) link:\doibase    10.1063/1.3478232 [ * * ,   ( ) ] link:\doibase    10.1038/nphoton.2011.16 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.83.063823 [ * * ,   ( ) ] @noop link:\doibase    10.1103/physrevlett.105.260502 [ * * ,   ( ) ] @noop * * ,   ( ) @noop ( ) ,   @noop _ _  ( ,  ,  ) @noop * * ,   ( ) link:\doibase    10.1103/physreva.79.013806 [ * * ,   ( ) ] link:\doibase 10.1080/09500340903159511 [ * * ,   ( ) ]