SEARCH · Engineering Papers

Results for “statistical mechanics 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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Universal Kardar-Parisi-Zhang Dynamics in Integrable Quantum Systems

Although the Bethe ansatz solution of the spin-1/2 Heisenberg model dates back nearly a century, the anomalous nature of its high-temperature transport dynamics has only recently been uncovered. Indeed, numerical and experimental observations have demonstrated that spin transport in this paradigmatic model falls into the Kardar-Parisi-Zhang (KPZ) universality class. This has inspired the significantly stronger conjecture that KPZ dynamics, in fact, occur in all integrable spin chains with non-Abelian symmetry. Here, in this letter, we provide extensive numerical evidence affirming this conjecture. Moreover, we observe that KPZ transport is even more generic, arising in both supersymmetric and periodically driven models. Motivated by recent advances in the realization of SU(N)-symmetric spin models in alkaline-earth-based optical lattice experiments, we propose and analyze a protocol to directly investigate the KPZ scaling function in such systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Spin crossovers and superdiffusion in the one-dimensional Hubbard model

We use tools from integrability and generalized hydrodynamics to study finite-temperature dynamics in the one-dimensional Hubbard model. First, we examine charge, spin, and energy transport away from half-filling and zero magnetization, focusing on the strong coupling regime where we identify a rich interplay of temperature and energy scales, with crossovers between distinct dynamical regimes. We identify an intermediate-temperature regime analogous to the spin-incoherent Luttinger liquid, where spin degrees of freedom are hot but charge degrees of freedom are at low temperature. We demonstrate that the spin Drude weight exhibits sharp features at the crossover between this regime and the low-temperature Luttinger liquid regime, which are absent in the charge and energy response, and rationalize this behavior in terms of the properties of Bethe ansatz quasiparticles.We then turn to the dynamics along special lines in the phase diagram corresponding to half-filling and/or zero magnetization where on general grounds we anticipate that the transport is subballistic but superdiffusive. We provide analytical and numerical evidence for Kardar-Parisi-Zhang (KPZ) dynamical scaling (with length and time scales related via x ~ t 2/3 ) along both lines and at the SO(4)-symmetric point where they intersect. Furthermore, our results suggest that both spin-coherence crossovers and KPZ scaling may be accessed in near-term experiments with optical lattice Hubbard emulators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Modeling single-molecule stretching experiments using statistical thermodynamics

Single-molecule stretching experiments are widely utilized within the fields of physics and chemistry to characterize the mechanics of individual bonds or molecules, as well as chemical reactions. Analytic relations describing these experiments are valuable, and these relations can be obtained through the statistical thermodynamics of idealized model systems representing the experiments. Since the specific thermodynamic ensembles manifested by the experiments affect the outcome, primarily for small molecules, the stretching device must be included in the idealized model system. Though the model for the stretched molecule might be exactly solvable, including the device in the model often prevents analytic solutions. In the limit of large or small device stiffness, the isometric or isotensional ensembles can provide effective approximations, but the device effects are missing. Here a dual set of asymptotically correct statistical thermodynamic theories are applied to develop accurate approximations for the full model system that includes both the molecule and the device. In conclusion, the asymptotic theories are first demonstrated to be accurate using the freely jointed chain model and then using molecular dynamics calculations of a single polyethylene chain.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Natural Fracture Pattern Development Using 3D Coupled Numerical Model: Mechanics and Diagenesis

In the earth's crust, while tectonic processes can mechanically create natural fractures, cementation can simultaneously grow on the surfaces of the fractures. This crystallization can affect the growth of the fractures. To more accurately address the modeling of natural fracture development, both mechanical deformation and chemical diagenesis should be carried out in an integrated manner. Here, we achieve this with a three-dimensional (3D) simulator that couples mechanical and chemical processes for natural fracture pattern development in layer confined situations. Tectonic loading over geologic time is mimicked by displacement driven boundary conditions in a 3D displacement discontinuity method (DDM) scheme. Quartz diagenesis competes with mechanical fracture opening, resulting in three different cementation conditions - open, bridging or fully infilled – which define fracture element stiffness that is used as feedback within the 3D-DDM code during propagation. An important milestone of this work is that this coupled structural diagenetic model can reproduce power-law aperture population statistics which cannot be generated using mechanics alone. More accurately estimated mechanical fracture apertures and their cementation condition is a critical step toward better permeability estimation and fluid path analysis for subsurface reservoirs.

58 GEOSCIENCES↗

DOE OSTI2024

A spatially-resolved model of neutron-irradiated tungsten coupling stochastic cluster dynamics and finite deformation plasticity

Structural materials used in nuclear reactors face severe degradation in mechanical properties, such as hardening and embrittlement. At the microscopic scale, this occurs due to creation and accumulation of irradiation-induced defects and their interaction with system dislocations. Although techniques exist which can model evolution of irradiation defects, for instance kinetic transport theory-based models, their interaction with mechanical deformation of the bulk material has not been investigated extensively. In this work, we demonstrate a novel spatially-resolved multiscale coupling between microscopic irradiation defect evolution, modeled using Stochastic Cluster Dynamics (SCD) and macroscopic mechanical deformation modeled using a finite-deformation plasticity model. SCD is used to determine the statistically averaged defect cluster spacing, dependent on operating conditions such as irradiation dose and temperature. This acts as an initial condition that governs the critical resolved shear stress of dislocation glide in the macroscopic plasticity model. This framework is used to predict mechanical behavior in post-mortem test of irradiated Tungsten samples, which has found its importance as structural material used in nuclear reactors. The results obtained using the coupled approach are in good agreement with experimental data of uniaxial tension tests. The model is able to capture the effect of temperature and irradiation dose on the material hardening. Two methods are proposed to estimate hardness – using Tabor's Law relating uniaxial yield stress to hardness and from flat-punch simulations. The results are in reasonable agreement with hardness data from micro-indentation experiments of irradiated Tungsten samples. Finally, the model is also able to reveal microstructural details such as spatial variation in defect density and local stress.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Size-Dependent Nucleation in Crystal Phase Transition from Machine Learning Metadynamics

In this Letter, we present a framework that combines machine learning potential (MLP) and metadynamics to investigate solid-solid phase transition. Here, based on the spectral descriptors and neural networks regression, we develop a scalable MLP model to warrant an accurate interpolation of the energy surface where two phases coexist. Applying it to the simulation of B4–B1 phase transition of GaN under 50 GPa with different model sizes, we observe sequential change of the phase transition mechanism from collective modes to nucleation and growths. When the size is at or below 128 000 atoms, the nucleation and growth appear to follow a preferred direction. At larger sizes, the nuclei occur at multiple sites simultaneously and grow to microstructures by passing the critical size. The observed change of the atomistic mechanism manifests the importance of statistical sampling with large system size in phase transition modeling.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Linear growth of circuit complexity from Brownian dynamics

How rapidly can a many-body quantum system generate randomness? Using path integral methods, we demonstrate that Brownian quantum systems have circuit complexity that grows linearly with time. In particular, we study Brownian clusters of N spins or fermions with time-dependent all-to-all interactions, and calculate the Frame Potential to characterize complexity growth in these models. In both cases the problem can be mapped to an effective statistical mechanics problem which we study using path integral methods. Within this framework it is straightforward to show that the kth Frame Potential comes within ϵ of the Haar value after a time of order t ~ kN + k log k + log ϵ –1 . Using a bound on the diamond norm, this implies that such circuits are capable of coming very close to a unitary k-design after a time of order t ~ kN. We also consider the same question for systems with a time-independent Hamiltonian and argue that a small amount of time-dependent randomness is sufficient to generate a k-design in linear time provided the underlying Hamiltonian is quantum chaotic. These models provide explicit examples of linear complexity growth that are analytically tractable and are directly applicable to practical applications calling for unitary k-designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Statistical Mechanics on Lattices (smol) v0.0.1

smol is a Python implementation of the Cluster Expansion Method and extensions of this methodology; which is used to fit applied lattice models from first principle calculations for subsequent use in Monte Carlo simulations for computing thermodynamic and statistical properties of atomic configuration. smol has a several notable advantages over other similar software packages. The first being its Python implementation which makes it easy to use even for users with little software and coding skills. Despite being implemented in Python, essential Monte Carlo routines are implemented in Cython such that performance is not compromised and is competitive even to similar software implemented in C/C++. Furthermore, the package has specific functionality for handling complex ionic materials (such as cathodes and electrolytes) that is not readily available in other packages. Lastly, the software has a flexible and modular design with the intention of making it fast and efficient to develop and extend the methodology.

Barroso-Luque, Luis↗

DOE OSTI2024

Degradation of a Cr-Mo steel by carbide precipitation over long-term service

Microstructure of the 40-year service-exposed 1Cr-0.5Mo steel extracted from pressure vessels through the boat sampling technique was characterized using x-ray diffraction, field-emission scanning electron microscopy, energy dispersive spectrometry, electron backscatter diffraction, focused ion beam, and transmission electron microscopy. The microstructure degradation during long-term service was identified. Statistical crystal plastic ity models were developed to correlate the microstructure and mechanical properties quantitatively. The steel is found significantly strengthened but slightly embrittled. Nanoscale Mo 2 C carbides have nucleated within the ferrite matrix over the long-term service, resulting in significant strengthening of the matrix. The coarser grain boundary Cr 23 C 6 carbides, on the other hand, act as the damage accumulation sites and lead to the degradation.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Progress toward hydro-equivalent ignition in OMEGA direct-drive DT-layered implosions

Considerable progress has been made in deuterium-tritium-layered implosion experiments on the OMEGA Laser System, bringing the prospects for thermonuclear ignition in direct-drive configurations with megajoule-class lasers closer to reality. Doing so has required navigating the balance between improved 1D performance and multidimensional stability. Using statistical modeling based on over 350 cryogenic implosions to identify various degradation mechanisms, and combined with multidimensional simulations and experimental techniques such as target offsets to combat residual flows, core conditions have repeatably been achieved that extrapolate to the burning-plasma state when scaled to 2.15 MJ of symmetric laser illumination. Using high implosion velocities (⁠> 450 km/s) and moderately high adiabats (⁠~5⁠), these experiments produced record-high scaled Lawson parameters in direct drive equal to 89 ± 2% of that required for ignition with expected yields of up to 1.5 ± 0.2 MJ. To improve these results still further, focused physics studies are performed to improve physics understanding and identify routes to even greater performance. Recent studies include investigations into the impact of mounting features, laser imprint, reduced fuel temperatures, and greater on-target intensities through subscale experiments. This manuscript gives a summary of the cryogenic direct-drive program on the OMEGA laser, including routes taken to achieve the current best performance, the status of recent focused physics investigations, and future designs—such as target solutions to laser imprint and reducing vapor density to increase convergence—that are expected lead to the demonstration of hydro-equivalent ignition on OMEGA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Phase diagram of generalized XY model using the tensor renormalization group

We use the higher-order tensor renormalization group method to study the two-dimensional generalized XY model that admits integer and half-integer vortices. This model is the deformation of the classical XY model and has a rich phase structure consisting of nematic, ferromagnetic, and disordered phases and three transition lines belonging to the Berezinskii-Kosterlitz-Thouless and Ising class. We explore the model for a wide range of temperatures, T , and the deformation parameter, Δ , and compute specific heat along with integer and half-integer magnetic susceptibility, finding both Berezinskii-Kosterlitz-Thouless-like and Ising-like transitions and the region where they meet. Published by the American Physical Society 2024

2-dimensional systems↗

DOE OSTI2022

Observation of a multiplicity dependence in the $\mathcal{p}_T$-differential charm baryon-to-meson ratios in proton–proton collisions at $\sqrt{s}$ = 13 TeV

The production of prompt D 0 , D$^+_s$, and Λ$^+_c$ hadrons, and their ratios, D$^+_s$/D 0 and Λ$^+_c$/D 0 , are measured in proton–proton collisions at $\sqrt{s}$ = 13 TeV at midrapidity (|y| < 0.5) with the ALICE detector at the LHC. The measurements are performed as a function of the charm-hadron transverse momentum (p T ) in intervals of charged-particle multiplicity, measured with two multiplicity estimators covering different pseudorapidity regions. While the strange to non-strange D$^+_s$/D 0 ratio indicates no significant multiplicity dependence, the baryon-to-meson p T -differential Λ$^+_c$/D 0 ratio shows a multiplicity-dependent enhancement, with a significance of 5.3σ for 1 < p T < 12 GeV/c, comparing the highest multiplicity interval with respect to the lowest one. The measurements are compared with a theoretical model that explains the multiplicity dependence by a canonical treatment of quantum charges in the statistical hadronisation approach, and with predictions from event generators that implement colour reconnection mechanisms beyond the leading colour approximation to model the hadronisation process. The Λ$^+_c$/D 0 ratios as a function of p T present a similar shape and magnitude as the Λ/K$^0_s$ ratios in comparable multiplicity intervals, suggesting a potential common mechanism for light- and charm-hadron formation, with analogous multiplicity dependence. The p T -integrated ratios, extrapolated down to p T = 0, do not show a significant dependence on multiplicity within the uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

A statistical mechanics approach to macroscopic limits of car-following traffic dynamics

Herein we study the derivation of macroscopic traffic models from car-following vehicle dynamics by means of hydrodynamic limits of an Enskog-type kinetic description. We consider the superposition of Follow-the-Leader (FTL) interactions and relaxation towards a traffic-dependent Optimal Velocity (OV) and we show that the resulting macroscopic models depend on the relative frequency between these two microscopic processes. If FTL interactions dominate then one gets an inhomogeneous Aw–Rascle–Zhang model, whose (pseudo) pressure and stability of the uniform flow are precisely defined by some features of the microscopic FTL and OV dynamics. Conversely, if the rate of OV relaxation is comparable to that of FTL interactions then one gets a Lighthill–Whitham–Richards model ruled only by the OV function. We further confirm these findings by means of numerical simulations of the particle system and the macroscopic models. Unlike other formally analogous results, our approach builds the macroscopic models as physical limits of particle dynamics rather than assessing the convergence of microscopic to macroscopic solutions under suitable numerical discretisations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Probabilistic Failure Criterion of SiC/SiC Composites Under Multiaxial Loading

Owing to its excellent mechanical properties and stability under high temperature and neutron irradiation conditions, SiC/SiC composites have emerged as a promising material for light water reactors (LWRs) in the development of accident-tolerant fuel (ATF) systems. Structural integrity and retention of hermeticity are two crucial requirements for SiC/SiC claddings during normal operations, and both of them are closely related to the proportional limit stress (PLS) of the material. Understanding the behavior of SiC/SiC composites under multiaxial stress states and developing a probabilistic approach for evaluating the structural vulnerability are of paramount importance for reliability-based analysis and design of SiC/SiC composite claddings. So far, there has been very limited effort towards experimental and analytical investigations of probabilistic failure of SiC/SiC claddings. This critical knowledge gap motivates this research. A probabilistic failure criterion for SiC/SiC composites under multi-axial loading is developed, and this criterion is incorporated into reliability analysis of the structural integrity of SiC/SiC fuel cladding. The research consists of two parts: 1) experimental investigation of multiaxial failure behavior of SiC/SiC composites, and 2) theoretical modeling of time-dependent probabilistic failure of SiC/SiC cladding. In the experimental investigation, the PLS is determined through the examination of stress-strain response, the acoustic emission measurement, as well as the X-ray computed tomography. The theoretical framework is derived by combin- ing the finite weakest-link statistical model and the subcritical damage growth model. This theoretical model captures the time-dependent failure mechanism of the material, which has a major consequence for predicting the lifetime distribution of the cladding. Meanwhile, the model also predicts that the failure statistics of the cladding depends strongly on the cladding length. The results of the multiaxial experiments reveal the level of statistical variation of the PLS of SiC/SiC materials under different stress states. The theoretical model provides a robust analytical tool for extrapolation of small-scale laboratory test results to the behavior of full-scale claddings. These findings establish a scientific foundation for the development of reliability-based design of SiC/SiC fuel claddings, which will play an essential role in improving the structural safety and integrity of LWRs.

42 ENGINEERING↗

DOE OSTI2024

Dynamics of magnetization at infinite temperature in a Heisenberg spin chain

Understanding universal aspects of quantum dynamics is an unresolved problem in statistical mechanics. In particular, the spin dynamics of the one-dimensional Heisenberg model were conjectured as to belong to the Kardar-Parisi-Zhang (KPZ) universality class based on the scaling of the infinite-temperature spin-spin correlation function. In a chain of 46 superconducting qubits, we studied the probability distribution of the magnetization transferred across the chain’s center, P M . The first two moments of P M show superdiffusive behavior, a hallmark of KPZ universality. However, the third and fourth moments ruled out the KPZ conjecture and allow for evaluating other theories. Our results highlight the importance of studying higher moments in determining dynamic universality classes and provide insights into universal behavior in quantum systems.

Science & Technology - Other Topics↗

DOE OSTI2021

Error Bounds for Dynamical Spectral Estimation

Dynamical spectral estimation is a well-established numerical approach for estimating eigenvalues and eigenfunctions of the Markov transition operator from trajectory data. Although the approach has been widely applied in biomolecular simulations, its error properties remain poorly understood. Here we analyze the error of a dynamical spectral estimation method called “the variational approach to conformational dynamics" (VAC). We bound the approximation error and estimation error for VAC estimates. Our analysis establishes VAC's convergence properties and suggests new strategies for tuning VAC to improve accuracy.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS↗

DOE OSTI2020

Enantioselective Allenoate-Claisen Rearrangement Using Chiral Phosphate Catalysts

Herein we report the first highly enantioselective allenoate-Claisen rearrangement using doubly axially chiral phosphate sodium salts as catalysts. This synthetic method provides access to β-amino acid derivatives with vicinal stereocenters in up to 95% ee. We also investigated the mechanism of enantioinduction by transition state (TS) computations with DFT as well as statistical modeling of the relationship between selectivity and the molecular features of both the catalyst and substrate. The mutual interactions of charge-separated regions in both the zwitterionic intermediate generated by reaction of an amine to the allenoate and the Na + -salt of the chiral phosphate leads to an orientation of the TS in the catalytic pocket that maximizes favorable noncovalent interactions. Crucial arene-arene interactions at the periphery of the catalyst lead to a differentiation of the TS diastereomers. Finally, these interactions were interrogated using DFT calculations and validated through statistical modeling of parameters describing noncovalent interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗