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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Engineered thermalization and cooling of quantum many-body systems

We develop a scheme for engineering genuine thermal states in analog quantum simulation platforms by coupling local degrees of freedom to driven, dissipative ancilla pseudospins. We demonstrate the scheme in a many-body quantum spin lattice simulation setting. A Born-Markov master equation describing the dynamics of the many-body system is developed, and we show that if the ancilla energies are periodically modulated, with a carefully chosen hierarchy of timescales, one can effectively thermalize the many-body system. Through analysis of the time-dependent dynamical generator, we determine the conditions under which the true thermal state is an approximate dynamical fixed point for general system Hamiltonians. Finally, we evaluate the thermalization protocol through numerical simulation and discuss prospects for implementation on current quantum simulation hardware.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploring nonlinear Rashba effect and spin Hall conductivity in Janus MXenes W 2 ⁢CO ⁢𝑋 (𝑋=S, Se, Te)

Rashba spin-orbit coupling (RSOC) facilitates spin manipulation without relying on an external magnetic field, opening up exciting possibilities for advanced spintronic devices. In this paper, we examine the effects of crystal momentum (𝑘) nonlinearity and anisotropy on the conventional Rashba effect, with a particular focus on their impact on the spin Hall conductivity (SHC) in a newly predicted family of 2D Janus materials, W 2 ⁢CO⁢𝑋 (𝑋 =S, Se, Te). Using first-principles density functional theory calculations, we confirm the dynamical and mechanical stability of the studied 2D materials. Strikingly, this materials family exhibits pronounced nonlinear Rashba spin splitting at the Γ point of Brillouin zone near the Fermi level, which cannot be adequately described by the linear-𝑘 Rashba model. Therefore, third-order momentum contributions (𝑘 3 ) must be incorporated into the Rashba Hamiltonian. Our analysis reveals that among the studied systems, W 2 ⁢COS exhibits the highest 𝑘 3 contribution of −45.9 eV Å 3 , despite having the lowest linear Rashba constant. Here, a detailed analysis of electronic structure reveals topologically nontrivial behavior in these 2D materials, yielding sizable SHC that is primarily governed by the nonlinear Rashba effect. Notably, these materials also exhibit large spin Hall angle (0.018–2.5 at E 𝐹 ), which is comparable to that of in bulk topological insulators like Bi 2 ⁢Se 3 and Bi 2 ⁢Te 3 , and surpassing those in narrow bandgap bulk semiconductors GeTe and SnTe, as well as heavy metals such as Pt. Sizable SHC, large spin Hall angles, and the ability to tune SHC via electric fields without altering the topological properties, rooted in the crystal field splitting, underscore the potential of these materials for spintronic applications.

Electronic structure↗
DOE OSTI2020

Effective spin-mixing conductance of topological-insulator/ferromagnet and heavy-metal/ferromagnet spin-orbit-coupled interfaces: A first-principles Floquet-nonequilibrium Green function approach

The spin-mixing conductance (SMC) is a key quantity determining efficiency of spin transport across interfaces. Thus, knowledge of its precise value is required for accurate measurement of parameters quantifying numerous effects in spintronics, such as spin-orbit torque, spin Hall magnetoresistance, spin Hall effect, and spin pumping. However, the standard expression for SMC, provided by the scattering theory in terms of the reflection probability amplitudes, is inapplicable when strong spin-orbit coupling (SOC) is present directly at the interface. This is the precisely the case of topological-insulator/ferromagnet and heavy-metal/ferromagnet interfaces of great contemporary interest. We introduce an approach where first-principles Hamiltonian of these interfaces, obtained from noncollinear density functional theory (ncDFT) calculations, is combined with charge-conserving Floquet-nonequilibrium-Green-function formalism to compute directly the pumped spin current $I^{S_z}_L$ into semi-infinite left lead of two-terminal heterostructures Cu/X/Co/Cu or Y/Co/Cu—where X = Bi 2 Se 3 and Y = Pt or W—due to microwave-driven steadily precessing magnetization of the Co layer. This allows us to extract an effective SMC as a prefactor in $I^{S_z}_L$ versus precession cone angle θ dependence, as long as it remains the same, $I^{S_z}$ ∝ sin 2 θ , as in the case where SOC is absent. By comparing calculations where SOC in switched off versus switched on in ncDFT calculations, we find that SOC consistently reduces the pumped spin current and therefore the effective SMC.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tidal effects for spinning particles

Expanding on the recent derivation of tidal actions for scalar particles, we present here the action for a tidally deformed spin-1 / 2 particle. Focusing on operators containing two powers of the Weyl tensor, we combine the Hilbert series with an on-shell amplitude basis to construct the tidal action. With the tidal action in hand, we compute the leading-post-Minkowskian tidal contributions to the spin-1/2–spin-1/2 amplitude, arising at $\mathcal{O}$( G 2 ). Our amplitudes provide evidence that the observed long range spin-universality for the scattering of two point particles extends to the scattering of tidally deformed objects. From the scattering amplitude we find the conservative two-body Hamiltonian, linear and angular impulses, eikonal phase, spin kick, and aligned-spin scattering angle. We present analogous results in the electromagnetic case along the way.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Relativistic Douglas–Kroll–Hess calculations of hyperfine interactions within first-principles multireference methods

A relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas–Kroll–Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods, including spin–orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field formalism with restricted active space state interaction for spin–orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativistic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs), such as Tb(III)Pc2 and Tb(II)(CpiPr5)2 without ligand truncation using well-converged basis sets. In particular, for the divalent SMM, which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).

Chemistry↗
DOE OSTI2022

Yang-Mills observables: from $\mathrm{KMOC}$ to eikonal through $\mathrm{EFT}$

We obtain a conservative Hamiltonian describing the interactions of two charged bodies in Yang-Mills through $\mathcal{O}$ (α 2 ) and to all orders in velocity. Our calculation extends a recently-introduced framework based on scattering amplitudes and effective field theory (EFT) to consider color-charged objects. These results are checked against the direct integration of the observables in the Kosower-Maybee-O’Connell (KMOC) formalism. At the order we consider we find that the linear and color impulses in a scattering event can be concisely described in terms of the eikonal phase, thus extending the domain of applicability of a formula originally proposed in the context of spinning particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Quadratic-in-spin Hamiltonian at $ \mathcal{O} $(G 2 ) from scattering amplitudes

We obtain the quadratic-in-spin terms of the conservative Hamiltonian describing the interactions of a binary of spinning bodies in General Relativity through O ( G 2 ) and to all orders in velocity. Our calculation extends a recently-introduced framework based on scattering amplitudes and effective field theory to consider non-minimal coupling of the spinning objects to gravity. At the order that we consider, we establish the validity of the formula proposed in [1] that relates the impulse and spin kick in a scattering event to the eikonal phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Inviscid damping of an elliptical vortex subject to an external strain flow

Inviscid spatial Landau damping is studied experimentally for the case of oscillatory motion of a two-dimensional vortex about its elliptical equilibrium in the presence of an applied strain flow. Here, the experiments are performed using electron plasmas in a Penning–Malmberg trap. They exploit the isomorphism between the two-dimensional Euler equations for an ideal fluid and the drift-Poisson equations for the plasma, where plasma density is the analog of vorticity. Perturbed elliptical vortex states are created using E x B strain flows, which are generated by applying voltages to electrodes surrounding the plasma. Measurements of spatial Landau damping (also called critical-layer damping) are in agreement with previous studies in the absence of an applied strain, where the damping is due to a resonance between the local fluid motion and the vortex oscillations. Interestingly, the damping rate does not change significantly over a wide range of applied strain rates. This can be accurately predicted from the initial vorticity profile, even though the resonant frequency is reduced substantially due to the applied strain. For higher amplitude perturbations, nonlinear trapping oscillations also exhibit behavior similar to the strain-free case. In principle, higher-order effects of the applied strain, such as separatrix crossing of peripheral vorticity and interactions with harmonics of the fundamental resonance, are expected to change the damping rate. However, this occurs only for conditions that are not realized in the experiments described here. Vortex-in-cell simulations are used to investigate the possible roles of these effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Stochastic quantum Krylov protocol with double-factorized Hamiltonians

Here we propose a class of randomized quantum Krylov diagonalization (rQKD) algorithms capable of solving the eigenstate estimation problem with modest quantum resource requirements. Compared to previous real-time evolution quantum Krylov subspace methods, our approach expresses the time evolution operator e –i$\widehat{H}$$\tau$ as a linear combination of unitaries and subsequently uses a stochastic sampling procedure to reduce circuit depth requirements. While our methodology applies to any Hamiltonian with fast-forwardable subcomponents, we focus on its application to the explicitly double-factorized electronic-structure Hamiltonian. To demonstrate the potential of the proposed rQKD algorithm on near-term quantum devices, we provide numerical benchmarks for a variety of molecular systems with circuit-based state-vector simulators including the effects of sampling noise, achieving ground-state energy errors of less than 1 kcal mol -1 with circuit depths orders of magnitude shallower than those required for low-rank deterministic Trotter-Suzuki decompositions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Hamiltonian learning using machine-learning models trained with continuous measurements

Here, we build upon recent work on the use of machine-learning models to estimate Hamiltonian parameters using continuous weak measurement of qubits as input. We consider two settings for the training of our model: (1) supervised learning, where the weak-measurement training record can be labeled with known Hamiltonian parameters, and (2) unsupervised learning, where no labels are available. The first has the advantage of not requiring an explicit representation of the quantum state, thus potentially scaling very favorably to a larger number of qubits. The second requires the implementation of a physical model to map the Hamiltonian parameters to a measurement record, which we implement using an integrator of the physical model with a recurrent neural network to provide a model-free correction at every time step to account for small effects not captured by the physical model. We test our construction on a system of two qubits and demonstrate accurate prediction of multiple physical parameters in both the supervised context and the unsupervised context. We demonstrate that the model benefits from larger training sets, establishing that it is “learning,” and we show robustness regarding errors in the assumed physical model by achieving accurate parameter estimation in the presence of unanticipated single-particle relaxation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

The role of three-dimensional shapes in the break-up of charged drops

The nonlinear dynamics of nonaxisymmetric inviscid charged conducting drops near the Rayleigh charge limit (R = 4) is investigated analytically. It is shown that only axisymmetric spheroid drops bifurcate from the sphere family when the charge is increased, that oblate spheroids at R greater than 4 are unstable to nonaxisymmetric disturbances governing drop breakup, and that prolate spheroids at R less than 4 are unstable only to axisymmetric disturbances tending to increase the length of the drop along its symmetry axis. The effects of external electric fields and rigid-body rotation are also analyzed, and the solutions for the amplitude equations at R just less than 4 (equivalent to the dynamical equations of the Henon-Heiles Hamiltonian) are explored.

Natarajan, R.↗
NASA NTRS1994

Modal decomposition of Hamiltonian variational equations

Over any finite arc of trajectory, the variational equations of a Hamiltonian system can be separated into 'normal' modes. This transformation is canonical, and the Lyapunov exponents over the trajectory arc occur as positive/negative pairs for conjugate modes, while the modal vectors remain unit vectors. This decomposition effectively solves the variational equations for any canonical, linear-dependent system. As an example, we study the Voyager I trajectory. In an interplanetary flyby, some of the modal variables increase by very large multiplicative factors, but this means that their conjugate modal variables decrease by those same very large multiplicative vectors. Maneuver strategies for this case are explored, and the minimum delta upsilon maneuver is found.

Wiesel, William E.↗
NASA NTRS2003

Spurious States-Free Solution of the k dot p Hamiltonian for Heterostructures

A method for eliminating spurious solution in the k dot p Hamiltonian has been proposed. Introduction of additional off-diagonal alpha k(exp 2) term converts spurious solution with large real wave vector to evanescent solution with large imaginary wave vector. This modification keeps the same effective masses at Gamma point and introduces negligible deviation from original nonparabolicity. A set of unphysical fast oscillation eigenfunctions in confined states of heterostructures are removed.

Kolokolov, Konstantin I.↗
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

Higher-order topological superconductivity in twisted bilayer graphene

In this work, we show that introducing spin-singlet or spin-triplet superconductivity into twisted bilayer graphene induces higher-order topological superconductivity. Multiple copies of C 2z ⁢T-protected Majorana Kramers pairs appear at corners on pairing domain walls. The topology originates from the anomaly analyzed in Song et al.—the absence of a lattice support—of the single-valley band structure of twisted bilayer graphene, which is protected by C 2z ⁢T and approximate particle-hole symmetry $\mathcal{P}$. We prove that any pairing (spin-singlet or spin-triplet) term preserving valley-U(1), spin-SU(2), time-reversal, C 2z ⁢T, and $\mathcal{P}$ must drive the system into a higher-order topological superconductor. Here spin-SU(2) is the global spin-SU(2) for the singlet pairing and is broken to U(1) for the triplet pairing. Using a Dirac Hamiltonian, we derive the corner modes and confirm with numerics. These corner states are stable even if $\mathcal{P}$ is weakly broken, which is true in experimental setups. Finally, we suggest experimental detection via the fractional Josephson effect in a TBG-TSC Josephson junction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1973

Rotation of solid bodies in the solar system

The effects of elastic distortion, nonprincipal axis rotation, precessing orbits, and internal dissipation on the rotation of a solid solar system body, which is in the gravitational field of an exterior body, are relatively easily analyzed by a Hamiltonian theory developed here. Examples of applications include the Chandler wobble, wobble of the moon, spin-orbit coupling, generalized Cassini laws, and tidal evolution.

Peale, S. J.↗
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

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

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗