SEARCH · Search NASA

Results for “Hamiltonian structures”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Ab initio many-fermion structure calculations on a quantum computer

To overcome the limitations of existing algorithms for solving self-bound quantum many-body problems—such as those encountered in nuclear and particle physics—that access only a restricted subset of energy levels and provide limited structural information, we introduce and demonstrate a novel quantum-classical approach capable of resolving the complete bound-state spectrum. This method also provides the total angular momentum 𝐽 associated with each eigenstate. Here, our approach is based on expressing the Hamiltonian in second-quantized form within a novel input model combined with a scan scheme, enabling broad applicability to configuration-interaction calculations across diverse fields. We apply this hybrid method to compute, for the first time, the bound-state spectrum together with corresponding 𝐽 values of 20 O using a realistic strong-interaction Hamiltonian. Our approach applies to hadron spectra and 𝐽 values solved in the relativistic basis light-front quantization approach.

Du, Weijie [Chinese Academy of Sciences (CAS), Lan↗

DOE OSTI2026

Theory of ab initio downfolding with arbitrary-range electron-phonon coupling

Ab initio downfolding describes the electronic structure of materials within a low-energy subspace, often around the Fermi level. Typically starting from mean-field calculations, this framework allows for the calculation of one- and two-electron interactions, and the parametrization of a many-body Hamiltonian representing the active space of interest. The subsequent solution of such Hamiltonians can provide insights into the physics of strongly correlated materials. While phonons can substantially screen electron-electron interactions, electron-phonon coupling has been commonly ignored within ab initio downfolding, and when considered, this is done only for short-range coupling. Here we propose a theory of ab initio downfolding that accounts for short- and long-range electron-phonon coupling on equal footing. Our practical computational implementation is readily compatible with current downfolding approaches. We apply our approach to polar materials MgO and GeTe, and we reveal the importance of both short-range and long-range electron-phonon coupling in determining the magnitude of electron-electron interactions. Our results show that in the static limit, phonons reduce the on-site repulsion between electrons by 40% for MgO and by 79% for GeTe. Our framework also predicts that overall attractive nearest-neighbor interactions arise between electrons in GeTe, consistent with superconductivity in this material.

Tubman, Norm M↗

DOE OSTI2024

Emulator-Based Bayesian Calibration of the CISNET Colorectal Cancer Models

Purpose To calibrate Cancer Intervention and Surveillance Modeling Network (CISNET)'s SimCRC, MISCAN-Colon, and CRC-SPIN simulation models of the natural history colorectal cancer (CRC) with an emulator-based Bayesian algorithm and internally validate the model-predicted outcomes to calibration targets.Methods We used Latin hypercube sampling to sample up to 50,000 parameter sets for each CISNET-CRC model and generated the corresponding outputs. We trained multilayer perceptron artificial neural networks (ANNs) as emulators using the input and output samples for each CISNET-CRC model. We selected ANN structures with corresponding hyperparameters (i.e., number of hidden layers, nodes, activation functions, epochs, and optimizer) that minimize the predicted mean square error on the validation sample. We implemented the ANN emulators in a probabilistic programming language and calibrated the input parameters with Hamiltonian Monte Carlo-based algorithms to obtain the joint posterior distributions of the CISNET-CRC models' parameters. We internally validated each calibrated emulator by comparing the model-predicted posterior outputs against the calibration targets.Results The optimal ANN for SimCRC had 4 hidden layers and 360 hidden nodes, MISCAN-Colon had 4 hidden layers and 114 hidden nodes, and CRC-SPIN had 1 hidden layer and 140 hidden nodes. The total time for training and calibrating the emulators was 7.3, 4.0, and 0.66 h for SimCRC, MISCAN-Colon, and CRC-SPIN, respectively. The mean of the model-predicted outputs fell within the 95% confidence intervals of the calibration targets in 98 of 110 for SimCRC, 65 of 93 for MISCAN, and 31 of 41 targets for CRC-SPIN.Conclusions Using ANN emulators is a practical solution to reduce the computational burden and complexity for Bayesian calibration of individual-level simulation models used for policy analysis, such as the CISNET CRC models. In this work, we present a step-by-step guide to constructing emulators for calibrating 3 realistic CRC individual-level models using a Bayesian approach.

artificial neural networks↗

DOE OSTI2024

Kondo effect and its destruction in heterobilayer transition metal dichalcogenides

Moiré structures, along with line-graph-based d -electron systems, represent a setting to realize flat bands. One form of the associated strong correlation physics is the Kondo effect. Here, we address the recently observed Kondo-driven heavy fermion state and its destruction in AB-stacked hetero-bilayer transition metal dichalcogenides, which can be controlled by the gate voltages. By studying an effective interacting Hamiltonian using the slave spin approach, we obtained a phase diagram with the total filling factor and the displacement field strength as the tunable parameters. In an extended range of the tunable displacement field, our numerical results show that the relative filling of the d orbital, which is associated with the highest moiré band from the MoTe 2 layer, is enforced to be ν d ≈ 1 by the interaction. This agrees with the experimental observation. We also argue that the observed high coherence temperature scale could be explained by the non-negligible bandwidth of the d orbital. Our results set the stage to address the amplified quantum fluctuations that the Kondo effect may produce in these structures and new regimes that the systems open up for Kondo-destruction quantum criticality. Published by the American Physical Society 2024

Xie, Fang (ORCID:000000021327767X)↗

DOE OSTI2022

Constraints on the two-dimensional pseudospin-$\frac{1}{2}$Mott insulator description of $\mathrm{Sr_2IrO_4}$

Sr 2 IrO 4 has often been described via a simple, one-band pseudospin-$\frac{1}{2}$ model subject to electron-electron interactions on a square lattice, fostering analogies with cuprate superconductors believed to be well described by a similar model. In this work we argue - based on a detailed study of the low-energy electronic structure by circularly polarized spin and angle-resolved photoemission spectroscopy combined with dynamical mean-field theory calculations - that a pseudospin-$\frac{1}{2}$ model fails to capture the full complexity of the system. We show instead that a realistic multiband Hubbard Hamiltonian, accounting for the full correlated t 2g manifold, provides a detailed description of the interplay between spin-orbital entanglement and electron-electron interactions and yields quantitative agreement with experiments. Our analysis establishes that the j 3/2 states make up a substantial percentage of the low-energy spectral weight, i.e., approximately 74% as determined from the integration of the j-resolved spectral function in the 0 to -1.64eV energy range. Therefore, the results in our work are of relevance not only to Ir-based materials but also more generally to multiorbital materials with closely spaced energy scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Trends of Neutron Skins and Radii of Mirror Nuclei from First Principles

The neutron skin of atomic nuclei impacts the structure of neutron-rich nuclei, the equation of state of nucleonic matter, and the size of neutron stars. Here we predict the neutron skin of selected light- and medium-mass nuclei using coupled-cluster theory and the auxiliary field diffusion Monte Carlo method with two- and three-nucleon forces from chiral effective field theory. We find a linear correlation between the neutron skin and the isospin asymmetry in agreement with the liquid-drop model and compare with data. We also extract the linear relationship that describes the difference between neutron and proton radii of mirror nuclei and quantify the effect of charge symmetry breaking terms in the nuclear Hamiltonian. In conclusion, our results for the mirror-difference charge radii and binding energies per nucleon agree with existing data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the lambda-exchange of free energy perturbation (FEP). Further, two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/lambda-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1998

The Opacity of TiO from a Coupled Electronic State Calculation Parameterized by ab initio and Experimental Data

We have carried out ab initio electronic structure calculations of the spin-orbit and rotation-orbit couplings among the 14 lowest electronic states of TiO and used them to predict ro-vibrational energy levels. We report on the qualitative results as well as our progress in optimizing our Hamiltonian parameters in order to improve agreement with experimental line positions,

Schwenke, David W.↗

NASA NTRS1998

The Opacity of TiO from a Coupled Electronic State Calculation Parameterized by ab initio and Experimental Data

Schwenke, David W.↗

DOE OSTI2020

Gamow-shell-model description of Li isotopes and their mirror partners

Background: Weakly bound and unbound nuclei close to particle drip lines are laboratories of new nuclear structure physics at the extremes of neutron/proton excess. The comprehensive description of these systems requires an open quantum system framework that is capable of treating resonant and nonresonant many-body states on equal footing. Purpose: In this work, we develop the complex-energy configuration interaction approach to describe binding energies and spectra of selected 5 ≤ $\textit{A}$ ≤ 11 nuclei. Method: We employ the complex-energy Gamow shell model (GSM) assuming a rigid 4 He core. The effective Hamiltonian, consisting of a core-nucleon Woods-Saxon potential and a simplified version of the Furutani-Horiuchi-Tamagaki interaction with the mass-dependent scaling, is optimized in the $\textit{sp}$ space. To diagonalize the Hamiltonian matrix, we employ the Davidson method and the Density Matrix Renormalization Group technique. Results: We report our optimized GSM Hamiltonian offers a good reproduction of binding energies and spectra with the root-mean-square (rms) deviation from experiment of 160 keV. Since the model performs well when used to predict known excitations that have not been included in the fit, it can serve as a reliable tool to describe poorly known states. A case in point is our prediction for the pair of unbound mirror nuclei 10 Li – 10 N in which a huge Thomas-Ehrman shift dramatically alters the pattern of low-energy excitations. Conclusion: The new model will enable comprehensive studies of structure and reactions aspects of light drip-line nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Evidence for triaxial shape coexistence in 74 Ge

The deformation properties of the low-lying states in 74 Ge have been investigated using multistep Coulomb excitation. The measurements were carried out with the advanced 𝛾-ray tracking array, GRETINA, and the CHICO2 particle detector. A comprehensive set of 𝐸⁢2 transition and diagonal matrix elements was deduced following an analysis with the semiclassical coupled-channels code GOSIA. The data were compared with results of calculations carried out within the framework of the generalized triaxial rotor model as well as with the configuration interaction shell model and the symmetric rotor model. Results from calculations with covariant density functional theory were used to construct a five-dimensional collective Hamiltonian for further comparisons with the data. Collectively, the calculations provide an accurate reproduction of the experimental matrix elements and further support an understanding in terms of the coexistence of two axially asymmetric shapes. In conclusion, this leads to an overall interpretation of the underlying structure of 74 Ge requiring triaxiality, as is also the case in the neighboring even-mass Ge isotopes.

59 ≤ A ≤ 89↗

DOE OSTI2023

Magnetic ground state and perturbations of the distorted kagome Ising metal TmAgGe

Here, we present the magnetic orders and excitations of the distorted kagome intermetallic magnet TmAgGe. Using neutron single crystal diffraction we identify the propagation vectors k = ($\frac{1}{2}$ 0 0) and k = (0 0 0) and determine the magnetic structures of the zero-field and magnetic field-induced phases for H along the $\mathcal{a}$ and [–110] crystal directions. We determine the experimental magnetic field- temperature (H, T)-phase diagram and reproduce it by Monte Carlo simulations of an effective spin exchange Hamiltonian for one distorted kagome layer. Our model includes a strong axial single-ion anisotropy and significantly smaller exchange couplings, which span up to the third-nearest neighbors within the layer. Single crystal inelastic neutron scattering (INS) measurements reveal an almost flat, only weakly dispersive mode around 7 meV that we use alongside bulk magnetization data to deduce the crystal-electric field (CEF) scheme for the Tm 3+ ions. Random phase approximation (RPA) calculations based on the determined CEF wave functions of the two lowest quasidoublets enable an estimation of the interlayer coupling that is compatible with the experimental INS spectra. No evidence for low-energy spin waves associated to the magnetic order was found, which is consistent with the strongly Ising nature of the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics

Since variational symplectic integrators for the guiding center was proposed, structure-preserving geometric algorithms have become an active research field in plasma physics. We found that the slow manifolds of the classical Pauli particle enable a family of structure-preserving geometric algorithms for guiding center dynamics with long-term stability and accuracy. This discovery overcomes the difficulty associated with the unstable parasitic modes for variational symplectic integrators when applied to the degenerate guiding center Lagrangian. It is a pleasant surprise that Pauli's Hamiltonian for electrons, which predated the Dirac equation and marks the beginning of particle physics, reappears in classical physics as an effective algorithm for solving an important plasma physics problem. This technique is applicable to other degenerate Lagrangians reduced from regular Lagrangians.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Strong-field physics in three-dimensional topological insulators

We investigate theoretically the strong-field regime of light-matter interactions in the topological-insulator class of quantum materials. In particular, we focus on the process of nonperturbative high-order harmonic generation from the paradigmatic three-dimensional topological insulator bismuth selenide (Bi 2 Se 3 ) subjected to intense midinfrared laser fields. We analyze the contributions from the spin-orbit-coupled bulk states and the topological surface bands separately and reveal a major difference in how their harmonic yields depend on the ellipticity of the laser field. Furthermore, bulk harmonics show a monotonic decrease in their yield as the ellipticity increases, in a manner reminiscent of high harmonic generation in gaseous media. However, the surface contribution exhibits a highly nontrivial dependence, culminating with a maximum for circularly polarized fields. We attribute the observed anomalous behavior to (i) the enhanced amplitude and the circular pattern of the interband dipole and the Berry connections in the vicinity of the Dirac point and (ii) the influence of the higher-order, hexagonal warping terms in the Hamiltonian, which are responsible for the hexagonal deformation of the energy surface at higher momenta. The latter are associated directly with spin-orbit-coupling parameters. Our results thus establish the sensitivity of strong-field-driven high harmonic emission to the topology of the band structure as well as to the manifestations of spin-orbit interaction.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Optimal Protocols in Quantum Annealing and Quantum Approximate Optimization Algorithm Problems

Quantum annealing (QA) and the quantum approximate optimization algorithm (QAOA) are two special cases of the following control problem: apply a combination of two Hamiltonians to minimize the energy of a quantum state. Which is more effective has remained unclear. Here we analytically apply the framework of optimal control theory to show that generically, given a fixed amount of time, the optimal procedure has the pulsed (or “bang-bang”) structure of QAOA at the beginning and end but can have a smooth annealing structure in between. This is in contrast to previous works which have suggested that bang-bang (i.e., QAOA) protocols are ideal. To support this theoretical work, we carry out simulations of various transverse field Ising models, demonstrating that bang-anneal-bang protocols are more common. Futher, the general features identified here provide guideposts for the nascent experimental implementations of quantum optimization algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Inelastic Neutron Scattering Study of Magnetic Exchange Pathways in MnS

We report an investigation of the magnetic structure and magnetic exchange pathways in MnS via neutron scattering methods, aided by density functional theory (DFT) modeling. The material has been confirmed to undergo antiferromagnetic (AFM) ordering at 152 K, with the magnetic structure representing AFM stacking of ferromagnetic (FM) (111) planes of Mn magnetic moments. Correspondingly, the magnetic structure is described by a propagation vector k = (1/2, 1/2, 1/2), with the volume of the magnetic unit cell being 8 times larger than the volume of the nuclear unit cell. Analysis of inelastic neutron scattering (INS) data collected on a powder sample of MnS revealed that the next-nearest-neighbor magnetic exchange constant (J 2 ) exceeds the nearest-neighbor exchange constant (J 1 ) by more than 3 times, while in the case of MnO, which exhibits the same nuclear and magnetic structures as MnS, the J 2 /J 1 ratio was reported to be below 1.5. Although for MnO the signs of both J 1 and J 2 indicated AFM exchange interactions, machine-learning INS data analysis in combination with DFT calculations suggests that the INS data collected on MnS are best described with J 1 < 0 and J 2 > 0, corresponding to FM and AFM exchange couplings, respectively. To achieve a satisfactory fit to the experimentally observed data, the Hamiltonian used to model the INS spectra also included the next-next-nearest-neighbor magnetic exchange constant (J 3 ). The best-fit model has been obtained with the values of the exchange constants J 1 = -0.27 meV, J 2 = 1.05 meV, and J 3 = -0.19 meV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Momentum conservation in current drive and alpha-channeling-mediated rotation drive

Alpha channeling uses waves to extract hot ash from a fusion plasma, transferring energy from the ash to the wave. It has been proposed that this process could create a radial electric field, efficiently driving E × B rotation. However, existing theories ignore the nonresonant particles, which play a critical role in enforcing momentum conservation in quasilinear theory. Because cross field charge transport and momentum conservation are fundamentally linked, this non-consistency throws the rotation drive into question. This paper has two main goals. First, we provide a pedantic and cohesive introduction to the recently developed simple, general, self-consistent quasilinear theory for electrostatic waves that explains the torques which allow for current drive parallel to the magnetic field, and charge extraction across it; a theory that has largely resolved the question of rotation drive by alpha channeling. We show how the theory reveals a fundamental difference between the reaction of nonresonant particles to plane waves that grow in time vs steady-state waves that have a nonuniform spatial structure, allowing rotation drive in the latter case while precluding it in the former, and we review the local and global conservation laws that lead to this result. Second, we provide two new results in support of the theory. First, we provide a novel two-particle Hamiltonian model that rigorously establishes the relationship between charge transport and momentum conservation. Second, we compare the new quasilinear theory to the oscillation-center theories of ponderomotive forces, showing how the latter often obscure the time-dependent nonresonant recoil, but ultimately lead to similar results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Magnetic properties of metastable honeycomb KCoAsO 4

We present comprehensive neutron scattering data on polycrystalline samples of a new metastable honeycomb material KCoAsO 4 . Below 𝑇 𝑁 = 14 K, the system orders into a zigzag antiferromagnetic state, with spins ordered into alternating ferromagnetic chains similar to the isostructural sister compound KNiAsO 4 . In the case of KCoAsO 4 , we find the moments are robustly canted out of plane closer to the crystallographic 𝑐 axis. A combination of weak interlayer coupling, lattice strain, and inhomogeneities lead to the coexistence of two magnetic ordering wave vectors 𝑘 1 = (1.5 0 0) and 𝑘 2 = (0.5 0 0.5), where the two structures differ only in their layer stacking. Inelastic data show the presence of a spin orbital mode at 24 meV, supporting a pseudospin $\tilde{𝑆}$ = 1/2 Kramer's doublet ground state of the Co 2+ ions. We model the low energy excitations using both a conventional XXZ Hamiltonian and generalized Kitaev Heisenberg Hamiltonian within a linear spin wave limit. While either model can qualitatively reproduce the observed spectra, the lack of fine features and observation of disorder prevent a clear-cut determination of the low energy Hamiltonian. In the case of an XXZ-type model, a large easy-axis anisotropy is necessary to reproduce the gapped spectra and canting of the magnetic moments. For the generalized Kitaev model, despite the large canting of the moments away from the honeycomb layers, we find a noticeable if nondominant Kitaev term persists. In conclusion, the contrast of the magnetic properties of KCoAsO 4 to other cobalt honeycombs highlights the sensitivity of the low energy magnetic properties of Co 2+ to fine details of its crystalline environment.

Honeycomb lattice↗