Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “excited state 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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An Efficient Three-Dimensional CFD-Based Numerical Wave Tank for a Wave Energy Converter in Extreme Irregular Waves: Preprint

A numerical wave tank approach for computational fluid dynamics (CFD) modelling of an extreme irregular seastate is presented. The technique couples a potential flow solution with a CFD solver for more efficient numerical predictions. This method has recently become attractive both for the research community and the industry working with offshore structures. The model is used to determine the response of a submerged pressure differential wave energy converter (WEC) in a fully nonlinear irregular wave condition using the high fidelity CFD code, STAR-CCM+. Potential flow based numerical models are commonly used to predict motions and performance of wave energy converters. Wave kinematics can deviate from potential flow predictions for extreme wave conditions; the excitation loads on an absorber can also be increasingly influenced by viscous effects, not captured in engineering level models. In these extreme conditions, a Reynolds-averaged Navier-Stokes CFD model can better predict motions and loads for a WEC. Long time series with varying random seed numbers can be used to identify singular extreme wave events from a stochastic irregular sea state. This approach simulates a more realistic wave series for a given sea state than a regular wave or a focused wave. However, it is computationally infeasible to run these long time series for three-dimensional (3D) CFD simulations. In this work, two-dimensional (2D) CFD simulations with a long domain allow the full development of an extreme nonlinear wave condition. The results are used to identify extreme events from a 50-year storm condition for the PacWave site off the coast of Oregon. A relatively short time window including this extreme event is then mapped to a 3D simulation using a user defined wave methodology. Convergence studies for domain length, wave forcing lengths, and time before the extreme event were conducted.

cfd↗
DOE OSTI2025

All-optical control of charge-trapping defects in rare-earth doped oxides

Abstract Charge-trapping defects in crystalline solids play important roles in applications ranging from microelectronics, optical storage, sensing and quantum technologies. On one hand, depleting trapped charges in the host matrix reduces charge noise and enhances coherence of solid-state quantum emitters. On the other hand, stable charge traps can enable high-density optical storage systems. Here we report all-optical control of charge-trapping defects via optical charge trapping (OCT) spectroscopy of a rare-earth ion doped oxide (Y 2 O 3 ). Charge trapping is realized by low intensity optical excitation in the 200–375 nm range. Charge detrapping or depletion is carried out by optically stimulated luminescence (OSL) under 532 nm stimulation. Using a Pr-doped Y 2 O 3 polycrystalline ceramic host matrix, we observe charging pathways via the inter-band optical absorption of Y 2 O 3 and via the 4f-5d transitions of Pr 3+ . We demonstrate effective control of the density of trapped charges within the Y 2 O 3 matrix at ambient environment. These results point to a viable method for controlling the local charge environment in rare-earth doped crystals via all-optical means, and pave the way for further development of efficient optical storage technologies with ultrahigh storage capacity, as well as for the localized control of quantum coherence in rare-earth doped solids.

França, Leonardo V. S. [Pritzker School of Molecul↗
DOE OSTI2021

Tracing the Mott-Hubbard transition in one-dimensional Hubbard models without Umklapp scattering

Here we apply the density-matrix renormalization group (DMRG) method to a one-dimensional Hubbard model that lacks Umklapp scattering and thus provides an ideal case to study the Mott-Hubbard transition analytically and numerically. The model has a linear dispersion and displays a metal-to-insulator transition when the Hubbard interaction U equals the bandwidth U c = W, where the single-particle gap opens linearly, Δ(U ≥ W) = U - W. The simple nature of the elementary excitations permits us to determine numerically with high accuracy the critical interaction strength and the gap function in the thermodynamic limit. The jump discontinuity of the momentum distribution $n_k$ at the Fermi wave number $k_F$ = 0 cannot be used to locate accurately U c from finite-size systems. However, the slope of $n_k$ at the band edges $k_B$ = ±π, reveals the formation of a single-particle bound state which can be used to determine U c reliably from $n_k$ using accurate finite-size data.

1-dimensional systems↗
DOE OSTI2021

Giant phonon anomalies in the proximate Kitaev quantum spin liquid α-RuCl 3

The Kitaev quantum spin liquid epitomizes an entangled topological state, for which two flavors of fractionalized low-energy excitations are predicted: the itinerant Majorana fermion and the Z 2 gauge flux. Here, it was proposed recently that fingerprints of fractional excitations are encoded in the phonon spectra of Kitaev quantum spin liquids through a novel fractional-excitation-phonon coupling. Here, we detect anomalous phonon effects in α-RuCl 3 using inelastic X-ray scattering with meV resolution. At high temperature, we discover interlaced optical phonons intercepting a transverse acoustic phonon between 3 and 7 meV. Upon decreasing temperature, the optical phonons display a large intensity enhancement near the Kitaev energy, J K ~8 meV, that coincides with a giant acoustic phonon softening near the Z 2 gauge flux energy scale. These phonon anomalies signify the coupling of phonon and Kitaev magnetic excitations in α-RuCl 3 and demonstrates a proof-of-principle method to detect anomalous excitations in topological quantum materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Energy conservation in real-time nuclear–electronic orbital Ehrenfest dynamics

Real-time nuclear–electronic orbital Ehrenfest (RT-NEO-Ehrenfest) dynamics methods provide a first-principles approach for describing nonadiabatic molecular processes with nuclear quantum effects. For an efficient description of proton transfer within RT-NEO-Ehrenfest dynamics, the basis function center associated with the quantum proton can be allowed to move classically. Here, this traveling proton basis (TPB) approach effectively captures proton quantum dynamics, although its energy conservation behavior is not yet fully satisfactory. Two recently proposed TPB approaches, in principle, conserve the extended energy, which includes both the system energy and the kinetic energy associated with the proton basis function center. Herein, a thermostatted TPB approach is proposed to improve the conservation of the system energy, excluding the kinetic energy associated with the proton basis function center. In this approach, the quantum proton dynamics are modulated by dynamically rescaling the proton momentum operator to maintain the system energy conservation. With the excited-state intramolecular proton transfer of o-hydroxybenzaldehyde as an example, this approach is shown to significantly improve the system energy conservation while preserving the accuracy of the quantum proton dynamics as achieved in the original TPB approach.

Ehrenfest dynamics↗
DOE OSTI2023

An Efficient Three-Dimensional CFD-Based Numerical Wave Tank for a Wave Energy Converter in Extreme Irregular Waves

A numerical wave tank approach for computational fluid dynamics (CFD) modelling of an extreme irregular seastate is presented. The technique couples a potential flow solution with a CFD solver for more efficient numerical predictions. This method has recently become attractive both for the research community and the industry working with offshore structures. The model is used to determine the response of a submerged pressure differential wave energy converter (WEC) in a fully nonlinear irregular wave condition using the high-fidelity CFD code, STAR-CCM+. Potential flow based numerical models are commonly used to predict motions and performance of wave energy converters. Wave kinematics can deviate from potential flow predictions for extreme wave conditions; the excitation loads on an absorber can also be increasingly influenced by viscous effects, not predicted by potential flow engineering level models. In these extreme conditions, a Reynolds-averaged Navier-Stokes CFD model can better predict motions and loads for a WEC. Long time series with varying random seed numbers can be used to identify singular extreme wave events from a stochastic irregular sea state. This approach simulates a more realistic wave series for a given sea state than a regular wave or a focused wave. However, it is computationally infeasible to run these long time series for three-dimensional (3D) CFD simulations. In this work, two-dimensional (2D) CFD simulations with a long domain allow the full development of an extreme nonlinear wave condition. The results are used to identify extreme events from a 50-year storm condition for the PacWave site off the coast of Oregon. A relatively short time window including this extreme event is then mapped to a 3D simulation using a user defined wave methodology. Convergence studies for domain length, wave forcing lengths, and time before the extreme event were conducted.

CFD↗
DOE OSTI2023

A Fast Butterfly-Compressed Hadamard–Babich Integrator for High-Frequency Helmholtz Equations in Inhomogeneous Media with Arbitrary Sources

Here we present a butterfly-compressed representation of the Hadamard-Babich (HB) ansatz for the Green's function of the high-frequency Helmholtz equation in smooth inhomogeneous media. For a computational domain discretized with Nv discretization cells, the proposed algorithm first solves and tabulates the phase and HB coefficients via eikonal and transport equations with observation points and point sources located at the Chebyshev nodes using a set of much coarser computation grids, and then butterfly compresses the resulting HB interactions from all Nv cell centers to each other. The overall CPU time and memory requirement scale as O(Nv log2 Nv) for any bounded two-dimensional (2D) domains with arbitrary excitation sources. A direct extension of this scheme to bounded 3D domains yields an O(Nv4/3) CPU complexity, which can be further reduced to quasi-linear complexities with proposed remedies. The scheme can also efficiently handle scattering problems involving inclusions in inhomogeneous media. Although the current construction of our HB integrator does not accommodate caustics, the resulting HB integrator itself can be applied to certain sources, such as concave-shaped sources, to produce caustic effects. Compared to finite-difference frequency domain methods, the proposed HB integrator is free of numerical dispersion and requires fewer discretization points per wavelength. As a result, it can solve wave propagation problems well beyond the capability of existing solvers. Remarkably, the proposed scheme can accurately model wave propagation in 2D domains with 640 wavelengths per direction and in 3D domains with 54 wavelengths per direction on a state-of-the-art supercomputer at Lawrence Berkeley National Laboratory.

Hadamard--Babich ansatz↗
DOE OSTI2023

Multidimensional Simulations of Core Convection

The cores of main sequence intermediate- and high-mass stars are convective. Mixing at the radiative–convective boundary, waves excited by the convection, and magnetic fields generated by convective dynamos all influence the main sequence and post-main sequence evolution of these stars. These effects must be understood to accurately model the structure and evolution of intermediate- and high-mass stars. Unfortunately, there are many challenges in simulating core convection due to the wide range of temporal and spatial scales, as well as many important physics effects. In this review, we describe the latest numerical strategies to address these challenges. We then describe the latest state-of-the-art simulations of core convection, summarizing their main findings. These simulations have led to important insights into many of the processes associated with core convection. Two outstanding problems with multidimensional simulations are, 1. it is not always straightforward to extrapolate from simulation parameters to the parameters of real stars; and 2. simulations using different methods sometimes appear to arrive at contradictory results. To address these issues, next generation simulations of core convection must address how their results depend on stellar luminosity, dimensionality, and turbulence intensity. Furthermore, code comparison projects will be essential to establish robust parameterizations that will become the new standard in stellar modeling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Fixed-phase diffusion Monte Carlo study of activation gap and skyrmion excitations of a ν = 1 system in the presence of charged impurities

The discrepancy between the theoretically calculated and experimentally measured activation gaps in quantum Hall effect has long been a puzzle. We revisit this issue in the context of the v = 1 quantum Hall state, while also incorporating the skyrmion physics. We find that the finite width and the Landau level mixing effects are not sufficient to explain the observed activation gap. We further show that the presence of charged impurities located adjacent to the quantum well can cause a significant reduction in the activation gap, while also causing a suppression of the skyrmion size.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Gas-phase synthesis of naphthalene through an unconventional thermal alkyne–alkene [2 + 2] cycloaddition mechanism

Exotic cycloaddition entrance channels were discovered for the bimolecular gas-phase reactions of the phenylethynyl radical (C 6 H 5 CC, X 2 A 1 ) with ethylene-d 4 (C 2 D 4 ) and propylene (C 3 H 6 ) as explored under single-collision conditions utilizing the crossed molecular beams technique combined with electronic structure and statistical calculations. Connecting the concepts of barrierless entrance channels, excited states, and facile non-photochemically activated cycloadditions, the reaction pathway features an unconventional thermal [2 + 2] cycloaddition forming a four-membered ring collision complex followed by multiple isomerizations prior to unimolecular decomposition via atomic hydrogen loss to (un)substituted naphthalenes—naphthalene-d 4 (C 10 H 4 D 4 ) and 1-/2-methylnaphthalene (C 11 H 10 ). The small energy gap between the singly-occupied a 1 highest occupied molecular orbital (HOMO) with a σ-character and the underlying doubly-occupied b 1 molecular orbital with a π-character allows a facile promotion of an electron. This in turn enables a versatile low-temperature reactivity of phenylethynyl, where the end-on and side-on barrierless approaches of ethylene are due to its interaction with the σ and π orbitals, respectively, thus suggesting this mechanism as a possible method for tuning substituents in polycyclic aromatic hydrocarbon (PAH) formation and highlighting its versatility as a probe of fundamental carbon chemistry via counterintuitive cycloaddition reactions under single-collision conditions.

Goettl, Shane J. [University of Hawaii at Manoa, H↗
DOE OSTI2025

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗
DOE OSTI2026

Computational schemes for the Magnus expansion of the in-medium similarity renormalization group

The in-medium similarity renormalization group (IMSRG) is a popular many-body method used for computations of nuclei. It solves the many-body Schrödinger equation through a continuous unitary transformation of the many-body Hamiltonian. The IMSRG transformation is typically truncated at the normal-ordered two-body level, the IMSRG(2), but recently several approaches have been developed to capture the effects of normal-ordered three-body operators, the IMSRG(3). In particular, a factorized approximation to the IMSRG(3) proposes to capture the leading effects of three-body operators at the same computational cost as the IMSRG(2) approximation. This approach often employs an approximate scheme for solving the IMSRG equations, the so-called hunter-gatherer scheme. In this work, I study the uncertainty associated with this scheme. I find that the hunter-gatherer scheme differs by up to 7MeV for ground-state energies and 0.5MeV for excitation energies from standard IMSRG(2) approaches. These differences are in some cases comparable to the expected size of IMSRG(3) corrections.

39 ≤ A ≤ 58↗
DOE OSTI2021

Electronic structure of reduced CeO 2 (111) surfaces interacting with hydrogen as revealed through electron energy loss spectroscopy in comparison with theoretical investigations

Based on both, ab-initio cluster calculations as well as periodic density functional based loss function calculations, we have assigned the origins of the valence electron excitation regime in electron energy loss spectra of well-ordered ceria films in (111) orientation at various states of reduction as well as after exposing the reduced films to hydrogen from the gas phase (Li et al, Angew. Chem. Int. Ed., 58 (2019) 14686–14693). The explicit calculation of intensity distributions using the dipole approximation allow us to draw conclusions about the nature of oxygen vacancies, which occur upon reduction and how those interact with hydrogen to form hydride species. Herein, we also reports a brief discussion of vibrational excitations of hydrogen loaded ceria and corroborates previous interpretations, based on inelastic neutron scattering.

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

How transparent is graphene? A surface science perspective on remote epitaxy

Remote epitaxy is the synthesis of a single crystalline film on a graphene-covered substrate, where the film adopts epitaxial registry to the substrate as if the graphene is transparent. Despite many exciting applications for flexible electronics, strain engineering, and heterogeneous integration, an understanding of the fundamental synthesis mechanisms remains elusive. Here we offer a perspective on the synthesis mechanisms, focusing on the foundational assumption of graphene transparency. We identify challenges for quantifying the strength of the remote substrate potential that permeates through graphene, and propose Fourier and beating analysis as a bias-free method for decomposing the lattice potential contributions from the substrate, from graphene, and from surface reconstructions, each at different frequencies. We highlight the importance of graphene-induced reconstructions on epitaxial templating, drawing comparison to moiré epitaxy. We highlight the role of the remote potential in tuning surface diffusion and adatom kinetics on graphene, which are crucial for navigating the competition between remote epitaxy and defect-seeded mechanisms like pinhole epitaxy. In light of this weak remote potential, we re-evaluate the current state-of-the-art experimental evidence, highlighting why it remains challenging to experimentally validate a ‘remote’ epitaxy mechanism that cannot be explained by alternatives, such as pinhole-seeded epitaxy or serial van der Waals epitaxy. We end with one experimental example that, to out knowledge, cannot be explained by competing mechanisms: a different long-range epitaxial relationship for GdPtSb films grown on graphene/sapphire, compared to direct epitaxy on sapphire. We suggest for future experiments that directly measure the remote potential and impact of tuneable growth kinetics.

epitaxy↗
DOE OSTI2026

Quantifying uncertainty in physics-based predictions of rare-isotope production cross sections via Bayesian-inspired model averaging across nuclear mass tables

Accurate prediction of fragmentation cross sections is essential for rare-isotope beam production, planning new-isotope searches, and designing experiments to study the most exotic regions of the nuclear chart. However, existing reaction models and phenomenological cross-section parametrizations often exhibit significant deviations over broad regions of mass and charge. In this work, a Bayesian-inspired model-averaging framework is developed to combine abrasion-ablation (AA) calculations based on multiple nuclear mass tables into a single statistically weighted estimate. For the calibrated systems, the model weights are assigned empirically according to the relative quality of fit to measured cross sections, thereby reducing systematic model bias while preserving the underlying physics content of the AA description. The weights are constrained using proton-rich fragmentation data for the 78 Kr and 124 Xe projectiles. The resulting parameter trends are then propagated to the 92 Mo and 144 Sm systems through a controlled scaling procedure. In the present implementation, the excitation-energy prescription is fixed, while the averaging is performed across nuclear-mass inputs; the framework provides both weighted cross sections and associated uncertainty estimates. Applied to proton-rich fragmentation, the present approach provides a practical basis for interpolation and limited extrapolation in regions relevant to rare-isotope production. The resulting predictions are used to assess the production of very proton-rich nuclei, and candidate new isotopes are discussed.

Bayesian methods↗
DOE OSTI2025

Evolution of Mass Spectrometers for High m / z Biological Ion Formation, Transmission, Analysis and Detection: A Personal Perspective

Mass spectrometry (MS) has become an essential tool in virtually all academic, pharmaceutical, and biopharmaceutical analytical laboratories. The specialized and bespoke area of MS research and application of high m / z ion (> m / z 6000 and high mass, >150 kDa) formation, transmission, analysis, and detection is a relatively new area of focus for MS that has seen dramatic acceleration in interest over the last two decades. Herein we delve into this exciting aspect of MS, discussing how MS instrumentation has been refined and evolved for native-MS analysis. We cover the early groundbreaking experiments showing high m / z ion formation, transmission, and preservation of protein structure in the gas phase. Additionally, we discuss specific instrument optimizations and modifications that have advanced high m / z ion generation, transmission, analysis, and detection, contributing to the research area known as gas-phase structural biology. Native-MS sample introduction methods, emerging technologies, and future perspectives are also examined. Finally, we share personal opinions, observations, and experiences that are new to the community or previously unpublished.

collisions↗
DOE OSTI2021

Photophysical properties of zero-dimensional perovskites studied by PBE0 and GW+BSE methods

Due to the exceptional photophysical properties and great potential applications in solid state lighting, zero-dimensional halide perovskites with broadband emitting have recently attracted much attention. In this work, to reveal the origin of the broadband absorption and emission spectra of 0D (C 4 N 2 H 14 Br) 4 SnBr x I 6-x (x = 3), the ground state and exciton dynamics properties have been studied by hybrid PBE0 functional and the GW+BSE method. Different self-trapped excitons (STEs) are obtained, and their emission energies by PBE0 are consistent well with experimental data. By considering the vibrational broadening of emission, our calculated full width at half maximum of emission accords well with the experiment. The absorption of free exciton and STEs is also calculated by considering the quasiparticle and excitonic effect. Our results show that these light-induced trapped states can be observed by experiments, as confirmed in other similar systems by the femtosecond transient absorption spectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Fundamental data for modeling electron-induced processes in plasma remediation of perfluoroalkyl substances

Plasma treatment of per- and polyfluoroalkyl substances (PFAS) contaminated water is a potentially energy efficient remediation method. In this treatment, an atmospheric pressure plasma interacts with surface-resident PFAS molecules. Developing a reaction mechanism and modeling of plasma–PFAS interactions requires fundamental data for electron–molecule reactions. In this paper, we present results of electron scattering calculations, potential energy landscapes and their implications for plasma modelling of a dielectric barrier discharge in PFAS contaminated gases, a first step towards modelling of plasma–water–PFAS intereactions. It is found that the plasma degradation of PFAS is dominated by dissociative electron attachment with the importance of other contributing processes varying depending on the molecule. All molecules posses a large number of shape resonances – transient negative ion states – from near-threshold up to ionization threshold. These states lie in the region of the most probable electron energies in the plasma (4–5 eV) and consequently are expected to further enhance the fragmentation dynamics in both dissociative attachment and dissociative excitation.

54 ENVIRONMENTAL SCIENCES↗