Search NASA⌕ Search

SEARCH · Search NASA

Results for “Effective Hamiltonian models”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Robust preparation of ground state phases under noisy imaginary time evolution

Nonunitary state preparation protocols such as imaginary time evolution (ITE) offer substantial advantages relative to unitary ones, including the ability to prepare certain long-range correlated states more efficiently. Here, we ask whether such protocols are also robust to noise arising due to coupling to the environment. We consider a nonunitary ITE “circuit” subjected to a variety of noise models and investigate whether the resulting steady state remains in the same phase as the target state of the ITE at finite noise strength. Taking the one-dimensional quantum Ising model as a concrete example, we find that the ground-state order and associated phase transition persist in the presence of noise, provided the noise does not explicitly break the symmetry that protects the phase transition. That is, the noise must possess the protecting symmetry in a weak (or average) form. Our analysis is facilitated by a mapping to an effective Hamiltonian picture in a doubled Hilbert space. We discuss possible implications of these findings for quantum simulation on noisy quantum hardware. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Active control of Alfvén eigenmodes by external magnetic perturbations with different spatial spectra

Alfvén eigenmodes have been suppressed and excited in tokamak plasmas by (just) modifying the poloidal spectra of externally applied static magnetic perturbations. This effect is observed experimentally when toroidal spectra of n = 2, n = 4 as well as a mixed spectrum of n = 2 and n = 4 is applied. Under the n = 2 magnetic perturbations, the modes are excited or suppressed by modifying the coil phasing between the upper and the lower set of coils. Regardless of the absolute rotation, an even parity for the n = 4 perturbation is observed to reduce the amplitude of the Alfvénic instabilities, while an odd parity amplifies it. To combine the stabilizing (and destabilizing) effect of n = 2 and n = 4, a mixed spectrum is applied, finding similar reduction (and amplification) trends. However, the impact on the mode amplitude is more subtle, due to the reduced coil current required for a mixed spectrum. The signal level on the fast-ion loss detector is sensitive to the applied poloidal spectrum, which is consistent with Hamiltonian full-orbit modelling of an edge resonant transport layer activated by the 3D perturbative fields. An internal redistribution of the fast-ion population is induced, modifying the phase-space gradients driving the Alfvénic instabilities, and ultimately determining their existence. The calculated edge resonant layers for both n = 2 and n = 4 toroidal spectra are consistent with the observed suppressed and excited phases. Moreover, hybrid kinetic-magnetohydrodynamic (MHD) simulations reveal that this edge resonant transport layer overlaps in phase-space with the population responsible for the fast-ion drive. The results presented here may help to control fast-ion driven Alfvénic instabilities in future burning plasmas with a significant fusion born alpha particle population.

Alfvén waves↗
DOE OSTI2022

Second-order magnetic responses in quantum magnets: Magnetization under ac magnetic fields

Here we investigate second-order magnetic responses of quantum magnets against ac magnetic fields. We focus on the case where the z component of the spin is conserved in the unperturbed Hamiltonian and the driving field is applied in the xy plane. We find that linearly polarized driving fields induce a second-harmonic response, while circularly polarized fields generate only a zero-frequency response, leading to a magnetization with a direction determined by the helicity. Employing an unbiased numerical method, we demonstrate the nonlinear magnetic effect driven by the circularly polarized field in the XXZ model and show that the magnitude of the magnetization can be predicted by the dynamical spin structure factor in the linear response regime.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Abelian SU(N) 1 chiral spin liquids on the square lattice

In the physics of the Fractional Quantum Hall (FQH) effect, a zoo of Abelian topological phases can be obtained by varying the magnetic field. Aiming to reach the same phenomenology in spin-like systems, in this work we propose a family of SU(N)-symmetric models in the fundamental representation, on the square lattice with short-range interactions restricted to triangular units, a natural generalization for arbitrary N of an SU(3) model studied previously where time-reversal symmetry is broken explicitly. Guided by the recent discovery of SU(2) 1 and SU(3) 1 chiral spin liquids (CSL) on similar models we search for topological SU(N) 1 CSL in some range of the Hamiltonian parameters via a combination of complementary numerical methods such as exact diagonaliza- tions (ED), infinite density matrix renormalization group (iDMRG) and infinite Projected Entangled Pair State (iPEPS). Extensive ED on small (periodic and open) clusters up to N = 10 and an innovative SU(N)-symmetric version of iDMRG to compute entanglement spectra on (infinitely-long) cylinders in all topological sectors pro- vide unambiguous signatures of the SU(N) 1 character of the chiral liquids. An SU(4)-symmetric chiral PEPS, constructed in a manner similar to its SU(2) and SU(3) analogs, is shown to give a good variational ansatz of the N = 4 ground state, with chiral edge modes originating from the PEPS holographic bulk-edge correspondence. Finally, we discuss the possible observation of such Abelian CSL in ultracold atom setups where the possibility of varying N provides a tuning parameter similar to the magnetic field in the physics of the FQH effect.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Ab initio benchmarks of neutrinoless double- β decay in light nuclei with a chiral Hamiltonian

Here, we report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-β (0νββ) decays in light nuclei with mass numbers ranging from A = 6 to A = 22. We use the transition operator derived from light-Majorana neutrino exchange and evaluate the NME with three different methods: two variants of in-medium similarity renormalization group (IMSRG) and importance-truncated no-core shell model (IT-NCSM). The same two plus three nucleon interaction from chiral effective field theory is employed, and both isospin-conserving (ΔT = 0) and isospin-changing (ΔT = 2) transitions are studied. We compare our resulting ground-state energies and NMEs to those of recent ab initio no-core shell model and coupled-cluster calculations, also with the same inputs. We show that the NMEs of ΔT = 0 transitions are in good agreement among all calculations, at the level of 10%. For ΔT = 2, relative deviations are more significant in some nuclei. The comparison with the exact IT-NCSM result allows us to analyze these cases in detail, and indicates the next steps toward improving the IMSRG-based approaches. The present study clearly demonstrates the power of consistent cross checks that are made possible by ab initio methodology. This capability is crucial for providing meaningful many-body uncertainties in the NMEs for the 0νββ decays in heavier candidate nuclei, where quasiexact benchmarks are not available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Dynamical model of φ meson photoproduction on the nucleon and 4 He

We investigate φ meson photoproduction on the nucleon and the 4 He targets within a dynamical model approach based on a Hamiltonian which describes the production mechanisms by the Pomeron-exchange, meson-exchanges, φ radiations, and nucleon resonance excitations mechanisms. The final φN interactions are included being described by the gluon-exchange, direct φN couplings, and the box-diagrams arising from the couplings with πN, ρN, KΛ, and KΣ channels. The parameters of the Hamiltonian are determined by the experimental data of γp → φp from the CLAS Collaboration. The resulting Hamiltonian is then used to predict the coherent φ-meson production on the 4 He targets by using the distorted-wave impulse approximation. For the proton target, the final φN rescattering effects, as required by the unitarity condition, are found to be very weak, which supports the earlier calculations in the literature. For the 4 He targets, the predicted differential cross sections are in good agreement with the data obtained by the LEPS Collaboration. Furthermore, the role of each mechanism in this reaction is discussed and predictions for a wide range of scattering angles are presented, which can be tested in future experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Towards a nonperturbative construction of the S-matrix

We present a nonperturbative recipe for directly computing the S-matrix in strongly-coupled QFTs. The method makes use of spectral data obtained in a Hamiltonian framework and can be applied to a wide range of theories, including potentially QCD. We demonstrate the utility of this prescription in the specific example of the 2+1d O(N) model at large N, using energy eigenstates computed with Hamiltonian truncation to reproduce the full 2 → 2 scattering amplitude for arbitrary (complex) center-of-mass energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Enhancement of charge density wave correlations in a Holstein model with an anharmonic phonon potential

The Holstein Hamiltonian describes itinerant electrons whose site density couples to local phonon degrees of freedom. In the single-site limit, at half filling, the electron-phonon coupling results in a double-well structure for the lattice displacement, favoring empty or doubly occupied sites. In two dimensions and on a bipartite lattice in d ≥ 2, intersite hopping causes these doubly occupied and empty sites to alternate in a charge density wave (CDW) pattern when the temperature is lowered. Because a discrete symmetry is broken, this occurs in a conventional second-order transition at finite T cdw . In this paper, we investigate the effect of changing the phonon potential energy to one with an intrinsic double-well structure even in the absence of electron-phonon coupling. While this aids in the initial process of pair formation, the implications for subsequent CDW order are nontrivial. One expects that, when the electron-phonon coupling is too strong, the double wells become deep and the polaron mass gets large, an effect which reduces T cdw . Here, we show here the existence of regions of parameter space where the double-well potential, while aiding local pair formation, does so in a way which also substantially enhances long-range CDW order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Engineering an effective three-spin Hamiltonian in trapped-ion systems for applications in quantum simulation

Trapped-ion quantum simulators, in analog and digital modes, are considered a primary candidate to achieve quantum advantage in quantum simulation and quantum computation. The underlying controlled ion–laser interactions induce all-to-all two-spin interactions via the collective modes of motion through Cirac–Zoller or Mølmer–Sørensen schemes, leading to effective two-spin Hamiltonians, as well as two-qubit entangling gates. In this work, the Mølmer–Sørensen scheme is extended to induce three-spin interactions via tailored first- and second-order spin–motion couplings. The scheme enables engineering single-, two-, and three-spin interactions, and can be tuned via an enhanced protocol to simulate purely three-spin dynamics. Analytical results for the effective evolution are presented, along with detailed numerical simulations of the full dynamics to support the accuracy and feasibility of the proposed scheme for near-term applications. With a focus on quantum simulation, the advantage of a direct analog implementation of three-spin dynamics is demonstrated via the example of matter-gauge interactions in the U(1) lattice gauge theory within the quantum link model. The mapping of degrees of freedom and strategies for scaling the three-spin scheme to larger systems, are detailed, along with a discussion of the expected outcome of the simulation of the quantum link model given realistic fidelities in the upcoming experiments. The applications of the three-spin scheme go beyond the lattice gauge theory example studied here and include studies of static and dynamical phase diagrams of strongly interacting condensed-matter systems modeled by two- and three-spin Hamiltonians.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Truncation effects in the charge representation of the O(2) model

The O(2) model in Euclidean space-time is the zero-gauge-coupling limit of the compact scalar quantum electrodynamics. In this work, we obtain a dual representation of it called the charge representation. We study the quantum phase transition in the charge representation with a truncation to “spin $\textit{S}$," where the quantum numbers have an absolute value less than or equal to $\textit{S}$. The charge representation preserves the gapless-to-gapped phase transition even for the smallest spin truncation $\textit{S}$ = 1. The phase transition for $\textit{S}$ = 1 is an infinite-order Gaussian transition with the same critical exponents $\textit{δ}$ and $\textit{η}$ as the Berezinskii-Kosterlitz-Thouless (BKT) transition, while there are true BKT transitions for $\textit{S}$ ≥ 2 . The essential singularity in the correlation length for $\textit{S}$ = 1 is different from that for $\textit{S}$ ≥ 2 . The exponential convergence of the phase-transition point is studied in both Lagrangian and Hamiltonian formulations. We discuss the effects of replacing the truncated $\hat{U}^±$ = exp ($± i\hat{θ}$) operators by the spin ladder operators $\hat{S}^±$ in the Hamiltonian. The marginal operators vanish at the Gaussian transition point for $\textit{S}$ = 1, which allows us to extract the $\textit{η}$ exponent with high accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Gaussian processes meet NeuralODEs: a Bayesian framework for learning the dynamics of partially observed systems from scarce and noisy data

We present a machine learning framework (GP-NODE) for Bayesian model discovery from partial, noisy and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling and Gaussian Process priors over the observed system states. This allows us to exploit temporal correlations in the observed data, and efficiently infer posterior distributions over plausible models with quantified uncertainty. The use of the Finnish Horseshoe as a sparsity-promoting prior for free model parameters also enables the discovery of parsimonious representations for the latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed GP-NODE method including predator–prey systems, systems biology and a 50-dimensional human motion dynamical system. This article is part of the theme issue ‘Data-driven prediction in dynamical systems’.

Science & Technology - Other Topics↗
DOE OSTI2022

Quantum-Accurate Multiscale Modeling of Shock Hugoniots, Ramp Compression Paths, Structural and Magnetic Phase Transitions, and Transport Properties in Highly Compressed Metals

Fully characterizing high energy density (HED) phenomena using pulsed power facilities (Z machine) and coherent light sources is possible only with complementary numerical modeling for design, diagnostic development, and data interpretation. The exercise of creating numerical tests, that match experimental conditions, builds critical insight that is crucial for the development of a strong fundamental understanding of the physics behind HED phenomena and for the design of next generation pulsed power facilities. The persistence of electron correlation in HED materials arising from Coulomb interactions and the Pauli exclusion principle is one of the greatest challenges for accurate numerical modeling and has hitherto impeded our ability to model HED phenomena across multiple length and time scales at sufficient accuracy. An exemplar is a ferromagnetic material like iron, while familiar and widely used, we lack a simulation capability to characterize the interplay of structure and magnetic effects that govern material strength, kinetics of phase transitions and other transport properties. Herein we construct and demonstrate the Molecular-Spin Dynamics (MSD) simulation capability for iron from ambient to earth core conditions, all software advances are open source and presently available for broad usage. These methods are multi-scale in nature, direct comparisons between high fidelity density functional theory (DFT) and linear-scaling MSD simulations is done throughout this work, with advancements made to MSD allowing for electronic structure changes being reflected in classical dynamics. Main takeaways for the project include insight into the role of magnetic spins on mechanical properties and thermal conductivity, development of accurate interatomic potentials paired with spin Hamiltonians, and characterization of the high pressure melt boundary that is of critical importance to planetary modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Ab initio multishell valence-space Hamiltonians and the island of inversion

In the shell-model framework, valence-space Hamiltonians connecting multiple major-oscillator shells are of key interest for investigating the physics of neutron-rich nuclei, which have been the subject of intense experimental activity for decades. Here we present an extension of the ab initio valence-space in-medium similarity renormalization group, which allows the derivation of such Hamiltonians nonperturbatively. Starting from initial two- and three-nucleon forces from chiral effective field theory, we then calculate properties of nuclei in the important island-of-inversion region above oxygen, so far unexplored with ab initio methods. In this work, our results in the neon and magnesium isotopes indicate the importance of neutron excitation from the $\textit{sd}$ to $\textit{pf}$ shells and ground states dominated by intruder configurations around $\textit{N}$ = 20, consistent with the conclusions from phenomenological studies. We also benchmark the excitation spectrum of 16 O with coupled-cluster theory, finding generally good agreement, and discuss implications for ground-state energies and charge radii in oxygen and calcium isotopes. Finally we outline the proper procedure for treating the longstanding issue of center-of-mass contamination, and show that with a particular choice of valence space, these spurious states can be removed successfully.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Emergent interaction-induced topology in Bose-Hubbard ladders

We investigate the quantum many-body dynamics of bosonic atoms hopping in a two-leg ladder with strong on-site contact interactions. We observe that when the atoms are prepared in a staggered pattern with pairs of atoms on every other rung, singlon defects, i.e., rungs with only one atom, can localize due to an emergent topological model, even though the underlying model in the absence of interactions admits only topologically trivial states. This emergent topological localization results from the formation of a zero-energy edge mode in an effective lattice formed by two adjacent chains with alternating strong and weak hoping links (Su-Schrieffer-Heeger chains) and opposite staggering which interface at the defect position. Our findings open the opportunity to dynamically generate nontrivial topological behaviors without the need for complex Hamiltonian engineering. Published by the American Physical Society 2025

Wellnitz, David (ORCID:0000000349788938)↗
DOE OSTI2023

Stability of Kuramoto networks subject to large and small fluctuations from heterogeneous and spatially correlated noise

Oscillatory networks subjected to noise are broadly used to model physical and technological systems. Due to their nonlinear coupling, such networks typically have multiple stable and unstable states that a network might visit due to noise. In this article, we focus on the assessment of fluctuations resulting from heterogeneous and spatially correlated noise inputs on Kuramoto model networks. We evaluate the typical, small fluctuations near synchronized states and connect the network variance to the overlap between stable modes of synchronization and the input noise covariance. Going beyond small to large fluctuations, we introduce the indicator mode approximation that projects the dynamics onto a single amplitude dimension. Such an approximation allows for estimating rates of fluctuations to saddle instabilities, resulting in phase slips between connected oscillators. Statistics for both regimes are quantified in terms of effective noise amplitudes that are compared and contrasted for several noise models. In conclusion, bridging the gap between small and large fluctuations, we show that a larger network variance does not necessarily lead to higher rates of large fluctuations.

42 ENGINEERING↗
DOE OSTI2024

Gravitational form factors and mechanical properties of quarks in protons: A basis light-front quantization approach

We compute the gravitational form factors (GFFs) and study their applications for the description of the mechanical properties such as the pressure, shear force distributions, and the mechanical radius of the proton from its light-front wave functions (LFWFs) based on basis light-front quantization (BLFQ). The LFWFs of the proton are given by the lowest eigenvector of a light-front effective Hamiltonian that incorporates a three-dimensional confining potential and a one-gluon exchange interaction with fixed coupling between the constituent quarks solved in the valence Fock sector. We find acceptable agreement between our BLFQ computations and the lattice QCD for the GFFs. Our D -term form factor also agrees well with the extracted data from the deeply virtual Compton scattering experiments at Jefferson Lab, and the results of different phenomenological models. The distributions of pressures and shear forces are similar to those from different models. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Floquet Engineering of Interactions and Entanglement in Periodically Driven Rydberg Chains

Neutral atom arrays driven into Rydberg states constitute a promising approach for realizing programmable quantum systems. Enabled by strong interactions associated with Rydberg blockade, they allow for simulation of complex spin models and quantum dynamics. We introduce a new Floquet engineering technique for systems in the blockade regime that provides control over novel forms of interactions and entanglement dynamics in such systems. Our approach is based on time-dependent control of Rydberg laser detuning and leverages perturbations around periodic many-body trajectories as resources for operator spreading. These time-evolved operators are utilized as a basis for engineering interactions in the effective Hamiltonian describing the stroboscopic evolution. As an example, we show how our method can be used to engineer strong spin exchange, consistent with the blockade, in a one-dimensional chain, enabling the exploration of gapless Luttinger liquid phases. In addition, we demonstrate that combining gapless excitations with Rydberg blockade can lead to dynamic generation of large-scale multipartite entanglement. Experimental feasibility and possible generalizations are discussed.

Floquet systems↗