Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STATISTICAL MECHANICS”

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 1,171 records · Page 65

On the Existence of the Plateau Emission in High-energy Gamma-Ray Burst Light Curves Observed by Fermi-LAT

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) shows long-lasting high-energy emission in many gamma-ray bursts (GRBs), similar to X-ray afterglows observed by the Neil Gehrels Swift Observatory (Swift). Some LAT light curves (LCs) show a late-time flattening reminiscent of X-ray plateaus. We explore the presence of plateaus in LAT temporally extended emission analyzing GRBs from the second Fermi-LAT GRB Catalog from 2008 to 2016 May with known redshifts, and check whether they follow closure relations corresponding to four distinct astrophysical environments predicted by the external forward shock model. We find that three LCs can be fit by the same phenomenological model used to fit X-ray plateaus and show tentative evidence for the existence of plateaus in their high-energy extended emission. The most favorable scenario is a slow-cooling regime, whereas the preferred density profile for each GRBs varies from a constant-density interstellar medium to an r –2 wind environment. We also compare the end time of the plateaus in γ-rays and X-rays using a statistical comparison with 222 Swift GRBs with plateaus and known redshifts from 2005 January to 2019 August. Within this comparison, the case of GRB 090510 shows an indication of chromaticity at the end time of the plateau. Finally, we update the 3D fundamental plane relation among the rest-frame end time of the plateau, its correspondent luminosity, and the peak prompt luminosity for 222 GRBs observed by Swift. We find that these three LAT GRBs follow this relation.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of Heavy Flavor Mesons and Flavor-Tagged Jets with the CMS Detector

The goal of this research program is to implement heavy flavor meson triggers in heavy-ion collisions for the Compact Muon Solenoid (CMS) experiment at the Large Hadron Collider (LHC) at CERN, including algorithm design, timing studies, offline validation, and online performance monitoring. The physics analyses which can be achieved by data from these new triggers is to address one of the most important questions in the field: parton flavor dependence of jet-quenching for the understanding of the transport properties of the Quark-Gluon Plasma. This program will allow CMS to collect the highest statistics heavy flavor meson and jet data ever recorded in heavy-ion colliders. The program includes two objectives: (1) Build and maintain the heavy flavor meson and jet triggers for heavy-ion collisions and deploy the trigger algorithms for 2015-2018 PbPb and pPb run at the LHC; (2) Perform heavy flavor meson and jet physics analyses, which can be used to study the parton flavor dependence of jet quenching, to extract the elastic energy loss coefficient of the QGP, and to test whether massive quarks also participate in collective expansion dynamics in heavy-ion collisions. With the heavy flavor physics trigger developed in this project, a competitive heavy flavor physics program in heavy-ion collisions has been established in CMS. This program allows studies of the fully reconstructed and flavor identified charm, beauty, and exotic hadrons that cover the widest transverse momentum range. The novel measurements supported by the award provide new constraints on the size of the flavor dependence of parton energy loss, the value of the in-medium charm quark diffusion coefficient, the mechanism of charm and beauty quark hadronization, and provide new insights to the nature of the X(3872) hadron.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Physics-informed Karhunen-Loeve and Neural Network Approximations for Solving Inverse Differential Equation Problems

Here we present the PI-CKL-NN method for parameter estimation in differential equation (DE) models given sparse measurements of the parameters and states. In the proposed approach, the space- or time-dependent parameters are approximated by Karhunen-Loeve (KL) expansions that are conditioned on the parameters’ measurements, and the states are approximated by deep neural networks (DNNs). The unknown weights in the KL expansions and DNNs are found my minimizing the cost function that enforces the measurements of the states the DE constraint. Regularization is achieved by adding the l2 norm of the conditional KL coefficients into the loss function. Our approach assumes that the parameter fields are correlated in space or time and enforces the statistical knowledge (the mean and the covariance function) in addition to the DE constraints and measurements as opposed to the physics-informed neural network (PINN) and other similar physics-informed machine learning methods where only DE constraints and data are used for parameter estimation. We use the PI-CKL-NN method for parameter estimation in an ordinary differential equation with an unknown time-dependent parameter and the one- and two-dimensional partial differential diffusion equations with unknown space-dependent diffusion coefficients. We also demonstrate that PI-CKL-NN is more accurate than the PINN method, especially when the observations of the parameters are very sparse

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Next-generation even-denominator fractional quantum Hall states of interacting composite fermions

The discovery of the fractional quantum Hall state (FQHS) in 1982 ushered a new era of research in many-body condensed matter physics. Among the numerous FQHSs, those observed at even-denominator Landau level filling factors are of particular interest as they may host quasiparticles obeying non-Abelian statistics and be of potential use in topological quantum computing. The even-denominator FQHSs, however, are scarce and have been observed predominantly in low-disorder two-dimensional (2D) systems when an excited electron Landau level is half filled. An example is the well-studied FQHS at filling factor ν = 5/2 which is believed to be a Bardeen-Cooper-Schrieffer-type, paired state of flux-particle composite fermions (CFs). Here, we report the observation of even-denominator FQHSs at ν = 3/10, 3/8, and 3/4 in the lowest Landau level of an ultrahigh-quality GaAs 2D hole system, evinced by deep minima in longitudinal resistance and developing quantized Hall plateaus. Quite remarkably, these states can be interpreted as even-denominator FQHSs of CFs, emerging from pairing of higher-order CFs when a CF Landau level, rather than an electron or a hole Landau level, is half-filled. Our results affirm enhanced interaction between CFs in a hole system with significant Landau level mixing and, more generally, the pairing of CFs as a valid mechanism for even-denominator FQHSs, and suggest the realization of FQHSs with non-Abelian anyons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental measurements of fluid flow in an 84-pin hexagonal rod bundle with spacer grid for a gas-cooled fast modular reactor

A 50 MW e helium-cooled fast modular reactor (FMR) is under development, and as part of the Department of Energy Integrated Research Project (IRP), Texas A&M University is conducting the thermal-hydraulic characterization of its baseline core configuration. Here, we experimentally investigated the axial and cross flow fields characteristics in the hexagonal fuel rod bundle composed of 84 rods, a central rod, and spacer grids at a Reynolds number of 12,000. A fully transparent experimental facility resembling of one unit of the fuel assembly was constructed. Time-resolved particle image velocimetry (TR-PIV) measurements were performed to characterize hydraulic behavior downstream the spacer grid. Velocity measurements were conducted in the axial and radial direction of the rod bundle. From the PIV velocity vector fields, the full-field flow statistics were computed for the mean velocity, vorticity, and Reynolds stresses. An analysis of the energy decay downstream of the spacer grid was conducted with calculations of secondary-flow intensities, turbulent kinetic energy, power spectrum, and spatial-temporal velocity cross correlations. The vorticity field was obtained and pairs of counter-rotating vortexes were identified using a 3D reconstruction of several measurement planes. The data from this experimental campaign will be used to validate numerical models based on Computational Fluid Dynamics and other codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simultaneous velocity and density measurements of fully developed Rayleigh-Taylor mixing

Here, the dynamics of molecular mixing and the energy transfer process in the Rayleigh-Taylor instability (RTI) are studied through the collection of simultaneous density-velocity measurements. These experiments provide simultaneous density-velocity field measurements, in contrast to previous point measurements. Statistically stationary experiments are performed in a “convective-type” gas tunnel facility, with density contrast achieved through the injection of helium into the bottom stream. Three experiments at Atwood number ≈ 0.1 are captured at three outer scale Reynolds numbers Re = 520, 2260, and 4050. Particle image velocimetry and laser induced fluorescence are employed simultaneously. Statistics of the density and velocity show self-similar collapse of RTI profiles at large Reynolds number Re >2000. Flat velocity profiles indicate homogeneous turbulence characteristics in the core of the mixing region. Significant anisotropy develops in the flow, with horizontal velocity fluctuations being only 60% of the vertical velocity fluctuations. The turbulent mass flux is found to be asymmetric about the centerline, with increased peak towards the spike. Measurements of the molecular mixing show that mixing is maximized at the core of the flow and increases with increased Reynolds number. The transport equation of density-specific-volume correlation b shows that it is mostly produced in the core of the mixing region, and that the spatial evolution of its profile is the result of transport by bulk motion of the bubble and spike. Energy transfer from gravitational potential energy to turbulent kinetic energy and viscous dissipation is observed to occur in the experiment with a ratio of dissipated energy to potential energy released of 38%. The analysis of the turbulent kinetic energy transport equation budget reveals that production is the primary mechanism towards the growth of turbulent kinetic energy in the core of the flow, and is asymmetrically slightly skewed towards the spike. However, it is through the transport that the strong advection at the edges of the mixing region is maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tracking atomic structure evolution during directed electron beam induced Si-atom motion in graphene via deep machine learning

Using electron beam manipulation, we enable deterministic motion of individual Si atoms in graphene along predefined trajectories. Structural evolution during the dopant motion was explored, providing information on changes of the Si atom neighborhood during atomic motion and providing statistical information of possible defect configurations. The combination of a Gaussian mixture model and principal component analysis applied to the deep learning-processed experimental data allowed disentangling of the atomic distortions for two different graphene sublattices. This approach demonstrates the potential of e-beam manipulation to create defect libraries of multiple realizations of the same defect and explore the potential of symmetry breaking physics. The rapid image analytics enabled via a deep learning network further empowers instrumentation for e-beam controlled atom-by-atom fabrication. Here, the analysis described in the paper can be reproduced via an interactive Jupyter notebook at https://git.io/JJ3Bx.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantifying Uncertainty in All-to-All Estimates of Space Object Conjunction Probabilities using U-Statistics

Predicting space object conjunctions is inherently probabilistic due to initial state and orbit model uncertainty. A commonly considered Monte Carlo estimator of the conjunction probability is the ’all-to-all’ estimator. Given independent random samples of the trajectories of both objects, the estimator is the percentage of all pairs of trajectories that result in a conjunction. Intuitively, the all-to-all estimator is the best possible estimator of the conjunction probability since it considers all pairs of Monte Carlo samples. However, its distribution is not available in closed-form, which limits its use in practice and makes this intuition difficult to make rigorous. In this paper, the all-to-all estimator is identified as a U-statistic, which implies that it has several favorable properties. Specifically, the estimator is the minimum variance unbiased estimator of the conjunction probability and is asymptotically Gaussian distributed. An approximate confidence interval for the conjunction probability is obtained from an estimate of the asymptotic Gaussian distribution. We show how to efficiently compute the confidence interval and demonstrate that the interval has the nominal coverage level. The confidence intervals are also seen to be narrower than those based on the commonly-used each-to-each estimator. Furthermore, the all-to-all estimator is shown to allow different Monte Carlo sample sizes, whereas the each-to-each estimator requires equal sample sizes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Role of Chemical Disorder in High Temperature Dislocation Glide in Refractory Multi-principal Element Alloys

Refractory multi-principal element alloys (RMPEAs) combine a chemically disordered lattice with a structurally ordered, single‐phase body‐centered cubic (bcc) crystal structure. Chemical fluctuations in these alloys give rise to significant energy barriers that impede dislocation motion. In this study, we use phase field dislocation dynamics to examine how spatial temperature fluctuations compete with randomness in energy barriers to affect the motion of long screw dislocations in three equi-atomic MoNbTa‐based RMPEAs: MoNbTa, MoNbTaW, and MoNbTaVW. Over a wide range of homologous temperatures (T h ≈ 0–0.6), we determined a screw dislocation ‘flow stress’ as the minimum applied stress to sustain continuous motion over a long excursion distance within a fixed timeframe. All three RMPEAs exhibit temperature dependent flow stress with three characteristic glide regimes: at low homologous temperatures, flow stress drops sharply, and screw glide remains planar and rectilinear; at intermediate homologous temperatures, the flow stress levels off as glide becomes planar but wavy; and at high homologous temperatures, flow stress plateaus as screws exhibit nonplanar, three‐dimensional motion. Glide kinetics and transition temperatures are controlled by the chemically induced fluctuations in the energy landscape. The rectilinear to wavy transition temperature is controlled by statistically weakest local barriers in the glide plane, whereas the wavy to 3D transition temperature is governed by statistically strongest local barriers. At low and intermediate homologous temperatures, the relative spread in energy barriers governs glide behavior by controlling the local kinetics of kink pair formation and kink pinning. At high homologous temperatures, the average barrier height governs the glide behavior by controlling the number of out-of-plane excursions during 3D glide. These findings reveal how random chemical fluctuations determine screw‐driven plasticity in RMPEAs, providing critical insight for the design of high temperature structural alloys.

Defects↗

Topology and many-body localization

We discuss the problem of localization in two dimensional electron systems in the quantum Hall (single Landau level) regime. After briefly summarizing the well-studied problem of Anderson localization in the non-interacting case, we concentrate on the problem of disorder induced many-body localization (MBL) in the presence of electron–electron interactions using numerical exact diagonalization and eigenvalue spacing statistics as a function of system size. We provide evidence showing that MBL is not attainable in a single Landau level with short range (white noise) disorder in the thermodynamic limit. We then study the interplay of topology and localization, by contrasting the behavior of topological and nontopological subbands arising from a single Landau level in two models — (i) a pair of extremely flat Hofstadter bands with an optimally chosen periodic potential, and (ii) a Landau level with a split-off nontopological impurity band. Both models provide convincing evidence for the strong effect of topology on the feasibility of many-body localization as well as slow dynamics starting from a nonequilibrium state with charge imbalance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mapping strain and structural heterogeneities around bubbles in amorphous ionically conductive Bi 2 O 3

While amorphous materials are often approximated to have a statistically homogeneous atomic structure, they frequently exhibit localized structural heterogeneity that challenges simplified models. This study uses 4D scanning transmission electron microscopy to investigate the strain and structural modifications around gas bubbles in amorphous Bi 2 O 3 induced by argon irradiation. We present a method for determining strain fields surrounding bubbles that can be used to measure the internal pressure of the gas. Compressive strain is observed around the cavities, with higher-order crystalline symmetries emerging near the cavity interfaces, suggesting paracrystalline ordering as a result of bubble coarsening. This ordering, along with a compressive strain gradient, indicates that gas bubbles induce significant localized changes in atomic packing. By analyzing strain fields with maximum compressive strains of 3%, we estimate a lower bound on the internal pressure of the bubbles at 2.5 GPa. These findings provide insight into the complex structural behavior of amorphous materials under stress, particularly in systems with gas inclusions, and offer new methods for probing the local atomic structure in disordered materials. Although considering structural heterogeneity in amorphous systems is non-trivial, these features have crucial impacts on material functionalities, such as mechanical strength, ionic conductivity, and electronic mobility.

36 MATERIALS SCIENCE↗

Oxidation Dynamics of Supported Catalytic Cu Clusters: Coupling to Fluxionality

Copper oxide nanoclusters have a wide range of catalytic applications, such as the selective oxidation of hydrocarbons. O 2 binding to the catalyst, activation, and release upon reagent oxidation are key events in these catalytic chemistries. Furthermore, these events are expected to be accompanied by significant structural changes of the Cu clusters, because O atoms integrate into the cluster, rather than bind to its surface. Topping the complexity of the problem, partially oxidized Cu clusters are known to exhibit strong fluxionality and feature diverse and interconverting structures and oxygen contents in conditions of oxidative dehydrogenation (ODH). Hence, a significant dynamic coupling between the “hot” O 2 molecule impacting the cluster at reaction temperatures and the cluster fluxionality can be expected. In this work, we focus on the dynamics of dioxygen integration into a partially oxidized Cu cluster supported on hydroxylated amorphous alumina–a system recently reported to be an exceptionally selective catalyst for cyclohexane ODH with very little CO produced, whose mechanistic underpinnings are of utmost interest. The statistics over a swarm of adsorption and scattering trajectories where O 2 hits various sites on the cluster at reaction temperature shows that the O 2 binding does not only follow the minimal energy paths. O 2 also rarely integrates into the cluster in a single step and instead first binds to a single Cu atom via either an η 1 -O 2 or an η 2 -O 2 mode. Surprisingly, this step often has a higher barrier than the subsequent O 2 integration and dissociation, which in turn take multiple steps and complete the oxidation process. Dynamic trajectories starting from the key transition state of integration of the adsorbed O 2 can also lead to different intermediate structures during or right after the dissociation, due to the energy released from the transition state and the thermal intracluster effects. From these activated O 2 chemisorbed structures, O 2 dissociation occurs with moderate barriers (~0.5 eV), producing multiple final oxidized Cu 4 O 4 states. Hence, a diversity of reaction profiles for the attack of supported Cu cluster by O 2 emerges due to the dynamic effects, with implications for mechanisms, kinetic models, and catalyst design principles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulating Z 2 lattice gauge theory on a quantum computer

The utility of quantum computers for simulating lattice gauge theories is currently limited by the noisiness of the physical hardware. Various quantum error mitigation strategies exist to reduce the statistical and systematic uncertainties in quantum simulations via improved algorithms and analysis strategies. We perform quantum simulations of Z 2 gauge theory with matter to study the efficacy and interplay of different error mitigation methods: readout error mitigation, randomized compiling, rescaling, and dynamical decoupling. We compute Minkowski correlation functions in this confining gauge theory and extract the mass of the lightest spin-1 state from fits to their time dependence. Quantum error mitigation extends the range of times over which our correlation function calculations are accurate by a factor of 6 and is therefore essential for obtaining reliable masses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dust pollution in China affected by different spatial and temporal types of El Niño

Dust is an important aerosol affecting air quality in China in the winter and spring seasons. Dust in China is potentially influenced by the interannual climate variability associated with El Niño. Here, the impacts of El Niño with different temporal and spatial types on dust pollution in boreal winter and spring in China and the potential mechanisms are investigated using a state-of-the-art Earth system model (E3SMv1). We find that the eastern Pacific (EP) and central Pacific (CP) El Niño both increase wintertime dust concentrations by 5–50 µg m –3 over central-eastern China. Due to a stronger wind and lower relative humidity, which favor dust emissions near sources, and a strengthened northwesterly and reduced precipitation, which are conducive to dust transport, dust concentrations during the CP El Niño are 5–20 µg m –3 higher in northern China than during the EP El Niño, although the changes are mostly insignificant. El Niño with a short duration (SD) increases boreal winter dust concentrations by 20–100 µg m –3 over northern China relative to the climatological mean, while there is a decrease of 5–50 µg m –3 during the long-duration (LD) El Niño, which is also related to the El Niño-induced changes in atmospheric circulation, precipitation, and relative humidity. In the following spring season, all types of El Niño events enhance dust over northern China, but only the increase during the LD El Niño is statistically significant, suggesting that the weaker intensity but longer duration of the LD El Niño events can significantly affect spring dust in China. Our results contribute to the current knowledge of the influence of El Niño on dust pollution, which has profound implications for air pollution control and dust storm prediction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Approximating the linear response of physical chaos

Abstract Parametric derivatives of statistics are highly desired quantities in prediction, design optimization and uncertainty quantification. In the presence of chaos, the rigorous computation of these quantities is certainly possible, but mathematically complicated and computationally expensive. Based on Ruelle’s formalism, this paper shows that the sophisticated linear response algorithm can be dramatically simplified in higher-dimensional systems featuring statistical homogeneity in the physical space. We argue that the contribution of the SRB (Sinai–Ruelle–Bowen) measure gradient, which is an integral yet the most cumbersome part of the full algorithm, is negligible if the objective function is appropriately aligned with unstable manifolds. This abstract condition could potentially be satisfied by a vast family of real-world chaotic systems, regardless of the physical meaning and mathematical form of the objective function and perturbed parameter. We demonstrate several numerical examples that support these conclusions and that present the use and performance of a simplified linear response algorithm. In the numerical experiments, we consider physical models described by differential equations, including Lorenz 96 and Kuramoto–Sivashinsky.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Unitary Subharmonic Response and Floquet Majorana Modes

Detection and manipulation of excitations with non-Abelian statistics, such as Majorana fermions, are essential for creating topological quantum computers. To this end, we show the connection between the existence of such localized particles and the phenomenon of unitary subharmonic response (SR) in periodically driven systems. In particular, starting from highly nonequilibrium initial states, the unpaired Majorana modes exhibit spin oscillations with twice the driving period, are localized, and can have exponentially long lifetimes in clean systems. While the lifetime of SR is limited in translationally invariant systems, we show in this work that disorder can be engineered to stabilize the subharmonic response of Majorana modes. A viable observation of this phenomenon can be achieved using modern multiqubit hardware, such as superconducting circuits and cold atomic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Transient hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation. II. Isomer-resolved unimolecular dynamics

Transient carbon-centered hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation are characterized by their time- and energy-resolved unimolecular dissociation dynamics to hydroxyl (OH) and cyclic ether products. Two distinct •QOOH isomers are examined with radical sites at a primary carbon of one of the methyl groups (β-Me) or a secondary carbon (β-Et) of the ethyl group. Energy-dependent unimolecular rates are obtained from the time-dependent appearance of OH products for the two isomers and compared with statistical microcanonical rates computed using RRKM theory, including heavy-atom tunneling, based on high-level theoretical calculations. A benchmark-corrected approach is utilized to compute high-accuracy stationary-point energies, most importantly, transition-state barriers, for the •QOOH_Me and •QOOH_Et isomers in isopentane oxidation, building on higher-level reference calculations for the oxidation of ethane (C 2 H 5 O 2 ) and propane (C 3 H 7 O 2 ), respectively. Here, the measured rates are compared with RRKM calculations incorporating the benchmark-corrected transition-state parameters, a vibrationally adiabatic multidimensional hindered-rotor treatment of key torsions, and quantum tunneling. Agreement between experiment and theory validates the statistical description and shows faster decay for •QOOH_Et due to its lower barrier. Both β-QOOH isomers decay almost exclusively to OH + cyclic ether products under the conditions studied.

Oxidation processes↗