Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Generalized 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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree–Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Furthermore, our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.

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

Anomalous structural recovery in the near glass transition range in a polymer glass: Data revisited in light of temperature variability in vacuum oven‐based experiments*

Abstract There is considerable interest in data reported by Cangialosi et al. [Cangialosi, D., Boucher, V. M., Alegría, A., & Colmenero, J. (2013). Physical Review Letters , 111 (9), 095701] that purportedly showed unusual irregular responses in the isothermal structural recovery behavior of glasses because the observation challenges the general view of the dynamics being smooth functions of time as observed in extensive dilatometric experiments by Kovacs [Kovacs, A. J. (1964). In Fortschritte der Hochpolymeren‐Forschung , 3 , 394–507. Springer, Berlin, Heidelberg.] in the 1960s. The reported data show an apparent two mechanism structural (enthalpy) recovery at aging temperatures ranging from 6 to 17 K below which is also controversial as it could not be reproduced in work from Koh and Simon [Koh, Y. P., & Simon, S. L. (2013). Macromolecules , 46(14), 5815–5821.] where a second plateau observed in the enthalpy loss curve was not reproduced in nearly 1‐year‐long experiments at 15 K below . Therefore, it is important to determine what might be possible reasons for the apparent non‐smooth relaxations. Here, we examine the possibility that the poor temperature control of a typical vacuum oven could explain the anomalous results. We used the Tool–Narayanaswamy–Moynihan (TNM) model of structural recovery to calculate the effect of typical vacuum oven temperature variation on the structural recovery of polystyrene and find that we can reproduce the reported experimental results. The issue of temperature control using a vacuum oven is discussed, and data from thermocouple measurements of temperature at various locations inside a vacuum oven are shown.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Generalized Ellipsometry Measurements of Crystalline Thin Film and Bulk Tin Oxide

Here, several bulk and thin-film crystals of SnO2 are grown and examined using generalized ellipsometry techniques. The bulk samples are grown using the chemical vapor transport technique and thin films of SnO2 are grown using the pulsed laser deposition technique. The bulk samples are examined using the two-modulator generalized ellipsometry microscope (2-MGEM) at normal incidence and the spectroscopic two-modulator generalized ellipsometer (2-MGE). The spectroscopic optical functions of tin oxide are then obtained using the 2-MGE from 1.46 to 6.2 eV. The material is highly birefringent, and the ordinary bandgap is less than the extraordinary band edge. 2-MGE measurements are also made on thin-film samples of crystalline tin oxide grown on sapphire and rutile, showing no cross polarization. Because of the complicated morphology of the tin oxide films grown on sapphire, the ellipsometry data are simulated using the Tauc–Lorentz model. Films grown on rutile had the optic axis perpendicular to the sample surface, but the film is strained, resulting in a more complicated ellipsometric spectrum. These films are modeled using the air/surface roughness/tin oxide/interface/rutile model, where the roughness and interface are modeled using the incomplete Beta function.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Solar Energy Research Institute for India and the United States (SERIIUS): Lessons and Results from a Binational Consortium

This book describes the development, functioning, and results of a successful binational program to promote significant scientific advances in Earth-abundant photovoltaics (PV) and concentrated solar power (CSP), advanced process/manufacturing technologies, multiscale modeling and reliability testing, and analysis of integrated solar energy systems. SERIIUS is a consortium between India and the United States dedicated to developing new solar technologies and assessing their potential impact in the two countries. The consortium consists of nearly 50 institutions including academia, national laboratories, and industry, with the goal of developing significant new technologies in all areas of solar deployment. In addition, the program focused on workforce development through graduate students, post-doctoral students, and an international exchange program. Particular emphasis was placed on the following efforts: creating disruptive technologies in PV and CSP through high-impact fundamental and applied research and development (R&D); identifying and quantifying the critical technical, economic, and policy issues for solar energy development and deployment in India; overcoming barriers to technology transfer by teaming research institutions and industry in an effective project structure; building a new platform for binational collaboration using a formalized R&D project structure, along with effective management, coordination, and decision processes; creating a sustainable network and workforce development program from which to build large collaborations and fostering a collaborative culture and outreach programs. This includes using existing and new methodologies for collaboration based on advanced electronic and web-based communication to facilitate functional international teams. The book summarizes the general lessons learned from these experiences.

14 SOLAR ENERGY↗
DOE OSTI2023

Holographic open quantum systems: toy models and analytic properties of thermal correlators

We present a unified picture of open quantum systems, the theory of a system probing a noisy thermal environment, distilling lessons learnt from previous holographic analyses. Our treatment is applicable both when the system is coupled to short-lived (Markovian), and long-lived (non-Markovian) environmental degrees of freedom. The thermal environment is modeled using an asymptotically AdS black hole, and the systems of interest are simple probe field theories. The effective stochastic dynamics of the system is governed by real-time thermal correlators, which we compute using the gravitational Schwinger-Keldysh (grSK) geometry. We describe the structure of arbitrary tree-level contact and exchange Witten diagrams in the grSK geometry. In particular, we argue, that all such diagrams reduce to integrals supported on a single copy of the exterior of the black hole. The integrand is obtained as a multiple discontinuity of a function comprising ingoing boundary-bulk propagators, monodromy functions which appear as radial Boltzmann weights, and vertex factors. These results allow us to deduce the analytic structure of real-time thermal n-point functions in holographic CFTs. We illustrate the general statements by a two-dimensional toy model, dual to fields in the BTZ background, which we argue captures many of the essential features of generic open holographic QFTs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Electrochemical Reduction Selectivity of Crotonaldehyde on Copper

Lignocellulosic waste is a potential feedstock for the generation of fuels and commodity chemicals, but existing conversion methods are too cost-intensive to be viable long-term solutions. Electrochemical reductions are promising for decentralized biomass valorization due to their modular scaling and capacity to run intermittently and without high temperatures or pressures. Using crotonaldehyde as a multi-functional model compound for the many partially unsaturated oxygenates found in processed biomass, we here demonstrate the production of butanal, butanol, butene, and butane (variously useful as commodity chemicals and major components of liquified petroleum gas) under ambient conditions by reductive bulk electrolysis with a copper mesh working electrode. We identify an optimum potential for reduced organic production under the reaction conditions and compare product distributions from reductions of intermediate species to further propose branching reaction pathways. Though butanal is typically the most abundant product from crotonaldehyde reduction, most of the butene and butane appear to result from a pathway involving initial reduction of the aldehyde group. We discuss evidence that selectivity is driven by interplay between crotonaldehyde reduction, local pH shifts due to the hydrogen evolution reaction, and changes in site reactivity and availability due to electrode fouling. This demonstration of model electrochemical biomass valorization also serves to inform further exploration into reduction of multi-functional molecules and electrochemical biomass processing in general.

biomass↗
DOE OSTI2022

Entropy–Preserving and Entropy–Stable Relaxation IMEX and Multirate Time–Stepping Methods

In this work, we propose entropy-preserving and entropy-stable partitioned Runge–Kutta (RK) methods. In particular, we extend the explicit relaxation Runge–Kutta methods to IMEX–RK methods and a class of explicit second-order multirate methods for stiff problems arising from scale-separable or grid-induced stiffness in a system. The proposed approaches not only mitigate system stiffness but also fully support entropy-preserving and entropy-stability properties at a discrete level. The key idea of the relaxation approach is to adjust the step completion with a relaxation parameter so that the time-adjusted solution satisfies the entropy condition at a discrete level. The relaxation parameter is computed by solving a scalar nonlinear equation at each timestep in general; however, as for a quadratic entropy function, we theoretically derive the explicit form of the relaxation parameter and numerically confirm that the relaxation parameter works the Burgers equation. Several numerical results for ordinary differential equations and the Burgers equation are presented to demonstrate the entropy-conserving/stable behavior of these methods. We also compare the relaxation approach and the incremental direction technique for the Burgers equation with and without a limiter in the presence of shocks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Angular momentum eigenstates of the isotropic 3-D harmonic oscillator: Phase-space distributions and coalescence probabilities

The isotropic 3-dimensional harmonic oscillator potential can serve as an approximate description of many systems in atomic, solid state, nuclear, and particle physics. In particular, the question of 2 particles binding (or coalescing) into angular momentum eigenstates in such a potential has interesting applications. Here we compute the probabilities for coalescence of two distinguishable, non-relativistic particles into such a bound state, where the initial particles are represented by generic wave packets of given average positions and momenta. We use a phase-space formulation and hence need the Wigner distribution functions of angular momentum eigenstates in isotropic 3-dimensional harmonic oscillators. These distribution functions have been discussed in the literature before but we utilize an alternative approach to obtain these functions. Along the way, we derive a general formula that expands angular momentum eigenstates in terms of products of 1-dimensional harmonic oscillator eigenstates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modified Energy Span Analysis Reveals Heterogeneous Catalytic Kinetics

Mechanistic modeling is a cornerstone of catalyst development generally conducted with microkinetic models or density functional theory-based energy profiles. We extend the energy span model of homogeneous catalysis to heterogeneous systems by introducing the modified energy span analysis (MESA) model by implementing collision theory and gas-phase concentration effects. We determine analytically turnover frequencies, coverages, rate-determining steps, apparent activation energies, and reaction orders in agreement with microkinetic and kinetic Monte Carlo simulations. The model applies to discrete catalyst state systems, including single-atom catalysts, homogeneous systems, and microporous materials. A generalizable rate expression is derived for reactor modeling or mechanistic insights from solely experimentally measured reaction orders. Here, we illustrate MESA on three published mechanisms and reveal unexpected phenomena absent in traditional energy span profiles.

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

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗
DOE OSTI2022

Frustrated Coulombic and Cation Size Effects on Nanoscale Boehmite Aggregation: A Tumbler Small- and Ultra-Small-Angle Neutron Scattering Study

To better understand the effects of solution chemistry on particle aggregation in the complex legacy tank wastes at the Hanford (WA) and Savannah River (SC) sites, we have performed a series of tumbler small- and ultra-small-angle neutron scattering experiments on 20 wt % solid slurries of nanoparticulate aluminum oxyhydroxide (boehmite) with M 1+ nitrates of various concentrations and radii. In this study, the solutes consisted of H, Li, Na, K, and Rb nitrates at 10 –5 , 10 –3 , 10 –1 , 2, and 4 molal (m) concentrations, as well as in pure H 2 O. Synthetic boehmite nanoparticles were used with a size range from ~20 to 30 nm. Tumbler cells were used to keep the solids from settling. Although particles initially form individual rhombohedral platelets, once placed in solution, they quickly form well-bonded stacks, primary aggregates, up to ~1500 Å long, and a second level of aggregates whose concentration and structure vary as a function of cation type and concentration. Aggregation generally increases with increased solute concentration and with cation radius up to a concentration somewhat above 10 –1 m, at which point the trend reverses. Primary aggregates become more rodlike and larger. The Kirkwood-like reversal probably reflects a change from Derjaguin–Landau–Verwey–Overbeek (DLVO)/Debye behavior controlled by surface chemistry to a frustrated Coulombic system controlled by the solution structure. These data suggest that an understanding of the effects of salt concentration and chemistry on nanoparticle aggregate structures provides useful physical insights into the microscopic origin of slurry rheology in the Hanford and Savannah River legacy wastes.

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

Numerical study of $\delta$-function current sheets arising from resonant magnetic perturbations

We report general three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially regions of stochastic field lines when islands overlap. Numerically resolving singular current sheets in the ideal magnetohydrodynamics (MHD) limit has been a significant challenge. This work presents numerical solutions of the Hahm–Kulsrud–Taylor (HKT) problem, which is a prototype for resonant singular current sheet formation. The HKT problem is solved by two codes: a Grad–Shafranov (GS) solver and the Stepped Pressure Equilibrium Code (SPEC) code. The GS solver has built-in nested flux surfaces with prescribed magnetic fluxes. The SPEC code implements multi-region relaxed magnetohydrodynamics (MRxMHD), whereby the solution relaxes to a Taylor state in each region while maintaining force balance across the interfaces between regions. As the number of regions increases, the MRxMHD solution appears to approach the ideal MHD solution assuming a continuum of nested flux surfaces. We demonstrate agreement between the numerical solutions obtained from the two codes through a convergence study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Cosmological angular trispectra and non-Gaussian covariance

Angular cosmological correlators are infamously difficult to compute due to the highly oscillatory nature of the projection integrals. Motivated by recent development on analytic approaches to cosmological perturbation theory, in this paper we present an efficient method for computing cosmological four-point correlations in angular space, generalizing previous works on lower-point functions. This builds on the FFTLog algorithm that approximates the matter power spectrum as a sum over power-law functions, which makes certain momentum integrals analytically solvable. The computational complexity is drastically reduced for correlators in a "separable" form—we define a suitable notion of separability for cosmological trispectra, and derive formulas for angular correlators of different separability classes. As an application of our formalism, we compute the angular galaxy trispectrum at tree level, with and without primordial non-Gaussianity. This includes effects of redshift space distortion and bias parameters up to cubic order. We also compute the non-Gaussian covariance of the angular matter power spectrum due to the connected four-point function, beyond the Limber approximation. Finally, we demonstrate that, in contrast to the standard lore, the Limber approximation can fail for the non-Gaussian covariance computation even for large multipoles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Meta-GGA performance in solids at almost GGA cost

A recent modification, r 2 SCAN, of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (GGA) exchange-correlation functional mostly eliminates numerical instabilities and attendant integration grid sensitivities exhibited by SCAN. Here, we show that the successful deorbitalization of SCAN to SCAN-L (SCAN with density Laplacian dependence) carries over directly to yield r 2 SCAN-L. A major benefit is that the high iteration counts that hindered the use of SCAN-L are eliminated in r 2 SCAN-L. It therefore is a computationally much faster meta-GGA than its orbital-dependent antecedent. In this work, validation data for molecular heats of formation, bond lengths, and vibration frequencies (G3/99X, T96-R, and T82-F test sets, respectively) and on lattice constants, and cohesive energies (for 55 solids) and bulk moduli (for 40 solids) are provided. In addition, we show that the overmagnetization of bcc Fe, hcp Co, and fcc Ni persists in r 2 SCAN but does not appear in r 2 SCAN-L. Distinct from SCAN, both r 2 SCAN and r 2 SCAN-L give the correct nonmagnetic ground state for bcc V.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Evolution of magnetic field induced ordering in the layered quantum Heisenberg triangular-lattice antiferromagnet Ba 3 CoSb 2 O 9

Quantum fluctuations in the effective spin-12 layered triangular-lattice quantum Heisenberg antiferromagnet Ba 3 CoSb 2 O 9 lift the classical degeneracy of the antiferromagnetic ground state in magnetic field, producing a series of novel spin structures for magnetic fields applied within the crystallographic ab plane, including a celebrated collinear “up-up-down” spin ordering with magnetization equal to 1/3 of the saturation magnetization over an extended field range. Theoretically unresolved, however, are the effects of interlayer antiferromagnetic coupling and transverse magnetic fields on the ground states of this system. Additional magnetic field induced phase transitions are theoretically expected and in some cases have been experimentally observed, but details regarding their number, location, and physical character appear inconsistent with the predictions of existing models. Conversely, an absence of experimental measurements as a function of magnetic-field orientation has left other key predictions of these models untested. To address these issues, in this study we have used specific heat, neutron diffraction, thermal conductivity, and magnetic torque measurements to map out the phase diagram as a function of magnetic field intensity and orientation relative to the crystallographic ab plane. For H||ab, we have discovered an additional magnetic field induced phase transition at low temperature and an unexpected tetracritical point in the high-field phase diagram, which coupled with the apparent second-order nature of the phase transitions eliminates several theoretically proposed spin structures for the high-field phases. Our calorimetric measurements as a function of magnetic field orientation are in general agreement with theory for field-orientation angles close to plane parallel (H||a) but diverge at angles near plane perpendicular; a predicted convergence of two phase boundaries at finite angle and a corresponding change in the order of the field induced phase transition are not observed experimentally. Our results emphasize the role of interlayer coupling in selecting and stabilizing field induced phases, provide guidance on the nature of the magnetic order in each phase, and reveal the need for new physics to account for the nature of magnetic ordering in this archetypal two-dimensional spin-12 triangular-lattice quantum Heisenberg antiferromagnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Pseudogauge dependence of the spin polarization and of the axial vortical effect

The mean spin polarization vector of spin 1/2 particles in a relativistic fluid at local thermal equilibrium (LTE) in different pseudogauges (PGs), i.e., with different choices for the decomposition of orbital and spin angular momentum, is obtained. The spin polarization obtained in the canonical PG differs from the one obtained in the Belinfante PG. It is found that this difference can not be written by replacing the thermal vorticity with the spin potential in the usual polarization formula. Other PG choices affect the contribution of thermal shear to the spin polarization. So, the choice of a PG affects the predictions for the polarization measured in heavy-ion collisions. In general, it is shown that the Wigner function of a noninteracting Dirac field at LTE is affected by the PG transformations. Explicit expressions of the mean axial current are also calculated and it is found that even the axial vortical effect conductivity depends on the PG.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Collins-Soper kernel in the QCD instanton vacuum

We outline a general framework for evaluating the nonperturbative soft functions in the quantum chromodynamics (QCD) instanton vacuum. In particular, from the soft function we derive the Collins-Soper (CS) kernel, which drives the rapidity evolution of the transverse-momentum-dependent parton distributions. The resulting CS kernel, when supplemented with the perturbative contribution, agrees well with recent lattice results and some phenomenological parametrizations. Moreover, our CS kernel depends logarithmically on the large quark transverse separation, providing a key constraint on its phenomenological parametrization. Finally, a lattice calculation can be directly compared to our generic results in Euclidean signature, thus providing a new approach for evalulating the soft function and extracting the CS kernel by analytical continuation.

QCD phenomenology↗
DOE OSTI2020

Giant piezoelectricity in oxide thin films with nanopillar structure

High-performance piezoelectric materials are critical components for electromechanical sensors and actuators. For more than 60 years, the main strategy for obtaining large piezoelectric response has been to construct multiphase boundaries, where nanoscale domains with local structural and polar heterogeneity are formed, by tuning complex chemical compositions. We used a different strategy to emulate such local heterogeneity by forming nanopillar regions in perovskite oxide thin films. We obtained a giant effective piezoelectric coefficient d 33 , f * of ~1098 picometers per volt with a high Curie temperature of ~450°C. Our lead-free composition of sodium-deficient sodium niobate contains only three elements (Na, Nb, and O). The formation of local heterogeneity with nanopillars in the perovskite structure could be the basis for a general approach to designing and optimizing various functional materials.

Science & Technology - Other Topics↗