Chirality manipulation of ultrafast phase switches in a correlated CDW-Weyl semimetal.

Light engineering of correlated states in topological materials provides a new avenue of achieving exotic topological phases inaccessible by conventional tuning methods. Here we demonstrate a light control of correlation gaps in a model charge-density-wave (CDW) and polaron insulator (TaSe4)2I recently predicted to be an axion insulator. Our ultrafast terahertz photocurrent spectroscopy reveals a two-step, non-thermal melting of polarons and electronic CDW gap via the fluence dependence of a longitudinal circular photogalvanic current. This helicity-dependent photocurrent reveals continuous ultrafast phase switches from the polaronic state to the CDW (axion) phase, and finally to a hidden Weyl phase as the pump fluence increases. Additional distinctive attributes aligning with the light-induced switches include: the mode-selective coupling of coherent phonons to the polaron and CDW modulation, and the emergence of a non-thermal chiral photocurrent above the pump threshold of CDW-related phonons. The demonstrated ultrafast chirality control of correlated topological states here holds large potentials for realizing axion electrodynamics and advancing quantum-computing applications.

Ground and excited states of even-numbered Hubbard ring at half-filling: comparison of the extended Gutzwiller approach with exact diagonalization.

It remains a great challenge in condensed matter physics to develop a method to treat strongly correlated many-body systems with balanced accuracy and efficiency. We introduce an extended Gutzwiller (EG) method incorporating a manifold technique, which builds an effective manifold of the many-body Hilbert space, to describe the ground-state (GS) and excited-state (ES) properties of strongly correlated electrons. We systematically apply an EG projector onto the GS and ES of a non-interacting system. Diagonalization of the true Hamiltonian within the manifold formed by the resulting EG wavefunctions gives the approximate GS and ES of the correlated system. To validate this technique, we implement it on even-numbered fermionic Hubbard rings at half-filling with periodic boundary conditions, and compare the results with the exact diagonalization (ED) method. The EG method is capable of generating high-quality GS and low-lying ES wavefunctions, as evidenced by the high overlaps of wavefunctions between the EG and ED methods. Favorable comparisons are also achieved for other quantities including the total energy, the double occupancy, the total spin and the staggered magnetization. With the capability of accessing the ESs, the EG method can capture the essential features of the one-electron removal spectral function that contains contributions from states deep in the excited spectrum. Finally, we provide an outlook on the application of this method on large extended systems.

A rotationally invariant approach based on Gutzwiller wave function for correlated electron systems.

We introduce a rotationally invariant approach combined with the Gutzwiller conjugate gradient minimization method to study correlated electron systems. In the approach, the Gutzwiller projector is parametrized based on the number of electrons occupying the onsite orbitals instead of the onsite configurations. The approach efficiently groups the onsite orbitals according to their symmetry and greatly reduces the computational complexity, which yields a speedup of20∼50×in the minimal basis energy calculation of dimers. The computationally efficient approach promotes more accurate calculations beyond the minimal basis that is inapplicable in the original approach. A large-basis energy calculation of F2demonstrates favorable agreements with standard quantum-chemical calculations Bytautaset al(2007J. Chem. Phys.127164317).

The Gutzwiller conjugate gradient minimization method for correlated electron systems.

We review our recent work on the Gutzwiller conjugate gradient minimization method, anab initioapproach developed for correlated electron systems. The complete formalism has been outlined that allows for a systematic understanding of the method, followed by a discussion of benchmark studies of dimers, one- and two-dimensional single-band Hubbard models. In the end, we present some preliminary results of multi-band Hubbard models and large-basis calculations of F2to illustrate our efforts to further reduce the computational complexity.

Coherent control of asymmetric spintronic terahertz emission from two-dimensional hybrid metal halides.

Next-generation terahertz (THz) sources demand lightweight, low-cost, defect-tolerant, and robust components with synergistic, tunable capabilities. However, a paucity of materials systems simultaneously possessing these desirable attributes and functionalities has made device realization difficult. Here we report the observation of asymmetric spintronic-THz radiation in Two-Dimensional Hybrid Metal Halides (2D-HMH) interfaced with a ferromagnetic metal, produced by ultrafast spin current under femtosecond laser excitation. The generated THz radiation exhibits an asymmetric intensity toward forward and backward emission direction whose directionality can be mutually controlled by the direction of applied magnetic field and linear polarization of the laser pulse. Our work demonstrates the capability for the coherent control of THz emission from 2D-HMHs, enabling their promising applications on the ultrafast timescale as solution-processed material candidates for future THz emitters.

A light-induced phononic symmetry switch and giant dissipationless topological photocurrent in ZrTe5.

Dissipationless currents from topologically protected states are promising for disorder-tolerant electronics and quantum computation. Here, we photogenerate giant anisotropic terahertz nonlinear currents with vanishing scattering, driven by laser-induced coherent phonons of broken inversion symmetry in a centrosymmetric Dirac material ZrTe5. Our work suggests that this phononic terahertz symmetry switching leads to formation of Weyl points, whose chirality manifests in a transverse, helicity-dependent current, orthogonal to the dynamical inversion symmetry breaking axis, via circular photogalvanic effect. The temperature-dependent topological photocurrent exhibits several distinct features: Berry curvature dominance, particle-hole reversal near conical points and chirality protection that is responsible for an exceptional ballistic transport length of ~10 µm. These results, together with first-principles modelling, indicate two pairs of Weyl points dynamically created by B1u phonons of broken inversion symmetry. Such phononic terahertz control breaks ground for coherent manipulation of Weyl nodes and robust quantum transport without application of static electric or magnetic fields.

Correlation matrix renormalization theory in multi-band lattice systems.

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) [1] for studying correlated-electron materials has shown good performance in molecular systems and a periodic hydrogen chain in comparison with various quantum chemistry and quantum Monte Carlo calculations [2]. This work gives a detailed formulation and computational code implementation of CMR in multi-band periodic lattice systems. This lattice CMRab initiotheory 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.

Efficient Step-Merged Quantum Imaginary Time Evolution Algorithm for Quantum Chemistry.

We develop a resource-efficient step-merged quantum imaginary time evolution approach (smQITE) to solve for the ground state of a Hamiltonian on quantum computers. This heuristic method features a fixed shallow quantum circuit depth along the state evolution path. We use this algorithm to determine the binding energy curves of a set of molecules, including H2, H4, H6, LiH, HF, H2O, and BeH2, and find highly accurate results. The required quantum resources of smQITE calculations can be further reduced by adopting the circuit form of the variational quantum eigensolver (VQE) technique, such as the unitary coupled cluster ansatz. We demonstrate that smQITE achieves a similar computational accuracy as VQE at the same fixed-circuit ansatz, without requiring a generally complicated high-dimensional nonconvex optimization. Finally, smQITE calculations are carried out on Rigetti quantum processing units, demonstrating that the approach is readily applicable on current noisy intermediate-scale quantum devices.

Origins of the odd optical observables in plutonium and americium tungstates.

A series of trivalent f-block tungstates, MW2O7(OH)(H2O) (M = La, Ce, Pr, Nd, and Pu) and AmWO4(OH), have been prepared in crystalline form using hydrothermal methods. Both structure types take the form of 3D networks where MW2O7(OH)(H2O) is assembled from infinite chains of distorted tungstate octahedra linked by isolated MO8 bicapped trigonal prisms; whereas AmWO4(OH) is constructed from edge-sharing AmO8 square antiprisms connected by distorted tungstate trigonal bipyramids. PuW2O7(OH)(H2O) crystallizes as red plates; an atypical color for a Pu(iii) compound. Optical absorption spectra acquired from single crystals show strong, broadband absorption in the visible region. A similar feature is observed for CeW2O7(OH)(H2O), but not for AmWO4(OH). Here we demonstrate that these significantly different optical properties do not stem directly from the 5f electrons, as in both systems the valence band has mostly O-2p character and the conduction band has mostly W-5d character. Furthermore, the quasi-particle gap is essentially unaffected by the 5f degrees of freedom. Despite this, our analysis demonstrates that the f-electron covalency effects are quite important and substantially different energetically in PuW2O7(OH)(H2O) and AmWO4(OH), indicating that the optical gap alone cannot be used to infer conclusions concerning the f electron contribution to the chemical bond in these systems.

First-principles calculation of correlated electron materials based on Gutzwiller wave function beyond Gutzwiller approximation.

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.

Benchmark of correlation matrix renormalization method in molecule calculations.

We report benchmark calculations of the correlation matrix renormalization (CMR) approach for 23 molecules in the well-established G2 molecule set. This subset represents molecules with spin-singlet ground state in a variety of chemical bonding and coordination environments. The QUAsi-atomic minimal basis-set orbitals (QUAMBOs) are used as local orbitals in both CMR and full configuration interaction (FCI) calculations for comparison. The results obtained from the calculations are also compared with available experimental data. It is shown that the CMR method produces binding and dissociation energy curves in good agreement with the QUAMBO-FCI calculations as well as experimental results. The CMR benchmark calculations yield a standard deviation of 0.09 Å for the equilibrium bond length and 0.018 Hartree/atom for the formation energy, with a gain of great computational efficiency which scales like Hartree-Fock method.

Magnetocrystalline anisotropy in cobalt based magnets: a choice of correlation parameters and the relativistic effects.

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.

Uncovering the Origin of Divergence in the CsM(CrO4)2 (M = La, Pr, Nd, Sm, Eu; Am) Family through Examination of the Chemical Bonding in a Molecular Cluster and by Band Structure Analysis.

A series of f-block chromates, CsM(CrO4)2 (M = La, Pr, Nd, Sm, Eu; Am), were prepared revealing notable differences between the AmIII derivatives and their lanthanide analogs. While all compounds form similar layered structures, the americium compound exhibits polymorphism and adopts both a structure isomorphous with the early lanthanides as well as one that possesses lower symmetry. Both polymorphs are dark red and possess band gaps that are smaller than the LnIII compounds. In order to probe the origin of these differences, the electronic structure of α-CsSm(CrO4)2, α-CsEu(CrO4)2, and α-CsAm(CrO4)2 were studied using both a molecular cluster approach featuring hybrid density functional theory and QTAIM analysis and by the periodic LDA+GA and LDA+DMFT methods. Notably, the covalent contributions to bonding by the f orbitals were found to be more than twice as large in the AmIII chromate than in the SmIII and EuIII compounds, and even larger in magnitude than the Am-5f spin-orbit splitting in this system. Our analysis indicates also that the Am-O covalency in α-CsAm(CrO4)2 is driven by the degeneracy of the 5f and 2p orbitals, and not by orbital overlap.

Ultrafast terahertz snapshots of excitonic Rydberg states and electronic coherence in an organometal halide perovskite.

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.

Slave Boson Theory of Orbital Differentiation with Crystal Field Effects: Application to UO_{2}.

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 differentiation within the 5f manifold, with Mott-localized Γ_{8} and extended Γ_{7} electrons.

Nature of the Insulating Ground State of the 5d Postperovskite CaIrO3.

The insulating ground state of the 5d transition metal oxide CaIrO3 has been classified as a Mott-type insulator. Based on a systematic density functional theory (DFT) study with local, semilocal, and hybrid exchange-correlation functionals, we reveal that the Ir t(2g) states exhibit large splittings and one-dimensional electronic states along the c axis due to a tetragonal crystal field. Our hybrid DFT calculation adequately describes the antiferromagnetic (AFM) order along the c direction via a superexchange interaction between Ir^{4+} spins. Furthermore, the spin-orbit coupling (SOC) hybridizes the t(2g) states to open an insulating gap. These results indicate that CaIrO_{3} can be represented as a spin-orbit Slater insulator, driven by the interplay between a long-range AFM order and the SOC. Such a Slater mechanism for the gap formation is also demonstrated by the DFT + dynamical mean field theory calculation, where the metal-insulator transition and the paramagnetic to AFM phase transition are concomitant with each other.

A coefficient average approximation towards Gutzwiller wavefunction formalism.

Gutzwiller wavefunction is a physically well-motivated trial wavefunction for describing correlated electron systems. In this work, a new approximation is introduced to facilitate the evaluation of the expectation value of any operator within the Gutzwiller wavefunction formalism. The basic idea is to make use of a specially designed average over Gutzwiller wavefunction coefficients expanded in the many-body Fock space to approximate the ratio of expectation values between a Gutzwiller wavefunction and its underlying noninteracting wavefunction. To check with the standard Gutzwiller approximation (GA), we test its performance on single band systems and find quite interesting properties. On finite systems, we noticed that it gives superior performance over GA, while on infinite systems it asymptotically approaches GA. Analytic analysis together with numerical tests are provided to support this claimed asymptotical behavior. Finally, possible improvements on the approximation and its generalization towards multiband systems are illustrated and discussed.

Principle of maximum entanglement entropy and local physics of strongly correlated materials.

We argue that, because of quantum entanglement, the local physics of strongly correlated materials at zero temperature is described in a very good approximation by a simple generalized Gibbs distribution, which depends on a relatively small number of local quantum thermodynamical potentials. We demonstrate that our statement is exact in certain limits and present numerical calculations of the iron compounds FeSe and FeTe and of the elemental cerium by employing the Gutzwiller approximation that strongly support our theory in general.

Strain effect on the adsorption, diffusion, and molecular dissociation of hydrogen on Mg (0001) surface.

The adsorption, diffusion, and molecular dissociation of hydrogen on the biaxially strained Mg (0001) surface have been systematically investigated by the first principle calculations based on density functional theory. When the strain changes from the compressive to tensile state, the adsorption energy of H atom linearly increases while its diffusion barrier linearly decreases oppositely. The dissociation barrier of H2 molecule linearly reduces in the tensile strain region. Through the chemical bonding analysis including the charge density difference, the projected density of states and the Mulliken population, the mechanism of the strain effect on the adsorption of H atom and the dissociation of H2 molecule has been elucidated by an s-p charge transfer model. With the reduction of the orbital overlap between the surface Mg atoms upon the lattice expansion, the charge transfers from p to s states of Mg atoms, which enhances the hybridization of H s and Mg s orbitals. Therefore, the bonding interaction of H with Mg surface is strengthened and then the atomic diffusion and molecular dissociation barriers of hydrogen decrease accordingly. Our works will be helpful to understand and to estimate the influence of the lattice deformation on the performance of Mg-containing hydrogen storage materials.

γ-α isostructural transition in cerium.

We present zero-temperature first-principles calculations of elemental cerium and we compute its pressure-volume phase diagram within a theoretical framework able to describe simultaneously both the α and the γ phases. A surprising result revealed by our study is the presence of a clear signature of the transition at zero temperature and that this signature can be observed if and only if the spin-orbit coupling is taken into account. Our calculations indicate that the transition line in the pressure-temperature phase diagram of this material has a low-T critical point at negative pressures, placed very close to zero temperature. This suggests that cerium is very close to being "quantum critical," in agreement with recent experiments.