Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hamiltonian”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

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

Stochastic exciton-scattering theory of optical line shapes: Renormalized many-body contributions

Spectral line shapes provide a window into the local environment coupled to a quantum transition in the condensed phase. In this paper, we build upon a stochastic model to account for non-stationary background processes produced by broad-band pulsed laser stimulation, as distinguished from those for stationary phonon bath. In particular, we consider the contribution of pair-fluctuations arising from the full bosonic many-body Hamiltonian within a mean-field approximation, treating the coupling to the system as a stochastic noise term. Herein, using the Itô transformation, we consider two limiting cases for our model, which lead to a connection between the observed spectral fluctuations and the spectral density of the environment. In the first case, we consider a Brownian environment and show that this produces spectral dynamics that relax to form dressed excitonic states and recover an Anderson–Kubo-like form for the spectral correlations. In the second case, we assume that the spectrum is Anderson–Kubo like and invert to determine the corresponding background. Using the Jensen inequality, we obtain an upper limit for the spectral density for the background. The results presented here provide the technical tools for applying the stochastic model to a broad range of problems.

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

Permutation invariant polynomial neural network based diabatic ansatz for the ( E + A ) × ( e + a ) Jahn–Teller and Pseudo-Jahn–Teller systems

In this work, the permutation invariant polynomial neural network (PIP-NN) approach is employed to construct a quasi-diabatic Hamiltonian for system with non-Abelian symmetries. It provides a flexible and compact NN-based diabatic ansatz from the related approach of Williams, Eisfeld, and co-workers. The example of $H^3_+$ is studied, which is an (E + A) × (e + a) Jahn–Teller and Pseudo-Jahn–Teller system. The PIP-NN diabatic ansatz is based on the symmetric polynomial expansion of Viel and Eisfeld, the coefficients of which are expressed with neural network functions that take permutation-invariant polynomials as input. This PIP-NN-based diabatic ansatz not only preserves the correct symmetry but also provides functional flexibility to accurately reproduce ab initio electronic structure data, thus resulting in excellent fits. Further, the adiabatic energies, energy gradients, and derivative couplings are well reproduced. A good description of the local topology of the conical intersection seam is also achieved. Therefore, this diabatic ansatz completes the PIP-NN based representation of DPEM with correct symmetries and will enable us to diabatize even more complicated systems with complex symmetries.

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

High-resolution CH stretch spectroscopy of jet-cooled cyclopentyl radical: First insights into equilibrium structure, out-of-plane puckering, and IVR dynamics

First, high-resolution sub-Doppler infrared spectroscopic results for cyclopentyl radical (C 5 H 9 ) are reported on the α-CH stretch fundamental with suppression of spectral congestion achieved by adiabatic cooling to T rot ≈ 19(4) K in a slit jet expansion. Surprisingly, cyclopentyl radical exhibits a rotationally assignable infrared spectrum, despite 3N – 6 = 36 vibrational modes and an upper vibrational state density (ρ ≈ 40–90 #/cm –1 ) in the critical regime (ρ ≈ 100 #/cm –1 ) necessary for onset of intramolecular vibrational relaxation (IVR) dynamics. Such high-resolution data for cyclopentyl radical permit detailed fits to a rigid-rotor asymmetric top Hamiltonian, initial structural information for ground and vibrationally excited states, and opportunities for detailed comparison with theoretical predictions. Specifically, high level ab initio calculations at the coupled-cluster singles, doubles, and perturbative triples (CCSD(T))/ANO0, 1 level are used to calculate an out-of-plane bending potential, which reveals a C 2 symmetry double minimum 1D energy surface over a C 2v transition state. Here, the inversion barrier [V barrier ≈ 3.7(1) kcal/mol] is much larger than the effective moment of inertia for out-of-plane bending, resulting in localization of the cyclopentyl wavefunction near its C 2 symmetry equilibrium geometry and tunneling splittings for the ground state too small (<1 MHz) to be resolved under sub-Doppler slit jet conditions. The persistence of fully resolved high-resolution infrared spectroscopy for such large cyclic polyatomic radicals at high vibrational state densities suggests a “deceleration” of IVR for a cycloalkane ring topology, much as low frequency torsion/methyl rotation degrees of freedom have demonstrated a corresponding “acceleration” of IVR processes in linear hydrocarbons.

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

Evolution of the marker distribution in gyrokinetic $δf$ particle-in-cell simulations

The evolution of the particle weight in a δf particle-in-cell simulation depends on the marker distribution that can evolve in a turbulent field due to turbulent diffusion. When Monte Carlo methods are used to implement the test particle collision operator, or when the particle motion is not strictly Hamiltonian in a collisionless simulation, the marker distribution will evolve along the particle trajectory and, in general, cannot be known exactly. A two-dimensional numerical marker distribution is proposed as an approximation. It is shown to be advantageous over other common methods for evaluating the marker distribution in long-time turbulence simulations. A generalized two-weight δf-method is proposed to mitigate the marker evolution problem.

, Monte Carlo methods↗
DOE OSTI2022

Helium droplet infrared spectroscopy of the butyl radicals

Butyl radicals (n-, s-, i-, and tert-butyl) are formed from the pyrolysis of stable precursors (1-pentyl nitrite, 2-methyl-1-butyl nitrite, isopentyl nitrite, and azo-tert-butane, respectively). The radicals are doped into a beam of liquid helium droplets and probed with infrared action spectroscopy from 2700 to 3125 cm -1 , allowing for a low temperature measurement of the CH stretching region. The presence of anharmonic resonance polyads in the 2800–3000 cm -1 region complicates its interpretation. To facilitate spectral assignment, the anharmonic resonances are modeled with two model Hamiltonian approaches that explicitly couple CH stretch fundamentals to HCH bend overtones and combinations: a VPT2+K normal mode model based on coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] quartic force fields and a semi-empirical local mode model. Both of these computational methods provide generally good agreement with the experimental spectra.

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

Quantum computing for fusion energy science applications

This is a review of recent research exploring and extending present-day quantum computing capabilities for fusion energy science applications. We begin with a brief tutorial on both ideal and open quantum dynamics, universal quantum computation, and quantum algorithms. Then, we explore the topic of using quantum computers to simulate both linear and nonlinear dynamics in greater detail. Because quantum computers can only efficiently perform linear operations on the quantum state, it is challenging to perform nonlinear operations that are generically required to describe the nonlinear differential equations of interest. In this work, we extend previous results on embedding nonlinear systems within linear systems by explicitly deriving the connection between the Koopman evolution operator, the Perron–Frobenius evolution operator, and the Koopman–von Neumann evolution (KvN) operator. We also explicitly derive the connection between the Koopman and Carleman approaches to embedding. Extension of the KvN framework to the complex-analytic setting relevant to Carleman embedding, and the proof that different choices of complex analytic reproducing kernel Hilbert spaces depend on the choice of Hilbert space metric are covered in the appendixes. Finally, we conclude with a review of recent quantum hardware implementations of algorithms on present-day quantum hardware platforms that may one day be accelerated through Hamiltonian simulation. We discuss the simulation of toy models of wave–particle interactions through the simulation of quantum maps and of wave–wave interactions important in nonlinear plasma dynamics.

Joseph, I. (ORCID:0000000255400840)↗
DOE OSTI2023

Accelerating the density-functional tight-binding method using graphical processing units

Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Herein two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1–6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1–2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters.

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

TeraChem protocol buffers ( TCPB ): Accelerating QM and QM/MM simulations with a client–server model

The routine use of electronic structures in many chemical simulation applications calls for efficient and easy ways to access electronic structure programs. Here, we describe how the graphics processing unit (GPU) accelerated electronic structure program TeraChem can be set up as an electronic structure server, to be easily accessed by third-party client programs. We exploit Google’s protocol buffer framework for data serialization and communication. The client interface, called TeraChem protocol buffers (TCPB), has been designed for ease of use and compatibility with multiple programming languages, such as C++, Fortran, and Python. To demonstrate the ease of coupling third-party programs with electronic structures using TCPB, we have incorporated the TCPB client into Amber for quantum mechanics/molecular mechanics (QM/MM) simulations. The TCPB interface saves time with GPU initialization and I/O operations, achieving a speedup of more than 2× compared to a prior file-based implementation for a QM region with ~250 basis functions. We demonstrate the practical application of TCPB by computing the free energy profile of p-hydroxybenzylidene-2,3-dimethylimidazolinone (p-HBDI - )—a model chromophore in green fluorescent proteins—on the first excited singlet state using Hamiltonian replica exchange for enhanced sampling. All calculations in this work have been performed with the non-commercial freely-available version of TeraChem, which is sufficient for many QM region sizes in common use.

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

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree–Fock reference is also), and orbital invariant. Here, our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H 2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

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

Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data

In this work, we present a machine learning approach to calculating electronic specific heat capacities for a variety of benchmark molecular systems. Our models are based on data from density matrix quantum Monte Carlo, which is a stochastic method that can calculate the electronic energy at finite temperature. As these energies typically have noise, numerical derivatives of the energy can be challenging to find reliably. In order to circumvent this problem, we use Gaussian process regression to model the energy and use analytical derivatives to produce the specific heat capacity. From there, we also calculate the entropy by numerical integration. We compare our results to cubic splines and finite differences in a variety of molecules in which Hamiltonians can be diagonalized exactly with full configuration interaction. We finally apply this method to look at larger molecules where exact diagonalization is not possible and make comparisons with more approximate ways to calculate the specific heat capacity and entropy.

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

Spectral localizer for line-gapped non-Hermitian systems

Short-ranged and line-gapped non-Hermitian Hamiltonians have strong topological invariants given by an index of an associated Fredholm operator. It is shown how these invariants can be accessed via the signature of a suitable spectral localizer. Here, this numerical technique is implemented in an example with relevance to the design of topological photonic systems, such as topological lasers.

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

Synergy of semiempirical models and machine learning in computational chemistry

Catalyzed by enormous success in the industrial sector, many research programs have been exploring data-driven, machine learning approaches. Performance can be poor when the model is extrapolated to new regions of chemical space, e.g., new bonding types, new many-body interactions. Another important limitation is the spatial locality assumption in model architecture, and this limitation cannot be overcome with larger or more diverse datasets. The outlined challenges are primarily associated with the lack of electronic structure information in surrogate models such as interatomic potentials. Given the fast development of machine learning and computational chemistry methods, we expect some limitations of surrogate models to be addressed in the near future; nevertheless spatial locality assumption will likely remain a limiting factor for their transferability. Here, we suggest focusing on an equally important effort—design of physics-informed models that leverage the domain knowledge and employ machine learning only as a corrective tool. In the context of material science, we will focus on semi-empirical quantum mechanics, using machine learning to predict corrections to the reduced-order Hamiltonian model parameters. The resulting models are broadly applicable, retain the speed of semiempirical chemistry, and frequently achieve accuracy on par with much more expensive ab initio calculations. These early results indicate that future work, in which machine learning and quantum chemistry methods are developed jointly, may provide the best of all worlds for chemistry applications that demand both high accuracy and high numerical efficiency.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modified Fermi’s golden rule rate expressions

Fermi’s golden rule (FGR) serves as the basis for many expressions of spectroscopic observables and quantum transition rates. The utility of FGR has been demonstrated through decades of experimental confirmation. However, there still remain important cases where the evaluation of a FGR rate is ambiguous or ill-defined. Examples are cases where the rate has divergent terms due to the sparsity in the density of final states or time dependent fluctuations of system Hamiltonians. Strictly speaking, assumptions of FGR are no longer valid for such cases. However, it is still possible to define modified FGR rate expressions that are useful as effective rates. The resulting modified FGR rate expressions resolve a long standing ambiguity often encountered in using FGR and offer more reliable ways to model general rate processes. Finally, simple model calculations illustrate the utility and implications of new rate expressions.

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

The Grassmann interpolation method for spin-unrestricted open-shell systems

The recently reported Grassmann interpolation (G-Int) method [J. A. Tan and K. U. Lao, J. Chem. Phys. 158, 051101 (2023)] has been extended to spin-unrestricted open-shell systems. In contrast to closed-shell systems, where G-Int has to be performed only once since the α and β density matrices are the same, spin-unrestricted open-shell systems require G-Int to be performed twice—one for the α spin and another for the β spin density matrix. In this work, we tested the performance of G-Int to the carbon monoxide radical cation CO●+ and nickelocene complex, which have the doublet and triple ground states, respectively. We found that the Frobenius norm errors associated with the interpolations for the α and β spin density matrices are comparable for a given molecular geometry. These G-Int density matrices, when used as an initial guess for a self-consistent field (SCF) calculation, outperform the conventional SCF guess schemes, such as the superposition of atomic densities, purified superposition of atomic densities, core Hamiltonian, and generalized Wolfsberg–Helmholtz approximation. Depending on the desired accuracy, these G-Int density matrices can be used to directly evaluate the SCF energy without performing SCF iterations. In addition, the spin-unrestricted G-Int density matrices have been used for the first time to directly calculate the atomic charges using the Mulliken and ChElPG population analysis.

Chemistry↗
DOE OSTI2023

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Four-wave mixing with anti-parity-time symmetry in hot 85 Rb vapor

Here, we report an experimental demonstration of anti-parity-time symmetric optical four-wave mixing in thermal rubidium vapor, where the propagation of probe and stokes fields in a double-Λ scheme is governed by a non-Hermitian Hamiltonian. We are particularly interested in studying quantum intensity correlations between the two fields near the exceptional point, taking into account loss and accompanied Langevin noise. Our experimental measurements of classical four-wave mixing gain and the associated two-mode relative-intensity squeezing are in reasonable agreement with the theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Orbital dependent complications for close vs well-separated electrons in diradicals

We investigate two limits in open-shell diradical systems: O3, in which the interesting orbitals are in close proximity to one another, and (C21H13)2, where there is a significant spatial separation between the two orbitals. In accord with earlier calculations, we find that standard density-functional approximations do not predict the open-shell character for the former case but uniformly predict the open-shell character for the latter case. We trace the qualitatively incorrect behavior in O3 predicted by these standard density functional approximations to self-interaction error and use the Fermi–Löwdin-orbital-self-interaction-corrected formalism to determine accurate triplet, closed-shell singlet, and open-shell broken-spin-symmetry electronic configurations. Analysis of the resulting many-electron overlap matrices allows us to unambiguously show that the broken-spin-symmetry configurations do not participate in the representation of the Ms = 0 triplet states and allows us to reliably extract the singlet–triplet splitting in O3 by analyzing the energy as a function of Fermi-orbital-descriptor permutations. The results of these analyses predict the percentage of open-shell character in O3, which agrees well with conventional wavefunction-based methods. While these techniques are expected to be required in cases near the Coulson–Fischer point, we find that they will be less necessary in diradical systems with well-separated electrons, such as (C21H13)2. Results based on energies from self-interaction-corrected generalized gradient, local density, and Hartree–Fock approximations and experimental results are in generally good agreement for O3. These results help form the basis for deriving extended Heisenberg-like Hamiltonians that are needed for descriptions of molecular magnets when there are competing low-energy electronic configurations.

Chemistry↗