Department of Physics, Anhui Normal University, Wuhu 241000, China
† Corresponding author. E-mail:
fengbf@mail.ahnu.edu.cn
1. IntroductionCold atomic and molecular species are extremely versatile systems to investigate a wealth of fascinating phenomena, ranging from precision measurement[1] and quantum-controlled chemistry,[2] to condensed-matter physics.[3] Ultracold atoms and molecules are usually confined in electromagnetic or magneto–optical traps, and collision cross sections of ultracold species change dramatically as the external field varies through a resonance. Tunable Feshbach resonances (FR) thus are powerful tools to control the interaction and scattering properties of ultracold species.[4] Magnetic FRs have been used to realize Bose–Einstein condensates of alkali metal atoms and create associated ultracold molecular species by coherently linking ultracold atoms.[5] In addition, FRs provide experimental access to interesting three-body Efimov states,[6] loosely bound molecules,[7] and strongly correlated many-body states of matter.[8]
The coupling between different internal states of the two-body collision system can cause interaction between different scattering channels. A Feshbach resonance[4, 7] occurs when a bound state in a closed-channel exists at the same energy as the dissociation energy threshold of collision pairs. Because the bound states and collision pairs have, in general, different electric or magnetic moments, an external field can be used to tune both their energies and corresponding energy differences, and thus tune the FRs. Magnetic FR has been the major tool in various experiments of ultracold atomic and molecular gases.[9, 10] A second way of tuning the FR is to apply a dc-electric field.[11–13] The electric field induces an electric dipole–dipole interaction between heteronuclear colliding atoms. This anisotropic interaction couples the states of different orbital angular momenta, and changes the width of FR resonance to some degree. A recent theoretical work[14] has demonstrated the possibility of inducing and modifying the position and width of narrow FRs by intense nonresonant radiation. The authors show that the resonance widths can be increased by three orders of magnitude up to a few Gauss. The use of combined electric and magnetic fields to control interatomic and intermolecular interactions has also been suggested. It has been demonstrated that dynamics of molecules in a regime of ultracold temperature may be sensitive to the magnitude of an applied field,[15–19] and the combination of electric and magnetic fields may be used to control both the position and width of FRs independently, leading to total control over the character of ultracold collisions.[20] In addition, the magnetic control of collision dynamics is limited to paramagnetic molecules; therefore, the application of electric fields may expand the scope of studies of correlative phenomena in ultracold polar gas.
Recently, we investigated the collisional dynamics of 3He atom and PH(
molecule in the presence of a magnetic field.[21] The work has its background of Helium buffer-gas cooling and magnetic trapping experiments, which was first performed successfully for CaH molecule and later for other molecules NH, CaF, and MnH, etc.[30] Our work shows that the PH molecule is a promising candidate for this type of experiment. In the present paper, we will explore the effects of an external electric field on magnetically induced FRs in the 3He–PH collision system. In Section 2 we briefly describe the methods we use to perform the bound state and scattering calculations. In Section 3, we discuss the tuning of magnetic FRs by a dc-electric field and shall conclude in Section 4.
2. Theoretical approachTo precisely estimate the position of a near zero-energy Feshbach resonance, we first compute the bound state energies of the collision complex as functions of external fields. The method used in this work is analogous to that described in our previous papers.[21–24] Here, we extend it to include the electric and magnetic field effect. The Hamiltonian for the atom–diatom complex has the following form:
![](cpb_26_1_013402/cpb_26_1_013402002.gif) | (1) |
where
μ is the reduced mass of the collisional system,
R is the distance between the center of mass of the diatomic molecule and the atom, and
is the orbital angular momentum for the collision. The
is the interaction potential which vanishes as
. The
denotes the Hamiltonian for a
molecule interacting with electric and magnetic fields and is expressed as
![](cpb_26_1_013402/cpb_26_1_013402008.gif) | (2) |
where
is the rotational constant of the molecule in the ground vibration state,
and
are the rotational and spin angular momenta,
and
are the spin–rotation and spin–spin interaction terms,
and
are the electric and magnetic fields, and
is the electric dipole moment of the molecule. We assume that the electric and magnetic field vectors are both oriented along the
Z space fixed axis.
In the presence of parallel electric and magnetic fields, the total angular momentum J and the total parity which is given by
are not the good quantum numbers, the only rigorously good quantum number is the total angular momentum projection
, where mn, ms, and ML are the projections of
,
, and
individually. There are several basis sets that could be used to expand the eigenfunctions of Eq. (1). We use an uncoupled basis set
.[25] The radial basis
are determined as numerical eigensolutions to a vibrational Schrödinger equation using the one-dimensional potential
cut at steady structure (
.[21] The matrix elements of the interaction with an electric field can be calculated as follows:[18]
![](cpb_26_1_013402/cpb_26_1_013402026.gif) | (3) |
The other angular matrix elements of operators in Eqs. (
1) and (
2) can be found in Ref. [
25] and the radial matrix elements are calculated by numerical integration. The energies of bound states and corresponding eigenfunctions of the complex are determined from solution of the corresponding secular Eq. (
1). The calculations are performed for each of given
MJ and electric and magnetic fields.
As for the scattering calculations, we adopt the method for atom–molecule collisions in combined electric and magnetic fields developed by Tscherbul and Krems.[17, 26, 27] Because the asymptotical Hamiltonian in Eq. (2) is not diagonal in the uncoupled angular basis set
,
then cannot be used to describe the asymptotic states of the dressed molecule in the presence of the field. We use instead a transformative basis set
,[28] which diagonalizes the monomer Hamiltonian (2) as
![](cpb_26_1_013402/cpb_26_1_013402030.gif) | (4) |
where
is the energy level of the dressed diatomic molecule in the external field. The field removes the degeneracy in
and each energy level
correlating with a unique value of
mj labeled by
, where
j is the total angular momentum of diatomic molecule,
. For given values of
MJ and
mj,
. The basis set describing the collision process is obtained by including the possible values of the quantum number
L for each value of
α. With the basis set, the close-coupling equations are
![](cpb_26_1_013402/cpb_26_1_013402037.gif) | (5) |
The cross sections for elastic and inelastic scattering are computed from the
S matrix as follows:
![](cpb_26_1_013402/cpb_26_1_013402038.gif) | (6) |
The coupled equations are propagated on radial grid points using Johnson's log-derivative algorithm.
[29] 3. Applications to the He–PH collisionsWe now focus on the electric-field modified Feshbach resonance in ultracold collisions of atom and molecule, taking the 3He–PH collision system as an example. The calculations are performed using the same potential–energy surface and parameters as presented in Ref. [21]. The positions of Feshbach resonances can be identified by the bound state energies of the collision complex as functions of an external field. Figure 1 shows the bound state energies of 3He–PH for
and
as a function of the magnetic field. The energies are relative to the ground state energy of PH
which is −0.088 cm
. When an external field is applied, each level at zero field splits into components that can be labeled with the approximate quantum numbers
, ±1, and L. As expected, the Zeeman splits are approximately linear for He–PH in the range of fields studied. This is because for a rotational n = 0 state, the coupling arise from the off-diagonal term of
and spin–spin interaction that mix the rotational excited levels, but the spacing (∼16.8 cm
between n = 0 and n = 1 levels for PH is larger than the Zeeman split tuned by a magnetic field. While an electric field is further added, the energies of bound states will shift slightly, but no new energy level emerges.
In Fig. 1, we also show the rotational n = 0 energies of a monomer PH molecule as functions of the magnetic field as dotted lines. When the collision pairs approach each other with near-zero kinetic collision energy, the total collision energy will be approximately equal to the bound state energy of complex, which is the energy constraint for the presentation of a Feshbach resonance. Another constraint is the same conserved quantity of collision channels as the relative bound states.[25] In the case of parallel electric and magnetic fields, MJ is the only conserved quantity. At low collision energies, the cross sections are dominated by incoming s-wave, L = 0, and
. It means that the
of a monomer molecule should be equal to the MJ of the bound state of complex. Based on this analysis, we can estimate the positions of FR shown in Fig. 1 by circles and labels of “Res”. A detailed scan across a small range of fields in the vicinity of crossings and careful extrapolation of the bound state energies with external fields will give more precise positions of FRs. The scattering calculations are then performed at surrounding values of electric and magnetic fields corresponding to these positions respectively.
3.1.
, Res3The FR at
, Res3 has been reported in our previous work[21] with only the magnetic field applied. In the ultracold collision, we want s-wave resonance, with L = 0 in the incoming channel. For 3He colliding with PH (
) at ultralow kinetic energy, only elastic scattering can occur. Scattering into the
channels is strongly suppressed by the centrifugal barriers. Fitting the scattering length, we obtain the resonance position
G, width
G and background scattering length
.
Figure 2 shows the magnetically tuned FR at three values of electric fields. It can be seen that electric fields can shift the position of FRs significantly. When the electric field increases, the position moves to the low magnetic-field direction (corresponding to a shift of the associated bound state to high energy). The width of resonance changes trivially and is still in order of
G. The magnetic field of this degree of accuracy is difficult to be produced in an experiment. The FR with a width of this order of magnitude is thus unlikely to be tuned by the magnetic field. From the energy point of view, a narrow FR width of magnetic field implies a narrow energy interval. The electric field energy density is customarily smaller than that of the magnetic field, so for a given energy interval, there is a relatively larger field interval in the electric field than the magnetic field. This fact can be used to tune the narrow magnetic FRs by the electric field. Figure 3 displays the cross section as a function of the electric-field. The magnetic field is fixed at 10917.1 G, which is in the vicinity of the magnetic FR. It can be seen that the electric field can also be used to achieve the control of the FR. The resonance width attains 0.14 kV/cm in the electric field scale. It thus is preferable that, for a narrow magnetic FR, we can first adjust the magnetic field to a position near the FR, and then finely tune the electric field to achieve resonance.
3.2.
, Res1, and Res2When only a magnetic field is applied, no FRs take place at the intersection points
, Res1 and
, Res2. This is because the magnetic field cannot couple the rotational states of the PH monomer molecule with different parities,[25] the incoming s-wave channel (L = 0, n = 0,
) has parity of −1, but the bound states of the complex relative to the resonances are (L = 3,
) with parity of +1. The lowest Zeeman energy level of the monomer molecule labeled with (n = 0,
) has a main ingredient of (n = 0,
,
). The electric field can produce coupling of the states of different parities (see Eq. (3)
sign) with
, and cause the dynamic mixing between (n = 0,
,
) and (n = 1,
,
) rotational states. Therefore, the electric field will induce FRs which cannot occur with only the magnetic field applied. Figure 4 shows the scattering length of FRs at
, Res1, and Res2 induced by an electric field of 30 kV/cm. The position of FR shifts to the high magnetic-field direction as the electric field increases, contrary to that at
, Res3. The position of the resonance shifted to the high magnetic field or low magnetic field direction, is determined according to whether the discrete bound state crosses resonance from above or below. The electric-induced FRs have first been demonstrated in heteronuclear mixtures of bi-alkali atomic gases.[11, 12, 15, 20] The anisotropic interaction between the instantaneous dipole moment of a heteronuclear collision complex with the external electric field, couples the states of different orbital angular momenta with
, and results in multiple resonances. The mechanism is thus different from that described here.
3.3.
, Res1In Fig. 1(b), the monomer molecule threshold state (n = 0,
) crosses two bound states of complex, one is (L = 3,
) and the other is (L = 4,
). Doing scattering calculations in the vicinity of two crossings, we only find a resonance that is marked with “Res1”, and do not find resonance at the other crossing point, possibly due to the suppression by a centrifugal energy barrier. This resonance has not been reported in our previous work.[21]
Figure 5 illustrates the cross section as functions of the magnetic field near the crossing point fixed the electric field at zero and 120 kV/cm, respectively. For He colliding with PH (
), elastic and inelastic scattering are possible. The incoming channel is (L = 0, n = 0,
), the outcoming channel may be either elastic scattering channels (L = 0–8, n = 0,
), or inelastic channels (L = 0–8, n = 0,
). In ultracold collision energy, the elastic channels with
make no significant contribution.
The following characteristics can be seen. (I) When the electric field is zero, the resonance phenomenon occurs only in the inelastic cross section, while the elastic cross section is essentially the same size of about 59.5 Å2. (II) The electric field induces the resonance in the elastic cross-section, and the resonance position shifts toward a low magnetic field with the electric field increasing. The resonance peak is antisymmetric. (III) The affect of the electric field on the resonance amplitude of the inelastic cross section is very obvious, for example, from about 1 Å2 at E = 0 to
Å2 at E = 120 kV/cm. (IV) The line shape of the inelastic cross section is symmetrical.
4. ConclusionWhen the collision energy of free atom and molecule coincides with that of a bound state of complex, the Feshbach resonance phenomenon occurs. In this work, we study the effect of the electric field on the magnetic Feshbach resonance in atom–molecule ultracold collision taking the He–PH system as an example. We begin with the examination of bound state energies of the collision complex with the variation of the electromagnetic field, and find the energy intersection points of the collision system and the bound state of the complex. This method is especially effective for determining the very narrow resonance position as in the case of most atom–molecule ultracold collisions. The close coupled scattering calculations are then performed by scanning the electromagnetic field around these cross-over points carefully. Results show as follows that: (i) the electric field can induce Feshbach resonance, which cannot occur when only a magnetic field is applied; (ii) the electric field can shift the position of the magnetic Feshbach resonance, and can change the amplitude of resonance significantly; (iii) the width of resonance in electric field scale is relatively larger than that in magnetic field scale. For narrow magnetic FR, the control by an electric field may be preferable.