Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HAMILTONIANS”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

ELSI — An open infrastructure for electronic structure solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn–Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. As a result, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics

Since variational symplectic integrators for the guiding center was proposed, structure-preserving geometric algorithms have become an active research field in plasma physics. We found that the slow manifolds of the classical Pauli particle enable a family of structure-preserving geometric algorithms for guiding center dynamics with long-term stability and accuracy. This discovery overcomes the difficulty associated with the unstable parasitic modes for variational symplectic integrators when applied to the degenerate guiding center Lagrangian. It is a pleasant surprise that Pauli's Hamiltonian for electrons, which predated the Dirac equation and marks the beginning of particle physics, reappears in classical physics as an effective algorithm for solving an important plasma physics problem. This technique is applicable to other degenerate Lagrangians reduced from regular Lagrangians.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Accelerated impurity solver for DMFT and its diagrammatic extensions

Here, we present ComCTQMC, a GPU accelerated quantum impurity solver. It uses the continuous-time quantum Monte Carlo (CTQMC) algorithm wherein the partition function is expanded in terms of the hybridisation function (CT-HYB). ComCTQMC supports both partition and worm-space measurements, and it uses improved estimators and the reduced density matrix to improve observable measurements whenever possible. ComCTQMC efficiently measures all one and two-particle Green's functions, all static observables which commute with the local Hamiltonian, and the occupation of each impurity orbital. ComCTQMC can solve complex-valued impurities with crystal fields that are hybridized to both fermionic and bosonic baths. Most importantly, ComCTQMC utilizes graphical processing units (GPUs), if available, to dramatically accelerate the CTQMC algorithm when the Hilbert space is sufficiently large. We demonstrate acceleration by a factor of over 600 (100) in a simulation of δ-Pu at 600 K with (without) crystal fields. In easier problems, the GPU offers less impressive acceleration or even decelerates the CTQMC. Here we describe the theory, algorithms, and structure used by ComCTQMC in order to achieve this set of features and level of acceleration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Calculating elements of matrix functions using divided differences

In this work, we introduce a method for calculating individual elements of matrix functions. Our technique makes use of a novel series expansion for the action of matrix functions on basis vectors that is memory efficient even for very large matrices. We showcase our approach by calculating the matrix elements of the exponential of a transverse-field Ising model and evaluating quantum transition amplitudes for large many-body Hamiltonians of sizes up to 2 64 x 2 64 on a single workstation. We also discuss the application of the method to matrix inverses. We relate and compare our method to the state-of-the-art and demonstrate its advantages. We also discuss practical applications of our method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

QuDPy: A Python-based tool for computing ultrafast non-linear optical responses

Nonlinear Optical Spectroscopy is a well-developed field with theoretical and experimental advances that have benefited multiple disciplines, including chemistry, biology, and physics. However, for the accurate interpretation of the corresponding multi-dimensional spectra, there is a need for precise quantum dynamical simulations based on model Hamiltonians. In this article, we present the initial release of our code, QuDPy (Quantum Dynamics in Python), which provides a robust numerical platform for performing quantum dynamics simulations based on model systems, including open quantum systems. Furthermore, a distinguishing feature of our approach is the ability to specify various high-order optical response pathways in the form of double-sided Feynman diagrams through a straightforward input syntax. This syntax outlines the time-ordering of ket-sided or bra-sided optical interactions acting on the time-evolving density matrix of the system. We utilize the quantum dynamics capabilities of QuTip to simulate the spectral response of complex systems, allowing us to compute virtually any $n$-th order optical response of the model system. To illustrate the utility of our approach, we provide a series of example calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

High-order finite element method for atomic structure calculations

We introduce featom, an open source code that implements a high-order finite element solver for the radial Schrödinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number and/or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. Here, to address the slow convergence of the $κ=±1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy (10 -8 Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schrödinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers (Z) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver dftatom. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples and tests, wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Advanced measurement techniques in quantum Monte Carlo: The permutation matrix representation approach

In a typical finite temperature quantum Monte Carlo (QMC) simulation, estimators for simple static observables such as specific heat and magnetization are known. With a great deal of system-specific manual labor, one can sometimes also derive more complicated non-local or even dynamic observable estimators. In contrast, we show that arbitrary static observables can be estimated within the permutation matrix representation (PMR) flavor for any Hamiltonian. We then generalize these results to general imaginary-time correlation functions and non-trivial integrated susceptibilities thereof. Finally, we demonstrate the practical versatility of our method by estimating various non-local, random observables for the transverse-field Ising model on a square lattice and a toy random model.

Permutation matrix representation↗
DOE OSTI2025

Quasi-optical beam tracing module development for millimeter-wave high-wavenumber collective scattering on the NSTX-U and EAST tokamaks

A Python3-based beam tracing code utilizing Quasi-Optics has been developed to track both incident and receiving beams in high-k collective millimeter wave scattering systems within magnetic fusion plasmas. In contrast to existing ray tracing codes that solely consider refraction, this beam tracing code incorporates diffraction phenomena, providing a more comprehensive calculation. Here, this enhanced capability allows for a more accurate calculation of the scattering volume and spatial resolution in high-k collective scattering systems, crucial for evaluating system performance and facilitating data analysis. Unlike Geometrical Optics, Quasi-Optics employs the complex eikonal method, representing a Gaussian beam as a collection of coupled rays to accurately preserve diffraction characteristics. The developed code is intended for application in NSTX-Upgrade and EAST high-k beam tracing analyses, targeting frequencies of 693 GHz and 270 GHz, respectively. The high-k system's primary objective is the observation of electron-scale instabilities. Employing a symplectic integrator, the code ensures numerical accuracy, assessed through the conservation of the Hamiltonian. With its precision and efficiency, the code facilitates rapid inter-shot analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Magnetothermal properties of Tm x Dy 1- x Al 2 ( x = 0.25, 0.50 and 0.75)

Here, we describe magnetic, thermal, and magnetocaloric properties of rare earth intermetallic compounds Tm x Dy 1-x Al 2 with 0.25, 0.5 and 0.75. Using model Hamiltonian we consider contributions of the crystalline electric field anisotropy in both Tm and Dy magnetic sublattices, disorder in exchange interactions among Tm-Tm, Dy-Dy and Tm-Dy magnetic ions, and the Zeeman effect. Employing earlier reported and new experimental measurements, we first determine a single free variable – the intersublattice magnetic exchange parameter – to properly model the temperature and magnetic field dependencies of heat capacity and magnetization, and then use the modeling results to explain the emergence of an anomalous spin reorientation transition and its influence on the magnetocaloric effect in the title compounds. Theoretical results agree with experimental data reasonably well.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Anisotropic magneto-resistivity and magnetocaloric effect in DyAl 2

Here, in this study, we experimentally and theoretically investigate the correlation between the anisotropic magnetic resistivity and magnetocaloric effect in DyAl 2 single crystals. DyAl 2 crystallizes in the C15 Laves phase structure with cubic symmetry. While some earlier studies revealed that the isothermal entropy change in dialuminides follows a similar trend as the electrical resistivity change as a function of temperature and magnetic field, the correlation between anisotropic behavior of these two physical parameters is yet to be explored to the best of our knowledge. To address this gap, we measured electrical resistivity of DyAl 2 single crystals along two different directions in zero and applied magnetic fields and compared the experimental results with our theoretical model. For the theoretical analysis, we employed a model Hamiltonian in the mean-field approximation, considering the exchange interaction, Zeeman effect, and crystalline electric field associated with the cubic symmetry. Our findings reveal a significant dependence of DyAl 2 resistivity on the direction of the applied magnetic field, in agreement with our theoretical results. These outcomes emphasize the interplay between magnetocaloric effect and magneto-resistivity, underscoring the potential of the magnetocaloric effect as a valuable tool for gaining deeper insights into basic physical properties including electron transport behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Learning functional priors and posteriors from data and physics

In this work, we develop a new Bayesian framework based on deep neural networks to be able to extrapolate in space-time using historical data and to quantify uncertainties arising from both noisy and gappy data in physical problems. Specifically, the proposed approach has two stages: (1) prior learning and (2) posterior estimation. At the first stage, we employ the physics-informed Generative Adversarial Networks (PI-GAN) to learn a functional prior either from a prescribed function distribution, e.g., Gaussian process, or from historical data and physics. At the second stage, we employ the Hamiltonian Monte Carlo (HMC) method to estimate the posterior in the latent space of PI-GANs. In addition, we use two different approaches to encode the physics: (1) automatic differentiation, used in the physicsinformed neural networks (PINNs) for scenarios with explicitly known partial differential equations (PDEs), and (2) operator regression using the deep operator network (DeepONet) for PDE-agnostic scenarios. We then test the proposed method for (1) meta-learning for one-dimensional regression, and forward/inverse PDE problems (combined with PINNs); (2) PDE-agnostic physical problems (combined with DeepONet), e.g., fractional diffusion as well as saturated stochastic (100-dimensional) flows in heterogeneous porous media; and (3) spatial-temporal regression problems, i.e., inference of a marine riser displacement field using experimental data from the Norwegian Deepwater Programme (NDP). The results demonstrate that the proposed approach can provide accurate predictions as well as uncertainty quantification given very limited scattered and noisy data, since historical data could be available to provide informative priors. In summary, the proposed method is capable of learning flexible functional priors, e.g., both Gaussian and non-Gaussian process, and can be readily extended to big data problems by enabling mini-batch training using stochastic HMC or normalizing flows since the latent space is generally characterized as low dimensional.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Structure preservation using discrete gradients in the Vlasov-Poisson-Landau system

We present a novel structure-preserving framework for solving the Vlasov-Poisson-Landau system of equations using a particle in cell (PIC) discretization combined with discrete gradient time integrators. The Vlasov-Poisson-Landau system is an accurate model for studying hot plasma dynamics at a kinetic scale where small-angle Coulomb collisions dominate. Our scheme guarantees conservation of mass, momentum and energy as well as preservation of the monotonicity of entropy production in both the time-continuous and discrete systems. We employ the conservative integrator for both the Hamiltonian Vlasov-Poisson equations and the dissipative Landau equation using the PETSc library (www.mcs.anl.gov/petsc) to showcase structure-preserving properties.

Discrete gradients↗
DOE OSTI2023

An atomistic theory of nucleation: Self-organization via non-equilibrium work and fluctuations

For this work, insights from non-equilibrium statistical mechanics, highlighting the role of work and fluctuations at the microscale, are applied toward the development of a fundamental, rigorous and purely atomistic theory of nucleation. Nanoscale fluctuations in order, density and heat influence the local nucleation rate by orders of magnitude, necessitating their inclusion through a modern approach. Coarse-graining over the underlying Hamiltonian dynamics allows derivation of a microscale expression for the nucleation rate in terms of a classical path integral over far from equilibrium trajectories and their associated work. Second law violating states at the microscale, as found from the dynamics of small critical nucleation clusters, contribute exponentially to the observable macroscale nucleation rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Magnetothermal properties of Ho 1-x Dy x Al 2 (x = 0, 0.05, 0.10, 0.15, 0.25 and 0.50) compounds

Magnetic and magnetocaloric properties of H o 1 - x D y x A l 2 compounds with x = 0 , 0.05 , 0.10 , 0.15 , 0.25 and 0.50 , modelled using a Hamiltonian that includes the exchange interactions between Ho-Dy, Ho-Ho and Dy-Dy ions in addition to the crystalline electric field and the Zeeman effects, have been compared with those determined experimentally. In order to reproduce experimentally observed global ferromagnetic ordering temperatures and spin reorientation transition temperatures as x Dy varies, the exchange interactions between Ho-Dy and Ho-Ho were set as free parameters and adjusted to match the experimental results. We demonstrate that heat capacity of polycrystalline materials in non-zero magnetic fields can be satisfactory reproduced by using the average of multiple magnetic field directions with respect to the crystallographic coordinate system, while reasonably good agreement between experimentally determined and theoretically predicted magnetocaloric effects can be achieved considering an average of only three field directions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Magnetic phase diagram of A 2 [FeCl 5 (H 2 O)] ( A = K, Rb, NH 4 )

Here, erythrosiderites with the formula A 2 FeX 5 ∙H 2 O, where A = Rb, K, and (NH 4 ) and X = Cl and Br are intriguing systems that possess various magnetic and electric phases, as well as multiferroic phases in which magnetism and ferroelectricity are coupled. In this report, we study the magnetic phase diagram of erythrosiderites as a function of superexchange interactions and magnetic anisotropies. To this end, we perform classical Monte Carlo simulations on magnetic Hamiltonians that contain five different superexchange interactions with single-ion anisotropies. Our phase diagram contains all magnetic ground states that have been experimentally observed in these materials. We argue that the ground states can be explained by varying the ratio of $\frac{J_4}{J_2}$. For $\frac{J_4}{J_2}$ > 0.95 a cycloidal spins structure is stabilized as observed in (NH 4 ) 2 FeCl 5 ∙H 2 O and otherwise a collinear spin structure is stabilized as observed in (K,Rb) 2 FeCl 5 ∙H 2 O. We also show that the difference in the single-ion anisotropy along a- and c- axes is essential to stabilize the intermediate state observed in (NH) 2 FeCl 5 ∙H 2 O.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ro-vibrational Spectroscopy of 2,5-Norbornadiene between 400 and 850 cm-1

Rotationally-resolved infrared absorption spectra of the two lowest-energy fundamental bands for each of the three infrared active symmetry species of 2,5-norbornadiene (C2v point group symmetry; asymmetry parameter ? = -0.2203), ?11 (a1), ?12(a1), ?28(b1), ?29(b1), ?38(b2), and ?39(b2), have been measured, assigned and fit to a semi-rigid rotor Hamiltonian to produce a set of spectroscopic constants for each band. The fit band centers are ?11 (a1), c-type transitions, 728.852357(3) cm-1; ?12 (a1), c-type transitions, 417.31604(1) cm-1; ?28 (b1), a-type transitions, 656.719741(3) cm-1; ?29 (b1), a-type transitions, 500.31425(1) cm-1; ?38 (b2), b-type transitions, 800.82243(2) cm-1; and ?39 (b2), b-type transitions, 539.23110(6) cm-1. Using ground state combination differences from the ?11 and ?28 bands and previously published microwave data, an improved set of ground state spectroscopic constants were fit for the parent isotopologue of the molecule. This research was conducted at the Pacific Northwest National Laboratory (PNNL), which is operated by Battelle Memorial Institute for the Department of Energy (DOE) under Contract No. DE-AC05-76RL01830.

Sams, R L.↗
DOE OSTI2023

High-resolution rovibrational spectroscopy of trans -formic acid in the v 1 OH stretching fundamental: Dark state coupling and evidence for charge delocalization dynamics

High-resolution infrared (IR) reduced-Doppler absorption spectra of jet-cooled gas phase trans-formic acid in the v 1 OH stretching fundamental region are reported for the first time, obtained by supersonically expanding trans- formic acid/Ar mixtures through a slit jet nozzle source and rotationally cooling to T rot ≈ 10.9(5) K, with ab- sorption signals recorded by high-resolution difference-frequency IR absorption spectroscopy. Two a/b-type rovibrational bands of comparable intensity, one ~10-fold weaker b-type band, and one ~6-fold weaker a-type band are observed, with vibrational band origins at 3570.493(5), 3566.793(5), 3560.032(9), and 3534.6869(2) cm –1 , respectively. Based on previous Raman jet spectroscopic work by Nejad and Sibert [A. Nejad, E.L. Sibert III, The Raman jet spectrum of trans-formic acid and its deuterated isotopologs: Combining theory and experi- ment to extend the vibrational database, J. Chem. Phys. 154(6) (2021) 064301.], these four rovibrational bands have been assigned to v 1 , (v 2 + v 7 ), (v 6 + 2v 7 + 2v 9 ), and 2v 3 , respectively. Specifically, two of the three upper dark states (2 1 7 1 (a') and 6 1 7 2 9 2 (a')) are close enough to the “bright” 1 1 (a') state to facilitate strong anharmonic resonance interactions, which results in intensity mixing into the two zero-order bands that would otherwise be “dark”. Furthermore, our high-resolution spectral analysis reveals that there are local rotational crossings be- tween these zero-order 1 1 and 2 1 7 1 states resulting in extra lines (i.e., some upper levels in the nominally v 1 band have majority zero-order 2 1 7 1 state character). This motivates development of a 3 coupled state (1 1 , 2 1 7 1 , and 6 1 7 2 9 2 ) picture to aid in the spectral analysis, which is able to match all 3 observed band origins and relative band intensities, as well as indicate the necessity of multistate (> 2) coupling. Though limited by range of J and Ka levels (J’ ≤ 9 and K a ’ ≤ 3) populated at supersonic jet temperatures, this work offers first precision spec- troscopic analysis of trans-formic acid in the v 1 OH stretch region, which should aid in assignment of the more complete yet highly congested room temperature FTIR spectra. Lastly, and in sharp contrast to the spectral complexity in the three predominantly b-type bands, the lone a-type 2v 3 rovibrational band at 3534.6869(2) cm –1 is well described by a simple, rigid asymmetric top Hamiltonian.

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

Additively manufactured refractory high-entropy alloys with superior radiation resistance

Refractory high-entropy alloys (RHEAs) are promising candidates for next-generation nuclear and high-temperature applications. Among many approaches to manufacture RHEAs, additive manufacturing (AM) represents the most recent and advanced metal manufacturing method which allows near-net-shape manufacturing to reduce material waste and post-processing time. However, performance of AM RHEAs under complex irradiation conditions remains largely unexplored. Here, in this study, we demonstrate for the first time the response of directed energy deposition (DED) AM quaternary RHEAs (HfTaVW, CrTaVW) subjected to sequential dual-beam ion irradiation, consisting of helium pre-implantation followed by high-dose heavy ion bombardment. Compositions of DED AM RHEAs were selected using Monte Carlo (MC) simulations based on a cluster expansion (CE) Hamiltonian parameterized by density functional theory (DFT). Post-irradiation microstructural characterization revealed that the AM RHEA maintained remarkable stability, with suppressed helium bubble growth and reduced defect accumulation compared to conventional alloys. Even at high doses (∼100 dpa), the alloy exhibited no void swelling, a low density of dislocation loops, and no evidence of severe degradation. These results highlight the intrinsic ability of AM-derived microstructures and multicomponent chemistry to synergistically mitigate irradiation effects. Our findings establish AM RHEAs as a class of materials with superior resistance to radiation damage under conditions relevant to advanced fusion and fission environments and demonstrate the importance of sequential ion beam studies in evaluating their long-term performance.

36 MATERIALS SCIENCE↗