Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large deviations”

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 37 records · Page 2

Physics-tailored machine learning reveals unexpected physics in dusty plasmas

Dusty plasma is a mixture of ions, electrons, and macroscopic charged particles that is commonly found in space and planetary environments. The particles interact through Coulomb forces mediated by the surrounding plasma, and as a result, the effective forces between particles can be nonconservative and nonreciprocal. Machine learning (ML) models are a promising route to learn these complex forces, yet their structure should match the underlying physical constraints to provide useful insight. Here, we demonstrate and experimentally validate an ML approach that incorporates physical intuition to infer force laws in a laboratory dusty plasma. Trained on 3D particle trajectories, the model accounts for inherent symmetries, nonidentical particles, and learns the effective nonreciprocal forces between particles with exquisite accuracy (R 2 > 0.99). We validate the model by inferring particle masses in two independent yet consistent ways. The model’s accuracy enables precise measurements of particle charge and screening length, identifying large deviations from common theoretical assumptions. Our ability to identify unknown physics from experimental data demonstrates how ML-powered approaches can guide new routes of scientific discovery in many-body systems. Furthermore, we anticipate our ML approach to be a starting point for inferring laws from dynamics in a wide range of many-body systems, from colloids to living organisms.

Science & Technology - Other Topics↗

Uniform modelling of the stellar density of thirteen tidal streams within the Galactic halo

ABSTRACT We present the results of fitting a flexible stellar stream density model to a collection of thirteen streams around the Milky Way, using photometric data from DES, DECaLS, and Pan-STARRS. We construct density maps for each stream and characterize their tracks on the sky, width, and distance modulus curves along the length of each stream. We use these measurements to compute lengths and total luminosities of streams and identify substructures. Several streams show prominent substructures, such as stream broadening, gaps, large deviations of stream tracks, and sharp changes in stream densities. Examining the group of streams as a population, as expected we find that streams with globular cluster progenitors are typically narrower than those with dwarf galaxy progenitors, with streams around 100 pc wide showing overlap between the two populations. We also note the average luminosity of globular cluster streams is significantly lower than the typical luminosity of intact globular clusters. The likely explanation is that observed globular cluster streams preferentially come from lower luminosity and lower density clusters. The stream measurements done in a uniform manner presented here will be helpful for more detailed stream studies such as identifying candidate stream members for spectroscopic follow up and stellar stream dynamical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗

Scaling functions for eigenstate entanglement crossovers in harmonic lattices

For quantum matter, eigenstate entanglement entropies obey an area law or log-area law at low energies and small subsystem sizes and cross over to volume laws for high energies and large subsystems. This transition is captured by crossover functions, which assume a universal scaling form in quantum critical regimes. We demonstrate this for the harmonic lattice model, which describes quantized lattice vibrations and is a regularization for free scalar field theories, modeling, e.g., spin-0 bosonic particles. In one dimension, the ground-state entanglement obeys a log-area law. For dimensions d ≥ 2, it displays area laws, even at criticality. The distribution of excited-state entanglement entropies is found to be sharply peaked around subsystem entropies of corresponding thermodynamic ensembles in accordance with the eigenstate thermalization hypothesis. Numerically, we determine crossover scaling functions for the quantum critical regime of the model and do a large-deviation analysis. We show how infrared singularities of the system can be handled and how to access the thermodynamic limit using a perturbative trick for the covariance matrix. Eigenstates for quasi-free bosonic systems are not Gaussian. Here, we resolve this problem by considering appropriate squeezed states instead. For these, entanglement entropies can be evaluated efficiently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Distinct universality classes of diffusive transport from full counting statistics

The hydrodynamic transport of local conserved densities furnishes an effective coarse-grained description of the dynamics of a many-body quantum system. However, the full quantum dynamics contains much more structure beyond the simplified hydrodynamic description. Here we show that systems with the same hydrodynamics can nevertheless belong to distinct dynamical universality classes, as revealed by new classes of experimental observables accessible in synthetic quantum systems, which can, for instance, measure simultaneous site-resolved snapshots of all of the particles in a system. Specifically, we study the full counting statistics of spin transport, whose first moment is related to linear-response transport, but the higher moments go beyond. We present an analytic theory of the full counting statistics of spin transport in various integrable and nonintegrable anisotropic one-dimensional spin models, including the XXZ spin chain. We find that spin transport, while diffusive on average, is governed by a distinct non-Gaussian dynamical universality class in the models considered. We consider a setup in which the left and right half of the chain are initially created at different magnetization densities, and consider the probability distribution of the magnetization transferred between the two half-chains. We derive a closed-form expression for the probability distribution of the magnetization transfer, in terms of random walks on the half-line. We show that this distribution strongly violates the large-deviation form expected for diffusive chaotic systems, and explain the physical origin of this violation. Here, we discuss the crossovers that occur as the initial state is brought closer to global equilibrium. Our predictions can directly be tested in experiments using quantum gas microscopes or superconducting qubit arrays.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Valence, charge transfer, and orbital-dependent correlation in bilayer nickelates Nd 3 Ni 2 O 7

We examine the bulk electronic structure of Nd 3 Ni 2 O 7 using Ni 2 p core-level hard x-ray photoemission spectroscopy combined with density functional theory + dynamical mean-field theory. Our results reveal a large deviation of the Ni 3 d occupation from the formal Ni 2.5 + valency, highlighting the importance of the charge transfer from oxygen ligands. We find that the dominant d 8 configuration is accompanied by nearly equal contributions from d 7 and d 9 states, exhibiting an unusual valence state among Ni-based oxides. Finally, we discuss the Ni d x 2 − y 2 and d z 2 orbital-dependent hybridization, correlation and local spin dynamics. Published by the American Physical Society 2025

Takegami, Daisuke (ORCID:0000000338321849)↗

Quasiparton distributions in massive QED2: Toward quantum computation

We analyze the quasiparton distributions of the lightest 𝜂′ meson in massive two-dimensional quantum electrodynamics (QED2) by exact diagonalization. The Hamiltonian and boost operators are mapped onto spin qubits in a spatial lattice with open boundary conditions. The lowest excited state in the exact diagonalization is shown to interpolate continuously between an anomalous 𝜂′ state at strong coupling, and a nonanomalous heavy meson at weak coupling, with a cusp at the critical point. The boosted 𝜂′ state follows relativistic kinematics but with large deviations in the luminal limit. The spatial quasiparton distribution function and amplitude for the 𝜂′ state are computed numerically for increasing rapidity both at strong and weak coupling, and compared to the exact light front results. The numerical results from the boosted form of the spatial parton distributions, compare fairly with the inverse Fourier transformation of the luminal parton distributions, derived in the lowest Fock space approximation. Our analysis points out some of the limitations facing the current lattice program for the parton distributions.

Lattice field theory↗

Thermodynamic Control of Activity Patterns in Cytoskeletal Networks

Biological materials, such as the actin cytoskeleton, exhibit remarkable structural adaptability to various external stimuli by consuming different amounts of energy. In this Letter, we use methods from large deviation theory to identify a thermodynamic control principle for structural transitions in a model cytoskeletal network. Specifically, we demonstrate that biasing the dynamics with respect to the work done by nonequilibrium components effectively renormalizes the interaction strength between such components, which can eventually result in a morphological transition. Further, our work demonstrates how a thermodynamic quantity can be used to renormalize effective interactions, which in turn can tune structure in a predictable manner, suggesting a thermodynamic principle for the control of cytoskeletal structure and dynamics.

59 BASIC BIOLOGICAL SCIENCES↗

Resilience of the surface code to error bursts

Quantum error correction works effectively only if the error rate of gate operations is sufficiently low. However, some rare physical mechanisms can cause a temporary increase in the error rate that affects many qubits; examples include ionizing radiation in superconducting hardware and large deviations in the global control of atomic systems. We refer to such rare transient spikes in the gate error rate as error bursts. In this work, we investigate the resilience of the rotated surface code to generic error bursts. We assume that, after appropriate mitigation strategies, the spike in the error rate lasts for only a single syndrome-extraction cycle; we also assume that the enhanced error rate is uniform across the code block. Under these assumptions, and for a circuit-level depolarizing noise model, we perform Monte Carlo simulations to determine the regime in burst error rate and background error rate for which the memory time becomes arbitrarily long as the code block size grows. Our results indicate that suitable hardware mitigation methods combined with standard decoding methods may suffice to protect against transient error bursts in the rotated surface code.

Quantum error correction↗

Variational Path Sampling of Rare Dynamical Events

This article reviews the concepts and methods of variational path sampling. These methods allow computational studies of rare events in systems driven arbitrarily far from equilibrium. Based upon a statistical mechanics of trajectory space and leveraging the theory of large deviations, they provide a perspective from which dynamical phenomena can be studied with the same types of ensemble reweighting ideas that have been used for static equilibrium properties. Applications to chemical, material, and biophysical systems are highlighted.

Singh, Aditya N↗

Isothermal Microcalorimetric Analysis of Li/CF x Battery Discharge

Lithium/carbon monofluoride (Li/CF x ) batteries generate substantial amounts of heat during discharge in part due to the large deviation of the loaded voltage from the measured and thermodynamically predicted open circuit voltage. Here, in this study, we further analyze this system by estimating the equilibrium voltage (V eq ) and its temperature dependence (dV eq /dT) over the entire discharge range. Based on these results, the ohmic and entropic heat contributions to the overall heat flow are calculated from experimental data. The ohmic heat flow is consistent with an electrochemical impedance spectroscopy model and electrode density measurements, indicating decreasing porosity of the electrode during discharge. Entropic heating increases during discharge as the cell reaction becomes increasingly entropically unfavorable. The total energy dissipated by the cell (electrical and thermal) remains similar over the entire discharge, evidence for a two-phase reaction with a constant, and rate-independent, ΔG rxn of about -465 kJ mol -1 . However, near the end of discharge (>90% DOD), ΔG rxn changes significantly indicating possible secondary reactions.

25 ENERGY STORAGE↗

Recent Updates in ETOE-2 and MC 2 -3

The MC 2 -3 code is a multigroup cross section generation code for fast reactor applications, developed by Argonne National Laboratory (ANL) under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, first, the cross section libraries for ENDF/B-VII.1 and ENDF/B-VIII.0 data were generated using the updated ETOE-2 code. Preliminary verification and validation tests of the ENDF/B-VII.1 MC 2 -3 library were performed with the selected fast reactor benchmark problems. Processing and verifying the libraries suggested that significant efforts would be required for thoroughly verifying the ENDF/B-VII.1 MC 2 -3 library and successfully processing the ENDF/B-VIII.0 MC 2 -3 library. Secondly, the cross section generation capability of MC 2 -3 was updated with the intermediate group lattice calculation and the equivalent Dancoff-factor cell (EDC) method to significantly improve the performance of a twodimensional assembly calculation using the method of characteristics (MOC). This effort was made useful in implementing and verifying the EDC method in the Griffin cross section API. Finally, we analyzed a SFR problem, for which unusually large deviations in core eigenvalues from Monte Carlo solutions were reported, by generating multigroup cross sections with MC 2 -3, performing core calculations with DIF3D-VARIANT, and analyzing cross section and eigenvalue results against Monte Carlo solutions. The analysis confirmed that the MC 2 -3/DIF3D solutions were in good agreement with Monte Carlo solutions, providing the correct process of accurately generating broad-group cross sections with MC 2 -3

97 MATHEMATICS AND COMPUTING↗

Evaluation of Adaptive Volt-VAR to Mitigate PV Impacts [Slides]

Distributed generation (DG) sources like photovoltaic (PV) systems with advanced inverters are able to perform grid-support functions, like autonomous Volt-VAR that attempts to mitigate voltage issues by injecting or consuming reactive power. However, the Volt-VAR function operates with VAR priority, meaning real power may be curtailed to provide additional reactive power support. Since some locations on the grid may be more prone to higher voltages than others, PV systems installed at those locations may be forced to curtail more power, adversely impacting the value of that PV system. Adaptive Volt-VAR (AVV) could be implemented as an alternative, whereby the Volt-VAR reference voltage changes over time, but this functionality has not been well-explored in the literature. In this work, the potential benefits and grid impacts of AVV were investigated using yearlong quasi-static time-series (QSTS) simulations. After testing a variety of allowable AVV settings, we found that even with aggressive settings AVV resulted in <0.01% real power curtailment and significantly reduced the reactive power support required from the PV inverter compared to conventional Volt-VAR but did not provide much mitigation for extreme voltage conditions. The reactive power support provided by AVV was injected to oppose large deviations in voltage (in either direction), indicating that it could be useful for other applications like reducing voltage flicker or minimizing interactions with other voltage regulating devices.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Microseismic Monitoring at the Farnsworth CO 2 -EOR Field

The Farnsworth Unit in northern Texas is a field site for studying geologic carbon storage during enhanced oil recovery (EOR) using CO 2 . Microseismic monitoring is essential for risk assessment by detecting fluid leakage and fractures. We analyzed borehole microseismic data acquired during CO 2 injection and migration, including data denoising, event detection, event location, magnitude estimation, moment tensor inversion, and stress field inversion. We detected and located two shallow clusters, which occurred during increasing injection pressure. The two shallow clusters were also featured by large b values and tensile cracking moment tensors that are obtained based on a newly developed moment tensor inversion method using single-borehole data. The inverted stress fields at the two clusters showed large deviations from the regional stress field. The results provide evidence for microseismic responses to CO 2 /fluid injection and migration.

02 PETROLEUM↗

KMT-2019-BLG-1715: Planetary Microlensing Event with Three Lens Masses and Two Source Stars

We investigate the gravitational microlensing event KMT-2019-BLG-1715, the light curve of which shows two short-term anomalies from a caustic-crossing binary-lensing light curve: one with a large deviation and the other with a small deviation. We identify five pairs of solutions, in which the anomalies are explained by adding an extra lens or source component in addition to the base binary-lens model. We resolve the degeneracies by applying a method in which the measured flux ratio between the first and second source stars is compared with the flux ratio deduced from the ratio of the source radii. Applying this method leaves a single pair of viable solutions, in both of which the major anomaly is generated by a planetary-mass third body of the lens, and the minor anomaly is generated by a faint second source. A Bayesian analysis indicates that the lens comprises three masses: a planet-mass object with ∼2.6 M {sub J} and binary stars of K and M dwarfs lying in the galactic disk. We point out the possibility that the lens is the blend, and this can be verified by conducting high-resolution follow-up imaging for the resolution of the lens from the source.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating wind profiles in a numerical weather prediction model with Doppler lidar

We use Doppler lidar wind profiles from six locations around the globe to evaluate the wind profile forecasts in the boundary layer generated by the operational global Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts (ECMWF). The six locations selected cover a variety of surfaces with different characteristics (rural, marine, mountainous urban, and coastal urban). We first validated the Doppler lidar observations at four locations by comparison with co-located radiosonde profiles to ensure that the Doppler lidar observations were of sufficient quality. The two observation types agree well, with the mean absolute error (MAE) in wind speed almost always less than 1 m s –1 . Large deviations in the wind direction were usually only seen for low wind speeds and are due to the wind direction uncertainty increasing rapidly as the wind speed tends to zero. Time–height composites of the wind evaluation with 1 h resolution were generated, and evaluation of the model winds showed that the IFS model performs best over marine and coastal locations, where the mean absolute wind vector error was usually less than 3 m s –1 at all heights within the boundary layer. Larger errors were seen in locations where the surface was more complex, especially in the wind direction. For example, in Granada, which is near a high mountain range, the IFS model failed to capture a commonly occurring mountain breeze, which is highly dependent on the sub-grid-size terrain features that are not resolved by the model. The uncertainty in the wind forecasts increased with forecast lead time, but no increase in the bias was seen. At one location, we conditionally performed the wind evaluation based on the presence or absence of a low-level jet diagnosed from the Doppler lidar observations. The model was able to reproduce the presence of the low-level jet, but the wind speed maximum was about 2 m s –1 lower than observed. This is attributed to the effective vertical resolution of the model being too coarse to create the strong gradients in wind speed observed. Our results show that Doppler lidar is a suitable instrument for evaluating the boundary layer wind profiles in atmospheric models.

54 ENVIRONMENTAL SCIENCES↗

Finite-time correlations boost large voltage angle fluctuations in electric power grids

Abstract Decarbonization in the energy sector has been accompanied by an increased penetration of new renewable energy sources in electric power systems. Such sources differ from traditional productions in that, first, they induce larger, undispatchable fluctuations in power generation and second, they lack inertia. Recent measurements have indeed reported long, non-Gaussian tails in the distribution of local voltage frequency data. Large frequency deviations may induce grid instabilities, leading in worst-case scenarios to cascading failures and large-scale blackouts. In this article, we investigate how correlated noise disturbances, characterized by the cumulants of their distribution, propagate through meshed, high-voltage power grids. For a single source of fluctuations, we show that long noise correlation times boost non-Gaussian voltage angle fluctuations so that they propagate similarly to Gaussian fluctuations over the entire network. However, they vanish faster, over short distances if the noise fluctuates rapidly. We furthermore demonstrate that a Berry–Esseen theorem leads to the vanishing of non-Gaussianities as the number of uncorrelated noise sources increases. Our predictions are corroborated by numerical simulations on realistic models of power grids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Testing quantum gravity using pulsed optomechanical systems

An interesting idea, dating back to Feynman [Report from Chapel Hill Conference, edited by C. M. DeWitt and D. Rickles (1957)], argues that quantum mechanics may break down for large masses if one entertains the possibility that gravity can be “classical,” thereby leading to predictions different from conventional low-energy quantum gravity. Despite the technical difficulty in testing such deviations, a large number of experimental proposals have been put forward due to the high level of fundamental interest. Here, we consider the Schrödinger-Newton (SN) theory and the correlated worldline (CWL) theory, and show that they can be distinguished from conventional quantum mechanics, as well as each other, by performing pulsed optomechanics experiments. For CWL specifically we develop a framework resembling the commonly used “Heisenberg-picture” treatment of coupled oscillators, allowing one to perform simple calculations for such systems without delving into the deeper path-integral formalism. We find that discriminating between the theories will be very difficult until experimental control over low frequency quantum optomechanical systems is pushed much further. Furthermore, the predicted departures of SN and CWL from quantum mechanics occur at the same scale, so both alternative models could in principle be probed by a single experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗