In this work, a chain of N identical two-level atoms coupled with a quantized electromagnetic field, initially prepared via a single-photon Fock state, is investigated. The N-particle state amplitude of the system is calculated for several space configurations of the atoms in the Weisskopf-Wigner approximation. It was shown that the space configuration of an atomic chain, the total number of atoms, and even the available volume for the field modes define the behavior of the system state amplitude with time. Applying the condition of ‘cyclic bonds’, presented in this work, to the elaborated theory allows to describe the system time evolution, practically, for any space configuration.
Keywords:Atoms; Fock state; Weisskopf-Wigner approximation
The collective absorption (emission) of photons by an ensemble of identical atoms ‘provides valuable insights into the many-body physics of photons and atoms’ (quoted from ). Taking into account the quantization of electromagnetic field, many fundamental and interesting properties of the coupled systems of atoms and field are revealed. For example, when the average distances between atoms are much less than the ‘resonant transition’ wavelength of emitted (absorbed) light, the cooperative coupling leads to a substantial radiative shift of the transition energy and significant change in decay rate of the ensemble state. The latter was revealed through the various theoretical (for example, some relatively modern researches in [2-5]) and experimental investigations (see starting, for example, from [6,7] to the modern applications like described in  and impressively effective experimental realizations as in ). Some peculiar behavior in spontaneous emission is proper even in a system of atoms which can have a relative distance larger than the emission wavelength (see, for instance, ), and initially, only one atom or one-photon state is excited, as discovered in this paper.
In the present paper, a system (chain) of N identical two-level non-interacting atoms, prepared ‘via a single-photon Fock state’ in the one- or two-mode resonator, is investigated. The main goal of the paper is to obtain the information about the state of electromagnetic field and atomic system (chain) in a Weisskopf-Wigner approximation (see  chapter 6, page 206 and some comments in ). The calculations of the state amplitudes of the atomic system are made for several approximations in resonator (cavity) characteristics and for several types of space configurations of the atoms in the chain.
In this work, we study the case, in which the distances between atoms are quite large, so that the average distances between atoms are greater or in the same order than the ‘resonant transition’ wavelength. Therefore, we prepare an ensemble of N two-level atoms initially in ground state, and a single mode of the radiation field is excited in a ‘Fock’ state (so called one-photon state). This is the case of a purely monochromatic wave with zero line width under the consideration. A laser output in single mode operation can approximate this situation due to its high degree of monochromaticity (small line width) for instance. The mode of electromagnetic field is specified completely by giving its wave vectors k0 with atomic transition frequency ω = c|k0| and its polarization j (j = 1, 2).
The main feature, differentiating our research from others in this domain, is the developed direct and consistent solution to the N-particle equations, describing the time evolution of the N atomic probability state amplitudes. Besides, in certain sense, we explained the nature of the widely used Weisskopf-Wigner approximation that was not found in the reviewed by us scientific literature.
The goal of this paper can be formulated as an attempt to propose an adapted and simple in practical use theory, for example in the highly applied nanoscale physics. The proposed theoretical material requires corresponding experimental verification. As an idea of an application, the model system can be realized on atomic (developing the method proposed in  for the nuclei of 57Fe in certain composites, but this time for a visible region), chains of trapped ions (like in ), and molecular structures for further developing such techniques like FRET (described for instance in ), atomic chains like carbyne loops (for example, ), and microhole array synthesized by femtosecond laser radiation (see , for an instance).
Let us first provide below some general theoretical premises. More detailed derivations of the corresponding mathematical model can be found in .
The equations of motion for the state amplitudes
We have assumed that the atomic energy levels have no linewidth, so that, only if , the atoms can be able to absorb a photon. Obviously, this is an unrealistic case since it is impossible to have a completely monochromatic wave. In addition, for the case of the Fock initial state, in which we measured the energy precisely of the mode, the average electric field will be zero. In the forth of the law of energy conservation, an emitted photon will correspond to the same frequency (we can say it will occur with a high probability after a quite long time interval if the system has a damping).
Therefore, consider a collection of N identical atoms, at positions r1,…,rα,…,rN, coupled to a one mode electromagnetic (EM) field. Each atom α = 1..N is assumed to have only the two states |a〉α and |b〉α, separated by energy . In the dipole approximation, the closed conservative system of identical atoms with the electromagnetic field in a cavity can be described by the Hamiltonian consisting of free atoms and electromagnetic field items with dipole-field coupling between the atoms and the electromagnetic field modes.
Inasmuch as at the initial time moment t = 0 all atoms α = 1..N of the ensemble are in the ground state |b〉α and EM field is in Fock state (that presents one photon with the wave vector k0), we look for a solution of the corresponding Schrödinger equation in the interacting picture in the following form:
with the initial conditions:
where is Kroneker’s delta symbol. if k = k0, and if k ≠ k0. βα (t) (α = 1..N) and γk,j (t) (j = 1, 2) are the αth atom excited state amplitude with the others in the ground states and excited Fock field state amplitude of the jth polarization with the wave vector k, accordingly.
Then, the corresponding Schrödinger equation in the interacting picture yields the following system of equations:
Here, θk,j is the angle between dipole transition vector ℘α (more accurately, non-diagonal dipole matrix element) and the jth unit polarization vector ek,j (j = 1,2 and ek,j · k= 0). V is the available by the system of atoms and field space volume. The frequencies νk correspond to the modes with the module of the wave vectors k equal to |k|.
Therefore, substituting Equation 4 into 3 and differentiating one more time, after applying the Weisskopf-Wigner approximation (details in ), we can derive the following system of evolution equations:
And the decay rates Dα in the approximation can be estimated by the formula:
The coefficient Dα (α = 1..N) describes the respective rate of decay for αth atom excited state. Note, that the ‘non-resonant’ items for the particle with distinguished from α indexes were disregarded in here in an assumption of quite large interatomic distances (see details in ).
Results and discussion
An atomic chain with cyclically distanced atoms
Next, we try to make the calculations, using here the particular case of space configuration for the system atoms field. Below, for simplicity, only one polarized mode (j = 1) of the resonant field modes is taken into account with the common parameters gα and ℘α for α = 1..N :
for |k| = k0. In other words, the space angle distribution for the components Φαδis disregarded here, assuming the direction of the transition dipole moment ℘α for any atom in the system coincides with the photon polarization in absorbing or emitting a resonant photon.
Then, from the system of Equation 7, in the case of a cavity with two resonant modes k = ± k0 and identical atoms with Dα ≡ D for α = 1..N, one derives that
Using the notation
Therefore, taking into account the initial conditions βδ (0) = 0 for α = 1..N, the solution of the above equation is as follows:
The field probability amplitudes can be obtained using the subsystem of Equation 4 of the full ‘conservative’ system of Equations 3 and 4. Therefore, substituting (15) and (17) into the Equation 4, and then taking into account the restrictions βα (0) = 0 for α = 1..N, we obtain that
Note, here, we neglected the possible space angle distribution for the direction of the resonant wave vector k.
Inasmuch as cos(k (rα - rδ)) = cos (krα) cos (krδ) + sin (krα) sin (krδ), then, after substitution of the found superpositions (15) and (17) into the initial Equation 12, we derive the following integrable differential equation:
Integrating the left and right sides of the equation above (22) over time yields
According to the definition of the functions Fc,s (t)
The solution of such linear first order differential equation, like (23), has the form:
The integration in the last expression can be performed, yielding
The initial condition βα (0) = 0, for α = 1..N, sets the coefficient C0 equals 0. The initial time derivative can be determined, for example, if the system of Equation 3 from the initial ‘conservative’ full system of Equations 3 and 4 is chosen as a basis at the time moment t = 0. Then, the initial condition for the field state amplitude γk (0) = 1, where k = k0, sets the time derivative to the following expression:
Now, the question arises how to choose correctly the coefficients C and C′. First of all, the choice has to satisfy the limitations on the probability amplitude, yielding the corresponding probability limited above by unit (the sum of all the modules squared of the introduced amplitudes equals unit probability). Secondly, the solution with the coefficients have to be consistent with the model decay (damping).
We observe that, formally, when the real part of the variable Ω is a negative quantity, that is Re (Ω) < 0, the introduced functions H and f have the following limits for quite long time intervals:
As for an open system, in our case, it should be expected for a quite long time interval the total electromagnetic energy of the atoms-field system to be emitted into the subsystem causing the state damping. Therefore, let us define the coefficients C and C′ in the following manner:
Then, after substitution into the expressions for the time limits, one derive the logical finale of the system evolution:
The possible space configurations of the atomic system, satisfying the condition of ‘circularity’, can be easily found. For example, the set s3a1 (the notation ‘s3a1’ is just introduced here): , , and kr3 = π. As an instance, it can also be the set s3a2: , , and . A space configuration for five atoms can be represented, for example, by such set like s5a1: , , , , and . Specifically, a combination, such like with i = 1..N, n > 2, and , can generate the necessary magnitudes of the characteristic system frequencies Ω2 and (that, actually, are the corresponding Rabi frequencies), comparable with the given magnitude of the decay coefficient D.
Below we depict the atomic system behavior in the several introduced above configurations. Note, that the cited thereby Rabi frequencies were calculated in the SI system of units with the following notations: ; the electric permittivity of free space ε0 ≈ 8.8542 × 10-12 F/m; the speed of light in free space c = 299792458 m/sec; resonant wavelength close to the D2-line of a sodium atom λD ≈ 589.29 × 10-9 m; corresponding circular (in radians per second) resonant frequency ; non-diagonal so called ‘transition’ dipole matrix element (in the same order as for the D2-line transition, that is about 1 Debye) ρex = 1 × 3.33564 × 10-30 C m. For instance, if the available for the system of atoms and field volume has the value equal to V = 0.001 m3, then .
Assume, for example, the available volume V = 10-13 m3 is somehow filled by the set s3a1 with D ≈ 107 rad/sec, initially coupled with one-photon Fock state. Then, , , and . The corresponding graphs for probability to find each atom in the excited state are shown in Figure 1.
Figure 1. Time evolution of |βα (t)|2. V = 10-13 m3. Atoms are arranged in the set s3a1 with D ≈ 107 rad/sec. The bold solid line represents the atom with the space phase kr1 = π/6, the dot line is for the space phase kr2 = 2π/3, and the thin solid line corresponds to kr3 = π.
Let us see what happens when the available volume is increased by one order. This yields V = 10-12 m3 with the same three atoms (D ≈ 107 rad/sec) of the configuration s3a1. Then, ; and . The corresponding graphs for each atom excited state probability are depicted in Figure 2.
Figure 2. Atom excited state probability |βα (t)|2. V = 10-12 m3. Atoms are arranged in the set s3a1 with D ≈ 107 rad/sec. The bold solid line represents the atom with the space phase kr1 = π/6, the dot line is for the space phase kr2 = 2π/3, and the thin solid line corresponds to kr3 = π.
Suppose now that the available volume is V = 10-13 m3, somehow filled by the set s5a1 with D ≈ 107 rad/sec initially coupled with one-photon Fock state. Then, ; , and . The corresponding graphs for each atom excited state probability are shown in Figure 3.
Figure 3. Atomic excitation probability |βα (t)|2 as a function of time. V = 10-13 m3. Atoms are arranged in the set s5a1 with D ≈ 107 rad/sec. The bold solid line represents the atom with the space phase kr1 = 2π/3, the dot line is for the space phase kr5 = 19π/6, and the thin solid line corresponds to kr3 = 5π/2.
And again, let us see what happens when the available volume is increased by one order. This yields V = 10-12 m3 with the same five atoms (D ≈ 107 rad/sec) of the configuration s5a1. Then, ; , and . The corresponding graph is in Figure 4. Note that the graphs Figures 3 and 4 of excited state probabilities are for the chosen three atoms with the following phases: , , and .
Figure 4. Probability |βα (t)|2. V = 10-12 m3. Atoms are arranged in the set s5a1 with D ≈ 107 rad/sec. The bold solid line represents the atom with the space phase kr1= 2π/3, the dot line is for the space phase kr5 = 19π/6, and the thin solid line corresponds to kr3 = 5π/2.
As it was supposed in the derivative of the differential equations with the damping items such like (12) (see the details in the work , the available volume V for the system of atoms and field defines the ‘available’ modes for the electromagnetic field. The value of volume V can determine one of the inequalities D < Ω2 and D > Ω2 ( and ), therefore defining the character of the system relaxation. Such fundamental system property was illustrated in the figures. It is interesting to note that increasing the system volume V, therefore increases the ‘available’ number of quantized field modes, the maximum probability to find an atom in its excited state decreases. Other interesting feature, shown in the proposed graphs, is the different character of relaxation for each excited atom. The latter depends, as shown here, on the space phase krα, where α = 1..N.
On this note, therefore, let our narration to come to the following conclusions, in short.
Thus, in this work, we investigated a chain of N identical two-level long distanced atoms prepared ‘via a single-photon Fock state’. The functional dependence of the atomic state amplitudes on a space configuration and time is derived in the Weiskopf-Wiegner approximation.
It was shown that in increasing the system volume V, the maximum value of probability to find an atom in its excited state decreases. The feature can be experimentally investigated at the proposed nanoscale limit for the space configuration of atoms.
Hence, the Weiskopf-Wiegner approximation was revealed through the provided application to the many-body system at the nanoscale limit for the atomic space phases. The found solution (30) cannot be counted as a particular one, or as a limit of such, for the initial systems of Equations 3 and 4 that represent only a closed conservative system of atoms and an electromagnetic field. Thus, we can say that the model described in this work, besides the atoms and the electromagnetic field, implicitly contains a third participant guaranteeing a total system relaxation with time. It is interesting to note here that the ‘complete’ decay of the system excitations was strongly imposed by the choice of the coefficients C (38) and C′ (39).
The methods, described in this work, of solving the system of linear differential equations can be applied even for more general situations when the boundary ‘circular’ conditions are not satisfied. In certain cases, the problem allows to extend the system by adding a subsystem of non-sufficient number of atoms (or an atom when the expression is actual) without influencing the main characteristics under our interest. Besides, it is interesting to investigate the dependence of value of the mentioned construction on the location of the origin of a coordinate system.
The authors declare that they have no competing interests.
SAS performed the calculations and analysis of the results. YSM analyzed the methods of investigation and drafted the manuscript. Both authors read and approved the final manuscript.
Texas AM University, Physics Department; 2008. http://repository.tamu.edu/bitstream/handle/1969.1/ETD-TAMU-3226/CHANG-DISSERTATION.pdf?sequence=1 webcite
Phys Rev Lett 1973, 30:309. Publisher Full Text
Phys Rep 1982, 93:301. Publisher Full Text
Nanoscale Res Lett 2013,.