We consider repulsively-interacting cold fermionic atoms loaded on an optical ladder lattice in a trapping potential. The density-matrix renormalization-group method is used to numerically calculate the ground state for systematically varied values o
f interaction U and spin imbalance p in the Hubbard model on the ladder. The system exhibits rich structures, where a fully spin polarized phase, spatially separated from other domains in the trapping potential, appears for large enough U and p. The phase-separated ferromagnetism can be captured as a real-space image of the energy gap between the ferromagnetic and other states arising from a combined effect of Nagaokas ferromagnetism extended to the ladder and the density dependence of the energy separation between competing states. We also predict how to maximize the ferromagnetic region.
In order to demonstrate that atomic Fermi gas is a good experimental reality in studying unsolved problems in frustrated interacting-spin systems, we numerically examine the Mott core state emerged by loading two-component atomic Fermi gases on trian
gular optical lattices. Consequently, we find that plateau like structures are observable in the Mott core polarization as a function of the population imbalance. These plateau states are caused by a flexibility that the surrounding metallic region absorbs the excess imbalance to keep the plateau states inside the Mott core. We also find spin patterns peculiar to the plateau states inside the Mott core.
We investigate effects of pseudo-spin population imbalance on Mott phases in 1D trapped two-component atomic Fermi gases loaded on optical lattices based on the repulsive Hubbard model in harmonic traps. By using the density matrix renormalization gr
oup method, we numerically calculate density profiles of each component and clarify the pseudo-spin magnetism. Consequently, we find that all the features from weakly imbalance to fully polarized cases are well described by S=1/2 antiferromagnetic Heisenberg chain under magnetic field. These results indicate that the Mott phases offer experimental stages for studying various interacting spin systems.
We study the interplay of disorder and correlation in the one-dimensional hole-doped Hubbard-model with disorder (Anderson-Hubbard model) by using the density-matrix renormalization group method. Concentrating on the doped-hole density profile, we fi
nd in a large $U/t$ regime that the clean system exhibits a simple fluid-like behavior whereas finite disorders create locally Mott regions which expand their area with increasing the disorder strength contrary to the ordinary sense. We propose that such an anomalous Mott phase formation assisted by disorder is observable in atomic Fermi gases by setup of the box shape trap.
In order to study an interplay of disorder, correlation, and spin imbalance on antiferromagnetism, we systematically explore the ground state of one-dimensional spin-imbalanced Anderson-Hubbard model by using the density-matrix renormalization group
method. We find that disorders localize the antiferromagnetic spin density wave induced by imbalanced fermions and the increase of the disorder magnitude shrinks the areas of the localized antiferromagnetized regions. Moreover, the antiferromagnetism finally disappears above a large disorder. These behaviors are observable in atomic Fermi gases loaded on optical lattices and disordered strongly-correlated chains under magnetic field.
We investigate effects of optical lattice potential in one- and two-dimensional two-component trapped Fermi gases with population imbalances. Using the exact diagonalization and the density matrix renormalization group methods complementarily, we cal
culate the atom density profile from the ground state many-body wavefunction as a function of attractive interaction strength for various population imbalances. The numerical results reveal that although a phase separation between the superfluid core and the shell cloud of excess atoms occurs as observed in experiments without the optical lattice, the population imbalance generally remains in the core region in contrast to the non-lattice cases. The essence of the numerical results in a strong attractive regime can be explained by an effective model composed of Cooper pairs and excess major fermions.
We suggest that box shape trap enables to observe intrinsic properties of the repulsive Hubbard model in a fixed doping in contrast to the harmonic trap bringing about spatial variations of atom density profiles. In order to predict atomic density pr
ofile under the box trap, we apply the directly-extended density-matrix renormalization group method to 4-leg repulsive Hubbard model with the open boundary condition. Consequently, we find that stripe formation is universal in a low hole doping range and the stripe sensitively changes its structure with variations of $U/t$ and the doping rate. A remarkable change is that a stripe formed by a hole pair turns to one by a bi-hole pair when entering a limited strong $U/t$ range. Furthermore, a systematic calculation reveals that the Hubbard model shows a change from the stripe to the Friedel like oscillation with increasing the doping rate.
We parallelize density-matrix renormalization group to directly extend it to 2-dimensional ($n$-leg) quantum lattice models. The parallelization is made mainly on the exact diagonalization for the superblock Hamiltonian since the part requires an eno
rmous memory space as the leg number $n$ increases. The superblock Hamiltonian is divided into three parts, and the correspondent superblock vector is transformed into a matrix, whose elements are uniformly distributed into processors. The parallel efficiency shows a high rate as the number of the states kept $m$ increases, and the eigenvalue converges within only a few sweeps in contrast to the multichain algorithm.
