Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lattice calculation”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Symmetrizing relativistic three-body partial wave amplitudes

S matrix principles and symmetries impose constraints on three-particle scattering amplitudes, which can be formulated as a class of integral equations for their partial wave projections. However, these amplitudes are typically expressed in an asymmetric basis, where one of the initial and final state particles is singled out, and all quantum numbers are defined relative to this spectator. In this work, we show how to construct symmetric partial wave amplitudes, which have been symmetrized over all possible spectator combinations, using their asymmetric counterparts and sets of recoupling coefficients. We derive these recoupling coefficients for arbitrary angular momentum and isospin for arbitrary systems of spinless particles with SU(2) flavor symmetry. We propose a simple intensity observable suitable for visualizing the structure of three-body dynamics in Dalitz distributions. Finally, we provide some numerical examples of Dalitz distributions relevant for future lattice QCD calculations of 3⁢𝜋 systems and provide numerical evidence that the symmetrization procedure is consistent with expected symmetries of Dalitz plots.

Lattice QCD↗
DOE OSTI2026

Implementing the finite-volume three-pion scattering formalism across all non-maximal isospins

We present a numerical exploration of the relativistic-field-theory (RFT) formalism for three pions with all possible values of non-maximal isospin, I πππ = 2, 1 and 0. Using the generic-isospin extension of the RFT formalism [1] and applying our open-source Python library to implement the framework, we predict a range of three-pion energies for illustrative values of the two-to-two scattering amplitudes for various finite-volume irreps also with non-zero total momentum P in the finite-volume frame. The results restrict attention to the case of a vanishing intrinsic three-body interaction so that the spectra can be understood as a baseline. In future lattice QCD calculations, deviations from these values will be translated into evidence for intrinsic three-body effects in the various scattering channels.

hadronic spectroscopy↗
DOE OSTI2025

Sequential multidimensional heteroepitaxy of chalcogen-sharing 3D ZnSe and 2D MoSe 2 with quasi van der Waals interface engineering

Two-dimensional (2D) materials are emerging as a promising platform for epitaxial growth, largely free from the constraints of lattice constant and thermal expansion coefficient mismatches. Among them, transition metal dichalcogenides (TMDs), known for their superior electrical properties, are ideal for ultrathin semiconductor applications. Their unique epitaxial characteristics enable seamless integration with 3D materials, facilitating the development of gate stacks and heterojunction devices. In this regard, developing a process for growing high-quality 3D epitaxial materials before and after the growth of 2D TMDs and understanding the 2D/3D interface are crucial. This study demonstrates the sequential growth of fully epitaxial ZnSe/MoSe 2 /ZnSe heterostructures using metal-organic chemical vapor deposition. ZnSe and MoSe 2 , sharing chalcogen elements, enable large-area quasi van der Waals epitaxy with sharp interfaces without intermediate phase. Multiscale analysis involving transmission electron microscopy and density functional theory calculation reveals lattice commensurability, van der Waals gaps, termination, and interfacial reconstruction. Understanding these interactions is crucial for advancing multidimensional integration of 2D and 3D materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electronic Structure Theory and Novel Materials

This grant supported research on electronic structure and materials theory, with focus on three main issues: (i) novel techniques to deal with correlation in the electronic ground-state, (ii) topological materials, (iii) the phase diagram of lattice spin models. Regarding (i), we applied to the homogeneous electron liquid an approach that we previously developed in the context of molecular systems. In this scheme the electronic occupation probabilities and the natural spin orbitals are used to construct an approximate two-body density matrix for the electronic ground-state. Regarding (ii) we used standard electronic structure methods based on density functional theory to model topological materials and interpret experimental observations. Finally, regarding (iii) we further developed a numerical approach to compute the renormalized couplings within real space renormalization group theory in the context of lattice spin models. The main findings were the following. (i) We found that with our approximate two-body density matrix, which works well for small molecules, is not sufficiently accurate for condensed phase systems. Missing a systematic way of improving on the adopted approximations, we decided not to pursue this approach. (ii) We performed two studies. In one, we investigated the influence of Te defects on the topological properties of a WTe2 monolayer, finding that while Te vacancies, even in modest concentration, destroy the topological character, Te adatoms do not, consistent with a recent experiment. In another study, we predicted Weyl semimetal character and strong anomalous Hall effect in the Heusler compensated ferrimagnet Ti2MnAl. (iii) We developed a new Monte Carlo method to do real space renormalization group calculations for lattice spin models. We subsequently extended the scheme to deal with lattice spin models in presence of quenched disorder, finding that the approach can distinguish systems with finite and strong disorder. In the finite disorder case, the method allows one to find with good approximation the critical coupling distribution and the critical exponents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Generalized parton distributions through universal moment parameterization: non-zero skewness case

We present the first global analysis of generalized parton distributions (GPDs) combing lattice quantum chromodynamics (QCD) calculations and experiment measurements including global parton distribution functions (PDFs), form factors (FFs) and deeply virtual Compton scattering (DVCS) measurements. Following the previous work where we parameterize GPDs in terms of their moments, we extend the framework to allow for the global analysis at non-zero skewness. Together with the constraints at zero skewness, we fit GPDs to global DVCS measurements from both the recent JLab and the earlier Hadron-Electron Ring Accelerator (HERA) experiments with two active quark flavors and leading order QCD evolution. With certain choices of empirical constraints, both sea and valence quark distributions are extracted with the combined inputs, and we present the quark distributions in the proton correspondingly. We also discuss how to extend the framework to accommodate more off-forward constraints beyond the small ξ expansion, especially the lattice calculated GPDs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Quark spin and orbital angular momentum from proton generalized parton distributions

Here, we calculate the leading-twist helicity-dependent generalized parton distributions (GPDs) of the proton at finite skewness in the Nambu–Jona-Lasinio (NJL) model of quantum chromodynamics (QCD). From these (and previously calculated helicity-independent GPDs) we obtain the spin decomposition of the proton, including predictions for quark intrinsic spin and orbital angular momentum. The inclusion of multiple species of diquarks is found to have a significant effect on the flavor decomposition, and resolving the internal structure of these dynamical diquark correlations proves essential for the mechanical stability of the proton. At a scale of Q 2 = 4 GeV 2 we find that the up and down quarks carry an intrinsic spin and orbital angular momentum of S u = 0.534, S d = -0.214, L u = -0.189, and L d = 0.210 , whereas the gluons have a total angular momentum of J g = 0.151 . The down quark is therefore found to carry almost no total angular momentum due to cancellations between spin and orbital contributions. Comparisons are made between these spin decomposition results and lattice QCD calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Solvent Isotherms and Structural Transitions in Nanoparticle Superlattice Assembly

Here, we introduce a Molecular Theory for Compressible Fluids (MOLT-CF) that enables us to compute free energies and other thermodynamic functions for nanoparticle superlattices with any solvent content, including the dry limit. Quantitative agreement is observed between MOLT-CF and united-atom molecular dynamics simulations performed to assess the reliability and precision of the theory. Among other predictions, MOLT-CF shows that the amount of solvent within the superlattice decreases approximately linearly with its vapor pressure and that in the late stages of drying, solvent-filled voids form at lattice interstitials. Applied to single-component superlattices, MOLT-CF predicts fcc-to-bcc Bain transitions for decreasing vapor pressure and for increasing ligand length, both in agreement with experimental results. We explore the stability of other single-component phases and show that the C14 Frank-Kasper phase, which has been reported in experiments, is not a global free-energy minimum. Implications for precise assembly and prediction of multicomponent nanoparticle systems are discussed.

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

Spectrum of Dirac eigenvalues across QCD chiral transition

Here, we show that the universal scaling properties of the chiral phase transition in quantum chromodynamics (QCD) at the macroscale are encoded in the microscopic energy levels of the quarks in the background of gluon fields, defined as eigenvalues of the Dirac operator. We introduce a novel relation between the cumulants of the chiral condensate, and the correlations among the energy levels of quarks, i.e., the eigenspectra of the massless QCD Dirac operator. This relation elucidates how the fluctuations of the chiral condensate arise from the correlations within the infrared part of the energy spectra of quarks, and naturally leads to a generalization of the Banks-Casher relation for the cumulants of the chiral condensate. Then, through (2+1)-flavor lattice QCD calculations using HISQ action with varying light quark masses around the chiral phase transition temperature, we demonstrate that the correlations among the infrared part of the Dirac eigenvalue spectra exhibit same universal scaling behaviors as expected of the cumulants of the chiral condensate. We find that these universal scaling behaviors extend up to the physical values of the up and down quark masses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

High-Field Magnetoelectric and Spin-Phonon Coupling in Multiferroic (NH 4 ) 2 [FeCl 5 ·(H 2 O)]

We combine high field polarization, magneto-infrared spectroscopy, and lattice dynamics calculations with prior magnetization to explore the properties of (NH 4 ) 2 [FeCl 5 ·(H 2 O)]-a type II molecular multiferroic in which the mixing between charge, structure, and magnetism is controlled by intermolecular hydrogen and halogen bonds. Electric polarization is sensitive to the series of field-induced spin reorientations, increasing linearly with the field and reaching a maximum before collapsing to zero across the quasi-collinear to collinear-sinusoidal reorientation due to the restoration of inversion symmetry. Magnetoelectric coupling is on the order of 1.2 ps/m for the P∥c, H∥c configuration between 5 and 25 T at 1.5 K. In this range, the coupling takes place via an orbital hybridization mechanism. Other forms of mixing are active in (NH 4 ) 2 [FeCl 5 ·(H 2 O)] as well. Magneto-infrared spectroscopy reveals that all of the vibrational modes below 600 cm –1 are sensitive to the field-induced transition to the fully saturated magnetic state at 30 T. We analyze these local lattice distortions and use frequency shifts to extract spin-phonon coupling constants for the Fe–O stretch, Fe–OH 2 rock, and NH 4 + libration. Inspection also reveals subtle symmetry breaking of the ammonium counterions across the ferroelectric transition. Here, the coexistence of such varied mixing processes in a platform with intermolecular hydrogen- and halogen-bonding opens the door to greater understanding of multiferroics and magnetoelectrics governed by through-space interactions.

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

Effects of rhodium doping on dislocation nucleation in a [001] SrNi 2 P 2 single crystal under spherical nanoindentation

Nanoindentation was performed on SrNi 2 P 2 single crystals and their dilute solid solutions, Sr(Ni 1–x Rh x ) 2 P 2 (x = 0 ~ 0.055), under [001] loading to investigate the effects of elemental doping on Young’s modulus and dislocation nucleation stress. The results show that Young’s modulus and the dislocation nucleation stress decrease as the Rh content increases, and their Rh content dependence also varies with the Rh content. Electron diffraction analysis and the anisotropic lattice distortion calculations revealed that a local structural transition from the orthorhombic superstructure to the tetragonal structure and local compressive residual stresses could be the main reasons for the decrease in Young’s modulus and dislocation nucleation stress, respectively. The overlap of stress–strain fields between Rh atoms and the associated local structural transition contribute to their Rh content dependence. Here, the results also indicate that elemental doping does not necessarily strengthen materials if the incipient plasticity is controlled by dislocation nucleation, unlike conventional solid solution strengthening.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Di-nucleons do not form bound states at heavy pion mass

We perform a high-statistics lattice QCD calculation of the low-energy two-nucleon scattering amplitudes. In order to address discrepancies in the literature, the calculation is performed at a heavy pion mass in the limit that the light quark masses are equal to the physical strange quark mass, $m_π= m_K \simeq 714 $ MeV. Using a state-of-the-art momentum space method, we rule out the presence of a bound di-nucleon in both the isospin 0 (deuteron) and 1 (di-neutron) channels, in contrast with many previous results that made use of compact hexaquark creation operators. In order to diagnose the discrepancy, we add such hexaquark interpolating operators to our basis and find that they do not affect the determination of the two-nucleon finite volume spectrum, and thus they do not couple to deeply bound di-nucleons that are missed by the momentum-space operators. Further, we perform a high-statistics calculation of the HAL QCD potential on the same gauge ensembles and find qualitative agreement with our main results. We conclude that two-nucleons do not form bound states at heavy pion masses and that previous identification of deeply bound di-nucleons must have arisen from a misidentification of the spectrum from off-diagonal elements of a correlation function.

FOS: Physical sciences↗
DOE OSTI2022

Rethinking the ill-posedness of the spectral function reconstruction — Why is it fundamentally hard and how Artificial Neural Networks can help

Reconstructing hadron spectral functions through Euclidean correlation functions are of the important missions in lattice QCD calculations. However, in a Källen–Lehmann (KL) spectral representation, the reconstruction is observed to be ill-posed in practice. It is usually ascribed to the fewer observation points compared to the number of points in the spectral function. In this paper, by solving the eigenvalue problem of continuous KL convolution, we show analytically that the ill-posedness of the inversion is fundamental and it exists even for continuous correlation functions. We discussed how to introduce regulators to alleviate the predicament, in which include the Artificial Neural Networks (ANNs) representations recently proposed by the Authors in another study. The uniqueness of solutions using ANNs representations is manifested analytically and validated numerically. Reconstructed spectral functions using different regularization schemes are also demonstrated, together with their eigen-mode decomposition. We observe that components with large eigenvalues can be reliably reconstructed by all methods, whereas those with low eigenvalues need to be constrained by regulators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Machine learning to alleviate Hubbard-model sign problems

Lattice Monte Carlo calculations of interacting systems on nonbipartite lattices exhibit an oscillatory imaginary phase known as the phase or sign problem, even at zero chemical potential. One method to alleviate the sign problem is to analytically continue the integration region of the state variables into the complex plane via holomorphic flow equations. For asymptotically large flow times, the state variables approach manifolds of constant imaginary phase known as Lefschetz thimbles. Furthermore, flowing such variables and calculating the ensuing Jacobian is a computationally demanding procedure. In this paper, we demonstrate that neural networks can be trained to parametrize suitable manifolds for this class of sign problem and drastically reduce the computational cost for different severely afflicted small volume systems. In particular, we apply our method to the Hubbard model on the triangle and tetrahedron, both of which are nonbipartite. At strong interaction strengths and modest temperatures, the tetrahedron suffers from a severe sign problem that cannot be overcome with standard reweighting techniques, while it quickly yields to our method. We benchmark our results with exact calculations and comment on future directions of this work.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Spin–phonon interactions and magnetoelectric coupling in Co 4 B 2 O 9 ( B = Nb, Ta)

In order to explore the consequences of spin–orbit coupling on spin–phonon interactions in a set of chemically similar mixed metal oxides, we measured the infrared vibrational properties of Co 4 B 2 O 9 (B = Nb, Ta) as a function of temperature and compared our findings with lattice dynamics calculations and several different models of spin–phonon coupling. Frequency vs temperature trends for the Co 2+ shearing mode near 150 cm –1 reveal significant shifts across the magnetic ordering temperature that are especially large in relative terms. Bringing these results together and accounting for noncollinearity, we obtain spin–phonon coupling constants of –3.4 and –4.3 cm –1 for Co 4 Nb 2 O 9 and the Ta analog, respectively. Analysis reveals that these coupling constants are derived from interlayer (rather than intralayer) exchange interactions and that the interlayer interactions contain competing antiferromagnetic and ferromagnetic contributions. At the same time, beyond-Heisenberg terms are minimized due to fortuitous symmetry considerations, different from most other 4d- and 5d-containing oxides. Comparison with other contemporary oxides shows that spin–phonon coupling in this family of materials is among the strongest ever reported, suggesting an origin for magnetoelectric coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Symmetry progression and band vs Mott character of CdPS 3 under pressure

Complex chalcogenides are renowned for their tunable electronic, magnetic, and optical properties under external stimuli. The MPX 3 family (M = Mn, Ni, Co, V; X = S, Se) is a platform for many exciting discoveries—especially under compression—although CdPS3 is thought to be different because the Cd center possesses a filled 4d shell, which precludes Mottness. Here, we combine synchrotron-based infrared absorbance and Raman scattering spectroscopies with diamond anvil cell techniques, complementary lattice dynamics calculations, and an analysis of the energy landscape to reveal a series of structural phase transitions in CdPS 3 . We find four distinct pressure-driven transitions, with low frequency modes detectable over the full 35 GPa range of our investigation. A group–subgroup analysis along with our first-principles calculations allows us to partially unravel the space group sequence. For instance, the first critical pressure is a monoclinic C2/m to trigonal $\overline{R}$3 transition at 10 GPa. Despite the softness and overall sensitivity to pressure, we do not locate an insulator-to-metal transition in this pressure range, indicating that the energy scale for gap closure is significantly higher than expected. We discuss these findings in terms of force-induced color change and Mott vs band character in this system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Multiple Lattice Instabilities and Complex Ground State in Cs 2 Ag Bi Br 6

Metal-halide perovskites (MHPs) are attracting considerable interest for optoelectronic applications, with Cs 2 Ag Bi Br 6 one of the main contenders among lead-free systems. Cs 2 Ag Bi Br 6 crystallizes in a nominally double-perovskite structure, but exhibits a soft lattice with large atomic fluctuations characteristic of MHPs. While crucial to understand electron-phonon and phonon-phonon couplings, the spatiotemporal correlations of these fluctuations remain largely unknown. Here, we reveal these correlations using comprehensive neutron and x-ray scattering measurements on Cs 2 Ag Bi Br 6 single crystals, complemented with first-principles simulations augmented with machine-learned neural-network potentials. We report the discovery of an unexpected complex modulated ground-state structure containing several hundred atoms, arising from a soft-phonon instability of the low-temperature tetragonal phase. Further, our experiments and simulations both reveal extensive correlated two-dimensional fluctuations of Br octahedra at finite temperature, arising from soft optic phonons that are strongly broadened by anhamonicity, reflecting very shallow potential wells. These results provide new insights into the atomic structure and fluctuations in MHPs, critical to understand and control their thermal and optoelectronic properties. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Lattice QCD Benchmark of Proton Helicity and Flavor-Dependent Unpolarized Transverse Momentum-Dependent Parton Distribution Functions at Physical Quark Masses

We present the first lattice QCD calculations of the isovector helicity transverse momentum-dependent parton distribution function (TMDPDF) and the flavor-dependent unpolarized TMDPDFs for up and down quarks in the proton. Our computations utilize domain-wall fermion discretization with physical quark masses. Employing Coulomb-gauge-fixed quark correlation functions within the large-momentum effective theory framework, we access nonperturbative transverse quark separations 𝑏 𝑇 up to approximately 1 fm, corresponding to transverse momenta as low as 200 MeV. Based on the quasi-TMD factorization theorem, we construct renormalization-group–invariant ratios that are equal to the corresponding light-cone TMDPDF ratios. At moderate 𝑥, our results reveal that the isovector helicity and unpolarized TMDPDFs exhibit nearly identical transverse structure up to a normalization factor, and the unpolarized distributions display only mild flavor dependence. These findings not only support key trends observed in recent global analyses but also provide robust, nonperturbative constraints that can distinguish between different parameterizations. This Letter establishes a first-principles benchmark for TMDPDFs, offering valuable input for ongoing and future experimental efforts to map the proton’s three-dimensional structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗