Two-component Fermi System Research Articles

Attraction from kinetic frustration in ladder systems

We analyze the formation of multiparticle bound states in ladders with frustrated kinetic energy in two-component bosonic and two-component fermionic systems. We focus on the regime of light doping relative to insulating states at half-filling, spin polarization close to 100%, and strong repulsive interactions. A special feature of these systems is that the binding energy scales with single-particle tunneling t rather than exchange interactions, since effective attraction arises from alleviating kinetic frustration. For two-component Fermi systems on a zigzag ladder we find a bound state between a hole and a flipped spin (magnon) with a binding energy that can be as large as 0.6t. We demonstrate that magnon-hole attraction leads to formation of clusters comprising several holes and magnons, and we expound on antiferromagentic correlations for the transverse spin components inside the clusters. We identify several many-body states that result from self-organization of multiparticle bound states, including a Luttinger liquid of hole-magnon pairs and a density wave state of two-hole–three-magnon composites. We establish a symmetry between the spectra of Bose and Fermi systems and use it to establish the existence of antibound states in two-component Bose mixtures with SU(2) symmetric repulsion on a zigzag ladder. We also consider Bose and Fermi systems on a square ladder with flux and demonstrate that both systems support bound states. We discuss experimental signatures of multiparticle bound states in both equilibrium and dynamical experiments. We point out intriguing connections between these systems and the quark bag model in QCD. Published by the American Physical Society 2024

Spatially random disorder in unitary fermion system in (4 − ϵ)-dimensions and effective action at finite temperature

Non-relativistic conformal field theory is significant to understand various aspects of an ultra-cold system. In this paper, we study a non-relativistic system of two-component fermions interacting with a complex boson with Yukawa-like interactions near d = 4-spatial dimensions in the presence of a quenched disorder. The homogeneous theory flows to an interacting fixed point describing a unitary fermion system. In the presence of the disorder, we find that the system has an interesting phase structure in the space of the coupling constants and exhibits an interacting disorder fixed point in ϵ-expansion. The correlation function obeys Lifshitz scaling behaviour at the disorder fixed point with the anisotropic exponent being z = 2 + γE. We also study the disorder system at finite temperature and compute the leading contribution to the 1PI effective action.

Density-Induced Hadron–Quark Crossover via the Formation of Cooper Triples

We discuss the hadron–quark crossover accompanied by the formation of Cooper triples (three-body counterpart of Cooper pairs) by analogy with the Bose–Einstein condensate to Bardeen–Cooper–Schrieffer crossover in two-component fermionic systems. Such a crossover is different from a phase transition, which often involves symmetry breaking. We calculate the in-medium three-body energy from the three-body T-matrix with a phenomenological three-body force characterizing a bound hadronic state in vacuum. With increasing density, the hadronic bound-state pole smoothly undergoes a crossover toward the Cooper triple phase where the in-medium three-body clusters coexist with the quark Fermi sea. The relation to the quarkyonic matter model can also be found in a natural manner.

Cosmological Evolution of a Statistical System of Degenerate Scalarly Charged Fermions with an Asymmetric Scalar Doublet. II. One-Component System of Doubly Charged Fermions$${}^{a}$$

Based on the previously formulated mathematical model of a statistical system with scalar interaction of fermions, a cosmological model based on a one-component statistical system of doubly scalar charged degenerate fermions interacting with an asymmetric scalar doublet-canonical and phantom scalar fields-is studied. The connection of the presented model with previously studied models based on one-component and two-component fermion systems is investigated. The asymptotic and limiting properties of the cosmological model are investigated, it is shown that among all models there is a class of models with a finite lifetime. The asymptotic behavior of models near the corresponding singularities is investigated, a qualitative analysis of the corresponding dynamical system is carried out, and numerical implementations of such models are constructed. Based on numerical integration, it is shown that in the presented model there can be transitions from a stable asymptotically vacuum state with a zero canonical field and a constant phantom field corresponding to the phase of cosmological compression to a symmetric state corresponding to the expansion phase. The time interval of the transition between phases is accompanied by oscillations of the canonical scalar field.

Effective interactions between nuclear clusters

The effective interactions between two nuclear clusters, $d+d$, $t+t$, and $\alpha+\alpha$, are investigated within a cluster model using local nucleon-nucleon~($NN$) forces. It is shown that the interaction in the spin-aligned $d+d$ system is repulsive for all inter-cluster distances, whereas the $\alpha+\alpha$ and spin-aligned $t+t$ systems are attractive at intermediate distances. The Pauli blocking between identical-nucleon pairs is responsible for the cluster-cluster repulsion and becomes dominant in the shallow binding limit. We demonstrate that two $d$-clusters could be bound if the $NN$ force has nonzero range and is strong enough to form a deeply bound $d$-cluster, or if the $NN$ force has both even-parity and odd-parity attraction. Effective dimer-dimer interactions for general quantum systems of two-component fermions are also discussed in heavy-light mass limit, where one component is much heavier than the other, and their relation to inter-cluster interactions in nuclear systems are discussed. Our findings provide a conceptual foundation for conclusions obtained numerically in the literature, that increasing the range or strength of the local part of the attractive nucleon-nucleon interaction results in a more attractive cluster-cluster interaction.

Superfluid condensate fraction and pairing wave function of the unitary Fermi gas

The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gases and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this work we present auxiliary-field lattice Monte Carlo simulations using a novel lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test we compute the ground state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy $E_{FG}$, we find $E_0/E_{FG}= 0.369(2), 0.372(2)$, using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is $\alpha = 0.43(2)$. We also extract the pairing wave function and find the pair correlation length to be $\zeta_pk_F = 1.8(3) \hbar$, where $k_F$ is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic P-wave and D-wave superfluids.

Order-Disorder Transition in the Half-Filled Two-Component Lattice Fermion Model with Nearest-Neighbor Repulsion

In this work, we study a model of asymmetric two-component lattice fermion system at half-filling, where particles of both species (distinguishable by, e.g., spin) can interact only via nearest-neighbor repulsion W. The analysis is performed on the Bethe lattice using the Hartree-Fock-type mean-field approximation, which is rigorous in the limit of high dimensions. At sufficiently low temperatures, both antiferromagnetic and charge orders (related to inhomogeneous distribution of particles with both spins) coexist in the system. We find that an order-disorder continuous transition occurs with increasing temperature. The transition temperature depends on the ratio t↑/t↓ of the hopping amplitudes of both fermion species (i.e., the asymmetry of the model). For fixed W, it is the biggest if one component is localized (e.g., t↑ = 0), and it decreases to its minimal value for the same hopping amplitudes (t↑ = t↓). Moreover, it increases with W for fixed t↑/t↓. Dependencies of order parameters with model parameters and temperature are also presented.

Comparative study of the condensate fraction of balanced and imbalanced uniform fermionic systems in the BEC-BCS crossover

In this paper, we have made a comparative study of the condensate fraction of balanced and imbalanced two-component uniform fermionic systems in the Bose–Einstein condensation-Bardeen–Cooper–Schrieffer (BEC-BCS) crossover. We have used extended BCS equations to calculate the condensate fraction of Cooper pairs as a function of chemical potential and gap energy in a balanced system using a path integral approach. For an imbalanced system, the same equation has been derived as a function of average chemical potential, imbalanced chemical potential, gap energy, and temperature. The condensate fraction in a balanced fermionic system varies with the interaction parameter, whereas its behavior in an imbalanced system changes with the interaction parameter and polarization, in different regimes. In particular, the effect of polarization on the condensate fraction has been studied in an imbalanced system for odd values of polarization. We have compared the condensate fraction for both systems analytically and graphically and observed that for an imbalanced system, it decreases with increasing polarization and almost vanishes for higher values of polarization in the BEC regime.

Leading order corrections to the free energy and phase separation in two-component Fermion systems

We study phase separation in a dilute two-component Fermi system with attractive interactions as a function of the coupling strength and the polarization or number density asymmetry between the two components. In weak and strong couplings with a finite number density asymmetry, phase separation is energetically more favorable. A heterogeneous phase containing a symmetric superfluid component and an asymmetric normal phase has lower energy than a homogeneous normal phase. We show that for a small number density asymmetry, taking into consideration the leading order corrections at order of the interaction parameter, phase separation is stable against the normal phase in the whole BCS range. We investigate the consequences of the consideration of the leading order corrections to the thermodynamic potentials of the normal and BCS phase on the Chandrasekhar–Clogston limit. We have also investigated the stability of a Bose–Fermi mixture in the far-BEC limit. We find that the molecular BEC is locally stable against an external magnetic field h, provided is smaller than the pairing gap .

Lattice calculations for two-component fermion systems with unequal masses: one dimension

We consider systems of two-component fermions with unequal masses and interacting via a short-range attractive potential. We discuss the case where the two-component fermions form a shallow dimer with large scattering length. The three-fermion and four-fermion systems with such properties are universal and characterized by the twofermion scattering length a ff and the ratio of the mass of spin-↑ fermion to the mass of spin-↓ fermion, m ↑ /m ↓ . In this study, using lattice effective field theory, we analyze fermion-dimer and dimer-dimer systems and calculate the universal fermion-dimer and dimer-dimer scattering lengths for various values of the mass ratio m ↑ /m ↓ . We find that these universal scattering lengths increase logarithmically with the mass ratio m ↑ /m ↓ .

Fulde-Ferrell-Larkin-Ovchinnikov pairing states of a polarized dipolar Fermi gas trapped in a one-dimensional optical lattice

We study the interplay between the long- and short-range interaction of a one-dimensional optical lattice system of two-component dipolar fermions by using the density matrix renormalization group method. The atomic density profile, pairing-pairing correlation function, and the compressibility are calculated in the ground state, from which we identify the parameter region of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) pairing state, half-metal (HM) state, FFLO-HM state, and the normal polarized state, and thus the phase diagram in the coordinates of the long- and short-range interaction strength. The effect of the long-range dipolar interaction on the FFLO state is discussed in details. We find that the long-range part of the dipole-dipole interaction does not sweep away the FFLO superconducting region that is driven by the short-range interaction in the Hubbard model, and thus the FFLO state survives in the wide parameter space of the long-range interaction, polarization and the filling.

Equation of state of non-relativistic matter from automated perturbation theory and complex Langevin

We calculate the pressure and density of polarized non-relativistic systems of two-component fermions coupled via a contact interaction at finite temperature. For the unpolarized one-dimensional system with an attractive interaction, we perform a thirdorder lattice perturbation theory calculation and assess its convergence by comparing with hybrid Monte Carlo. In that regime, we also demonstrate agreement with real Langevin. For the repulsive unpolarized one-dimensional system, where there is a so-called complex phase problem, we present lattice perturbation theory as well as complex Langevin calculations. For our studies, we employ a Hubbard-Stratonovich transformation to decouple the interaction and automate the application of Wick’s theorem for perturbative calculations, which generates the diagrammatic expansion at any order. We find excellent agreement between the results from our perturbative calculations and stochastic studies in the weakly interacting regime. In addition, we show predictions for the strong coupling regime as well as for the polarized one-dimensional system. Finally, we show a first estimate for the equation of state in three dimensions where we focus on the polarized unitary Fermi gas.

Inhomogeneity growth in two-component fermionic systems

The dynamics of fermionic many-body systems is investigated in the framework of Boltzmann-Langevin (BL) stochastic one-body approaches. Within the recently introduced BLOB model, we examine the interplay between mean-field effects and two-body correlations, of stochastic nature, for nuclear matter at moderate temperature and in several density conditions, corresponding to stable or mechanically unstable situations. Numerical results are compared to analytic expectations for the fluctuation amplitude of isoscalar and isovector densities, probing the link to the properties of the employed effective interaction, namely symmetry energy (for isovector modes) and incompressibility (for isoscalar modes). For unstable systems, clusterization is observed. The associated features are compared to analytical results for the typical length and time scales characterizing the growth of unstable modes in nuclear matter and for the isotopic variance of the emerging fragments. We show that the BLOB model is generally better suited than simplified approaches previously introduced to solve the BL equation, and it is therefore more advantageous in applications to open systems, like heavy ion collisions.

Multi-regulator functional renormalization group for many-fermion systems

We propose a method of multi-regulator functional renormalization group (MR-FRG) which is a novel formulation of functional renormalization group with multiple infrared (IR) regulators. It is applied to a two-component fermionic system with an attractive contact interaction to study crossover phenomena between the Bardeen–Cooper–Schrieffer (BCS) phase and the Bose–Einstein condensation (BEC) phase. To control both the fermionic one-particle excitations and the bosonic collective excitations, IR regulators are introduced, one for the fermionic two-point function and another for the four-fermion vertex. It is shown that the Nozières–Schmitt-Rink (NSR) theory, which is successful to capture qualitative features of the BCS–BEC crossover, can be derived from MR–FRG. Some aspects of MR-FRG to go beyond the NSR theory are also discussed.

Population imbalance in the extended Fermi-Hubbard model

We study the interplay between population imbalance in a two-component fermionic system and nearest-neighbor interaction using matrix product states method. Our analysis reveals the existence of a new type of Fulde-Ferrell-Larkin-Ovchinnikov phase in the presence of competing interactions. Furthermore, we find distinct evidence for the presence of hidden order in the system. We present an effective model to understand the emergent oscillations in the string correlations due to the imbalance, and show how they can become an efficient tool to investigate systems with imbalance.

Magnetism in Strongly Interacting One-Dimensional Quantum Mixtures.

We consider two species of bosons in one dimension near the Tonks-Girardeau limit of infinite interactions. For the case of equal masses and equal intraspecies interactions, the system can be mapped to a S=1/2 XXZ Heisenberg spin chain, thus allowing one to access different magnetic phases. Using a powerful ansatz developed for the two-component Fermi system, we elucidate the evolution from few to many particles for the experimentally relevant case of an external harmonic confinement. In the few-body limit, we already find clear evidence of both ferromagnetic and antiferromagnetic spin correlations as the ratio of intraspecies and interspecies interactions is varied. Furthermore, we observe the rapid emergence of symmetry-broken magnetic ground states as the particle number is increased. We therefore demonstrate that systems containing only a few bosons are an ideal setting in which to realize the highly sought-after itinerant ferromagnetic phase.

Competing superconducting instabilities in the one-dimensionalp-band degenerate cold fermionic system

The zero-temperature phase diagram of $p$-orbital two-component fermionic system loaded into a one-dimensional optical lattice is mapped out by means of analytical and numerical techniques. It is shown that the $p$-band model away from half-filling hosts various competing superconducting phases for attractive and repulsive interactions. At quarter filling, we analyze the possible formation of incompressible Mott phases and in particular for repulsive interactions, we find the occurrence of a Mott transition with the formation of fully gapped bond-ordering waves.

Hyperspherical treatment of strongly-interacting few-fermion systems in one dimension

We examine a one-dimensional two-component fermionic system in a trap, assuming that all particles have the same mass and interact through a strong repulsive zero-range force. First we show how a simple system of three strongly interacting particles in a harmonic trap can be treated using the hyperspherical formalism. Next we discuss the behavior of the energy for the N-body system.

Sarma phase in relativistic and non-relativistic systems

We investigate the stability of the Sarma phase in two-component fermion systems in three spatial dimensions. For this purpose we compare strongly-correlated systems with either relativistic or non-relativistic dispersion relation: relativistic quarks and mesons at finite isospin density and spin-imbalanced ultracold Fermi gases. Using a Functional Renormalization Group approach, we resolve fluctuation effects onto the corresponding phase diagrams beyond the mean-field approximation. We find that fluctuations induce a second-order phase transition at zero temperature, and thus a Sarma phase, in the relativistic setup for large isospin chemical potential. This motivates the investigation of the cold atoms setup with comparable mean-field phase structure, where the Sarma phase could then be realized in experiment. However, for the non-relativistic system we find the stability region of the Sarma phase to be smaller than the one predicted from mean-field theory. It is limited to the BEC side of the phase diagram, and the unitary Fermi gas does not support a Sarma phase at zero temperature. Finally, we propose an ultracold quantum gas with four fermion species that has a good chance to realize a zero-temperature Sarma phase.

Pairing effects on spinodal decomposition of asymmetric nuclear matter

We present an analysis framed in the general context of two-component fermionic systems subjected to pairing correlations. The study is conducted for unstable asymmetric nuclear matter at low temperature, along the clusterization process driven by spinodal instabilities. It is shown that, especially around the transition temperature from the superfluid to the normal phase, pairing correlations may have non-negligible effects on the isotopic features of the clusterized low-density matter, which could be of interest also in the astrophysical context.