Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional renormalization group”

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 73 records · Page 4

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spatially anisotropic 𝑆=1 square-lattice antiferromagnet with single-ion anisotropy realized in a Ni(II) pyrazine-𝑛,𝑛′-dioxide coordination polymer

The Ni(NCS) 2 (pyzdo) 2 coordination polymer is found to be an S = 1 spatially anisotropic square lattice with easy-axis single-ion anisotropy. This conclusion is based upon considering in concert the experimental probes x-ray diffraction, magnetic susceptibility, magnetic-field-dependent heat capacity, muon-spin relaxation, neutron diffraction, neutron spectroscopy, and pulsed-field magnetization. Long-range antiferromagnetic (AFM) order develops at T N = 18.5 K. Although the samples are polycrystalline, there is an observable spin-flop transition and saturation of the magnetization at ≈80 T. Linear spin-wave theory yields spatially anisotropic exchanges within an AFM square lattice, J x =0.235 meV, J y =2.014 meV, and an easy-axis single-ion anisotropy D=−1.622 meV (after renormalization). The anisotropy of the exchanges is supported by density functional theory.

2-dimensional systems↗

Bypassing the Energy Functional in Density Functional Theory: Direct Calculation of Electronic Energies from Conditional Probability Densities

Density functional calculations can fail for want of an accurate exchange-correlation approximation. The energy can instead be extracted from a sequence of density functional calculations of conditional probabilities (CP DFT). Simple CP approximations yield usefully accurate results for two-electron ions, the hydrogen dimer, and the uniform gas at all temperatures. CP DFT has no self-interaction error for one electron, and correctly dissociates H2, both major challenges. For warm dense matter, classical CP DFT calculations can overcome the convergence problems of Kohn-Sham DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Origin of Insulating Ferromagnetism in Iron Oxychalcogenide Ce 2 O 2 FeSe 2

An insulating ferromagnetic (FM) phase exists in the quasi-one-dimensional iron oxychalcogenide Ce 2 O 2 FeSe 2 , but its origin is unknown. To understand the FM mechanism, here a systematic investigation of this material is provided, analyzing the competition between ferromagnetic and antiferromagnetic tendencies and the interplay of hoppings, Coulomb interactions, Hund’s coupling, and crystal-field splittings. Overall, our intuitive analysis based on second-order perturbation theory shows that large entanglements between doubly occupied and half filled orbitals play a key role in stabilizing the FM order in Ce 2 O 2 FeSe 2 . In addition, via many-body computational techniques applied to a multiorbital Hubbard model, the phase diagram confirms the proposed FM mechanism.

1-dimensional spin chains↗

Monte Carlo renormalization-group investigation of the two-dimensional O(4) sigma model

An improved Monte Carlo renormalization-group method is used to determine the beta function of the two-dimensional O(4) sigma model. While for (inverse) couplings beta = greater than about 2.2 agreement is obtained with asymptotic scaling according to asymptotic freedom, deviations from it are obtained at smaller couplings. They are, however, consistent with the behavior of the correlation length, indicating 'scaling' according to the full beta function. These results contradict recent claims that the model has a critical point at finite coupling.

Heller, Urs M.↗

Two-loop radiative jet function for exclusive B-meson and Higgs decays

The rare radiative B -meson decay \( {B}^{-}\to \gamma {\mathrm{\ell}}^{-}\overline{v} \) and the radiative Higgs-boson decay h → γγ mediated by light-quark loops both receive large logarithmic corrections in QCD, which can be resummed using factorization theorems derived in soft-collinear effective theory. In these factorization theorems the same radiative jet function appears, which is a central object in the study of factorization beyond the leading order in scale ratios. We calculate this function at two-loop order both in momentum space and in a dual space, where its renormalization-group evolution takes on a simpler form. We also derive the two-loop anomalous dimension of the jet function and present the exact solution to its evolution equation at two-loop order. Another important outcome of our analysis is the explicit form of the two-loop anomalous dimension of the B -meson light-cone distribution amplitude in momentum space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

c -functions in higher-derivative flows across dimensions

In the context of gravitational theories describing renormalization group flows across dimensions via AdS/CFT, we study the role of higher-derivative corrections to Einstein gravity. We use the Null Energy Condition to derive monotonicity properties of candidate holographic central charges formed by combinations of metric functions. We also implement an entropic approach to the characterization of the four-derivative flows using the Jacobson-Myers functional and demonstrate, under reasonable conditions, monotonicity of certain terms in the entanglement entropy via the appropriate generalization of the Ryu-Takayanagi prescription. In particular, we show that any flow from a higher dimensional theory to a holographic CFT 2 satisfies a type of monotonicity. We also uncover direct relations between NEC-motivated and entropic central charges.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

3D Heisenberg universality in the van der Waals antiferromagnet NiPS 3

Van der Waals (vdW) magnetic materials are comprised of layers of atomically thin sheets, making them ideal platforms for studying magnetism at the two-dimensional (2D) limit. These materials are at the center of a host of novel types of experiments, however, there are notably few pathways to directly probe their magnetic structure. We confirm the magnetic order within a single crystal of NiPS 3 and show it can be accessed with resonant elastic X-ray diffraction along the edge of the vdW planes in a carefully grown crystal by detecting structurally forbidden resonant magnetic X-ray scattering. We find the magnetic order parameter has a critical exponent of β ~ 0.36, indicating that the magnetism of these vdW crystals is more adequately characterized by the three-dimensional (3D) Heisenberg universality class. We verify these findings with first-principles density functional theory, Monte-Carlo simulations, and density matrix renormalization group calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spectroscopic signatures of next-nearest-neighbor hopping in the charge and spin dynamics of doped one-dimensional antiferromagnets

In this work, we study the impact of next-nearest-neighbor (nnn) hopping on the low-energy collective excitations of strongly correlated doped antiferromagnetic cuprate spin chains. Specifically, we use exact diagonalization and the density matrix renormalization group method to study the single-particle spectral function, the dynamical spin and charge structure factors, and the Cu L-edge resonant inelastic x-ray scattering (RIXS) intensity of the doped t-t'-J model for a set of t' values. We find evidence for the breakdown of spin-charge separation as |t'| increases and identify its fingerprints in the dynamical response functions. The inclusion of nnn hopping couples the spinon and holon excitations, resulting in the formation of a spin-polaron, where a ferromagnetic spin polarization cloud dresses the doped carrier. The spin-polaron manifests itself as additional spectral weight in the dynamical correlation functions, which appear simultaneously in the spin- and charge-sensitive channels. We also demonstrate that RIXS can provide a unique view of the spin-polaron, due to its sensitivity to both the spin and charge degrees of freedom.

1D spin chains↗

Ferromagnetism and doublon localization in a Wannier-Hubbard chain

We derive a “Wannier-Hubbard” model consisting of an array of overlapping atomic orbitals interacting via a local Coulomb interaction. Transforming to an orthogonal Wannier basis set, the resulting Hamiltonian displays long-range hopping and interactions, with new terms such as correlated hopping and ferromagnetic direct exchange, among others. We numerically study the one-dimensional version of the model at half-filling using the density matrix renormalization group method, unveiling a rich phase diagram as a function of the interaction U and the overlap s with metallic and ferromagnetic phases separated by a ferrimagnetic region. Our results indicate a path toward understanding emergent phases under pressure and beyond standard model Hamiltonians.

36 MATERIALS SCIENCE↗

Charge Radii Measurements of Exotic Tin Isotopes in the Proximity of 𝑁 =50 and 𝑁 =82

We report nuclear charge radii for the isotopes 104–134 Sn , measured using two different collinear laser spectroscopy techniques at ISOLDE-CERN. These measurements clarify the archlike trend in charge radii along the isotopic chain and reveal an odd-even staggering that is more pronounced near the 𝑁 =50 and 𝑁 =82 shell closures. The observed local trends are well described by both nuclear density functional theory and valence space in-medium similarity renormalization group calculations. Both theories predict appreciable contributions from beyond-mean-field correlations to the charge radii of the neutron-deficient tin isotopes. The models, however, fall short of reproducing the magnitude of the known 𝐵⁡(𝐸⁢2) transition probabilities, highlighting the remaining challenges in achieving a unified description of both ground-state properties and collective phenomena.

laser spectroscopy↗

Charge radii measurements of exotic tin isotopes in the proximity of $N=50$ and $N=82$

We report nuclear charge radii for the isotopes $^{104-134}$Sn, measured using two different collinear laser spectroscopy techniques at ISOLDE-CERN. These measurements clarify the archlike trend in charge radii along the isotopic chain and reveal an odd-even staggering that is more pronounced near the $N=50$ and $N=82$ shell closures. The observed local trends are well described by both nuclear density functional theory and valence space in-medium similarity renormalization group calculations. Both theories predict appreciable contributions from beyond-mean-field correlations to the charge radii of the neutron-deficient tin isotopes. The models fall short, however, of reproducing the magnitude of the known $B(E2)$ transition probabilities, highlighting the remaining challenges in achieving a unified description of both ground-state properties and collective phenomena.

FOS: Physical sciences↗

237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes

A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

CFT and lattice correlators near an RG domain wall between minimal models

Conformal interfaces separating two conformal field theories (CFTs) provide maps between different CFTs, and naturally exist in nature as domain walls between different phases. One particularly interesting construction of a conformal interface is the renormalization group (RG) domain wall between CFTs. For a given Virasoro minimal model \mathcal{M}_{k+3,k+2} ℳ k + 3 , k + 2 , an RG domain wall can be generated by a specific deformation which triggers an RG flow towards its adjacent Virasoro minimal model \mathcal{M}_{k+2,k+1} ℳ k + 2 , k + 1 with the deformation turned on over part of the space. An algebraic construction of this domain wall was proposed by Gaiotto in [J. High Energy Phys. 12, 103 (2012)]. In this paper, we will provide a study of this RG domain wall for the minimal case k=2 k = 2 , which can be thought of as a nonperturbative check of the construction. In this case the wall is separating the Tricritical Ising Model (TIM) CFT and the Ising Model (IM) CFT. We will check the analytical results of correlation functions from the RG brane construction with the numerical density matrix renormalization group (DMRG) calculation using a lattice model proposed in [arXiv.1206.1332, Science 344, 280 (2014)], and find a perfect agreement. We comment on possible experimental realizations of this RG domain wall.

Cogburn, Cameron V.↗

Exploring Hilbert space on a budget: Novel benchmark set and performance metric for testing electronic structure methods in the regime of strong correlation

This work explores the ability of classical electronic structure methods to efficiently represent (compress) the information content of full configuration interaction (FCI) wave functions. We introduce a benchmark set of four hydrogen model systems of different dimensionalities and distinctive electronic structures: a 1D chain, a 1D ring, a 2D triangular lattice, and a 3D close-packed pyramid. To assess the ability of a computational method to produce accurate and compact wave functions, we introduce the accuracy volume, a metric that measures the number of variational parameters necessary to achieve a target energy error. Using this metric and the hydrogen models, we examine the performance of three classical deterministic methods: (i) selected configuration interaction (sCI) realized both via an a posteriori (ap-sCI) and variational selection of the most important determinants, (ii) an a posteriori singular value decomposition (SVD) of the FCI tensor (SVD-FCI), and (iii) the matrix product state representation obtained via the density matrix renormalization group (DMRG). We find that the DMRG generally gives the most efficient wave function representation for all systems, particularly in the 1D chain with a localized basis. For the 2D and 3D systems, all methods (except DMRG) perform best with a delocalized basis, and the efficiency of sCI and SVD-FCI is closer to that of DMRG. For larger analogs of the models, the DMRG consistently requires the fewest parameters but still scales exponentially in 2D and 3D systems, and the performance of SVD-FCI is essentially equivalent to that of ap-sCI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗