Search NASA⌕ Search

SEARCH · Search NASA

Results for “Harmonic Constants”

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

Quantum and anharmonic effects in non-adiabatic transition state theory

Quantitative descriptions of non-adiabatic transition rates at intermediate temperatures are challenging due to the simultaneous importance of quantum and anharmonic effects. Here, in this paper, the interplay between quantum effects—for motion across or along the seam of crossing—and anharmonicity in the seam potential is considered within the weak coupling limit. The well-known expression for quantized 1-D motion across the seam (i.e., tunneling) in the linear terms approximation is derived in the thermal domain using the Lagrangian formalism, which is then applied to the case when tunneling is distributed along the seam of crossing (treating motion along the seam classically). For high-frequency quantum modes, a vibrationally adiabatic (VA) approach is developed that introduces to the non-adiabatic rate constant a factor associated with high-frequency wavefunction overlap; this approach treats the high-frequency motion along the seam quantum mechanically. To test these methodologies, the reaction N 2 O ↔ N 2 + O( 3 P) was chosen. CCSD(T)-F12b/cc-pVTZ-F12 explorations of the 3 A'- 1 A' seam of N 2 O revealed that seam anharmonicity has a strong effect on the rate constant (a factor of ~20 at 2000 K). Several quantum effects were found to be significant at intermediate/lower temperatures, including the quantum N–N vibration that was coupled with seam anharmonicity using the VA approach. Finally, a 1-D approximation to non-adiabatic instanton theory is presented to estimate the validity limit of the linear terms model at low temperatures (⁠~ 250 K for N 2 O). We recommend that the assumptions built into many statistical theories for non-adiabatic reactions—harmonic behavior, classical motion, linear terms, and weak coupling—should be verified on a case-by-case basis.

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

Imposing equilibrium on experimental 3-D stress fields using Hodge decomposition and FFT-based optimization

Here, we present a methodology to impose micromechanical constraints, i.e. stress equilibrium at grain and sub-grain scale, to an arbitrary (non-equilibrated) voxelized stress field obtained, for example, by means of synchrotron X-ray diffraction techniques. The method consists in finding the equilibrated stress field closest (in L 2 -norm sense) to the measured non-equilibrated stress field, via the solution of an optimization problem. The extraction of the divergence-free (equilibrated) part of a general (non-equilibrated) field is performed using the Hodge decomposition of a symmetric matrix field, which is the generalization of the Helmholtz decomposition of a vector field into the sum of an irrotational field and a solenoidal field. The combination of: a) the Euler–Lagrange equations that solve the optimization problem, and b) the Hodge decomposition, gives a differential expression that contains the bi-harmonic operator and two times the curl operator acting on the experimental stress field. These high-order derivatives can be efficiently performed in Fourier space. The method is applied to filter the non-equilibrated parts of a synthetic piecewise constant stress fields with a known ground truth, and stress fields in Gum Metal, a beta-Ti-based alloy measured in-situ using Diffraction Contrast Tomography (DCT). In both cases, the largest corrections were obtained near grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ab initio calculations of the rate of carrier trapping and release at dopant sites in NaI: Tl beyond the harmonic approximation

We present ab initio calculations of electron capture and release coefficients as a function of temperature for Tl dopants in the widely used scintillator NaI, which is a soft ionic crystal. The modelled capture and release events occur by transitions mediated by multiple phonon absorption and emission between states with significantly different local geometries around the trapping site. We demonstrate that such transitions are not well-described by the normal harmonic approximation to the nuclear dynamics. We avoid the harmonic approximation by numerically solving the vibrational Schrödinger equation for the motion along/of a single effective phonon coordinate. The localized trapped state is not correctly described by semilocal density functionals, so we employ a hybrid functional tuned to reproduce the band gap and lattice constant of the host material. Our calculations, which combine plane-wave and cluster electronic structure methods in a novel way, are in reasonable agreement with available experimental data for detrapping and predict an unusual temperature dependence for trapping rates.

Prange, Micah P.↗
DOE OSTI2026

Low-cost calculation and analysis of 2D IR spectra of model diiron trinitrosyl complexes in the NO stretch region with vibrational perturbation theory

Two-dimensional infrared spectroscopy offers unique capabilities for probing vibrational coupling in complex metal–ligand systems. In this paper, we combine two-dimensional infrared spectroscopy with vibrational perturbation theory to investigate vibrational coupling in a diiron trinitrosyl complex across three stable redox states. Although these systems are challenging for electronic structure methods, we demonstrate that key features of experimental 2D IR spectra can be accurately reproduced using reduced-dimensional anharmonic calculations with a small harmonic frequency scaling. Analysis reveals that N–O stretching modes maintain high locality across all redox states, with coupling patterns that directly reflect variations in Fe–N bond strength. Using curvilinear coordinate analysis, we demonstrate these differences result from systematic changes in cubic anharmonic force constants rather than mode delocalization. Our results establish N–O stretches as sensitive probes of metal–ligand bonding strength, expanding the toolkit for studying biologically relevant nitrosyl complexes.

Moran, Hayden A. [Texas A&M University, College St↗
DOE OSTI2025

Crystal Growth of Quaternary Rare Earth Selenosilicates by Using the Flux-Assisted Boron Chalcogen Mixture Method: Investigation of Their Magnetic and Optical Properties

A series of rare earth magnesium selenosilicates, RE3Mg 0.5 SiSe 7 (RE = Ce, Pr, Nd, Sm, Gd, Tb, Dy) were obtained as single crystal using the flux assisted boron chalcogen mixture (BCM) method. The structures of the crystals were determined by single-crystal X-ray diffraction. The RE3Mg 0.5 SiSe 7 series crystallizes in the hexagonal crystal system in the space group P6 3 . Polycrystalline powders were synthesized to perform physical property measurements. Magnetic measurements over the 2–300 K temperature range reveal that Ce 3 Mg 0.5 SiSe 7 and Gd 3 Mg 0.5 SiSe 7 exhibit paramagnetic behavior with negative Weiss constants (θ W = −14.50, θ W = −6.13 K, respectively). The optical properties of RE 3 Mg 0.5 SiSe 7 (RE = Ce, Pr, Nd, Sm, Gd) were measured by ultraviolet–visible (UV–vis) diffuse reflectance. Density functional theory (DFT) electronic structure calculations were performed. Furthermore, a second harmonic generation measurement was performed on a polycrystalline powder of Ce 3 Mg 0.5 SiSe 7 and was found to be SHG active with an efficiency of 0.11 times the standard potassium dihydrogen phosphate (KDP).

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

Energy Exascale Earth System Model v2.0.1

First patch release of v2.0.0 Changes since v2.0.0 [Important change] Fix ocean threading bug seen in debug cases on Chrysalis with Intel 20.0.4. Was introduced around time of v2.0.0 tag. Does not change v2.0.0 answers on Chrysalis because those didn't use threading or debugging. [EAM] Add semi-lagrangian tracer transport for theta-l (F90 and C++), add new algorithm for finding tropopause, add DSCREAM to allow v2 and SCREAM settings in same code such as adjust_ps [EAM-MMF] 60L default, allow C++ back end of RRTMGP (EAM too). [EAMxx] add nu-top functionality, fix forcing functor, add ttype9 and dcmip2012 tests 2.1, 2.2, and 3 HOMME: remove obsolete remap algs, option to specify dynamics alg indep of tracer, new sponge layer, add imex tests [ELM] Add topography-based subgrid (topounits), add FATES-ELM Nitro., Phos. and CH4 coupling, add land-use ts for NARRM, add lulc for SSP3 RCP7, Fix nutrient fertilization exp test and carbon isotope flux, Fix xactive lnd dry deposition, add lake water storage option, fix plant hydraulics 2d params, fix carbon budget calc, fix soil nutrient conc. bug, fix mosart dam bug, add test for new ELM, MOSART features, fix bug in O3 dry dep stomatal resistances, fix plant hydraulics restart BFB error, update mkmapdata. [MOSART] fix bug for reading the latitude from an unstructured input file, fix oversat in bubble test. [MPAS-ocean] Add CFC11, CFC12 tracers, add 2D spherical transport tests, fix del4 tracer mixing, add MARBL ocean tracer mixing, modify harmonic analysis options, add GPU port of vmix routines, fix calc of ML-averaged BV freq. [MPAS-seaice] Change extents of initial polar disks for oRRS18to6v3 grid, fix ice BGC with MARBL, update spherical test cases, fix DON coupling, Remove Cf from sea ice constants. [MPAS-landice] add CRYO1850-4xCO2 compset [CIME] add GCP, ANL GCE, Spock, Perlmutter, deprecate config_compilers.xml, fix and clean-up cmake macros, fix slurm bindings, refactor CIME internal testing, cleanup SCORPIO perf data, allow position independent compset naming, [also] update v2 benchmarking suite, extend e3sm_prod with throughput and memory checks

E3SM Project, DOE↗
DOE OSTI2020

A predictive strain-gradient model with no undetermined constants or length scales

A general meso–scale (GM) crystal plasticity (CP) model was developed that accounts for lower-order (strain hardening) and higher-order (internal stress) effects of geometrically necessary dislocations (GNDs). It is predictive: no arbitrary parameters or length scales were invoked and no ad hoc numerical techniques were employed. It uses general stress field equations for GND content and a novel harmonization technique to enforce consistency of elastic long-range singular defect fields with applied elastic-plastic fields. The model facilitates implementation in commercial finite element programs without requiring special elements, special boundary conditions, or access to element shape functions. GM simulations confirmed, with improved accuracy, previously published predictions of the Hall-Petch effect, Bauschinger effect, and anelasticity. Previously unpredicted phenomena were also predicted: anelasticity and hysteresis for single Ta crystals and strain-hardening stagnation. The internal stresses (higher-order effect) dominate at large length scales, while at small length scales, the GND density hardening (lower-order effect) dominates. As a result, GM predicts that strain heterogeneity and consequent GND internal stresses are important factors in anelasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fluctuation–dissipation relation for a quantum Brownian oscillator in a parametrically squeezed thermal field

In this paper we study the nonequilibrium evolution of a quantum Brownian oscillator, modeling the internal degree of freedom of a harmonic atom or an Unruh–DeWitt detector, coupled to a nonequilibrium and nonstationary quantum field bath and inquire whether a fluctuation–dissipation relation (FDR) can exist after/if it approaches equilibration. This is a nontrivial issue because a squeezed field bath cannot reach equilibration and yet, as this work shows, the system oscillator indeed can, which is a necessary condition for FDRs. We discuss three different settings: (A) The bath field essentially remains in a squeezed thermal state throughout, whose squeeze parameter is a mode- and time-independent constant. This situation is often encountered in quantum optics and quantum thermodynamics. (B) The bath field is initially in a thermal state, but is subjected to a parametric process leading to mode- and time-dependent squeezing. This scenario is encountered in cosmology and dynamical Casimir effects. The squeezing in the bath in both types of processes will affect the oscillator’s nonequilibrium evolution. We show that at late times it approaches equilibration and this stationarity condition warrants the existence of a FDR. The trait of squeezing is marked by the oscillator’s effective equilibrium temperature, and the proportionality factor in the FDR is only related to the stationary component of the noise kernel of the bath field. Setting (C) is more subtle: A finite system–bath coupling strength can set the oscillator in a squeezed state even though the bath field is stationary and does not engage in any parametric process. The squeezing of the system in this case is in general time-dependent but becomes constant when the internal dynamics is fully relaxed. We begin with comments on the broad range of physical processes involving squeezed thermal baths and end with some remarks on the significance of FDRs in capturing the essence of quantum backreaction in nonequilibrium and stochastic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Colloquium : Physical properties of group-IV monochalcogenide monolayers

The state-of-the-art knowledge of ferroelectric and ferroelastic group-IV monochalcogenide monolayers is surveyed. These semiconductors feature remarkable structural and mechanical properties, such as a switchable in-plane spontaneous polarization, soft elastic constants, structural degeneracies, and thermally driven two-dimensional structural transformations. Additionally, these 2D materials display selective valley excitations, valley Hall effects, and persistent spin helix behavior. After a description of their Raman spectra, a discussion of optical properties arising from their lack of centrosymmetry (such as an unusually strong second-harmonic intensity, large bulk photovoltaic effects, photostriction, and tunable exciton binding energies) is provided as well. Here, the physical properties observed in these materials originate from (correlate with) their intrinsic and switchable electric polarization, and the physical behavior hereby reviewed could be of use in nonvolatile memory, valleytronic, spintronic, and optoelectronic devices: these 2D multiferroics enrich and diversify the 2D material toolbox.

2D ferroelectrics↗
DOE OSTI2024

Structure Sensitive Reaction Kinetics of Chiral Molecules on Intrinsically Chiral Surfaces

Enantiospecific heterogeneous catalysis utilizes chiral surfaces to resolve enantiomers via structure sensitive surface chemistry. The catalyst design challenge is the identification of chiral surface structures that maximize enantiospecificity. Herein, we develop data driven models for the enantiospecificity of tartaric acid reactions on chiral Cu(hkl) R&S surfaces. Measurements of enantiospecific rate constants were obtained by using curved Cu(hkl) R&S surfaces that enable kinetic measurements on hundreds of chiral surface orientations. One model uses feature vectors derived from generalized coordination numbers to capture the local structure around Cu atoms exposed by the Cu(hkl) R&S surfaces. The second model introduces the use of chiral cubic harmonic functions to capture the symmetry constraints of the face-centered cubic Cu structure. The model using 58 generalized coordination numbers has a fitting error similar to that of the model using only 5 cubic harmonic functions. The two models predict maxima in the enantiospecificity on surfaces with very similar surface orientations. The models developed in this work are applicable for any enantiospecific reaction happening on any chiral material with a cubic lattice structure, opening the way to understanding the surface structure sensitivity of the enantiospecific reaction kinetics.

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

Spatio-temporal electric field distributions in an atmospheric plasma jet impinging on a microchannel array surface

The electric field distribution in the ionization waves (IWs) propagating over a microchannel array dielectric surface, with the channels either empty or filled with distilled water, is measured by ps electric field induced second harmonic generation. The surface IW is initiated by the atmospheric pressure N 2 –Ar plasma jet impinging on the surface vertically and powered by ns pulse discharge bursts. The results show that the electric field inside the microchannels, specifically its horizontal component, is enhanced by up to a factor of 2. The field enhancement region is localized within the channels. The vertical electric field inside the channels lags in time compared to the field measured at the ridges, indicating the transient reversal of the IW propagation direction across the channels (toward the jet). This is consistent with the phase-locked plasma emission images and confirmed by the kinetic modeling predictions, which show that the IW ‘jumps’ over the empty channels and propagates into the channels only after the jump between the adjacent ridges. When the channels are filled with water, the wave speed increases by up to 50%, due to the higher effective dielectric constant of the surface. No evidence of a significant electric field enhancement near the dielectric surface (ceramic or water) has been detected, within the spatial resolution of the present diagnostic, ~100 μ m.

plasmas and liquids↗
DOE OSTI2024

Measuring signatures in photon angular spectra to distinguish nonlinear Compton scattering models

The collision of a high-energy electron beam with a laser pulse may be used to study radiation reaction and nonlinear Compton scattering among many other processes in strong-field quantum electrodynamics. Predictions from simulation and theory for these interactions rely on a number of approximations and assumptions that have not been experimentally tested. Here, experimentally measurable signatures are identified that might be able to distinguish between radiation reaction models, i.e., classical or quantum, or between the local constant field and local monochromatic approximations used to calculate the properties of the nonlinear Compton process. These signatures are considered through Monte Carlo simulations of various experimental conditions that are relevant to today's laser facilities. Potential detection schemes for measuring the signatures are proposed. We find that single-photon counting of keV photons to resolve harmonics and scintillator-based detection of MeV photons may allow us to validate nonlinear Compton scattering models and radiation reaction models respectively. This will require electron beams with divergence angles less than 2 mrad and less than 20% energy spread.

Russell, Brandon K↗
DOE OSTI2020

ABINIT: Overview and focus on selected capabilities

ABINIT is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe–Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the “temperature-dependent effective potential” approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which ABINIT relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The ABINIT DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library libpaw. ABINIT has strong links with many other software projects that are briefly mentioned.

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

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

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

Islands in the fluid: islands are common in cosmology

We discuss the possibility of entanglement islands in cosmological spacetimes with a general perfect fluid with an equation of state w. We find that flat universes with time-symmetric slices where the Hubble parameter vanishes always have islands on that slice. We then move away from such slices, considering still universes with a general perfect fluid. Under the local thermal equilibrium assumption, the comoving entropy density s c is constant. As a result, the conditions for an island become an inequality between the energy density (or Hubble parameter) and the temperature at some time of normalization. The consequences are that islands can exist for practically all fluids that are not radiation, i.e. w ≠ 1/3. We also discuss the ramifications of our results for universes with spatial curvature. Finally, we show that islands occur in the Simple Harmonic Universe model which has no classical singularity at the background level, in contrast to all previous examples where islands occurred only in space-times with singularities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Axion magnetic resonance: A novel enhancement in axion-photon conversion

We identify a new resonance, axion magnetic resonance, that can greatly enhance the conversion rate between axions and photons. A series of axion search experiments rely on converting them into photons inside a constant magnetic field background. A common bottleneck of such experiments is the conversion amplitude being suppressed by the axion mass when m a ≳ 10 − 4 eV . We point out that a spatial or temporal variation in the magnetic field can cancel the difference between the photon dispersion relation and that of the axion, hence greatly enhancing the conversion probability. We demonstrate that the enhancement can be achieved by both a helical magnetic field profile and a harmonic oscillation of the magnitude. Our approach can extend the projected ALPS II reach in the axion-photon coupling ( g a γ ) by two orders of magnitude at m a = 10 − 3 eV with moderate assumptions (see .) Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dynamic Shaping of Grid Response of Multi-Machine Multi-Inverter Systems Through Grid-Forming IBRs: Preprint

We consider the problem of controlling the frequency response of weakly-coupled multi-machine multi-inverter low-inertia power systems via grid-forming inverter-based resources (IBRs). In contrast to existing methods, our approach relies on dividing the larger system into multiple strongly-coupled subsystems, without ignoring either the underlying network or approximating the subsystem response as an aggregate harmonic mean model. Rather, through a structured clustering and recursive dynamic shaping approach, the frequency response of the overall system to load perturbations is shaped appropriately. We demonstrate the proposed approach for a three-node triangular configuration and a small-scale radial network. Furthermore, for small-scale radial microgrids, we demonstrate the ability of IBRs to tune the effective transfer functions of synchronous machines. This enables us to relax the uniform turbine time-constant assumptions and widen the scope of existing synchronization results for proportionally heterogeneous machines.

frequency-shaping control↗
DOE OSTI2025

Development of a wide bandwidth heterodyne dispersion interferometer for electron density measurement of atmospheric pressure plasmas

One of the challenges of electron density measurements of an atmospheric pressure plasma (APP) with a laser interferometer is the significant and unwanted phase shift caused by changes in the neutral gas density. These unwanted phase shifts can be mitigated and plasma density measured using an interferometer architecture called a dispersion interferometer (DI). A DI is composed of two nonlinear orientation patterned GaAs crystals for frequency doubling and measures the phase shift induced by a plasma in the interference signal between two second-harmonic beams. Measurement of plasma dynamics or a short plasma pulse in less than a millisecond is enabled with heterodyne detection in a DI using an acousto-optic cell with a frequency of 40 MHz. This heterodyne DI (HDI) is targeted to measure APPs in an electron density range of 10 20 –10 24 m −3 . Finally, the HDI achieves a line-integrated density resolution of 2 × 10 15 m −2 (a phase resolution of 0.005°) with a time constant of 1 μs using ensemble averaging techniques.

atmospheric pressure plasma↗