# Speed of the Goldstone Sound Mode of an Atomic Fermi Gas Loaded on a Square Optical Lattice with a Non-Abelian Gauge Field in the Presence of a Zeeman Field.

1. IntroductionIt is well known that topological superfluids are new states of matter that can be observed in two-dimensional atomic Fermi gases with strong Rashba spin-orbit coupling (SOC), conventional s-wave pairing, and out-of-plane magnetic field (Zeeman field) which breaks time-reversal symmetry. In what follows, we study the speed of sound of a two-component pseudospin-[up arrow], [down arrow] fermionic gas loaded on a square optical lattice with a non-Abelian gauge field A = ([alpha][[sigma].sub.y], - [alpha][[sigma].sub.x]) in the presence of a Zeeman field, where a is independently tunable parameter, and [[sigma].sub.i], i = x, y, z are the Pauli matrices. The external non-Abelian gauge field effectively creates a SOC in the gas. The pseudospin of atoms can couple with the Zeeman field and the orbital degrees of freedom of atoms, which gives rise to a topological quantum phase transition (TQPT) between gapped (topological trivial) and gapless (topological nontrivial) superfluid states. In other words, the spectrum of the Hamiltonian becomes gapless at some set of points in the momentum space. Due to absence of symmetry breaking across the TQPTs, it remains experimentally challenging to detect these phase transitions. It is experimentally possible to detect TQPTs by measuring the dynamic response of the bulk under an external force [1], by observing the edge modes [2, 3], or the momentum distributions of atoms across the phase transition boundary using a noise-correlation imaging [4] and momentum-resolved spectroscopy [5, 6].

Recent numerical calculations of the sound speed of atomic Fermi gases with s-wave attraction, synthetic Rashba SOC, and out-of-plane Zeeman field in two- and three-dimensional free space [7-9] show sharp changes of the slope of the Goldstone sound mode across the phase transition points, and this can be used to detect TQPTs. Our goal is to study whether a similar behavior of the sound speed occurs in the lattice case if the SOC is created by the non-Abelian gauge field. We found that the sharp change in the slope does not exist in the presence of lattice geometry.

It is widely accepted among the cold-atom community that the system of a two-component Fermi gas with s-wave attraction in square optical lattice with an external non-Abelian gauge field (or a synthetic Rashba SOC) and out-of-plane Zeeman field can be described by a negative-U Hubbard model [10-26]. We restrict the discussion to the case of atoms confined to the lowest-energy band (single-band model). We assume that there are M = [N.sub.[up arrow]] + [N.sub.[down arrow]] atoms distributed along N sites, and the filling factor f = [f.sub.[up arrow]] + [N.sub.[down arrow]] = M/N is smaller than unity. The corresponding polarization is P = ([f.sub.[up arrow]] - [N.sub.[down arrow]])/f. The Hamiltonian for an uniform system is [mathematical expression not reproducible], where

[mathematical expression not reproducible]. (1)

The out-of-plane Zeeman field is described by the term Hz:

[mathematical expression not reproducible]. (2)

The Hubbard part of the Hamiltonian is

[mathematical expression not reproducible]. (3)

Here, [mu] is the chemical potential, and [[??].sub.i,[sigma]] = [[psi].sup.[dagger].sub.i,[sigma]] is the density operator on site i. The Fermi operator [[psi].sup.[dagger].sub.i,[sigma]] = (1/N)[[summation].sub.k] exp(-1k.[r.sub.i])[[psi].sup.[dagger].sub.k,[sigma]] ([[psi].sub.i,[sigma]] = (1/N) [[summation].sub.k] exp(-k.[r.sub.i])[[psi].sub.[dagger].sub.k,[sigma]]) creates (destroys) a fermion on the lattice site [r.sub.i] with pseudospin projection [sigma]. The sum means sum [[summation].sub.<ij>] over the nearest-neighbor sites of the 2D lattice. [t.sup.[sigma][sigma]'.sub.i,j] = [t.sub.i,j] exp [- l[[theta].sup.[sigma][sigma]'.sub.i,j] describes the hoping of the atoms between site j and site i, and [mathematical expression not reproducible] A.dr is accumulated via the Peierls substitution phase factor. The strength of the on-site Hubbard interaction is U > 0, which corresponds to attractive interaction.

Hamiltonian (1) in a momentum representation assumes the form [mathematical expression not reproducible]. Here, [[psi].sup.[dagger].sub.k] = ([[psi].sup.[dagger].sub.k,[up arrow]], [[psi].sup.[dagger].sub.k,[down arrow]]), [xi](k) = [member of](k) - [mu], [member of](k) is the tight-binding energy, [[sigma].sub.0] is the identity matrix, and [J.sub.x,y](k) characterize the SOC.

In the case of a square lattice the vectors k = ([k.sub.x], [k.sub.y]) (in units [a.sup.-1]) are restricted to be within the Brillouin zone, i.e., a square defined by the following four points ([+ or -][pi]/a, [+ or -][pi]/a). The tight-binding energy and [J.sub.x,y] are given by

[mathematical expression not reproducible], (4)

where is the strength of the nearest-neighbor hopping (in what follows, the lattice spacing and spacing are set to unity).

The opening or closing of the gap can be achieved by varying the tunable parameter a along with the Zeeman field. We first assume that the opening or closing of the gap is achieved by varying the parameter, while by changing the out-of-plane Zeeman field the polarization is kept fixed. We also discuss the second possibility, when one can close or open the gap by changing the Zeeman field while keeping the parameter [alpha] fixed. In this case the polarization changes during the transition from gapless to gapped superfluid states.

The single-particle and collective excitations of the system under consideration manifest themselves as poles of the single-particle and two-particle Green's functions, but the corresponding expressions for the Green's functions cannot be evaluated exactly because the interaction part of the Hubbard Hamiltonian is quartic in the fermion fields. The simplest way to solve this problem is to apply the so-called mean field decoupling of the quartic interaction. According to the standard mean field theory on the Hubbard Hamiltonian, the on-site interaction is decoupled (up to Hartree-Fock correction terms) as

[mathematical expression not reproducible], (5)

where the order parameter [DELTA] = -U<[[psi].sub.i,[down arrow]] [[psi].sub.i,[up arrow]]> is a real number.

Instead of applying the mean field decoupling, we shall transform the quartic terms to a quadratic form by introducing a four-component boson field which mediates the interaction of fermions in the same manner as in quantum electrodynamics, where the photons mediate the interaction of electric charges. Green's functions are thermodynamic averages of the [[??].sub.u]-ordered products of field operators. The standard procedure for calculating Green's functions is to apply Wick's theorem. This enables us to evaluate the [[??].sub.u]-ordered products of field operators as perturbation expansions involving only wholly contracted field operators. These expansions can be summed formally to yield different equations of Green's functions. The main disadvantage of this procedure is that the validity of the equations must be verified diagram by diagram. For this reason, we shall use the method of Legendre transforms [27] of the generating functional for connected Green's functions to derive the Schwinger-Dyson [28, 29] (SD) equation [G.sup.-1] = [G.sup.(0)-1] - [summation] for the poles of the single-particle Green's function, G, and the Bethe-Salpeter [30, 31] (BS) equation [[K.sup.(0)-1] - I][PSI] = 0 for the poles of the poles of the two-particle Green's function, K. Here, [G.sup.(0)] is the free single-particle propagator, [summation] is the fermion self-energy, I is the BS kernel, and the two-particle free propagator [K.sup.(0)] = GG is a product of two fully dressed single-particle Green's functions. The kernel of the BS equation is defined as a sum of the direct interaction, [I.sub.d] = [delta][Z.sup.F]/[delta]G, and the exchange interaction [I.sub.exc] = [delta][[summation].sup.H]/[delta]G, where [[summation].sup.F] and [[summation].sup.H] are the Fock and the Hartree parts of the fermion self-energy [summation]. Since the fermion self-energy depends on the two-particle Green's function, the positions of both poles must to be obtained by solving the SD and BS equations self-consistently. Instead of solving self-consistently the SD and BS equations, we shall employ the generalized random phase approximation (GRPA). In this approximation, the single-particle excitations are obtained in the mean field approximation, while the collective modes are obtained by solving the BS equation in which single-particle Green's functions are calculated in Hartree-Fock approximation, and the BS kernel is obtained by summing ladder and bubble diagrams.

2. Functional-Integral Formalism

2.1. Field-Theoretical Approach to the Hubbard Model. The functional-integral formulation of the Hubbard model requires the representation of the Hubbard interaction [mathematical expression not reproducible] in terms of squares of fermion operators. This can be done by employing a certain Hubbard-Stratonovich transformation. This field-theoretical approach to the Hubbard Hamiltonian has already been used to describe the collective modes of ultracold [sup.6][Li-.sup.40]K mixture in a square optical lattice (see [32]). Hamiltonian (1) in the present study has two more interactions. It is possible to include both the SOC and the Zeeman field into initial definition of the free fermion Green's function, and therefore most of the general equations in [32] remain unchanged. In particular, the general expressions for the SD and BS equations are the same as those without the SOC and the Zeeman field, but with more complicated fermion self-energy and BS kernel.

As it is well known, Green's functions in the functional-integral approach are defined by means of the so-called generating functional with sources for the boson and fermion fields, but the corresponding functional integrals cannot be evaluated exactly because the interaction part of the Hubbard Hamiltonian is quartic in the Grassmann fermion fields. However, we can transform the quartic terms to a quadratic form by introducing a model system which consists of a four-component boson field [A.sub.[alpha]](z) interacting with fermion fields [mathematical expression not reproducible] and [mathematical expression not reproducible], where

[mathematical expression not reproducible]. (6)

Here [alpha] = 1, 2, 3, 4, x = ([r.sub.i], u), y = ([r.sub.j], u'), and z = ([r.sub.k], v) are composite variables, where [r.sub.i], [r.sub.j], and [r.sub.k] are the lattice site vectors. According to imaginary-time (Matsubara) formalism the variables 0 [less than or equal to] v, u, u' [less than or equal to] [beta] = [([k.sub.B]T).sup.-1]; [k.sub.B] is the Boltzmann constant ([??] = [k.sub.B] = 1).

The field operators (6) allow us to define generalized single-particle Green's function by using a tensor product of these two matrices. The corresponding Green function, represented by a 4x4 matrix, includes all possible thermodynamic averages:

[mathematical expression not reproducible]. (7)

The action of the above-mentioned model system is assumed to be of the following form S = [S.sup.(F).sub.0] + [S.sup.(B).sub.0] + [S.sup.(F-B)], where [mathematical expression not reproducible]. The fermion action [s.sup.(F).sub.0] corresponds to the Hamiltonian [[??].sub.0]. The corresponding inverse Green function of free fermions [[??].sup.(0)-1] (y; x) is given by the following 4x4 matrix:

[mathematical expression not reproducible]. (8)

The symbol [mathematical expression not reproducible] is used to denote [[beta].sup.-1] [[summation].sub.m] (for fermion fields [[omega].sub.m] = (2[pi]/[beta])(m + 1/2); m = 0, [+ or -]1, [+ or -]2, ...).

The Fourier transform of noninteracting Green's function is

[mathematical expression not reproducible]. (9)

After including the Zeeman field and the SOC interaction into free fermion Green's function (9), one can perform the same steps as in [32] to come up with the following single-particle Green function:

[mathematical expression not reproducible], (10)

where [mathematical expression not reproducible], where [DELTA] is a constant in space. Physically, it describes a superfluid state of Cooper pairs with zero momentum when the pairing is only between atoms with equal and opposite momenta.

2.2. Mean Field Approximation. The Fourier transform of zero-temperature single-particle Green's function (10) in the mean field approximation is

[mathematical expression not reproducible]. (11)

The matrix elements [mathematical expression not reproducible] of [[??].sub.MF] (k, w) are given by

[mathematical expression not reproducible], (12)

where the corresponding expressions for [mathematical expression not reproducible], and [mathematical expression not reproducible] can be obtained by inverting matrix (11) and setting T [right arrow] 0.

For a fixed filing factor f = [f.sub.[up arrow]] + [f.sub.[down arrow]] and at a given Zeeman field h, the chemical potential [mu] and the gap [DELTA] have to be determined by solving the mean field number and gap equations:

[mathematical expression not reproducible], (13)

[mathematical expression not reproducible]. (14)

Here, S(k) = [J.sup.2.sub.x](k) + [J.sup.2.sub.y](k), F([omega]) is the Fermi- Dirac distribution function, and the following notations have been introduced:

[mathematical expression not reproducible], (15)

[mathematical expression not reproducible]. (16)

The number equation (13) follows from the definitions of the filing factors [mathematical expression not reproducible] through corresponding mean field Green's functions. It is known that the gap equation can be obtained by minimizing the mean field Helmholtz free energy with respect to [DELTA].

As we mentioned before, in order to keep the polarization fixed, one has to adjust the corresponding Zeeman field. This means that [mu], [DELTA], and h have to be calculated by solving the number and the gap equations, along with the following equation:

[mathematical expression not reproducible]. (17)

Another interesting feature is that because of the SOC term in Hamiltonian (1) the pairing field contains both singlet and triplet component. The singlet [[PHI].sub.[down arrow][up arrow]](k) = - [[PHI].sub.[up arrow][down arrow]] = <[[psi].sub.k[down arrow]] [[psi].sub.k[up arrow]]) and triplet [[PHI].sub.[up arrow][up arrow]](k) = <[[psi].sub.k[up arrow]] [[psi].sub.k[down arrow]]> amplitudes, obtained by means of Green's function elements [G.sup.MF.sub.23], [G.sup.MF.sub.13], and [G.sup.MF.sub.24], are

[mathematical expression not reproducible]. (18)

With the help of the pairing amplitudes, one can calculate the total condensate fraction [f.sub.c] = [f.sub.s] + [f.sub.tr], where [f.sub.s] and [f.sub.tr] are the singlet and the triplet contributions, correspondingly. At zero temperature, we obtain

[mathematical expression not reproducible], (19)

[mathematical expression not reproducible]. (20)

2.3. The Bethe-Salpeter Equation for the Collective Excitations in the Generalized Random Phase Approximation. The spectrum of the collective modes will be obtained by solving the BS equation in the GRPA. As we have already mentioned, the kernel of the BS equation is a sum of the direct [I.sub.d] = [delta][[summation].sup.F]/[delta]G and exchange [I.sub.exc] = [partial derivative][[summation].sub.H]/[delta]G interactions, written as derivatives of the Fock and the Hartree parts of the self-energy. Thus, in the GRPA the corresponding equation for the BS amplitude [mathematical expression not reproducible] is given by (see [32])

[mathematical expression not reproducible], (21)

where [mathematical expression not reproducible] and [mathematical expression not reproducible] are the direct and the exchange interactions, correspondingly. The two-particle propagator [K.sup.(0)] in the GRPA is defined as follows:

[mathematical expression not reproducible]. (22)

The BS equation (21) can be written in the matrix form as [mathematical expression not reproducible], where I is the unit matrix, the matrix [??] is a 16x16 matrix, and the transposed matrix of [??] is given by

[mathematical expression not reproducible]. (23)

At a given vector Q from the Brillouin zone, the collective excitation spectrum [omega](Q) is obtained by the condition [mathematical expression not reproducible] to have a nontrivial solution of the BS equation. By applying simple matrix algebra, the above secular 16 x 16 determinant can be simplified to a 10 x 10 BS determinant det [absolute value of BS([omega](Q), Q)] =0, where the BS matrix is given by

[mathematical expression not reproducible]. (24)

The elements [A.sub.i,j], [B.sub.i,j], [C.sub.i,j], and [D.sub.i,j] of the blocks [A.sub.2x2] ([omega], Q), [B.sub.2x8]([omega], Q), [C.sub.8x2]([omega], Q), and [D.sub.8x8]([omega], Q) are given in the Appendix in terms of the propagators [mathematical expression not reproducible] ([omega], Q).

In contrast to our functional-integral formalism, one can use the Hubbard-Stratonovich transformation to introduce the energy gap as an order parameter field. In this approach, one can integrate out the fermion fields and to arrive at an effective action. The next steps are to consider the state which corresponds to the saddle point of the effective action and to write the effective action as a series in powers of the fluctuations and their derivatives. The exact result can be obtained by explicitly calculating the terms up to second order in the fluctuations and their derivatives. This approximation is called the Gaussian approximation. Within this approximation, the collective excitation spectrum [omega](Q) is defined by the 2x2 Gaussian secular determinant

det |[A.sub.2x2] ([omega](Q), Q)] = 0, (25)

while within the BS formalism the secular determinant is more complicated (([omega](Q), Q) dependence is understood):

det|[A.sub.2x2] - [B.sub.2x8] x [D.sup.-1.sub.8x8] x [C.sub.8x2]] = 0 (26)

3. Numerical Results

There are various mean field quantities of physical interest, such as the chemical potential, the pairing gap, the singlet and triplet pairing amplitudes, and the singlet and triplet condensate fractions. We focus on the zero-temperature case assuming two different filling factors of f = 0.4 and f = 0.6. The strength of the attractive interaction is U/t = 4, and for both filling factors the polarization is fixed to be P = 0.48.

In Figure 1, we plot the chemical potential [mu], the gap [DELTA], and the corresponding Zeeman field h, as functions of the SOC parameter [alpha], obtained by solving the mean field equations (13)-(17). In Figure 2 we have shown the singlet and the triplet condensate fractions, as functions of the SOC parameter [alpha], obtained by solving the mean field equations along with definitions (19) and (20). It can be seen that at the value [alpha] = [pi]/2, where [xi](k) = -[mu] and S(k) = 4[t.sup.2] [[sin.sup.2]([k.sub.x]) + [sin.sup.2]([k.sub.y])], the chemical potential has a minimum, while the gap and the Zeeman field reach their maximum values. As expected, the singlet and triplet condensate fraction-graphs have maximum at the same [alpha] as the gap.

Next, we have calculated the speed of sound as a function of the SOC parameter [alpha]. The sound speed c is proportional to the slope u of the Goldstone sound mode c = uat/[??]. The slope u has been numerically obtained using the collective-mode dispersion [omega](Q) for small wave vectors in [Q.sub.x]-direction.

It is clear that the single-particle dispersion, given by (16), becomes precisely zero at some special points [k.sub.0] in momentum space satisfying the condition S([k.sub.0]) = 0 when h = [square root of [([member of]([k.sub.0]) - [mu]).sup.2] + [DELTA]2]. The condition S([k.sub.0]) = 0 is met at the following points: (0, 0), (0, [+ or -][pi]), ([+ or -][pi], 0), and ([+ or -][pi], [+ or -][pi]). For U/t = 4, P = 0.48, and a filling factor f = 0.4, we have two points at which the single-particle dispersion [omega]([k.sub.0,1]) = 0. The first one is at the center of the Brillouin zone [k.sub.0] = (0, 0), and [omega]([k.sub.0]) = 0 for [alpha]/[pi] between 0.39 and 0.40. At the second point, [k.sub.1] = ([pi], [pi]), we found also the existence of opening and closing of a gap in momentum space when 0.61 < [alpha]/[pi] < 0.62. Therefore, at both points one should have a TQPT between superfluid states with different topologies. During the TQPT the different phases of matter are not characterized by an order parameter, but rather an integer number, the Chern number, describing the system as a whole. This is because the different phases across the TQPT are characterized by different topologies of Fermi surface instead of being classified by different symmetries.

In Figure 3, we have shown the slope calculated by using the Gaussian and the BS secular determinants (25) and (26), correspondingly. It is worth mentioning the significant difference between the slope obtained within the Gaussian approximation and the slope obtained by the BS secular determinant. The arrows show the points where the opening and closing of a gap take place. As can be seen, the sharp change in the slope does not exist in the presence of a non-Abelian gauge field.

The boundaries for the TQPT in [7-9] were accessed by varying the Zeeman filed at a constant strength of the synthetic Rashba SOC. It was found that on a topological trivial side, where h < [h.sub.c] = [square root of [([member of]([k.sub.0]) - [mu]).sup.2] + [DELTA]2], the slope is suppressed by the Zeeman field, while in the nontrivial side, where h > [h.sub.c] = [square root of [([member of]([k.sub.0]) - [mu]).sup.2] + [DELTA]2], the slope is enhanced by the Zeeman field. As mentioned in [7, 8], this is quite unusual since we generally expect that the superfluidity should be suppressed by the Zeeman field.

We have fixed the non-Abelian parameter [alpha] = 0.38[pi], and we have accessed the opening and the closing of a gap in momentum space at the point [k.sub.0] = (0, 0) by varying the Zeeman filed. To do this, first we used (13)-(14) to calculate numerically the chemical potential and the gap as functions of the Zeeman field. The results are shown in Figure 4. The TQPT takes place for [h.sub.c]/t = 1.39 and [h.sub.c]/t = 1.49. The insert shows the behavior of the slope close to [h.sub.c]/t = 1.49. We found that there is no sharp change in the slope across the two TQPT points.

4. Discussion

In this paper, we are concerned with the question which naturally arises here as to whether there is a possibility of detecting TQPTs by monitoring the behavior of slope of the sound mode if a non-Abelian gauge field is used instead of a synthetic Rashba SOC. The answer is not trivial because the non-Abelian gauge field is taken into account via the Peierls substitution. As a result, the tight-binding energy [member of](k) = -2t cos([alpha])[cos([k.sub.x]) + cos([k.sub.y])] does depend on the parameter [alpha]. The second difference is that the nondiagonal part of [[theta].sup.[sigma][sigma]'.sub.i,j], i.e., [sigma] [not equal to] [sigma]', describes the spin-flipped hopping between site i and site j. In the case of a synthetic Rashba SOC we do not have spin-flipped hopping, and the hopping to the nearest-neighbor site does not depend on the strength of Rashba SOC.

In summary, we have derived the BS equation in the GRPA for the collective excitation energy of a Fermi gas loaded on a square optical lattice with a non-Abelian gauge field in the presence of a Zeeman field. To the best of our knowledge, there is no other calculation of the speed of sound in a lattice geometry as a function of the non-Abelian parameter with which we can make comparisons.

According to our numerical calculations, there is no sharp change of the slope of the Goldstone sound mode across the phase transition point. It is found that the Gaussian approximation significantly overestimates the speed of sound of the Goldstone mode compared to the value obtained within the BS formalism.

Appendix

Blocks of the BS Matrix

The two-particle propagator in BS equation (21) is defined as

[mathematical expression not reproducible]. (A.1)

The elements [A.sub.i,j], [B.sub.i,j], [C.sub.i,j] and [D.sub.i,j] of the blocks [A.sub.2x2]([omega], Q), [B.sub.2x8] ([omega], Q), [C.sub.8x2] ([omega], Q), and [D.sub.8x8] ([omega], Q) are given in terms of the propagators [mathematical expression not reproducible] as follows:

[mathematical expression not reproducible], (A.2)

[mathematical expression not reproducible], (A.3)

[mathematical expression not reproducible], (A.4)

[mathematical expression not reproducible]. (A.5)

https://doi.org/10.1155/2018/6809265

Conflicts of Interest

The authors declare that they have no conflicts of interest. Acknowledgments

Israel Chavez Villalpando gratefully acknowledges funding from CONACyT Grant no. 291001 and UNAM-DGAPA PAPIIT IN102417.

References

[1] H. Pino, E. Alba, J. Taron, J. J. Garcia-Ripoll, and N. Barberan, "Hall response of interacting bosonic atoms in strong gauge fields: from condensed to fractional-quantum-Hall states," Physical Review A: Atomic, Molecular and Optical Physics, vol. 87, no. 5, Article ID 053611, 2013.

[2] N. Goldman, J. Beugnon, and F. Gerbier, "Detecting chiral edge states in the hofstadter optical lattice," Physical Review Letters, vol. 108, no. 25, Article ID 255303, 2012.

[3] N. Goldman, J. Beugnon, and F. Gerbier, "Identifying topological edge states in 2D optical lattices using light scattering," The European Physical Journal Special Topics, vol. 217, no. 1, pp. 135-152, 2013.

[4] M. Greiner, C. A. Regal, J. T. Stewart, and D. S. Jin, "Probing pair-correlated fermionic atoms through correlations in atom shot noise," Physical Review Letters, vol. 94, no. 11, Article ID 110401, 2005.

[5] L. Tarruell, D. Greif, T. Uehlinger, G. Jotzu, and T. Esslinger, "Creating, moving and merging Dirac points with a Fermi gas in a tunable honeycomb lattice," Nature, vol. 483, no. 7389, pp. 302-305, 2012.

[6] Z. Fu, L. Huang, Z. Meng et al., "Radio-frequency spectroscopy of a strongly interacting spin-orbit-coupled Fermi gas," Physical Review A: Atomic, Molecular and Optical Physics, vol. 87, no. 5, Article ID 053619, 2013.

[7] L. He and X.-G. Huang, "Nonperturbative effects on the ferromagnetic transition in repulsive Fermi gases," Physical Review A: Atomic, Molecular and Optical Physics, vol. 85, no. 4, Article ID 043618, 2012.

[8] L. He and X.-G. Huang, "Superfluidity and collective modes in Rashba spin-orbit coupled Fermi gases," Annals of Physics, vol. 337, pp. 163-207, 2013.

[9] Y. Xu, R.-L. Chu, and C. Zhang, "Anisotropic weyl fermions from the quasiparticle excitation spectrum of a 3D fulde-ferrell superfluid," Physical Review Letters, vol. 112, no. 13, Article ID 136402, 2014.

[10] 1.1. Satija, D. C. Dakin, J. Y. Vaishnav, and C. W. Clark, "Physics of a two-dimensional electron gas with cold atoms in non-Abelian gauge potentials," Physical Review A: Atomic, Molecular and Optical Physics, vol. 77, no. 4, Article ID 043410, 2008.

[11] A. Bermudez, N. Goldman, A. Kubashiak, M. Lewenstein, and M. A. Martin-Delgado, "Non-abelian optical lattices: anomalous quantum hall effect and dirac fermions," Physical Review Letters, vol. 103, Article ID 035301, 2009.

[12] A. Bermudez, N. Goldman, A. Kubasiak, M. Lewenstein, and M. A. Martin-Delgado, "Topological phase transitions in the non-Abelian honeycomb lattice," New Journal of Physics, vol. 12, Article ID 033041, 2010.

[13] N. Goldman, A. Kubasiak, P. Gaspard, and M. Lewenstein, "Ultracold atomic gases in non-Abelian gauge potentials: the case of constant Wilson loop," Physical Review A: Atomic, Molecular and Optical Physics, vol 79, no. 2, Article ID 023624, 2009.

[14] A. Kubasiak, P. Massignan, and M. Lewenstein, "Topological superfluids on a lattice with non-Abelian gauge fields," EPL (Europhysics Letters), vol 92, no. 4, Article ID 46004, 2010.

[15] L. Wang and L. Fu, "Interaction-induced merging of Dirac points in non-Abelian optical lattices," Physical Review A: Atomic, Molecular and Optical Physics, vol 87, no. 5, Article ID 053612, 2013.

[16] M. Iskin, "Superfluid phases of ultracold Fermi gases on a checkerboard superlattice," Physical Review A: Atomic, Molecular and Optical Physics, vol 88, no. 5, Article ID 013631, 2013.

[17] M. Burrello, I. C. Fulga, E. Alba, L. Lepori, and A. Trombettoni, "Topological phase transitions driven by non-Abelian gauge potentials in optical square lattices," Physical Review A: Atomic, Molecular and Optical Physics, vol. 88, no. 5, Article ID 053619, 2013.

[18] Q. Sun, G.-B. Zhu, W.-M. Liu, and A.-C. Ji, "Spin-orbit coupling effects on the superfluidity of a Fermi gas in an optical lattice," Physical Review A: Atomic, Molecular and Optical Physics, vol. 88, no. 6, Article ID 063637, 2013.

[19] Y. Zhang and C. Zhang, "Bose-einstein condensates in spin-orbit-coupled optical lattices: flat bands and superfluidity," Physical Review A, vol. 87, Article ID 023611, 2013.

[20] F. Lin, C. Zhang, and V. W. Scarola, "Emergent kinetics and fractionalized charge in 1D spin-orbit coupled flatband optical lattices," Physical Review Letters, vol. 112, no. 11, Article ID 110404, 2014.

[21] D. Toniolo and J. Linder, "Superfluidity breakdown and multiple roton gaps in spin-orbit-coupled Bose-Einstein condensates in an optical lattice," Physical Review A: Atomic, Molecular and Optical Physics, vol. 89, no. 6, Article ID 061605(R), 2014.

[22] C. Qu, M. Gong, and C. Zhang, "Fulde-Ferrell-Larkin-Ovchinnikov or Majorana superfluids: the fate of fermionic cold atoms in spin-orbit-coupled optical lattices," Physical Review A, vol. 89, Article ID 053618, 2014.

[23] Y. Xu, C. Qu, M. Gong, and C. Zhang, "Competing superfluid orders in spin-orbit-coupled fermionic cold-atom optical lattices," Physical Review A, vol. 89, Article ID 013607, 2014.

[24] H. Zhai, "Degenerate quantum gases with spin-orbit coupling: a review," Reports on Progress in Physics, vol. 78, no. 2, Article ID 026001, 25 pages, 2015.

[25] H.-K. Tang, X. Yang, J. Sun, and H.-Q. Lin, "Berezinskii-Kosterlitz-Thoules phase transition of spin-orbit coupled Fermi gas in optical lattice," Europhysics Letters, vol. 107, Article ID 40003, 2014.

[26] M. Iskin, "Topological superfluids on a square optical lattice with non-Abelian gauge fields: effects of next-nearest-neighbor hopping in the BCS-BEC evolution," Physical Review A: Atomic, Molecular and Optical Physics, vol. 93, Article ID 013608, 2016.

[27] C. Itzykson and J. B. Zuber, Quantum Field Theory, McGraw-Hill, New York, NY, USA, 1980.

[28] J. Schwinger, "The theory of quantizied fields," Physical Review A: Atomic, Molecular and Optical Physics, vol. 82, p. 914, 1951.

[29] F. J. Dyson, "The S matrix in quantum electrodynamics," Physical Review A: Atomic, Molecular and Optical Physics, vol. 75, no. 11, pp. 1736-1755, 1949.

[30] H. A. Bethe and E. E. Salpeter, "A relativistic equation for bound-state problems," Physical Review A: Atomic, Molecular and Optical Physics, vol. 82, p. 309, 1951.

[31] H. A. Bethe and E. E. Salpeter, "A relativistic equation for bound-state problems," Physical Review A: Atomic, Molecular and Optical Physics, vol. 84, p. 1232, 1951.

[32] Z. G. Koinov, "Collective modes of an ultracold [sup.6][Li-.sup.40] K mixture in an optical lattice," Advances in Condensed Matter Physics, vol. 2015, Article ID 952852, 12 pages, 2015.

Zlatko Koinov (iD) (1) and Israel Chavez Villalpando (2)

(1) Department of Physics and Astronomy, University of Texas at San Antonio, San Antonio, TX 78249, USA

(2) Instituto de Investigaciones en Materiales, Universidad Nacional Autonoma de Mexico, Apartado Postal 20-364, 01000 Mexico Distrito Federal, Mexico

Correspondence should be addressed to Zlatko Koinov; zlatko.koinov@utsa.edu

Received 15 December 2017; Revised 24 May 2018; Accepted 27 May 2018; Published 2 July 2018

Academic Editor: Oleg Derzhko

Caption: Figure 1: Chemical potential pairing gap A, and the Zeeman field h of a Fermi gas in a square optical lattice subject to a non-Abelian gauge field A = ([alpha][[sigma].sub.y], -[alpha][[sigma].sub.x]). The graphs are obtained by solving the mean field equations (13)-(17) using two different filling factors: f = 0.4 and f = 0.6. The on-site attractive strength is U = 4t, and the polarization is P = 0.48.

Caption: Figure 2: The singlet [f.sub.s] and the triplet [f.sub.tr] contributions to the total condensed fraction f = [f.sub.s] + [f.sub.tr] of a Fermi gas in a square optical lattice subject to a non-Abelian gauge field A = ([alpha][[sigma].sub.y], - [alpha][[sigma].sub.x]). The graphs are obtained by solving the mean field equations along with definitions (19) and (20). The system parameters are the same as in Figure 1.

Caption: Figure 3: The slope of the Goldstone sound mode of a Fermi gas in a square optical lattice subject to a non-Abelian gauge field A = ([alpha][[sigma].sub.y], -[alpha][[sigma].sub.x]). The graphs are obtained by applying the Gaussian and the BS secular determinants. The system parameters are the same as in Figure 1. The arrows show the values of a for which the opening and closing of a gap take place.

Caption: Figure 4: Chemical potential [mu], pairing gap [DELTA], and the slope of the Goldstone sound mode of a Fermi gas in a square optical lattice subject to a non-Abelian gauge field with [alpha] = 0.38[pi] as a function of the Zeeman field h. The chemical potential and the gap are obtained by solving the mean field number and gap equations (13)-(14). The arrows show the values of the Zeeman field for which the opening and closing of a gap take place. The insert shows the slope close to [h.sub.c]/t = 1.49.

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Koinov, Zlatko; Villalpando, Israel Chavez |

Publication: | Advances in Condensed Matter Physics |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Jan 1, 2018 |

Words: | 5475 |

Previous Article: | Cytotoxicity Evaluation of Ammonia-Modified Graphene Oxide Particles in Lung Cancer Cells and Embryonic Stem Cells. |

Next Article: | Highly Efficient Broadband Light Absorber Based on Nonuniform Hyperbolic Metamaterial Film. |

Topics: |