Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correlation energy”

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 145 records · Page 8

Cohesion and excitations of diamond-structure silicon by quantum Monte Carlo: Benchmarks and control of systematic biases

We have carried out quantum Monte Carlo (QMC) calculations of silicon crystal focusing on the accuracy and systematic biases that affect the electronic structure characteristics. The results show that 64 and 216 atom supercells provide an excellent consistency for extrapolated energies per atom in the thermodynamic limit for ground, excited, and ionized states. We have calculated the ground state cohesion energy with both systematic and statistical errors below ≈ 0.05 eV. The ground state exhibits a fixed-node error of only 1.3(2)% of the correlation energy, suggesting an unusually high accuracy of the corresponding single-reference trial wave function. We obtain a very good agreement between optical and quasiparticle gaps that affirms the marginal impact of excitonic effects. Our most accurate results for band gaps differ from the experiments by about 0.2 eV. This difference is assigned to a combination of residual finite-size and fixed-node errors. In conclusion, we have estimated the crystal Fermi level referenced to vacuum that enabled us to calculate the edges of valence and conduction bands in agreement with experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Photonuclear cross sections for 197 Au: An update on the gold standard

Cross sections for the 197 Au(γ, n) reaction are broadly used in nuclear physics as a standard for normalizing photonuclear reaction cross-section data at photon beam energies above approximately 8 MeV. In this paper, we report cross-section measurements for the 197 Au(γ, n) 196 Au g+m1 reaction at beam energies from 13 to 31 MeV. Our measurements provide the first cross-section data for this reaction at beam energies above 20 MeV, enabling the use of this reaction as a cross-section standard up to 30 MeV. Also, this work provides first cross-section measurements for the 197 Au(γ, n) 196 Au m2 reaction. In addition, we measured cross-section data for the 197 Au(γ, 3n) 194 Au reaction, which can be used as a cross-section standard above about 25 MeV. These measurements were performed using a new target activation method that is based on the angle-energy correlation of the laser Compton- scattered photon beams at the High Intensity Gamma-ray Source (HIγS). The technique enables measuring photonuclear reaction cross-sections at several discrete beam energies concurrently via a single irradiation on a stack of different targets. Measurements were carried out by irradiating a stack of concentric-ring targets consisting of Au, TiO 2 , Zn, Os, and Au (in order of the γ-ray beam direction). Our data for the 197 Au(γ,n) 196 Au g+m1 reaction in the energy range of 13 to 20 MeV are in good agreement with existing ones measured using monoenergetic γ-ray beams, but differ from data acquired using a bremsstrahlung γ-ray beam. Also, above 18 MeV, our data for the 197 Au(γ,n) 196 Au g+m1 and 197 Au(γ,n) 196 Au m2 reactions differ significantly from the most recent TENDL and JENDL evaluations, suggesting a need to update these data libraries. Furthermore, the TENDL evaluation and existing data are consistent with our data for the 197 Au(γ,3n) reaction, but differ significantly from the JENDL evaluation above 26 MeV.

190 ≤ A ≤ 219↗

Fully thermal meta-GGA exchange correlation free-energy density functional

The application of density functional theory to materials in the warm dense matter regime has motivated the development of exchange-correlation functionals which incorporate proper, explicit temperature dependence. Previous work has yielded fully-thermal exchange-correlation free energy functionals at the local density approximation (LDA) and generalized gradient approximation (GGA) levels of refinement. Recently an additive thermal correction scheme was utilized to construct a meta-GGA exchange-correlation (XC) functional in which thermal effects are treated at the GGA level. Here, the f TSCAN free-energy XC functional presented here includes thermal effects through the meta-GGA level in the context of the SCAN (strongly constrained and appropriately normed) ground-state functional. The f TSCAN functional provides generality while achieving similar performance to a thermal GGA functional at high temperatures, e.g. pressures within 1% of path integral Monte Carlo simulations of warm dense hydrogen, and a significant improvement over ground-state functionals. At low temperatures, f TSCAN demonstrates improvements in accuracy relative to lower-level and deorbitalized functionals, indicating that calculations using f TSCAN may be expected to perform well across experimentally relevant densities and pressures.

36 MATERIALS SCIENCE↗

Many-Body Basis Set Amelioration Method for Incremental Full Configuration Interaction

Incremental full configuration interaction (iFCI) is a polynomial-cost electronic structure method that systematically approaches the FCI limit by employing the method of increments to solve the Schrödinger equation through a many-body expansion. This article introduces the many-body basis set amelioration (MBBSA) method, which is designed to allow iFCI to be applicable to larger atomic orbital basis sets. MBBSA uses a series of inexpensive iFCI calculations to approximate the correlation energy that would be found using a more expensive, highly accurate iFCI calculation. Here, when compared to standard iFCI computations on smaller molecules in triple-zeta and larger basis sets, MBBSA provides approximations to the total and relative energies within chemical accuracy. MBBSA exhibits a reduced cost of between 60-92% when compared to standard iFCI calculations, with larger systems experiencing the largest benefit. Tests of MBBSA on two reactions that involve highly correlated systems, the automerization of cyclobutadiene and a Criegee intermediate reaction, show that MBBSA has practical utility for studying realistic chemistries.

Basis sets↗

Abrupt structural transition in exotic molybdenum isotopes unveils an isospin-symmetric island of inversion

Like electrons in atoms, protons and neutrons in nuclei occupy orbitals in a shell structure with energy gaps at magic numbers. Radioactive-beam experiments revealed the disappearance of magic numbers in some neutron-rich isotopes. In these nuclei, configurations involving particles excited across the shell gap gain correlation energy, becoming the ground state. Neutron-rich regions of the nuclear chart that exhibit this property are known as “Islands of Inversion”. Here we present the lifetime measurement of the first 2 + states in 84 Mo (N = Z) and 86 Mo (N = Z + 2) revealing an unexpected sharp structural change between them defining the edge of the region of deformation around 80 Zr. Similarly to the neutron-rich N = 40 Island of Inversion near 64 Cr where cross-shell excitations dominate, we identify this region as an Island of Inversion with symmetrical proton and neutron excitations that we term “Isospin-Symmetric Island of Inversion”. Three-nucleon forces are suggested to drive Mo isotope structural changes.

experimental nuclear physics↗

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter ($α$) was set to recover the exact dissociation limit of H 2 in a minimal basis set. Alternatively $α$ can be viewed as a regularization parameter, where each value of $α$ represents a valid variant of BW-s2, which we denote as BW-s2($α$). In this work, we semi-empirically optimize $α$ for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the ab initio parameterization of BW-s2 and MP2. We demonstrate that the optimal $α$ parameter ($α$ = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2($α$) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2($α$) is iterative $\mathcal{O}$($N^5$), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fermion Pair Dynamics in Open Quantum Systems

Three-body decay is a rare decay mode observed in a handful of unbound rare isotopes. The angular and energy correlations between emitted nucleons are of particular interest, as they provide invaluable information on the interplay between structure and reaction aspects of the nuclear open quantum system. To study the mechanism of two-nucleon emission, we developed a time-dependent approach that allows us to probe emitted nucleons at long times and large distances. We successfully benchmarked the new method against the Green’s function approach and applied it to low-energy two-proton and two-neutron decays. In particular, we studied the interplay between initial-state nucleon-nucleon correlations and final-state interaction. We demonstrated that the time evolution of the two-nucleon wave function is strongly impacted by the diproton/dineutron dynamics and that the correlations between emitted nucleons provide invaluable information on the dinucleon structure in the initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Isotopic mass effects of tritium-fueled high-performance TFTR supershots

An increase in total stored energy correlated with the addition of tritium fuel was observed in supershots during the TFTR DT campaign. This supershot regime had strikingly high, centrally peaked ion and electron temperatures, and the largest neutron emission rates observed in TFTR. This paper presents a study of the causes of this increase in stored energy in supershots. Twenty-six supershots have been recently reanalyzed with the TRANSP plasma analysis code. Early TRANSP simulations did not accurately match the measured magnitude and time evolutions of the neutron emission rates. This mismatch is attributed to neglecting apparent increases of trace amounts of heavy impurities during neutral beam injection. The new TRANSP runs were tuned to accommodate this and match the measured global neutron emission rates. These new runs also had improved fidelity in predicting the time histories and radial dependencies of measured DT neutron emission rates. That in turn adds confidence in the simulated thermal deuterium and tritium density profiles that are needed for calculating the average hydrogenic atomic mass profiles. Six subsets of these supershots had well matched toroidal field B tor , plasma current I p , flux geometry, and total injected neutral beam power. The mix of D and T beam ions was varied for different discharges. The magnitude of the increase of the thermal ion energy W i with added tritium was relatively small, and the total W tot increase is dominated by the increase in fast beam ions with T. Analyses at times before the occurrences of deleterious MHD instabilities yielded scaling of W tot with the volume-average isotopic mass consistent with previous publications. The relative fraction of fast energy ions is expected to be small in practical tokamak reactors. Thus the increase in stored energy W tot observed in TFTR supershots does not appear likely to be significantly helpful for producing useful fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Calculation of electron-ion temperature equilibration rates and friction coefficients in plasmas and liquid metals using quantum molecular dynamics

In this work, we discuss a method to calculate with quantum molecular dynamics simulations the rate of energy exchanges between electrons and ions in two-temperature plasmas, liquid metals, and hot solids. Promising results from this method were recently reported for various materials and physical conditions. Like other ab initio calculations, the approach offers a very useful comparison with the experimental measurements and permits an extension into conditions not covered by the experiments. The energy relaxation rate is related to the friction coefficients felt by individual ions due to their nonadiabatic interactions with electrons. Each coefficient satisfies a Kubo relation given by the time integral of the autocorrelation function of the interaction force between an ion and the electrons. These Kubo relations are evaluated using the output of quantum molecular dynamics calculations in which electrons are treated in the framework of finite-temperature density functional theory. The calculation presents difficulties that are unlike those encountered with the Kubo formulas for the electrical and thermal conductivities. In particular, the widely used Kubo-Greenwood approximation is inapplicable here. Indeed, the friction coefficients and the energy relaxation rate diverge in this approximation since it does not properly account for the electronic screening of electron-ion interactions. The inclusion of screening effects considerably complicates the calculations. We discuss the physically motivated approximations we applied to deal with these complications in order to investigate a widest range of materials and physical conditions. Unlike the standard method used for the electronic conductivities, the Kubo formulas are evaluated directly in the time domain and not in the energy domain, which spares one from needing to introduce an extraneous undetermined numerical parameter to account for the discrete character of the numerical density of states. We highlight interesting properties of the energy relaxation rate not shared by other electronic properties, in particular its self-averaging character. We then present a detailed parametric and convergence study with the numerical parameters, including the system size, the number of bands and k-points, and the physical approximations for the dielectric function and the exchange-correlation energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An electronic origin to the oscillatory segregation behavior in Ni-Cr and other BCC defects in FCC metals

Ni-Cr alloys exhibit oscillatory segregation behaviors near low index surfaces, in which the preferred segregation species changes from Ni in the first layer to Cr in the second layer. In many dilute-alloy systems, this oscillatory pattern is attributed to the elastic release of stresses in the local lattice around the segregating solute or impurity atom. These stresses are mostly thought to originate from mismatches in the atomic size of the solute and host atoms. In Ni-Cr alloys, however, an appreciable mismatch in atomic size is not present, leading to questions about the origins of the oscillatory behavior in this alloy. Therefore, using density functional theory, we have modeled the segregation of a single Cr atom in the and surfaces of FCC Ni, an alloy which exhibits this oscillatory behavior. Using Bader charge analysis, we show that the negative energy correlates directly with the amount of charge on the Cr atom. As Ni atoms strip valence charge from the Cr, the Cr contracts slightly in size. The greatest contraction and highest positive charge for the Cr occurs when it is in the second layer of the surface where the system exhibits the oscillating negative segregation energy. We then find that this behavior persists in other alloy systems (Ag-Nb, Cu-Cr, Pt-Nb, and Pt-V), which exhibit similar atomic radii and electronegativity differences between host and solute to Ni-Cr. These represent alloys in which the host metal exhibits an FCC ground-state structure while the solute metal exhibits a BCC ground-state structure.

36 MATERIALS SCIENCE↗

Predicting a Key Catalyst-Performance Descriptor for Supported Metal Nanoparticles: Metal Chemical Potential

The chemical potential of metal atoms in supported nanoparticles provides a convenient descriptor of their performance as heterogeneous catalysts that captures the effects of particle size, support and alloying. Based on microcalorimetric measurements, the chemical potential is shown in this work to be predictable as a function of monometallic particle diameter and the adhesion energy of the particle to the support, and, for oxide supports, this adhesion energy correlates predictably with metal oxophilicity. These correlations provide predictions of metal chemical potential that can enable catalyst design. They also suggest an improvement in the Gibbs-Thomson relation for free-standing nanoparticles, whether metals or molecular solids.

metal/oxide interface↗

Electrostatics, Charge Transfer, and the Nature of the Halide–Water Hydrogen Bond

Binary halide–water complexes X – (H 2 O) are examined by means of symmetry-adapted perturbation theory, using charge-constrained promolecular reference densities to extract a meaningful charge-transfer component from the induction energy. As is known, the X – (H 2 O) potential energy surface (for X = F, Cl, Br, or I) is characterized by symmetric left and right hydrogen bonds separated by a C 2v -symmetric saddle point, with a tunneling barrier height that is <2 kcal/mol except in the case of F – (H 2 O). Our analysis demonstrates that the charge-transfer energy is correspondingly small (<2 kcal/mol except for X = F), considerably smaller than the electrostatic interaction energy. Nevertheless, charge transfer plays a crucial role determining the conformational preferences of X – (H 2 O) and provides a driving force for the formation of quasi-linear X··· H–O hydrogen bonds. Here, charge-transfer energies correlate well with measured O–H vibrational redshifts for the halide–water complexes and also for OH – (H 2 O) and NO 2 – (H 2 O), providing some indication of a general mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A formalism for extracting track functions from jet measurements

The continued success of the jet substructure program will require widespread use of tracking information to enable increasingly precise measurements of a broader class of observables. The recent reformulation of jet substructure in terms of energy correlators has simplified the incorporation of universal non-perturbative matrix elements, so called “track functions”, in jet substructure calculations. These advances make it timely to understand how these universal non-perturbative functions can be extracted from hadron collider data, which is complicated by the use jet algorithms. In this paper we introduce a new class of jet functions, which we call (semi-inclusive) track jet functions, which describe measurements of the track energy fraction in identified jets. These track jet functions can be matched onto the universal track functions, with perturbatively calculable matching coefficients that incorporate the jet algorithm dependence. We perform this matching, and present phenomenological results for the charged energy fraction in jets at the LHC and EIC/HERA at collinear next-to-leading logarithmic accuracy. We show that higher moments of the charged energy fraction directly exhibit non-linear Lorentzian renormalization group flows, allowing the study of these flows with collider data. Our factorization theorem enables the extraction of universal track functions from jet measurements, opening the door to their use for a precision jet substructure program.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Intense laser interaction with micro-bars

Intense laser fields interact very differently with micrometric rough surfaces than with flat objects. The interaction features high laser energy absorption and increased emission of MeV electrons, ions, and of hard x-rays. In this work, we irradiated isolated, translationally-symmetric objects in the form of micrometric Au bars. The interaction resulted in the emission of two forward-directed electron jets having a small opening angle, a narrow energy spread in the MeV range, and a positive angle to energy correlation. Our numerical simulations show that following ionization, those electrons that are pulled into vacuum near the object’s edge, remain in-phase with the laser pulse for long enough so that the Lorentz force they experience drive them around the object’s edge. After these electrons pass the object, they form attosecond duration bunches and interact with the laser field over large distances in vacuum in confined volumes that trap and accelerate them within a narrow range of momentum. The selectivity in energy of the interaction, its directionality, and the preservation of the attosecond duration of the electron bunches over large distances, offer new means for designing future laser-based light sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron Transfer Theory Elucidates the Hidden Role Played by Triethylamine and Triethanolamine during Photocatalysis

Triethylamine (TEA) and triethanolamine (TEOA) are renowned, in part, for their ability to reductively quench excited states by outer-sphere electron transfer with vast and still growing applications as sacrificial electron donors for photocatalytic systems. Upon amine oxidation, the resulting TEA •+ and TEOA •+ radical cations undergo proton transfer (or hydrogen atom transfer), resulting in the formation of a chemical reductant that has an α-carbon centered radical adjacent to the nitrogen center (TEA • and TEOA • ). In this contribution, we have electrochemically and spectroscopically characterized a set of electron acceptors which, upon accepting an electron, are a series of photocatalysts, [ReCl(R 1 R 2 -bpy)(CO) 3 ] •– , where R 1 and R 2 are electron-donating and electron-withdrawing groups in the 4,4′- and 5,5′-positions on the bipyridyl ligand. Here, we substantiated the formation of the electron donors, TE(O)A • , by spin trapping electron paramagnetic resonance spectroscopy, where TE(O)A • reacts with 2,4,6-tri-tert-butylnitrosobenzene to generate N-centered and O-centered radical adducts. Having established the chemical behaviors of the electron acceptors and donors individually, the electron transfer rate constants were determined across a 1.43 V range in driving force. The redox potential of TEA • was benchmarked to within ±80 mV on an absolute scale in V vs Fc + /Fc in CH 3 CN by using an empirical rate vs free-energy correlation, electron transfer theory, and density functional theory calculations. The equilibrium potentials for TEA • and TEOA • were determined to be −1.98 V and −1.76 V, respectively. Based on the kinetic and thermochemical analysis presented for TEA • and TEOA • , these transient radicals can be broadly considered strong homogeneous chemical reductants within the wider context of photoredox potentials. Thus, this work clarifies a frequently unnoticed secondary function for these sacrificial electron donors during photocatalysis and rationalizes the possibility of a one-photon/two-electron conversion process that is dependent on the free-energy exchange between TE(O)A • and photocatalysts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Linearization of an Electron Beam’s Longitudinal Phase Space Using a Hollow-Channel Plasma

The removal of undesired beam nonlinear energy chirp (time-energy correlation) or linearization of the beam longitudinal phase space (LPS) is crucial for high-brightness linac-based scientific applications, such as x-ray free-electron lasers. In this paper, we propose that a low-density hollow channel plasma can be used as a near-ideal passive linearizer to significantly linearize the beam LPS and, at the same time, preserve the beam emittance. Physically, the passage of the beam through the hollow plasma channel excites a strong quasicosinoidal longitudinal wakefield that acts to mitigate the beam nonlinear energy chirp by superimposing a reverse chirp on the beam. The theoretical analyses and large-scale three-dimensional start-to-end simulations confirm that the beam longitudinal phase space can be almost completely linearized without noticeable beam emittance growth. In conclusion, application of such a near-ideal linearizer may significantly improve the performance of numerous accelerator-based applications.

43 PARTICLE ACCELERATORS↗

A Quantum Mechanical MP2 Study of the Electronic Effect of Nonplanarity on the Carbon Pyramidalization of Fullerene C 60

Among C 60 ’s diverse functionalities, its potential application in CO 2 sequestration has gained increasing interest. However, the processes involved are sensitive to the molecule’s electronic structure, aspects of which remain debated and require greater precision. To address this, we performed structural optimization of fullerene C 60 using the QM MP2/6–31G* method. The nonplanarity of the optimized icosahedron is characterized by two types of dihedral angles: 138° and 143°. The 120 dihedrals of 138° occur between two hexagons intersecting at C–C bonds of 1.42 Å, while the 60 dihedrals of 143° are observed between hexagons and pentagons at C–C bonds of 1.47 Å. NBO analysis reveals less pyramidal sp 1.78 hybridization for carbons at the 1.42 Å bonds and more pyramidal sp 2.13 hybridization for the 1.47 Å bonds. Electrostatic potential charges range from -0.04 a.u. to 0.04 a.u. on the carbon atoms. Second-order perturbation analysis indicates that delocalization interactions in the C–C bonds of 1.42 Å (143.70 kcal/mol) and 1.47 Å (34.98 kcal/mol) are 22% and 38% higher, respectively, than those in benzene. MP2/Def2SVP calculations yield a correlation energy of 13.49 kcal/mol per electron for C 60 , slightly higher than the 11.68 kcal/mol for benzene. However, the results from HOMO-LUMO calculations should be interpreted with caution. This study may assist in the rational design of fullerene C 60 derivatives for CO 2 reduction systems.

36 MATERIALS SCIENCE↗

Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA–ellipticine and buckycatcher–fullerene

In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA–ellipticine intercalation complex, and the buckycatcher–C60 complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller−Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA–ellipticine, and buckycatcher–C60, respectively. Such sizable discrepancy above “chemical accuracy” for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree–Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.

Ballesteros, Francisco (ORCID:0000000325732307)↗