nonlinearities are important in optical systems since they enable us to manipulate light by light .
especially , nonlinearities sensitive to single photons are necessary for quantum optical applications such as optical quantum computation @xcite , entanglement generation @xcite , and quantum non - demolition measurements @xcite . in particular , strong nonlinearities are indispensable for quantum computation , since conditional operations at the single photon level are required to realize elementary gate operations .
however , conventional optical media require very high photon densities to obtain significant nonlinear effects @xcite .
single atom cavity quantumelectrodynamics may offer a possible solution @xcite .
in particular , a cavity can focus all the light in an input beam on a single atom , effectively creating a situation where a single atom is coupled to a one - dimensional light field @xcite .
such a one - dimensional atom is an ideal device for the non - linear manipulation of few photon pulses .
it is therefore of great interest to investigate the response of a one - dimensional atom to one- and two - photon input pulses .
there are many theories for such atom - cavity systems @xcite .
however these theories usually eliminate the quantum state of the field outside the atom - cavity system and therefore the analysis of the precise spatiotemporal coherence of the input and output photon pulses can not be achieved . in our previous work ,
we have therefore presented a theory for one- and two - photon pulses that includes the propagation to and from the system in a bad - cavity regime @xcite , based on a one - dimensional model of light - atom interaction @xcite . in these studies
, we focused on the changes of the photon - photon correlations in time @xcite and in frequency @xcite .    in this paper
, we study the effect of the non - linear interaction with the one - dimensional atom on the shape of one- and two - photon wavefunctions using gaussian input pulses of different pulse lengths . for this purpose
, we first introduce an analysis of the interaction in terms of absorption and transmission processes in section [ sec : level2 ] . the characteristic time scale defining the response of the atom
is given by the dipole relaxation time @xmath0 . in section [ sec : level3 ] , we then analyze the interaction processes for a long pulse length of @xmath1 in terms of pulse amplitude and pulse delay time . from these results
, we can estimate that a maximal non - linear effect should be obtained for a pulse length equal to the dipole relaxation time of the atom . in section [ sec :
level4 ] , we therefore analyze the corresponding interaction processes for a short pulse length of @xmath0 , confirming the dominance of the non - linear effect throughout the short two - photon pulse . the results of the analysis are summarized and discussed in section [ sec : level5 ] .
the model of a single two - level atom in one - dimensional free space is shown in figure .
[ fig : figure1](a ) . in this model ,
an incoming one - photon or two - photon wavepacket interacts with the atom locally , and the output wavepacket is then emitted into the output field .
the physical realization of this model can be implemented by a two - level atom coupled with a single mode of a one - sided cavity in the bad cavity regime @xcite .
the cavity geometry is shown in figure .
[ fig : figure1 ] ( b ) .
the input field of the one - dimensional free space shown in figure .
[ fig : figure1 ] ( a ) corresponds to the input of the one - sided cavity in figure .  [ fig : figure1 ] ( b ) . likewise , the output field in figure .
[ fig : figure1 ] ( a ) corresponds to the output of the one - sided cavity in figure .
[ fig : figure1 ] ( b ) . here
, we assume that the spontaneous emission rate through the cavity mode is much larger than the spontaneous emission rate through non - cavity modes . in terms of the conventional cavity
quantum electrodynamics parameters , this regime is characterized by @xmath2 , where @xmath3 is the cavity damping rate through the left mirror , @xmath4 is the dipole coupling between the atom and the cavity mode and @xmath5 is the rate of spontaneous emission into the non - cavity modes .
since the cavity damping rate @xmath3 is much faster than the dipole coupling @xmath4 , the method of adiabatic elimination can be applied to the time evolution of the cavity field @xcite .
this means that the interaction between the atom and the outside field through the cavity field can be expressed by an effective dipole relaxation rate @xmath6 . the dipole relaxation rate @xmath7 describes the dipole damping caused by emissions through the left mirror of the cavity , and the corresponding rate of spontaneous emission through the cavity is equal to @xmath8 @xcite .
in our case , we assume that the rate of spontaneous emission into the non - cavity modes @xmath5 is negligible ( @xmath9 ) .
nearly all emissions from the atom can then be confined to the cavity and 2@xmath7 is the total spontaneous emission rate of the excited atom in the cavity . in present cavity designs , this can be realized by exploiting the purcell effect .
for example , in the case of turchette et al.s experiment @xcite , the cavity parameters indicate that about 80% of the sponaneous emission from the atom is emitted through the single cavity mode .
( 0,0 ) ( -10,160)(a ) ( 270,160)(b )    ( a ) schematic of a two - level atom in one - dimensional free space and ( b ) physical realization of the one - dimensional atom using an atom - cavity system .
an incoming one - photon or two - photon wavepacket interacts with the atom locally , and the output wavepacket is then emitted into the output field .
the input field of the one - dimensional free space shown in ( a ) corresponds to the input of the one - sided cavity in ( b ) .
likewise , the output field of the one - dimensional free space shown in ( a ) corresponds to the output of the one - sided cavity in ( b).,title="fig:",width=340 ]   ( a ) schematic of a two - level atom in one - dimensional free space and ( b ) physical realization of the one - dimensional atom using an atom - cavity system .
an incoming one - photon or two - photon wavepacket interacts with the atom locally , and the output wavepacket is then emitted into the output field .
the input field of the one - dimensional free space shown in ( a ) corresponds to the input of the one - sided cavity in ( b ) .
likewise , the output field of the one - dimensional free space shown in ( a ) corresponds to the output of the one - sided cavity in ( b).,title="fig:",width=245 ]      in the following , we investigate the effect of the one - dimensional atom on one and two - photon input pulses of different pulse length . for this purpose
, we have to formulate the quantum states for a given input pulse shape for both one photon and two photon cases .
if the real space representation for the one - photon state reads , @xmath10 the two - photon state can be expressed as the product state @xmath11 the light pulses characterized by the above input states are propagating at the velocity of light in the input field of the schematic shown in figure .
[ fig : figure1](a ) . to simplify the expression , the spatial coordinates , @xmath12 and @xmath13 ,
therefore represent a coordinate system moving at the velocity of light .
the indices 1 and 2 of @xmath12 and @xmath13 in eq .
( [ eq : spatialinput ] ) identify the two particles .
the spatial features of the probability amplitudes given by the wavefunction @xmath14 correspond to the spatial features of the pulse mode a. thus , the wavefunction @xmath15 characterizes the single mode two - photon wavefunction , where the wavefunctions of both particles overlap perfectly . in our previous work @xcite
, we derived the unitary time evolution operator in the hilbert spaces for the one photon component and the two photon component by solving the schrdinger equation of the interaction between a one - dimensional field and a single atom .
the output one - photon state in the far field can then be obtained by convoluting the one - photon input wavefunction @xmath16 with the matrix description of the unitary operator in real space , @xmath17 @xmath18    the matrix elements @xmath19 are given by @xmath20 likewise , the output two - photon state in the far field can be obtained by a linear convolution , @xmath21 @xmath22 the matrix elements @xmath23 are given by @xmath24 where the contribution of the non - linear photon - photon interaction is given by @xmath25 .
\label{eq : nonlineffect}\end{aligned}\ ] ] the output wavefunctions @xmath26 and @xmath27 describe the spatial feature of the output states in the far field of the atom shown in figure .
[ fig : figure1 ] ( a ) .
note that the spatiotemporal features of the input one - photon and two - photon wavefunctions are generally not preserved in the output due to the atomic response described by the matrix elements @xmath28 and @xmath29 .
the matrix element of the time evolution @xmath28 in eq.([eq : onemat ] ) can be expanded in terms of two interaction processes , @xmath30 where the two processes are ( i ) single photon transmission without absorption , given by @xmath31 and ( ii ) single photon reemission after absorption , given by @xmath32 the output one - photon wavefunction @xmath26 can thus be interpreted as the result of quantum interference of the two processes , @xmath33 where the components are given by @xmath34    in the same way , the matrix element of the time evolution @xmath35 can be expanded in terms of three interaction processes .
process ( 1 ) is the transmission of both photons without absorption .
process ( 2 ) is the transmission of one photon without absorption and the absorption and reemission of the other photon .
process ( 3 ) is the absorption and reemission of both photons . in the following
, we will refer to process ( 1 ) as the transmission process , to process ( 2 ) as the one photon absorption process , and to process ( 3 ) as the two photon absorption process .
the transformation of the two - photon wavefunction can then be written as @xmath36 the components of @xmath35 read @xmath37 note that the non - linear contribution @xmath38 only occurs in the two photon absorption process .
this is because the non - linearity is due to the saturation of the two - level atom , which prevents the absorption of two photons at the same time .
likewise , the output wavefunction @xmath39 can be expanded as @xmath40 where the @xmath41 are given by @xmath42 this analysis will be convenient to understand the origin of the statistical properties of the output photons . in the next section ,
these results are applied to the cases of gaussian one - photon and two - photon input wavepackets .
the dependence of the coherent spatiotemporal correlations between the output photons on the input pulse length is analyzed numerically using the results in the previous subsection . for our analysis
, we use gaussian one - photon and two - photon input wavefunctions , @xmath43 the pulse length @xmath44 is defined as twice the dispersion of the gaussian distribution of the photons in time .
note that the response of the atom is governed by the time scale @xmath0 , therefore the atom responds differently to long pulses and short pulses . in the following
, the term long pulse refers to a pulse whose pulse length @xmath44 is much longer than the dipole relaxation time @xmath0 .
likewise , the term short pulse refers to a pulse whose pulse length @xmath44 is about equal to the dipole relaxation time @xmath0 . in the next section ,
we analyze the output wavefunctions for long one - photon and two - photon pulses in order to identify the characteristic quantum interference effects .
[ fig : figure2 ] ( a ) shows an example of a one - photon long pulse ( broken line ) and the corresponding output wavefunction ( solid line ) . in this example , we have chosen an input pulse length of @xmath45 which is 10 times greater than the dipole relaxation time @xmath0 .
[ t][t]@xmath46 [ b][b]@xmath47 [ b][b]@xmath48    ( 0,0 ) ( 0,200)(a )    gaussian input and output one - photon wave packets for an input wave packet length of @xmath45 . the horizontal axes @xmath49 represents the spatial coordinate of one photon . in ( a )
, the broken line corresponds to the input wavefunction @xmath16 and the solid line corresponds to the output wavefunction @xmath26 . in ( b ) , the dotted line corresponds to the component of the process ( i ) : single photon transmission without absorption and the thin line corresponds to the component of the process ( ii ) : single photon reemission after absorption . ]       ( 0,0 ) ( 0,200)(b )    gaussian input and output one - photon wave packets for an input wave packet length of @xmath45 .
the horizontal axes @xmath49 represents the spatial coordinate of one photon . in ( a )
, the broken line corresponds to the input wavefunction @xmath16 and the solid line corresponds to the output wavefunction @xmath26 . in ( b ) , the dotted line corresponds to the component of the process ( i ) : single photon transmission without absorption and the thin line corresponds to the component of the process ( ii ) : single photon reemission after absorption . ]
the pulse shape of the output one - photon pulse is almost identical with the input pulse except for the @xmath50 phase flip in the output amplitudes compared to the input amplitudes , @xmath51 .
furthermore the minimum of the output amplitude is delayed compared to the maximum of the input amplitudes .
the delay time is about @xmath52 .
let us now analyze the origin of these features .    figure .
[ fig : figure2 ] ( b ) shows the components corresponding to different light - atom interaction processes as given by eqs.([eq : oneoutputcomponents ] ) .
the dotted line is the component of single photon transmission without absorption , @xmath53 ( see eqs.([eq : oneoutputcomponents ] ) ) .
the thin line is the component of single photon absorption , @xmath54 . since the absorption process convolutes the input amplitude with an exponential function according to eqs .
( [ eq : oneoutputcomponents ] ) , this component is not a precise gaussian anymore .
however , in the long pulse limit , the convolution function is so much narrower than the gaussian that the change in pulse shape is negligibly small .
the component @xmath54 can then be approximated as @xmath55 as this approximation indicates , the amplitude of the absorption process is about two times as large as the amplitude of the transmission process .
it is then possible to understand the @xmath50 phase flip in the output as a result of interference , @xmath56 where the negative amplitude of the absorption process is twice as high as the positive amplitude of the transmission process .
however this approximation can not explain the delay of the peak of the output wavefunction in figure .  2 ( a ) .
in the long pulse limit , the delay time of @xmath0 in @xmath54 can be included in the approximation given in eq.([eq : appabs ] ) by expanding @xmath54 as follows , @xmath57 the delay time of @xmath26 can then be obtained from the interference of @xmath58 and @xmath54 , @xmath59 the delay time of the output wavefunction is therefore equal to @xmath52 , as seen in figure .
[ fig : figure2 ] ( a ) .
let us now consider the output wavefunction for the long two photon pulse .
the output wavefunction can be analyzed in the same way as the one photon case , using the three interaction processes discussed in section [ sec : level2].3 . however , there are now two photon coordinates , @xmath12 and @xmath13 , and , as shall be seen in the following , the distance @xmath60 between the two photons significantly influences the two photon absorption process .
note also that the bosonic symmetry of the indistinguishable photons requires that @xmath61 .    figure .
[ fig : figure3 ] ( a ) shows the contour plot of a long two - photon input pulse @xmath15 .
the probability amplitude increases from black to white shading . as before
, the input pulse has a length of @xmath45 . figure .
[ fig : figure3 ] ( b ) shows the cross - section of the input wavefunction at @xmath62 for time differences of @xmath63 between the two photons .
note that the horizontal axis gives the average position of the two photons .
the position of the peaks then defines the average delay of both photons at a given distance of @xmath64 between the two photons .
[ t][t]@xmath65 [ b][b]@xmath66 [ t][t]@xmath67 [ ] [ ] @xmath68    [ tr][tr]20 [ tr][tr]10 [ tr][tr]0 [ tr][tr]-10 [ tr][tr]-20    ( 0,0 ) ( -5,200)(a )    input and output wavefunction for a two - photon input pulse with a pulse length of @xmath45 .
( a ) is a contour plot of the input wavefunction @xmath15 .
( b ) is the cross section of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath70 ( broken line ) and @xmath71 ( dotted line ) . the horizontal axis @xmath72 in ( b ) represents the average position of the two photons .
likewise , ( c ) is a contour plot of the output wavefunction @xmath39 and ( d ) shows cross sections of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath70 ( broken line ) , and @xmath71 ( dotted line ) . ]       ( 0,0 ) ( -20,200)(b )    input and output wavefunction for a two - photon input pulse with a pulse length of @xmath45 .
( a ) is a contour plot of the input wavefunction @xmath15 .
( b ) is the cross section of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath70 ( broken line ) and @xmath71 ( dotted line ) . the horizontal axis @xmath72 in ( b ) represents the average position of the two photons .
likewise , ( c ) is a contour plot of the output wavefunction @xmath39 and ( d ) shows cross sections of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath70 ( broken line ) , and @xmath71 ( dotted line ) . ]
[ t][t]@xmath67 [ ] [ ] @xmath73    [ tr][tr]20 [ tr][tr]10 [ tr][tr]0 [ tr][tr]-10 [ tr][tr]-20    ( 0,0 ) ( -5,200)(c )    input and output wavefunction for a two - photon input pulse with a pulse length of @xmath45 .
( a ) is a contour plot of the input wavefunction @xmath15 .
( b ) is the cross section of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath70 ( broken line ) and @xmath71 ( dotted line ) . the horizontal axis @xmath72 in ( b ) represents the average position of the two photons .
likewise , ( c ) is a contour plot of the output wavefunction @xmath39 and ( d ) shows cross sections of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath70 ( broken line ) , and @xmath71 ( dotted line ) . ]       ( 0,0 ) ( -20,200)(d )    input and output wavefunction for a two - photon input pulse with a pulse length of @xmath45 .
( a ) is a contour plot of the input wavefunction @xmath15 .
( b ) is the cross section of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath70 ( broken line ) and @xmath71 ( dotted line ) . the horizontal axis @xmath72 in ( b ) represents the average position of the two photons .
likewise , ( c ) is a contour plot of the output wavefunction @xmath39 and ( d ) shows cross sections of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath70 ( broken line ) , and @xmath71 ( dotted line ) . ]    ( 0,0 ) ( -5,35)(e ) ( -5,-15)-0.35 ( 190,-15)0.3 ( 110,-15)0    input and output wavefunction for a two - photon input pulse with a pulse length of @xmath45 .
( a ) is a contour plot of the input wavefunction @xmath15 .
( b ) is the cross section of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath70 ( broken line ) and @xmath71 ( dotted line ) . the horizontal axis @xmath72 in ( b ) represents the average position of the two photons . likewise , ( c ) is a contour plot of the output wavefunction @xmath39 and ( d ) shows cross sections of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath70 ( broken line ) , and @xmath71 ( dotted line).,height=37 ]    figure .  [
fig : figure3 ] ( c ) and figure .
[ fig : figure3 ] ( d ) show the output wavefunction @xmath39 .
a remarkable new feature of this output wavefunction is the negative amplitude region around @xmath74 .
the output amplitude changes from negative to positive values at about @xmath75 . when going further away from @xmath74 , the spatial features of the output wavefunction eventually become similar to that of the input wavefunction .    in order to illustrate these spatial features in detail , figure .
[ fig : figure3 ] ( d ) shows the cross - sections of the output wavefunction at @xmath62 for @xmath76 .
the negative amplitude of the cross section at a time difference of @xmath69 ( solid line ) is about 3 times as large as the amplitude of the input wavefunction at @xmath69 shown in figure .
[ fig : figure3 ] ( b ) ( solid line ) .
moreover , the peak of the output at @xmath69 is delayed by about @xmath77 compared to the peak of the input wavefunction at the same time difference .
the cross - section at a time difference of @xmath70 ( broken line ) is nearly zero .
this means that the probability density of detecting one photon at a time difference of @xmath70 after detecting the other photon is vanishingly small at all times .
the cross - section at @xmath71 ( dotted line ) is similar to the cross - section of the input wavefunction at the same time difference between the photons .
however the output at @xmath71 is delayed compared to the input wavefunction by about @xmath52 , which is equal to the delay time observed in the one photon case .
let us now consider why such features appear
. the two photon case can be understood as a superposition of three processes in the same way that the one photon case can be understood as a superposition of two processes .
figure.[fig : figure4 ] shows the contour plots of the three components given by eq.([eq : outputcomponents ] ) .
[ fig : figure4 ] ( a ) corresponds to the transmission process , figure .
[ fig : figure4 ] ( b ) corresponds to the one photon absorption process , and figure .
[ fig : figure4 ] ( c ) corresponds to the two photon absorption process , respectively .
the component of the transmission process ( figure .
[ fig : figure4 ] ( a ) ) is identical with the gaussian input wavefunction , just as in the one photon case , @xmath78 . likewise , the component of the one photon absorption process ( figure .
[ fig : figure4 ] ( b ) ) , although not a precise gaussian , can be approximated by the products of single photon transmission and single photon absorption , @xmath79 this approximation indicates that the negative amplitude of the one photon absorption process in figure .  [ fig : figure4 ] ( b ) is about 4 times as large as the positive amplitude of the transmission process in figure .
[ fig : figure4 ] ( a ) .
moreover , the wavepacket of the one photon absorption process is delayed by @xmath80 compared to the transmission process .    the component of the two photon absorption process ( figure .
[ fig : figure4 ] ( c ) ) includes the non - linear interaction between the two photons . for @xmath81
, it resembles the shifted gaussian expected for separate one photon absorption processes ( see eq.([eq : appabs2 ] ) ) , @xmath82 however , the saturation of the atom completely suppresses this process for @xmath83 .
closer inspection shows that the transition between @xmath83 and @xmath81 can be approximated by an exponential function of @xmath60 , @xmath84 the combination of the two photon absorption process is therefore directly dependent on the time difference @xmath85 between the two photons .
[ t][t]@xmath65 [ b][b]@xmath66 [ t][t]@xmath12 [ t][t]@xmath12 [ tr][tr]20 [ tr][tr]10 [ tr][tr]0 [ tr][tr]-10 [ tr][tr]-20    ( 0,0 ) ( -5,200)(a )     contour plots of the two photon processes @xmath86 .
( a ) corresponds to the transmission process , ( b ) corresponds to the one photon absorption process , ( c ) corresponds to the two photon absorption process . ]       ( 0,0 ) ( -5,200)(b )     contour plots of the two photon processes @xmath86 .
( a ) corresponds to the transmission process , ( b ) corresponds to the one photon absorption process , ( c ) corresponds to the two photon absorption process . ]    ( 0,0 ) ( -5,200)(c )     contour plots of the two photon processes @xmath86 .
( a ) corresponds to the transmission process , ( b ) corresponds to the one photon absorption process , ( c ) corresponds to the two photon absorption process . ]
+    ( 0,0 ) ( -10,35)(d ) ( -10,-10)-0.35 ( 190,-10)0.3 ( 110,-10)0     contour plots of the two photon processes @xmath86 . (
a ) corresponds to the transmission process , ( b ) corresponds to the one photon absorption process , ( c ) corresponds to the two photon absorption process.,height=37 ]    in order to illustrate the dependence of the output component on the relative distance between two photons in detail , figure .
[ fig : figure5 ] shows the cross - sections at @xmath62 of the output amplitude @xmath39 ( solid lines ) and the components of the transmission process @xmath87 ( dotted lines ) , the one photon absorption process @xmath88 ( thin lines ) , and the two photon absorption process @xmath89 ( broken lines ) .
the time differences are @xmath90 .
[ t][t]@xmath67 [ t][t]@xmath67 [ t][t]@xmath67 [ t][t]@xmath67 [ b][b ] [ b][b ] [ b][b ] [ b][b ]    ( 0,0 ) ( -5,215)(a )     cross - sections of the output amplitude @xmath39 and its components @xmath86 at time differences @xmath85 of ( a ) @xmath69 , ( b ) @xmath70 and ( c ) @xmath71 .
the solid lines correspond to the total output amplitude @xmath39 , the dotted lines correspond to the transmission amplitude @xmath87 , the thin lines correspond to the one photon absorption amplitude @xmath88 , and the broken lines correspond to the two photon absorption amplitude @xmath89 . ]
( 0,0 ) ( -5,215)(b )     cross - sections of the output amplitude @xmath39 and its components @xmath86 at time differences @xmath85 of ( a ) @xmath69 , ( b ) @xmath70 and ( c ) @xmath71 .
the solid lines correspond to the total output amplitude @xmath39 , the dotted lines correspond to the transmission amplitude @xmath87 , the thin lines correspond to the one photon absorption amplitude @xmath88 , and the broken lines correspond to the two photon absorption amplitude @xmath89 . ]
( 0,0 ) ( -5,215)(c )     cross - sections of the output amplitude @xmath39 and its components @xmath86 at time differences @xmath85 of ( a ) @xmath69 , ( b ) @xmath70 and ( c ) @xmath71 .
the solid lines correspond to the total output amplitude @xmath39 , the dotted lines correspond to the transmission amplitude @xmath87 , the thin lines correspond to the one photon absorption amplitude @xmath88 , and the broken lines correspond to the two photon absorption amplitude @xmath89 . ]    for all three time differences @xmath91 , the amplitude of the component of the one photon absorption process ( thin lines ) is about 4 times as large as the amplitude of the corresponding component of the transmission process ( dotted lines ) . on the other hand , the ratio of the amplitude of the component of the two photon absorption process ( broken lines ) to the amplitude of the corresponding component of the transmission process ( dotted lines ) changes depending on the time difference @xmath91 between the two photons as discussed in the following .    figure .
[ fig : figure5 ] ( a ) shows the cross sections of the output components at @xmath69 . in this case , the amplitude of the component of the two photon absorption process , @xmath92 , is exactly zero ( see eq.([eq : appprocess3 ] ) ) .
the output amplitude is then given by the sum of the transmission amplitude @xmath78 and the one photon absorption amplitude @xmath93 , @xmath94 the amplitude of the output at the time difference of @xmath69 ( solid line in figure .
[ fig : figure5 ] ( a ) ) is therefore about 3 times as large as the corresponding amplitude of the input .    figure .
[ fig : figure5 ] ( b ) shows the cross sections of the output components at @xmath70 . at this special delay time , the factor of @xmath95 describing the effect of the nonlinearity in eq.([eq : appprocess3 ] ) is approximately equal to @xmath96 .
therefore , the two photon absorption component in the total output can be approximated by @xmath97 .
the output wavefunction is then given by @xmath98 at a delay time of @xmath70 , the positive amplitudes of the transmission process and the two photon absorption process thus compensate the negative amplitude of the one photon absorption process , resulting in a total output amplitude of nearly zero .    figure .
[ fig : figure5 ] ( c ) shows the cross sections of the output components at @xmath71 . in this case
, the amplitude of the component of the two photon absorption process , @xmath92 , can be approximated as @xmath99 since the time difference @xmath91 is sufficiently larger than @xmath0 .
the output amplitude is then given by @xmath100 the amplitude of the output at the time difference of @xmath71 ( solid line in figure .
[ fig : figure5 ] ( c ) ) is thus almost identical with the amplitude of the input at the same time difference .    to explain the delay of the output wavefunction
, we can include the delay times in the arguments of the approximations given by eq.([eq : appprocess2 ] ) and eq.([eq : appprocess3 ] ) . for @xmath69 , the output wavefunction
can then be approximated as @xmath101 the amplitude of the output wavefunction at a time difference of @xmath69 is thus delayed by about @xmath77 compared to the amplitude of the input wavefunction at the same time difference . in the case of @xmath70 ,
the output wavefunction is nearly zero , and no delay time can be defined . for @xmath71 , the two photon component @xmath89 of the output wavefunction can be approximated using the limit of @xmath81 given in eq.([eq : shiftedgaussian ] ) . since in this limit
, the two photons are absorbed and reemitted independently of each other , the output wavefunction can then be factorized into one photon processes ( see eq.([eq : apponeoutput ] ) ) , @xmath102 the delay of the peak of the output wavefunction for @xmath81 is thus equal to the delay for the one photon output wavefunction discussed in section [ sec : level3].1 . in summary
, the delay of the negative amplitude at @xmath74 is equal to @xmath103 of the delay of the positive amplitude at @xmath81 .
this means that the effect of the non - linear interaction between the photons at @xmath74 reduces the total delay of the output wavefunction as compared to the delay caused by the uncorrelated transmission and absorption processes at @xmath81 .
as shown in the long pulse case , the atomic nonlinearity is characterized by a negative amplitude of the output wavefunction around @xmath74 . in the long pulse case ,
this nonlinear effect is limited to a region of @xmath104 . in order to realize a non - linear phase change throughout the pulse , it is necessary to shorten the pulse to a pulse length comparable to @xmath0 . in the following ,
we therefore discuss the response of the atom to a pulse of length @xmath105 .      the line plot in figure .
[ fig : figure6 ] ( a ) shows the short one - photon pulse ( broken line ) and the corresponding output wavefunction ( solid line ) .
the length of the input wavepacket is @xmath105 .
[ t][t]@xmath46 [ b][b]@xmath106 [ b][b]@xmath48    ( 0,0 ) ( -5,200)(a )    gaussian input and output one - photon wavepackets for an input wavepacket length of @xmath105 . in ( a )
, the broken line corresponds to the input wavefunction @xmath16 and the solid line corresponds to the output wavefunction @xmath107 . in ( b ) , the dotted line corresponds to the component of the process ( i ) : single photon transmission without absorption and the thin line corresponds to the component of the process ( ii ) : single photon reemission after absorption . ]       ( 0,0 ) ( -5,200)(b )    gaussian input and output one - photon wavepackets for an input wavepacket length of @xmath105 . in ( a )
, the broken line corresponds to the input wavefunction @xmath16 and the solid line corresponds to the output wavefunction @xmath107 . in ( b ) , the dotted line corresponds to the component of the process ( i ) : single photon transmission without absorption and the thin line corresponds to the component of the process ( ii ) : single photon reemission after absorption . ]
the pulse shape of the output pulse is now quite different from the input pulse . specifically , the output amplitude is still positive for @xmath108 and only changes to negative values in the latter part of the pulse .
let us now analyze the origin of these features .    figure .
[ fig : figure6 ] ( b ) shows the components corresponding to different light - atom interaction processes as given by eqs.([eq : oneoutputcomponents ] ) .
the dotted line is the component of single photon transmission without absorption , @xmath58 .
the thin line is that of single photon reemission after absorption , @xmath54 .
as before , @xmath109 according to eqs.([eq : oneoutputcomponents ] ) . on the other hand , the absorption component @xmath54 ( thin line ) is given by the convolution with the exponential function of the dipole relaxation , given by eq.([eq : absgenerater ] ) and eqs.([eq : oneoutputcomponents ] ) .
the peak of this component is shifted by about @xmath110 compared to the component of the transmission process .
furthermore , the peak of the amplitude is only about 1.3 times as large as the corresponding amplitude of the component of the transmission process .
therefore , the front part of the output wavefunction is dominated by the positive component of the transmission process @xmath58 , even though the major part of the output wavefunction is still dominated by the negative amplitude component of the absorption process @xmath54 . the large delay between the components and the reduced amplitude of the dominant negative absorption component
thus cause the asymmetry of the output wavefunction seen in figure .
[ fig : figure6 ] ( a ) .
let us now consider the output wavefunction for the short two photon pulse . figure .
[ fig : figure7 ] ( a ) shows the contour plot of the short two - photon input pulse @xmath15 .
the probability amplitude increases from black to white shading . as before
, the input pulse has a length of @xmath105 .
[ t][t]@xmath65 [ b][b]@xmath66 [ t][t]@xmath67 [ b][b]@xmath111    [ tr][tr]2 [ tr][tr]1 [ tr][tr]0 [ tr][tr]-1 [ tr][tr]-2 [ tr][tr]-3 [ tr][tr]-4    ( 0,0 ) ( -5,200)(a )     input and output wavefunction for a two - photon input pulse with a pulse length of @xmath105 .
( a ) is a contour plot of the input wavefunction @xmath15 .
( b ) is the cross - section of the contour line plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath112 ( broken line ) and @xmath113 ( dotted line ) . the horizontal axis @xmath72 in ( b ) is the average position of the two photons . likewise ,
( c ) is a contour plot of the output wavefunction @xmath39 and ( d ) shows cross sections of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath112 ( broken line ) and @xmath113 ( dotted line ) . ]
( 0,0 ) ( -20,200)(b )     input and output wavefunction for a two - photon input pulse with a pulse length of @xmath105 .
( a ) is a contour plot of the input wavefunction @xmath15 .
( b ) is the cross - section of the contour line plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath112 ( broken line ) and @xmath113 ( dotted line ) .
the horizontal axis @xmath72 in ( b ) is the average position of the two photons .
likewise , ( c ) is a contour plot of the output wavefunction @xmath39 and ( d ) shows cross sections of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath112 ( broken line ) and @xmath113 ( dotted line ) . ]
[ t][t]@xmath65 [ b][b]@xmath73    [ tr][tr]2 [ tr][tr]1 [ tr][tr]0 [ tr][tr]-1 [ tr][tr]-2 [ tr][tr]-3 [ tr][tr]-4    ( 0,0 ) ( -5,200)(c )     input and output wavefunction for a two - photon input pulse with a pulse length of @xmath105 .
( a ) is a contour plot of the input wavefunction @xmath15 .
( b ) is the cross - section of the contour line plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath112 ( broken line ) and @xmath113 ( dotted line ) . the horizontal axis @xmath72 in ( b ) is the average position of the two photons .
likewise , ( c ) is a contour plot of the output wavefunction @xmath39 and ( d ) shows cross sections of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath112 ( broken line ) and @xmath113 ( dotted line ) . ]       ( 0,0 ) ( -20,200)(d )     input and output wavefunction for a two - photon input pulse with a pulse length of @xmath105 .
( a ) is a contour plot of the input wavefunction @xmath15 .
( b ) is the cross - section of the contour line plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath112 ( broken line ) and @xmath113 ( dotted line ) . the horizontal axis @xmath72 in ( b ) is the average position of the two photons .
likewise , ( c ) is a contour plot of the output wavefunction @xmath39 and ( d ) shows cross sections of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath112 ( broken line ) and @xmath113 ( dotted line ) . ]    ( 0,0 ) ( -15,25)(e ) ( -5,-10)-2.0 ( 185,-10)1.5 ( 120,-10)0     input and output wavefunction for a two - photon input pulse with a pulse length of @xmath105 .
( a ) is a contour plot of the input wavefunction @xmath15 .
( b ) is the cross - section of the contour line plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath112 ( broken line ) and @xmath113 ( dotted line ) . the horizontal axis @xmath72 in ( b ) is the average position of the two photons .
likewise , ( c ) is a contour plot of the output wavefunction @xmath39 and ( d ) shows cross sections of the contour plot at @xmath62 , where the time differences are @xmath69 ( solid line ) , @xmath112 ( broken line ) and @xmath113 ( dotted line).,height=37 ]    figure .
[ fig : figure7 ] ( b ) shows the cross - section of the input wavefunction at @xmath62 for time differences of @xmath69 , @xmath112 , and @xmath113 between the two photons . figure .
[ fig : figure7 ] ( c ) and figure .
[ fig : figure7 ] ( d ) show the output wavefunction @xmath39 .
the output amplitude is now negative everywhere .
however , the pulse shape of the output shown in figure .
[ fig : figure7 ] ( c ) is quite different from the shape of the input pulse shown in figure .  [
fig : figure7 ] ( a ) . in order to illustrate the spatial features of the output wavefunction in detail , figure .
[ fig : figure7 ] ( d ) shows the cross - sections of the output wavefunction for time differences of @xmath69 , @xmath112 , and @xmath113 between the two photons .
note that the shape of the cross - sections of the output in figure .
[ fig : figure7 ] ( d ) is similar to that of the input cross - sections in figure .
[ fig : figure7 ] ( b ) .
however , the amplitudes are quite different . specifically , the output amplitude at a time difference of @xmath69 ( figure .
[ fig : figure7 ] ( d ) , solid line ) is larger than the input amplitude at the same time difference ( figure .
[ fig : figure7 ] ( b ) , solid line ) , while the output amplitude at a time difference of @xmath113 ( figure .
[ fig : figure7 ] ( d ) , dotted line ) is smaller than the input amplitude at the same time difference ( figure .
[ fig : figure7 ] ( b ) , dotted line ) .
this indicates that the main difference in pulse shape between the input and output pulses is that the output amplitude decreases faster than the input amplitude between @xmath69 and @xmath114 .
on the other hand , the delay time of the peak of the output amplitude does not change much .
its value is about @xmath115 for all time differences @xmath91 between the photons .
let us now consider why such features appear . figure .
[ fig : figure8 ] shows the contour plot of the three components given by eq.([eq : outputcomponents ] ) . figure .
[ fig : figure8 ] ( a ) corresponds to the transmission process , figure .  [
fig : figure8 ] ( b ) corresponds to the one photon absorption process and figure .
[ fig : figure8 ] ( c ) corresponds to the two photon absorption process .
as before , the component of the transmission process ( figure .
[ fig : figure8 ] ( a ) ) is identical with the input wavefunction , @xmath116 .
likewise , the component of the one photon absorption process ( figure .  [ fig : figure8 ] ( b ) ) can be described as @xmath117 since the single photon transmission @xmath58 is positive and the single photon absorption @xmath54 is negative , the total amplitude of this component is negative for all @xmath118 .
the component of the two photon absorption process ( figure .
[ fig : figure8 ] ( c ) ) includes the non - linear interaction between the two photons .
the positive amplitude of this component is therefore always smaller than the product of the single photon absorption amplitudes , @xmath119 .
specifically , the two photon absorption component @xmath89 is exactly zero at @xmath74 and increases rapidly up to a maximal amplitude at about @xmath120 . for larger photon distances , @xmath121 ,
@xmath89 decreases as a consequence of the decrease in the amplitudes of the photon absorption processes , @xmath119 .
note however that , at its maximal amplitudes near @xmath120 , the two photon absorption amplitude @xmath89 is still significantly lower than the product of the single photon absorption amplitudes , @xmath119 .
this means that the non - linear interaction of the photons significantly affects the whole output pulse in the case of an input pulse length of @xmath105 .
[ t][t]@xmath65 [ b][b]@xmath66 [ tr][tr]2 [ tr][tr]1 [ tr][tr]0 [ tr][tr]-1 [ tr][tr]-2 [ tr][tr]-3 [ tr][tr]-4    ( 0,0 ) ( -5,200)(a )     contour plots of the two photon processes @xmath122 .
( a ) corresponds to the transmission process , ( b ) corresponds to the one photon absorption process , ( c ) corresponds to the two photon absorption process . ]       ( 0,0 ) ( -5,200)(b )     contour plots of the two photon processes @xmath122 .
( a ) corresponds to the transmission process , ( b ) corresponds to the one photon absorption process , ( c ) corresponds to the two photon absorption process . ]
( 0,0 ) ( -5,200)(c )     contour plots of the two photon processes @xmath122 .
( a ) corresponds to the transmission process , ( b ) corresponds to the one photon absorption process , ( c ) corresponds to the two photon absorption process . ]     +    ( 0,0 ) ( -10,35)(d ) ( -10,-10)-0.2 ( 190,-10)1.5 ( 120,-10)0     contour plots of the two photon processes @xmath122 . ( a ) corresponds to the transmission process , ( b ) corresponds to the one photon absorption process , ( c ) corresponds to the two photon absorption process.,height=37 ]    in order to illustrate the dependence of the output component on the relative distance between two photons in detail , figure .
[ fig : figure9 ] shows the cross - sections at @xmath62 of the output amplitude @xmath39 ( solid line ) and the components of the transmission process @xmath123 ( dotted line ) , the one photon absorption process @xmath88 ( thin line ) , and the two photon absorption process @xmath89 ( broken line ) . the time differences are @xmath69 in ( a ) , @xmath112 in ( b ) , and @xmath113 in ( c ) . in all figures , the component of the one photon absorption process ( thin line ) is always negative and significantly greater than the component of the transmission process ( dotted line ) .
furthermore , the ratio of the peak of the one photon absorption amplitude ( thin line ) to the peak of the corresponding transmission amplitude does not change very much depending on the time difference @xmath91 between the two photons .
on the other hand , the ratio of the peak of the two photon absorption amplitude ( broken line ) to the peak of the corresponding transmission amplitude ( dotted line ) changes significantly depending on the time difference @xmath91 between the two photons as discussed in the following .
[ fig : figure9 ] ( a ) shows the cross - section of the output components at @xmath69 . in this case , the two photon absorption amplitude ( broken line ) is exactly zero .
the output amplitude is then given by the sum of the transmission amplitude ( dotted line ) and the one photon absorption amplitude ( thin line ) .
since the one photon absorption amplitude is about @xmath124 times as large as the transmission amplitude , the total output pulse ( solid line ) is negative and about @xmath125 times as large as the transmission amplitude .    figure .
[ fig : figure9 ] ( b ) shows the cross - sections of the output components at @xmath112 .
the one photon absorption amplitude ( thin line ) and the transmission amplitude ( dotted line ) are not much different from @xmath69 ( figure .
[ fig : figure9 ] ( a ) ) .
however , the two photon absorption amplitude ( broken line ) is now about half as high as the transmission amplitude ( dotted line ) .
therefore , the output amplitude is reduced to about the same height as the transmission amplitude .
[ fig : figure9 ] ( c ) shows the cross - sections of the output components at @xmath113 .
the two photon absorption amplitude ( broken line ) and the transmission amplitude ( dotted line ) are now significantly lower than for @xmath69 ( figure .
[ fig : figure9 ] ( a ) ) .
however , the negative one photon absorption amplitude is still about @xmath126 times as large as the transmission amplitude .
the two photon absorption amplitude ( broken line ) has increased to about @xmath125 times the transmission amplitude , leaving a total negative output amplitude ( solid line ) of about @xmath127 times the transmission amplitude .
thus the output amplitude is reduced to very low values , even though the components themselves are still very large .
the fast decrease of the total output amplitude shown in figure .
[ fig : figure7 ] ( c ) and [ fig : figure7 ] ( d ) can then be understood as a consequence of the rapid increase of the two photon absorption amplitude with @xmath91 .
interestingly , the delay times of the pulse components are not sufficiently high to cause the same change in the sign of the amplitude observed for the one photon case ( figure .
[ fig : figure6 ] ( a ) ) .
no part of the two - photon output wavefunction is dominated by either the transmission or the two photon absorption amplitude .
this is mostly because the two photon absorption amplitude is suppressed by the non - linearity .
in fact , it can be seen in figure .
[ fig : figure9 ] ( c ) that the nonlinearity even reduces the delay of the two photon absorption amplitude from a delay time of @xmath110 for @xmath119 to only about @xmath128 for @xmath89 at @xmath113 . therefore , the change of the pulse shape in the two photon pulses is actually less drastic than the change of pulse shape observed in the one photon case .
[ t][t]@xmath67 [ ] [ ] [ ] [ ] [ ] [ ]    ( 0,0 ) ( -5,210)(a )     cross - sections of the output amplitude @xmath39 and its components @xmath86 at time differences @xmath85 of ( a ) @xmath69 , ( b ) @xmath112 and ( c ) @xmath113 .
the solid lines correspond to the total output amplitude @xmath39 , the dotted lines correspond to the transmission amplitude @xmath87 , the thin lines correspond to the one photon absorption amplitude @xmath88 , and the broken lines correspond to the two photon absorption amplitude @xmath89 . ]
( 0,0 ) ( -5,210)(b )     cross - sections of the output amplitude @xmath39 and its components @xmath86 at time differences @xmath85 of ( a ) @xmath69 , ( b ) @xmath112 and ( c ) @xmath113 .
the solid lines correspond to the total output amplitude @xmath39 , the dotted lines correspond to the transmission amplitude @xmath87 , the thin lines correspond to the one photon absorption amplitude @xmath88 , and the broken lines correspond to the two photon absorption amplitude @xmath89 . ]     +    ( 0,0 ) ( -5,210)(c )     cross - sections of the output amplitude @xmath39 and its components @xmath86 at time differences @xmath85 of ( a ) @xmath69 , ( b ) @xmath112 and ( c ) @xmath113 .
the solid lines correspond to the total output amplitude @xmath39 , the dotted lines correspond to the transmission amplitude @xmath87 , the thin lines correspond to the one photon absorption amplitude @xmath88 , and the broken lines correspond to the two photon absorption amplitude @xmath89 . ]
we have studied the effect of the strong nonlinearity of a one - dimensional atom on one- and two - photon input wavepakets .
the response of the atom can be characterized by transmission and absorption processes . for long pulses ,
the absorption components dominate in the output .
however , in the two photon case , the two photon absorption component is suppressed around @xmath74 by the saturation of the two - level atom .
the long two photon pulse therefore has a non - linear region of negative amplitude dominated by the one photon absorption for @xmath129 , and a linear region of positive amplitude dominated by the two photon absorption for @xmath130 . for a short pulse length of @xmath105 , the two photon absorption is strongly suppressed throughout the whole output pulse .
the one photon absorption process then dominates everywhere , and the amplitude of the short two - photon output pulse is completely negative .
it is therefore possible to obtain a very strong non - linear effect by adjusting the pulse length @xmath44 to the dipole relaxation time of the atom .    besides the total amplitude
, we have also characterized the change in delay time and pulse shape . for long pulses ,
the delay of the non - linear region is reduced to @xmath103 of the delay of the liner region . in the long pulse case
, the interaction of the two photons thus significantly reduces the average delay time caused by the absorption and reemission processes at the atom and causes a faster reemission into the output field . in the case of a short one photon pulse ,
the delay between the transmission and the absorption component is so great that the pulse shape is changed significantly .
in particular , the single photon output amplitude has both negative and positive regions . on the other hand , the output pulse shape for a short two photon pulse
is changed much less due to the reduced delay caused by the non - linearity .
it is therefore possible to obtain a negative output amplitude throughout the short output pulse .
this result indicates that the atomic nonlinearity is particularly strong at a pulse length of @xmath105 .
the analysis of the output in terms of components corresponding to different interaction processes thus allow us to understand the origin of the non - linear photon - photon interaction in great detail .
these insights may be very useful for the development of novel non - linear optical devices operating at the few photon level .
this work was partly supported by the program `` research and development on quantum communication technology '' of the ministry of public management , home affairs , posts and telecommunications of japan .
99 nielsen ,  m. , chuang ,  i. ( 2000 ) .
_ quantum computation and quantum information .
_ cambrige : cambridege university press .
hofmann ,  h. , kojima ,  k. , takeuchi ,  s. , sasaki ,  k. ( 2003 ) .
entanglement and four - wave mixing effects in the dissipation - free nonlinear interaction of two photons at a single atom .
a , 68 _ , 043813 .
milburn ,  g. , walls ,  d. ( 1983 ) .
quantum nondemolition measurements via quadratic coupling .
rev .  a , 28 _ , pp .  2065 - 2070 ; n.  imoto , h.  a.  haus , y.  yamamoto . ( 1985 ) .
quantum nondemolition measurement of the photon number via the optical kerr effect .
rev .  a , 32 _ , pp .
2287 - 2292 .
mills ,  d. ( 1991 ) . _ nonlinear optics ( basic concepts ) . _
new york : springer - verlag .
bermann ,  p. ( 1994 ) .
_ cavity quantum electrodynamics supplement 2 to advances in atomic , molecular , and optical physics .
_ san diego : academic .
turchette ,  q. , hood ,  c. , lange ,  w. , mabuchi ,  h. , kimble ,  h. ( 1995 ) .
measurement of conditional phase shifts for quantum logic .
lett . , 75 _ , pp .  4710 - 4713 .
hofmann ,  h. , kojima ,  k. , takeuchi ,  s. , sasaki ,  k. ( 2003 ) .
optimized phase switching using a single - atom nonlinearity _ j.  opt .
b : quantum .
semiclass .  opt .
, 5 _ , pp .  218 - 221 .
kojima ,  k. , hofmann ,  h. , takeuchi ,  s. , sasaki ,  k. ( 2003 ) .
nonlinear interaction of two photons with a one - dimensional atom : spatiotemporal quantum coherence in the emitted field . _ phys .  rev .
a , 68 _ , 013803 .
hofmann ,  h. , mahler ,  g. ( 1995 ) .
measurement models for time - resolved spectroscopy : a comment .
_ quantum .
semiclass .
opt . 7 _ , pp .  489 - 497 .
rice ,  p. , carmichael ,  h. ( 1988 ) .
single - atom cavity - enhanced absorption i : photon statistics in the bad - cavity limit .
_ ieee j.  quantum  electron . , 24 _ , pp
.  1351 - 1365 .
truchette ,  q. , thompson ,  r. , kimble ,  h. ( 1992 ) .
one - dimensional atoms .
b , 60 _ , pp .  s1-s10 .