Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION METHODS”

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 OSTI2021

Coexistence and Coupling of Multiple Charge Orderings and Spin States in Hexagonal Ferrite

The coupling between charge and spin orderings in strongly correlated systems plays a crucial role in fundamental physics and device applications. As a candidate of multiferroic materials, LuFe 2 O 4 with a nominal Fe 2.5+ valence state has the potential for strong charge–spin interactions; however, these interactions have not been fully understood until now. Here, combining complementary characterization methods with theoretical calculations, two types of charge orderings with distinct magnetic properties are revealed. The ground states of LuFe 2 O 4 are decided by the parallel/antiparallel coupling of both charge and spin orderings in the adjacent FeO double layers. Whereas the ferroelectric charge ordering remains ferrimagnetic below 230 K, the antiferroelectric ordering undergoes antiferromagnetic–ferrimagnetic–paramagnetic transitions from 2 K to room temperature. This study demonstrates the unique aspects of strong spin-charge coupling within LuFe 2 O 4 . Our results shed light on the coexistence and competing nature of orderings in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Low Thermal Conductivity in Heteroanionic Materials with Layers of Homoleptic Polyhedra

Although BiAgOSe, an analogue of a well-studied thermoelectric material BiCuOSe, is thermodynamically stable, its synthesis is complicated by the low driving force of formation from the stable binary and ternary intermediates. Here we have developed a “subtraction strategy” to suppress byproducts and produce pure phase BiAgOSe using hydrothermal methods. Electronic structure calculations and optical characterization show that BiAgOSe is an indirect bandgap semiconductor with a bandgap of 0.95 eV. The prepared sample exhibits lower lattice thermal conductivities (0.61 W·m –1 ·K –1 at room temperature and 0.35 W·m –1 ·K –1 at 650 K) than BiCuOSe. Further, lattice dynamical simulations and variable temperature diffraction measurements demonstrate that the low lattice thermal conductivity arises from both the low sound velocity and high phonon–phonon scattering rates in BiAgOSe. These in turn result primarily from the soft Ag–Se bonds in the edge-sharing AgSe 4 tetrahedra and large sublattice mismatch between the quasi-two-dimensional [Bi 2 O 2 ] 2+ and [Ag 2 Se 2 ] 2– layers. These results highlight the advantages of manipulating the chemistry of homoleptic polyhedra in heteroanionic compounds for electronic structure and phonon transport control.

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

Quantifying the Influences of Epoxide Binding in Epoxide/CO 2 Ring Opening Copolymerization Catalysis

Understanding and predicting the effect of epoxide structure on the rate of polymerization in epoxide/CO 2 ring opening copolymerization catalysis is a long-standing challenge. Here, a known highly active Co(III)K(I) catalyst is used to investigate the influences of six different epoxides' binding strengths on their rates of copolymerization. Since calculations and experiments indicate that studying the catalytically relevant Co(III)−epoxide adduct directly is experimentally challenging, epoxide−catalyst binding interactions are quantified using a Co(II)K(I) complex to model the key catalytic intermediate. Epoxide−catalyst coordination is investigated using UV−vis spectroscopy titrations which provide fast and effective determination of association or binding constants. The epoxide−catalyst equilibrium constants show a clear exponential correlation with copolymerization rates and a new catalyst performance linear free energy relationship is revealed. Epoxides exhibiting stronger catalyst binding constants show higher copolymerization rates. The structure−activity correlation is consistent with the polymerization kinetics, mechanism and DFT calculations. Both the methods to investigate epoxide−catalyst coordination and the linear free energy relationship are shown to apply to the series of six epoxides and a second Co(III)K(I) catalyst. These structure− performance relationships are likely applicable to other transition metal catalysts and should expedite future epoxide and catalyst selection to make useful poly(carbonate) materials.

catalysts↗
DOE OSTI2021

ESI-MS Identification of the Cationic Phosphine-Ligated Gold Clusters Au1-Au22: Insight into the Gold-Ligand Ratio and Abundance of Larger Clusters

Triphenylphosphine (PPh3)-ligated gold clusters offer promising potential applications due to their relative ease of synthesis and usefulness in forming advanced cluster architectures. While previous studies reported cationic PPh3-ligated gold clusters with core sizes of Au1 - Au4, Au6 - Au11, and Au¬13 - Au14, there has not been definitive identification by mass spectrometry of larger clusters in the Au12 - Au25¬ range. Herein, we survey a polydisperse solution of cationic PPh3-ligated gold clusters using high mass-resolution (M/?M = 60,000) electrospray ionization mass spectrometry (ESI-MS). To improve the sensitivity and mass resolution of larger clusters for unambiguous identification, we increased the number of scan averages and reduced the range of mass collection windows to 200 m/z, thereby mitigating potential mass and ion abundance bias resulting from smaller “building block” gold clusters and other solution components present in higher abundance. In addition to the previously reported clusters, we identified several new species including Au5(PPh3)5+, Au12(PPh3)9HCl2+, Au15(PPh3)9Cl2+, Au16(PPh3)10Cl22+, Au17(PPh3)113+, Au18(PPh3)102+, Au19(PPh3)10Cl2+, Au20(PPh3)12H33+, Au21(PPh3)10Cl2+, and Au22(PPh3)10Cl22+, indicating that a full range of clusters between Au1 - Au22 may be observed in a single polydisperse solution. Considering all of the observed clusters, our findings provide evidence that the “magic number” icosahedral Au13 may be the transition point in cluster growth between smaller clusters, exhibiting a 1:1 gold-to-ligand ratio, and larger clusters, wherein subsequent gold atoms are added to the core without an equal number of accompanying ligands. Our method demonstrates that reducing the range of m/z collection windows and increasing the number of scan averages can improve instrument sensitivity for cationic gold clusters and enable a more complete survey of polydisperse solutions, thereby providing new insights to guide and validate the results of other characterization methods and theoretical calculations. This work was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. MH acknowledges support from the DOE Science Undergraduate Laboratory Internship (SULI) program. HH acknowledges support from the DOE Office of Workforce Development for Teachers and Scientist (WDTS) under the Visiting Faculty Program (VFP). This work was performed using EMSL, a national scientific user facility sponsored by the DOE's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is a multiprogram national laboratory operated for DOE by Battelle.

Hewitt, Michael↗
DOE OSTI2023

Electronic structure and energetics of a heterodimeric BChl g ′/Chl a ′ special pair generated by exposure of Heliomicrobium modesticaldum to dioxygen

Heliobacteria are anoxygenic phototrophs that have a Type I homodimeric reaction center containing bacteriochlorophyll g (BChl g). Previous experimental studies have shown that in the presence of light and dioxygen, BChl g is converted into 8 1 -OH-chlorophyll a F (hereafter Chl a F ), with an accompanying loss of light-driven charge separation. These studies suggest that the reaction center only loses the ability to transfer electrons once both BChl g′ molecules of the P 800 special pair have been converted to Chl a F ′. Here, the present work confirms that the partially converted BChl g′/Chl a F ′ special pair remains functional in samples exposed to dioxygen by demonstrating its presence using hyperfine couplings obtained from Q-band 1 H ENDOR, 2D 14 N HYSCORE and DFT methods. The DFT calculations of the BChl g′/BChl g′ homodimeric primary donor, which are based on the recently published X-ray crystal structure, predict that the unpaired electron spin is equally delocalized over both BChl g′ molecules and provide an excellent match to the experimental hyperfine couplings of the anaerobic samples. Exposure to dioxygen leads to substantial changes in the hyperfine interactions, indicative of greater localization of the unpaired electron spin. The measured hyperfine couplings are reproduced in the DFT calculations by replacing one of the BChl g′ molecules of the primary donor with a Chl a F ′ molecule. The calculations reveal that the spin density becomes localized on BChl g′ in the heterodimeric primary donor. Time-dependent DFT calculations demonstrate that conversion of either or both of the accessory BChl g molecules and/or one of the BChl g′ molecules of P800 to Chl a F ′ results in minor effects on the energy of the charge-separated states. In contrast, if both of the BChl g′ molecules of P 800 are converted a large increase in the energy of the charge-separated state occurs. This suggests that the reaction center remains functional when only one half of the dimer is converted, however, conversion of both halves of the P 800 dimer leads to loss of function

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

Molecular dynamics simulations of inelastic x-ray scattering from shocked copper

By taking the spatial and temporal Fourier transforms of the coordinates of the atoms in molecular dynamics simulations conducted using an embedded-atom-method potential, we calculate the inelastic scattering of x rays from copper single crystals shocked along [001] to pressures of up to 70 GPa. Above the Hugoniot elastic limit, we find that the copious stacking faults generated at the shock front introduce strong quasi-elastic scattering (QES) that competes with the inelastic scattering signal, which remains discernible within the first Brillouin zone; for specific directions in reciprocal space outside the first zone, the QES dominates the inelastic signal overwhelmingly. The synthetic scattering spectra we generate from our Fourier transforms suggest that energy resolutions of order 10 meV would be required to distinguish inelastic from quasi-elastic scattering within the first Brillouin zone of shock-loaded copper. We further note that high-resolution inelastic scattering also affords the possibility of directly measuring particle velocities via the Doppler shift. These simulations are of relevance to future planned inelastic scattering experiments at x-ray Free Electron Laser facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid–solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Self-interaction error overbinds water clusters but cancels in structural energy differences

Here, we gauge the importance of self-interaction errors in density functional approximations (DFAs) for the case of water clusters. To this end, we used the Fermi–Löwdin orbital self-interaction correction method (FLOSIC) to calculate the binding energy of clusters of up to eight water molecules. Three representative DFAs of the local, generalized gradient, and metageneralized gradient families [i.e., local density approximation (LDA), Perdew–Burke–Ernzerhof (PBE), and strongly constrained and appropriately normed (SCAN)] were used. We find that the overbinding of the water clusters in these approximations is not a density-driven error. We show that, while removing self-interaction error does not alter the energetic ordering of the different water isomers with respect to the uncorrected DFAs, the resulting binding energies are corrected toward accurate reference values from higher-level calculations. In particular, self-interaction–corrected SCAN not only retains the correct energetic ordering for water hexamers but also reduces the mean error in the hexamer binding energies to less than 14 meV/ H 2 O from about 42 meV/ H 2 O for SCAN. By decomposing the total binding energy into many-body components, we find that large errors in the two-body interaction in SCAN are significantly reduced by self-interaction corrections. Higher-order many-body errors are small in both SCAN and self-interaction–corrected SCAN. These findings suggest that orbital-by-orbital removal of self-interaction combined with a proper DFA can lead to improved descriptions of water complexes.

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

Influence of nuclear quantum effects on the electronic properties of amorphous carbon

We carry out quantum simulations to study the physical properties of diamond-like amorphous carbon by coupling first-principles molecular dynamics with a quantum thermostat, and we analyze multiple samples representative of different defective sites present in the disordered network. We show that quantum vibronic coupling is critical in determining the electronic properties of the system, in particular its electronic and mobility gaps, while it has a moderate influence on the structural properties. We find that despite localized electronic states near the Fermi level, the quantum nature of the nuclear motion leads to a renormalization of the electronic gap surprisingly similar to that found in crystalline diamond. We also discuss the notable influence of nuclear quantum effects on band-like and variable-hopping mechanisms contributing to electrical conduction. Our calculations indicate that methods often used to evaluate electron–phonon coupling in ordered solids are inaccurate to study the electronic and transport properties of amorphous semiconductors composed of light atoms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Three-Dimensional Full-Core BEAVRS Using OpenMOC with Transport Equivalence

Using an optimized implementation of the three-dimensional (3D) method of characteristics for neutron transport, along with a novel equivalence method for transport calculations that was designed to correct self-shielding errors from neglecting the angular dependence of resonant group absorption, a 3D full-core light water reactor hybrid stochastic-deterministic eigenvalue calculation was achieved. This paper presents the optimizations developed and compares the transport solutions obtained. For the statepoint, run times near 10 000 CPU hours are achieved—improving on previous works by an order of magnitude—with near 1% error on pin fission to 238 U capture ratios and a few dozen pcms on the eigenvalue.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Observation of d(t,n) α Neutrons Following d(d,p)t Reactions in a Deuterium Gas Cell: An Attempt to Repeat Ruhlig’s 1938 Observation of Secondary Reactions

In order to benchmark methods used to calculate reaction-in-flight fusion reactions in inertial confinement fusion and address issues related to the first claimed observation of d(t,n)α reactions in 1938, secondary d(t,n)α reactions have been observed following d(d,p)t reactions in deuterium gas. A pulsed 200-nA, 2.2-MeV deuterium beam from the Triangle Universities Nuclear Laboratory FN tandem accelerator was injected into a cylindrical multiatmosphere deuterium gas target. The incident beam traversed along the target cylinder’s 3-cm symmetry axis after its passage through a Havar entrance foil. Two different Havar foil thicknesses were used to obtain 1.5- and 0.6-MeV deuteron beams entering the deuterium cell. The cylinder’s radius was 2 cm to allow for d(d,p)t tritons emitted perpendicular to the beam to range out in the deuterium gas. The neutron emission from the cell was observed via its time of flight to a liquid scintillator placed at various angles to the beam direction, at a distance of 243 cm. Pulse-shape-discrimination techniques were used to separate neutron and gamma-ray signals seen in the liquid scintillator. The observed probability of ~2 × 10 –4 for inducing secondary d(t,n)α fusion in the gas cell per d(d,p)t reaction is consistent with theoretical expectations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Disorder and Itinerant Magnetism in Full Heusler Pd 2 TiIn

We report electronic and magnetic properties of full Heusler Pd 2 TiIn based on first principles calculations. This compound has been variously characterized as magnetic or non-magnetic. We use first principles calculations with accurate methods to reexamine this issue. We find that ideal ordered Heusler Pd 2 TiIn remains non-magnetic, in accord with prior work. Furthermore, we do find that it is possible to explain the magnetism seen in experiments through disorder and in particular we find that site disorder can lead to moment formation in this compound. In addition, we find an alternative low energy cubic crystal structure, which will be of interest to explore experimentally.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First on-line commissioning experiments at the St. Benedict facility

Nuclear beta decays provide an excellent probe of fundamental symmetries due to their mediation by the weak interaction. In particular, precise measurements of these decays provide constraints on the unitarity of the Cabbibo-Kobayashi-Maskawa (CKM) quark-mixing matrix. While superallowed pure Fermi decays currently set the most precise limits, the alternative suite of superallowed mixed mirror decays has been ill-studied. These nuclei can provide an important consistency check of calculation and measurement methods employed for the pure Fermi decays, more critically needed now in the wake of a 2.4σ deviation from unitarity of the CKM matrix. In order to remedy the gap in data for mirror decays, the Superallowed Transition Beta-Neutrino Decay Ion Coincidence Trap (St. Benedict) facility is being commissioned at the University of Notre Dame’s Nuclear Science Laboratory (NSL). In this paper, we present first results of the commissioning of the St. Benedict facility on-line at the TwinSol radioactive beam facility. The results of initial commissioning experiments involving the St. Benedict gas catcher, RF carpet, RFQ ion guide and RFQ cooler-buncher will be presented.

beta decay↗
DOE OSTI2023

Cosmic variance and the inhomogeneous UV luminosity function of galaxies during reionization

When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from z ≃ 6–15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same z. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by z = 6 (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, ‘rolling-over’ at M UV ≳ -17. CoDa II’s UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at z = 6–8, including the faint end, except for the faintest data point at z = 6, based on one galaxy at M UV = -12.5. According to CoDa II, the probability of observing the latter is $\sim 5~{{\ \rm per\ cent}}$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. Here, we find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2–4, primarily by underestimating the variance in halo abundance when the sample volume is small.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Linear polarization fractions of Fulcher- α fluorescence in electron collisions with H 2

We apply the molecular convergent close-coupling method to the calculation of linear polarization fractions for Fulcher-band fluorescence following electron-impact excitation of the H 2 d 3 Π u state. Furthermore, the results exhibit the opposite threshold behavior compared to the only previous calculations, but are in agreement with the most recent measurements for the Q(1), R(1), and Q(3) transitions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Checking nonflow assumptions and results via PHENIX published correlations in p+p, p+Au, d+Au, and 3 He+Au at √ s NN=200 GeV

Recently the PHENIX Collaboration has made available two-particle correlation Fourier coefficients for multiple detector combinations in minimum bias p+p and 0–5% central p+Au, d+Au, and 3 He+Au collisions at √ s NN = 200 GeV [Phys. Rev. C 105, 024901 (2022)]. Using these coefficients for three sets of two-particle correlations, azimuthal anisotropy coefficients v 2 and v 3 are extracted for midrapidity charged hadrons as a function of transverse momentum. In this paper, we use the available coefficients to explore various nonflow hypotheses as well as to compare the results with theoretical model calculations. The nonflow methods fail basic closure tests with ampt and pythia/angantyr, particularly when including correlations with particles in the low multiplicity light-projectile going direction. In data, the nonflow adjusted v 2 results are modestly lower in p+Au and the adjusted v 3 results are more significantly higher in p+Au and d+Au. However, the resulting higher values for the ratio v 3 /v 2 in p+Au at RHIC compared to p+Pb at the LHC is additional evidence for a significant overcorrection. Incorporating these additional checks, the conclusion that these flow coefficients are dominated by initial geometry coupled with final-state interactions (e.g., hydrodynamic expansion of quark-gluon plasma) remains true, and explanations based solely on initial-state glasma are ruled out. The detailed balance between intrinsic and fluctuation-driven geometry and the exact role of weakly versus strongly coupled prehydrodynamic evolution remains an open question for triangular flow, requiring further theoretical and experimental investigation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Uncertainty quantification for neutrino opacities in core-collapse supernovae and neutron star mergers

We perform an extensive study of the correlations between the neutrino-nucleon inverse mean free paths (IMFPs) and the underlying equation of states (EoSs). Strong interaction uncertainties in the neutrino mean free path are investigated in different density regimes. The nucleon effective mass, the nucleon chemical potentials, and the residual interactions in the medium play an important role in determining neutrino-nucleon interactions in a density-dependent manner. Here we study how the above quantities are constrained by an EoS consistent with (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. We then study the uncertainties of both the charged current and the neutral current neutrino-nucleon inverse mean free paths due to the variation of these quantities, using the Hartree-Fock+random phase approximation method. Finally, we calculate the Pearson correlation coefficients between (i) the EoS-based quantities and the EoS-based quantities; (ii) the EoS-based quantities and the IMFPs; (iii) the IMFPs and the IMFPs. We find a strong impact of residual interactions on neutrino opacity in the spin and spin-isospin channels, which are not well constrained by current nuclear modelings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Predicting the neutrinoless double-β-decay matrix element of 136 Xe using a statistical approach

Calculation of the nuclear matrix elements (NMEs) for double-β decay is of paramount importance for guiding experiments and for analyzing and interpreting the experimental data, especially for the search of the neutrinoless double β decay mode (0v⁢ββ). However, there are currently still large differences between the NME values calculated by different methods, hence a quantification of their uncertainties is very much required. Here, in this paper, we propose a statistical analysis of 0⁢v⁢ββ NME for the 136 Xe isotope, based on the interacting shell model, but using three independent effective Hamiltonians, emphasizing the range of the NMEs' most probable values and its correlations with observables that can be obtained from the existing nuclear data. Consequently, we propose a common probability distribution function for the 0⁢vββ NME, which has a range of (1.55–2.65) at 90% confidence level, with a mean value of 1.99 and a standard deviation of 0.37.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗