We introduce a computational approach to study porphyrin-like transition metal complexes, bridging density functional theory and exact many-body techniques, such as the density matrix renormalization group (DMRG). We first derive a multi-orbital Anderson impurity Hamiltonian starting from first principles considerations that qualitatively reproduce generalized gradient approximation (GGA)+U results when ignoring inter-orbital Coulomb repulsion U ′ and Hund exchange J. An exact canonical transformation is used to reduce the dimensionality of the problem and make it amenable to DMRG calculations, including all many-body terms (both intra- and inter-orbital), which are treated in a numerically exact way. We apply this technique to FeN 4 centers in graphene and show that the inclusion of these terms has dramatic effects: as the iron orbitals become single occupied due to the Coulomb repulsion, the inter-orbital interaction further reduces the occupation, yielding a non-monotonic behavior of the magnetic moment as a function of the interactions, with maximum polarization only in a small window at intermediate values of the parameters. Furthermore, U ′ changes the relative position of the peaks in the density of states, particularly on the iron d z 2 orbital, which is expected to affect the binding of ligands greatly.
We introduce a computational approach to study porphyrin-like transition metal complexes, bridging density functional theory and exact many-body techniques, such as the density matrix renormalization group (DMRG). We first derive a multi-orbital Anderson impurity Hamiltonian starting from first principles considerations that qualitatively reproduce generalized gradient approximation (GGA)+U results when ignoring inter-orbital Coulomb repulsion U ′ and Hund exchange J. An exact canonical transformation is used to reduce the dimensionality of the problem and make it amenable to DMRG calculations, including all many-body terms (both intra- and inter-orbital), which are treated in a numerically exact way. We apply this technique to FeN 4 centers in graphene and show that the inclusion of these terms has dramatic effects: as the iron orbitals become single occupied due to the Coulomb repulsion, the inter-orbital interaction further reduces the occupation, yielding a non-monotonic behavior of the magnetic moment as a function of the interactions, with maximum polarization only in a small window at intermediate values of the parameters. Furthermore, U ′ changes the relative position of the peaks in the density of states, particularly on the iron d z 2 orbital, which is expected to affect the binding of ligands greatly.
We introduce a block Lanczos (BL) recursive technique to construct quasi-one-dimensional models, suitable for density-matrix renormalization group (DMRG) calculations, from single- as well as multiple-impurity Anderson models in any spatial dimensions. This new scheme, named BL-DMRG method, allows us to calculate not only local but also spatially dependent static and dynamical quantities of the ground state for general Anderson impurity models without losing elaborate geometrical information of the lattice. We show that the BL-DMRG method can be easily extended to treat a multiorbital Anderson impurity model where not only inter- and intraorbital Coulomb interactions but also Hund's coupling and pair hopping interactions are included. We also show that the symmetry adapted BL bases can be utilized, when it is appropriate, to reduce the computational cost. As a demonstration, we apply the BL-DMRG method to three different models for graphene with a structural defect and with a single hydrogen or fluorine absorbed, where a single Anderson impurity is coupled to conduction electrons in the honeycomb lattice. These models include (i) a single adatom on the honeycomb lattice, (ii) a substitutional impurity in the honeycomb lattice, and (iii) an effective model for a single carbon vacancy in graphene. Our analysis of the local dynamical magnetic susceptibility and the local density of states at the impurity site reveals that, for the particle-hole symmetric case at half-filling of electron density, the ground state of model (i) behaves as an isolated magnetic impurity with no Kondo screening, while the ground state of the other two models forms a spin-singlet state where the impurity moment is screened by the conduction electrons. We also calculate the real-space dependence of the spin-spin correlation functions between the impurity site and the conduction sites for these three models. Our results clearly show that, reflecting the presence or absence of unscreened magnetic moment at the impurity site, the spin-spin correlation functions decay as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}3}$, differently from the noninteracting limit ($\ensuremath{\propto}$ ${r}^{\ensuremath{-}2}$), for model (i) and as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}4}$, exactly the same as the noninteracting limit, for models (ii) and (iii) in the asymptotic r, where r is the distance between the impurity site and the conduction site. Finally, based on our results, we shed light on recent experiments on graphene where the formation of local magnetic moments as well as the Kondo-like behavior have been observed.
We present a completely unbiased and controlled numerical method to solve quantum impurity problems in d-dimensional lattices. This approach is based on a canonical transformation, of the Lanczos form, where the complete lattice Hamiltonian is exactly mapped onto an equivalent one-dimensional system, in the same spirit as Wilson's numerical renormalization, and Haydock's recursion method. We introduce many-body interactions in the form of a Kondo or Anderson impurity and we solve the low-dimensional problem using the density matrix renormalization group. The technique is particularly suited to study systems that are inhomogeneous, and/or have a boundary. The resulting dimensional reduction translates into a reduction of the scaling of the entanglement entropy by a factor ${L}^{d\ensuremath{-}1}$, where L is the linear dimension of the original d-dimensional lattice. This allows one to calculate the ground state of a magnetic impurity attached to an $L\ifmmode\times\else\texttimes\fi{}L$ square lattice and an $L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L$ cubic lattice with L up to 140 sites. We also study the localized edge states in graphene nanoribbons by attaching a magnetic impurity to the edge or the center of the system. For armchair metallic nanoribbons we find a slow decay of the spin correlations as a consequence of the delocalized metallic states. In the case of zigzag ribbons, the decay of the spin correlations depends on the position of the impurity. If the impurity is situated in the bulk of the ribbon, the decay is slow as in the metallic case. On the other hand, if the adatom is attached to the edge, the decay is fast, within few sites of the impurity, as a consequence of the localized edge states, and the short correlation length. The mapping can be combined with ab initio band structure calculations to model the system, and to understand correlation effects in quantum impurity problems starting from first principles.
Using the density-matrix renormalization group method and quantum Monte Carlo simulation, we studied numerically the energy spectrum and thermodynamics of the quantum Heisenberg spin model for narrow graphene nanoribbons and their derivatives with periodically embedded heteroatoms. For several nanoribbon structures we found macroscopic ground state spin, gapless lowest excitation spectra and intermediate magnetization plateaus at low temperatures. We also studied the lowest energy states of frustrated systems formed by triangular graphitic clusters connected by bridged ions of transition metals. On the base of many-body perturbation theory and the exact diagonalization method, we showed the possibility of spin switching for this model due to the change the corresponding coupling parameters.
We investigated the honeycomb lattice using a two-orbital (α and β) Hubbard model and many-body exact calculations. The spin gap (Δ s ), charge gap (Δ c ), intrasite ( γ i i α β = ⟨ S → i α · S → i β ⟩ ), intersite ( γ i j α β = ⟨ S → i α · S → j β ⟩ ), spin–spin correlations, and orbital magnetic moments μ iα are analyzed for different band-filling ( n ), which denotes the number of electrons per site. Special attention was set in one-quarter-filled ( n = 1) and half-filled ( n = 2) bands. The results are presented by varying Coulomb repulsion ( U), Hund’s coupling ( J H ), external magnetic field ( B ), and orbital energy splitting (Δϵ). For the one-quarter-filled band, an intrasite magnetic disorder ( γ i i α β ≈ 0 ) begins to form for large values of U / t, which is interpreted as a possible signature of a resonating valence bond state. Full orbital polarization ( γ i i α β = 0.25 ) with intersite antiferromagnetic (AFM) order ( γ i j α β = − 0.25 ) was observed in the half-filled band for U / t ≫ 1, indicating that the two-orbital Hubbard model can be mapped to an AFM Heisenberg model with S = 1. The increase of J H is conductive to insulating and metallic phases for half- and one-quarter-filled bands, respectively. For a fractional n, the honeycomb lattice can exhibit spin-glass-like behavior, predominating as a disordered magnetic system characterized by random values of γ i j α β . The magnetic field induced a reduction in Δ c, which became zero at the critical external magnetic field. Twisted bilayer graphene (TBG) exhibits approximately flat electronic bands at magic angles and interesting correlated phenomena, such as unconventional superconductivity and correlated insulating states around n = 1 and n = 2, respectively. The dominance of electron–electron interactions in TBG makes the Hubbard model, in various versions (extended or multi-orbital), offer a practical, structure-scale account for investigating many-body effects in large moiré unit cells, where first-principles density-functional theory calculations are computationally challenging.
We propose a quantum gate architecture that allows for the systematic control of the effective exchange interactions between magnetic impurities embedded in nanoscale graphene flakes connected by a gated bridge. The entanglement between the magnetic moment and the edge states of the fragments is used to electrostatically tune the exchange interaction from ferro to antiferromagnetic by merely changing the bridge's carrier density. By characterizing the effects of size and coupling parameters, we explore different operation regimes of this device by means of exact calculations with the density matrix renormalization group. We analyze the results utilizing a simplified model that accounts for the main many-body mechanisms. Finally, we discuss how to use arrays of these devices to build quantum simulators for quantum many-body Hamiltonians.
The applicability of ab initio multireference wavefunction-based methods to the study of magnetic complexes has been restricted by the quickly rising active-space requirements of oligonuclear systems and dinuclear complexes with S > 1 spin centers. Ab initio density matrix renormalization group (DMRG) methods built upon an efficient parameterization of the correlation network enable the use of much larger active spaces, and therefore may offer a way forward. Here, we apply DMRG-CASSCF to the dinuclear complexes [Fe2OCl6](2-) and [Cr2O(NH3)10](4+). After developing the methodology through systematic basis set and DMRG M testing, we explore the effects of extended active spaces that are beyond the limit of conventional methods. We find that DMRG-CASSCF with active spaces including the metal d orbitals, occupied bridging-ligand orbitals, and their virtual double shells already capture a major portion of the dynamic correlation effects, accurately reproducing the experimental magnetic coupling constant (J) of [Fe2OCl6](2-) with (16e,26o), and considerably improving the smaller active space results for [Cr2O(NH3)10](4+) with (12e,32o). For comparison, we perform conventional MRCI+Q calculations and find the J values to be consistent with those from DMRG-CASSCF. In contrast to previous studies, the higher spin states of the two systems show similar deviations from the Heisenberg spectrum, regardless of the computational method.