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

DeepONet-grid-UQ: A trustworthy deep operator framework for predicting the power grid’s post-fault trajectories

This paper proposes a novel data-driven method for the reliable prediction of the power grid’s post-fault trajectories, i.e., the power grid’s dynamic response after a disturbance or fault. Here, the proposed method is based on the recently proposed concept of Deep Operator Networks (DeepONets). Unlike traditional neural networks that learn to approximate functions, DeepONets are designed to approximate nonlinear operators, i.e., mappings between infinite-dimensional spaces. Under this operator framework, we design a novel and efficient DeepONet that (i) takes as inputs the trajectories collected before and during the fault and (ii) outputs the predicted post-fault trajectories. In addition, we endow our method with the much-needed ability to balance efficiency with reliable/trustworthy predictions via uncertainty quantification. To this end, we propose and compare two novel methods that enable quantifying the predictive uncertainty. First, we propose a Bayesian DeepONet (B-DeepONet) that uses stochastic gradient Hamiltonian Monte-Carlo to sample from the posterior distribution of the DeepONet trainable parameters. Then, we design a Probabilistic DeepONet (Prob-DeepONet) that uses a probabilistic training strategy to enable quantifying uncertainty at virtually no extra computational cost. Finally, we validate the proposed methods’ predictive power and uncertainty quantification capability using the New York-New England power grid model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Six-dimensional matching of intense beam with linear accelerating structure

Beam matching is a common technique that is routinely employed in accelerator design with the aim of minimizing beam losses and preservation of beam brightness. Despite being widely used, a full theoretical understanding of beam matching in 6D remains elusive. In this work, we present an analytical treatment of 6D beam matching of a high-intensity beam onto an RF structure. We begin our analysis within the framework of a linear model, and apply the averaging method to a set of 3D beam envelope equations. Accordingly, we obtain a matched solution that is comprised of smoothed envelopes and periodic terms, describing envelope oscillations with the period of the focusing structure. We then consider the nonlinear regime, where the beam size is comparable with the separatrix size. Stating with a Hamiltonian analysis in 6D phase space, we attain a self-consistent beam profile and show that it is significantly different from the commonly used ellipsoidal shape. Subsequently, we analyze the special case of an equilibrium with equal space charge depression between all degrees of freedom. Comparison of beam dynamics for equipartitioned, equal space charge depression, and equal emittances beams is given. Finally, we present experimental results on beam matching in the LANSCE linac.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

On the connection between perturbation theory and new semiclassical expansion in quantum mechanics

Here, it is shown that for the one-dimensional anharmonic oscillator with potential V(x) = ax 2 + bgx 3 + ... = $\frac{1}{g^2}$ $\hat{V}$ (gx), as well as for the radial oscillator V(r) = $\frac{1}{g^2}$ $\hat{V}$ (gr) and for the perturbed Coulomb problem V(r) = $\frac{a}{r}$ + βgr + ... = g $\tilde{V}$ (gr), the Perturbation Theory in powers of the coupling constant g (weak coupling regime) and the semiclassical expansion in powers of $\hbar$ 1/2 for the energies coincide. This is related to the fact that the dynamics developed in two spaces: x (r)-space and gx (gr)-space, lead to the same energy spectra. The equations which govern dynamics in these two spaces, the Riccati-Bloch equation and the Generalized Bloch equation, respectively, are presented. It is shown that the perturbation theory for the logarithmic derivative of the wavefunction in - space leads to (true) semiclassical expansion in powers of $\hbar$ 1/2 ; for the one-dimensional case this corresponds to the flucton calculus for the density matrix in the path integral formalism in Euclidean (imaginary) time proposed by one of the authors. Matching the perturbation theory in powers of g and the semiclassical expansion in powers of $\hbar$ 1/2 for the wavefunction leads to a highly accurate local approximation in the entire coordinate space, its expectation value for the Hamiltonian provides a prescription for the summation of the perturbative (trans)-series.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Spin crossover transition driven by pressure: Barocaloric applications

This article describes a mean-field theoretical model for Spin-Crossover (SCO) materials and explores its implications. It is based on a simple Hamiltonian that yields the high spin molar fraction as a function of temperature and pressure, as well as a temperature–pressure phase diagram for the SCO transition. In order to test the model, it was compared with the giant Barocaloric Effect (BCE) of the SCO material [FeL 2 ][B 4 ] 2 . Furthermore, we found that optical phonons are responsible for 92% of the total barocaloric entropy change. DFT calculations successively indicates that, as expected, the majority of this effect can be traced to low frequencies modes of vibration (400 cm -1 ), associated to the Fe coordination.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Structure formation in turbulence as an instability of effective quantum plasma

Structure formation in turbulence can be understood as an instability of “plasma” formed by fluctuations serving as effective particles. These “particles” are quantumlike in the sense that their wavelengths are non-negligible compared to the sizes of background coherent structures. The corresponding “kinetic equation” describes the Wigner matrix of the turbulent field, and the coherent structures serve as collective fields. This formalism is usually applied to manifestly quantumlike or scalar waves. Here, we show how to systematically extend it to more complex systems using compressible Navier–Stokes turbulence as an example. In this case, the fluctuation Hamiltonian is a five-dimensional matrix operator and diverse modulational modes are present. As an illustration, we calculate these modes for a sinusoidal shear flow and find two modulational instabilities. One of them is specific to supersonic flows, and the other one is a Kelvin–Helmholtz-type instability that is a generalization of the known zonostrophic instability. Our calculations are readily extendable to other types of turbulence, for example, magnetohydrodynamic turbulence in plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fibonacci fast convergence for neutrino oscillations in matter

Understanding neutrino oscillations in matter requires a non-trivial diagonalization of the Hamiltonian. As the exact solution is very complicated, many approximation schemes have been pursued. Here we show that one scheme, systematically applying rotations to change to a better basis, converges exponentially fast wherein the rate of convergence follows the Fibonacci sequence. We find that the convergence rate of this procedure depends very sensitively on the initial choices of the rotations as well as the mechanism of selecting the pivots. We then apply this scheme for neutrino oscillations in matter and discover that the optimal convergence rate is found using the following simple strategy: first apply the vacuum (2-3) rotation and then use the largest off-diagonal element as the pivot for each of the following rotations. The Fibonacci convergence rate presented here may be extendable to systems beyond neutrino oscillations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Basis light-front quantization approach to photon

We solve for the light-front wave functions (LFWFs) of the physical photon from the eigenvectors of the light-front quantum electrodynamics (QED) Hamiltonian with the aim to determine its bare photon and electron-positron Fock components. We then employ the resulting LFWFs to compute the transverse momentum dependent parton distributions (TMDs) and the generalized parton distributions (GPDs) at zero skewness of the photon. The TMDs are found to be in excellent agreement with the lowest-order perturbative results calculated using the electron-positron quantum fluctuation of the photon. The GPDs are also consistent with the perturbative calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Transverse structure of the pion beyond leading twist with basis light-front quantization

We investigate the twist-3transverse-momentum-dependent parton distribution functions (TMDs) of the pion with basis light-front quantization. The twist-3 TMDs are not independent and can be decomposed into twist-2and genuine twist-3 terms from the equations of motion (EOM). We compute the TMDs from the resulting light-front wave functions obtained by diagonalizing the light-front QCD Hamiltonian, determined for pion’s constituent quark-antiquark and quark-antiquark-gluon Fock sectors, with three-dimensional confinement. We also obtain the twist-3 parton distribution functions (PDFs) and show that they preserve the sum rule, which affirms the robustness of our approach. This is the first time that theoretical predictions are made for subleading twist structures of the pion containing interference terms between two light-front Fock sectors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗