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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Strong influence of nonmagnetic ligands on the momentum-dependent spin splitting in antiferromagnets

Here, recent studies show the nonrelativistic antiferromagnetic ordering could generate momentum-dependent spin splitting analogous to the Rashba effect but free from the requirement of relativistic spin-orbit coupling. Whereas the classification of such compounds can be illustrated by different spin-splitting prototypes (SSTs) from symmetry analysis and density-functional-theory calculations, the huge variation in chemical bonding and structures of these diverse compounds possibly clouds the issue of how much of the variation in spin splitting can be traced back to the symmetry-defined characteristics, rather to the underlining chemical and structural diversity. The alternative model Hamiltonian approaches do not confront the issues of chemical and structural complexity but often consider only the magnetic sublattice, dealing with the all-important effects of the nonmagnetic ligands via renormalizing the interactions between the magnetic sites. To this end, we constructed a DFT model Hamiltonian that allows us to study SSTs at constant chemistry while retaining the realistic atomic-scale structure including ligands. This is accomplished by using a single, universal magnetic skeletal lattice (Ni 2+ ions in rocksalt NiO) and designing small displacements of the nonmagnetic (oxygen) sublattice which produce, by design, the different SST magnetic symmetries. We show that (i) even similar crystal structures having very similar band structures can lead to contrasting behavior of spin splitting vs momentum, and (ii) even subtle deformations of the nonmagnetic ligand sublattice could cause a giant spin splitting in AFM-induced SST. This is a paradigm shift relative to the convention of modeling magnets without considering the nonmagnetic ligand that mediates indirect magnetic interaction (e.g., superexchange).

36 MATERIALS SCIENCE↗
DOE OSTI2026

SCF Framework, HF Stability, and RPA Correlation for Jordan–Wigner-Transformed Spin Hamiltonians on Arbitrary Coupling Topologies

Mapping spins to fermions via the Jordan–Wigner (JW) transformation can render mean-field (Hartree–Fock, HF) descriptions effective for strongly correlated spin systems. As established in recent work, the application of such approaches is not limited by the nonlocal structure of JW strings or by site ordering because string operators can be absorbed into Thouless rotations of a Slater determinant, and the variational optimization of a unitary Lie-algebraic similarity transformation removes any ordering dependence. Leveraging these ideas, we develop a self-consistent field (SCF) scheme that expresses the mean-field energy as a functional of the single-particle density matrix, providing an alternative to gradient-based optimization of Thouless parameters. We derive the analytical orbital Hessian to diagnose HF stability and compute the ground-state correlation energy through the random-phase approximation (RPA). Benchmark results for the XXZ and J 1 –J 2 model on one- and two-dimensional lattices demonstrate that RPA significantly improves mean-field accuracy.

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

Sub Coulomb barrier d+ 208 Pb scattering in the time-dependent basis function approach

We employ the non-perturbative time-dependent basis function (tBF) approach to study the scattering of the deuteron on 208 Pb below the Coulomb barrier. We obtain the bound and discretized scattering states of the projectile, which form the basis representation of the tBF approach, by diagonalizing a realistic Hamiltonian in a large harmonic oscillator basis. We find that the higher-order inelastic scattering effects are noticeable for sub barrier scatterings with the tBF method. We have successfully reproduced experimental sub Coulomb barrier elastic cross section ratios with the tBF approach by considering only the electric dipole (E1) component of the Coulomb interaction between the projectile and the target during scatterings. Here, we find that the correction of the polarization potential to the Rutherford trajectory is dominant in reproducing the data at very low bombarding energies, whereas the role of internal transitions of the deuteron projectile induced by the E1 interaction during the scattering becomes increasingly significant at higher bombarding energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Deformed Fredkin model for the ν = 5 / 2 Moore-Read state on thin cylinders

We propose a frustration-free model for the Moore-Read quantum Hall state on sufficiently thin cylinders with circumferences ≲ 7 magnetic lengths. While the Moore-Read Hamiltonian involves complicated long-range interactions between triplets of electrons in a Landau level, our effective model is a simpler one-dimensional chain of qubits with deformed Fredkin gates. We show that the ground state of the Fredkin model has high overlap with the Moore-Read wave function and accurately reproduces the latter's entanglement properties. Moreover, we demonstrate that the model captures the dynamical response of the Moore-Read state to a geometric quench, induced by suddenly changing the anisotropy of the system. We elucidate the underlying mechanism of the quench dynamics and show that it coincides with the linearized bimetric field theory. The minimal model introduced here can be directly implemented as a first step towards quantum simulation of the Moore-Read state, as we demonstrate by deriving an efficient circuit approximation to the ground state and implementing it on an IBM quantum processor. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Control optimization for parametric Hamiltonians by pulse reconstruction

Optimal control techniques provide a means to tailor the control pulses required to generate customized quantum gates, which helps to improve the resilience of quantum simulations to gate errors and device noise. However, the significant amount of (classical) computation required to generate customized gates can quickly undermine the effectiveness of this approach, especially when pulse optimization needs to be iterated. We propose a method to reduce the computational time required to generate the control pulse for a Hamiltonian that is parametrically dependent on a time-varying quantity. We use simple interpolation schemes to accurately reconstruct the control pulses from a set of pulses obtained in advance for a discrete set of predetermined parameter values. We obtain a reconstruction with very high fidelity and a significant reduction in computational effort. We report the results of the application of the proposed method to device-level quantum simulations of the unitary (real) time evolution of two interacting neutrons based on superconducting qubits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2021

Quantum Monte Carlo studies of a trimer scaling function with microscopic two- and three-body interactions

Here, we present an energy scaling function to predict, in a specific range, the energy of bosonic trimers with large scattering lengths and finite range interactions, which is validated by quantum Monte Carlo calculations using microscopic Hamiltonians with two- and three-body potentials. The proposed scaling function depends on the scattering length, effective range, and a reference energy, which we chose as the trimer energy at unitarity. We obtained the scaling function as a limit cycle from the solution of the renormalized zero-range model with effective range corrections. We proposed a simple parametrization of the energy scaling function. Besides the intrinsic interest in theoretical and experimental investigations, this scaling function allows one to probe Efimov physics with only the trimer ground states, which may open opportunities to identify Efimov trimers whenever access to excited states is limited.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Temperature dependence of magnetic anisotropy and magnetoelasticity from classical spin-lattice calculations

Here we present a classical molecular-spin dynamics (MSD) methodology that enables accurate computations of the temperature dependence of the magnetocrystalline anisotropy as well as magnetoelastic properties of magnetic materials. The nonmagnetic interactions are accounted for by a spectral neighbor analysis potential (SNAP) machine-learned interatomic potential, whereas the magnetoelastic contributions are accounted for using a combination of an extended Heisenberg Hamiltonian and a Néel pair interaction model, representing both the exchange interaction and spin-orbit-coupling effects, respectively. All magnetoelastic potential components are parameterized using a combination of first-principles and experimental data. Our framework is applied to the α phase of iron. Initial testing of our MSD model is done using a 0 K parametrization of the Néel interaction model. After this, we examine how individual Néel parameters impact the $B$ 1 and $B$ 2 magnetostrictive coefficients using a moment-independent δ sensitivity analysis. The results from this study are then used to initialize a genetic algorithm optimization which explores the Néel parameter phase space and tries to minimize the error in the B 1 and B 2 magnetostrictive coefficients in the range of 0–1200 K. Our results show that while both the 0 K and genetic algorithm optimized parametrization provide good experimental agreement for $B$ 1 and $B$ 2 , only the genetic algorithm optimized results can capture the second peak in the $B$ 1 magnetostrictive coefficient which occurs near approximately 800 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Computing Edge States without Hard Truncation

We present a numerical method which accurately computes the discrete spectrum and associated bound states of semi-infinite Hamiltonians which model electronic “edge” states localized at boundaries of one- and two-dimensional crystalline materials. The problem is nontrivial since arbitrarily large finite “hard” (Dirichlet) truncations of the Hamiltonian in the infinite bulk direction tend to produce spurious bound states partially supported at the truncation. Our method, which overcomes this difficulty, is to compute the Green's function of the semi-infinite Hamiltonian by imposing an appropriate boundary condition in the bulk direction; then, the spectral data is recovered via Riesz projection. We demonstrate our method's effectiveness by studies of edge states at a graphene zig-zag edge in the presence of defects modeled both by a discrete tight-binding model and a continuum PDE model under finite difference discretization. Our method may also be used to study states localized at domain wall-type edges in one- and two-dimensional materials where the edge Hamiltonian is infinite in both directions; we demonstrate this for the case of a tight-binding model of distinct honeycomb structures joined along a zig-zag edge. Here, we expect our method to be useful for designing novel devices based on precise wave-guiding by edge states.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Coupled Cluster Study of the Heats of Formation of UF 6 and the Uranium Oxyhalides, UO 2 X 2 (X = F, Cl, Br, I, and At)

The atomization enthalpies of the U(VI) species UF6 and the uranium oxyhalides UO 2 X 2 (X=F, Cl, Br, I, At) were calculated using a composite relativistic Feller-Peterson-Dixon (FPD) approach based on scalar relativistic DKH3-CCSD(T) with extrapolations to the CBS limit. The inherent multideterminantal nature of the U atom was mitigated by utilizing the singly charged atomic cation in all calculations with correction back to the neutral asymptote via the accurate ionization energy of the U atom. The effects of SO coupling were recovered using full 4-component CCSD(T) with contributions due to the Gaunt Hamiltonian calculated using Dirac-Hartree-Fock. The final atomization enthalpy for UF 6 (752.2 kcal/mol) was within 2.5 kcal/mol of the experimental value, but unfortunately the latter carries a ±2.4 kcal/mol uncertainty that is predominantly due to the experimental uncertainty in the formation enthalpy of U atom. The analogous value for UO 2 F 2 (607.6 kcal/mol) was in nearly exact agreement with experiment, but the latter has a stated experimental uncertainty of ±4.3 kcal/mol. The FPD atomization enthalpy for UO 2 Cl 2 (540.4 kcal/mol) was within the experimental error limits of ±5.5 kcal/mol. FPD atomization energies for the non-U-containing molecules (used for reaction enthalpies) H 2 O and HX (X=F, Cl, Br, I, At) were within at most 0.3 kcal/mol of their experimental values where available. The FPD atomization enthalpies, together with FPD reaction enthalpies for two different reactions, were used to determine heats of formation for all species of this work with estimated uncertainties of ±4 kcal/mol. The calculated heat of formation for UF 6 (-511.0 kcal/mol) is within 2.5 kcal/mol of the accurately-known (±0.45 kcal/mol) experimental value.

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

Toward more accurate adiabatic connection approach for multireference wavefunctions

A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the “fixed-RDM” approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.

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

Optimizing qubit control pulses for state preparation

In the burgeoning field of quantum computing, the precise design and optimization of quantum pulses are essential for enhancing qubit operation fidelity. This study focuses on refining the pulse engineering techniques for superconducting qubits, employing a detailed analysis of square and Gaussian pulse envelopes under various approximation schemes. We evaluated the effects of coherent errors induced by naive pulse designs. Furthermore, we identified the sources of these errors in the Hamiltonian model’s approximation level. We mitigated these errors through adjustments to the external driving frequency and pulse durations, thus implementing a pulse scheme with stroboscopic error reduction. Our results demonstrate that these refined pulse strategies improve performance and reduce coherent errors. Moreover, the techniques developed herein are applicable across different quantum architectures, such as ion-trap, atomic, and photonic systems.

Chirp modulation↗
DOE OSTI2024

Hamiltonian switching control of noisy bipartite qubit systems

Abstract We develop a Hamiltonian switching ansatz for bipartite control that is inspired by the quantum approximate optimization algorithm, to mitigate environmental noise on qubits. We demonstrate the control for a central spin coupled to bath spins via isotropic Heisenberg interactions, and then make physical applications to the protection of quantum gates performed on superconducting transmon qubits coupling to environmental two-level-systems (TLSs) through dipole-dipole interactions, as well as on such qubits coupled to both TLSs and a Lindblad bath. The control field is classical and acts only on the system qubits. We use reinforcement learning with policy gradient to optimize the Hamiltonian switching control protocols, using a fidelity objective for specific target quantum gates. We use this approach to demonstrate effective suppression of both coherent and dissipative noise, with numerical studies achieving target gate implementations with fidelities over 0.9999 (four nines) in the majority of our test cases and showing improvement beyond this to values of 0.999 999 999 (nine nines) upon a subsequent optimization by GRadient Ascent Pulse Engineering (GRAPE). We analyze how the control depth, total evolution time, number of environmental TLS, and choice of optimization method affect the fidelity achieved by the optimal protocols and reveal some critical behaviors of bipartite control of quantum gates.

Physics↗
DOE OSTI2026

Toward scalable quantum computations of atomic nuclei

We solve the nuclear two-body and three-body bound states via quantum simulations of pionless effective field theory on a lattice in position space. While the employed lattice remains small, the usage of local Hamiltonians including two- and three-body forces ensures that the number of Pauli terms scales linearly with increasing numbers of lattice sites. We use an adaptive ansatz grown from unitary coupled cluster theory to parametrize the ground states of the deuteron and 3 He, compute their corresponding energies, and analyze the scaling of the required computational resources. Our quantum simulations reproduce exact benchmarks for 2 H and 3 He within 100 keV, requiring at most 30 layers in the ansatz and thus resulting in modest circuit depths. Additionally, we find the number of shots required to reach a given precision scales linearly in the lattice size and more mildly in the system size. Furthermore, based on the agreement with exact benchmarks and mild scaling, we conclude that this can be an efficient, scalable approach for quantum computations of nuclear ground states, particularly to prepare initial states for quantum phase estimation or other filtering algorithms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Management and Operation of the Lawrence Livermore National Laboratory (Final Report)

Optimize existing characterization and control methods by developing methods for rapid synthesis of multi-qubit control for error mitigation and complex gate design using optimal control theory. The Subcontractor shall develop Hamiltonian learning methods that can be used for characterizing super-conducting quantum devices, in particular develop strategies for effective probing of such devices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Nitrogen Vacancy Center in Diamond for the Stress and Field Sensing Applications

The nitrogen-vacancy (NV) center in a nanodiamond (ND) crystal is a promising material for quantum information processing, sensing, and computing applications. It is one of the best candidate materials for quantum sensing and metrology expected to work at elevated temperatures and pressures conditions. We computationally show the effect of strain on the defect band edges and band gaps in the NV center diamond. A low energy Hamiltonian is developed for the ±1 spin manifold at the ground state. We show the quantum sensing device is a few orders of magnitude superior in sensing than the traditional optical sensing devices. We also discuss experimental results from the optically detected magnetic resonance (ODMR) and the spin relaxometry for the field sensing applications. The presentation concludes by providing a model for free spins detection of the rare earth ions.

field sensing applications↗
DOE OSTI2023

Modeling of standing wave RF cavities for tracking through multi-pass energy recovery linac

Short bunches, high current and multiple linac pass are all characteristics of Energy Recovery Linacs (ERLs), which may result in collective effects. They in turn, may affect the beam, degrading its quality, or even yield to instabilities causing a beam loss. To study and mitigate these effects one needs a numerical simulation code, that can take into account both the collective effects, as well as, particular ERL features, such as a multi-turn design that does not reach a steady state or the multiple passages of the beam through Radio-Frequency (RF) cavities at different energies. CODAL [1], a code developed by SOLEIL in collaboration with IJCLab, enables such studies. It is a 6 dimensional (6D) tracking code applying 'kicks' based on the integration of the local Hamiltonian for each element of the lattice. It is also capable of simulating space charge, wakefields and coherent synchrotron radiation. However, to correctly take into account the ERL dynamics, an upgrade had to be made to include the effect of a standing wave RF cavity in 6D. In this paper, we will concentrate on the implementation and benchmarking (with DESY’s tracking code ASTRA [2]) of both the longitudinal and the transverse models (by J.B. Rosenzweig and L. Serafini [3]), which we use to carry out tracking of fully analytical 6D RF cavity.

Bogacz, A.↗