Search

The demand for disorder-tolerant quantum logic and spin electronics can be met by generating and controlling dissipationless spin currents protected by topology. Dirac fermions with helical spin-locking surface transport offer a way of achieving such a goal. Yet, surface-bulk coupling can lead to strong Dirac electron scattering with bulk carriers and phonons as well as impurities, assisted by such dissipative channel, which results in “topological breakdown”. Here, we demonstrate that coherent lattice vibrations periodically driven by a single-cycle terahertz (THz) pulse can significantly suppress such dissipative channel in topological insulators. This is achieved by reducing the phase space in the bulk available for Dirac fermion scattering into during coherent lattice oscillations in Bi2Se3. This light-induced suppression manifests as a remarkable transition exclusively in surface transport, absent for bulk, above the THz electric fields for driving coherent phonons, which prolongs the surface transport lifetime. These results, together with simulations, identify the critical role of spin–orbit coupling for the “phase space contraction” mechanism that suppresses the surface-bulk coupling. Imposing vibrational quantum coherence into topological states of matter may become a universal light control principle for reinforcing the symmetry-protected helical transport.

The choice that a solid system “makes” when adopting a crystal structure (stable or metastable) is ultimately governed by the interactions between electrons forming chemical bonds. Here we analyze six prototypical binary transition metal compounds and shed light on the connection between Mott physics and the behavior of the energy as a function of the spatial arrangement of the atoms in these systems. Remarkably, we find that the main qualitative features of this complex behavior in the Mott phase of these systems can be traced back to the fact that the strong d-electron correlations influence substantially the charge transfer mechanism, which, in turn, controls the electrostatic interactions. This result advances our understanding of the influence of strong correlations on the crystal structure, opens a new avenue for extending structure prediction methodologies to strongly correlated materials, and paves the way for predicting and studying metastability and polymorphism in these systems.

We analyze several methods of obtaining the accurate relativistic total energy (TE) variations using traditional perturbation theories (PTs) and proposed coupling constant integration (CCI) methods. For this purpose, we perform benchmark calculations within the density functional theory taking the spin-orbit coupling (SOC) and its derivative as a perturbation. The TE change due to SOC addition obtained from both PTs and CCI is shown to reach the accuracy of fully self-consistent TE calculations. Similar accuracy is also obtained even for the magnetocrystalline anisotropy energy (MAE). The real advantage of the proposed methods is to use PTs and CCI methods in those electronic structure methods where accurate total energies currently cannot be obtained with required accuracy. Correspondingly, we demonstrate the applicability of suggested methods for calculations of MAE in different magnetic materials using a dynamic mean-field method. All suggested PTs and CCI methods also provide convenient site, orbital, and spin decompositions of the TE variation, creating a powerful way to analyze microscopic physics in strongly correlated materials.

We propose a scheme of decomposition of the total relativistic energy in solids to intra- and interatomic contributions. The method is based on a site variation of such fundamental constant as the speed of light. As a practical illustration of the method, we tested such decomposition in the case of a spin-orbit interaction variation for the decomposition of the magnetic anisotropy energy (MAE) in CoPt. We further studied the ${\ensuremath{\alpha}}^{\ensuremath{''}}\text{\ensuremath{-}}{\mathrm{Fe}}_{16}{\mathrm{N}}_{2}$ magnet doped by Bi, Sb, Co, and Pt atoms. It was found that the addition of Pt atoms can enhance the MAE by as much as five times while Bi and Sb substitutions double the total MAE. Using the proposed technique, we demonstrate the spatial distribution of these enhancements. Our studies also suggest that Sb, Pt, and Co substitutions could be synthesized by experiments.

An appropriate treatment of electronic correlation effects plays an important role in accurate descriptions of physical and chemical properties of real materials. The recently proposed Correlation Matrix Renormalization theory with Sum Rule correction (CMR) for studying correlated electron materails has shown good performance in molecular systems and a periodic Hydrogen chain in comparison with various quantum chemistry and quantum Monte Carlo calculations. This work gives a detailed formulation and computational code implementation of CMR in multi-band periodic lattice systems. This lattice CMR ab initio theory is highly efficient, has no material specific adjustable parameters, and has no double counting issues faced by the hybrid approaches like LDA+U, DFT+DMFT and DFT+GA type theories. Benchmark studies on materials with s and p orbitals in this study show that CMR in its current implementation consistently performs well for these systems as the electron correlation increases from the bonding region to the bond breaking region.

A cardinal obstacle to performing quantum-mechanical simulations of strongly-correlated matter is that, with the theoretical tools presently available, sufficiently-accurate computations are often too expensive to be ever feasible. Here we design a computational framework combining quantum-embedding (QE) methods with machine learning. This allows us to bypass altogether the most computationally-expensive components of QE algorithms, making their overall cost comparable to bare Density Functional Theory (DFT). We perform benchmark calculations of a series of actinide systems, where our method describes accurately the correlation effects, reducing by orders of magnitude the computational cost. We argue that, by producing a larger-scale set of training data, it will be possible to apply our method to systems with arbitrary stoichiometries and crystal structures, paving the way to virtually infinite applications in condensed matter physics, chemistry and materials science., Main text (6 pages, 4 figures ) and supplemental material (5 pages)

We introduce Gutzwiller conjugate gradient minimization (GCGM) theory, an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. GCGM uses the Gutzwiller wave function but does not use the commonly adopted Gutzwiller approximation (GA), which is a major source of inaccuracy. Instead, the theory uses an approximation that is based on the occupation probability of the on-site configurations, rather than approximations that decouple the site-site correlations as used in the GA. We test the theory in the one-dimensional and two-dimensional Hubbard models at various electron densities and find that GCGM reproduces energies and double occupancies in reasonable agreement with benchmark data at a very small computational cost.

A grand challenge underlies the entire field of topology-enabled quantum logic and information science: how to establish topological control principles driven by quantum coherence and understand the time dependence of such periodic driving. Here we demonstrate a few-cycle THz-pulse-induced phase transition in a Dirac semimetal ZrTe_{5} that is periodically driven by vibrational coherence due to excitation of the lowest Raman active mode. Above a critical THz-pump field threshold, there emerges a long-lived metastable phase, approximately 100 ps, with unique Raman phonon-assisted topological switching dynamics absent for optical pumping. The switching also manifests itself by distinct features: nonthermal spectral shape, relaxation slowing near the Lifshitz transition where the critical Dirac point occurs, and diminishing signals at the same temperature that the Berry-curvature-induced anomalous Hall effect magnetoresistance vanishes. These results, together with first-principles modeling, identify a mode-selective Raman coupling that drives the system from strong to weak topological insulators with a Dirac semimetal phase established at a critical atomic displacement controlled by the phonon coherent pumping. Harnessing of vibrational coherence can be extended to steer symmetry-breaking transitions, i.e., Dirac to Weyl ones, with implications for THz topological quantum gate and error correction applications.

We discover hidden Rashba fine structure in ${\mathrm{CH}}_{3}{\mathrm{NH}}_{3}{\mathrm{PbI}}_{3}$ and demonstrate its quantum control by vibrational coherence through symmetry-selective vibronic (electron-phonon) coupling. Above a critical threshold of a single-cycle terahertz pump field, a Raman phonon mode distinctly modulates the middle excitonic states with persistent coherence for more than ten times longer than the ones on two sides that predominately couple to infrared phonons. These vibronic quantum beats, together with first-principles modeling of phonon periodically modulated Rashba parameters, identify a threefold excitonic fine structure splitting, i.e., optically forbidden, degenerate dark states in between two bright ones with a narrow, $\ensuremath{\sim}3\text{ }\text{ }\text{nm}$ splitting. Harnessing of vibronic quantum coherence and symmetry inspires light-perovskite quantum control and sub-THz-cycle ``Rashba engineering'' of spin-split bands for ultimate multifunction device.

The coherence of collective modes, such as phonons and polarons, and their modulation of electronic states is long sought in complex systems, which is a crosscutting issue in photovoltaics and quantum electronics. In photovoltaic cells and lasers based on metal halide perovskites, the presence of polarons, i.e., photocarriers dressed by the macroscopic motion of charged lattice, assisted by terahertz (THz) longitudinal-optical (LO) phonons, has been intensely studied yet is still debated. This may be key for explaining the remarkable properties of the perovskite materials, e.g., defect tolerance, long charge lifetimes, and diffusion lengths. Here we use the intense single-cycle THz pulse with peak electric field up to ${E}_{\mathrm{THz}}=1000$ kV/cm to drive coherent polaronic band-edge oscillations at room temperature in ${\mathrm{CH}}_{3}{\mathrm{NH}}_{3}{\mathrm{PbI}}_{3}\phantom{\rule{4pt}{0ex}}({\mathrm{MAPbI}}_{3})$. We reveal the oscillatory behavior is dominated by a specific quantized lattice vibration mode at ${\ensuremath{\omega}}_{\mathrm{LO}}\ensuremath{\sim}4\phantom{\rule{0.16em}{0ex}}\mathrm{THz}$, which is both dipole and momentum forbidden. THz-driven coherent polaron dynamics exhibits distinguishing features: room temperature coherent oscillations at ${\ensuremath{\omega}}_{\mathrm{LO}}$ longer than 1 ps in both single crystals and thin films, mode-selective modulation of different band-edge states assisted by electron-phonon interaction, and dynamic mode splitting at low temperature due to entropy and anharmonicity of organic cations. Our results demonstrate intense THz-driven coherent band-edge modulation is a powerful probe of electron-lattice coupling phenomena and polaronic quantum control in perovskites.

We recently proposed the correlation matrix renormalisation (CMR) theory to efficiently and accurately calculate ground state total energy of molecular systems, based on the Gutzwiller variational wavefunction (GWF) to treat the electronic correlation effects. To help reduce numerical complications and better adapt the CMR to infinite lattice systems, we need to further refine the way to minimise the error originated from the approximations in the theory. This conference proceeding reports our recent progress on this key issue, namely, we obtained a simple analytical functional form for the one-electron renormalisation factors, and introduced a novel sum-rule correction for a more accurate description of the intersite electron correlations. Benchmark calculations are performed on a set of molecules to show the reasonable accuracy of the method.

We propose an approach that is under the framework of Gutzwiller wave function but goes beyond the commonly adopted Gutzwiller approximation to improve the accuracy and flexibility in treating the correlation effects. Detailed formalism is described for a dimer which is straightforwardly generalized later to more complicated periodic bulk systems. The accuracy of the approach is demonstrated by evaluating the potential energy curves of spin-singlet N2 dimer, spin-triplet O2 dimer, and 1D hydrogen chain. The computational workload of the approach can be easily handled by efficient parallel computing.

We present detailed benchmark ground-state calculations of the one- and two-dimensional Hubbard model utilizing the cluster extensions of the rotationally invariant slave-boson mean-field theory and the density matrix embedding theory. Our analysis shows that the overall accuracy and the performance of these two methods are very similar. Furthermore, we propose a unified computational framework that allows us to implement both of these techniques on the same footing. This provides us with a different line of interpretation and paves the ways for developing systematically distinct generalizations of these complementary approaches.

We report terahertz (THz) electrodynamics of a moderately clean A15 superconductor (SC) following ultrafast excitation to manipulate quasiparticle (QP) transport. In the Martensitic normal state, we observe a photo enhancement in the THz conductivity using optical pulses, while the opposite is observed for the THz pump. This demonstrates wavelength-selective nonthermal control of conductivity distinct from sample heating. The photo enhancement persists up to an additional critical temperature, above the SC one, from a competing electronic order. In the SC state, the fluence dependence of pair-breaking kinetics together with an analytic model provides an implication for a ``one photon to one Cooper pair'' nonresonant energy transfer during the 35-fs laser pulse; i.e., the fitted photon energy $\ensuremath{\hbar}{\ensuremath{\omega}}_{}$ absorption to create QPs set by $2{\mathrm{\ensuremath{\Delta}}}_{SC}/\ensuremath{\hbar}\ensuremath{\omega}=0.33%$. This is more than one order of magnitude smaller than in previously studied BCS SCs, which we attribute to strong electron-phonon coupling and possible influence of phonon condensation.

We recently proposed the correlation matrix renormalization (CMR) theory to treat the electronic correlation effects [Phys. Rev. B 2014, 89, 045131 and Sci. Rep. 2015, 5, 13478] in ground state total energy calculations of molecular systems using the Gutzwiller variational wave function (GWF). By adopting a number of approximations, the computational effort of the CMR can be reduced to a level similar to Hartree-Fock calculations. This paper reports our recent progress in minimizing the error originating from some of these approximations. We introduce a novel sum-rule correction to obtain a more accurate description of the intersite electron correlation effects in total energy calculations. Benchmark calculations are performed on a set of molecules to show the reasonable accuracy of the method.

The dependence of the magnetocrystalline anisotropy energy (MAE) in MCo5 (M = Y, La, Ce, Gd) and CoPt on the Coulomb correlations and strength of spin orbit (SO) interaction within the GGA + U scheme is investigated. A range of parameters suitable for the satisfactory description of key magnetic properties is determined. We show that for a large variation of SO interaction the MAE in these materials can be well described by the traditional second order perturbation theory. We also show that in these materials the MAE can be both proportional and negatively proportional to the orbital moment anisotropy (OMA) of Co atoms. Dependence of relativistic effects on Coulomb correlations, applicability of the second order perturbation theory for the description of MAE, and effective screening of the SO interaction in these systems are discussed using a generalized virial theorem. Such determined sets of parameters of Coulomb correlations can be used in much needed large scale atomistic simulations.

Developing accurate and computationally efficient methods to calculate the electronic structure and total energy of correlated-electron materials has been a very challenging task in condensed matter physics and materials science. Recently, we have developed a correlation matrix renormalization (CMR) method which does not assume any empirical Coulomb interaction $U$ parameters and does not have double counting problems in the ground-state total energy calculation. The CMR method has been demonstrated to be accurate in describing both the bonding and bond breaking behaviors of molecules. In this study, we extend the CMR method to the treatment of electron correlations in periodic solid systems. Using a linear hydrogen chain as a benchmark system, we show that the results from the CMR method compare very well with those obtained recently by accurate quantum Monte Carlo (QMC) calculations. We also study the equation of states of three-dimensional crystalline phases of atomic hydrogen. We show that the results from the CMR method agree much better with the available QMC data in comparison with those from density functional theory and Hartree-Fock calculations.

We develop a unified theoretical picture for excitations in Mott systems, portraying both the heavy quasiparticle excitations and the Hubbard bands as features of an emergent Fermi liquid state formed in an extended Hilbert space, which is non-perturbatively connected to the physical system. This observation sheds light on the fact that even the incoherent excitations in strongly correlated matter often display a well defined Bloch character, with pronounced momentum dispersion. Furthermore, it indicates that the Mott point can be viewed as a topological transition, where the number of distinct dispersing bands displays a sudden change at the critical point. Our results, obtained from an appropriate variational principle, display also remarkable quantitative accuracy. This opens an exciting avenue for fast realistic modeling of strongly correlated materials., 5 pages, 3 figures

We investigate the normal state of the superconducting compound ${\mathrm{PuCoGa}}_{5}$ using the combination of density functional theory (DFT) and dynamical mean-field theory (DMFT), with the continuous time quantum Monte Carlo (CTQMC) and the vertex-corrected one-crossing approximation (OCA) as the impurity solvers. Our DFT+DMFT (CTQMC) calculations suggest a strong tendency of Pu-$5f$ orbitals to differentiate at low temperatures. The renormalized $5{f}_{5/2}$ states exhibit a Fermi-liquid behavior whereas one electron in the $5{f}_{7/2}$ states is at the edge of a Mott localization. We find that the orbital differentiation is manifested as the removing of $5{f}_{7/2}$ spectral weight from the Fermi level relative to DFT. We corroborate these conclusions with DFT+DMFT (OCA) calculations which demonstrate that $5{f}_{5/2}$ electrons have a much larger Kondo scale than the $5{f}_{7/2}$.

How photoexcitations evolve into Coulomb-bound electron and hole pairs, called excitons, and unbound charge carriers is a key cross-cutting issue in photovoltaics and optoelectronics. Until now, the initial quantum dynamics following photoexcitation remains elusive in the hybrid perovskite system. Here we reveal excitonic Rydberg states with distinct formation pathways by observing the multiple resonant, internal quantum transitions using ultrafast terahertz quasi-particle transport. Nonequilibrium emergent states evolve with a complex co-existence of excitons, carriers and phonons, where a delayed buildup of excitons under on- and off-resonant pumping conditions allows us to distinguish between the loss of electronic coherence and hot state cooling processes. The nearly ∼1 ps dephasing time, efficient electron scattering with discrete terahertz phonons and intermediate binding energy of ∼13.5 meV in perovskites are distinct from conventional photovoltaic semiconductors. In addition to providing implications for coherent energy conversion, these are potentially relevant to the development of light-harvesting and electron-transport devices., The generation of bound electron and hole pairs—excitons—is a key process in photovoltaic technologies, yet it is challenging to follow their initial dynamics. Here, Luo et al. probe the Rydberg eigenstates that characterize the excitonic transport and coherent conversion in a perovskite material.

We derive an exact operatorial reformulation of the rotational invariant slave boson method and we apply it to describe the orbital differentiation in strongly correlated electron systems starting from first principles. The approach enables us to treat strong electron correlations, spin-orbit coupling and crystal field splittings on the same footing by exploiting the gauge invariance of the mean-field equations. We apply our theory to the archetypical nuclear fuel UO$_2$, and show that the ground state of this system displays a pronounced orbital differention within the $5f$ manifold, with Mott localized $\Gamma_8$ and extended $\Gamma_7$ electrons., Comment: Accepted for publication in Physical Review Letters (5 pages, 1 figure)

We develop a variational scheme called "Gutzwiller renormalization group" (GRG), which enables us to calculate the ground state of Anderson impurity models (AIM) with arbitrary numerical precision. Our method can exploit the low-entanglement property of the ground state in combination with the framework of the Gutzwiller wavefunction, and suggests that the ground state of the AIM has a very simple structure, which can be represented very accurately in terms of a surprisingly small number of variational parameters. We perform benchmark calculations of the single-band AIM that validate our theory and indicate that the GRG might enable us to study complex systems beyond the reach of the other methods presently available and pave the way to interesting generalizations, e.g., to nonequilibrium transport in nanostructures., Comment: 5 pages, 3 figures

Using resistivity, heat-capacity, thermal-expansion, and susceptibility measurements we study the normal-state behavior of ${\mathrm{KFe}}_{2}{\mathrm{As}}_{2}$. Both the Sommerfeld coefficient ($\ensuremath{\gamma}\ensuremath{\approx}103\text{ }\text{ }\mathrm{mJ}\text{ }{\mathrm{mol}}^{\ensuremath{-}1}\text{ }{\mathrm{K}}^{\ensuremath{-}2}$) and the Pauli susceptibility ($\ensuremath{\chi}\ensuremath{\approx}4\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}4}$) are strongly enhanced, which confirm the existence of heavy quasiparticles inferred from previous de Haas--van Alphen and angle-resolved photoemission spectroscopy experiments. We discuss this large enhancement using a Gutzwiller slave-boson mean-field calculation, which shows the proximity of ${\mathrm{KFe}}_{2}{\mathrm{As}}_{2}$ to an orbital-selective Mott transition. The temperature dependence of the magnetic susceptibility and the thermal expansion provide strong experimental evidence for the existence of a coherence-incoherence crossover, similar to what is found in heavy fermion and ruthenate compounds, due to Hund's coupling between orbitals.