We show that self-localized ground states can be created in the spin-balanced gas of fermions with repulsion between the spin components, whose strength grows from the center to periphery, in combination with the harmonic-oscillator (HO) trapping potential acting in one or two transverse directions. We also consider the ground state in the noninteracting Fermi gas under the action of the spatially growing tightness of the one- or two-dimensional (1D or 2D) HO confinement. These settings are considered in the framework of the Thomas-Fermi--von Weizs\"acker (TF--vW) density functional. It is found that the vW correction to the simple TF approximation (the gradient term) is nearly negligible in all situations. The properties of the ground state under the action of the 2D and 1D HO confinement with the tightness growing in the transverse directions are investigated too for the Bose-Einstein condensate with the self-repulsive nonlinearity.
We show that self-localized ground states can be created in the spin-balanced gas of fermions with repulsion between the spin components, whose strength grows from the center to periphery, in combination with the harmonic-oscillator (HO) trapping potential acting in one or two transverse directions. We also consider the ground state in the noninteracting Fermi gas under the action of the spatially growing tightness of the one- or two-dimensional (1D or 2D) HO confinement. These settings are considered in the framework of the Thomas-Fermi--von Weizs"acker (TF--vW) density functional. It is found that the vW correction to the simple TF approximation (the gradient term) is nearly negligible in all situations. The properties of the ground state under the action of the 2D and 1D HO confinement with the tightness growing in the transverse directions are investigated too for the Bose-Einstein condensate with the self-repulsive nonlinearity.
We investigate separations of trapped balanced two-component atomic Fermi gases with repulsive contact interaction. Candidates for ground-state densities are obtained from the imaginary-time evolution of a nonlinear pseudo-Schr\"odinger equation in three dimensions, rather than from the cumbersome variational equations of the underlying energy density functional. With the employed hydrodynamical approach, gradient corrections to the Thomas-Fermi approximation are conveniently included and are shown to be vital for reliable density profiles. We provide critical repulsion strengths that mark the onset of phase transitions in a harmonic trap. We present transitions from identical density profiles of the two fermion species towards isotropic and anisotropic separations for various confinements, including harmonic and double-well-type traps. Our proposed method is suited for arbitrary trap geometries and can be straightforwardly extended to study dynamics in the light of ongoing experiments on degenerate Fermi gases.
We derive an effective nonpolynomial Schr\"odinger equation (NPSE) for self-repulsive or attractive BEC in the nearly one-dimensional cigar-shaped trap, with the transverse confining frequency periodically modulated along the axial direction. In addition to the usual linear cigar-shaped trap, where the periodic modulation emulates the action of an optical lattice (OL), the model may be also relevant to toroidal traps, where an ordinary OL cannot be created. For either sign of the nonlinearity, extended and localized states are found, in the numerical form [using both the effective NPSE and the full three-dimensional (3D) Gross-Pitaevskii equation] and by means of the variational approximation (VA). The latter is applied to construct ground-state solitons and predict the collapse threshold in the case of self-attraction. It is shown that numerical solutions provided by the one-dimensional NPSE are always very close to full 3D solutions, and the VA yields quite reasonable results too. The transition from delocalized states to gap solitons, in the first finite bandgap of the linear spectrum, is examined in detail, for the repulsive and attractive nonlinearities alike.
We investigate the ground-state properties of two-component Bose gases confined in one-dimensional harmonic traps in the scheme of density-functional theory. The density-functional calculations employ a Bethe-ansatz-based local-density approximation for the correlation energy, which accounts for the correlation effect properly from the weakly interacting regime to the strongly interacting regime. For the binary Bose mixture with spin-independent interaction, the homogeneous reference system is exactly solvable by the Bethe-ansatz method. Within the local-density approximation, we determine the density distribution of each component and study its evolution from Bose distributions to Fermi-like distribution with the increase in interaction. For the binary mixture of Tonks-Girardeau gases with a tunable interspecies repulsion, with a generalized Bose-Fermi transformation we show that the Bose mixture can be mapped into a two-component Fermi gas, which corresponds to exact soluble Yang-Gaudin model for the homogeneous system. Based on the ground-state energy function of the Yang-Gaudin model, the ground-state density distributions are calculated for various interspecies interactions. It is shown that with the increase in interspecies interaction, the system exhibits composite-fermionization crossover.
We consider a strongly interacting one-dimensional (1D) Bose-Fermi mixture confined in a harmonic trap. It consists of a Tonks-Girardeau (TG) gas (1D Bose gas with repulsive hard-core interactions) and of a noninteracting Fermi gas (1D spin-aligned Fermi gas), both species interacting through hard-core repulsive interactions. Using a generalized Bose-Fermi mapping, we determine the one-body density matrices, exact particle density profiles, momentum distributions, and behavior of the mixture under 1D expansion when opening the trap. In real space, bosons and fermions do not display any phase separation: the respective density profiles extend over the same region and they both present a number of peaks equal to the total number of particles in the trap. In momentum space the bosonic component has the typical narrow TG profile, while the fermionic component shows a broad distribution with fermionic oscillations at small momenta. Due to the large boson-fermion repulsive interactions, both the bosonic and the fermionic momentum distributions decay as $C{p}^{\ensuremath{-}4}$ at large momenta, like in the case of a pure bosonic TG gas. The coefficient C is related to the two-body density matrix and to the bosonic concentration in the mixture. When opening the trap, both momentum distributions ``fermionize'' under expansion and turn into that of a Fermi gas with a particle number equal to the total number of particles in the mixture.
We investigate the ground-state properties of ultracold two-component Fermi gases in the presence of a transverse harmonic potential, focusing on the strongly interacting regime in which pairs of fermions form tightly bound molecules. Using the fixed-node diffusion Monte Carlo method, we calculate the equation of state and density profiles for the full fermionic system, which allows us to address the importance of finite-range corrections arising from the internal fermionic structure of the dimers. We interpret the results in terms of a molecular Bose gas in quasi-two-dimensional confinement and compare them with theoretical predictions for a weakly interacting two-dimensional Bose gas, identifying the range of validity of mean-field and beyond-mean-field descriptions. We also develop an analytical theory for the transverse density profile, capturing its broadening with increasing interaction strength. This work provides a benchmark for an effective bosonic description of strongly bound fermionic dimers and offers insights into confinement-induced dimensional effects.
We develop a method for creating two- and one-dimensional (2D and 1D) self-trapped modes in binary spin-orbit-coupled Bose–Einstein condensates with the contact repulsive interaction, whose local strength grows sufficiently rapidly from the center to the periphery. In particular, an exact semi-vortex (SV) solution is found for the anti-Gaussian radial modulation profile. The exact modes are included in the numerically produced family of SV solitons. Other families, in the form of mixed modes (MMs), as well as excited states of SVs and MMs, are also produced. Although the excited states are unstable in all previously studied models, they are partially stable in the present one. In the 1D version of the system, exact solutions for the counterpart of SVs, namely, semi-dipole solitons, are also found. Families of semi-dipoles, as well as the 1D version of MMs, are produced numerically.