arxiv: v1 [cond-mat.quant-gas] 29 Aug 2012

Transcription

1 Radio-freuency spectroscopy of wealy bound molecules in spin-orbit coupled atomic Fermi gases arxiv:8.58v cond-mat.uant-gas 9 Aug Hui Hu, Han Pu, Jing Zhang, Shi-Guo Peng, and Xia-Ji Liu ARC Centre of Excellence for Quantum-Atom Optics and Centre for Atom Optics and Ultrafast Spectroscopy, Swinburne University of Technology, Melbourne, Australia Department of Physics and Astronomy and Rice Quantum Institute, Rice University, Houston, TX 775, USA State Key Laboratory of Quantum Optics and Quantum Optics Devices, Institute of Opto-Electronics, Shanxi University, Taiyuan 6, P. R. China State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 7, P. R. China (Dated: May, 8) We investigate theoretically radio-freuency spectroscopy of wealy bound molecules in an ultracold spin-orbit-coupled atomic Fermi gas. We consider two cases with either eual Rashba and Dresselhaus coupling or pure Rashba coupling. The former system has been realized very recently at Shanxi University Wang et al., arxiv:.887 and MIT Cheu et al., arxiv:5.8. We predict realistic radio-freuency signals for revealing the uniue properties of anisotropic molecules formed by spin-orbit coupling. PACS numbers:.75.ss,.75.hh, 5..F, d I. INTRODUCTION The coupling between the spin of electrons to their orbital motion, the so-called spin-orbit coupling, lies at the heart of a variety of intriguing phenomena in diverse fields of physics. It is responsible for the wellnown fine structure of atomic spectra in atomic physics, as well as the recently discovered topological state of matter in solid-state physics, such as topological insulators and superconductors,. For electrons, spinorbit coupling is a relativistic effect and in general not strong. Most recently, in a milestone experiment at the National Institute of Standards and Technology (NIST), synthetic spin-orbit coupling was created and detected in an atomic Bose-Einstein condensate (BEC) of 87 Rb atoms. Using the same experimental techniue, noninteracting spin-orbit-coupled Fermi gases of K atoms and 6 Li atoms have also been realized, respectively, at Shanxi University and at Massachusetts Institute of Technology (MIT) 5. These experiments have paved an entirely new way to investigate the celebrated effects of spin-orbit coupling. Owing to the high-controllability of ultracold atoms in atomic species, interactions, confining geometry and purity, the advantage of using synthetic spin-orbit coupling is apparent: (i) The strength of spin-orbit coupling between ultracold atoms can be made very strong, much stronger than that in solids; (ii) New bosonic topological states that have no analogy in solid-state systems may be created; (iii) Ultracold atoms are able to realize topological superfluids that are yet to be observed in the solid state; (iv) Strongly correlated topological states can be readily realized, whose understanding remains a grand Electronic address: xiajiliu@swin.edu.au challenge. At the moment, there has been a flood of theoretical wor on synthetic spin-orbit coupling in BECs 6 8 and atomic Fermi gases 9, addressing particularly new exotic superfluid phases arising from spinorbit coupling 8, 9,,. In this paper, we investigate theoretically momentumresolved radio-freuency (rf) spectroscopy of an interacting two-component atomic Fermi gas with spin-orbit coupling. The interatomic interactions can be easily manipulated using a Feshbach resonance in K or 6 Li atoms. It is nown that wealy bound molecules with anisotropic mass and anisotropic wave function (in momentum space) may be formed due to spin-orbit coupling 9,,. Here, we aim to predict observable rf signals of these anisotropic molecules. Our calculation is based on the Fermi s golden rule for the bound-free transition of a stationary molecule. We consider two inds of spinorbit coupling: () the eual Rashba and Dresselhaus coupling, λ x σ y, which has been realized experimentally at Shanxi University and MIT 5, and () the pure Rashba coupling, λ( y σ x x σ y ),which is yet to be realized. Here σ x and σ y are the Pauli matrices, x and y are momenta and λ is the strength of spin-orbit coupling. The latter case with pure Rashba spin-orbit coupling is of particular theoretical interest, since molecules induced by spin-orbit coupling exist even for a negative s-wave scattering length above Feshbach resonances 9,,. The paper is structured as follows. In the next section, we discuss briefly the rf spectroscopy and the Fermi s golden rule for the calculation of rf transfer strength. In Sec. III, we consider the experimental case of eual Rashba and Dresselhaus spin-orbit coupling. We introduce first the model Hamiltonian and the single-particle and two-particle wave functions. We then derive an analytic expression for momentum-resolved rf transfer strength following the Fermi s golden rule. It can be written explicitly in terms of the two-body wave function.

2 We discuss in detail the distinct features of momentumresolved rf spectroscopy in the presence of spin-orbit coupling, with the use of realistic experimental parameters. In Sec. IV, we consider an atomic Fermi gas with pure Rashba spin-orbit coupling, a system anticipated to be realized in the near future. Finally, in Sec. V, we conclude and mae some final remars. II. RADIO-FREQUENCY SPECTROSCOPY AND THE FERMI S GOLDEN RULE Radio-freuency spectroscopy, including momentumresolved rf-spectroscopy, is a powerful tool to characterize interacting many-body systems. It has been widely used to study fermionic pairing in a two-component atomic Fermi gas near Feshbach resonances when it crosses from a Bardeen-Cooper-Schrieffer (BCS) superfluid of wealy interacting Cooper pairs into a BEC of tightly bound molecules 6. Most recently, it has also been used to detect new uasiparticles nown as repulsive polarons 7, 8, which occur when impurity fermionic particles interact repulsively with a fermionic environment. The underlying mechanics of rf-spectroscopy is simple. For an atomic Fermi gas with two hyperfine states, denoted as = and =, the rf field drives transitions between one of the hyperfine states (i.e., ) and an empty hyperfine state which lies above it by an energy ω due to the magnetic field splitting in bare atomic hyperfine levels. The Hamiltonian for rf-coupling may be written as, V rf = V ˆ dr ψ (r)ψ (r)+ψ (r)ψ (r), () ψ (r) (ψ (r)) is the field operator which creates an atom in ( ) at the position r and, V is the strength of the rf drive and is related to the Rabi freuency ω R with V = ω R /. For the rf-spectroscopy of wealy bound molecules that is of interest in this wor, a molecule is initially at rest in the bound state Φ B with energy E = ǫ B. Here ǫ B is the binding energy of the molecules. A radio-freuency photon with energyω will brea the molecule and transfer one of the atoms to the third state. In the case that there is no interaction between the state and the spin-up and spin-down states, the final state Φ f involves a free atom in the third state and a remaining atom in the system. According to the Fermi s golden rule, the rf strength of breaing molecules and transferring atoms is proportional to the Franc-Condon factor, F (ω) = Φ f V rf Φ B δ ω ω E f E, () the Dirac delta function ensures energy conservation and E f is the energy of the final state. The integrated Franc-Condon factor over freuency should be + unity, F (ω)dω =, if we can find a complete set of final states for rf transition. Hereafter, without any confusion we shall ignore the energy splitting in the bare atomic hyperfine levels and set ω =. To calculate the Franc-Condon factor E. (), it is crucial to understand the initial two-particle bound state Φ B and the final two-particle state Φ f. III. EQUAL RASHBA AND DRESSELHAUS SPIN-ORBIT COUPLING Let us first consider a spin-orbit-coupled atomic Fermi gas realized recently at Shanxi University and at MIT 5. In these two experiments, the spin-orbit coupling is induced by the spatial dependence of two counter propagating Raman laser beams that couple the two spin states of the system. Near Feshbach resonances, the system may be described by a model Hamiltonian H = H + H int, H = σ ˆ ˆ drψ σ ˆ (r) m ψ σ(r)+ ( dr ψ (r) ΩR eirx ) ψ (r)+h.c. () is the single-particle Hamiltonian and H int = U ˆ drψ (r)ψ (r)ψ (r)ψ (r) () is the interaction Hamiltonian describing the contact interaction between two spin states. Here, ψ σ (r) is the creation field operator for atoms in the spin-state σ, ˆ i is the momentum operator, Ω R is the coupling strength of Raman beams, R = π/λ R is determined by the wave length λ R of two Raman lasers and therefore R is the momentum transfer during the twophoton Raman process. The interaction strength is denoted by the bare interaction parameter U. It should be regularized in terms of the s-wave scattering length a s, i.e., /U = m/(π a s ) m/( ). To solve the many-body Hamiltonian E. (), it is useful to remove the spatial dependence of the Raman coupling term, by introducing the following new field operators ψ σ via ψ (r) = e +irx ψ (r), (5) ψ (r) = e irx ψ (r). (6) With the new field operators ψσ, the single-particle Hamiltonian then becomes, ) H = ˆ (ˆ± R e x dr ψ σ(r) ψ σ (r) m σ + Ω ˆ R dr ψ (r) ψ (r)+h.c., (7)

3 in the first term on the right hand side of the euation we tae + for spin-up atoms and for spindown atoms. The form of the interaction Hamiltonian is invariant, H int = U ˆ dr ψ (r) ψ (r) ψ (r) ψ (r). (8) However, the rf Hamiltonian acuires an effective momentum transfer R e x, V rf = V ˆ dr e irx ψ (r) ψ (r)+h.c.. (9) For later reference, we shall rewrite the rf Hamiltonian in terms of field operators in the momentum space, ( ) V rf = V c c Re x, + H.c., () ψ (r) c ei r and ψ (r) c e i r. Hereafter, we shall denote c and c σ as the field operators (in the momentum space) for atoms in the third state and in the spin-state σ, respectively. A. Single-particle solution Using the Pauli matrices, the single-particle Hamiltonian taes the form, ˆ H = dr ψ (r), ψ (r) ( R +) m ψ (r) +hσ x +λ x σ z, () ψ (r) for convenience we have defined the spin-orbit coupling constant denoted as λ R /m and an effective Zeeman field h Ω R /. This Hamiltonian is euivalent to the one with eual Rashba and Dresselhaus spin-orbit coupling, λ x σ y. To see this, let us tae the second transformation and introduce new field operators Ψ σ (r) via ψ (r) = Ψ (r) iψ (r), () ψ (r) = Ψ (r)+iψ (r). () Using these new field operators, the single-particle Hamiltonian now taes the form, ˆ H = drψ (r),ψ (r) ( R +) m Ψ (r) +λ x σ y +hσ z, () Ψ (r) which is precisely the Hamiltonian with eual Rashba and Dresselhaus spin-orbit coupling. The single-particle Hamiltonian E. () can be diagonalized to yield two eigenvalues E ± = R m + m ± h +λ x. (5) Here ± stands for the two helicity branches. The singleparticle eigenstates or the field operators in the helicity basis tae the form c + = +c cosθ +c sinθ, (6) c = c sinθ +c cosθ, (7) θ = arctan( h +λ x λ x)/h > (8) is an angle determined by h and x. Note that, ( cos θ = + ( sin θ = λ x h +λ x λ x h +λ x ) ), (9). () Note also that the minimum energy of the single-particle energy dispersion is given by E min = R m mλ h mλ = h mλ, () if h < mλ /.

4 B. The initial two-particle bound state Φ B In the presence of spin-orbit coupling, the wave function of the initial two-body bound state has both spin singlet and triplet components 9,,. The wave function at zero center-of-mass momentum, Φ B, may be written as, Φ B = ψ ()c c +ψ ()c c +ψ ()c c +ψ ()c c vac, () C c and c are creation field operators of spin-up and spin-down atoms with momentum and C ψ () + ψ () + ψ () + ψ () is the normalization factor. From the Schrödinger euation (H +H int ) Φ B = E Φ B, we can straightforwardly derive the following euations for coefficients ψ σσ appearing in the above two-body wave function : E ( R ) ψ () = + U m + m +λ x E ( R m + E ( R m λ x m + m E ( R m + m ψ ( ) ψ ( )+hψ ()+hψ (), () ) ψ () = U ψ ( ) ψ ( )+hψ ()+hψ (), () ) ψ () = hψ ()+hψ (), (5) ) ψ () = hψ ()+hψ (), (6) E = ǫ B < is the energy of the two-body bound state. Let us introduce A ǫ B ( R /m+ /m) < and different spin components of the wavefunctions, ψ s () = ψ () ψ (), (7) ψ a () = ψ ()+ψ (). (8) It is easy to see that, h ψ () = ψ a (), (9) A h ψ () = ψ a (), () A ψ a () = λ x A h + ψ s (), () A +h and A λ x A h /A ψ s () = U ψ s ( ). () As reuired by the symmetry of fermionic system, the spin-singlet wave function ψ s () is an even function of the momentum, i.e., ψ s ( ) = ψ s (), and the spin-triplet wave functions are odd functions, satisfying ψ a ( ) = ψ a (), ψ ( ) = ψ () and ψ ( ) = ψ (). The un-normalized wavefunction ψ s () = A λ x /(A h /A ) is given by, h λ x ψ s () = h +λ x + A A A (h +λ x ). () Using E. () and un-normalized wave function ψ s (), the bound-state energy E or the binding energy ǫ B is determined by U ψ s() =, or more explicitly, m π ψ s ()+ m a s =. () Here we have replaced the bare interaction strength U by the s-wave scattering length a s using the standard regularization scheme mentioned earlier. The normalization factor of the total two-body wave function is given by, C = ψ s () + λ x (A h) + λ x (A +h) C. The final two-particle state Φ f. (5) Let us consider now the final state Φ f. For this purpose, it is useful to calculate V rf Φ B = V c + c Re x, Φ B (6)

5 5 and then determine possible final states. It can be readily seen that, V rf Φ B = V C c + Re x, { } ψ () ψ ( )c +ψ () ψ ( )c vac. (7) Rewriting ψ and ψ in terms of ψ s and ψ a as shown in Es. (7) and (8), and exploiting the parity of the wave functions, we obtain a general result valid for any type of spin-orbit coupling, V rf Φ B = C V { } c + Re x, ψ s ()+ψ a ()c + ψ ()c vac. (8) To proceed, we need to rewrite the field operators c and c in terms of creation operators in the helicity basis. For the case of eual Rashba and Dresselhaus spin-orbit coupling, using Es. (6) and (7), we find that c = cosθ c + sinθ c, (9) c = sinθ c + +cosθ c. () Thus, we obtain V rf Φ B = C V c + Re x, s + c + s c vac, () s + = ψ s ()+ψ a ()cosθ + ψ ()sinθ, () s = ψ s ()+ψ a ()sinθ ψ ()cosθ. () E. () can be interpreted as follows. The rf photon breas a stationary molecule and transfers a spin-down atom to the third state. We have two possible final states: () we may have two atoms in the third state and the upper helicity state, respectively, with a possibility of s + /C; and () we may also have a possibility of s /C for having two atoms in the third state and the lower helicity state, respectively. D. Momentum-resolved rf spectroscopy Taing into account these two final states and using the Fermi s golden rule, we end up with the following expression for the Franc-Condon factor, F (ω) = ( s C +δ ω E ) ( + +s δ ω E ), () E ± ǫ B + ( R +) m ± h +λ x + (+ R e x ). (5) m The two Dirac delta functions in E. () are due to energy conservation. For example, the energy of the initial state (of the stationary molecule) is E = ǫ B, while the energy of the final state is (+ R e x ) /(m) for the free atom in the third state and (R + )/(m)+ h +λ x for the remaining atom in the upper branch. Therefore, the rf energy ω reuired to have such a transfer is given by E +, as shown by the first Dirac delta function. It is easy to + chec that the Franc-Condon factor is integrated to unity, F (ω) =. Experimentally, in addition to measuring the total number of atoms transferred to the third state, which is proportional to F(ω), we may also resolve the transferred number of atoms for a given momentum or wave-vector x. Such

6 a momentum-resolved rf spectroscopy has already been implemented for a non-interacting spin-orbit coupled atomic Fermi gas at Shanxi University and at MIT. Accordingly, we may define a momentum-resolved Franc-Condon factor, F ( x,ω) = ( s + C δ ω E ( + )+s δ ω E ), (6) the summation are now over the wave-vector ( y, z ) and we have defined x x + R by shifting the wave-vector x by an amount R. This shift is due to the gauge transformation used in Es. (5) and (6). With the help of the two Dirac delta functions, the summation over may be done analytically. We finally arrive at, Θ(x) is the step function and,± = m F ( x,ω) = m s 8π C + Θ (,+ ) +s Θ (, ), (7) ( ω ǫ B ) ( R + x + x R ± m ) h +λ x. (8) 6 It is understood that we will use = ( x,,+ ) in the calculation of s + and = ( x,, ) in s. We may immediately realize from the above expression that the momentum-resolved Franc-Condon factor is an asymmetric function of x, due to the coexistence of spin-singlet and spin-triplet wave functions in the initial two-body bound state. Moreover, the contribution from two final states or two branches should manifest themselves in the different freuency domain in rf spectra. As we shall see below, these features give us clear signals of anisotropic bound molecules formed by spin-orbit coupling. On the other hand, from E. (7), it is readily seen that once the momentum-resolved rf spectroscopy is measured with high resolution, it is possible to determine precisely s + and s and then re-construct the two-body wave function of spin-orbit bound molecules. E. Numerical results and discussions For eual Rashba and Dresselhaus spin-orbit coupling, the bound molecular state exists only on the BEC side of Feshbach resonances with a positive s-wave scattering length, a s >. Thus, it is convenient to tae the characteristic binding energy E B = /(ma s ) as the unit for energy and freuency. For wave-vector, we use R = mλ/ as the unit. The strength of spin-orbit coupling may be measured by the ratio E λ E B = mλa s, (9) we have defined the characteristic spin-orbit energy E λ mλ / = R /m. Note that, the spin-orbit coupling is also controlled by the effective Zeeman field h = Ω R /. In particular, in the limit of zero Zeeman field Ω R =, there is no spin-orbit coupling term as shown in the original Hamiltonian E. (), although the characteristic spin-orbit energy E λ. Using R and E B as the units for wave-vector and energy, we can write a set of dimensionless euations for the binding energy ǫ B = E, normalization factor C, Franc-Condon factor F (ω) and the momentum-resolved Franc-Condon factor F (ω, x ). We then solve them for given parameters E λ /E B and h/e λ. In accord with the normalization condition F (ω) =, the units for F (ω) and F ( x,ω) + are taen to be /E B and /(E B R ), respectively. F( )... E /E B =.5... E = E B. 6 h =.5E. 5 6 /E B Figure : (color online) Franc-Condon factor of wealy bound molecules formed by eual Rashba and Dresselhaus spin-orbit coupling, in units of E B. Here we tae h = E λ/ or Ω R = R/m and set E λ /E B =.5,, and. The result without spin-orbit coupling is plotted by the thin dashed line. The Inset shows the different contribution from the two final states at E λ /E B =. The one with a remaining atom in the lower (upper) helicity branch is plotted by the dashed (dot-dashed) line. Fig. displays the Franc-Condon factor as a function of the rf freuency at h/e λ =.5 and at several ratios of E λ /E B as indicated. For comparison, we show also the rf line-shape without spin-orbit coupling, F (ω) = (/π) ω E B /ω, by the thin dashed line. In the presence of spin-orbit coupling, the existence of two possible final states is clearly revealed by the two peas in the rf spectra. This is highlighted in the inset for

7 7 E λ /E B =, the contribution from the two possible final states is plotted separately. The main rf response is from the final state with the remaining atom staying in the lower helicity branch, i.e., the second term in the Franc-Condon factor E. (). The two pea positions may be roughly estimated from E. (8) for the threshold freuency ω c± of two branches, ( R + x + ) x R ω c± = ǫ B + m ± h +λ x min (5) With increasing spin-orbit coupling, the low-freuency pea becomes more and more pronounced and shifts slightly towards lower energy. In contrast, the highfreuency pea has a rapid blue-shift. /E B /E B 6 5 (a) E /E B = (c) E /E B = x / R (b) E /E B = (d) E /E B = x / R Figure : (color online) Linear contour plot of momentumresolved Franc-Condon factor, in units of (E B R). Here we tae h = E λ / and consider E λ /E B =.,.5,, and. F( x, ).. /E =..5 h =.5E, E = E B x / R Figure : (color online) Energy distribution curve of the momentum-resolved Franc-Condon factor, in units of (E B R). We consider several values of the rf freuency ω as indicated, under given parameters h = E λ / and E λ = E B... Fig. presents the corresponding momentum-resolved Franc-Condon factor. We find a strong asymmetric distribution as a function of the momentum x. In particular, the contribution from two final states are well separated in different freuency domains and therefore should be easily observed experimentally. The asymmetric distribution of F ( x,ω) is mostly evident in energy distribution curve, as shown in Fig., we plot F ( x,ω) as a function of x at several given freuencies ω. In the experiment, each of these energy distribution curves can be obtained by a single-shot measurement. F( ) h/e = /E B h=e - x / R. 6 /E B Figure : (color online) Zeeman-field dependence of the Franc-Condon factor at E λ = E B. Here we vary the effective Zeeman fields, h/e λ =.,.6,.8 and.. The inset shows the momentum-resolved Franc-Condon factor at h = E λ. We finally discuss the effect of the effective Zeeman field h = Ω R /. Fig. shows how the line-shape of Franc-Condon factor evolves as a function of the Zeeman field at E λ /E B =. In general, the larger Zeeman field the stronger spin-orbit coupling. Therefore, as the same as in Fig., the increase in Zeeman field leads to a pronounced pea at about the binding energy. There is a red-shift in the pea position as the binding energy becomes smaller as the Zeeman field increases. As anticipated, the larger the Zeeman field, the more asymmetric F ( x,ω) becomes. In the inset, we show as an example the contour plot of F ( x,ω) at h/e λ =. IV. RASHBA SPIN-ORBIT COUPLING We now turn to the case with pure Rashba spin-orbit coupling, λ( y σ x x σ y ), which may be realized experimentally in the near future. The single-particle Hamiltonian may be written as, ˆ H = dr ψ (r),ψ (r) S ψ (r), (5) ψ (r) the matrix ( S = R +) /(m) iλ( x i y ) iλ( x +i y ) ( R +) /(m). (5) Here λ is the coupling strength of Rashba spin-orbit coupling, R mλ/, and we have added a constant term

8 8 R /(m) to mae the minimum single-particle energy zero, i.e., E min =. Here ϕ arg( x, y ) is the azimuthal angle of the wavevector in the x y plane. A. Single-particle solution We diagonalize the matrix S to obtain two helicity eigenvalues, E ± = ( R +) ±λ, (5) m x +y and ± stands for the two helicity branches. For later reference, the field operators in the original spin basis and in the helicity basis are related by, c = ) (c + +ieiϕ c, (5) c = ) (ie iϕ c + +c. (55) B. The initial two-particle bound state Φ B In the case of Rashba spin-orbit coupling, the two-body wave function can still be written in the same form as in the previous case, i.e., E. (), and the Schrödinger euation leads to, A ψ () = + U ψ ( ) ψ ( ) λ( y i x )ψ ()+λ( y +i x )ψ (), (56) A ψ () = U ψ ( ) ψ ( )+λ( y i x )ψ () λ( y +i x )ψ (), (57) A ψ () = λ( y +i x )ψ ()+λ( y +i x )ψ (), (58) A ψ () = +λ( y i x )ψ () λ( y i x )ψ (), (59) A E ( R /m+ /m) <. It is easy to show that ψ a () = and A λ ψ s () = U ψ s ( ). (6) A Thus, we obtain the (un-normalized) wavefunction: ψ s () = +, (6) E E + E E and the euation for the energy E, m π = / / + + m a s E E + E E. (6) The spin-triplet wave functions ψ () and ψ () are given by, ψ () = ie iϕ λ ψ s (), (6) E ǫ ψ () = ie +iϕ λ ψ s (), (6) E ǫ ǫ /(m). The normalization factor for the two-body wave function is therefore, C = ψ s () + λ (E ǫ ) C. The final two-particle state Φ f. (65) To obtain the final state, we consider again V rf Φ B. In the present case, we assume that the rf Hamiltonian is given by, V rf = V ( ) c c +c c. (66) Following the same procedure as in the case of eual Rashba and Dresselhaus coupling, it is straightforward to show that,

9 9 V rf Φ B = C V c + ψ s ()c + ψ ()c vac. (67) Using Es. (5) and (55) to rewrite c and c in terms of c + and c, we obtain, V rf Φ B = C V c + c + c + + ieiϕ c vac. (68) E E + E E Therefore, we have again two final states, differing in the helicity branch that the remaining atom stays. The remaining atom stays in the upper branch with probability (C) (E E + ), and in the lower branch with probability (C) (E E ). D. Momentum-resolved rf spectroscopy Using the Fermi s golden rule, we have immediately the Franc-Condon factor, F (ω) = δ(ω E + /) C (ǫ B +E + ) + δ(ω E /) (ǫ B +E ), (69) E ± ǫ B + R m + m ±λ. (7) For Rashba spin-orbit coupling, it is reasonable to define the following momentum-resolved Franc-Condon factor, F (,ω) = δ(ω E + /) C (ǫ B +E + ) + δ(ω E /) (ǫ B +E ), (7) z we have summed over the momentum z. Integrating over z with the help of the two Dirac delta functions, we find that, F (,ω) = m Θ ( ) z+ Θ ( ) z 6π C (ω + R /m+λ ) + z+ (ω + R /m λ ), (7) z z± = m ( ω ǫ B ) ) ( R + ± R. (7) It is easy to see that the threshold freuencies for the two final states are given by, ω c+ = ǫ B + R m, (7) ω c = ǫ B + R m, (75) which differ by an amount of R /(m) = E λ/. Near ω c, we find approximately that F (ω) Θ(ω ω c )/ω. Thus, the lineshape near the threshold is similar to that of a two-dimensional (D) Ferm gas 6. This similarity is related to the fact that at low energy a D Fermi gas with Rashba spin-orbit coupling has exactly the same density of states as a D Fermi gas. E. Numerical results and discussions For the pure Rashba spin-orbit coupling, the molecular bound state exists for arbitrary s-wave scattering length a s 9,,. We shall tae R = mλ/ and E λ mλ / as the units for wave-vector and energy, respectively. With these units, the dimensionless interaction strength is given by /(mλa s ). The spin-orbit

10 effect should be mostly significant on the BCS side with /(mλa s ) <, the bound state cannot exist without spin-orbit coupling. F( ) /(m a s ). Figure 5: (color online) Franc-Condon factor of wealy bound molecules formed by Rashba spin-orbit coupling, in units of E λ. Here we tae /(mλa s) = (dashed line), (solid line), and (dot-dashed line). In the deep BCS limit, /(mλa s), the Franc-Condon factor peas sharply at ω E λ / and becomes a delta-lie distribution. Fig. 5 shows the Franc-Condon factor at three different interaction strengths /(mλa s ) =,, and +. The strong response in the BCS regime (a s < ) or in the unitary limit (a s ± ) is an unambiguous signal of the existence of Rashba molecules. In particular, the rf line-shape in the BCS regime shows a sharp pea at about ω E λ / and decays very fast at high freuency. In Fig. 6, we present the corresponding momentumresolved Franc-Condon factor F (,ω), in the form of contour plots. We can see clearly the different response from the two final states. The momentum-resolved rf spectroscopy is particularly useful to identify the contribution from the final state that has a remaining atom in the upper branch, which, being integrated over, becomes too wea to be resolved in the total rf spectroscopy. Finally, we report in Fig. 7 energy distribution curves of F (,ω) in the unitary limit /(mλa s ) =. We find two sharp peas in each energy distribution curve, arising from the two final states. When measured experimentally, these sharp peas would become much broader owing to the finite experimental energy resolution. /E V. CONCLUSIONS In conclusions, we have investigated theoretically the radio-freuency spectroscopy of wealy bound molecules in a spin-orbit coupled atomic Fermi gas. The wave function of these molecules is greatly affected by spinorbit coupling and has both spin-singlet and spin-triplet components. As a result, the line-shape of the total radio-freuency spectroscopy is ualitatively different from that of the conventional molecules at the BEC- /E /E /E (a) /(m a s ) = (b) /(m a s ) = (c) /(m a s ) = / R. Figure 6: (color online) Contour plot of momentum-resolved Franc-Condon factor of wealy bound molecules formed by Rashba spin-orbit coupling, in units of (E λ R). The intensity increases from blue to red in a logarithmic scale. We consider /(mλa s) = (a), (b), and (c). BCS crossover without spin-orbit coupling. In addition, the momentum-resolved radio-freuency becomes highly asymmetric as a function of the momentum. These features are easily observable in current experiments with spin-orbit coupled Fermi gases of K atoms and 6 Li atoms. On the other hand, from the highresolution momentum-resolved radio-freuency, we may re-construct the two-body wave function of the bound molecules. We consider so far the molecular response in the radiofreuency spectroscopy. Our results should be uantiatively reliable in the deep BEC limit with negligible number of atoms, i.e., in the interaction parameter regime with /( F a s ) >. However, in a real experiment, in order to maximize the spin-orbit effect, it is better to wor closer to Feshbach resonances, i.e., /( F a s ).5. Under this situation, the spin-orbit coupled Fermi gas consists of both atoms and wealy bound molecules, which may strongly interact with each other. Our prediction for the molecular response is still ualitatively valid, with the understanding that there would be an additional pronounced atomic response in the rf spectra. A more indepth investigation of radio-freuency spectroscopy re-.

11 F(, ).8.6 /E = / R Figure 7: (color online) Momentum-resolved Franc-Condon factor of Rashba molecules at /(mλa s) =, shown in the form of energy distribution curves at several rf freuencies as indicated..5 uires complicated many-body calculations beyond our simple two-body picture pursued in the present wor. Acnowledgments We would lie to than Zeng-Qiang Yu and Hui Zhai for useful discussions. HH and XJL are supported by the ARC Discovery Projects (Grant No. DP985 and No. DP9867) and the National Basic Research Program of China (NFRP-China, Grant No. CB95). HP is supported by the NSF and the Welch Foundation (Grant No. C-669). JZ is supported by the NSFC Project for Excellent Research Team (Grant No. 66) and the NFRP-China (Grant No. CB96). SGP is supported by the NSFC-China (Grant No. ) and the NFRP-China (Grant No. CB96). X.-L. Qi and S.-C. Zhang, Physics Today 6, (). M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 8, 5 (). Y.-J. Lin, K. Jiménez-García, and I. B. Spielman, Nature (London) 7, 8 (). P. Wang, Z.-Q. Yu, Z. Fu, J. Miao, L. Huang, S. Chai, H. Zhai, and J. Zhang, Phys. Rev. Lett. 9, 95 (). 5 L. W. Cheu, A. T. Sommer, Z. Hadzibabic, T. Yefsah, W. S. Bar, and M. W. Zwierlein, Phys. Rev. Lett. 9, 95 (). 6 T. D. Stanescu, B. Anderson, and V. Galitsi, Phys. Rev. A 78, 66 (8). 7 J. Larson and E. Sjövist, Phys. Rev. A 79, 67 (9). 8 C. Wang, C. Gao, G.-M. Jian, and H. Zhai, Phys. Rev. Lett. 5, 6 (). 9 C. Wu, I. Mondragon-Shem, and X.-F. Zhou, Chin. Phys. Lett. 8, 97 (). T.-L. Ho and S. Zhang, Phys. Rev. Lett. 7, 5 (). X. Q. Xu and J. H. Han, Phys. Rev. Lett. 7, (). S. Sinha, R. Nath, and L. Santos, Phys. Rev. Lett. 7, 7 (). H. Hu, B. Ramachandhran, H. Pu, and X.-J. Liu, Phys. Rev. Lett. 8, (). H. Hu and X.-J. Liu, Phys. Rev. A 85, 69 (). 5 B. Ramachandhran, B. Opanchu, X.-J. Liu, H. Pu, P. D. Drummond, and H. Hu, Phys. Rev. A 85, 66 (). 6 R. Barnett, S. Powell, T. Gra, M. Lewenstein, and S. Das Sarma, Phys. Rev. A 85, 65 (). 7 Y. Deng, J. Cheng, H. Jing, C.-P. Sun, and S. Yi, Phys. Rev. Lett. 8, 5 (). 8 Q. Zhu, C. Zhang, and B. Wu, eprint arxiv: J. P. Vyasanaere, V. B. Shenoy, Phys. Rev. B 8, 955 (). M. Isin and A. L. Subas, Phys. Rev. Lett. 7, 5 (). S. L. Zhu, L. B. Shao, Z. D. Wang, and L. M. Duan, Phys. Rev. Lett. 6, (). Z. Q. Yu and H. Zhai, Phys. Rev. Lett. 7, 955 (). H. Hu, L. Jiang, X.-J. Liu, and H. Pu, Phys. Rev. Lett. 7, 95 (). L. Jiang, X.-J. Liu, H. Hu, and H. Pu, Phys. Rev. A 8, 668 (). 5 M. Gong, S. Tewari, and C. Zhang, Phys. Rev. Lett. 7, 95 (). 6 L. Han and C. A. R. Sá de Melo, Phys. Rev. A 85, 66(R) (). 7 X.-J. Liu, L. Jiang, H. Pu, and H. Hu, Phys. Rev. A 85, 6 (). 8 X. Cui, Phys. Rev. A 85, 75 (). 9 X.-J. Liu and H. Hu, Phys. Rev. A 85, 6 (). L. He and X.-G. Huang, Phys. Rev. Lett. 8, 5 (). W. Yi and W. Zhang, arxiv:.676. C. Chin, R. Grimm, P. S. Julienne, and E. Tiesinga, Rev. Mod. Phys. 8, 5 (). C. Chin and P. S. Julienne, Phys. Rev. A 7, 7 (5). C. Chin, M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, J. Hecer Denschlag, and R. Grimm, Science 5, 8 (). 5 C. H. Schunc, Y. Shin, A. Schirotze, and W. Ketterle, Nature (London) 5, 79 (8). 6 Y. Zhang, W. Ong, I. Araelyan, J. E. Thomas, Phys. Rev. Lett. 8, 5 (). 7 C. Kohstall, M. Zaccanti, M. Jag, A. Trenwalder, P. Massignan, G. M. Bruun, F. Schrec and R. Grimm, Nature (London) 85, 65 (). 8 M. Koschorrec, D. Pertot, E. Vogt, B. Fröhlich, M. Feld and M. Köhl, Nature (London) 85, 69 ().