Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density Matrix”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Block-Type Antiferromagnetism in Single Chain Quasi-One-Dimensional K 3 ⁢Fe 2 ⁢Se 4

One-dimensional (1D) structures provide a unique platform to study the correlated quantum interactions and phase transitions such as unconventional magnetism and superconducting states. Here, we report that iron chalcogenide K 3 ⁢Fe 2 ⁢Se 4 exhibits an unusual block-type canted antiferromagnetic (AFM) order with a clear single chain quasi-1D structure, which is structurally different from the two-leg ladder BaFe 2 ⁢Se 3 , through both experimental measurements and density matrix renormalization group (DMRG) calculations. The narrow bandgap semiconductor K 3 ⁢Fe 2 ⁢Se 4 has a quasi-1D edge-shared FeSe4 tetrahedra chain structure and orders antiferromagnetically below 110 K. The magnetic moments couple antiferromagnetically along the quasi-1D chain direction of the 𝑏 axis and form an up-down-down-up (↑−↓−↓−↑)–like spin structure with a commensurate propagation vector 𝒌=⁢(0,0,0), where block-type spin ↑−↑ or ↓−↓ coupling are between the longer Fe-Fe bonds of the quasi-1D chain. DMRG results show that block antiferromagnetic state is stable in K 3 ⁢Fe 2 ⁢Se 4 and reveal that the block-ordered arrangement of Fe 2.5+ ions spins arise from the competition between ferromagnetic and AFM interaction in the presence of strong electronic correlation. Our research results not only report the discovery of a clear block-type canted antiferromagnetic structure in a real quasi-1D chain material but also provide a theoretical approach to understand the block-type antiferromagnetism in quasi-1D iron chalcogenides.

antiferromagnetism↗
DOE OSTI2021

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Theory of optical transitions in π-conjugated macrocycles

We describe a theoretical and computational investigation of the optical properties of π-conjugated macrocycles. Since the low-energy excitations of these systems are Frenkel excitons that couple to high-frequency dispersionless phonons, we employ the quantized Frenkel-Holstein model and solve it via the density matrix renormalization group (DMRG) method. First we consider optical emission from perfectly circular systems. Owing to optical selection rules, such systems radiate via two mechanisms: (i) within the Condon approximation, by thermally induced emission from the optically allowed j = ± 1 states and (ii) beyond the Condon approximation, by emission from the j = 0 state via coupling with a totally non-symmetric phonon (namely, the Herzberg-Teller effect). Using perturbation theory, we derive an expression for the Herzberg-Teller correction and show via DMRG calculations that this expression soon fails as ħ ω/J and the size of the macrocycle increase. Next, we consider the role of broken symmetry caused by torsional disorder. In this case the quantum number j no longer labels eigenstates of angular momentum, but instead labels localized local exciton groundstates (LEGSs) or quasi-extended states (QEESs). As for linear polymers, LEGSs define chromophores, with the higher energy QEESs being extended over numerous LEGSs. Within the Condon approximation (i.e., neglecting the Herzberg-Teller correction) we show that increased disorder increases the emissive optical intensity, because all the LEGSs are optically active. We next consider the combined role of broken symmetry and curvature, by explicitly evaluating the Herzberg-Teller correction in disordered systems via the DMRG method. The Herzberg-Teller correction is most evident in the emission intensity ratio, I{sub 00}/I{sub 01}. In the Condon approximation I{sub 00}/I{sub 01} is a constant function of curvature, whereas in practice it vanishes for closed rings and only approaches a constant in the limit of vanishing curvature. We calculate the optical spectra of a model system, cyclo-poly(para-phenylene ethynylene), for different amounts of torsional disorder within and beyond the Condon approximation. We show how broken symmetry and the Herzberg-Teller effect explain the spectral features. The Herzberg-Teller correction to the 0–1 emission vibronic peak is always significant. Finally, we note the qualitative similarities between the optical properties of conformationally disordered linear polymers and macrocycles in the limit of sufficiently large disorder, because in both cases they are determined by the optical properties of curved chromophores.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Golden rule kinetics of transfer reactions in condensed phase: The microscopic model of electron transfer reactions in disordered solid matrices

The algorithm for a theoretical calculation of transfer reaction rates for light quantum particles (i.e., the electron and H-atom transfers) in non-polar solid matrices is formulated and justified. The mechanism postulated involves a local mode (an either intra- or inter-molecular one) serving as a mediator which accomplishes the energy exchange between the reacting high-frequency quantum mode and the phonon modes belonging to the environment. This approach uses as a background the Fermi golden rule beyond the usually applied spin-boson approximation. The dynamical treatment rests on the one-dimensional version of the standard quantum relaxation equation for the reduced density matrix, which describes the frequency fluctuation spectrum for the local mode under consideration. The temperature dependence of a reaction rate is controlled by the dimensionless parameter ξ{sub 0}=ℏω{sub 0}/k{sub B}T where ω{sub 0} is the frequency of the local mode and T is the temperature. The realization of the computational scheme is different for the high/intermediate (ξ{sub 0} < 1 − 3) and for low (ξ{sub 0}≫ 1) temperature ranges. For the first (quasi-classical) kinetic regime, the Redfield approximation to the solution of the relaxation equation proved to be sufficient and efficient in practical applications. The study of the essentially quantum-mechanical low-temperature kinetic regime in its asymptotic limit requires the implementation of the exact relaxation equation. The coherent mechanism providing a non-vanishing reaction rate has been revealed when T→ 0. An accurate computational methodology for the cross-over kinetic regime needs a further elaboration. The original model of the hopping mechanism for electronic conduction in photosensitive organic materials is considered, based on the above techniques. The electron transfer (ET) in active centers of such systems proceeds via local intra- and intermolecular modes. The active modes, as a rule, operate beyond the kinetic regimes, which are usually postulated in the existing theories of the ET. Our alternative dynamic ET model for local modes immersed in the continuum harmonic medium is formulated for both classical and quantum regimes, and accounts explicitly for the mode/medium interaction. The kinetics of the energy exchange between the local ET subsystem and the surrounding environment essentially determine the total ET rate. The efficient computer code for rate computations is elaborated on. The computations are available for a wide range of system parameters, such as the temperature, external field, local mode frequency, and characteristics of mode/medium interaction. The relation of the present approach to the Marcus ET theory and to the quantum-statistical reaction rate theory [V. G. Levich and R. R. Dogonadze, Dokl. Akad. Nauk SSSR, Ser. Fiz. Khim. 124, 213 (1959); J. Ulstrup, Charge Transfer in Condensed Media (Springer, Berlin, 1979); M. Bixon and J. Jortner, Adv. Chem. Phys. 106, 35 (1999)] underlying it is discussed and illustrated by the results of computations for practically important target systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

The origin of linear scaling Fock matrix calculation with density prescreening

A theorem was proven, which reads that the number of nonzero two-electron integrals scales linearly with respect to the number of basis functions for large molecular systems. This permits to show that linear scaling property of the Fock matrix calculation with using density prescreening arises due to linear scaling properties of the number of nonzero two-electron integrals and the number of leading matrix elements of density matrix. This property is reinforced by employing the density prescreening technique. The use of the density difference prescreening further improves the linear scaling property of the Fock matrix calculation method. As a result, the linear scaling regime of the Fock matrix calculation can begin from the number of basis functions of 2000–3000 in dependence on the basis function type in molecular calculations. It was also shown that the conventional algorithm of Fock matrix calculation from stored nonzero two-electron integrals with density prescreening possesses linear scaling property.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Imaging of transition charge densities involving carbon core excitations by all X-ray sum-frequency generation

X-ray diffraction signals from the time-evolving molecular charge density induced by selective core excitation of chemically inequivalent carbon atoms are calculated. A narrowband X-ray pulse selectively excites the carbon K-edge of the –CH 3 or –CH 2 F groups in fluoroethane (CH 3 –CH 2 F). Each excitation creates a distinct core coherence which depends on the character of the electronic transition. Direct propagation of the reduced single-electron density matrix, using real-time time-dependent density functional theory, provides the time-evolving charge density following interactions with external fields. Here, the interplay between partially filled valence molecular orbitals upon core excitation induces characteristic femtosecond charge migration which depends on the core–valence coherence, and is monitored by the sum-frequency generation diffraction signal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Gamow-shell-model description of Li isotopes and their mirror partners

Background: Weakly bound and unbound nuclei close to particle drip lines are laboratories of new nuclear structure physics at the extremes of neutron/proton excess. The comprehensive description of these systems requires an open quantum system framework that is capable of treating resonant and nonresonant many-body states on equal footing. Purpose: In this work, we develop the complex-energy configuration interaction approach to describe binding energies and spectra of selected 5 ≤ $\textit{A}$ ≤ 11 nuclei. Method: We employ the complex-energy Gamow shell model (GSM) assuming a rigid 4 He core. The effective Hamiltonian, consisting of a core-nucleon Woods-Saxon potential and a simplified version of the Furutani-Horiuchi-Tamagaki interaction with the mass-dependent scaling, is optimized in the $\textit{sp}$ space. To diagonalize the Hamiltonian matrix, we employ the Davidson method and the Density Matrix Renormalization Group technique. Results: We report our optimized GSM Hamiltonian offers a good reproduction of binding energies and spectra with the root-mean-square (rms) deviation from experiment of 160 keV. Since the model performs well when used to predict known excitations that have not been included in the fit, it can serve as a reliable tool to describe poorly known states. A case in point is our prediction for the pair of unbound mirror nuclei 10 Li – 10 N in which a huge Thomas-Ehrman shift dramatically alters the pattern of low-energy excitations. Conclusion: The new model will enable comprehensive studies of structure and reactions aspects of light drip-line nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Hilbert space multireference coupled cluster tailored by matrix product states

In the past decade, the quantum chemical version of the density matrix renormalization group method has established itself as the method of choice for strongly correlated molecular systems. However, despite its favorable scaling, in practice, it is not suitable for computations of dynamic correlation. Several approaches to include that in post-DMRG methods exist; in our group, we focused on the tailored coupled cluster (TCC) approach. This method works well in many situations; however, in exactly degenerate cases (with two or more determinants of equal weight), it exhibits a bias toward the reference determinant representing the Fermi vacuum. Although sometimes it is possible to use a compensation scheme to avoid this bias for energy differences, it is certainly a drawback. In order to overcome this bias of the TCC method, we have developed a Hilbert-space multireference version of tailored CC, which can treat several determinants on an equal footing. We have implemented and compared the performance of three Hilbert-space multireference coupled cluster (MRCC) variants—the state universal one and the Brillouin–Wigner and Mukherjee’s state specific ones. We have assessed these approaches on the cyclobutadiene and tetramethyleneethane molecules, which are both diradicals with exactly degenerate determinants at a certain geometry. We have also investigated the sensitivity of the results on the orbital rotation of the highest occupied and lowest unoccupied molecular orbital (HOMO–LUMO) pair, as it is well known that Hilbert-space MRCC methods are not invariant to such transformations.

Chemistry↗
DOE OSTI2016

FatMan vs. LittleBoy: Scaling Up Linear Algebraic Operations in Scale-Out Data Platforms

Linear algebraic operations such as matrix manipulations form the kernel of many machine learning and other crucial algorithms. Scaling up as well as scaling out such algorithms are highly desirable to enable efficient processing over millions of data points. To this end, we present a matrix manipulation approach to effectively scale-up each node in a scale-out data parallel platform such as Apache Spark. Specifically, we enable hardware acceleration for matrix multiplications in a distributed Spark setup without user intervention. Our approach supports both dense and sparse distributed matrices, and provides flexible control of acceleration by matrix density. We demonstrate the benefit of our approach for generalized matrix multiplication operations over large matrices with up to four billion elements. To connect the effectiveness of our approach with machine learning applications, we performed Gramian matrix computation via generalized matrix multiplications. Our experiments show that our approach achieves more than 2× performance speed-up, and up to 96.1% computation improvement, compared to a state of the art Spark MLlib for dense matrices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

FatMan vs. LittleBoy: Scaling up Linear Algebraic Operations in Scale-out Data Platforms

Linear algebraic operations such as matrix manipulations form the kernel of many machine learning and other crucial algorithms. Scaling up as well as scaling out such algorithms are highly desirable to enable efficient processing over millions of data points. To this end, we present a matrix manipulation approach to effectively scale-up each node in a scale-out data parallel platform such as Apache Spark. Specifically, we enable hardware acceleration for matrix multiplications in a distributed Spark setup without user intervention. Our approach supports both dense and sparse distributed matrices, and provides flexible control of acceleration by matrix density. We demonstrate the benefit of our approach for generalized matrix multiplication operations over large matrices with up to four billion elements. To connect the effectiveness of our approach with machine learning applications, we performed Gramian matrix computation via generalized matrix multiplications. Our experiments show that our approach achieves more than 2x performance speed-up, and up to 96.1% computation improvement, compared to a state of the art Spark MLlib for dense matrices.

Lim, Seung-Hwan↗
DOE OSTI2024

Probing Quantum Telecloning on Superconducting Quantum Processors

Quantum information cannot be perfectly cloned, but approximate copies of quantum information can be generated. Quantum telecloning combines approximate quantum cloning, more typically referred to as quantum cloning, and quantum teleportation. Quantum telecloning allows approximate copies of quantum information to be constructed by separate parties, using the classical results of a Bell measurement made on a prepared quantum telecloning state. Quantum telecloning can be implemented as a circuit on quantum computers using a classical coprocessor to compute classical feedforward instructions using if statements based on the results of a midcircuit Bell measurement in real time. We present universal symmetric optimal 1 → M telecloning circuits and experimentally demonstrate these quantum telecloning circuits for M = 2 up to M = 10, natively executed with real-time classical control systems on IBM Quantum superconducting processors, known as dynamic circuits. We perform the cloning procedure on many different message states across the Bloch sphere, on seven IBM Quantum processors, optionally using the error suppression technique X–X sequence digital dynamical decoupling. Two circuit optimizations are utilized: one that removes ancilla qubits for M = 2,3 , and one that reduces the total number of gates in the circuit but still uses ancilla qubits. Parallel single-qubit tomography with maximum likelihood estimation density matrix reconstruction is used in order to compute the mixed-state density matrices of the clone qubits, and clone quality is measured using quantum fidelity. These results present one of the largest and most comprehensive noisy intermediate-scale quantum computer experimental analyses on (single qubit) quantum telecloning to date. The clone fidelity sharply decreases to 0.5 for M > 5 , but for M = 2 , we are able to achieve a mean clone fidelity of up to 0.79 using dynamical decoupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

DMRG on Top of Plane-Wave Kohn–Sham Orbitals: A Case Study of Defected Boron Nitride

In this paper, we analyze the numerical aspects of the inherent multi reference density matrix renormalization group (DMRG) calculations on top of the periodic Kohn-Sham density functional theory using the complete active space approach. The potential of the framework is illustrated by studying hexagonal boron nitride nanoflakes embedding a charged single boron vacancy point defect by revealing a vertical energy spectrum with a prominent multireference character. We investigate the consistency of the DMRG energy spectrum from the perspective of sample size, basis size, and active space selection protocol. Results obtained from standard quantum chemical atom-centered basis calculations and plane-wave based counterparts show excellent agreement. Furthermore, we also discuss the spectrum of the periodic sheet which is in good agreement with extrapolated data of finite clusters. These results pave the way toward applying the DMRG method in extended correlated solid-state systems, such as point defect qubit in wide band gap semiconductors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

AGR-5/6/7 Fuel Fabrication Report

The U.S. Department of Energy Office of Nuclear Energy (DOE NE) and the Idaho National Laboratory (INL) Advanced Reactor Technologies (ART) Advanced Gas Reactor (AGR) Fuel Development and Qualification program (referred to as AGR Fuel program hereafter) are pursuing qualification of tristructural isotropic (TRISO) coated particle fuel for use in high temperature gas cooled reactors (HTGRs). The AGR Fuel program was established to provide a fuel qualification data set in support of the licensing and operation of an HTGR. BWX Technologies Nuclear Operations Group (BWXT-NOG) was subcontracted to fabricate the fuel for the AGR program. Several investments and innovations were realized in preparation to fabricate fuel for the AGR-5/6/7 irradiation experiments that brought fuel fabrication fully out of the laboratory and into engineering-scale operations. These included: • Increased the kernel fabrication line capacity and uniformity • Upgraded ancillary support equipment and processes for the tristructural isotropic (TRISO) coating furnace • Demonstrated efficient production of the matrix precursor powder by dry jet milling of co mingled components • Demonstrated an engineering-scale method for quick and efficient overcoating TRISO particles with the matrix precursor • Demonstrated an automated, multi cavity compacting system with a volumetric feed system • Demonstrated a combined-cycle thermal treatment furnace These changes increased production rates of some of these processes by an order of magnitude or more while eliminating the use of flammable solvents, multiple grinding and sorting operations, and the weighing out of individual die charges. Three fuel kernel lots were fabricated for production of the fuel for AGR-5/6/7. The initial lot (J52R-16-39316) was certified to fuel specifications but was not used because the kernels had a high fraction of internal fissures that caused an unacceptable fraction of the kernels to fragment when charged to the coating furnace where the TRISO coating would be deposited. Fragmented kernels increased the dispersed uranium in the particles and produced a worrisome fraction of dimpled particles with an elevated probability of in-pile failure. After some efforts to identify the cause of the fissure formation, two additional lots were produced with much lower fissure fractions, J52R-16-69317 and 69318. The latter kernel lot was a backup to the first and was not needed. Multiple kernel batches were composited to form each of the lots so as to simulate a commercial-scale operation where kernel batches would also be composited. Multiple TRISO coating runs were performed and the product characterized so that several could be composited into a TRISO lot. TRISO lot J52R-16-98005 conformed to all fuel specifications except the mean outer pyrocarbon (OPyC) layer thickness was thinner than specified. Furthermore, it was determined that the TRISO lot had a dispersed uranium fraction (DUF) that might result in the compacts not meeting the DUF specification. A review of the role of the OPyC layer and consequences of the DUF by the Technical Coordination Team and INL concluded that the fuel was acceptable for use in the AGR-5/6/7 irradiation experiment. The TRISO particles were overcoated with the matrix precursor that had been produced in a jet-milling operation. The overcoating was performed in equipment originally designed to coat pharmaceuticals. The overcoating process performed well; producing highly spherical and uniform overcoats requiring little upgrading and no recycle or rework. TRISO particles were overcoated with the matrix precursor to achieve nominal volumetric packing fractions (PFs) of TRISO particles of 25% and 40% for the irradiation experiments. The 40% PF compacts occupy the first and fifth test capsule in the test train while the inner three capsules are loaded with 25% PF compacts. The resinated graphite matrix precursor powder was a derivative of the German A3-27 matrix formulation, which differs from previous AGR irradiation campaigns that used an A3-3 formulation. Jet milling of the matrix powder precursor produced a finer mean graphite particle size than the milling operations used for the A3-3 matrix powder precursor. Changes made in the matrix formula and equipment yielded compacts with significantly higher matrix density than was attained in previous AGR irradiation campaigns. The changes in the matrix formulation and the means of milling the powders also complicated resolution of the three fuel compact defect fractions, DUF, exposed kernel fraction (EKF), and the silicon carbide defect fraction (SDF). Characterization data from BWXT-NOG had some anomalous results, so samples of the fuel compacts and overcoated TRISO particles were also analyzed by Oak Ridge National Laboratory (ORNL) to ensure that the defect fractions were accurately characterized.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Prediction of orbital-selective Mott phases and block magnetic states in the quasi-one-dimensional iron chain Ce 2 O 2 FeSe 2 under hole and electron doping

The recent detailed study of quasi-one-dimensional iron-based ladders, with the 3d iron electronic density n=6, has unveiled surprises, such as orbital-selective phases. However, similar studies for n=6 iron chains are still rare. Here a three-orbital electronic Hubbard model was constructed to study the magnetic and electronic properties of the quasi-one-dimensional n=6 iron chain Ce 2 O 2 FeSe 2 , with focus on the effect of doping. Specifically, introducing the Hubbard U and Hund J H couplings and studying the model via the density matrix renormalization group, we report the ground-state phase diagram varying the electronic density away from n=6. For the realistic Hund coupling J H /U=1/4, several electronic phases were obtained, including a metal, orbital-selective Mott, and Mott insulating phases. Doping away from the parent phase, the competition of many tendencies leads to a variety of magnetic states, such as ferromagnetism, as well as several antiferromagnetic and magnetic “block” phases. In the hole-doping region, two different interesting orbital-selective Mott phases were found: OSMP1 (with one localized orbital and two itinerant orbitals) and OSMP2 (with two localized orbitals and one itinerant orbital). Moreover, charge disproportionation phenomena were found in special doping regions. We argue that our predictions can be tested by simple modifications in the original chemical formula of Ce 2 O 2 FeSe 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Tip‐Enhanced Imaging and Control of Infrared Strong Light‐Matter Interaction

Optical antenna resonators enable control of light‐matter interactions on the nano‐scale via electron–photon hybrid states in strong coupling. Specifically, mid‐infrared (MIR) nano‐antennas coupled to saturable intersubband transitions in multi‐quantum‐well (MQW) semiconductor heterostructures allow for the coupling strength to be tuned through antenna resonance and field intensity. Here, in this study, tip‐enhanced nano‐scale variation of antenna‐MQW coupling across the antenna is demonstrated, with a spatially‐dependent coupling strength $g_{\textrm{aq}}$ varying from 73 (strong coupling) to 24 cm -1 (weak coupling). This behavior is modeled based on the spatially dependent local constructive and destructive interference between tip and antenna fields. Using a quantum‐mechanical density‐matrix model of the MQW system with its designed values of transition dipole moment, doping density, and population decay time, the picosecond IR pulse coupling to intersubband transitions and the associated tip induced strong‐field saturation effects are described. These results present a new regime of nonlinear IR light‐matter control based on the dynamic manipulation of quantum hybrid states on the nanoscale and in the infrared, with a perspective regarding extension to molecular vibrations.

Wang, Yueying↗
DOE OSTI2017

Magnetic field dependent dynamics and field-driven metal-to-insulator transition of the half-filled Hubbard model: A DMFT+DMRG study

Here, we study the magnetic field-driven metal-to-insulator transition in half-filled Hubbard model on the Bethe lattice, using the dynamical mean-field theory by solving the quantum impurity problem with density-matrix renormalization group algorithm. The method enables us to obtain a high-resolution spectral densities in the presence of a magnetic field. It is found that the Kondo resonance at the Fermi level splits at relatively high magnetic field: the spin-up and -down components move away from the Fermi level and finally form a spin-polarized band insulator. By calculating the magnetization and spin susceptibility, we clarify that an applied magnetic field drives a transition from a paramagnetic metallic phase to a band insulating phase. In the weak interaction regime, the nature of the transition is continuous and captured by the Stoner's description, while in the strong interaction regime the transition is very likely to be metamagnetic, evidenced by the hysteresis curve. Furthermore, we determine the phase boundary by tracking the kink in the magnetic susceptibility, and the steplike change of the entanglement entropy and the entanglement gap closing. Interestingly, the phase boundaries determined from these two different ways are largely consistent with each other.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spin–Orbit Matrix Elements for a Combined Spin-Flip and IP/EA approach

We introduce a practical approach for computing the Breit–Pauli spin–orbit matrix elements of multiconfigurational systems with both spin and spatial degeneracies based on our recently developed RAS-nSF-IP/EA method (Houck, S. E.; et al. J. Chem. Theory Comput. 2019, 15, 2278). The spin–orbit matrix elements over all the multiplet components are computed using a single one-particle reduced density matrix as a result of the Wigner–Eckart theorem. A mean field spin–orbit approximation was used to account for the two-electron contributions. Basis set dependence as well as the effect of including additional excitations is presented. The effect of correlating the core and semicore orbitals is also examined. Surprisingly accurate results are obtained for spin–orbit coupling constants, despite the fact that the efficient wave function approximations we research neglect the bulk of dynamical correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗