Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pair distribution functions”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effects of anisotropy and disorder on the superconducting properties of niobium

We report results for the superconducting transition temperature and anisotropic energy gap for pure niobium based on Eliashberg’s equations and electron and phonon band structures computed from density functional theory. The electronic band structure is used to construct the Fermi surface and calculate the Fermi velocity at each point on the Fermi surface. The phonon bands are in excellent agreement with inelastic neutron scattering data. The corresponding phonon density of states and electron–phonon coupling define the electron–phonon spectral function, α 2 F(p, p'; ω), and the corresponding electron–phonon pairing interaction, which is the basis for computing the superconducting properties. The electron–phonon spectral function is in good agreement with existing tunneling spectroscopy data except for the spectral weight of the longitudinal phonon peak at ℏω LO = 23 meV. We obtain an electron–phonon coupling constant of λ = 1.057, renormalized Coulomb interaction μ* = 0.218, and transition temperature T c = 9.33 K. The corresponding strong-coupling gap at T = 0 is modestly enhanced, Δ 0 = 1.55 meV, compared to the weak-coupling BCS value $Δ^{\text {wc}}_{0}$ = 1.78 k B T c = 1.43 meV. The superconducting gap function exhibits substantial anisotropy on the Fermi surface. We analyze the distribution of gap anisotropy and compute the suppression of the superconducting transition temperature using a self-consistent T-matrix theory for quasiparticle-impurity scattering to describe niobium doped with non-magnetic impurities. We compare these results with experimental results on niobium SRF cavities doped with nitrogen impurities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Observation of $t\overline{t}$ production in the lepton+jets and dilepton channels in p+Pb collisions at $\sqrt{s_{\textrm{NN}}}$ = 8.16 TeV with the ATLAS detector

This paper reports the observation of top-quark pair production in proton-lead collisions in the ATLAS experiment at the Large Hadron Collider. The measurement is performed using 165 nb -1 of p+Pb data collected at $\sqrt{^sNN}$ = 8.16 TeV in 2016. Events are categorised in two analysis channels, consisting of either events with exactly one lepton (electron or muon) and at least four jets, or events with two opposite-charge leptons and at least two jets. In both channels at least one b-tagged jet is also required. Top-quark pair production is observed with a significance over five standard deviations in each channel. The top-quark pair production cross-section is measured to be σ $t\overline{t}$ =58.1 ± 2.0(stat.)$^{+4.8}_{-4.4}$(syst.) nb, with a total uncertainty of 9%. In addition, the nuclear modification factor is measured to be R pA = 1.090 ± 0.039(stat.)$^{+0.094}_{-0.087}$(syst.). The measurements are found to be in good agreement with theory predictions involving nuclear parton distribution functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

CO 2 Capture Characteristics of Hyperbranched Poly(alkylene imine): A Molecular Dynamics Simulation Approach

This study explores the CO 2 capture characteristics of hyperbranched poly(ethylenimine) (HB-PEI) and poly- (propyleneimine) (HB-PPI) through molecular dynamics simulations using density functional theory-calibrated force fields. Key features such as density, free volume, glass transition temperature, CO 2 /H 2 O distribution, and molecular diffusion are systematically investigated to elucidate structure−function relationships under dry and hydrated conditions. HB-PEI demonstrates a slightly higher density and lower free volume compared to HB-PPI yet shows superior CO 2 capture due to the high amine concentration. Glass transition analysis indicates a higher thermal mobility in HBPEI, enhancing the CO 2 diffusivity. Pair correlation and coordination analyses confirm a stronger affinity of CO 2 with primary and secondary amines, particularly in hydrated environments where water competes with CO 2 for binding sites. Despite its more compact structure, HB-PEI outperformed HB-PPI in CO 2 and H 2 O transport, as confirmed by higher diffusion coefficients across all hydration levels. These findings highlight a critical balance among polymer architecture, amine accessibility, and hydration in designing next-generation solid amine sorbents for efficient direct air capture applications.

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

A universal and constant rate of gene content change traces pangenome flux to LUCA

Abstract Prokaryotic genomes constantly undergo gene flux via lateral gene transfer, generating a pangenome structure consisting of a conserved core genome surrounded by a more variable accessory genome shell. Over time, flux generates change in genome content. Here, we measure and compare the rate of genome flux for 5655 prokaryotic genomes as a function of amino acid sequence divergence in 36 universally distributed proteins of the informational core (IC). We find a clock of gene content change. The long-term average rate of gene content flux is remarkably constant across all higher prokaryotic taxa sampled, whereby the size of the accessory genome—the proportion of the genome harboring gene content difference for genome pairs—varies across taxa. The proportion of species-level accessory genes per genome, varies from 0% (Chlamydia) to 30%–33% (Alphaproteobacteria, Gammaproteobacteria, and Clostridia). A clock-like rate of gene content change across all prokaryotic taxa sampled suggest that pangenome structure is a general feature of prokaryotic genomes and that it has been in existence since the divergence of bacteria and archaea.

Microbiology↗
DOE OSTI2021

Multimodal Bayesian registration of noisy functions using Hamiltonian Monte Carlo

Functional data registration is a necessary processing step for many applications. The observed data can be inherently noisy, often due to measurement error or natural process uncertainty; which most functional alignment methods cannot handle. A pair of functions can also have multiple optimal alignment solutions, which is not addressed in current literature. In this paper, a flexible Bayesian approach to functional alignment is presented, which appropriately accounts for noise in the data without any pre-smoothing required. Additionally, by running parallel MCMC chains, the method can account for multiple optimal alignments via the multi-modal posterior distribution of the warping functions. To most efficiently sample the warping functions, the approach relies on a modification of the standard Hamiltonian Monte Carlo to be well-defined on the infinite-dimensional Hilbert space. In this work, this flexible Bayesian alignment method is applied to both simulated data and real data sets to show its efficiency in handling noisy functions and successfully accounting for multiple optimal alignments in the posterior; characterizing the uncertainty surrounding the warping functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Designing a quantum-accurate machine-learning potential to enable large-scale simulations of deuterium under shock

Large-scale molecular dynamics of deuterium under shock can elucidate kinetic processes vital to the target design in inertial confinement fusion and high-energy-density experiments. However, modeling the complex evolution of this material from an insulating molecular state at ambient pressure to an ionized, atomic fluid under strong shock is beyond the capability of simple pair and even bond order potentials. We thus train a quantum-accurate and broadly transferable machine-learning interatomic potential for deuterium using the Chebyshev Interaction Model for Efficient Simulations framework. We show that due to an improved description of the molecular-to-atomic transition, our model is able to better reproduce the ab initio equation of state, radial distribution functions, and principal Hugoniot than bond order potentials. This represents an important step toward large-scale quantum-accurate and nonequilibrium simulations of complicated systems under dynamic changes including phase transitions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fission Product Evaluation (Quarterly Report)

The LLNL team has developed a theoretical framework to compute the fission fragment charge and mass distributions of odd-mass nuclei, which can be formed in neutron-induced reactions on even-even nuclei. For example, n+ 238 U leads to the compound nucleus 239 U which has 147 neutrons. Odd-mass and odd-odd nuclei are spin-polarized: the spin of the nucleus is not equal to 0 as in even-even nuclei (a consequence of pairing correlations) but is in fact set by a single particle. In practice, that single particle can occupy different quantum-mechanical configurations, each with different spin. This affects the deformation properties of the nucleus and, as a result, the distribution of fission fragments. These distributions require two inputs for the calculation: a potential energy surface (= the total energy of the nucleus as a function of the deformation of the nuclear shape, often parametrized with 2 degrees of freedom) and the collective inertia tensor associated with each point of the PES. In an odd-mass or odd-odd nucleus, each quantity depends on the configuration of the single particle. We computed of the order of 20 to 30 different such configurations to extract the PES for spins K=1/2, 3/2, 5/2, 7/2, and 9/2 of 239 U. We then computed the primary yields for each such PES, as shown in Fig. 1 below. These yields were weighed with the spin distribution associated with the reaction 238 U+n and the resulting distribution used as input to FREYA to compute the independent yields shown in Fig. 2. These yields are then used as inputs to the FIER code to extract the cumulative yields and compare with LLNL-LANL measurements at TUNL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Statistical characterization of twin transmission across grain boundaries in magnesium

Here, a detailed statistical analysis of twin transmission (TT) across grain boundaries is performed in rolled, commercial purity magnesium compressed along the rolling direction. EBSD images are acquired from two different cuts: a section containing the rolling (RD) and normal (ND) directions; and a section at 35° to the ND which contains transverse direction (TD). An automated twinning analysis software, METIS, is used to obtain the statistical correlations between deformation twins and other microstructural features by analyzing EBSD microstructures comprising thousands of grains and twins. This detailed statistical analysis reveals that the TT propensity is sensitive to the grain boundary (GB) misorientation angle but not to the GB misorientation axis. Specifically, TT propensity decreases with increasing GB misorientation angle; however, the decreasing trend is not monotonic. Further, the detailed analysis of TT events combining macroscopic Schmid factor and the geometric measures (m’: accounts for the alignment of both twin plane normal and shear directions; m”: accounts for the alignment of only twin shear directions; and m”’: accounts for the alignment of only twin plane normal) helps in understanding and identifying the process of twin-pair formation, i.e., co-nucleation versus transmission, also to identify the role of local stresses at GBs induced by inter-grain interaction versus twinning shear transformation processes. In addition, Bayesian inference is used to draw statistically meaningful conclusions as to the likelihood of transmission given misorientation angles. Lastly, twin chain frequency as a function of chain length is predicted using the probability distribution obtained from the Bayesian inference and compared with the actual data from the EBSD microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Information content of the differences in the charge radii of mirror nuclei

Differences in the charge radii of mirror nuclei have been recently suggested to contain information on the slope of the symmetry energy L. To test this hypothesis, we perform statistical correlation analysis using quantified energy density functionals that are consistent with our previous knowledge on global nuclear observables such as binding energies and charge radii. We conclude that the difference in charge radii between a mirror pair, $ΔR$$^{mir}_{ch}$, is an inferior isovector indicator compared to other observables, such at the neutron skin or electric dipole polarizability α D . In particular, this quantity correlates poorly with both the neutron skin and L. Here, we demonstrate that $ΔR$$^{mir}_{ch}$ is influenced by pairing correlations in the presence of low-lying proton continuum in the proton-rich mirror-partner nucleus. Considering the large theoretical uncertainties on $ΔR$$^{mir}_{ch}$, we conclude that the precise data on mirror charge radii cannot provide a stringent constraint on L.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Solid earth as a recycling systems and the lateral growth of Precambrian North America

If plotted on mass vs time diagrams, geologic entities (for example, continental and oceanic crust, sediments, and mineral resources) display an exponential (power law) relationship, with entity per unit time increasing toward the present. This relationship is consistent with the concept of recycling and can be simulated mathematically. The approach is based on the plate tectonic theory and considers area-age or mass-age distributions of crystalline basement and sediments for major global tectonic realms. Each tectonic realm is characterized by a specific lifespan, which is an inverse function of its recycling rate. The estimated average half-area of half-mass ages are given. The corresponding parameters for continental crust are 690 Ma for K/Ar, and approximately 1200 Ma for Rb/St and U-Th/Pb dating pairs. Tectonic diversity preserved in the geologic record is therefore a function of time, with oceanic tectonic realms, because of their rapid recycling, underrepresented in the rocks older than approximately 300 Ma. The Sm/Nd isotopic systematic of sediments suggest that, for a near steady-state post-Archean sedimentary mass, recycling is approximately 90 + or - 5 percent cannibalistic. This yields an estimated upper limit on crust-mantle exchange via sediment subduction of approximately 1.1 + or - 0.5 x 10 g a(sup -1) considerably less than demanded by isotopic constraints. The discrepancy may indicate the existence of additional loci, such as orogenic belts, for significant crust-mantle interaction.

Veizer, Jan↗
DOE OSTI2024

Coherent γ *-nucleus scattering and coherent nuclear states

In the context of a McLerran-Venugopalan (MV) model for a large nucleus, coherent scattering of a virtual photon on that nucleus is evaluated in the A − = 0 gauge, the gauge appropriate for the target nucleus. The evaluation of the scattering in A − = 0 gauge is very intricate compared to the usual A + = 0 gauge evaluation natural for the scattering process, but has the advantage of directly giving the scattering in terms of a partonic description of the nucleus. In the limit where a tagged forward jet puts the dipole-nucleus scattering in the saturation regime the coherent reactions are equal to the inelastic reactions. In terms of the nuclear wave function the coherent reactions come from color singlet and zero total transverse momentum quark-antiquark pairs in the wave function and in the saturation regime the nuclear wave function is a coherent state for these pairs. In the saturation region half of all quarks (or antiquarks) come from zero momentum and color charge pairs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Hydrogen Localization and Cluster Formation in α-Zr from First-Principles Investigations

We report alpha-zirconium (α-Zr) is a hydrogen atom host material used in nuclear energy and semiconducting applications due to Zr's excellent properties such as corrosion resistance, low thermal neutrons cross-section, and thermodynamically favorable interstitial sequestering of H atoms. Herein we use density functional theory simulations to systematically analyze the initial stages of H insertion in the Zr matrix, and determine optimal hydrogen localization, distribution, and clustering tendencies from a thermodynamics and kinetics perspective. Our results suggest that H occupation of tetrahedral sites is strongly favored, consistent with previous experimental and theoretical studies. Further, a cooperative effect in H localization is shown where H atoms favor clustering. The H pair is found to be the most stable by occupying tetrahedral sites with pairs at second nearest-neighbor positions (TTn2). Among the four H atom clusters, the most favorable configuration has two TTn2 pairs forming "box " or "zig-zag" shaped-cluster with the pairs separated by 3.03 or 3.25 Å, respectively. Nudged elastic band simulations suggest that the formation of TTn2 is facile and is limited by the small diffusion barrier of the monomer of less than 0.1 eV. Further, the formation routes for the 4H atom clusters are also presented. The current study can serve as a foundation to systematically model H atom clustering in metals and alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Search for the critical point of strongly-interacting matter in $^{40}$Ar + $^{45}$Sc collisions at 150A Ge V /c using scaled factorial moments of protons

The critical point of dense, strongly interacting matter is searched for at the CERN SPS in $^{40}$Ar + $^{45}$Sc collisions at 150A Ge V /c. The dependence of second-order scaled factorial moments of proton multiplicity distribution on the number of subdivisions of transverse momentum space is measured. The intermittency analysis is performed using both transverse momentum and cumulative transverse momentum. For the first time, statistically independent data sets are used for each subdivision number. The obtained results do not indicate any statistically significant intermittency pattern. An upper limit on the fraction of correlated proton pairs and the power of the correlation function is obtained based on a comparison with the Power-law Model developed for this purpose.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Anisotropic effective redshift and evolving clustering amplitude

A typical galaxy survey geometry results in galaxy pairs of different separation and angle to the line-of-sight having different distributions in redshift and consequently a different effective redshift. However, clustering measurements are analysed assuming that the clustering is representative of that at a single effective redshift. Here, we investigate the impact of variations in the galaxy-pair effective redshift on the large-scale clustering measured in galaxy surveys. We find that galaxy surveys spanning a large redshift range have different effective redshifts as a function of both pair separation and angle. Furthermore, when considering tracers whose clustering amplitude evolves strongly with redshift, this combination can result in an additional scale-dependent clustering anisotropy. We demonstrate the size of this effect on the eBOSS DR16 Quasar sample and show that, while the impact on monopole is negligible, neglecting this effect can result in a tilt of ~ 4% and ~40% in quadrupole and hexadecapole, respectively, on scales ~1 h -1 Gpc. We discuss strategies to mitigate this effect when making

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗
NASA NTRS1994

An outer gap model of high-energy emission from rotation-powered pulsars

We describe a refined calculation of high-energy emission from rotation-powered pulsars based on the outer gap model of Cheng, Ho, & Ruderman (1986 a, b). In this calculation, vacuum gaps form in regions near the speed-of-light cylinder of the pulsar magnetosphere along the boundary between the closed and open field line zones. We have improved upon previous efforts to model the spectra from these pulsars (e.g., Cheng et al. 1986b; Ho 1989) by following the variation in particle production and radiation properties with position in the outer gap. Curvature, synchotron, and inverse-Compton scattering fluxes vary significantly over the gap, and their interactions via photon-photon pair production build up the radiating charge populations at varying rates. We have also incorporated an approximate treatment of the transport particle and photon fluxes between gap emission zones. These effects, along with improved computations of the particle and photon distribution, provide very important modifications of the model gamma-ray flux. In paticular, we attempt to make specific predictions of pulse profile shapes and spectral variations as a function of pulse phase and suggest further extensions to the model which may provide accurate computations of the observed high-energy emissions.

Chiang, James↗
DOE OSTI2022

Next-to-leading power endpoint factorization and resummation for off-diagonal “gluon” thrust

The lack of convergence of the convolution integrals appearing in next-to-leading-power (NLP) factorization theorems prevents the applications of existing methods to resum power-suppressed large logarithmic corrections in collider physics. We consider thrust distribution in the two-jet region for the flavour-nonsinglet off-diagonal contribution, where a gluon-initiated jet recoils against a quark-antiquark pair, which is power-suppressed. With the help of operatorial endpoint factorization conditions, we obtain a factorization formula, where the individual terms are free from endpoint divergences in convolutions and can be expressed in terms of renormalized hard, soft and collinear functions in four dimensions. This allows us to perform the first resummation of the endpoint-divergent SCETI observables at the leading logarithmic accuracy using exclusively renormalization-group methods. The presented approach relies on universal properties of the soft and collinear limits and may serve as a paradigm for the systematic NLP resummation for other 1 → 2 and 2 → 1 collider physics processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗