Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective 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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Perturbation theory for the logarithm of a positive operator

In various contexts in mathematical physics, such as out-of-equilibrium physics and the asymptotic information theory of many-body quantum systems, one needs to compute the logarithm of a positive unbounded operator. Examples include the von Neumann entropy of a density matrix and the flow of operators with the modular Hamiltonian in the Tomita-Takesaki theory. Often, one encounters the situation where the operator under consideration, which we denote by ∆, can be related by a perturbative series to another operator ∆ 0 , whose logarithm is known. We set up a perturbation theory for the logarithm log ∆. It turns out that the terms in the series possess a remarkable algebraic structure, which enables us to write them in the form of nested commutators plus some “contact terms”.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Quantum Neural Networks: Issues, Training, and Applications

Our work in the field aims at explaining the limitations and expressive power of Quantum Machine Learning models, as well as finding feasible training algorithms that could be implemented in near-term Quantum Computers. The promise of Quantum Machine Learning is that by incorporating quantum effects, such as entanglement, into machine learning models researchers can improve model performance and understand more complex datasets. This pledge is particularly pronounced in the design of Quantum neural networks (QNNs), a promising framework for creating quantum algorithms, that promise to outperform classical models by combining the speedups of quantum computation with the widespread successes of deep learning. We show that applying this approach alone to quantum deep learning is problematic given that an excess of entanglement between the hidden and visible layers can destroy the predictive power of our QNN models. We address the barren plateau problem by suggesting the use of a generative, unbounded, nonlinear loss function with simple gradients. The loss function quantifies how much the quantum states generated by the QNNs differ from the data and the goal during training is to minimize it. Finally, we showcase how to use generative training to construct a "classical-quantum" neural network to accurately interpolate between the ground states of a Molecular Hamiltonian, a central question in Quantum Chemistry.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

A computational procedure for the dynamics of flexible beams within multibody systems

The dynamic analysis of three dimensional elastic beams which experience large rotational and large deformational motions are examined. The beam motion is modeled using an inertial reference for the translational displacements and a body-fixed reference for the rotational quantities. Finite strain rod theories are then defined in conjunction with the beam kinematic description which accounts for the effects of stretching, bending, torsion, and transverse shear deformations. A convected coordinate representation of the Cauchy stress tensor and a conjugate strain definition is introduced to model the beam deformation. To treat the beam dynamics, a two-stage modification of the central difference algorithm is presented to integrate the translational coordinates and the angular velocity vector. The angular orientation is then obtained from the application of an implicit integration algorithm to the Euler parameter/angular velocity kinematical relation. The combined developments of the objective internal force computation with the dynamic solution procedures result in the computational preservation of total energy for undamped systems. The present methodology is also extended to model the dynamics of deployment/retrieval of the flexible members. A moving spatial grid corresponding to the configuration of a deployed rigid beam is employed as a reference for the dynamic variables. A transient integration scheme which accurately accounts for the deforming spatial grid is derived from a space-time finite element discretization of a Hamiltonian variational statement. The computational results of this general deforming finite element beam formulation are compared to reported results for a planar inverse-spaghetti problem.

Downer, Janice Diane↗
DOE OSTI2023

Partial fillings of the bosonic $E$ 8 quantum Hall state

Here, we study bosonic topological phases constructed from electrons. In addition to a bulk excitation energy gap, these bosonic phases also have a fermion energy gap, below which all local excitations in the bulk and on the edge are even combinations of electrons. We focus on chiral phases, in which all low-energy edge excitations move in the same direction, that arise from the short-range entangled $E_8$ quantum Hall state, the bosonic analog of the filled lowest Landau level of electrons. The $E_8$ edge-state theory features an $E_8$ Kac-Moody symmetry that can be decomposed into ${\cal G}_A \times {\cal G}_B$ subalgebras, such as $SU(3) \times E_6$, $SO(M) \times SO(16-M)$, and $G_2 \times F_4$. (Here, $\{SO(M) \}$, $\{SU(N)\}$, and $\{E_8, G_2, F_4 \}$ denote orthogonal, unitary, and exceptional Lie algebras.) Using these symmetry decompositions, we construct exactly solvable coupled-wire model Hamiltonians for families of long-range entangled ${\cal G}_A$ or ${\cal G}_B$ bosonic fractional quantum Hall states that "partially fill" the $E_8$ state and are pairwise related by a generalized particle-hole symmetry. These long-range entangled states feature either Abelian or non-Abelian topological order. Some support the emergence of non-local Dirac and Majorana fermions, Ising anyons, metaplectic anyons, Fibonacci anyons, as well as deconfined $\mathbb{Z}_2$ gauge fluxes and charges.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

The Singlet–Triplet Gap of Cyclobutadiene: The CIPSI-Driven CC( P ; Q ) Study

An accurate determination of singlet−triplet gaps in biradicals, including cyclobutadiene in the automerization barrier region where one has to balance the substantial nondynamical many-electron correlation effects characterizing the singlet ground state with the predominantly dynamical correlations of the lowest-energy triplet, remains a challenge for many quantum chemistry methods. High-level coupled-cluster (CC) approaches, such as the CC method with a full treatment of singly, doubly, and triply excited clusters (CCSDT), are often capable of providing reliable results, but routine application of such methods is hindered by their high computational costs. We have recently proposed a practical alternative to converging the CCSDT energetics at small fractions of the computational effort, even when electron correlations become stronger and connected triply excited clusters are larger and nonperturbative, by merging the CC(P;Q) moment expansions with the selected configuration interaction methodology abbreviated as CIPSI. We demonstrate that one can accurately approximate the highly accurate CCSDT potential surfaces characterizing the lowest singlet and triplet states of cyclobutadiene along the automerization coordinate and the gap between them using tiny fractions of triply excited cluster amplitudes identified with the help of relatively inexpensive CIPSI Hamiltonian diagonalizations.

Basis sets↗
DOE OSTI2021

Magnetic excitations of the hybrid multiferroic (ND 4 ) 2 FeCl 5 ·D 2 O

We report a comprehensive inelastic neutron scattering study of the hybrid molecule-based multiferroic compound ( ND 4 ) 2 FeCl 5 · D 2 O in the zero-field incommensurate cycloidal phase and the high-field quasicollinear phase. The spontaneous electric polarization changes its direction concurrently with the field-induced magnetic transition, from mostly aligned with the crystallographic a axis to the c axis. To account for such a change in polarization direction, the underlying multiferroic mechanism was proposed to switch from the spin-current model induced via the inverse Dzyaloshinskii-Moriya interaction to the p - d hybridization model. We perform a detailed analysis of the inelastic neutron data of ( ND 4 ) 2 FeCl 5 · D 2 O using linear spin-wave theory to quantify magnetic interaction strengths and investigate the possible impact of different multiferroic mechanisms on the magnetic couplings. Overall our result reveals that the spin dynamics of both multiferroic phases can be well described by a Heisenberg Hamiltonian with easy-plane anisotropy. We do not find notable differences between the optimal model parameters of the two phases. The hierarchy of exchange couplings and the balance among frustrated interactions remain the same between two phases, suggesting that magnetic interactions in ( ND 4 ) 2 FeCl 5 · D 2 O are much more robust than the electric polarization in response to delicate reorganizations of the electronic degrees of freedom in an applied magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Study on the magnetothermal properties of the D y 1 - x T b x A l 2 series of compounds

Here, in this work, we developed a theoretical model Hamiltonian, in the mean field approximation, to describe the magnetic and magnetocaloric behavior of the series of compounds Dy 1-x Tb x Al 2 (x = 0.00, 0.15, 0.25, 0.30, 0.35, 0.40 and 0.75). We adjusted the exchange parameters λ DyDy , λ TbTb , and λ DyTb to obtain the spin reorientation temperatures (T SR ) and the critical temperature (T C ) for each compound in the series. The results obtained by the Hamiltonian model agree satisfactorily with the experimental results. The heat capacity curves with and without an applied magnetic field, adiabatic temperature variation and isothermal entropy variation were modeled and compared with the experimental data. As the experimental results show, our model was also able to reproduce the change in the spin reorientation process: a first order spin reorientation transition appears for concentrations x = 0.15, 0.25, 0.30, 0.35, and no spin reorientation transitions after x = 0.40.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Magnetic properties of a staggered S = 1 chain with an alternating single-ion anisotropy direction

Materials composed of spin-1 antiferromagnetic (AFM) chains are known to adopt complex ground states that are sensitive to the single-ion-anisotropy (SIA) energy ( D ), and intrachain ( J 0 ) and interchain ( J 1 , 2 ′ ) exchange energy scales. While theoretical and experimental studies have extended this model to include various other energy scales, the effect of the lack of a common SIA axis is not well explored. Here we investigate the magnetic properties of Ni ( pyrimidine ) ( H 2 O ) 2 ( NO 3 ) 2 , a chain compound where the tilting of Ni octahedra leads to a twofold alternation of the easy-axis directions along the chain. Muon-spin relaxation measurements indicate a transition to long-range order at T N = 2.3 K and the magnetic structure is initially determined to be antiferromagnetic and collinear using elastic neutron diffraction experiments. Inelastic neutron scattering measurements were used to find J 0 = 5.107 ( 7 ) K , D = 2.79 ( 1 ) K , J 1 ′ = 0.00 ( 5 ) K , J 2 ′ = 0.18 ( 3 ) K , and a rhombic anisotropy energy E = 0.19 ( 9 ) K . Mean-field modeling reveals that the ground state structure hosts spin canting of ϕ ≈ 6 . 5 ∘ , which is not detectable above the noise floor of the elastic neutron diffraction data. Monte Carlo simulation of the powder-averaged magnetization, M ( H ) , is then used to confirm these Hamiltonian parameters, while single-crystal M ( H ) simulations provide insight into features observed in the data. Published by the American Physical Society 2025

Vaidya, S. (ORCID:000000033975497X)↗
DOE OSTI2025

Classical Benchmarks for Variational Quantum Eigensolver Simulations of the Hubbard Model

Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to 32 qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

High field magnetization of FePS 3

High field magnetization measurements in pulsed fields up to 65 T have been performed on FePS 3 , which is nominally a good example of a two-dimensional Ising-like antiferromagnet on a honeycomb lattice. Measurements with the field parallel to the moment direction confirm the presence of two first-order transitions above 35 T, to M/M sat = 1/2 and M/M sat = 1, respectively, at 4 K. The measurements are in contradiction with published estimates for the magnetic exchange parameters, but the contradiction can be resolved by allowing for anisotropic exchange parameters in the Hamiltonian. Furthermore, the magnetization with the field perpendicular to the moment direction is anisotropic, with no transitions observed for fields along the a axis while a cascade of first-order transitions is observed for fields above 50 T along the b axis, the latter case also showing a strong degradation of the sample after repeated pulses. The results indicate a strong magnetolattice coupling in FePS 3 . Temperature-dependent measurements hint at a possible tricritical point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Characterization and thermometry of dissipatively stabilized steady states

In this work we study the properties of dissipatively stabilized steady states of noisy quantum algorithms, exploring the extent to which they can be well approximated as thermal distributions, and proposing methods to extract the effective temperature T. We study an algorithm called the relaxational quantum eigensolver (RQE), which is one of a family of algorithms that attempt to find ground states and balance error in noisy quantum devices. In RQE, we weakly couple a second register of auxiliary ‘shadow’ qubits to the primary system in Trotterized evolution, thus engineering an approximate zero-temperature bath by periodically resetting the auxiliary qubits during the algorithm’s runtime. Balancing the infinite temperature bath of random gate error, RQE returns states with an average energy equal to a constant fraction of the ground state. We probe the steady states of this algorithm for a range of base error rates, using several methods for estimating both T and deviations from thermal behavior. In particular, we both confirm that the steady states of these systems are often well-approximated by thermal distributions, and show that the same resources used for cooling can be adopted for thermometry, yielding a fairly reliable measure of the temperature. These methods could be readily implemented in near-term quantum hardware, and for stabilizing and probing Hamiltonians where simulating approximate thermal states is hard for classical computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Spin pumping from antiferromagnetic insulator spin-orbit-proximitized by adjacent heavy metal: a first-principles Floquet-nonequilibrium Green function study

Motivated by the recent experiment on spin pumping from sub-THz radiation-driven uniaxial antiferromagnetic insulator (AFI) MnF 2 into heavy metal (HM) Pt hosting strong spin-orbit (SO) coupling, we compute and compare pumped spin currents in Cu/MnF 2 /Cu and Pt/MnF 2 /Cu heterostructures. Recent theories of spin pumping by AFI have relied on simplistic Hamiltonians (such as tight-binding) and the scattering approach to quantum transport yielding the so-called interfacial spin mixing conductance (SMC), but the concept of SMC ceases to be applicable when SO coupling is present directly at the interface. In contrast, we use a more general first-principles quantum transport approach which combines noncollinear density functional theory with Floquet-nonequilibrium Green's functions in order to take into account: SO-proximitized AFI as a new type of quantum material, different from isolated AFI and brought about by AFI hybridization with adjacent HM layer; strong SO coupling at the interface; and evanescent wavefunctions penetrating from Pt or Cu into AFI layer to make its interfacial region conducting rather than insulating as in the isolated AFI. The DC component of pumped spin current $I_\mathrm{DC}^{S_z}$ vs. precession cone angle $\theta_{\boldsymbol{l}}$ of the Néel vector l of AFI does not follow putative $I^{S_z}_\mathrm{DC} \propto \sin^2 \theta_{\boldsymbol{l}}$, except for very small angles $\theta_{\boldsymbol{l}} \lesssim 10^\circ$ for which we define an effective SMC from the prefactor and find that it doubles from MnF2/Cu to MnF2/Pt interface. In addition, the angular dependence $I^{S_z}_\mathrm{DC}(\theta_{\boldsymbol{l}})$ differs for opposite directions of precession of the Néel vector, leading to twice as large SMC for the right-handed than for the left-handed chirality of the precession mode.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Variational Self-Consistent Theory for Beam-Loaded Cavities

A new variational theory is presented for beam loading in microwave cavities. The beam{field interaction is formulated as a dynamical interaction whose stationarity according to Hamilton’s principle will naturally lead to steady-state solutions that indicate how a cavity’s resonant frequency, Q and optimal coupling coefficient will detune as a result of the beam loading. A driven cavity Lagrangian is derived from first principles, including the effects of cavity wall losses, input power and beam interaction. The general formulation is applied to a typical klystron input cavity to predict the appropriate detuning parameters required to maximize the gain (or modulation depth) in the average Lorentz factor boost, (Δ γ). Numerical examples are presented, showing agreement with the general detuning trends previously observed in the literature. The developed formulation carries several advantages for beam-loaded cavity structures. It provides a self-consistent model for the dynamical (nonlinear) beam{field interaction, a procedure for maximizing gain under beam-loading conditions, and a useful set of parameters to guide cavity-shape optimization during the design of beam-loaded systems. Enhanced clarity of the physical picture underlying the problem seems to be gained using this approach, allowing straightforward inclusion or exclusion of different field configurations in the calculation and expressing the final results in terms of measurable quantities. Two field configurations are discussed for the klystron input cavity, using finite magnetic confinement or no confinement at all. Formulating the problem in a language that is directly accessible to the powerful techniques found in Hamiltonian dynamics and canonical transformations may potentially carry an additional advantage in terms of analytical computational gains, under suitable conditions

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Massless and massive photons within basis light-front quantization

We compute the structure function, transverse momentum dependent parton distributions (TMDs), and generalized parton distributions (GPDs) for the photon from the light-front quantum electrodynamics (QED) Hamiltonian, determined for its constituent bare photon and electron-positron Fock components within the basis light-front quantization framework. While the structure function and GPDs align well with expectations, some oscillatory behavior in the TMDs emerges as a finite-basis artifact. We also compute the GPD and the corresponding structure function for the longitudinal polarization mode of the massive photon. We further study the photon mass dependence of the calculated observables. After performing nonperturbative renormalization, we obtain a good quality description of the perturbative QED properties of the photon, as well as good agreement with the experimental data for the massless photon’s structure function.

Ab initio calculations↗
DOE OSTI2020

Molecular structure determination: Equilibrium structure of pyrimidine (m-C 4 H 4 N 2 ) from rotational spectroscopy (reSE) and high-level ab initio calculation (re) agree within the uncertainty of experimental measurement

The pure rotational spectrum of pyrimidine (m-C 4 H 4 N 2 ), the meta-substituted dinitrogen analog of benzene, has been studied in the millimeter-wave region from 235 GHz to 360 GHz. The rotational spectrum of the ground vibrational state has been assigned and fit to yield accurate rotational and distortion constants. Over 1700 distinct transitions were identified for the normal isotopologue in its ground vibrational state and least-squares fit to a partial sextic S-reduced Hamiltonian. Transitions for all four singly substituted 13 C and 15 N isotopologues were observed at natural abundance and were likewise fit. Deuterium-enriched samples of pyrimidine were synthesized, giving access to all eleven possible deuterium-substituted isotopologues, ten of which were previously unreported. Experimental values of rotational constants and computed values of vibration–rotation interaction constants and electron-mass corrections were used to determine semi-experimental equilibrium structures (r e SE ) of pyrimidine. The r e SE structure obtained using coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] corrections shows exceptional agreement with the re structure computed at the CCSD(T)/cc-pCV5Z level (≤0.0002 Å in bond distance and ≤0.03° in bond angle). Of the various computational methods examined, CCSD(T)/cc-pCV5Z is the only method for which the computed value of each geometric parameter lies within the statistical experimental uncertainty (2σ) of the corresponding semi-experimental coordinate. The exceptionally high accuracy and precision of the structure determination is a consequence of the large number of isotopologues measured, the precision and extent of the experimental frequency measurements, and the sophisticated theoretical treatment of the effects of vibration–rotation coupling and electron mass. Taken together, these demanding experimental and computational studies establish the capabilities of modern structural analysis for a prototypical monocyclic aromatic compound.

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

Conservative Spin-Magnitude Change in Orbital Evolution in General Relativity

We show that physical scattering observables for compact spinning objects in general relativity can depend on additional degrees of freedom in the spin tensor beyond those described by the spin vector alone. The impulse, spin kick, and leading-order waveforms exhibit such a nontrivial dependence. A signal of this additional structure is the change in the magnitude of the spin vector under conservative Hamiltonian evolution, similar to our previous studies in electrodynamics. These additional degrees of freedom describe dynamical mass multipoles of compact objects and decouple for black holes. We also show that the conservative impulse, spin kick, and change of the additional degrees of freedom are encoded in the eikonal phase.

Classical black holes↗
DOE OSTI2022

On the S-matrix of Ising field theory in two dimensions

We explore the analytic structure of the non-perturbative S-matrix in arguably the simplest family of massive non-integrable quantum field theories: the Ising field theory (IFT) in two dimensions, which may be viewed as the Ising CFT deformed by its two relevant operators, or equivalently, the scaling limit of the Ising model in a magnetic field. Our strategy is that of collider physics: we employ Hamiltonian truncation method (TFFSA) to extract the scattering phase of the lightest particles in the elastic regime, and combine it with S-matrix bootstrap methods based on unitarity and analyticity assumptions to determine the analytic continuation of the 2 → 2 S-matrix element to the complex s-plane. Focusing primarily on the “high temperature” regime in which the IFT interpolates between that of a weakly coupled massive fermion and the E 8 affine Toda theory, we will numerically determine 3-particle amplitudes, follow the evolution of poles and certain resonances of the S-matrix, and exclude the possibility of unknown wide resonances up to reasonably high energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗