Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative 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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Perturbative QCD contribution to transverse single spin asymmetries in the Drell-Yan process and SIDIS

In a previous publication [Benić , ], we have computed the perturbative QCD contribution to transverse single spin asymmetries (SSAs) in semi-inclusive deep inelastic scattering (SIDIS) involving the g T ( x ) distribution. In this paper, we first present a more efficient derivation of the asymmetries which is applicable to both transverse and longitudinal SSAs and correct some inconsistencies in our previous calculation. We then adapt the method to compute transverse SSAs in the Drell-Yan process proportional to g T ( x ) and its gluonic counterpart and discuss the crossing symmetry between the results for SIDIS and the Drell-Yan process. Finally, we present numerical results for various asymmetries measurable at the EIC, RHIC, COMPASS, and Fermilab (SpinQuest), including also part of the genuine twist-three corrections to g T ( x ) from a recent global analysis. We find that the asymmetries can reach percent-level magnitude, if the kinematics predominantly probes the large- x (valence) region of the polarized proton, but remain at subpercent levels otherwise. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Roadmap on methods and software for electronic structure based simulations in chemistry and materials

This Roadmap article provides a succinct, comprehensive overview of the state of electronic structure methods and software for molecular and materials simulations. Seventeen distinct sections collect insights by 51 leading scientists in the field. Each contribution addresses the status of a particular area, as well as current challenges and anticipated future advances, with a particular eye towards software related aspects and providing key references for further reading. Foundational sections cover density functional theory and its implementation in real-world simulation frameworks, Green's function based many-body perturbation theory, wave-function based and stochastic electronic structure approaches, relativistic effects and semiempirical electronic structure theory approaches. Subsequent sections cover nuclear quantum effects, real-time propagation of the electronic structure, challenges for computational spectroscopy simulations, and exploration of complex potential energy surfaces. The final sections summarize practical aspects, including computational workflows for complex simulation tasks, the impact of current and future high-performance computing architectures, software engineering practices, education and training to maintain and broaden the community, as well as the status of and needs for electronic structure based modeling from the vantage point of industry environments. Overall, the field of electronic structure software and method development continues to unlock immense opportunities for future scientific discovery, based on the growing ability of computations to reveal complex phenomena, processes and properties that are determined by the make-up of matter at the atomic scale, with high precision.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Overview of the recent experimental research on the J-TEXT tokamak

The J-TEXT capability is enhanced compared to two years ago with several upgrades of its diagnostics and the increase of electron cyclotron resonance heating (ECRH) power to 1 MW. With the application of electron cyclotron wave (ECW), the ECW assisted plasma startup is achieved; the tearing mode is suppressed; the toroidal injection of 300 kW ECW drives around 24 kA current; fast electrons are generated with toroidal injected ECW and the runaway current conversion efficiency increases with ECRH power. The mode coupling between 2/1 and 3/1 modes are extensively studied. The coupled 2/1 and 3/1 modes usually lead to major disruption. Their coupling can be either suppressed or avoided by external resonant magnetic perturbation fields and hence avoids the major disruption. It is also found that the 2/1 threshold of external field is significantly reduced by a pre-excited 3/1 mode, which can be either a locked island or an external kink mode. The disruption control is studied by developing prediction methods capable of cross tokamak application and by new mitigation methods, such as the biased electrode or electromagnetic pellet injector. The high-density operation and related disruptions are studied from various aspects. Approaching the density limit, the collapse of the edge shear layer is observed and such collapse can be prevented by applying edge biasing, leading to an increased density limit. The density limit is also observed to increase, if the plasma is operated in the poloidal divertor configuration or the plasma purity is increased by increasing the pre-filled gas pressure or ECRH power during the start-up phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method

Recent advances in selected configuration interaction methods have made them competitive with the most accurate techniques available and, hence, creating an increasingly powerful tool for solving quantum Hamiltonians. In this work, we build on recent advances from the adaptive sampling configuration interaction (ASCI) algorithm. We show that a useful paradigm for generating efficient selected CI/exact diagonalization algorithms is driven by fast sorting algorithms, much in the same way iterative diagonalization is based on the paradigm of matrix vector multiplication. We present several new algorithms for all parts of performing a selected CI, which includes new ASCI search, dynamic bit masking, fast orbital rotations, fast diagonal matrix elements, and residue arrays. The ASCI search algorithm can be used in several different modes, which includes an integral driven search and a coefficient driven search. The algorithms presented here are fast and scalable, and we find that because they are built on fast sorting algorithms they are more efficient than all other approaches we considered. After introducing these techniques, we present ASCI results applied to a large range of systems and basis sets to demonstrate the types of simulations that can be practically treated at the full-CI level with modern methods and hardware, presenting double- and triple-ζ benchmark data for the G1 data set. The largest of these calculations is Si$_2$H$_6$ which is a simulation of 34 electrons in 152 orbitals. We also present some preliminary results for fast deterministic perturbation theory simulations that use hash functions to maintain high efficiency for treating large basis sets.

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

Adsorption separation of heavier isotope gases in subnanometer carbon pores

Isotopes of heavier gases including carbon ( 13 C/ 14 C), nitrogen ( 13 N), and oxygen ( 18 O) are highly important because they can be substituted for naturally occurring atoms without significantly perturbing the biochemical properties of the radiolabelled parent molecules. These labelled molecules are employed in clinical radiopharmaceuticals, in studies of brain disease and as imaging probes for advanced medical imaging techniques such as positron-emission tomography (PET). Established distillation-based isotope gas separation methods have a separation factor ( S ) below 1.05 and incur very high operating costs due to high energy consumption and long processing times, highlighting the need for new separation technologies. Here, we show a rapid and highly selective adsorption-based separation of 18 O 2 from 16 O 2 with S above 60 using nanoporous adsorbents operating near the boiling point of methane (112 K), which is accessible through cryogenic liquefied-natural-gas technology. A collective-nuclear-quantum effect difference between the ordered 18 O 2 and 16 O 2 molecular assemblies confined in subnanometer pores can explain the observed equilibrium separation and is applicable to other isotopic gases.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Veridical data science

Building and expanding on principles of statistics, machine learning, and scientific inquiry, we propose the predictability, computability, and stability (PCS) framework for veridical data science. Our framework, composed of both a workflow and documentation, aims to provide responsible, reliable, reproducible, and transparent results across the data science life cycle. The PCS workflow uses predictability as a reality check and considers the importance of computation in data collection/storage and algorithm design. It augments predictability and computability with an overarching stability principle. Stability expands on statistical uncertainty considerations to assess how human judgment calls impact data results through data and model/algorithm perturbations. As part of the PCS workflow, we develop PCS inference procedures, namely PCS perturbation intervals and PCS hypothesis testing, to investigate the stability of data results relative to problem formulation, data cleaning, modeling decisions, and interpretations. We illustrate PCS inference through neuroscience and genomics projects of our own and others. Moreover, we demonstrate its favorable performance over existing methods in terms of receiver operating characteristic (ROC) curves in high-dimensional, sparse linear model simulations, including a wide range of misspecified models. Finally, we propose PCS documentation based on R Markdown or Jupyter Notebook, with publicly available, reproducible codes and narratives to back up human choices made throughout an analysis. The PCS workflow and documentation are demonstrated in a genomics case study available on Zenodo.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Observing Profiles of Derived Kinematic Field Quantities Using a Network of Profiling Sites

Abstract Observations of thermodynamic and kinematic parameters associated with derivatives of the thermodynamics and wind fields, namely, advection, vorticity, divergence, and deformation, can be obtained by applying Green’s theorem to a network of observing sites. The five nodes that comprise the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) profiling network, spaced 50–80 km apart, are used to obtain measurements of these parameters over a finite region. To demonstrate the applicability of this technique at this location, it is first applied to gridded model output from the High-Resolution Rapid Refresh (HRRR) numerical weather prediction model, using profiles from the locations of ARM network sites, so that values calculated from this method can be directly compared to finite difference calculations. Good agreement is found between both approaches as well as between the model and values calculated from the observations. Uncertainties for the observations are obtained via a Monte Carlo process in which the profiles are randomly perturbed in accordance with their known error characteristics. The existing size of the ARM network is well suited to capturing these parameters, with strong correlations to model values and smaller uncertainties than a more closely spaced network, yet it is small enough that it avoids the tendency for advection to go to zero over a large area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Thermal Effects on Far-Field Distributed Acoustic Strain-Rate Sensors

Summary Fiber-optic cables cemented outside of the casing of an unconventional well measure crosswell strain changes during fracturing of neighboring wells with low-frequency distributed acoustic sensing (LF-DAS). As a hydraulic fracture intersects an observation well instrumented with fiber-optic cables, the fracture fluid injected at ambient temperatures can cool a section of the sensing fiber. Often, LF-DAS and distributed temperature sensing (DTS) cables are run in tandem, enabling the detection of such cooling events. The increasing use of LF-DAS for characterizing unconventional hydraulic fracture completions demands an investigation of the effects of temperature on the measured strain response by LF-DAS. Researchers have demonstrated that LF-DAS can be used to extract the temporal derivative of temperature for use as a differential-temperature-gradient sensor. However, differential-temperature-gradient sensing is predicated on the ability to filter strain components out of the optical signal. In this work, beginning with an equation for optical phase shift of LF-DAS signals, a model relating strain, temperature, and optical phase shift is explicitly developed. The formula provides insights into the relative strength of strain and temperature effects on the phase shift. The uncertainty in the strain-rate measurements due to thermal effects is estimated. The relationship can also be used to quantify uncertainties in differential-temperature-gradient sensors due to strain perturbations. Additionally, a workflow is presented to simulate the LF-DAS response accounting for both strain and temperature effects. Hydraulic fracture geometries are generated with a 3D fracture simulator for a multistage unconventional completion. The fracture width distributions are imported by a displacement discontinuity method (DDM) program to compute the strain rates along an observation well. An analytic model is used to approximate the temperature in the fracture. Using the derived formulae for optical phase shift, the model outputs are then used to compute the LF-DAS response at a fiber-optic cable, enabling the generation of waterfall plots including both strain and thermal effects. The model results suggest that before, during, and immediately following a fracture intersecting a well instrumented with fiber, the strain on the fiber drives the LF-DAS signal. However, at later times, as completion fluid cools the observation well, the temperature component of the LF-DAS signal can be equal to or exceed the strain component. The modeled results are compared to a published field case in an attempt to enhance the interpretation of LF-DAS waterfall plots. Finally, we propose a sensing configuration to identify the events when “wet fractures” (fractures with fluids) intersect the observation well.

Engineering↗
DOE OSTI2022

Taylor wave solution for a general equation of state

This document describes a solution procedure for calculating the Taylor wave behind an unsupported Chapman–Jouguet (CJ) detonation in planar, cylindrical, and spherical geometries given a general equation of state. The resulting semi-analytic solution can be utilized to examine new equation of state models for detonation products and during the verification of hydrodynamic codes. The governing partial differential equations are reduced to ordinary differential equations in both characteristic and self-similar forms. The first-order systems corresponding to each geometry are amenable to solution numerically using commonly available methods. A difficulty arises at the CJ point in radial coordinates where the similarity equations become singular. Two separate strategies are proposed to integrate the first-order system. The first one uses an asymptotic approximation near the CJ point that can be used to perturb the boundary conditions. The second one applies a change of variables which removes the singularity at the expense of an additional equation to be integrated. A test problem is provided for the Davis products equation of state to illustrate the qualitative features of the Taylor wave in each geometric configuration and compared with a Lagrangian hydrodynamics research code. A Python code listing gives an implementation using the SciPy library to assists users in generating the results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Catalog of Candidate Double and Lensed Quasars from Gaia and WISE Data

Making use of strong correlations between closely separated multiple or double sources and photometric and astrometric metadata in Gaia Early Data Release 3 (EDR3), we generate a catalog of candidate double- and multiply imaged lensed quasars and active galactic nuclei (AGNs), comprising 3140 systems. It includes two partially overlapping parts: a sample of distant (redshifts mostly greater than 1) sources with perturbed data; and systems that have been resolved into separate components by Gaia at separations less than 2''. For the first part, which is roughly one-third of the published catalog, we synthesized 0.617 million redshifts using multiple machine-learning prediction and classification methods, using independent photometric and astrometric data from Gaia EDR3 and the Wide-field Infrared Survey Explorer, with accurate spectroscopic redshifts from the Sloan Digital Sky Survey (SDSS) as a training set. Using these synthetic redshifts, we estimate a 4.9% rate of interlopers with spectroscopic redshifts below 1 in this part of the catalog. Unresolved candidate double and dual AGNs and quasars are selected as sources with a marginally high BP/RP excess factor (phot_bp_rp_excess_factor), which is sensitive to source extent, limiting our search to high-redshift quasars. For the second part of the catalog, additional filters on measured parallax and near-neighbor statistics are applied to diminish the propagation of the remaining stellar contaminants. The estimated rate of the positives (double or multiple sources) is 98%, and the estimated rate of dual (physically related) quasars is greater than 54%. A few dozen serendipitously found objects of interest are discussed in more detail, including known and new lensed images, planetary nebulae, young IR stars of peculiar morphology, and quasars with catastrophic redshift errors in SDSS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Extracting scattering amplitudes for arbitrary two-particle systems with one-particle left-hand cuts via lattice QCD

We derive a general formalism that relates the spectrum of two-particle systems in a finite volume to physical scattering amplitudes, taking into account the presence of any left-hand branch cuts due to single-particle exchanges. The method first relates the finite-volume spectrum to an infinite-volume short-range quantity, denoted ${\mathcal{M}}_0$, and then relates the latter to the physical scattering amplitudes via known integral equations. The derivation of both relations is performed using all-orders perturbation theory and is exact up to neglected exponentially suppressed volume dependence. The relations hold for arbitrary two-particle systems with any number of coupled channels, non-identical and non-degenerate particles, and any intrinsic spin.

algorithms↗
DOE OSTI2026

C–C Bond Formation during Electrochemical CO 2 Reduction on Pristine Cu(100) Unlikely to Involve Adsorbed CO at Any Potential

Formation of hydrocarbons containing two or more carbon atoms (C 2+ ) during heterogeneous electrochemical CO and CO 2 reduction (ECOR and ECO 2 R) only occurs, among pure metals, on Cu electrodes. Moreover, the activity and selectivity is facet dependent, with Cu(100) generally preferentially forming ethylene over methane. Previously, we found via quantum-mechanics-based modeling that, unlike standard density functional theory, more accurate correlated wavefunction methods predict that non-electroactive coupling pathways involving two adsorbed COs (*CO) or a *CO and a *COH to form C–C bonds on Cu(100) are kinetically inhibited, with the former also thermodynamically unfavorable. Here, we extend that embedded complete active space second order perturbation theory (ECASPT2) study, further showing that electrochemical coupling of two *COs to form an anionic dimer [OC*–*CO] (1+δ)– , followed by protonation to form [OC*–*COH] δ− , is not kinetically competitive with the reduction of *CO to *COH at relevant ECO/CO 2 R potentials. Our simulations therefore suggest that the ability of Cu(100) to electrochemically synthesize C 2+ molecules from CO and CO 2 is unlikely to be via *CO, at least on pristine Cu(100). Instead, hydrogenated CO species (*COH, *CH x OH, or *CH x ) are most likely to be the key intermediates in C–C bond formation.

Martirez, John Mark P. [Princeton Plasma Physics L↗
DOE OSTI2024

Excited state electronic structure of dimethyl disulfide involved in photodissociation at ∼200 nm

Dimethyl disulfide (DMDS), one of the smallest organic molecules with an S–S bond, serves as a model system for understanding photofragmentation in polypeptides and proteins. Prior studies of DMDS photodissociation excited at ∼266 nm and ∼248 nm have elucidated the mechanisms of S–S and C–S bond cleavage, which involve the lowest excited electronic states S 1 and S 2 . Far less is known about the dissociation mechanisms and electronic structure of relevant excited states of DMDS excited at ∼200 nm. Herein we present calculations of the electronic structure and properties of electronic states S 1 –S 6 accessed when DMDS is excited at ∼200 nm. Our analysis includes a comparison of theoretical and experimental UV spectra, as well as theoretically predicted one-dimensional cuts through the singlet and triplet potential energy surfaces along the S–S and C–S bond dissociation coordinates. Finally, we present calculations of spin–orbit coupling constants at the Franck–Condon geometry to assess the likelihood of ultrafast intersystem crossing. We show that choosing an accurate yet computationally efficient electronic structure method for calculating the S 0 –S 6 potential energy surfaces along relevant dissociation coordinates is challenging due to excited states with doubly excited character and/or mixed Rydberg-valence character. Our findings demonstrate that the extended multi-state complete active space second-order perturbation theory (XMS-CASPT2) balances this computational efficiency and accuracy, as it captures both the Rydberg character of states in the Franck–Condon region and multiconfigurational character toward the bond-dissociation limits. In conclusion, we compare the performance of XMS-CASPT2 to a new variant of equation of motion coupled cluster theory with single, double, and perturbative triple corrections, EOM-CCSD(T)(a)*, finding that EOM-CCSD(T)(a)* significantly improves the treatment of doubly excited states compared to EOM-CCSD, but struggles to quantitatively capture asymptotic energies along bond dissociation coordinates for these states.

Rishi, Varun↗
DOE OSTI2023

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Particle Swarm Optimisation for group structure optimization for radiotherapy shielding

Neutron transport simulations are ubiquitous in nuclear engineering because they allow one to model experimental systems and render a model platform for easy perturbation of experimental designs. In addition, simulations allow one to gain experimental insight without actually having to go through the trouble of building a physical experiment. Neutron transport simulations can be stochastic or deterministic based. Stochastic neutron transport simulations are typically simulated using the Monte Carlo method and yield very accurate solutions but are computationally expensive, while deterministic methods are typically faster but can be less accurate. Here we focus on optimizing the accuracy of deterministic neutron transport simulations for radiotherapy simulations. Deterministic neutron transport requires discretization of angle, energy, and space to appropriately analyze the system one is trying to model. Discretization of energy is challenging because of the highly variable neutron flux at certain neutron energies. Improper discretization of energy in the transport model can lead to erroneous results and therefore inaccurate interpretations of the solution. In this study, we evaluate Particle Swarm Optimization (PSO) as a mechanism for selecting optimal group structures for radiotherapy shielding. We tested the particle swarm optimization algorithm on radiotherapy shielding problems using Los Alamos National Laboratory's (LANL) main deterministic transport code PARTISN. Results show that the optimized energy group structures generated from the optimization algorithm outperformed LANL's standard energy group structures, and therefore demonstrate utility in using PSO to expedite computation times due to the increased accuracy obtained with a smaller but optimized group structure. (authors)

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Importance of charge self-consistency in first-principles description of strongly correlated systems

First-principles approaches have been successful in solving many-body Hamiltonians for real materials to an extent when correlations are weak or moderate. As the electronic correlations become stronger often embedding methods based on first-principles approaches are used to better treat the correlations by solving a suitably chosen many-body Hamiltonian with a higher level theory. The success of such embedding theories, often referred to as second-principles, is commonly measured by the quality of self-energy Σ which is either a function of energy or momentum or both. However, Σ should, in principle, also modify the electronic eigenfunctions and thus change the real space charge distribution. While such practices are not prevalent, some works that use embedding techniques do take into account these effects. In such cases, choice of partitioning, of the parameters defining the correlated Hamiltonian, of double-counting corrections, and the adequacy of low-level Hamiltonian hosting the correlated subspace hinder a systematic and unambiguous understanding of such effects. Further, for a large variety of correlated systems, strong correlations are largely confined to the charge sector. Then an adequate nonlocal low-order theory is important, and the high-order local correlations embedding contributes become redundant. Here we study the impact of charge self-consistency within two example cases, TiSe 2 and CrBr 3 , and show how real space charge re-distribution due to correlation effects taken into account within a first-principles Green’s function-based many-body perturbative approach is key in driving qualitative changes to the final electronic structure of these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Non-local slicing approaches for NNLO QCD in MCFM

We present the implementation of several processes at Next-to-Next-to Leading Order (NNLO) accuracy in QCD in the parton-level Monte Carlo program MCFM. The processes treated are pp → H, W ± , Z, W ± H, ZH, W ± γ, Zγ and γγ and, for the first time in the code, W + W - , W ± Z and ZZ. Decays of the unstable bosons are fully included, resulting in a flexible fully differential Monte Carlo code. The NNLO corrections have been calculated using two non-local slicing approaches, isolating the doubly unresolved region by cutting on the zero-jettiness, $\mathscr{T}_0$, or on q T , the transverse momentum of the colour singlet final-state particles. We find that for most, but not all processes the q T slicing method leads to smaller power corrections for equal computational burden.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗