Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Towards a Quintic Ginzburg-Landau Description of the (2,7) Minimal Model

We discuss dimensional continuation of the massless scalar field theory with the 𝑖⁢𝜙 5 interaction term. It preserves the so-called 𝒫⁢𝒯 symmetry, which acts by 𝜙 →−𝜙 accompanied by 𝑖 →−𝑖. Below its upper critical dimension 10/3, this theory has interacting infrared fixed points. We argue that the fixed point in 𝑑 = 2 describes the nonunitary minimal conformal model 𝑀⁡(2,7). We identify the operators 𝜙 and 𝜙 2 with the Virasoro primaries 𝜙 1,2 and 𝜙 1,3 , respectively, and 𝑖⁢𝜙 3 with a quasiprimary operator, which is a Virasoro descendant of 𝜙 1,3 . Our identifications appear to be consistent with the operator product expansions and with considerations based on integrability. Using constrained Padé extrapolations, we provide estimates of the critical exponents in 𝑑 = 3. We also comment on possible lattice descriptions of 𝑀⁡(2,7) and discuss RG flows to and from this conformal field theory (CFT). Finally, we conjecture that the minimal models 𝑀⁡(2,2⁢𝑛 +1) are described by the massless scalar field theories with the 𝑖⁢𝜙 2⁢𝑛−1 interaction terms.

Conformal field theory↗
DOE OSTI2021

A dissolution model of alite coupling surface topography and ions transport under different hydrodynamics conditions at microscale

Portland cement is the most produced material in the world. The hydration process of cement consists of a group of complex chemical reactions. In order to investigate the mechanism of cement hydration, it is vital to study the hydration of each phase separately. An integrated model is proposed in this paper to simulate the dissolution of alite under different hydrodynamic conditions at microscale, coupling Kinetic Monte Carlo model (KMC), Lattice Boltzmann method (LBM) and diffusion boundary layer (DBL). The dissolution of alite is initialised with KMC. Two Multiple-relaxation-time (MRT) LB models are used to simulate the fluid flow and transport of ions, respectively. For solid-liquid interface, DBL is adapted to calculate the concentration gradient and dissolution flux. The model is validated with experiment from literature. The simulation results show good agreements with the results published in the literature.

36 MATERIALS SCIENCE↗
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 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 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 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↗
DOE OSTI2024

String-based parametrization of nucleon GPDs at any skewness: A comparison to lattice QCD

We introduce a string-based parametrization for nucleon quark and gluon generalized parton distributions (GPDs) valid at all skewness values. The conformal moments of the GPDs are expressed as sums of the spin-j nucleon A-form factor and the skewness-dependent spin-j nucleon D-form factor. This representation, which fulfills the polynomiality condition (due to Lorentz invariance) and does not rely on model-specific assumptions, is derived from t-channel string exchanges in anti-de Sitter spaces. The spin-j nucleon D-form factor is closely related to the spin-j nucleon A-form factor. We use the Mellin moments from empirical parton distributions to model the spin-j nucleon A-form factors. Using only five Regge slope parameters, fixed from the electromagnetic and gravitational form factors, our string-based parametrization generates accurate singlet, nonsinglet, isovector, and flavor-separated nucleon quark GPDs, along with symmetric nucleon gluon GPDs from their Mellin-Barnes integral representations. Our isovector nucleon quark GPD is in agreement with existing lattice data. Our string-based parametrization should facilitate the empirical extraction and global analysis of nucleon GPDs in exclusive processes, bypassing the deconvolution challenge.

Electron-ion collisions↗
DOE OSTI2026

Advanced measurement techniques in quantum Monte Carlo: The permutation matrix representation approach

In a typical finite temperature quantum Monte Carlo (QMC) simulation, estimators for simple static observables such as specific heat and magnetization are known. With a great deal of system-specific manual labor, one can sometimes also derive more complicated non-local or even dynamic observable estimators. In contrast, we show that arbitrary static observables can be estimated within the permutation matrix representation (PMR) flavor for any Hamiltonian. We then generalize these results to general imaginary-time correlation functions and non-trivial integrated susceptibilities thereof. Finally, we demonstrate the practical versatility of our method by estimating various non-local, random observables for the transverse-field Ising model on a square lattice and a toy random model.

Permutation matrix representation↗
DOE OSTI2023

Mesoscopic critical fluctuations

Here, we investigate the magnetic fluctuations in a mesoscopic critical region formed at the interface due to smooth time-independent spatial variations of a control parameter around its critical value. In the proximity of the spatial critical point, the order parameter fluctuations exhibit a mesoscopic nature, characterized by their significant size compared to the lattice constant, while gradually decaying away from the critical region. To explain this phenomenon, we present a minimal model that effectively captures this behavior and demonstrates its connection to the integrable Painlevé-II equation governing the local order parameter. By leveraging the well-established mathematical properties of this equation, we gain valuable insights into the nonlinear susceptibilities exhibited within this region.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Improving modular bootstrap bounds with integrality

We propose methods that efficiently impose integrality — i.e., the condition that the coefficients of characters in the partition function must be integers — into numerical modular bootstrap. We demonstrate the method with a number of examples where it can be used to strengthen modular bootstrap results. First, we show that, with a mild extra assumption, imposing integrality improves the bound on the maximal allowed gap in dimensions of operators in theories with a U(1) c symmetry at c = 3, and reduces it to the value saturated by the SU(4) 1 WZW model point of c = 3 Narain lattices moduli space. Second, we show that our method can be used to eliminate all but a discrete set of points saturating the bound from previous Virasoro modular bootstrap results. Finally, when central charge is close to 1, we can slightly improve the upper bound on the scaling dimension gap.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Photoinduced Melting of V 4 O 7 Correlated State

The compound V 4 O 7 is one of the Magnéli phase (V n O 2n - 1 , n = 3, 4, …, 9) correlated vanadium oxides with distinct intriguing electronic and structural properties. The possibility to manipulate the phase state of V 4 O 7 on an ultrafast time scale by light makes this material promising for potential applications in photonics, optoelectronics, quantum, and neuromorphic circuit design. In this work, the ultrafast spectroscopy of V 4 O 7 reveals the second-order nature of the photoinduced insulator-to-metal transition, emphasizing electronic and lattice contributions. The findings reveal the influence of the laser excitation level and temperature on these dynamics, providing a comprehensive understanding of V 4 O 7 structural changes and response to external stimuli. The phenomenological model based on the Landau–Ginzburg formalism provides a robust framework for explaining the photoinduced transition dynamics, showing a detailed picture of the light interaction with the electronic and lattice subsystems. This integrated approach significantly enhances the understanding of V 4 O 7 complex behavior upon photoexcitation, opening new possibilities for developing new optoelectronic devices and noninvasive optical control of the phase transition pathways in vanadates.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Development of Numerical Model of Metal Foam with PCM for the Estimation of Effective Thermal Conductivity

Global warming due to climate change is a threat to humankind. Nuclear energy is one of the promising solutions to reduce fossil fuel usage. Nuclear energy can handle the base load, compensating for the volatility of renewable energy. If nuclear energy could achieve load following capability, the combination with renewable energy would be more suitable. Thermal energy storage (TES) is one of the options for enabling load following of nuclear reactors. The TES makes it possible to store surplus nuclear thermal energy and release it later as needed. In Idaho National Laboratory (INL), a new concept of latent heat TES integrated with high-temperature heat pipe has been proposed and is under development, which is called Heat pipe-Integrated Thermal Battery (HITB). HITB exchanges thermal energy between the reactor system and TES via heat pipe. The heat transferred to TES medium, made of phase change material (PCM), stores energy as sensible heat and/or latent heat. As PCM typically has poor thermal conductivity, however, various heat transfer enhancement techniques are required to achieve a rapid charging cycle. There are many techniques to enhance the heat transfer ability of TES medium such as disk, fin, and metal foam. Among them, metal foam is an appropriate option to enhance the heat transfer because it maximizes the heat transfer area through metal wicks. Metal foam is a lightweight metal structure that has a high porosity of over 0.9. The typical materials for metal foam are Aluminum, Copper, Nickel, and Silicon Carbide (SiC). Metal foam not only enhances heat transfer via conduction but also increases contact surface area. In the HITB design , the metal foam is being considered as one of the options to enhance the heat transfer of TES medium (PCM) [1]. To predict the enhanced thermal performance of TES, one should properly estimate the effective thermal conductivity of metal foam combined with PCM material or calculate heat transfer in distributed model. There are many experimental works that provides effective thermal conductivity of metal foam with various PCM [2,3]. Also, many theoretical models were developed based on the unit cell model of metal foam [4,5]. With a distributed model, on the other hand, detail heat transfer characteristics between metal foam and PCM material can be analyzed considering the geometry or buoyancy effect. However, due to the complex geometry of metal foam pores, the computational cost for three-dimensional modeling highly increases. Therefore, if metal foam structure can be modeled in simple and repetitive design, the computational cost would decrease Among the various metal foam models [2], lattice model is one of the simple and extendable design. The porosity and pores per inch (PPI) can be characterized by the size and spatial distance of lattice structure. If the three-dimensional metal foam model consists of lattice structure could properly estimate the heat transfer, which is characterized by effective thermal conductivity, it would be a good option to assess the thermal performance of metal foam with PCM. In this study, a three-dimensional numerical model was developed to simulate conductive heat transfer between metal foam and PCM. The three-dimensional lattice structure of square pillars was selected as a basic structure of the metal foam. The calculation result was characterized by the effective thermal conductivity of the whole domain. A sensitivity study was conducted for mesh size, domain size, and PPI to check whether the calculation result gives a converged result or not. Lastly, the effective thermal conductivity from the lattice model was compared with existing experimental data to validate the model result

25 ENERGY STORAGE↗