The ultracold Rydberg atoms have been extensively studied due to their strong and controllable long-rang interactions. The interaction-induced dipole blockade effect has been proposed as an ideal candidate for realizing a scalable quantum logic gate applying in quantum information processing   . The dipole blockade effects based on ultracold Rydberg atoms have been experimentally studied for van der Waals interaction case by changing density or the principal quantum number of Rydberg state   , and for the dipole-dipole interaction case by tuning an electric field  . The ultracold Rydberg atoms move only a few percent of typical interatomic spacing on the time scale of experimental interest; in other words, the ultracold Rydberg atoms have smaller thermal kinetic energy relative to their average interaction energies, so ultracold Rydberg atoms provide an alternative platform to study interaction induced collision dynamics and evolution process. Some groups experimentally investigated collision and ionization mechanism of a pair of Rydberg atoms initially excited on the attractive potential for dipole-dipole interaction case  , and van der Waals interaction case  , which will give rise to the faster collisions and even avalanche ionization forming ultracold plasmas  .
Blackbody radiation (BBR) is known to strongly affect the populations of Rydberg atoms by transitions from initial state to nearby dipole allowed states or BBR-induced direct ionization  . Amthor et al. showed that the Rydberg atoms’ collisional ionization was also observed on longer time scales when atoms are excited to the repulsive van der Waals potential, which is attributed to the BBR redistribution to the effective attractive interaction potential  . And the initial thermal motion of the ultracold Rydberg atoms can significantly shorten the collision time  . For a pair of Cs nP atoms, penning ionization and the formation of an ultracold plasma mostly occur on attractive potential for states with n > 42  . So the long-range potential character of Rydberg atoms determines the collisional and ionization dynamic mechanism. However, above works on collision and ionization induced by Rydberg atomic interactions were mainly focused on the spectroscopic measurements based on the pulsed field ionization techniques, which have proven to be a powerful measuring method for Rydberg atoms and utilized widely. Based on the analysis of the decay curve provides a new method for studying Rydberg atoms interactions and collisional ionization process, which is significantly affected by the Rydberg atomic strong interaction and a weak background DC electric field.
In this paper, we in detail analyze the collisional loss mechanisms induced by Rydberg atoms interactions for different Cs nP states using the theoretical model used before. The validity of this theoretical model has been proved in our previous article by comparing the theoretical and experimental values of collisional loss rate coefficient of Cs 63S state  , but it is different from this paper: First, for Cs nS Rydberg states, the interaction potentials are repulsive, and we only study the effect of repulsive interactions on collisional loss, and consider simultaneously the BBR effect. For the Cs nP states, the interaction potential only for n < 42 is repulsive, but for n > 42 is attractive potential. In this paper, we can simultaneously study the collisional loss of nP states for attractive and repulsive case, and compare the effect on collisional loss rate under the same condition, especially the role of BBR in this process for two cases; in addition, the different spontaneous radiation rates, BBR transition rate and their relative ratio for nS and nP Rydberg states will also lead to different collisional loss mechanisms. Second, we investigate the important role of initial van der Waals interaction in total collision process for attractive and repulsive case, in detail analyze the effect on the total collisional loss rate, and improve the analytical formula for both cases. This content has not been studied in previous papers.
2. Theoretical and Analysis
The collisional characteristics of cold atoms are conventionally determined from the analysis of the decay curve. We can obtain the analytic solution of the time evolution of the Rydberg atoms number 
Here, α is the loss rate coefficient including the spontaneous radiation rate GSpon and blackbody radiation rate GBBR. The blackbody radiation rate include not only the blackbody radiation transition rates to all nearby Rydberg states through absorption and stimulated emission, but also the direct blackbody radiation photoionization rate. β(βRy-G) is the loss rate coefficient caused by collisions between ultracold Rydberg-Rydberg (Rydberg-ground) atoms; βRy-e is the loss rate coefficient caused by collisions between ultracold Rydberg atoms and electrons. N is the atomic number; subscript Ry (G) denotes Rydberg (ground) state; wRy, wG and we are the waist radius of Rydberg atoms, cold atoms sample and ionized electrons, respectively.
The value of dipole elements between Rydberg-ground atoms is two orders of magnitude smaller than that of between Rydberg atoms, and decrease with increasing of n, so the collisional loss rate induced by the weak interactions between Rydberg-ground atoms has negligible effect on the total collisional loss rate. Although the large geometrical cross section of Rydberg atoms will lead to the large quenching collisional loss for higher Rydberg states, but has weak contribution to the total collisional loss rate for n < 50  , so the Equation (2) can be rewritten as
In order to investigate the collisional loss induced by Rydberg atomic interactions, we detailed analyze the relation between collisional loss rate coefficient and interactions. The collisional loss rate coefficient can be written as , the sRy and vRy are the collisional cross section and relative collisional velocity of ultracold Rydberg atoms, both depend on their van der Waals interaction or dipole-dipole interaction. The collisional cross section induced by interactions given by 
where mCs is Cs atom mass; vc ≈ (8kBT/pmCs)1/2 (kB and T are the Boltzmann constant and the temperature of cold atoms) is the average velocity of cold atoms; Cn is the interaction coefficient, which is given by theoretical calculations  ; f(3) = 4 and f(6) = 15π/8 are f(n) values for van der Waals interaction and dipole-dipole interaction case. For a pair of Rydberg atoms separated by a distance R, they will experience a force and be accelerated toward (away from) each other if they are initially excited on attractive (repulsive) potentials. The force is the gradient of the interaction potential, the relative accelerations caused by van der Waals interaction force and dipole-dipole interaction force can be written as avdW ≈ 6C6/(mCsR7) and aDD ≈ 3C3/(mCsR4), respectively. vesc is a minimum escape velocity of Rydberg atomic from center to edge of average interaction region under corresponding acceleration, i.e., vesc = (2ar)1/2, the average interaction region is approximately equal to the Rydberg excitation volume with the average radius r, in actual experiment, the overlap between trapping laser beam and the Rydberg excitation laser beam defines an effective Rydberg excitation volume. We also obtain the relative collisional velocity between Rydberg atoms under corresponding acceleration, vRy = (2aR)1/2, the detailed analysis process has been described in our previous work  . Combining the above relations, the corresponding analytic solution of the collisional loss rate coefficient βvdW and βDD induced by van der Waals interaction and dipole-dipole interaction can be written as
Due to the existence of BBR, some Rydberg atoms in initial nP state will be gradually transferred to nearby dipole-allowed n’D and nS states, so the initial weak van der Waals interaction between Rydberg atoms also gradually evolve into strong dipole-dipole interaction between initial nP state and nearby dipole-allowed states, these strong attractive dipole-dipole interactions will give rise to faster collision and even ionization between Rydberg atoms. In order to obtain the effective collisional loss rate coefficient including the contributions of all larger dipole-allowed states, we use an effective coefficient instead of C3 in Equation (6), where μi is dipole matrix elements between initial nP state and nearby dipole-allowed states, pi is the transition probability from the initial nP Rydberg atoms to these dipole-allowed states, which is proportional to the strength of transition dipole moments. The effective collisional loss rate coefficient βDD and βvdW of nP3/2 at 3 μm atoms separation as a function of n are shown in Figure 1. The corresponding parameters for these and following calculations including r = 220 μm and T = 100 μK. The collisional loss rate coefficients induced by the interactions between nP3/2 and some nearby larger dipole-allowed states for nP3/2 − nS1/2, nP3/2 − (n + 1)S1/2 and nP3/2 − (n − 1)D5/2 are also shown in Figure 1. It clear seen that the effective collisional loss rate coefficient induced by dipole-dipole interactions are much larger than that of the van der Waals interactions at the same state, which indicate that the strong dipole-dipole interactions between initial nP and dipole-allowed states
Figure 1. The variation of βvdW, the effective βDD of nP3/2 and with some nearby larger dipole-allowed states as a function of n at R = 3 μm.
are the dominant collisional loss mechanism. Because the values of dipole matrix element between initial nP and nearby dipole-allowed states increase with increasing of n, and BBR transitions rates gradually equal to even exceed the corresponding spontaneous decay rates with n increasing, which lead to the effective collisional loss rate coefficient slightly nonlinear increase.
In Figure 1, we just separately analyzed the variation of the collisional loss rate coefficient with n for van der Waals interactions and dipole-dipole interactions case, we did not investigate the influence of atomic separation variation induced by the initial Rydberg atomic interactions on collisional loss rate coefficient, and we also require obtain the total collisional loss rates coefficient. For the Cs nP state, the interaction potentials are repulsive potential for n < 42, while n > 42 for the potential are attractive  . For n < 42 sates, many experiments observe the appearance of ionized ions on typical 10 μs delay time t   if ultracold Rydberg atoms initially excited on the repulsive van der Waals potential, in other words, this delay time is the sum of Rydberg atomic transition time form initial repulsive nP state to dipole-allowed states, and their collision time under attractive dipole-dipole interactions between nP state and dipole-allowed states. The interatomic separation will increase under their repulsive van der Waals forces, which maybe affect the latter interaction strength and collisional distance between nP state and dipole-allowed states. For n > 42 states, the interatomic separation will decrease under their initial attractive van der Waals forces, and even sharp decrease because the initial weak van der Waals interactions forces gradually evolve into the strong dipole-dipole interactions forces, the average collision time under their initial attractive van der Waals forces equal to, even shorter than the time of the BBR redistribution, which maybe lead to the contributions of interactions between nP state and dipole-allowed states gradually diminished, so we must investigate the effect of atoms separation variation on collisional loss rate under initial van der Waals interactions.
In order to obtain the moving distance of atoms compare to the initial atoms separation, we need to consider the relative acceleration of Rydberg atoms pair under initial interaction force and initial velocity of ultracold Rydberg atoms, which is given by ΔR = vct + 1/2at2. Figure 2 shows the variation of ΔR with time for different Rydberg states with 3 μm initial atomic separation, the initial Rydberg atomic velocity approximately equal to average velocity of cold atoms vc, which is about 18 cm/s at cold atoms temperature of 100 μK. It clear seen that the ΔR rapidly increase with delay time and even much greater than their initial atoms separation for higher n. For lower Rydberg states, the weak interactions give rise to relatively small distance even if long delay time. But the initial velocity of ultracold Rydberg atoms plays an important role, which will give rise to the relative moving distance up to about 1 μm at typical 5 μs delay time, the similar behaviors are experimental investigated for Rb Rydberg atoms  . When we only consider the Rydberg atomic initial velocity effect at 5 μs delay time, the variation of the βvdW with n are shown in Figure 3, the value of C6 is the average value of all molecular states. The Rydberg atomic separation after 5 μs delay time will gradually increase to about 4 μm for n < 42 and gradually decrease to about 2 μm for n > 42 because their different interaction characters. In contrast to the results of the fixed initial atoms separation R = 3 μm, we can see that the collisional loss coefficient βvdW increases with the decreasing of R for n > 42, and have slight decrease for n < 42 after the same change of R, which are shown in Figure 3 with green triangles and red circles, respectively. For n > 42, the interatomic separation will rapidly decrease due to their attractive van der Waals interactions, which also lead to larger collisional loss rate coefficient βDD compare to the initial atoms separation if the BBR has been occurred. Analogous to the
Figure 2. The variation of ΔR with delay time for different Rydberg states.
Figure 3. βvdW as a function of n for different effects contributions, compared to fixed initial atomic separation of 3 μm for all states. The green triangles and red circles only include the contribution of initial Rydberg atomic velocity; the blue triangles consider the contributions of initial Rydberg atomic velocity and van der Waals interactions.
above results and analysis, the variation of the βDD with n due to the change of interatomic separation induced by initial velocity are also shown in Figure 4, which have similar results with the van der Waals interactions case.
When we simultaneously consider the initial Rydberg atomic velocity and the van der Waals interactions at the same delay time, because the relative acceleration between ultracold Rydberg atoms will rapidly increase with n scaling as n11, the relative moving distances of ultracold Rydberg atoms for high states are much more than that of the low states at the same delay time, which will lead to more obviously changes of interactions with n. The variation of the βvdW considering the above two effects for n > 42 are shown in Figure 3 with blue triangles. We can see that the collisional loss rate coefficient obviously increase due to faster decrease of the interatomic separation under their initial attractive van der Waals interactions, compare to the case of only considering the initial Rydberg atomic velocity contribution. The similar variation of the βDD for n > 42 under the effect of the above two factors are also shown in Figure 4 with red circles. Due to slightly increase of atomic separation induced by the weak repulsive van der Waals interactions n < 42 at this short and even longer delay time, which will give rise to the negligible variation of the collisional loss rate coefficient, so we does not give these corresponding results in Figure 3 and Figure 4.
The total collisional loss rates induced by ultracold Rydberg-Rydberg interactions not only include the attractive (repulsive) interactions between initial nP states, but also the attractive dipole-dipole interactions between nP and nearby dipole-allowed states. Whether Rydberg atomic pairs initially excited on repulsive or attractive potential, the BBR redistribution effect always exists and must
Figure 4. βDD as a function of n for different effects contributions, compared to fixed initial atomic separation of 3 μm for all states. The triangles only include the contribution of initial Rydberg atomic velocity; the red circles consider the contributions of initial Rydberg atomic velocity and van der Waals interactions.
be taken into account, especially for the repulsive potential case. In order to comprehensive consider all Rydberg interactions contributions to collisional loss rate, we define a ratio such that F = GBBR/α, which is the transition ratio from initial nP state to nearby dipole-allowed states induced by BBR, compare to the total decay rates, the corresponding values of GBBR and α are obtained by theoretical calculations from reference  . For n > 42 attractive potential case, Rydberg atoms will accelerate toward each other, collide and even ionize under the initial attractive van der Waals interactions, the average collision time at R = 3 μm is about 8 μs for n = 45 by  , which is much shorter than the corresponding spontaneous emission lifetime and the effective lifetime at room temperature of 300 K  . This relative long collisional time will lead to most atoms been involved in the collisional process under initial attractive van der Waals interactions before BBR occurrence, and then gradually experiences the attractive dipole-dipole interactions between nP and nearby dipole-allowed states. So the total collisional loss rate coefficient of ensemble atoms system for n > 42 can be written as
where ηM is enhancement factor of many-body effects because Rydberg atoms in many atoms system maybe collide each other several times  , we consider many-body effects with ηM = 10 in the following calculations; The van der Waals interactions will introduce gradually increasing contribution for the total collisional loss rate coefficient, so the attractive van der Waals interactions between initial nP states for n > 42 have non-negligible effect on the total collisional loss compared to the strong attractive dipole-dipole interactions between initial nP and nearby dipole-allowed states.
For n < 42 repulsive potential case, Rydberg atoms will move away from each other under the initial repulsive van der Waals interactions, the collisional ionization does not occur. The collisional probability of ultracold Rydberg atoms initially excited on attractive potential is two orders of magnitude higher than that of repulsive potential case even consider the BBR redistribution effect  , so we only need to consider the attractive dipole-dipole interactions between initial nP state and nearby dipole-allowed states for n < 42, the corresponding total collisional loss rate coefficient can be written as
At 3 μm initial atoms separation and typical 5 μs delay time, the total collisional loss rate coefficient βT of nP3/2 versus n is shown in Figure 5 with blank blocks. It is clear that the total collisional loss rate obviously enhance with n because of collision effect induced by increasing attractive interactions for n > 42 compare to n < 42 case. The total collisional loss rates for n > 42 include not only the contribution of attractive interactions between initial nP state, but also the contribution of attractive interactions between nP and nearby dipole-allowed states. While the increasing slope of curves for n > 42 states mainly attributed to the important role of the initial van der Waals interaction, the results for n > 42 states without considering the atomic separation effect induced by the initial van der Waals interaction are also shown in Figure 5 with red circles. It indicates that ultracold Rydberg atoms for n > 42 states gradually experience increasing attractive van der Waals interactions, and even stronger dipole-dipole interaction between nP states with decreasing of atomic separation, and simultaneously
Figure 5. The interaction induced total collisional loss rate coefficient βT of nP3/2 with and without BBR.
enhance the interactions between initial nP and nearby dipole-allowed states at the longer delay time that the BBR occurred. Because the BBR has occurred before collision induced by attractive van der Waals interactions and stronger dipole-dipole interaction between nP states, the above analysis is suitable to larger Rydberg atomic separation or the lower Rydberg density case. For high Rydberg density, the smaller initial Rydberg atomic separation will lead to strong attractive van der Waals interactions between initial nP state and the rapid decrease of Rydberg atomic separation, then the attractive van der Waals interactions will quickly evolve into stronger dipole-dipole interactions before BBR transition occurrence, we do not consider the contribution of interactions between initial nP and nearby dipole-allowed states, because the Rydberg atoms have been collision and ionization under initial attractive van der Waals interactions and latter stronger dipole-dipole interactions between initial nP state at shorter delay time than that of BBR transition occurrence.
The attractive interaction case have higher ionization probabilities than that of the repulsive case due to their different collisional mechanisms, in actual experiment, the collision ionization characteristics have significantly difference for n > 42 and n < 42 at longer delay time  . As an example, we estimate the Cs 40P collisional ionization rate using our above model, and compare with the measured results from the reference  . For experimental results of 40P, the density of produced ions varies linearly with initial Rydberg atoms density, the ratio is about 2%. First, the ionization maybe come from BBR photoionization and collisions between the captured electrons and Rydberg atoms  , the calculated corresponding ionization rates of 40P are about 660 s−1  and 600s−1  respectively, which are far smaller than the experimental observed value of 2400 s−1. Second, we consider the interaction induced the collisional ionization rate, the initial atomic separation value R is about 3 μm for the Rydberg density
r = 2.5 × 109 cm−3 by , which will give rise to about 10% of the
Rydberg atoms close together  , the theoretical value of the total collisional loss rate coefficient for 40P is about 1.4 × 10−5 cm3/s using our model, and the corresponding collisional loss rate is about 3.5 × 104 s−1. We can extract the ionization rate is about 700 s−1 with ionization ratio of 2%, even for higher Rydberg density of 4.0 × 109 cm−3, the ionization rate is only about 1120 s−1, which is also smaller than the experimental value. For n < 42, we only consider the contribution of the initial Rydberg atomic velocity to the atomic separation, and think that the weak interactions between Rydberg atoms make atoms separation almost unchanged before BBR transition occurrence. When some close atoms pairs exist at the beginning, they will repeatedly approach each other on long delay time, but the atoms redistributed to nearby dipole-allowed states by BBR will have the chance to collide with these repelled and coming closer atoms, the smaller interatomic distance will give rise to faster collisions due to their stronger dipole-dipole interactions  , which could contribute to an increased collisional loss rate, so the van der Waals interaction play an important role in collision process for repulsive case. If we consider the effect of initial van der Waals interaction and the initial atoms velocity on interatomic distance, for initial Rydberg density of 4.0 × 109 cm3/s, the distance between 40P Rydberg atoms under initial interatomic interaction will increase 2 μm on 10 μs delay time, and a fraction of these repelled atoms may closer to redistributed atoms in nearby dipole-allowed states, compared to their initial atomic separation, which will lead to the faster collisions and ionizations under their stronger dipole-dipole interactions, the calculated ionization rate is about 2240 s−1 with ionization ratio of 2%. In additional, the stray fields from ions can transfer more initial Rydberg atoms to nearby attractive dipole-coupled Rydberg state, and tune their interaction strengths by potential curve determining atomic collisional properties; the present of electric field also can accelerate the impact ionization of Rydberg-electron.
We in detail analyze the different collisional loss mechanisms of Cs nP Rydberg atoms, present simple analytical formula to estimate these collisional loss rate coefficients for the van der Waals and dipole-dipole case, respectively, estimate Cs 40P collisional ionization rate using our above model, and compare it with the previous measured results, which have the important guiding significance for analyzing the experimental results.
This work was supported by the Fundamental Research Fund Project of NIM (Grant No. 27-AKY1705) and the Exploration and Innovation Project of NIM (Grant No. 27-AKYCX1604).