Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEfficient simulation of moire materials using the density matrix renormalization group
104
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present an infinite density-matrix renormalization group (DMRG) study of an interacting continuum model of twisted bilayer graphene (tBLG) near the magic angle. Because of the long-range Coulomb interaction and the large number of orbital degrees of freedom, tBLG is difficult to study with standard DMRG techniques -- even constructing and storing the Hamiltonian already poses a major challenge. To overcome these difficulties, we use a recently developed compression procedure to obtain a matrix product operator representation of the interacting tBLG Hamiltonian which we show is both efficient and accurate even when including the spin, valley and orbital degrees of freedom. To benchmark our approach, we focus mainly on the spinless, single-valley version of the problem where, at half-filling, we find that the ground state is a nematic semimetal. Remarkably, we find that the ground state is essentially a k-space Slater determinant, so that Hartree-Fock and DMRG give virtually identical results for this problem. Our results show that the effects of long-range interactions in magic angle graphene can be efficiently simulated with DMRG, and opens up a new route for numerically studying strong correlation physics in spinful, two-valley tBLG, and other moire materials, in future work.
rate research
Read More
We introduce the transcorrelated Density Matrix Renormalization Group (tcDMRG) theory for the efficient approximation of the energy for strongly correlated systems. tcDMRG encodes the wave function as a product of a fixed Jastrow or Gutzwiller correlator and a matrix product state. The latter is optimized by applying the imaginary-time variant of time-dependent (TD) DMRG to the non-Hermitian transcorrelated Hamiltonian. We demonstrate the efficiency of tcDMRG at the example of the two-dimensional Fermi-Hubbard Hamiltonian, a notoriously difficult target for the DMRG algorithm, for different sizes, occupation numbers, and interaction strengths. We demonstrate fast energy convergence of tcDMRG, which indicates that tcDMRG could increase the efficiency of standard DMRG beyond quasi-monodimensional systems and provides a generally powerful approach toward the dynamic correlation problem of DMRG.
In some cases the state of a quantum system with a large number of subsystems can be approximated efficiently by the density matrix renormalization group, which makes use of redundancies in the description of the state. Here we show that the achievable efficiency can be much better when performing density matrix renormalization group calculations in the Heisenberg picture, as only the observable of interest but not the entire state is considered. In some non-trivial cases, this approach can even be exact for finite bond dimensions.
A distributed-memory parallelization strategy for the density matrix renormalization group is proposed for cases where correlation functions are required. This new strategy has substantial improvements with respect to previous works. A scalability analysis shows an overall serial fraction of 9.4% and an efficiency of around 60% considering up to eight nodes. Sources of possible parallel slowdown are pointed out and solutions to circumvent these issues are brought forward in order to achieve a better performance.
We investigate the application of the Density Matrix Renormalization Group (DMRG) to the Hubbard model in momentum-space. We treat the one-dimensional models with dispersion relations corresponding to nearest-neighbor hopping and $1/r$ hopping and the two-dimensional model with isotropic nearest-neighbor hopping. By comparing with the exact solutions for both one-dimensional models and with exact diagonalization in two dimensions, we first investigate the convergence of the ground-state energy. We find variational convergence of the energy with the number of states kept for all models and parameter sets. In contrast to the real-space algorithm, the accuracy becomes rapidly worse with increasing interaction and is not significantly better at half filling. We compare the results for different dispersion relations at fixed interaction strength over bandwidth and find that extending the range of the hopping in one dimension has little effect, but that changing the dimensionality from one to two leads to lower accuracy at weak to moderate interaction strength. In the one-dimensional models at half-filling, we also investigate the behavior of the single-particle gap, the dispersion of spinon excitations, and the momentum distribution function. For the single-particle gap, we find that proper extrapolation in the number of states kept is important. For the spinon dispersion, we find that good agreement with the exact forms can be achieved at weak coupling if the large momentum-dependent finite-size effects are taken into account for nearest-neighbor hopping. For the momentum distribution, we compare with various weak-coupling and strong-coupling approximations and discuss the importance of finite-size effects as well as the accuracy of the DMRG.
We present a method for computing resonant inelastic x-ray scattering (RIXS) spectra in one-dimensional systems using the density matrix renormalization group (DMRG) method. By using DMRG to address the problem, we shift the computational bottleneck from the memory requirements associated with exact diagonalization (ED) calculations to the computational time associated with the DMRG algorithm. This approach is then used to obtain RIXS spectra on cluster sizes well beyond state-of-the-art ED techniques. Using this new procedure, we compute the low-energy magnetic excitations observed in Cu $L$-edge RIXS for the challenging corner shared CuO$_4$ chains, both for large multi-orbital clusters and downfolded $t$-$J$ chains. We are able to directly compare results obtained from both models defined in clusters with identical momentum resolution. In the strong coupling limit, we find that the downfolded $t$-$J$ model captures the features of the magnetic excitations probed by RIXS after a uniform scaling of the spectra is taken into account.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
