The minimally entangled typical thermal states algorithm is applied to fermionic systems using the Krylov-space approach to evolve the system in imaginary time. The convergence of local observables is studied in a tight-binding system with a site-dep
endent potential. The temperature dependence of the superconducting correlations of the attractive Hubbard model is analyzed on chains, showing an exponential decay with distance and exponents proportional to the temperature at low temperatures, as expected. In addition, the non-local parity correlator is calculated at finite temperature. Other possible applications of the minimally entangled typical thermal states algorithm to fermionic systems are also discussed.
A detailed description of the time-step-targetting time evolution method within the DMRG algorithm is presented. The focus of this publication is on the implementation of the algorithm, and on its generic application. The case of one-site excitations
within a Hubbard model is analyzed as a test for the algorithm, using open chains and two-leg ladder geometries. The accuracy of the procedure in the case of the recently discussed holon-doublon photo excitations of Mott insulators is also analyzed. Performance and parallelization issues are discussed. In addition, the full open-source code is provided as supplementary material.
In the Density Matrix Renormalization Group (DMRG) algorithm, Hamiltonian symmetries play an important role. Using symmetries, the matrix representation of the Hamiltonian can be blocked. Diagonalizing each matrix block is more efficient than diagona
lizing the original matrix. This paper explains how the the DMRG++ code has been extended to handle the non-local SU(2) symmetry in a model independent way. Improvements in CPU times compared to runs with only local symmetries are discussed for the one-orbital Hubbard model, and for a two-orbital Hubbard model for iron-based superconductors. The computational bottleneck of the algorithm and the use of shared memory parallelization are also addressed.
A quasi-classical theoretical description of polarization and relaxation of nuclear spins in a quantum dot with one resident electron is developed for arbitrary mechanisms of electron spin polarization. The dependence of the electron-nuclear spin dyn
amics on the correlation time $tau_c$ of electron spin precession, with frequency $Omega$, in the nuclear hyperfine field is analyzed. It is demonstrated that the highest nuclear polarization is achieved for a correlation time close to the period of electron spin precession in the nuclear field. For these and larger correlation times, the indirect hyperfine field, which acts on nuclear spins, also reaches a maximum. This maximum is of the order of the dipole-dipole magnetic field that nuclei create on each other. This value is non-zero even if the average electron polarization vanishes. It is shown that the transition from short correlation time to $Omegatau_c>1$ does not affect the general structure of the equation for nuclear spin temperature and nuclear polarization in the Knight field, but changes the values of parameters, which now become functions of $Omegatau_c$. For correlation times larger than the precession time of nuclei in the electron hyperfine field, it is found that three thermodynamic potentials ($chi$, $bm{xi}$, $varsigma$) characterize the polarized electron-nuclear spin system. The values of these potentials are calculated assuming a sharp transition from short to long correlation times, and the relaxation mechanisms of these potentials are discussed. The relaxation of the nuclear spin potential is simulated numerically showing that high nuclear polarization decreases relaxation rate.
The purpose of this paper is (i) to present a generic and fully functional implementation of the density-matrix renormalization group (DMRG) algorithm, and (ii) to describe how to write additional strongly-correlated electron models and geometries by
using templated classes. Besides considering general models and geometries, the code implements Hamiltonian symmetries in a generic way and parallelization over symmetry-related matrix blocks.
The one-particle spectral function of a state formed by superconducting (SC) clusters is studied via Monte Carlo techniques. The clusters have similar SC amplitudes but randomly distributed phases. This state is stabilized by the competition with ant
i-ferromagnetism, after quenched disorder is introduced. Fermi arcs between the critical temperature Tc and the cluster formation temperature scale T* are observed, similarly as in the pseudo-gap state of the cuprates. The arcs originate at metallic regions in between the neighboring clusters that present large SC phase differences.
