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

Failure of molecular dynamics to provide appropriate structures for quantum mechanical description of the aqueous chloride ion charge-transfer-to-solvent ultraviolet spectrum

The lowest band in the charge-transfer-to-solvent ultraviolet absorption spectrum of aqueous chloride ion is studied by experiment and computation. Interestingly, the experiments indicate that at concentrations up to at least 0.25 M, where calculations indicate ion pairing to be significant, there is no notable effect of ionic strength on the spectrum. The experimental spectra are fitted to aid comparison with computations. Classical molecular dynamic simulations are carried out on dilute aqueous Cl - , Na + , and NaCl, producing radial distribution functions in reasonable agreement with experiment and, for NaCl, clear evidence of ion pairing. Clusters are extracted from the simulations for quantum mechanical excited state calculations. Accurate ab initio coupled-cluster benchmark calculations on a small number of representative clusters are carried out and used to identify and validate an efficient protocol based on time-dependent density functional theory. The latter is used to carry out quantum mechanical calculations on thousands of clusters. The resulting computed spectrum is in excellent agreement with experiment for the peak position, with little influence from ion pairing, but is in qualitative disagreement on the width, being only about half as wide. It is concluded that simulation by classical molecular dynamics fails to provide an adequate variety of structures to explain the experimental CTTS spectrum of aqueous Cl - .

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

The three-point energy correlator in the coplanar limit

Energy correlators are a type of observables that measure how energy is distributed across multiple detectors as a function of the angles between pairs of detectors. In this paper, we study the three-point energy correlator (EEEC) at lepton colliders in the three-particle near-to-plane (coplanar) limit. The leading-power contribution in this limit is governed by the three-jet (trijet) configuration. We introduce a new approach by projecting the EEEC onto the volume of the parallelepiped formed by the unit vectors aligned with three detected final-state particles. Analogous to the back-to-back limit of the two-point energy correlator probing the dijet configuration, the small-volume limit of the EEEC probes the trijet configuration. We derive a transverse momentum dependent (TMD) based factorization theorem that captures the soft and collinear logarithms in the coplanar limit, which enables us to achieve the next-to-next-to-next-to-leading logarithm (N 3 LL) resummation. To our knowledge, this is the first N 3 LL result for a trijet event shape. Additionally, we demonstrate that a similar factorization theorem can be applied to the fully differential EEEC in the three-particle coplanar limit, which provides a clean environment for studying different coplanar trijet shapes.

Effective Field Theories of QCD↗
DOE OSTI2024

Machine Learning Potentials with the Iterative Boltzmann Inversion: Training to Experiment

Methodologies for training machine learning potentials (MLPs) with quantum-mechanical simulation data have recently seen tremendous progress. Experimental data have a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. Here, we investigate a training procedure based on iterative Boltzmann inversion that produces a pair potential correction to an existing MLP using equilibrium radial distribution function data. By applying these corrections to an MLP for pure aluminum based on density functional theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require autodifferentiating through a molecular dynamics solver and does not make assumptions about the MLP architecture. Our results suggest a practical framework for incorporating experimental data into machine learning models to improve the accuracy of molecular dynamics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pursuit of paired dijet resonances in the Run 2 dataset with ATLAS

New particles with large masses that decay into hadronically interacting particles are predicted by many models of physics beyond the Standard Model. A search for a massive resonance that decays into pairs of dijet resonances is performed using 140 fb –1 of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded by the ATLAS detector during Run 2 of the Large Hadron Collider. Resonances are searched for in the invariant mass of the tetrajet system, and in the average invariant mass of the pair of dijet systems. A data-driven background estimate is obtained by fitting the tetrajet and dijet invariant mass distributions with a four-parameter dijet function and a search for local excesses from resonant production of dijet pairs is performed. No significant excess of events beyond the Standard Model expectation is observed, and upper limits are set on the production cross sections of new physics scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Beam Commissioning Result of Polarized Target at SpinQuest

SpinQuest is a fixed-target experiment at Fermilab to measure the Drell-Yan process using transversely polarized NH$_3$ and ND$_3$ targets and unpolarized 120-GeV proton beam. In the Drell-Yan process, a quark in one scattering hadron and an anti-quark in the other hadron annihilate into a virtual photon and then decay into a muon (lepton) pair. The angular distribution of final-state muon pairs with respect to the target polarization is sensitive to the Sivers function of light antiquarks in the nucleon, which is one of the eight Transverse Momentum Dependent (TMD) parton distribution functions. The intensity of the proton beam is as high as $10^{12}$ protons/second, in order to accumulate the required statistics of Drell-Yan events. The polarized target system at SpinQuest has been carefully designed to accommodate such a high beam intensity. It is equipped with an evaporation refrigerator with a cooling power of around 5 W at 1 K. The target temperature is maintained at 1 K even under heat load caused by the proton beam and the DNP microwave. The target system and the spectrometer were commissioned with the proton beam from May through July 2024. The experimental setup of SpinQuest and the performance of the target system during the beam commissioning is reported.

Nakano, Kenichi [Virginia U.]↗
NASA NTRS1976

A mechanism for the production of pulsar radio radiation

Recent theoretical studies have suggested that avalanches of electron-positron pairs are produced by spark discharges near the magnetic poles of pulsars. This pair production mechanism may result in a double-peaked particle distribution function. We investigate collective effects caused by a high-energy electron-positron beam traversing a low-energy electron-positron plasma for an astrophysically interesting range of physical parameters. The principal result of this paper is to identify transverse unstable waves whose phase propagation may occur at any arbitrary angle (other than zero) with respect to the magnetic field but whose group velocity is along the magnetic field and approximately the speed of light. It is suggested that the observed radio pulses from pulsars may be caused by the focusing of these transverse waves along the open field lines extending from a pulsar's polar cap to the light cylinder.

Hardee, P. E.↗
NASA NTRS1993

Thermophysical properties of simple liquid metals: A brief review of theory

In this paper, we review the current theory of the thermophysical properties of simple liquid metals. The emphasis is on thermodynamic properties, but we also briefly discuss the nonequilibrium properties of liquid metals. We begin by defining a 'simple liquid metal' as one in which the valence electrons interact only weakly with the ionic cores, so that the interaction can be treated by perturbation theory. We then write down the equilibrium Hamiltonian of a liquid metal as a sum of five terms: the bare ion-ion interaction, the electron-electron interaction, the bare electron-ion interaction, and the kinetic energies of electrons and ions. Since the electron-ion interaction can be treated by perturbation, the electronic part contributes in two ways to the Helmholtz free energy: it gives a density-dependent term which is independent of the arrangement of ions, and it acts to screen the ion-ion interaction, giving rise to effective ion-ion pair potentials which are density-dependent, in general. After sketching the form of a typical pair potential, we briefly enumerate some methods for calculating the ionic distribution function and hence the Helmholtz free energy of the liquid: monte Carlo simulations, molecular dynamics simulations, and thermodynamic perturbation theory. The final result is a general expression for the Helmholtz free energy of the liquid metal. It can be used to calculate a wide range of thermodynamic properties of simple metal liquids, which we enumerate. They include not only a range of thermodynamic coefficients of both metals and alloys, but also many aspects of the phase diagram, including freezing curves of pure elements and phase diagrams of liquid alloys (including liquidus and solidus curves). We briefly mention some key discoveries resulting from previous applications of this method, and point out that the same methods work for other materials not normally considered to be liquid metals (such as colloidal suspensions, in which the suspended microspheres behave like ions screened by the salt solution in which they are suspended). We conclude with a brief discussion of some non-equilibrium (i.e., transport) properties which can be treated by an extension of these methods. These include electrical resistivity, thermal conductivity, viscosity, atomic self-diffusion coefficients, concentration diffusion coefficients in alloys, surface tension and thermal emissivity. Finally, we briefly mention two methods by which the theory might be extended to non-simple liquid metals: these are empirical techniques (i.e., empirical two- and three-body potentials), and numerical many-body approaches. Both may be potentially applicable to extremely complex systems, such as nonstoichiometric liquid semiconductor alloys.

Stroud, David↗
DOE OSTI2024

Measurements of soil protist richness and community composition are influenced by primer pair, annealing temperature, and bioinformatics choices

ABSTRACT Protists are a diverse and understudied group of microbial eukaryotic organisms especially in terrestrial environments. Advances in molecular methods are increasing our understanding of the distribution and functions of these creatures; however, there is a vast array of choices researchers make including barcoding genes, primer pairs, PCR settings, and bioinformatic options that can impact the outcome of protist community surveys. Here, we tested four commonly used primer pairs targeting the V4 and V9 regions of the 18S rRNA gene using different PCR annealing temperatures and processed the sequences with different bioinformatic parameters in 10 diverse soils to evaluate how primer pair, amplification parameters, and bioinformatic choices influence the composition and richness of protist and non-protist taxa using Illumina sequencing. Our results showed that annealing temperature influenced sequencing depth and protist taxon richness for most primer pairs, and that merging forward and reverse sequencing reads for the V4 primer pairs dramatically reduced the number of sequences and taxon richness of protists. The data sets of primers that targeted the same 18S rRNA gene region (e.g., V4 or V9) had similar protist community compositions; however, data sets from primers targeting the V4 18S rRNA gene region detected a greater number of protist taxa compared to those prepared with primers targeting the V9 18S rRNA region. There was limited overlap of protist taxa between data sets targeting the two different gene regions (80/549 taxa). Together, we show that laboratory and bioinformatic choices can substantially affect the results and conclusions about protist diversity and community composition using metabarcoding. IMPORTANCE Ecosystem functioning is driven by the activity and interactions of the microbial community, in both aquatic and terrestrial environments. Protists are a group of highly diverse, mostly unicellular microbes whose identity and roles in terrestrial ecosystem ecology have been largely ignored until recently. This study highlights the importance of choices researchers make, such as primer pair, on the results and conclusions about protist diversity and community composition in soils. In order to better understand the roles protist taxa play in terrestrial ecosystems, biases in methodological and analytical choices should be understood and acknowledged.

Biotechnology & Applied Microbiology↗
NASA NTRS1990

H + O3 Fourier-transform infrared emission and laser absorption studies of OH(X2Pi) radical - An experimental dipole moment function and state-to-state Einstein A coefficients

FTIR emission/absorption spectroscopy is used to measure the relative intensities of 88 pairs of rovibrational transitions of OH(X2Pi) distributed over 16 vibrational bands. The experimental technique used to obtain the Einstein A ratios is discussed. The dipole moment function which follows from the intensity ratios along with Einstein A coefficients calculated from mu(r) is presented.

Nelson, David D., Jr.↗
DOE OSTI2024

Elucidating the Impact of Cis – Trans Organic Structure Directing Agent Isomer Ratios on the Aluminum Distribution Within SSZ-39

Despite their widespread use, the mechanisms governing the synthesis of zeolite catalysts are still poorly understood. A notable example of this problem is the uncertainty surrounding the influence of synthesis conditions on the placement of Al atoms in the zeolite framework which determines the active sites available for catalytic species. In this work, the role of the cis to trans isomer ratio of the OSDA N,N-dimethyl-3-5-dimethylpiperidinium on the energetics of 26 distinct Al pair distributions in SSZ-39 is examined both in the presence and absence of Na using density functional theory calculations. The initial orientation of the OSDA was found to have a significant impact on the final energies present, necessitating the screening of a large number of initial orientations with force field calculations and single point DFT calculations. Ground state energies were found to vary significantly with the ratio of cis to trans OSDAs with a Boltzmann distribution revealing the most likely Al pair distributions shift from sharing the same 8 membered rings to sharing the same double six membered rings to having no shared subunits as one increases the amount of cis OSDA present within the framework. The presence of Na was found to favor Al pair distributions where both Als occupied the same 6-membered ring. When an implicit solvent model was used to evaluate ground state energies the ideal Na sites shifted from 6-membered rings to empty SSZ-39 cages while OSDA positions and orientations remained largely the same. To provide insight on how kinetic factors may influence Al distributions, formation energies we calculated for connected double six membered rings. Further, these formation energies revealed a preference for Al pairs to occupy the same 4-membered ring which indicates kinetic and thermodynamic control may lead to different Al distributions in SSZ-39.

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

Dynamic-local-field approximation for the quantum solids

A local-molecular-field description for the ground-state properties of the quantum solids is presented. The dynamical behavior of atoms contributing to the local field, which acts on an arbitrary pair of test particles, is incorporated by decoupling the pair correlations between these field atoms. The energy, pressure, compressibility, single-particle-distribution function, and the rms atomic deviations about the equilibrium lattice sites are calculated for H2, He-3, and He-4 over the volume range from 5 to 24.5 cu cm/mole. The results are in close agreement with existing Monte Carlo calculations wherever comparisons are possible. At very high pressure, the results agree with simplified descriptions which depend on negligible overlap of the system wave function between neighboring lattice sites.

Etters, R. D.↗
DOE OSTI2026

Thermophysical Properties of Liquid Tritium: A Path Integral Monte Carlo Study

Here, we present worm-algorithm, path integral Monte Carlo simulations of bulk liquid tritium. The simulations are benchmarked against empirically known thermophysical properties of liquid deuterium and liquid tritium. Results for the pair correlation function, chemical potential, isothermal compressibility, isochoric heat capacity, and single-particle momentum distributions are reported. Given the benchmark comparisons, our predictions of liquid tritium properties are expected to be accurate to within a few percent. Our simulations unambiguously demonstrate the significance of nuclear quantum effects to the properties of liquid tritium. In particular, under saturated vapor pressure, the average molecular kinetic energy of the liquid is found to be more than 60% higher than the value expected from the classical equipartition theorem.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Efficient Active Learning for Gaussian Process Classification by Error Reduction

Active learning sequentially selects the best instance for labeling by optimizing an acquisition function to enhance data/label efficiency. The selection can be either from a discrete instance set (pool-based scenario) or a continuous instance space (query synthesis scenario). In this work, we study both active learning scenarios for Gaussian Process Classification (GPC). The existing active learning strategies that maximize the Estimated Error Reduction (EER) aim at reducing the classification error after training with the new acquired instance in a onestep-look-ahead manner. The computation of EER-based acquisition functions is typically prohibitive as it requires retraining the GPC with every new query. Moreover, as the EER is not smooth, it can not be combined with gradient-based optimization techniques to efficiently explore the continuous instance space for query synthesis. To overcome these critical limitations, we develop computationally efficient algorithms for EER-based active learning with GPC. Further, we derive the joint predictive distribution of label pairs as a one-dimensional integral, as a result of which the computation of the acquisition function avoids retraining the GPC for each query, remarkably reducing the computational overhead. We also derive the gradient chain rule to efficiently calculate the gradient of the acquisition function, which leads to the first query synthesis active learning algorithm implementing EER-based strategies. Our experiments clearly demonstrate the computational efficiency of the proposed algorithms. We also benchmark our algorithms on both synthetic and real-world datasets, which show superior performance in terms of sampling efficiency compared to the existing state-of-the-art algorithms.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1957

Theoretical Calculation of the Power Spectra of the Rolling and Yawing Moments on a Wing in Random Turbulence

The correlation functions and power spectra of the rolling and yawing moments on an airplane wing due to the three components of continuous random turbulence are calculated. The rolling moments to the longitudinal (horizontal) and normal (vertical) components depend on the spanwise distributions of instantaneous gust intensity, which are taken into account by using the inherent properties of symmetry of isotropic turbulence. The results consist of expressions for correlation functions or spectra of the rolling moment in terms of the point correlation functions of the two components of turbulence. Specific numerical calculations are made for a pair of correlation functions given by simple analytic expressions which fit available experimental data quite well. Calculations are made for four lift distributions. Comparison is made with the results of previous analyses which assumed random turbulence along the flight path and linear variations of gust velocity across the span.

Eggleston, John M↗
NASA NTRS1992

Bivariate extreme value distributions

In certain engineering applications, such as those occurring in the analyses of ascent structural loads for the Space Transportation System (STS), some of the load variables have a lower bound of zero. Thus, the need for practical models of bivariate extreme value probability distribution functions with lower limits was identified. We discuss the Gumbel models and present practical forms of bivariate extreme probability distributions of Weibull and Frechet types with two parameters. Bivariate extreme value probability distribution functions can be expressed in terms of the marginal extremel distributions and a 'dependence' function subject to certain analytical conditions. Properties of such bivariate extreme distributions, sums and differences of paired extremals, as well as the corresponding forms of conditional distributions, are discussed. Practical estimation techniques are also given.

Elshamy, M.↗
DOE OSTI2020

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

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

Permissiveness and competition within and between Neurospora crassa syncytia

A multinucleate syncytium is a common growth form in filamentous fungi. Comprehensive functions of the syncytial state remain unknown, but it likely allows for a wide range of adaptations to enable filamentous fungi to coordinate growth, reproduction, responses to the environment, and to distribute nuclear and cytoplasmic elements across a colony. Indeed, the underlying mechanistic details of how syncytia regulate cellular and molecular processes spatiotemporally across a colony are largely unexplored. Here, we implemented a strategy to analyze the relative fitness of different nuclear populations in syncytia of Neurospora crassa , including nuclei with loss-of-function mutations in essential genes, based on production of multinucleate asexual spores using flow cytometry of pairings between strains with differentially fluorescently tagged nuclear histones. The distribution of homokaryotic and heterokaryotic asexual spores in pairings was assessed between different auxotrophic and morphological mutants, as well as with strains that were defective in somatic cell fusion or were heterokaryon incompatible. Mutant nuclei were compartmentalized into both homokaryotic and heterokaryotic asexual spores, a type of bet hedging for maintenance and evolution of mutational events, despite disadvantages to the syncytium. However, in pairings between strains that were blocked in somatic cell fusion or were heterokaryon incompatible, we observed a “winner-takes-all” phenotype, where asexual spores originating from paired strains were predominantly one genotype. These data indicate that syncytial fungal cells are permissive and tolerate a wide array of nuclear functionality, but that cells/colonies that are unable to cooperate via syncytia formation actively compete for resources.

59 BASIC BIOLOGICAL SCIENCES↗