# Quantum phases of a two-dimensional dipolar Fermi gas

## Abstract

We examine the superfluid and collapse instabilities of a quasi two-dimensional gas of dipolar fermions aligned by an orientable external field. It is shown that the interplay between the anisotropy of the dipole-dipole interaction, the geometry of the system, and the -wave symmetry of the superfluid order parameter means that the effective interaction for pairing can be made very large without the system collapsing. This leads to a broad region in the phase diagram where the system forms a stable superfluid. Analyzing the superfluid transition at finite temperatures, we calculate the Berezinskii–Kosterlitz–Thouless temperature as a function of the dipole angle.

###### pacs:

03.75.Ss,03.75.HhTrapped ultracold gases are increasingly being used to simulate solid-state systems, where clear experimental signatures of theoretical predictions are often lacking Bloch (). A limitation of these gases is that the interactions are typically -wave Reviews () whereas order parameters in solid state systems exhibit richer - and -wave symmetries. Recent progress in the production of trapped, cold dipolar gases Lahaye (); Koch (); Ospelkaus () promises to change this since the dipole-dipole interaction is long-range and anisotropic. Importantly, the interaction in fermionic heteronuclear Feshbach molecules can be large, with electric dipole moments on the order of a Debye Kotochigova (); Ospelkaus (); Weber (). This opens the door to experimentally reaching the superfluid phase. One complicating feature is that the gases tend to become unstable when interactions are sufficiently strong and attractive.

In this letter, we propose confining a gas of fermionic dipoles of mass in the 2D plane by a harmonic trapping potential . For where is the Fermi energy, determined by the 2D density , the system is effectively 2D. The dipoles are aligned by an external field (see Fig. 1), subtending an angle with respect to the – plane. The interaction between two dipole moments separated by is given by . Here, for electric dipoles and is the angle between and . The strength of the interaction is parametrized by the dimensionless ratio . Since we consider identical fermions at low temperatures , additional short range interactions are suppressed. We show that in this configuration, the effective interaction for pairing has -wave symmetry and can be made large without the system collapsing. This makes the 2D system of dipoles a promising candidate to study quantum phase transitions and pairing with unconventional symmetry.

A key point is that the effective dipole-dipole interaction in the – plane can be tuned by changing the angle . This gives rise to several interesting effects. For , the interaction is repulsive and isotropic in the – plane and the dipoles are predicted to undergo a quantum phase transition to a crystalline phase for at Buchler (). As is decreased, the interaction becomes anisotropic in the – plane and for , an attractive sliver appears along the axis as illustrated in Fig. (1). This gives rise to two competing phenomena: superfluidity and collapse Koch (). The resulting phase diagram is shown in Fig. (2) and is the main focus of this letter.

We first examine the critical angle below which the attractive part of the interaction overcomes the Fermi pressure and the gas is unstable towards collapse. The instability is identified with a negative value of the inverse compressibility . We analyze the stability of the dipolar Fermi gas using the normal phase energy density with

(1) |

the kinetic energy density. We allow for the possibility that the Fermi surface can deform Miyakawa () due to the anisotropy of the interaction, by writing the Fermi distribution function as . The direct and exchange energies are ()

(2) |

and

(3) |

with , the dipolar field operator, and the size of the gas. In the superfluid phase, there is a contribution to the energy arising from the condensation energy which is small for . We ignore this contribution in the following.

In the 2D limit , all dipoles reside in the lowest harmonic oscillator level in the -direction. We then have with

(4) |

Here, is the Bessel function of first order, , and . Evaluating the integrals in (2) and (3), we find

(5) |

where

(6) |

is a dimensionless function of the Fermi surface deformation parameter and the dipole angle , and . Equations (1), (5), and (6) give the energy as a function of , , and . For a given coupling strength and angle , is found by minimizing . In general, this has to be done numerically. However, for small Fermi surface deformation we can expand (6) in the small parameter . The resulting integrals are straightforward and using we obtain

(7) |

Note that as expected.

Since , it is now straightforward to evaluate the compressibility . The resulting region of the collapse instability is shown in Fig. 2. There is a large region where the system is stable even though the interaction is strong () and partly attractive []. This is most easily understood by ignoring the Fermi surface deformation (), which yields for . In this limit, the interaction energy dominates and the gas collapses when the dipole-dipole interaction is attractive in more than half the – plane, i.e. for . Fermi surface deformation effects are significant however; in Fig. 3 we show how deviates substantially from along the critical line for collapse.

The stability of the 2D system for a strong, attractive interaction makes it a promising candidate to observe superfluid pairing with unconventional symmetry. To investigate this, we solve the effective 2D BCS gap equation at , derived from the usual 3D BCS pairing Hamiltonian using . This ansatz results in a 2D gap equation in terms of an effective interaction , where with is the lowest harmonic oscillator wavefunction for the reduced mass . Using the rescaled momentum to describe pairing about the deformed Fermi surface, the 2D gap equation becomes

(8) |

with , , and the chemical potential. is the Fourier transform

(9) |

of , expressed in terms of the scaled momenta ; i.e. . Note that only components contribute to the Fourier series since the order parameter is antisymmetric. The Pauli exclusion principle therefore cancels the short range divergence in the dipole-dipole interaction, making the Fourier transform (9) finite even in the 2D limit.

Using the expansions and in (8), where is the angle between and the axis, we find numerically that for and all angles . This shows that to a very good approximation the gap is -wave, , and we can replace with . Using this in (8), it reduces to

(10) |

where .

The dimensionless effective pairing interaction corresponds to the dipole-dipole interaction averaged over the deformed Fermi surface, weighted by the -wave symmetry of the order parameter. After a straightforward but lengthy calculation, we find

(11) | |||||

Here, and are dimensionless functions given by

(12) |

and for and ; for one simply swaps and . and are the complete elliptic integrals of first and second kind respectively and . We have . For , one has showing that the 2D effective interaction has a high energy cut-off for . The last line in (11) is exact for .

We now analyze the gap equation (10) as the dipole angle is adjusted. Since for all , the gap is zero unless . The resulting superfluid region found by evaluating the sign of is shown in Fig. 2. Numerically, we find that the Fermi surface deformation has little effect on the effective interaction for pairing as can be seen in Fig. 2. Ignoring Fermi surface deformation effects, we find from (11) that the system is superfluid for . The critical angle for superfluidity is therefore to a good approximation independent of the interaction strength and is purely determined by geometry.

Crucially, as decreases from , the system becomes unstable toward pairing before the gas collapses, i.e. , and there is a significant region in the phase diagram where the system is superfluid yet stable. Because the -wave pair wavefunction is predominantly oriented along the axis where the interaction is maximally attractive, the effective pairing interaction (11) is stronger than the “bare” 2D interaction in (5) that determines the collapse instability. For this reason, even for strong interactions, there remains a window where the system is a stable superfluid.

To examine the strength of the pairing interaction further, we make use of the fact that pairing occurs primarily at the Fermi surface, and in Fig. 3 we plot from (11) along the critical line for collapse. This gives the largest possible attractive pairing interaction, before the system collapses. We see that the effective pairing interaction increases monotonically with and can become very large. Thus, one can produce a strongly paired gas without the system collapsing. This should be compared with the 3D trapped case where recent results indicate that the system is superfluid and stable only in a narrow region in phase space where the pairing is relatively weak Miyakawa (); Baranov04 (). Solving the gap equation in the weak coupling regime yields for (i.e. ) with an ultraviolet cut-off. This shows that the superfluid phase transition is infinite order in the sense that for all .

At finite temperatures, long-range order is destroyed by phase fluctuations in 2D, and . The superfluid density remains finite, however, describing a Berezinskii-Kosterlitz-Thouless (BKT) superfluid of bound vortex-antivortex pairs KTB (). The critical temperature for this phase is the temperature at which the free energy of a single unbound vortex vanishes ChaikinLubensky (),

(13) |

Here, is the average of the diagonal components of the superfluid mass density tensor . In estimating the energy of a single vortex of radius in a 2D box with sides of length , we have assumed that the contribution from the off-diagonal component due to the anisotropy of the vortex velocity field Yi06 () is small. Generalizing the usual expression for Leggett () to allow for the effects of Fermi surface deformation, one finds

(14) |

As shown in the inset of Fig. 4, rotational symmetry is broken by the external field and . The component of the superfluid density tensor is suppressed since the quasiparticle spectrum is gapless in the direction perpendicular to the field.

Since depends on , (13) has to be solved self-consistently. We determine from (14) using a -dependent gap calculated by including the usual Fermi functions in the gap equation. The critical temperature determined this way is plotted in Fig. (4) for . We also plot the mean field transition temperature obtained from the gap equation. Phase fluctuations suppress below . For weak coupling, . Since and , we have . This suggests that one can cross the critical temperature by adiabatically expanding the gas keeping constant so that decreases Petrov (). For stronger interactions, quickly approaches its upper bound .

Let us consider the experimental requirements for observing the effects discussed in this letter. Recently, a gas of K-Rb polar molecules was created with a density of cm at a temperature Ospelkaus (). These molecules have dipole moments of Debye in their vibrational ground state. Writing , where Debye is the dipole moment in Debyes, m is the interparticle spacing measured in , and is the mass of the dipoles in atomic mass units, the experimental parameters of Ref. Ospelkaus () give . Also, the creation of 2D systems has already been achieved for bosons Hadzibabic (). The observation of the effects discussed in this letter is therefore within experimental reach once further cooling has been achieved.

In conclusion, we studied the quantum phases of a dipolar Fermi gas in 2D aligned by an external field. We demonstrated that by partially orienting the dipoles into the 2D plane, they can experience a strong -wave pairing attraction without the system collapsing. This makes the system a promising candidate to study quantum phase transitions and pairing with unconventional symmetry. We also analyzed the BKT transition to the normal state at non-zero .

Stimulating discussions with S. Stringari are acknowledged.

### References

- I. Bloch, Nature Phys. 1, 23 (2005).
- See G. Pupillo et al., e-print: arXiv:0805.1896; C. Menotti et al., e-print: arXiv:0711.3422.
- T. Lahaye et al., Nature 448, 672 (2007).
- T. Koch et al., Nature Phys. 4, 218 (2008).
- S. Ospelkaus et al., Nature Phys. 4, 622 (2008); K. K. Ni et al., e-print: arXiv:0808:2963.
- S. Kotochigova, P. S. Julienne, and E. Tiesinga, Phys. Rev. A 68, 022501 (2003).
- C. Weber et al., e-print: arXiv:0808:4077.
- H. P. Büchler et al., Phys. Rev. Lett. 98, 060404 (2007); G. E. Astrakharchik et al., ibid. 98, 060405 (2007). These calculations were for bosons but quantum statistics are not relevant in the crystalline phase.
- T. Miyakawa, T. Sogo, and H. Pu, Phys. Rev. A 77, 061603(R) (2008).
- M. A. Baranov, L . Dobrek and M. Lewenstein, Phys. Rev. Lett. 92, 250403 (2004).
- V. L. Berezinskii, JETP 34, 610 (1972); J. M. Kosterlitz and D. J. Thouless, J. Phys. C 7, 1047 (1974).
- P. M. Chaikin and T. C. Lubensky, Principles of Condensed Matter Physics (Cambridge University Press, Cambridge, 1995).
- S. Yi and H. Pu, Phys. Rev. A 73, 061602(R) (2006).
- Corresponding to half the value for the normal fluid density in a two-component superfluid; see e.g. A. J. Leggett, Rev. Mod. Phys. 47, 331 (1975).
- D. S. Petrov, M. A. Baranov, and G. V. Shlyapnikov, Phys. Rev. A 67, 031601(R) (2003).
- Z. Hadzibabic et al., Nature 441, 1118 (2006); P. Krüger et al. Phys. Rev. Lett. 99, 040402 (2007).