Search

YMn$_6$Sn$_6$ consists of two types of Mn-based kagome planes stacked along $c$-axis having a complex magnetic interactions. We report a spin reconstruction in YMn$_6$Sn$_6$ from ferromagnet (FM) into a combination of two incommensurate spin spirals (SSs) originating from two different type of Mn kagome planes driven by frustrated magnetic exchanges along the $c$-axis with inclusion of Hubbard U. The pitch angle and wave vector of the incommensurate SSs are $\sim$89.3$^{\circ}$ and $\sim$ (0 0 0.248) respectively which are in excellent agreement with experiment. We employed an effective model Hamiltonian constructed out of exchange interactions to capture experimentally observed non-equivalent nature of the two incommensurate SSs which also explain FM-SS crossover due to antiferromagntic spin exchange with correlation. We further report the existence of topological magnon with spin-orbit coupling in incommensurate SS phase of YMn$_6$Sn$_6$ by calculating the topological invariants and Berry curvature profile. The location of Dirac magnon in energy landscape at 73 meV matches with another experimental report. We demonstrate the accuracy of our results by highlighting experimental features in YMn$_6$Sn$_6$.

The theoretical calculation of the interaction of electromagnetic radiation with matter remains a challenging problem for contemporary $ab$ $initio$ electronic structure methods, in particular for x-ray spectroscopies. This is not only due to the strong interaction between the core-hole and the photo-excited electron, but also due to the elusive multiplet effects that arise from the Coulomb interaction among the valence electrons. In this work we report a method based on density-functional theory in conjunction with multiplet ligand-field theory to investigate various core-level spectroscopies, in particular x-ray absorption spectroscopy (XAS) and x-ray magnetic circular dichroism (XMCD). The developed computational scheme is applied to the $L_{2,3}$ XAS edges of magnetite (Fe$_3$O$_4$) as well as cobalt ferrite (CoFe$_2$O$_4$) and nickel ferrite (NiFe$_2$O$_4$) and the corresponding XMCD spectra. The results are in overall good agreement with experimental observations, both regarding the XAS $L_2$/$L_3$ branching ratio, the peak positions as well as the relative intensities. The agreement between theory and experiment is equally good for XAS and the XMCD spectra, for all studied systems. The results are analyzed in terms of $e_g$ and $t_{2g}$ orbitals contributions and the importance of optimizing the Slater parameters. The analysis also highlights the strong effect of the $2p$-$3d$ interaction in x-ray spectroscopy., Comment: 18 pages, 15 figures

PdFe/Ir(111) has attracted tremendous attention for next-generation spintronics devices due to existence of magnetic skyrmions with the external magnetic field. Our density functional theoretical calculations in combination with spin dynamics simulation suggest that the spin spiral phase in fcc stacked PdFe/Ir(111) flips into the skyrmion lattice phase around B$_{ext} \sim$ 6 T. This leads to the microscopic understanding of the thermodynamic and kinetic behaviours affected by the intrinsic spin-lattice couplings (SLCs) in this skyrmion material for magneto-mechanical properties. Here we calculate fully relativistic SLC parameters from first principle computations and investigate the effect of SLC on dynamical magnetic interactions in skyrmion multilayers PdFe/Ir(111). The exchange interactions arising from next nearest-neighbors (NN) in this material are highly frustrated and responsible for enhancing skyrmion stability. We report the larger spin-lattice effect on both dynamical Heisenberg exchanges and Dzyaloshinskii-Moriya interactions for next NN compared to NN which is in contrast with recently observed spin-lattice effect in bulk bcc Fe and CrI$_3$ monolayer. Based on our analysis, we find that the effective measures of SLCs in fcc (hcp) stacking of PdFe/Ir(111) are $\sim 2.71 ( \sim 2.36)$ and $\sim 14.71 ( \sim21.89)$ times stronger for NN and next NN respectively, compared to bcc Fe. The linear regime of displacement for SLC parameters is $\leq$ 0.02 {\AA} which is 0.72\% of the lattice constant for PdFe/Ir(111). The microscopic understanding of SLCs provided by our current study could help in designing spintronic devices based on thermodynamic properties of skyrmion multilayers., Comment: 8(main text)+4(appendix) pages and 5(main text)+4(appendix) figures

Chromium triiodide is an intrinsically magnetic van der Waals material down to the single-layer limit. Here, we provide a first-principles description of finite-temperature magnetic and spectral properties of monolayer (ML) CrI$_3$ based on fully charge self-consistent density functional theory (DFT) combined with dynamical mean-field theory, revealing a formation of local moments on Cr from strong local Coulomb interactions. We show that the presence of local dynamical correlations leads to a modification of the electronic structure of ferromagnetically ordered CrI$_3$. In contrast to conventional DFT+$U$ calculations, we find that the top of the valence band in ML CrI$_3$ demonstrates essentially different orbital character for minority and majority spin states, which is closer to the standard DFT results. This leads to a strong spin polarization of the optical conductivity upon hole doping, which could be verified experimentally., Comment: Final version; 9 pages, 4 figures

The Dzyaloshinskii-Moriya (DM) interaction, as well as symmetric anisotropic exchange, are important ingredients for stabilizing topologically non-trivial magnetic textures, such as, e.g., skyrmions, merons and hopfions. These types of textures are currently in focus from a fundamental science perspective and they are also discussed in the context of future spintronics information technology. While the theoretical understanding of the Heisenberg exchange interactions is well developed, it is still a challenge to access, from first principles theory, the DM interaction as well as the symmetric anisotropic exchange, which both require a fully-relativistic treatment of the electronic structure, in magnetic systems where substantial electron-electron correlations are present. Here, we present results of a theoretical framework which allows to compute these interactions in any given system and demonstrate its performance for several selected cases, for both bulk and low-dimensional systems. We address several representative cases, including the bulk systems CoPt and FePt, the B20 compounds MnSi and FeGe as well as the low-dimensional transition metal bilayers Co/Pt(111) and Mn/W(001). The effect of electron-electron correlations is analyzed using dynamical mean-field theory on the level of the spin-polarized $T$-matrix + fluctuating exchange (SPTF) approximation, as regards the strength and character of the isotropic (Heisenberg) and anisotropic (DM) interactions in relation to the underlying electronic structure. Our method can be combined with more advanced techniques for treating correlations, e.g., quantum Monte Carlo and exact diagonalization methods for the impurity solver of dynamical mean-field theory. We find that correlation-induced changes of the DM interaction can be rather significant, with up to five-fold modifications in the most distinctive case., Comment: 20 pages, 22 figures

With progress towards more compact quantum computing architectures, fundamental questions regarding the entanglement of indistinguishable particles need to be addressed. In a solid state device, this quest is naturally connected to the quantum correlations of electrons. Here, we investigate the entanglement between electrons, focusing on the entanglement of modes, the entanglement of particles and the effect of particle-number superselection rules. We elucidate the formation of mode and particle entanglement in strongly correlated materials and show that both represent important resources in quantum information tasks such as quantum teleportation. To this end, we qualitatively and quantitatively analyze the entanglement in three electronic teleportation schemes: (i) quantum teleportation within a molecule on graphene, (ii) a nitrogen-vacancy center and (iii) a quantum dot array.

Whether SmB6 is a topological insulator remains hotly debated. Our density functional theory plus dynamical mean field theory calculations are in excellent agreement with a large range of experiments, from the 4f5.5 intermediate valency to x-ray and photoemission spectra. Using the pole extended (PE) Hamiltonian, which fully captures the self-energy, we show that SmB6 is a strongly correlated topological insulator, albeit not a Kondo insulator. The PE Hamiltonian is proved to be topological (in)equivalent to the topological Hamiltonian for (non-)local self-energies. The topological surface states are analyzed, addressing conflicting interpretations of photoemission data., Comment: The Supplemental Material starts at page 7

The interaction of light with solids gives rise to new bosonic quasiparticles, with the exciton being---undoubtedly---the most famous of these polaritons. While excitons are the generic polaritons of semiconductors, we show that for strongly correlated systems another polariton is prevalent---originating from the dominant antiferromagnetic or charge density wave fluctuations in these systems. As these are usually associated with a wave vector ${\mathbf k}= (\pi,\pi,\ldots)$ or close to it, we propose to call the derived polaritons $\pi$-tons. These $\pi$-tons yield the leading vertex correction to the optical conductivity in all correlated models studied: the Hubbard, the extended Hubbard model, the Falicov-Kimball, and the Pariser-Parr-Pople model, both in the insulating and in the metallic phase., Comment: 6 pages, 3 figures + Supplemental Material (6 pages, 8 figures)

We study a simple model system for the conjugated $\pi$-bonds in benzene, the Pariser-Parr-Pople (PPP) model, within the parquet approximation (PA), exemplifying the prospects of the PA for molecules. Advantages of the PA are the polynomial scaling with the number of orbitals, and the natural calculation of one- and two-particle spectral functions as well as of response and correlation functions. We find large differences in the electronic correlations in the PPP model compared to a Hubbard model with only local interactions. The quasiparticle renormalization (or mass enhancement) is much weaker in the PPP than in the Hubbard model, but the static part of the self-energy enhances the band gap of the former. Furthermore, the vertex corrections to the optical conductivity are much more important in the PPP model. Because non-local interactions strongly alter the self-energy, we conclude that the PA is more suitable for calculating conjugated $\pi$-bonds in molecules than single site dynamical mean-field theory., Comment: 14 pages, 7 figures

We have performed a systematic study of the emergence of meta-stable states in density functional theory plus Hubbard U (DFT+U ) simulations of NiO, CoO, FeO. Particular attention is given to the spin-polarization of the exchange-correlation functional and the double counting term, and the role of the spin-orbit coupling. The method of occupation matrix control is extended to use constrained random density matrices to map out the local minima in the total energy landscape. The extended scheme, random density matrix control, is successfully benchmarked against UO2, one of the most investigated systems in the field. When applied to the transition metal oxides it yields several meta- stable states which are well-characterized by their local spin and orbital moments. We find that the addition of spin-orbit coupling helps the simulations to converge to the global high-spin energy minimum. The random density matrix control scheme combined with LDA+U yields accurate magnetic moments for all the studied AFM transition metal oxides., Comment: 14 figures

Two-particle generalized susceptibilities and their irreducible vertex functions play a prominent role in the quantum many-body theory for correlated electron systems. They act as basic building blocks in the parquet formalism which provides a flexible scheme for the calculation of spectral and response functions. The irreducible vertices themselves have recently attracted increased attention as unexpected divergences in these functions have been identified. Remarkably, such singularities appear already for one of the simplest strongly interacting systems: the atomic limit of the half-filled Hubbard model (Hubbard atom). In this paper, we calculate the analytical expressions for all two-particle irreducible vertex functions of the Hubbard atom in all scattering channels as well as the fully irreducible two-particle vertices. We discuss their divergences and classify them by the eigenvalues and eigenvectors of the corresponding generalized susceptibilities. In order to establish a connection to the recently found multivaluedness of the exact self-energy functional $\Sigma[G]$, we show that already an approximation akin to iterated perturbation theory is sufficient to capture, qualitatively, the divergent structure of the vertex functions. Finally, we show that the localized divergences in the disordered binary mixture model are directly linked to a minimum in the single-particle Matsubara Green's function., Comment: 16 pages, 3 figures

The ab initio extension of the dynamical vertex approximation (D$\Gamma$A) method allows for realistic materials calculations that include non-local correlations beyond $GW$ and dynamical mean-field theory. Here, we discuss the AbinitioD$\Gamma$A algorithm, its implementation and usage in detail, and make the program package available to the scientific community.

X-ray absorption spectroscopy measured at the $L$-edge of transition metals (TMs) is a powerful element-selective tool providing direct information about the correlation effects in the $3d$ states. The theoretical modeling of the $2p\rightarrow3d$ excitation processes remains to be challenging for contemporary \textit{ab initio} electronic structure techniques, due to strong core-hole and multiplet effects influencing the spectra. In this work we present a realization of the method combining the density-functional theory with multiplet ligand field theory, proposed in Haverkort et al. (https://link.aps.org/doi/10.1103/PhysRevB.85.165113), Phys. Rev. B 85, 165113 (2012). In this approach a single-impurity Anderson model (SIAM) is constructed, with almost all parameters obtained from first principles, and then solved to obtain the spectra. In our implementation we adopt the language of the dynamical mean-field theory and utilize the local density of states and the hybridization function, projected onto TM $3d$ states, in order to construct the SIAM. The developed computational scheme is applied to calculate the $L$-edge spectra for several TM monoxides. A very good agreement between the theory and experiment is found for all studied systems. The effect of core-hole relaxation, hybridization discretization, possible extensions of the method as well as its limitations are discussed., Comment: 8 figures

Diagrammatic extensions of dynamical mean field theory (DMFT) such as the dynamical vertex approximation (D$\Gamma$A) allow us to include non-local correlations beyond DMFT on all length scales and proved their worth for model calculations. Here, we develop and implement an AbinitioD$\Gamma$A approach for electronic structure calculations of materials. Starting point is the two-particle irreducible vertex in the two particle-hole channels which is approximated by the bare non-local Coulomb interaction and all local vertex corrections. From this we calculate the full non-local vertex and the non-local self-energy through the Bethe-Salpeter equation. The AbinitioD$\Gamma$A approach naturally generates all local DMFT correlations and all non-local $GW$ contributions, but also further non-local correlations beyond: mixed terms of the former two and non-local spin fluctuations. We apply this new methodology to the prototypical correlated metal SrVO$_3$., Comment: 13 pages, 6 figures

We present an efficient implementation of the parquet formalism which respects the asymptotic structure of the vertex functions at both single- and two-particle levels in momentum- and frequency-space. We identify the two-particle reducible vertex as the core function which is essential for the construction of the other vertex functions. This observation stimulates us to consider a two-level parameter-reduction for this function to simplify the solution of the parquet equations. The resulting functions, which depend on fewer arguments, are coined "kernel functions". With the use of the "kernel functions", the open boundary of various vertex functions in the Matsubara-frequency space can be faithfully satisfied. We justify our implementation by accurately reproducing the dynamical mean-field theory results from momentum-independent parquet calculations. The high-frequency asymptotics of the single-particle self-energy and the two-particle vertex are correctly reproduced, which turns out to be essential for the self-consistent determination of the parquet solutions. The current implementation is also feasible for the dynamical vertex approximation., Comment: 12 pages, 11 figures

A proper theoretical description of electronic structure of the 3d orbitals in the metal centers of functional metalorganics is a challenging problem. In this letter, we apply density functional theory and an exact diagonalization method in a many body approach to study the ground state electronic configuration of an iron porphyrin (FeP) molecule. Our study reveals that dynamical correlation effects are important, and FeP is a potential candidate for realizing a spin crossover due to a subtle balance of crystal field effects, on-site Coulomb repulsion and hybridization between the Fe d-orbitals and ligand N p-states. The mechanism of switching between two close lying electronic configurations of Fe-d orbitals is shown. We discuss the generality of the suggested approach and the possibility to properly describe the electronic structure and related low energy physics of the whole class of correlated metal centered organometallic molecules.

We propose to calculate spectral functions of quantum impurity models using the Time Evolving Block Decimation (TEBD) for Matrix Product States. The resolution of the spectral function is improved by a so-called linear prediction approach. We apply the method as an impurity solver within the Dynamical Mean Field Theory (DMFT) for the single- and two-band Hubbard model on the Bethe lattice. For the single-band model we observe sharp features at the inner edges of the Hubbard bands. A finite size scaling shows that they remain present in the thermodynamic limit. We analyze the real time-dependence of the double occupation after adding a single electron and observe oscillations at the same energy as the sharp feature in the Hubbard band, indicating a long-lived coherent superposition of states that correspond to the Kondo peak and the side peaks. For a two-band Hubbard model we observe an even richer structure in the Hubbard bands, which cannot be related to a multiplet structure of the impurity, in addition to sharp excitations at the band edges of a type similar to the single-band case., Comment: 14 figures, 12 + pages including appendix. New Fig. 4b, Fig. 6, Fig.10, Fig.11 and Fig.A4

We improve a recently developed expansion technique for calculating real frequency spectral functions of any one-dimensional model with short-range interactions, by postprocessing computed Chebyshev moments with linear prediction. This can be achieved at virtually no cost and, in sharp contrast to existing methods based on the dampening of the moments, improves the spectral resolution rather than lowering it. We validate the method for the exactly solvable resonating level model and the single impurity Anderson model. It is capable of resolving sharp Kondo resonances, as well as peaks within the Hubbard bands when employed as an impurity solver for dynamical mean-field theory (DMFT). Our method works at zero temperature and allows for arbitrary discretization of the bath spectrum. It achieves similar precision as the dynamical density matrix renormalization group (DDMRG), at lower cost. We also propose an alternative expansion, of 1-exp(-tau H) instead of the usual H, which opens the possibility of using established methods for the time evolution of matrix product states to calculate spectral functions directly., Comment: 13 pages, 9 figures

Here we present a study of the entanglement in the electronic structure of the late transition metal monoxides - MnO, FeO, CoO, and NiO - obtained by means of density-functional theory in the local density approximation combined with dynamical mean-field theory (LDA+DMFT). The impurity problem is solved through Exact Diagonalization (ED), which grants full access to the thermally mixed many-body ground state density operator. The quality of the electronic structure is affirmed through a direct comparison between the calculated electronic excitation spectrum and photoemission experiments. Our treatment allows for a quantitative investigation of the entanglement in the electronic structure. Two main sources of entanglement are explicitly resolved through the use of a fidelity based geometrical entanglement measure, and additional information is gained from a complementary entropic entanglement measure. We show that the interplay of crystal field effects and Coulomb interaction causes the entanglement in CoO to take a particularly intricate form., Comment: Minor changes. Journal reference added

We generalize the adiabatic approximation to the case of open quantum systems, in the joint limit of slow change and weak open system disturbances. We show that the approximation is ``physically reasonable'' as under wide conditions it leads to a completely positive evolution, if the original master equation can be written on a time-dependent Lindblad form. We demonstrate the approximation for a non-Abelian holonomic implementation of the Hadamard gate, disturbed by a decoherence process. We compare the resulting approximate evolution with numerical simulations of the exact equation., Comment: New material added, references added and updated, journal reference added