Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lattice Integrable Models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Concentration fluctuations in miscible polymer blends: Influence of temperature and chain rigidity

In contrast to binary mixtures of small molecule fluids, homogeneous polymer blends exhibit relatively large concentration fluctuations that can strongly affect the transport properties of these complex fluids over wide ranges of temperatures and compositions. The spatial scale and intensity of these compositional fluctuations are studied by applying Kirkwood-Buff theory to model blends of linear semiflexible polymer chains with upper critical solution temperatures. The requisite quantities for determining the Kirkwood-Buff integrals are generated from the lattice cluster theory for the thermodynamics of the blend and from the generalization of the random phase approximation to compressible polymer mixtures. We explore how the scale and intensity of composition fluctuations in binary blends vary with the reduced temperature τ ≡ (T − T{sub c})/T (where T{sub c} is the critical temperature) and with the asymmetry in the rigidities of the components. Knowledge of these variations is crucial for understanding the dynamics of materials fabricated from polymer blends, and evidence supporting these expectations is briefly discussed.

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

Efficient field-theoretic simulation of polymer solutions

We present several developments that facilitate the efficient field-theoretic simulation of polymers by complex Langevin sampling. A regularization scheme using finite Gaussian excluded volume interactions is used to derive a polymer solution model that appears free of ultraviolet divergences and hence is well-suited for lattice-discretized field theoretic simulation. We show that such models can exhibit ultraviolet sensitivity, a numerical pathology that dramatically increases sampling error in the continuum lattice limit, and further show that this pathology can be eliminated by appropriate model reformulation by variable transformation. We present an exponential time differencing algorithm for integrating complex Langevin equations for field theoretic simulation, and show that the algorithm exhibits excellent accuracy and stability properties for our regularized polymer model. These developments collectively enable substantially more efficient field-theoretic simulation of polymers, and illustrate the importance of simultaneously addressing analytical and numerical pathologies when implementing such computations.

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

Nonlinear Optics at UMER: Lessons Learned in Simulation

Design of accelerator lattices with nonlinear optics to suppress transverse resonances is a novel approach and may be crucial for enabling low-loss high-intensity beam transport. Large amplitude-dependent tune spreads, driven by nonlinear field inserts, damp resonant response to driving terms. This presentation will focus on simulations of the UMER lattice operated as a quasi-integrable system (1 invariant of transverse motion) with a single strong octupole insert. We will discuss the evolution of simulation models, including the observation of losses associated with the original operating point near a fourth-order resonance. Other operating points farther from this resonance are considered and shown to be more promising.

Ruisard, Kiersten↗
DOE OSTI2025

Utility of a hybrid approach to the hadronic vacuum polarization contribution to the muon anomalous magnetic moment

An accurate calculation of the leading-order hadronic vacuum polarization (LOHVP) contribution to the anomalous magnetic moment of the muon ( a μ ) is key to determining whether a discrepancy, suggesting new physics, exists between the Standard Model and experimental results. This calculation can be expressed as an integral over Euclidean time of a current-current correlator G ( t ) , where G ( t ) can be calculated using lattice QCD or, with dispersion relations, from experimental data for e + e − → hadrons . The BMW/DMZ collaboration recently presented a hybrid approach in which G ( t ) is calculated using lattice QCD for most of the contributing t range, but using experimental data for the largest t (lowest energy) region. Here we study the advantages of varying the position t = t 1 separating lattice QCD from data-driven contributions. The total LOHVP contribution should be independent of t 1 , providing both a test of the experimental input and the robustness of the hybrid approach. We use this criterion and a correlated fit to show that Fermilab/HPQCD/MILC lattice QCD results from 2019 strongly favor the CMD-3 cross section data for e + e − → π + π − over a combination of earlier experimental results for this channel. Further, the resulting total LOHVP contribution obtained is consistent with the result obtained by BMW/DMZ, and supports the scenario in which there is no significant discrepancy between the experimental value for a μ and that expected in the Standard Model. We then discuss how improved lattice results in this hybrid approach could provide a more accurate total LOHVP across a wider range of t 1 values with an uncertainty that is smaller than that from either lattice QCD or data-driven approaches on their own. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗
DOE OSTI2016

Motion of a distinguishable Impurity in the Bose gas: Arrested expansion without a lattice and impurity snaking

We consider the real-time dynamics of an initially localized distinguishable impurity injected into the ground state of the Lieb-Liniger model. Focusing on the case where integrability is preserved, we numerically compute the time evolution of the impurity density operator in regimes far from analytically tractable limits. We find that the injected impurity undergoes a stuttering motion as it moves and expands. For an initially stationary impurity, the interaction-driven formation of a quasibound state with a hole in the background gas leads to arrested expansion—a period of quasistationary behavior. In conclusion, when the impurity is injected with a finite center-of-mass momentum, the impurity moves through the background gas in a snaking manner, arising from a quantum Newton’s cradlelike scenario where momentum is exchanged back and forth between the impurity and the background gas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

First-principles thermodynamic assessment of Sr-containing secondary phase formation in strontium-substituted lanthanum manganites for solid oxide cell applications

Sr-secondary phase formation is a potentially significant degradation mode with direct impact upon solid-oxide cell (SOC) commercial viability. A first-principles based thermodynamic study was performed for La 1−x Sr x MnO 3±δ (LSM) perovskites to assess their stability against formation of different Sr-secondary phases, including SrO, SrCrO 4 , SrSO 4 , SrCO 3 , and Sr(OH) 2 , for SOC applications. The Sr-secondary phase formation reaction free energies were determined via a thermodynamic model by combining ab initio lattice dynamics calculations for the solid phases and ab initio thermodynamic data for the gas phases. The current approach expands the previously reported thermodynamic modeling studies by integrating first-principles based point defect equilibria into the thermodynamic analysis. The modeling results obtained using this new approach indicate an increased tendency to form the SrO oxide upon decreasing the oxygen partial pressure. Additionally, the enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with the contaminant gas species as a function of temperature and pressure.

defect thermodynamics modeling↗
DOE OSTI2018

Importance sampling large deviations in nonequilibrium steady states. I

Large deviation functions hold information on the stability and response of systems driven into nonequilibrium steady states and in such a way are similar to free energies for systems at equilibrium. As with equilibrium free energies, evaluating large deviation functions numerically for all but the simplest systems is difficult because by construction they depend on exponentially rare events. In this first paper of a series, we assess different trajectory-based sampling methods capable of computing large deviation functions of time integrated observables within nonequilibrium steady states. We highlight some convergence criteria and best practices using a number of different models, including a biased Brownian walker, a driven lattice gas, and a model of self-assembly. We show how two popular methods for sampling trajectory ensembles, transition path sampling and diffusion Monte Carlo, suffer from exponentially diverging correlations in trajectory space as a function of the bias parameter when estimating large deviation functions. Improving the efficiencies of these algorithms requires introducing guiding functions for the trajectories.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Structure of the Ndc80 complex and its interactions at the yeast kinetochore–microtubule interface

The conserved Ndc80 kinetochore complex, Ndc80c, is the principal link between mitotic spindle microtubules and centromere-associated proteins. We used AlphaFold 2 (AF2) to obtain predictions of the Ndc80 ‘loop’ structure and of the Ndc80 : Nuf2 globular head domains that interact with the Dam1 subunit of the heterodecameric DASH/Dam1 complex (Dam1c). The predictions guided design of crystallizable constructs, with structures close to the predicted ones. The Ndc80 ‘loop’ is a stiff, α-helical ‘switchback’ structure; AF2 predictions and positions of preferential cleavage sites indicate that flexibility within the long Ndc80c rod occurs instead at a hinge closer to the globular head. Conserved stretches of the Dam1 C terminus bind Ndc80c such that phosphorylation of Dam1 serine residues 257, 265 and 292 by the mitotic kinase Ipl1/Aurora B can release this contact during error correction of mis-attached kinetochores. We integrate the structural results presented here into our current molecular model of the kinetochore–microtubule interface. The model illustrates how multiple interactions between Ndc80c, DASH/Dam1c and the microtubule lattice stabilize kinetochore attachments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

First-Principles Thermodynamic Assessments of Sr-Containing Secondary Phase Formation in La1-xSrxMnO3±δ Perovskites for Solid Oxide Cell Applications

Sr-secondary phase formation is a potentially significant degradation mode threatening solid-oxide cell (SOC) commercial viability. A first-principles thermodynamic study was performed for rhombohedral perovskite (La1-xSrx) MnO3±δ (LSM) to assess its stability against Sr secondary phase formation in SOC applications. In this work, the Sr secondary phase formation reaction free energies were determined by combining ab initio lattice dynamics calculations for the solid phases and an ab initio thermodynamics approach for the gas phases. Furthermore, this approach goes beyond previous thermodynamic modeling studies by integrating first-principles based point-defect equilibria into the analyses. The modeling results indicate an increased tendency to form SrO oxide from LSM upon decreasing the oxygen partial pressure. Additionally, enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with contaminant gas species as a function of temperature and gas pressure.

Defect and phase stability↗
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 OSTI2025

Deciphering the Scattering of Mechanically Driven Polymers Using Deep Learning

Here, we present a deep learning approach for analyzing two-dimensional scattering data of semiflexible polymers under external forces. In our framework, scattering functions are compressed into a three-dimensional latent space using a Variational Autoencoder (VAE), and two converter networks establish a bidirectional mapping between the polymer parameters (bending modulus, stretching force, and steady shear) and the scattering functions. The training data are generated using off-lattice Monte Carlo simulations to avoid the orientational bias inherent in lattice models, ensuring robust sampling of polymer conformations. The feasibility of this bidirectional mapping is demonstrated by the organized distribution of polymer parameters in the latent space. By integrating the converter networks with the VAE, we obtain a generator that produces scattering functions from given polymer parameters and an inferrer that directly extracts polymer parameters from scattering data. While the generator can be utilized in a traditional least-squares fitting procedure, the inferrer produces comparable results in a single pass and operates 3 orders of magnitude faster. This approach offers a scalable automated tool for polymer scattering analysis and provides a promising foundation for extending the method to other scattering models, experimental validation, and the study of time-dependent scattering data.

Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2013

Chiral symmetry and fermion doubling in the zero-mode Landau levels of massless Dirac fermions with disorder

The effect of disorder on the Landau levels of massless Dirac fermions is examined for the cases with and without the fermion doubling. To tune the doubling a tight-binding model having a complex transfer integral is adopted to shift the energies of two Dirac cones, which is theoretically proposed earlier and realizable in cold atoms in an optical lattice. In the absence of the fermion doubling, the n = 0 Landau level is shown to exhibit an anomalous sharpness even if the disorder is uncorrelated in space (i.e., large K-K' scattering). This anomaly occurs when the disorder respects the chiral symmetry of the Dirac cone.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Phenomenology of Semileptonic B-Meson Decays with Form Factors from Lattice QCD

We study the exclusive semileptonic B-meson decays B→K(π)ℓ + ℓ - , B→K(π)νν¯, and B→πτν, computing observables in the Standard model using the recent lattice-QCD results for the underlying form factors from the Fermilab Lattice and MILC Collaborations. These processes provide theoretically clean windows into physics beyond the Standard Model because the hadronic uncertainties are now under good control. The resulting partially-integrated branching fractions for B→πμ + μ - and B→Kμ + μ - outside the charmonium resonance region are 1-2σ higher than the LHCb Collaboration's recent measurements, where the theoretical and experimental errors are commensurate. The combined tension is 1.7σ. Combining the Standard-Model rates with LHCb's measurements yields values for the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements |V td |=7.45(69)×10 -3 , |V ts |=35.7(1.5)×10 -3 , and |V td /V ts |=0.201(20), which are compatible with the values obtained from neutral B (s) -meson oscillations and have competitive uncertainties. Alternatively, taking the CKM matrix elements from unitarity, we constrain new-physics contributions at the electroweak scale. Furthermore, the constraints on the Wilson coefficients Re(C 9 ) and Re(C 10 ) from B→πμ + μ - and B→Kμ + μ - are competitive with those from B→K*μ + μ - , and display a 2.0σ tension with the Standard Model. Our predictions for B→K(π)νν¯ and B→πτν are close to the current experimental limits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Diffraction-geometry refinement in the DIALS framework

Rapid data collection and modern computing resources provide the opportunity to revisit the task of optimizing the model of diffraction geometry prior to integration. A comprehensive description is given of new software that builds upon established methods by performing a single global refinement procedure, utilizing a smoothly varying model of the crystal lattice where appropriate. This global refinement technique extends to multiple data sets, providing useful constraints to handle the problem of correlated parameters, particularly for small wedges of data. Examples of advanced uses of the software are given and the design is explained in detail, with particular emphasis on the flexibility and extensibility it entails.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Unveiling the mechanism of remote epitaxy of crystalline semiconductors on 2D materials-coated substrates

Abstract Remote epitaxy has opened novel opportunities for advanced manufacturing and heterogeneous integration of two-dimensional (2D) materials and conventional (3D) materials. The lattice transparency as the fundamental principle of remote epitaxy has been studied and challenged by recent observations defying the concept. Understanding remote epitaxy requires an integrated approach of theoretical modeling and experimental validation at multi-scales because the phenomenon includes remote interactions of atoms across an atomically thin material and a few van der Waals gaps. The roles of atomically thin 2D material for the nucleation and growth of a 3D material have not been integrated into a framework of remote epitaxy research. Here, we summarize studies of remote epitaxy mechanisms with a comparison to other epitaxy techniques. In the end, we suggest the crucial topics of remote epitaxy research for basic science and applications. Graphical Abstract

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology for pin lattice models is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. The user can alternatively provide their own geometry and power profile instead of inheriting it from DIF3D. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and a sufficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1985

Adaptive control of large space structures using recursive lattice filters

The use of recursive lattice filters for identification and adaptive control of large space structures is studied. Lattice filters were used to identify the structural dynamics model of the flexible structures. This identification model is then used for adaptive control. Before the identified model and control laws are integrated, the identified model is passed through a series of validation procedures and only when the model passes these validation procedures is control engaged. This type of validation scheme prevents instability when the overall loop is closed. Another important area of research, namely that of robust controller synthesis, was investigated using frequency domain multivariable controller synthesis methods. The method uses the Linear Quadratic Guassian/Loop Transfer Recovery (LQG/LTR) approach to ensure stability against unmodeled higher frequency modes and achieves the desired performance.

Sundararajan, N.↗