Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Zero point energy”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Topological and magnetic properties of a noncollinear spin state on a honeycomb lattice in a magnetic field

Here, this paper studies the topological and magnetic properties of a noncollinear spin state on a honeycomb lattice that evolves from coplanar to ferromagnetic with a magnetic field applied along the z axis. The coplanar state is stabilized by nearest-neighbor ferromagnetic interactions, single-ion anisotropy along z, and DzyaloshinskiiMoriya interactions between next-nearest-neighbor sites. Below the critical field H$_c$ that aligns the spins, the magnetic unit cell contains six sites and the spin dynamics contains six magnon modes. Although the classical energy is degenerate with respect to the twist angle φ between nearest-neighbor spins, the dependence of the free energy on φ at low temperatures is dominated by the magnon zero-point energy, which contains extremum at φ = πl/3 for integer l. The only unique ground states GS(φ) have l = 0 or 1. For H < H$'_c$, the zero-point energy has minima at even l and the ground state is GS(0); for H$'_c$ < H < H$'_c$, the zero-point energy has minima at odd l and the ground state is GS(π/3). In GS(0), the magnon density of states exhibits five distinct topological phases with increasing field associated with the opening and closing of energy gaps between two or three magnonic bands. While the Berry curvature vanishes for the coplanar φ = 0 phase in zero field, the Berry curvature and Chern numbers exhibit signatures of the five topological phases below H$'_c$. Whereas the Berry curvature and Chern number are sensitive to changes in the magnon density of states within GS(π/3), the inelastic spectrum S(k,ω) is sensitive to changes in the intensity of the magnon modes in the different magnetic phases GS(0) and GS(π/3) rather than the five topological phases within GS(π/3).

Fishman, Randy S. [Oak Ridge National Laboratory (↗
DOE OSTI2025

Quantum Sieving for Isotopic Separations of Gases Using Porous Materials─30 Years of Progress

Here, this paper reviews theoretical and experimental efforts to establish microporous materials that can be used to separate isotopologues of molecular gases such as H 2 and D 2 . Emphasis is placed on use of simplified models to highlight the quantum phenomena that make these separations possible. In equilibrium adsorption, differences in zero point energy in the adsorbed states of molecules favor binding of heaver isotopologues (e.g. D 2 relative to H 2 ). Experimental data showing this effect was reported as early as 1933, but the theoretical work of Beenakkers et al. in 1995 [Chem. Phys. Lett. 232 (1995) 379-382] spurred modern efforts to develop materials with strong so-called quantum sieving. In porous materials where the transition state for molecular diffusion is more strongly confined than the energy minima associated with equilibrium adsorption, differences in zero point energies between these two sites can lead to isotopic differences in molecular diffusivities. This effect favors diffusion of heavier species (e.g. D 2 ) relative to lighter species (e.g. H 2 ). This so-called kinetic quantum sieving has been observed experimentally in porous carbons and in porous organic cage materials. We show that quantum tunneling, which favors hopping of lighter species across energy barriers, diminishes the strength of kinetic quantum sieving but that it appears to make only a small contribution to the net molecular diffusivities in many porous materials of practical interest.

Sholl, David S. [Oak Ridge National Laboratory (OR↗
DOE OSTI2025

The reaction NH 2 + CH 2 O: Kinetic measurements and computational studies

The reaction of amidogen with formaldehyde is relevant to astrochemistry and couples nitrogen and carbon chemistry in flames. The first measurements of the temperature dependence of the rate constant k were made, over 409–643 K, and are summarized as k = 7.1 × 10 −12 exp(−18.2 kJ mol −1 /RT) cm 3 molecule −1 s −1 with a confidence interval of ±20%. These results address a controversy over the presence of a significant barrier in the addition/elimination pathway leading to the production of formamide between theoretical models that include the zero-point energy of all modes at the transition state and a model that excludes transitional modes. The latter would lead to a negligible barrier, which is inconsistent with the experiments. Contrary to earlier claims, CCSDT(Q)-corrected energies combined with transition state theory yield quantitative accord with the measurements.

Astrochemistry↗
DOE OSTI2026

Unimolecular Dynamics of Partially Deuterated Hydroperoxyalkyl Intermediates (•QOOD) in Cyclohexane and Cyclopentane Oxidation

The hydroperoxyalkyl radicals (•QOOH) formed in cyclohexane and cyclopentane oxidation with carbon radical center at the β site relative to the hydroperoxy group have recently been observed through their infrared fingerprint and time- and energy-resolved unimolecular dissociation dynamics to hydroxyl (OH) radical and cyclic ether products. Partial deuteration shifts the IR transitions associated with OH/OD motion to lower frequencies, providing a new window on the IR spectroscopy and unimolecular dynamics of these transient intermediates. Specifically, the overtone OD stretch (2ν OD ) has been identified at 5215.0 cm –1 and 5216.0 cm –1 for the β-QOOD intermediates in cyclohexane and cyclopentane oxidation, respectively. Unimolecular decay rates of the •QOOD intermediates at these energies have been obtained from the time-resolved appearance of OD products. The experimental rates are compared with statistical microcanonical rates evaluated using RRKM theory, including heavy-atom tunneling associated with simultaneous O–O bond elongation and C–O–O angle contraction along the reaction pathway. The experimental rates provide a further test of the previously determined transition state (TS) barriers, which were computed utilizing a benchmarking approach that builds on higher-level reference calculations for the smaller ethane oxidation system. Here, the experimental rate measurements for both •QOOD intermediates agree well with the computed rates after accounting for small changes in zero-point energy and a minor empirical adjustment of the TS barriers in both systems.

Ethers↗
DOE OSTI2025

Reliable p K a Prediction through Efficient Incorporation of Anharmonicity within the Nuclear–Electronic Orbital Framework

Accurate pK a prediction is critical for understanding chemical reactivity and molecular properties across a wide range of applications. Computational approaches usually invoke a harmonic treatment of the vibrational modes for zero-point energies, as well as thermal and entropic contributions. Herein, we present a general protocol for relative pK a prediction that incorporates the significant anharmonic effects using nuclear–electronic orbital (NEO) theory. This protocol is validated against experimental data for a range of molecules in acetonitrile, including protonated nitrogen bases, nitrophenols, anilines, and diamines, as well as cobalt electrocatalysts. For simple acids, the NEO approach offers only a slight improvement over conventional density functional theory with the standard harmonic vibrational treatment, whereas for hydrogen-bonded acids, the NEO approach offers more significantly improved performance at a comparable computational cost. This accessible methodology provides a practical route for accurate pKa prediction in challenging systems and is extendable to related thermodynamic properties such as hydricities and proton-coupled redox potentials.

Density functional theory↗
DOE OSTI2024

Exploring direct photodetachment and photodissociation–photodetachment dynamics of platinum iodide anions (PtI n - , n = 2–5) using cryogenic photoelectron spectroscopy

The direct photodetachment and two-photon photodissociation–photodetachment processes of a series of PtI n - (n = 2–5) anions were systematically studied using cryogenic anion photoelectron spectroscopy and first-principles electronic structure calculations. The adiabatic/vertical detachment energies (ADEs/VDEs) of these anions were determined from their 193 nm photoelectron (PE) spectra, i.e., 3.54/3.63, 4.04/4.09, 4.33/4.36, and 4.37/4.41 eV for n = 2–5, respectively, and well reproduced by B3LYP-D3(BJ)/aug-cc-pVTZ-pp calculations. As the coordination number increases, the electron affinity (EA) of PtI n • (n = 2–5) neutrals (equivalent to the corresponding anion’s ADE) gradually increases, exceeding the EA of Cl at n = 3 and exhibiting superhalogen characteristics for n ≥ 3. Meanwhile, the ground state transition contributed from detaching electrons in the highest occupied molecular orbital gradually evolves from the central metal Pt to the iodine ligands. For the PtI 3 - anion, besides one-photon direct detachment, four distinct two-photon photodissociation–photodetachment channels were identified, and the competition between them was discussed.

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

The Remarkable I 2 O 3 Molecule: A New View from Theory

Atmospheric iodine chemistry has garnered increasing attention as a result of increased iodine emissions. A key subset of this chemistry involves iodine oxides (I 2 O 2–5 ), which serve as precursors to particle formation. Among these, I 2 O 3 is the simplest iodine oxide involved in particle formation, but it has remained undetected in the atmosphere. Previous theoretical studies have characterized this peculiar molecule, primarily using energies to refine geometries obtained at low levels of theory. Due to the reemerging interest in I 2 O 3 , this study presents geometries optimized at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory─marking the first instance, to the best of our knowledge, where this system has been studied exclusively with CCSD(T). Harmonic vibrational frequencies were computed at the same level of theory. Final energetics were obtained using the very high level CCSDT(Q) method with basis sets up to quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) and extrapolated to the CBS limit to yield CCSDT(Q)/CBS//CCSD(T)/aug-cc-pwCVTZ-PP energies. These energies include harmonic zero-point vibrational energy corrections and scalar relativistic energy corrections. Additionally, this study discovers new isomers along the I 2 O 3 potential energy surface, a novel contribution to the field. The performance of different computational methods and DFT functionals commonly used in atmospheric chemistry is also assessed relative to high-level theoretical methods.

basis sets↗
DOE OSTI2025

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗
DOE OSTI2024

A windowed mean trajectory approximation for condensed phase dynamics

We propose a trajectory-based quasi-classical method for approximating dynamics in condensed phase systems. Building upon the previously developed optimized mean trajectory approximation that has been used to compute linear and nonlinear spectra, we borrow some ideas from filtering trajectory methods to obtain a novel semiclassical method for the dynamical propagation of density matrices. This new approximation is tested rigorously against standard multistate electronic models, spin-boson models, and models of the Fenna–Matthews–Olson complex. For dissipative systems, the current method is significantly better or as good as many other semiclassical methods available, especially at low temperatures and for off-diagonal density matrix elements, whereas for scattering models, the current method bears similar limitations as mean-field propagation schemes. All results are tested against the numerically exact hierarchical equations of motion method. In conclusion, the new method shows excellent agreement across various parameter regimes with numerically exact results, highlighting the robustness and accuracy of our approach.

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

The DESI DR1 Peculiar Velocity Survey: Global Zero-point and H 0 Constraints

The Dark Energy Spectroscopic Instrument (DESI) in its first Data Release (DR1) already provides more than 100,000 galaxies with relative distance measurements. The primary purpose of this paper is to perform the calibration of the zero-point for the DESI Fundamental Plane and Tully–Fisher relations, which allows us to measure the Hubble constant, H 0 . This sample has a lower statistical uncertainty than any previously used to measure H 0 , and we investigate the systematic uncertainties in absolute calibration that could limit the accuracy of that measurement. We improve upon the DESI Early Data Release Fundamental Plane H 0 measurement by (a) using a group catalog to increase the number of calibrator galaxies and (b) investigating alternative calibrators in the nearby Universe. Our baseline measurement calibrates to the SH0ES/Pantheon+ type Ia supernovae, and finds H 0 = 73.7 ± 0.06 (stat.) ± 1.1 (syst.) km s −1 Mpc −1 . Calibrating to surface brightness fluctuation distances yields a similar H 0 . We explore measurements using other calibrators, but these are currently less precise since the overlap with DESI peculiar velocity tracers is much smaller. In future data releases with an even larger peculiar velocity sample, we plan to calibrate directly to Cepheids and the tip of the red giant branch, which will enable the uncertainty to decrease towards a percent-level measurement of H 0 . This will provide an alternative to supernovae as the Hubble flow sample for H 0 measurements.

Carr, Anthony [Korea Astronomy and Space Science I↗
DOE OSTI2026

Chirality reversal at finite magnetic impurity strength and local signatures of a topological phase transition

Here, we study the honeycomb lattice with a single magnetic impurity modeled by adding imaginary next-nearest-neighbor hopping 𝑖⁢ℎ on a single hexagon. This Haldane defect gives a topological mass term to the gapless Dirac cones and generates chirality. For a small density of defects, Neehus et al. [Phys. Rev. Lett. 135, 126604 (2025)] found that the system's chirality reverses at a critical ℎ 𝑐 ≈ 0.95 associated with an unexpected tricritical point of Dirac fermions at zero defect density. We investigate this zero-density limit by analyzing a single defect and computing two experimentally relevant measures of chirality: (1) orbital magnetization via local Chern marker, a bulk probe of all occupied states; and (2) electronic currents of low-energy states. Both probes show a chirality reversal at a critical ℎ 𝑐 ≈ 0.9–1.0. Motivated by this consistency, we propose a defect-scale toy model whose low-energy states reverse their chirality at ℎ$^{'}_{c}$ ≈ 0.87. Remarkably, the same pair of zero-energy bound states also generates the critical point ℎ 𝑐 in the full impurity projected T-matrix. Our results show how the chirality reversal produced by an impurity can be observed either in local probes or in the global topology, and suggest a possible role of the microscopic defect structure at the critical point.

Chern insulators↗
DOE OSTI2025

Anharmonic phonons, thermal expansion, and nuclear quantum effects in Zn

The anharmonic behavior of phonons and thermal expansion of hexagonal zinc were studied from 15 to 690 K by inelastic neutron scattering (INS) and ab initio simulations. Phonon spectra were measured for Q-points covering the full Brillouin zone, giving the phonon density of states (DOS) and dispersions along high-symmetry directions. The dispersions were sharp at 15 K, but diffuse intensity was observed at energies above them. The dispersions broadened with temperature T and the diffuse intensity became stronger. This diffuse intensity appeared in all INS measurements and simulations, except for classical molecular dynamics at 15 K. The temperature-dependent effective potential (TDEP) method, which included the nuclear quantum effect from zero-point vibrational dynamics, was used to calculate the free energy and thermal expansion. For T < 100 K nuclear quantum effects were important for obtaining the correct negative thermal expansion, and path integral molecular dynamics (PIMD) was particularly effective for obtaining the negative thermal expansion in the basal plane. A Heisenberg-Langevin model for interacting phonons coupled to a thermal bath was able to reproduce the shape and intensity of the diffuse spectral features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

GridSTIX

SF-25-112 Grid-STIX is a comprehensive extension of the STIX (Structured Threat Information Expression) 2.1 ontology specifically designed for electrical grid cybersecurity applications. This ontology provides a standardized, machine-readable framework for modeling grid assets, operational technology devices, threats, vulnerabilities, supply chain risks, and security relationships in electrical power systems. ## Key Features - **Comprehensive Grid Coverage**: Physical assets, OT devices, grid components, sensors, and energy storage systems - **Zero Trust Architecture**: Policy decision points, enforcement points, trust brokers, and continuous monitoring - **AMI Infrastructure**: Advanced metering networks, head-end systems, mesh gateways, and MDM systems - **Advanced Security Modeling**: Attack patterns, vulnerabilities, mitigations, and supply chain risks - **Critical Grid Relationships**: Power flow, protection, control, and synchronization relationships - **Supply Chain Security**: Supplier modeling, country of origin tracking, and risk assessment - **Protocol Support**: DNP3, Modbus, IEC 61850, IEC 60870-5-104, OPC-UA, and IEEE standards - **Python Code Generation**: Automated STIX-compliant Python class generation from ontologies - **Interactive Visualization**: Enhanced HTML network graphs with grid-specific categorization - **STIX 2.1 Compliance**: Full compatibility with STIX threat intelligence ecosystem

Blakely, Benjamin [Argonne National Laboratory (AN↗
DOE OSTI2024

Exotic edge states of C 3 high-fold fermions in honeycomb lattices

A generalization of the graphene honeycomb model to the case where each site in the honeycomb lattice contains a n -fold degenerate set of eigenstates of the C 3 symmetry has been recently proposed to describe several systems, including triangulene crystals and photonic lattices. These generalized honeycomb models are defined by ( n a , n b ) , the number of C 3 eigenstates in the a and b sites of the unit cell, resulting in n a + n b bands. Thus, the (1,1) case gives the coventional honeycomb model that describes the two low-energy bands in graphene. Generalizations, such as (2,1), (2,2), and (3,3) display several nontrivial features, such as coexisting graphenelike Dirac cones with flat bands, both at zero and finite energy, as well as robust degeneracy points where a flat band and a parabolic band meet at the Γ point. Here we explore the edge states of this class of crystals, using as reference triangulene crystals, and we find several types of edge states absent in the conventional (1,1) honeycomb case, associated to the nontrivial features of the two-dimensional bands of the high-fold case. First, we find dispersive edge states associated to the finite-energy flat bands, that occur both at the armchair and zigzag termination. Second, in the case of noncentrosymmetric triangulene crystals that lead to a S = 1 Dirac band, we have a bonding-antibonding pair of dispersive edge states, localized in the same edge so that their energy splitting is reduced as their localization increases, opposite to the conventional behavior of pairs of states localized in opposite edges. Third, for the (3,3) case, that hosts a gap separating a pair of flat conduction and valence bands, we find nondispersive edge states with E = 0 in all edge terminations. Published by the American Physical Society 2024

Madail, L. (ORCID:0000000208457748)↗
DOE OSTI2025

The Unitarity-Limit Expansion for Two Nucleons with Perturbative Pions: Digest and Ideas

Theorists love nontrivial fixed points. In the Unitarity Limit, the NN 𝑆-wave binding energies are zero, the scattering lengths infinite, Physics is universal, i.e. insensitive to details of the interactions, and observables display richer symmetries, namely invariance under both scaling and Wigner’s combined SU(4) transformation of spin and isospin. In “Pionless” EFT, both are explicitly but weakly broken and hence perturbative in the Unitarity Window (phase shifts 45° ≲ δ(k) ≲ 135°, i.e. momenta k≈mπ). This Unitarity Expansion provides strong hints that Nuclear Physics resides indeed in a sweet spot: bound weakly enough to be insensitive to the details of the nuclear interaction; and therefore interacting strongly enough that the NN scattering lengths are perturbatively close to the Unitarity Limit. In this paradigm change, NN details are less important than NNN interactions to explain the complexity and patterns of the nuclear chart. This presentation is a digest of the first quantitative exploration of corrections to this picture when pions are included [1] (see there for a more comprehensive list of references). Since the pion mass and decay constant introduce dimensionful scales in the NN system, they explicitly break the symmetries of the Unitarity fixed point. In χEFT, these symmetries must therefore be hidden and instead be classified as emergent.

Griesshammer, Harald W. [The George Washington Uni↗
DOE OSTI2025

Superfluid weight cross-over and critical temperature enhancement in singular flat bands

Nonanalytic Bloch eigenstates at isolated band degeneracy points exhibit singular behavior in the quantum metric. Here, a description of superfluid weight for zero-energy flat bands in proximity to other high-energy bands is presented, where they together form a singular band gap system. When the singular band gap closes, the geometric and conventional contributions to the superfluid weight as a function of the superconducting gap exhibit different cross-over behaviors. The scaling behavior of superfluid weight with the band gap is studied in detail, and the effect on the Berezinskii–Kosterlitz–Thouless transition temperature is explored. It is found that tuning the singular band gap provides a unique mechanism for enhancing the supercurrent and critical temperature of two-dimensional superconductors.

Science & Technology - Other Topics↗
DOE OSTI2025

Spontaneous Localization at a Potential Saddle Point from Edge State Reconstruction in a Quantum Hall Point Contact

Quantum point contacts (QPCs) are an essential component in mesoscopic devices. Here, in this work, we study the transmission of quantum Hall edge modes through a gate-defined QPC in monolayer graphene. We observe resonant tunneling peaks and a nonlinear conductance pattern characteristic of Coulomb-blockaded localized states. The in-plane electric polarizability reveals the states are localized at a classically unstable electrostatic saddle point. We explain this unexpected finding within a self-consistent Thomas-Fermi model, finding that localization of a zero-dimensional state at the saddle point is favored whenever the applied confinement potential is sufficiently soft compared to the Coulomb energy. Our results provide a direct demonstration of Coulomb-driven reconstruction at the boundary of a quantum Hall system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Long-term decarbonization impacts on residential energy security across income groups and US states

Abstract The impact of a transition to a net-zero economy on the residential energy sector across diverse income groups in the US remains uncertain. Here, we employ an integrated human-Earth system model, incorporating an expanded set of ten income groups in the residential energy sector, to examine the distributional impacts of long-term decarbonization scenarios on residential energy security at the state level through 2050. We use multiple metrics of energy security, including energy burden, energy satiation gap, and the distribution of energy service across income groups. Our findings show that the net-zero decarbonization scenarios affect residential energy security differently across income groups, with low-to-mid-income groups experiencing larger negative impacts on the dimensions studied here. Comparatively, climate change impact on residential energy security is minor through 2050 based on our model outcomes. Specifically, the net-zero decarbonization scenarios lead to increased energy burden across all income groups and states in 2050, where the lowest (highest) income group in each state shows an average of 0.6 (0.2) percentage point increase in energy burden, relative to the business-as-usual in 2050. The distribution of energy service consumption across income groups is also slightly more skewed under these scenarios. As incomes grow across all deciles in the future, residential energy security generally improves through 2050. Targeted interventions could mitigate the disproportionate impacts that some groups could incur under a transition to a net-zero economy.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗