Search NASA⌕ Search

SEARCH · Search NASA

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 55 records · Page 3

Adiabatic approximation and Aharonov-Casher bands in twisted homobilayer transition metal dichalcogenides

Topological flat moiré bands with nearly ideal quantum geometry have been identified in homobilayer transition metal dichalcogenide moiré superlattices, and are thought to be crucial for understanding the fractional Chern insulating states recently observed therein. Previous study proposed viewing the system using an adiabatic approximation that replaces the position dependence of the layer spinor with a nonuniform periodic effective magnetic field. When the local zero-point kinetic energy of this magnetic field cancels identically against that of an effective Zeeman energy, a Bloch-band version of Aharonov-Casher zero-energy modes, which we refer to as Aharonov-Casher band, emerges leading to ideal quantum geometry. Here, we critically examine the validity of the adiabatic approximation and identify the parameter regimes under which Aharonov-Casher bands emerge. We show that the adiabatic approximation is accurate for a wide range of parameters including those realized in experiments. Furthermore, we show that while the cancellation leading to the emergence of Aharonov-Casher bands is generally not possible beyond the leading Fourier harmonic, the leading harmonic is the dominant term in the Fourier expansions of the zero-point kinetic energy and Zeeman energy. Finally, the leading harmonic expansion accurately captures the trend of the bandwidth and quantum geometry, although it may fail to quantitatively reproduce more detailed information about the bands such as the Berry curvature distribution.

36 MATERIALS SCIENCE↗

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. Here, this study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Convergent ab initio analysis of the multi-channel HOBr + H reaction

High-level potential energy surfaces for three reactions of hypobromous acid with atomic hydrogen were computed at the CCSDTQ/CBS//CCSDT(Q)/complete basis set level of theory. Focal point analysis was utilized to extrapolate energies and gradients for energetics and optimizations, respectively. The H attack at Br and subsequent Br–O cleavage were found to proceed barrierlessly. The slightly submerged transition state lies −0.2 kcal mol−1 lower in energy than the reactants and produces OH and HBr. The two other studied reaction paths are the radical substitution to produce H2O and Br with a 4.0 kcal mol−1 barrier and the abstraction at hydrogen to produce BrO and H2 with an 11.2 kcal mol−1 barrier. The final product energies lie −37.2, −67.9, and −7.3 kcal mol−1 lower in energy than reactants, HOBr + H, for the sets of products OH + HBr, H2O + Br, and H2 + BrO, respectively. Additive corrections computed for the final energetics, particularly the zero-point vibrational energies and spin–orbit corrections, significantly impacted the final stationary point energies, with corrections up to 6.2 kcal mol−1.

Chemistry↗

Analysis of the classical trajectory treatment of photon dynamics for polaritonic phenomena

Simulating photon dynamics in strong light–matter coupling situations via classical trajectories is proving to be powerful and practical. Here, in this work, we analyze the performance of the approach through the lens of the exact factorization approach. Since the exact factorization enables a rigorous definition of the potentials driving the photonic motion, it allows us to identify that the underestimation of photon number and intensities observed in earlier work is primarily due to an inadequate accounting of light–matter correlation in the classical Ehrenfest force rather than errors from treating the photons quasiclassically per se. The latter becomes problematic when the number of photons per mode begins to exceed a half.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transformation of an Embedded Five-Membered Ring in Polycyclic Aromatic Hydrocarbons via the Hydrogen-Abstraction–Acetylene-Addition Mechanism: A Theoretical Study

Five-membered rings are constituents of many polycyclic aromatic hydrocarbons (PAHs) and their presence on the edges of large PAHs has been repeatedly observed experimentally. However, modern kinetic combustion models often do not consider the growth of PAHs through the transformation of the five-membered rings. Here, in connection with the above, we carried out a theoretical study of the mechanism of Hydrogen-Abstraction-aCetylene-Addition (HACA) transformation of an embedded five-membered ring on the armchair PAH edge to a six-membered ring, considering cyclopenta[d,e,f]phenanthrene (4,5-methylenephenanthrene) as a prototype system for this process. The potential energy surface for the reactions of cyclopenta[d,e,f]phenanthrenyl radicals produced by direct H abstractions from cyclopenta[d,e,f]phenanthrene with acetylene has been compiled at the G3(MP2,CC)//B3LYP/6-311G(d,p) level of theory including zero-point vibrational energy corrections. The computed energies and molecular parameters were then used to solve the Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) in order to calculate reaction rate at various pressures and temperatures, which were fitted to modified Arrhenius equation for further kinetic modeling. The results show that the HACA transformation of the embedded five-membered ring to a six-membered ring is possible, albeit slow. The most viable reaction mechanism involves the R2 + C 2 H 2 reaction, where the acetylene molecules adds to a σ-radical in the six-membered ring adjacent to the five-membered ring via a low entrance barrier. The predominant product of R2 + C 2 H 2 is predicted to be 3 ethynyl 4H cyclopenta[def]phenanthrene Pr5 via immediate H elimination from the initial addition complex. Next, Pr5 undergoes H-assisted isomerization to 4aH-pentaleno[4,3,2,1-cdef]phenanthrene Pr4, and the latter adds an H atom eventually forming the 1-pyrenylmethyl radical Pr3: R2 + C 2 H 2 ⇌ 3 ethynyl 4H cyclopenta[def]phenanthrene (Pr5) + H ⇌ 4aH-pentaleno[4,3,2,1-cdef]phenanthrene (Pr4) + H Pr5 + H ⇌ Pr4 + H Pr4 + H → 1-pyrenylmethyl (Pr3) This HACA sequence may be competitive with the methyl radical addition to the R1 radical formed by H abstraction from the CH 2 group in the five-membered ring of cyclopenta[d,e,f]phenanthrene, which provides a pathway to pyrene following two H atom losses. Relative contributions of the two mechanisms of the five- to six-membered ring transformation would strongly depend on the branching ratios of the R1 and R2 radicals produced by the H abstractions and the available concentration of C 2 H 2 vs. CH 3 and hence differ in different flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The water trimer reaction OH + (H 2 O) 3 → (H 2 O) 2 OH + H 2 O

All important stationary points on the potential energy surface (PES) for the reaction OH + (H 2 O) 3 → (H 2 O) 2 OH + H 2 O have been fully optimized using the “gold standard” CCSD(T) method with the large Dunning correlation-consistent cc-pVQZ basis sets. Three types of pathways were found. For the pathway without hydrogen abstraction, the barrier height of the transition state (TS1) is predicted to lie 5.9 kcal mol -1 below the reactants. The two major complexes (H 2 O) 3 ∙∙∙OH (CP1 and CP2a) are found to lie 6.3 and 11.0 kcal mol -1 , respectively, below the reactants [OH + (H 2 O) 3 ]. For one of the H-abstraction pathways the lowest classical barrier height is predicted to be much higher, 6.1 kcal mol -1 (TS2a) above the reactants. For the other H-abstraction pathway the barrier height is even higher, 15.0 (TS3) kcal mol -1 . Vibrational frequencies and the zero-point vibrational energies connected to the PES are also reported. Fianlly, the energy barriers for the H-abstraction pathways are compared with those for the OH + (H 2 O) 2 and OH + H 2 O reactions, and the effects of the third water on the energetics are usually minor (0.2 kcal mol -1 ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring the Tl 2H 2 potential energy surface: A comparative analysis with group 13 systems and experiment

Abstract Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet TlH potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug‐cc‐pwCVQZ‐PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug‐cc‐pwCVTZ‐PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple‐zeta cardinality (aug‐cc‐pwCV5Z‐PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero‐point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 MH (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist.

Chemistry↗

Energetic and Spectroscopic Properties of Astrophysically Relevant MgC 4 H Radicals Using High-Level Ab Initio Calculations

Considering the importance of magnesium-bearing hydrocarbon molecules (MgC n H; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC 4 H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2 Σ + ), is the global minimum geometry on the MgC 4 H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ~ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ 3 -magnesacycloprop-2-ene-2,3-diyl (2, 2 A 1 , μ = 2.54 D), is 146 kJ mol –1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anisotropic surface state in a topological semimetal candidate Ta 3 SiTe 6

Using high-resolution angle-resolved photoemission spectroscopy, we studied the surface state of Ta 3 SiTe 6 , a candidate of topological semimetal protected by nonsymmorphic symmetry. Through photon-energy dependent measurements, the surface state near the Fermi level was detected. By determining the band dispersions, we found that the surface state is topological trivial. Around the surface Brillouin zone center, the surface bands disperse linear-like, forming a band crossing (zero energy gap) point at about 0.2 eV below the Fermi level. In Ta 3 SiTe 6 's rectangular surface Brillouin zone, the surface bands only cross the Fermi level along one axis, forming two hole-like Fermi pockets near surface Brillouin zone boundaries. Although the observed surface state is not topological, we proposed that its anisotropic Fermi surface topography could be potentially useful for future applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

Two-atom-thin topological crystalline insulators lacking out of plane inversion symmetry

A two-dimensional topological crystalline insulator (TCI) with a single unit cell (u.c.) thickness is demonstrated here. To that end, one first shows that tetragonal (C 4 in-plane) symmetry is not a necessary condition for the creation of zero-energy metallic surface states on TCI slabs of finite-thicknesses, because zero-energy states persist even as all the in-plane rotational symmetries—furnishing topological protection—are completely removed. In other words, zero-energy levels on the model are not due to (nor are they protected by) topology. Furthermore, effective two-fold energy degeneracies taking place at few discrete k-points away from zero energy in the bulk Hamiltonian—that are topologically protected—persist at the u.c. thickness limit. Here, the chiral nature of the bulk TCI Hamiltonian permits creating a $2\times 2$ square Hamiltonian, whose topological properties remarkably hold invariant at both the bulk and at the single u.c. thickness limits. The identical topological characterization for bulk and u.c.-thick phases is further guaranteed by a calculation involving Pfaffians. This way, a two-atom-thick TCI is deployed hereby, in a demonstration of a topological phase that holds both in the bulk, and in two dimensions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Many-Body Perspective of Nuclear Quantum Effects in Aqueous Clusters

Nuclear quantum effects play an important role in the structure and thermodynamics of aqueous systems. Here, by performing a many-body expansion with nuclear-electronic orbital (NEO) theory, we show that proton quantization can give rise to significant energetic contributions for many-body interactions spanning several molecules in single-point energy calculations of water clusters. Although zero-point motion produces a large increase in energy at the one-body level, nuclear quantum effects serve to stabilize higher-order molecular interactions. These results are significant because they demonstrate that nuclear quantum effects play a nontrivial role in many-body interactions of aqueous systems. Our approach also provides a pathway for incorporating nuclear quantum effects into water potential energy surfaces. The NEO approach is advantageous for many-body expansion analyses because it includes nuclear quantum effects directly in the energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Topological Mott transition in a two band model of spinless fermions with on-site Coulomb repulsion

Highlights: • The effect of the on-site Coulomb repulsion and the spin–orbit Rashba coupling lead to the 2D topological semimetal. • The ground state phase diagram is calculated. • The topological semimetal phase is characterized by the zero energy Majorana states. In the framework of mean field approach, we study topological Mott transition in a two band model of spinless fermions on a square lattice at half filling. We consider the combined effect of the on-site Coulomb repulsion and the spin–orbit Rashba coupling. The ground state phase diagram is calculated as a function of the strength of the spin–orbit Rashba coupling and the Coulomb repulsion. The spin–orbit Rashba coupling leads to a distinct phase of matter, the topological semimetal. We study a new type of phase transition between the non-topological insulator and topological semimetal states. Topological phase state is characterized by the zero energy Majorana states, which there are in defined region of the wave vectors and are localized at the boundaries of the sample. The region of existence of the zero energy Majorana states tends to zero at the point of the Mott phase transition. The zero energy Majorana states are dispersionless (they can be considered as flat bands), the Chern number and Hall conductance are equal to zero (note in two dimensional model).

36 MATERIALS SCIENCE↗

Digital pulse analysis for fast neutron recoil spectroscopy with a 4 He scintillation detector

Helium-4-based fast neutron scintillation detectors are an attractive alternative to pulse-shape discrimination-capable organic scintillators for fast neutron detection and spectroscopy, as the response of the detectors to gamma rays is intrinsically limited to low energy deposition. Consequently, the neutron recoil distribution can be measured with these detectors without the need for pulse shape analysis. In this work, the response of an Arktis S670 4 He scintillation detector to D-D, D-T, and 252 Cf neutrons was measured. The detector has a unique construction and readout mechanism, with multiple output channels observing the same scintillation event, and an analysis method was developed to aggregate the outputs from all channels into a single list. The D-D and D-T neutron responses were used to perform a two-point energy calibration, which yielded a near-zero intercept, suggesting that the 4 He scintillation medium behaves linearly to a higher energy than previously reported, and that a two-point calibration is sufficient for nuclear recoil energies below 9 MeV. As a result, the detector was measured to have 16.7-ns FWHM time resolution when using the developed custom analysis, a reduction of 4.9 ns when compared to the conventional pulse analysis.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

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↗

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↗