Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GROUND STATES”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Strong-coupling expansion of multi-band interacting models: Mapping onto the transverse-field J 1 - J 2 Ising model

We investigate a class of two-dimensional two-band microscopic models in which the inter-band repulsive interactions play the dominant role. We first demonstrate three different schemes of constraining the ratios between the three types of inter-band interactions – density-density, spin exchange, and pair-hopping – that render the model free of the fermionic sign-problem for any filling and, consequently, amenable to efficient Quantum Monte Carlo simulations. We then study the behavior of these sign-problem-free models in the strong-coupling regime. In the cases where spin-rotational invariance is preserved or lowered to a planar symmetry, the strong-coupling ground state is a quantum paramagnet. However, in the case where there is only a residual Ising symmetry, the strong-coupling expansion maps onto the transverse-field J 1 - J 2 Ising model, whose pseudospins are associated with local inter-band magnetic order. We show that by varying the band structure parameters within a reasonable range of values, a variety of ground states and quantum critical points can be accessed in the strong-coupling regime, some of which are not realized in the weak-coupling regime. Here, we compare these results with the case of the single-band Hubbard model, where only intra-band repulsion is present, and whose strong-coupling behavior is captured by a simple Heisenberg model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Antiferromagnetic order in Co-doped Fe 5 GeTe 2 probed by resonant magnetic x-ray scattering

The quasi-two-dimensional van der Waals magnet Fe 5-δ GeTe 2 has emerged as a promising platform for electronic and spintronic functionalities at room temperature, owing to its large ferromagnetic ordering temperature T C ≈ 315 K. Interestingly, by cobalt (Co) substitution of iron in F5GT, i.e., (Fe 1-x Co x ) 5-δ GeTe 2 (Co-F5GT), not only can its magnetic transition temperature be further enhanced, but the magnetic and structural ground states can also be tuned. Specifically, an antiferromagnetic (AFM) order is induced beyond the Co doping level x ≥ 0.4. Here, in this work, we investigate the magnetic properties of a Co-F5GT single crystal at x = 0.45(1), by utilizing the element-specific, resonant magnetic x-ray scattering technique. Our study reveals an A-type, Ising-like AFM ground state, with a transition temperature T N ≈ 340 K. In addition, our work unveils an important contribution from Co magnetic moments to the magnetic order. The application of the in-plane magnetic fields gradually polarizes the spin moments along the field direction, but without inducing incommensurate spin texture(s).

36 MATERIALS SCIENCE↗
NASA NTRS1991

A conformal oxidation-resistant, plasma-polymerized coating

A comparative study was made of the surface recession (etching) of thin films of plasma polymerized tetrafluoro ethylene (PPTFE), polytetrafluoro ethylene (PTFE), and ion-beam sputter deposited polytetrafluoro ethylene (SPTFE) exposed to ground-state atomic oxygen downstream from a nonequilibrium radio-frequency O2 plasma. At 22 C, the etch rates for PTFE, SPTFE, and PPTFE were in the ratio of 8.7:1.8:1.0. A thin, conformal coating of PPTFE (etch rate of 0.3 nm/h at 22 C) was found to protect an underlying cast film of a reactive polymer, cis-1,4 polybutadiene, against ground-state atomic oxygen attack for the time required to fully etch away the PPTFE coating. From ESCA analysis, PTFE exhibited only minor surface oxidation (uptake of 0.5 atom percent O) upon etching, its F/C ratio decreasing slightly from 2.00 to 1.97; PPTFE exhibited considerable surface oxidation (uptake of 5.9 atom percent O) intermediate between those of PTFE and PPTFE, with a decrease in F/C ratio from 1.73 to 1.67. A plasma-polymerized fluorocarbon coating such as PPTFE might be useful for space applications to protect polymers that are vulnerable to oxidation or degradation by oxygen atoms.

Golub, Morton A.↗
DOE OSTI2025

Quenching of single-particle strength inferred from nucleon-removal transfer reactions on 15 C

The difference in the proton and neutron separation energies (ΔS) of the weakly bound 15 C ground state is -19.86 MeV, an extreme value. Data from intermediate-energy heavy-ion induced (HI-induced) knockout reactions on nuclei spanning -20 ≲ ΔS ≲ +20 MeV, suggest that the degree to which single-particle strength is quenched, R s , has a negative correlation with ΔS, decreasing from unity around -20 MeV to around 0.2 at +20 MeV. For the 15 C ground state (R s = 0.96 (4) in HI-induced knockout), contrasting results have recently been obtained via the neutron-adding transfer reaction, which reveal a value of R s = 0.64 (15), similar to the value observed at modest and more extreme values of ΔS with reaction probes other than HI knockout. In order to explore the any potential differences between adding and removing processes in transfer reactions at extreme ΔS, single-neutron removal transfer reactions on 15 C were performed at 7.1 MeV/u in inverse kinematics. The removal of a valence neutron in 2s 1/2 orbit using both (p, d) and (d, t) reactions shows consistent quenching factors and agrees with those from the neutron-adding reaction. The present results, which can be compared with neutron knockout reaction, suggest that correlations, represented by the quenching factor, show limited dependence on neutron-proton asymmetry under the most extreme asymmetry conditions so far achieved in transfer reactions.

quenching factor↗
DOE OSTI2020

Optimal Orbital Selection for Full Configuration Interaction (OptOrbFCI): Pursuing the Basis Set Limit under a Budget

Full configuration interaction (FCI) solvers are limited to small basis sets due to their expensive computational costs. An optimal orbital selection for FCI (OptOrbFCI) is proposed to boost the power of existing FCI solvers to pursue the basis set limit under a computational budget. The optimization problem coincides with that of the complete active space SCF method (CASSCF), while OptOrbFCI is algorithmically quite different. OptOrbFCI effectively finds an optimal rotation matrix via solving a constrained optimization problem directly to compress the orbitals of large basis sets to one with a manageable size, conducts FCI calculations only on rotated orbital sets, and produces a variational ground-state energy and its wave function. Coupled with coordinate descent full configuration interaction (CDFCI), we demonstrate the efficiency and accuracy of the method on the carbon dimer and nitrogen dimer under basis sets up to cc-pV5Z. We also benchmark the binding curve of the nitrogen dimer under the cc-pVQZ basis set with 28 selected orbitals, which provide consistently lower ground-state energies than the FCI results under the cc-pVDZ basis set. Furthermore, the dissociation energy in this case is found to be of higher accuracy.

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

Grand canonical inverse design of multicomponent colloidal crystals

Inverse design methods are powerful computational approaches for creating colloidal systems which self-assemble into a target morphology by reverse engineering the Hamiltonian of the system. Despite this, these optimization procedures tend to yield Hamiltonians which are too complex to be experimentally realized. An alternative route to complex structures involves the use of several different components, however, conventional inverse design methods do not explicitly account for the possibility of phase separation into compositionally distinct structures. We present an inverse design scheme for multicomponent colloidal systems by combining active learning with a method to directly compute their ground state phase diagrams. This explicitly accounts for phase separation and can locate stable regions of Hamiltonian parameter space which grid-based surveys are prone to miss. Using this we design low-density, binary structures with Lennard-Jones-like pairwise interactions that are simpler than in the single component case and potentially realizable in an experimental setting. This reinforces the concept that ground states of simple, multicomponent systems might be rich with previously unappreciated diversity, enabling the assembly of non-trivial structures with only few simple components instead of a single complex one.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Peierls transition, ferroelectricity, and spin-singlet formation in monolayer VOI 2

Using ab initio density functional theory and single-orbital Hubbard model calculations via the density matrix renormalization group method, we systematically studied the monolayer VOI 2 with a 3d 1 electronic configuration. Our phonon calculations indicate that the orthorhombic Pmm2 FE-II phase is the most likely ground state, involving a ferroelectric (FE) distortion along the a axis and V-V dimerization along the b axis. Specifically, the “pseudo Jahn-Teller” effect caused by the coupling between empty V (d xz/yz and d 3z 2 -r 2 ) and O 2p states is proposed as the mechanism that stabilizes the FE distortion from the paraelectric phase. Moreover, the half-filled metallic d xy band displays a Peierls instability along the b axis, inducing a V-V dimerization. Additionally, we found very short-range antiferromagnetic coupling along the V-V chain due to the formation of nearly decoupled spin singlets in the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Structural, Electrical, and Electromechanical Properties of Inverse Hybrid Perovskites from First-Principles: The Case of (CH 3 NH 3 ) 3 OI

Inverse-hybrid perovskites (IHPs) with large polarization have recently been predicted from first-principles computations. Here, we use one representative from the IHP class of materials, (CH 3 NH 3 ) 3 OI (MA3OI), to propose a route to the first-principles prediction of structural and electrical properties, such as polarization, polarization reversibility, and the associated coercive field for hybrid organic–inorganic perovskites. The route relies on the construction of the polarization reversal path that models experimental measurements. Such a path was found to play an important role in the ground-state search as well as in the identification of competing structural variations. The latter is believed to be the origin of the structural disorder that is characteristic of hybrid organic–inorganic perovskites. The application of such an approach to MA3OI leads to the prediction of several structural variations that are expected to result in a structurally disordered phase above 766 K and of the polar ground state with a polarization of 25.3 μC/cm 2 that is reversible with the application of an electric field. The upper estimate for the coercive field associated with homogeneous polarization reversal is 6.9 GV/m. The piezoelectric constants of MA3OI are predicted to be an order of magnitude smaller in comparison with a prototypical inorganic ferroelectric PbTiO 3 ; however, the low symmetry of the MA3OI structure yields finite values for all components of the piezoelectric tensor. The polarization in MA3OI is tunable by the epitaxial strain (11.5% change under 5% epitaxial strain), although less so as compared with PbTiO 3 .

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

Sudden Collapse of Magnetic Order in Oxygen-Deficient Nickelate Films

Antiferromagnetic order is a common and robust ground state in the parent (undoped) phase of several strongly correlated electron systems. The progressive weakening of antiferromagnetic correlations upon doping paves the way for a variety of emergent many-electron phenomena including unconventional superconductivity, colossal magnetoresistance, and collective charge-spin-orbital ordering. In this study, we explored the use of oxygen stoichiometry as an alternative pathway to modify the coupled magnetic and electronic ground state in the family of rare earth nickelates ( RENiO 3 - x ). Using a combination of x-ray spectroscopy and resonant soft x-ray magnetic scattering, we find that, while oxygen vacancies rapidly alter the electronic configuration within the Ni and O orbital manifolds, antiferromagnetic order is remarkably robust to substantial levels of carrier doping, only to suddenly collapse beyond 0.21 e - / Ni without an accompanying structural transition. Our work demonstrates that ordered magnetism in RENiO 3 - x is mostly insensitive to carrier doping up to significant levels unseen in other transition-metal oxides. The sudden collapse of ordered magnetism upon oxygen removal may provide a new mechanism for solid-state magnetoionic switching and new applications in antiferromagnetic spintronics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Experimental evaluation of the stabilization of the COT orbitals by 4f orbitals in COT 2 Ce using a Hubbard model

Significant orbital mixing is rare in lanthanide complexes because of the limited radial extent of the 4f orbitals, which results in a generally small stabilization due to 4f orbital interactions. Nevertheless, even a small amount of additional stabilization could enhance lanthanide separations. One lanthanide complex in which orbital mixing has been extensively studied both experimentally and computationally is cerocene, COT 2 Ce, where COT is cyclooctatetraene. This compound has a singlet ground state with a low-lying, triplet excited state. Previous fluorescence studies on trimethylsilyl-substituted cerocenes indicate the triplet state is 0.4 eV higher in energy than the singlet state. In addition, computational studies predict that the triplet is 0.3 to 1 eV higher in energy than the singlet. The synthesis of highly pure COT 2 Ce by Walter and Andersen allowed its physical properties to be accurately measured. Using these measurements, we evaluate the stabilization of the 4f orbitals using two, independent approaches. A Hubbard model is used to evaluate the stabilization of the ground state due to orbital mixing. This stabilization, which is also the singlet–triplet gap, is –0.29 eV using this model. Furthermore, this gap was also from the temperature independent paramagnetism of COT 2 Ce, which yielded a value of –0.32 eV.

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

Superconductivity with T c ≈ 7 K under pressure for Cu- and Au-doped BaFe 2 As 2

It is noteworthy that chemical substitution of BaFe 2 As 2 (122) with the noble elements Cu and Au gives superconductivity with a maximum T c ≈ 3 K, while Ag substitution (Ag-122) stays antiferromagnetic. For Ba(Fe 1–x TM x ) 2 As 2 , TM = Cu, Au, or Ag, and by doping an amount of x = 0.04, a-lattice parameter slightly increases (0.4%) for all TM dopants, while c-lattice decreases (–0.2%) for TM = Cu, barely moves (0.05%) for Au, and increases (0.2%) for Ag. Despite the naive expectation that the noble elements of group 11 should affect the quantum properties of 122 similarly, they produce significant differences extending to the character of the ground state. For the Ag-122 crystal, evidence of only a filamentary superconductivity is noted with pressure. However, for Au and Cu doping (x ≈ 0.03) we find a substantial improvement in the superconductivity, with T c increasing to 7 K and 7.5 K, respectively, under 20 kbar of pressure. As with the ambient pressure results, the identity of the dopant therefore has a substantial impact on the ground state properties. Further, density functional theory calculations corroborate these results and find evidence of strong electronic scattering for Au and Ag dopants, while Cu is comparatively less disruptive to the 122 electronic structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Constructing Hubbard models for the hydrogen chain using sliced-basis density matrix renormalization group

Sliced-basis DMRG(sb-DMRG) is used to simulate a chain of hydrogen atoms and to construct low-energy effective Hubbard-like models. The downfolding procedure first involves a change of basis to a set of atom-centered Wannier functions constructed from the natural orbitals of the exact DMRG one-particle density matrix. The Wannier function model is then reduced to a fewer-parameter Hubbard-like model, whose parameters are determined by minimizing the expectation value of the Wannier Hamiltonian in the ground state of the Hubbard Hamiltonian. This indirect variational procedure not only yields compact and simple models for the hydrogen chain, but also allows us to explore the importance of constraints in the effective Hamiltonian, such as the restricting the range of the single-particle hopping and two-particle interactions, and to assess the reliability of more conventional downfolding. The entanglement entropy for a model's ground state, cut in the middle, is an important property determining the ability of DMRG and tensor networks to simulate the model, and here we study its variation with the range of the interactions. Counterintuitively, we find that shorter ranged interactions often have larger entanglement.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Lowest ionization potentials of Al2

Potential curves for the lowest two electronic states, X 2Sigma(+)g and A 2Pi(u), of Al2(+) were computed using complete active space SCF/multireference CI wave functions and large Gaussian basis sets. The lowest observable vertical ionization potential, to Al2(+) X 2Sigma(+)g, of the Al2 X 3 Pi(u) ground state is calculated to occur around 6.1 eV, in excellent agreement with the experimental range of 6.0 to 6.42 eV obtained in recent cluster ionization studies by Cox and co-workers. The second vertical ionization potential, to Al2(+) A 2Pi(u), occurs near 6.4 eV, also within the experimental range. The adiabatic IP of 5.90 eV is in good agreement with the value of 5.8 to 6.1 eV deduced by Hanley and co-workers from the difference in thresholds between collision induced dissociation processes of Al3(+). The computed IP values are somewhat larger than those deduced from branching ratios in cluster fragmentation experiments by Jarrold and co-workers. The observation of an ionization threshold below 6.42 eV is shown to be incompatible with an Al2 ground electronic state assignment of 3Sigma(-)g, but the separation between the two lowest states of Al2 is so small that it is likely that both are populated in the experiments, so that this does not provide unambiguous support for the recent theoretical assignment of the ground state as 3Pi(u).

Bauschlicher, Charles W., Jr.↗
DOE OSTI2020

Iron(III)–Oxo Cluster Chemistry with Dimethylarsinate Ligands: Structures, Magnetic Properties, and Computational Studies

A program has been initiated to develop Fe III /oxo cluster chemistry with the ‘pseudo-carboxylate’ ligand dimethylarsinate (Me 2 AsO 2 - ) for comparison with the well investigated Fe III /oxo/carboxylate cluster area. The synthesis and characterization of three polynuclear Fe III complexes are reported, [Fe 12 O 4 (O 2 C t Bu) 8 (O 2 AsMe 2 ) 17 (H 2 O) 3 ]Cl 3 (1), Na 2 [Fe 12 Na 2 O 4 (O 2 AsMe 2 ) 20 (NO 3 ) 6 (Me 2 AsO 2 H) 2 (H 2 O) 4 ](NO 3 ) 6 (2) and [Fe 3 (O 2 AsMe 2 ) 6 (Me 2 AsO 2 H) 2 (hqn) 2 ](NO 3 ) (3), where hqn is 8-hydroxyquinoline. The Fe 12 core of 1 is a type never previously encountered in Fe III carboxylate chemistry, consisting of two Fe 6 units each of which comprises two {Fe 3 (μ 3 -O 2- )} units bridged by three Me 2 AsO 2 - groups and linked into an Fe 12 loop structure by two anti-anti η 1 :η 1 :μ Me 2 AsO 2 - groups, a bridging mode extremely rare with carboxylates. 2 also consists of two Fe 6 units, differing in their ligation from those in 1, and this time linked together into a linear structure by a central {Na 2 (NO 3 ) 2 } bridging unit. 3 is a linear Fe 3 complex with no monoatomic bridges between Fe III ions, a very rare situation in Fe III chemistry with any ligands, and unprecedented in Fe carboxylate chemistry. The distinct differences observed in arsinate vs carboxylate ligation modes are rationalized largely based on the greater basicity of the former vs the latter. Variable-temperature dc and ac magnetic susceptibility data reveal all Fe 2 pairwise interactions to be antiferromagnetic. For 1 and 2, the different J ij couplings were estimated by use of a magnetostructural correlation for high nuclearity Fe III -oxo clusters and by density functional theory calculations using broken symmetry methods, allowing identification of their relative spin vector alignments and thus rationalization of their S = 0 ground states. The J ij values were then used as input values to give excellent fits of the experimental χM T vs T data. For 3, the fits of the experimental χM T vs T data to the Van Vleck equation or with PHI gave a very weak J 12 = -0.8(1) cm -1 (H = –2JŠ i ·Š j convention) between adjacent Fe III ions, and an S = 5/2 ground state. Furthermore, these initial Fe III arsinate complexes also provide structural parameters that help validate literature assignments of arsinate binding modes to iron oxide/hydroxide minerals as part of environmental concerns of using arsenic-containing herbicides in agriculture.

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

Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C 36 Fullerene and an Iron Porphyrin Model Complex

We present three distinct examples where phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) can be reliably performed with a single-determinant trial wave function with essential symmetry breaking. Essential symmetry breaking was first introduced by Lee and Head-Gordon [ Phys. Chem. Chem. Phys. 2019, 21, 4763-4778, 10.1039/C8CP07613H]. We utilized essential complex and time-reversal symmetry breaking with ph-AFQMC to compute the triplet–singlet energy gap in the TS12 set. We found statistically better performance of ph-AFQMC with complex-restricted orbitals than with spin-unrestricted orbitals. We then showed the utilization of essential spin symmetry breaking when computing the singlet–triplet gap of a known biradicaloid, C 36 . ph-AFQMC with spin-unrestricted Hartree–Fock (ph-AFQMC+UHF) fails catastrophically even with spin-projection and predicts no biradicaloid character. With approximate Brueckner orbitals obtained from regularized orbital-optimized second-order Møller–Plesset perturbation theory (κ-OOMP2), ph-AFQMC quantitatively captures strong biradicaloid character of C 36 . Lastly, we applied ph-AFQMC to the computation of the quintet–triplet gap in a model iron porphyrin complex where brute-force methods with a small active space fail to capture the triplet ground state. We show unambiguously that neither triplet nor quintet is strongly correlated using UHF, κ-OOMP2, and coupled-cluster with singles and doubles (CCSD) performed on UHF and κ-OOMP2 orbitals. There is no essential symmetry breaking in this problem. By virtue of this, we were able to perform UHF+ph-AFQMC reliably with a cc-pVTZ basis set and predicted a triplet ground state for this model geometry. The largest ph-AFQMC in this work correlated 186 electrons in 956 orbitals. Our work highlights the utility, scalability, and accuracy of ph-AFQMC with a single-determinant trial wave function with essential symmetry breaking for systems mainly dominated by dynamical correlation with little static correlation.

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

Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo

The ability to accurately predict the finite temperature properties and phase diagrams of realistic quantum solids is central to uncovering new phases and engineering materials with novel properties ripe for device applications. Nonetheless, there remain comparatively few many-body techniques capable of elucidating the finite temperature physics of solids from first principles. Here, we take a significant step toward developing such a technique by generalizing our previous, exact fully ab initio finite temperature Auxiliary Field Quantum Monte Carlo (FT-AFQMC) method to model periodic solids and employing it to uncover the finite temperature physics of periodic hydrogen chains. Our chains’ unit cells consist of 10 hydrogen atoms modeled in a minimal basis, and we sample 5 k-points from the first Brillouin zone to arrive at a supercell consisting of 50 orbitals and 50 electrons. Based upon our calculations of these chains’ many-body energies, free energies, entropies, heat capacities, double and natural occupancies, and charge and spin correlation functions, we outline their metal–insulator and magnetic ordering as a function of both H–H bond distance and temperature. At low temperatures approaching the ground state, we observe both metal–insulator and ferromagnetic–antiferromagnetic crossovers at bond lengths between 0.5 and 0.75 Å. We then demonstrate how this low-temperature ordering evolves into a metallic phase with decreasing magnetic order at higher temperatures. In order to contextualize our results, we compare the features we observe to those previously seen in one-dimensional, half-filled Hubbard models at finite temperature and in ground state hydrogen chains. Interestingly, we identify signatures of the Pomeranchuk effect in hydrogen chains for the first time and show that spin and charge excitations that typically arise at distinct temperatures in the Hubbard model are indistinguishably coupled in these systems. Beyond qualitatively revealing the many-body phase behavior of hydrogen chains in a numerically exact manner without invoking the phaseless approximation, our efforts shed light on the further theoretical developments that will be required to construct the phase diagrams of the more complex transition metal, lanthanide, and actinide solids of longstanding interest to physicists.

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

Towards accurate orbital-free simulations: A generalized gradient approximation for the noninteracting free energy density functional

For orbital-free ab initio molecular dynamics, especially on systems in extreme thermodynamic conditions, we provide the first pseudo-potential-adapted generalized gradient approximation (GGA) functional for the non-interacting free energy. This is achieved by systematic finite-temperature extension of our recent LKT ground state non-interacting kinetic energy GGA functional (Phys. Rev. B 98, 041111(R) (2018)). We test the performance of the new functional first via static lattice calculations on crystalline aluminum and silicon. Then we compare deuterium equation of state results against both path-integral Monte Carlo and conventional (orbital-dependent) Kohn-Sham results. The new functional, denoted LKTF, outperforms the previous best semi-local free energy functional, VT84F (Phys. Rev. B 88, 161108(R) (2013)), and provides modestly faster simulations. Finally, we discuss subtleties of identification of kinetic and entropic contributions to non-interacting free-energy functionals obtained by extension from ground state orbital-free kinetic energy functionals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Accuracy of ghost rotationally invariant slave-boson and dynamical mean field theory as a function of the impurity-model bath size

Here, we compare the accuracy of the ghost rotationally invariant slave-boson (g-RISB) theory and dynamical mean field theory (DMFT) on the single-band Hubbard model, as a function of the number of bath sites in the embedding impurity Hamiltonian. Our benchmark calculations confirm that the accuracy of g-RISB can be systematically improved by increasing the number of bath sites, similar to DMFT. With a few bath sites, we observe that g-RISB is systematically more accurate than DMFT for the ground-state observables. On the other hand, the relative accuracy of these methods is generally comparable for the quasiparticle weight and the spectral function. As expected, we observe that g-RISB satisfies the variational principle in infinite dimensions, as the total energy decreases monotonically towards the exact value as a function of the number of bath sites, suggesting that the g-RISB wave function may approach the exact ground state in infinite dimensions. Our results suggest that the g-RISB is a promising method for first-principles simulations of strongly correlated matter, which can capture the behavior of both static and dynamical observables, at a relatively low computational cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗