Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microcanonical ensemble”

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.

15 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The endpoint of partial deconfinement

We study the matrix quantum mechanics of two free hermitian N × N matrices subject to a singlet constraint in the microcanonical ensemble. This is the simplest example of a theory that at large N has a confinement/deconfinement transition. In the microcanonical ensemble, it also exhibits partial deconfinement with a Hagedorn density of states. We argue that the entropy of these configurations, based on a combinatorial counting of Young diagrams, are dominated by Young diagrams that have the VKLS shape. When the shape gets to the maximal depth allowed for a Young diagram of SU(N), namely N, we argue that the system stops exhibiting the Hagedorn behavior. The number of boxes (energy) at the transition is N 2 /4, independent of the charge of the state.

1/N Expansion↗
DOE OSTI2025

Parton distribution functions from scalar light-front parton gas model

Here, we propose an application of a microcanonical ensemble with light-front kinematics to model the phase-space distribution of relativistic constituents of a bound state. These constituents denoted by partons are treated as classical spin-zero particles confined inside the bound state with inter-parton collisions as their only interaction. The microcanonical molecular dynamics ensemble is applied to obtain the phase-space distribution of such a thermodynamic system. We sample this phase-space distribution using Monte Carlo algorithms to obtain the parton distribution functions (PDFs) in scenarios with 3, 4, and 5 identical partons. In addition PDFs when a selected number of massless partons are mixed with 3 massive partons are also presented.

Microcanonical ensemble↗
DOE OSTI2024

Stochastic symplectic reduced-order modeling for model-form uncertainty quantification in molecular dynamics simulations in various statistical ensembles

Here, this work focuses on the representation of model-form uncertainties in molecular dynamics simulations in various statistical ensembles. In prior contributions, the modeling of such uncertainties was formalized and applied to quantify the impact of, and the error generated by, pair-potential selection in the microcanonical ensemble (NVE). In this work, we extend this formulation and present a linear-subspace reduced-order model for the canonical (NVT) and isobaric (NPT) ensembles. The symplectic reduced-order basis is randomized on the tangent space of the Stiefel manifold to provide topological relationships and capture model-form uncertainty. Using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), we assess the relevance of these stochastic reduced-order atomistic models on canonical problems involving a Lennard-Jones fluid and an argon crystal melt.

42 ENGINEERING↗
DOE OSTI2025

Mind the crosscap: $τ$-scaling in non-orientable gravity and time-reversal-invariant systems

Spectral statistics of quantum chaotic systems are governed by random matrix universality. In many cases of interest, time-reversal symmetry selects the Gaussian Orthogonal Ensemble (GOE) as the relevant universality class. In holographic CFTs, this is mirrored by the presence of non-orientable geometries in the dual gravitational path integral. In this work, we analyze general properties of these matrix models and their gravitational counterparts. First, we develop a formalism to express the universal level statistics in the canonical ensemble for arbitrary spectral curves, leading to a topological expansion with finite radius of convergence in the late-time $τ$-scaling limit. Then, we focus on topological gravity and study topological recursion on the moduli space of non-orientable surfaces. We find that the Weil-Petersson volumes display non-analytic behaviour multiplying polynomials in the boundary lengths. The volumes give rise to wormholes with late-time divergences, in contrast with the orientable case, which is finite. We identify systematic cancellations among WP volumes implied by the consistency and finiteness of the $τ$-scaling limit. In particular, the cancellation of late-time divergences requires a nontrivial genus resummation. Working in the gravitational microcanonical ensemble, we derive and resum all orders of the topological expansion matching the GOE matrix model in the high-energy regime.

Chaotic Dynamics (nlin.CD)↗
DOE OSTI2021

SYK wormhole formation in real time

We study the real time formation of the ground state of two coupled SYK models. This is a highly entangled state which is close to the thermofield double state and can be viewed as a wormhole. We start from a high temperature state, we let it cool by coupling to a cold bath. We numerically solve for the large N dynamics. Our main result is that the system forms a wormhole by going through a region with negative specific heat, taking time that is independent of N. The dynamics is smooth everywhere and it seems to follow the equilibrium thermodynamic configurations of the microcanonical ensemble. Also we comment on the relation between this coupled SYK model and Jackiw-Teitelboim gravity theory with bulk fields.

1/N expansion↗
DOE OSTI2023

Large N algebras and generalized entropy

We construct a Type II ∞ von Neumann algebra that describes the large N physics of single-trace operators in AdS/CFT in the microcanonical ensemble, where there is no need to include perturbative 1/N corrections. Using only the extrapolate dictionary, we show that the entropy of semiclassical states on this algebra is holographically dual to the generalized entropy of the black hole bifurcation surface. From a boundary perspective, this constitutes a derivation of a special case of the QES prescription without any use of Euclidean gravity or replicas; from a purely bulk perspective, it is a derivation of the quantum-corrected Bekenstein-Hawking formula as the entropy of an explicit algebra in the G → 0 limit of Lorentzian effective field theory quantum gravity. In a limit where a black hole is first allowed to equilibrate and then is later potentially re-excited, we show that the generalized second law is a direct consequence of the monotonicity of the entropy of algebras under trace-preserving inclusions. Finally, by considering excitations that are separated by more than a scrambling time we construct a “free product” von Neumann algebra that describes the semiclassical physics of long wormholes supported by shocks. We compute Rényi entropies for this algebra and show that they are equal to a sum over saddles associated to quantum extremal surfaces in the wormhole. Surprisingly, however, the saddles associated to “bulge” quantum extremal surfaces contribute with a negative sign.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Large-scale molecular dynamics simulations of bubble collapse in water: Effects of system size, water model, and nitrogen

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 10 6 and 4.5 × 10 9 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L –1 . Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models’ differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models’ thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. Finally, we also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state.

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

Path integral factorization and the gravitational effective action

Abstract We discuss the factorization and continuity properties of fields in the Euclidean gravitational path integral with higher dimension operators constructed from powers of the Riemann tensor. We construct the boundary terms corresponding to the microcanonical ensemble and show that the saddle point approximation to the path integral with a quasilocal energy constraint generally yields a saddle point with discontinuous temperature. This extends a previous result for the Euclidean Schwarzschild–de Sitter geometry in Einstein gravity and shows that it is robust against at least some types of quantum corrections from heavy fields. As an application, we compute the entropy of SdS in D = 4 using the BTZ method. Our result matches the entropy calculated using Wald’s formula.

Astronomy & Astrophysics↗
DOE OSTI2024

Fluctuations in the entropy of Hawking radiation

We use the gravitational path integral (GPI) to compute the fluctuations of the Hawking radiation entropy around the Page curve in a two-dimensional model introduced by Penington Before the Page time, we find that δ S = e − S / 2 , where S is the black hole entropy. This result agrees with the Haar-averaged entropy fluctuations in a bipartite system. After the Page time, we find that δ S ∼ e − S , up to a prefactor that depends logarithmically on the width of the microcanonical energy window. This is not symmetric under exchange of subsystem sizes and so does not agree with the Haar average for a subsystem of fixed Hilbert space dimension. The discrepancy can be attributed to the fact that the black hole Hilbert space dimension is not fixed by the state preparation: even in a microcanonical ensemble with a top-hat smearing function, the GPI yields an additive fluctuation in the number of black hole states. This result, and the fact that the Page curve computed by the GPI is smooth, all point towards an ensemble interpretation of the GPI. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Neutrino many-body flavor evolution: The full Hamiltonian

We study neutrino flavor evolution in the quantum many-body approach using the full neutrino-neutrino Hamiltonian, including the usually neglected terms that mediate nonforward scattering processes. Working in the occupation number representation with plane waves as single-particle states, we explore the time evolution of simple initial states with up to N = 10 neutrinos. We discuss the time evolution of the Loschmidt echo, one body flavor and kinetic observables, and the one-body entanglement entropy. For the small systems considered, we observe “thermalization” of both flavor and momentum degrees of freedom on comparable time scales, with results converging towards expectation values computed within a microcanonical ensemble. We also observe that the inclusion of nonforward processes generates a faster flavor evolution compared to the one induced by the truncated (forward) Hamiltonian. Published by the American Physical Society 2024

Cirigliano, Vincenzo (ORCID:000000029056754X)↗
DOE OSTI2024

Improved Precision Scaling for Simulating Coupled Quantum-Classical Dynamics

In this paper, we present a superpolynomial improvement in the precision scaling of quantum simulations for coupled quantum-classical systems. Such systems are found in, e.g., molecular-dynamics simulations within the Born-Oppenheimer approximation. By employing a framework based on the Koopman–von Neumann formulation of classical mechanics, we express the Liouville equation of motion as unitary dynamics and utilize phase kickback from a dynamical quantum simulation to calculate the quantum forces acting on classical particles. This approach allows us to simulate the dynamics of these classical particles without the overheads associated with measuring gradients and solving the equations of motion on a classical computer, resulting in a superpolynomial advantage at the price of increased space complexity. We demonstrate that these simulations can be performed in both microcanonical and canonical ensembles, enabling the estimation of thermodynamic properties from the prepared probability density. Published by the American Physical Society 2024

Simon, Sophia (ORCID:0000000332612338)↗
DOE OSTI2024

Quantifying Quantum Chaos through Microcanonical Distributions of Entanglement

A characteristic feature of “quantum chaotic” systems is that their eigenspectra and eigenstates display universal statistical properties described by random matrix theory (RMT). However, eigenstates of local systems also encode structure beyond RMT. To capture this feature, we introduce a framework that allows us to compare the properties of eigenstates in local systems with those of pure random states. In particular, our framework defines a notion of distance between quantum state ensembles that utilizes the Kullback-Leibler divergence to compare the microcanonical distribution of entanglement entropy (EE) of eigenstates with a reference RMT distribution generated by pure random states (with appropriate constraints). This notion gives rise to a quantitative metric for quantum chaos that not only accounts for averages of the distributions but also higher moments. The differences in moments are compared on a highly resolved scale set by the standard deviation of the RMT distribution, which is exponentially small in system size. As a result, the metric can distinguish between chaotic and integrable behaviors and, in addition, quantify and compare the of chaos (in terms of proximity to RMT behavior) between two systems that are assumed to be chaotic. We implement our framework in local, minimally structured, Floquet random circuits, as well as a canonical family of many-body Hamiltonians, the mixed-field Ising model (MFIM). Importantly, for Hamiltonian systems, we find that the reference random distribution must be appropriately constrained to incorporate the effect of energy conservation in order to describe the ensemble properties of midspectrum eigenstates. The metric captures deviations from RMT across all models and parameters, including those that have been previously identified as strongly chaotic, and for which other diagnostics of chaos such as level spacing statistics look strongly thermal. In Floquet circuits, the dominant source of deviations is the second moment of the distribution, and this persists for all system sizes. For the MFIM, we find significant variation of the KL divergence in parameter space. Notably, we find a small region where deviations from RMT are minimized, suggesting that “maximally chaotic” Hamiltonians may exist in fine-tuned pockets of parameter space. Published by the American Physical Society 2024

Physics↗
DOE OSTI2023

Wormholes from heavy operator statistics in AdS/CFT

We construct higher dimensional Euclidean AdS wormhole solutions that reproduce the statistical description of the correlation functions of an ensemble of heavy CFT operators. We consider an operator which effectively backreacts on the geometry in the form of a thin shell of dust particles. Assuming dynamical chaos in the form of the ETH ansatz, we demonstrate that the semiclassical path integral provides an effective statistical description of the microscopic features of the thin shell operator in the CFT. The Euclidean wormhole solutions provide microcanonical saddlepoint contributions to the cumulants of the correlation functions over the ensemble of operators. We finally elaborate on the role of these wormholes in the context of non-perturbative violations of bulk global symmetries in AdS/CFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Coupled fission fragment angular momenta

Here, nuclear fission produces fragments whose spins are coupled to the relative angular motion via angular momentum conservation. It is shown how ensembles of such spins can readily be obtained by either direct microcanonical sampling or by sampling of the associated normal modes of rotation. The resulting distribution of the spin-spin opening angle is illustrated in various three- and two-dimensional scenarios and it is demonstrated how recent mutually conflicting model calculations can be well reproduced with different assumptions about the scission geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Calculation of mode lifetimes in weakly anharmonic solids using self-consistent ensemble eigenstates of the Liouvillian

Recent schemes for the efficient calculation of vibrational mode frequencies and lifetimes in weakly anharmonic systems such as solid crystal structures have shown promise, with simple relationships being observed between low order moments of the Liouvillian and the observed lifetime. While only shown empirically, two parameters seemed sufficient to determine the lifetime to high accuracy over a large range of temperatures for a few simple low degree of freedom models and a Lennard–Jones solid. In this paper we analytically reproduce the numerically observed relationship using a few generally applicable assumptions about the density of states. This advance is made by studying the autocorrelation based on an “ensemble eigenstate” of the Liouvillian. We propose the use of this ensemble eigenstate as the basis for autocorrelations to simplify calculations in richer systems. Finally, an analysis is also performed of the relationship between canonical and microcanonical autocorrelation functions and a long standing error in the long time behavior of such functions originally presented by Sen, Sinkovits, and Chakravarti is corrected.

36 MATERIALS SCIENCE↗