*Nat Commun ; 12(1): 5744, 2021 Sep 30.*

##### RESUMO

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.

*Nat Mater ; 20(3): 329-334, 2021 Mar.*

##### RESUMO

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.

*J Phys Condens Matter ; 2020 Oct 05.*

##### RESUMO

An appropriate treatment of electronic correlation eï¬ects 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 eï¬cient, has no material speciï¬c 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.

*J Chem Theory Comput ; 16(10): 6256-6266, 2020 Oct 13.*

##### RESUMO

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.

*Chem Sci ; 10(26): 6508-6518, 2019 Jul 14.*

##### RESUMO

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.

*J Phys Condens Matter ; 31(33): 335601, 2019 Aug 21.*

##### RESUMO

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.

*J Phys Condens Matter ; 31(19): 195902, 2019 May 15.*

##### RESUMO

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.

*J Phys Condens Matter ; 30(19): 195801, 2018 May 16.*

##### RESUMO

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.

*J Am Chem Soc ; 140(5): 1674-1685, 2018 02 07.*

##### RESUMO

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.

*Nat Commun ; 8: 15565, 2017 06 01.*

##### RESUMO

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.

*Phys Rev Lett ; 118(12): 126401, 2017 Mar 24.*

##### RESUMO

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.

*Phys Rev Lett ; 115(9): 096401, 2015 Aug 28.*

##### RESUMO

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.

*J Phys Condens Matter ; 27(24): 245604, 2015 Jun 24.*

##### RESUMO

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.

*Phys Rev Lett ; 113(3): 036402, 2014 Jul 18.*

##### RESUMO

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.

*J Chem Phys ; 139(22): 224702, 2013 Dec 14.*

##### RESUMO

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.

*Phys Rev Lett ; 111(19): 196801, 2013 Nov 08.*

##### RESUMO

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.

*Nanoscale ; 4(8): 2567-70, 2012 Apr 21.*

##### RESUMO

Based on density-functional calculation and genetic algorithm structure search, we propose a series of 16-coordinate core-shell clusters: M@Li(16)(M = Ca, Sr, Ba, Ti, Zr, Hf). A tetrahedral (T(d)) structure with an outer shell of 16 lithium atoms and one enclosed heavy atom is found to be the global minimum in the structural exploration of BaLi(16) based on genetic algorithm. This structure also has lower energy compared to the other isomers we employed in all the MLi(16) clusters. In this structure, the atoms are bonded together by metallic bonds with alkali (IA) and alkaline-earth (IIA) metal atoms. Their corresponding first electronic shells are closed with significant energy gaps because their total numbers of valence electrons fulfil the 18-electron rule. Such a combination could be extended to 20-electron systems by enclosing IVB elements. With simple valence electrons and highly symmetric structures, superatomic molecular orbitals are identified in all of the T(d) clusters.