Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated data”

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,351 records · Page 75

Autotuning of Double-Dot Devices In Situ with Machine Learning

The current practice of manually tuning quantum dots (QDs) for qubit operation is a relatively time-consuming procedure that is inherently impractical for scaling up and applications. In this work, we report on the in situ implementation of a recently proposed autotuning protocol that combines machine learning (ML) with an optimization routine to navigate the parameter space. In particular, we show that a ML algorithm trained using exclusively simulated data to quantitatively classify the state of a double-QD device can be used to replace human heuristics in the tuning of gate voltages in real devices. Here, we demonstrate active feedback of a functional double-dot device operated at millikelvin temperatures and discuss success rates as a function of the initial conditions and the device performance. Modifications to the training network, fitness function, and optimizer are discussed as a path toward further improvement in the success rate when starting both near and far detuned from the target double-dot range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improved Vertical Carrier Transport for Green III-Nitride LEDs Using ( In , Ga ) N Alloy Quantum Barriers

We report on experimental and simulation-based results using ( In , Ga ) N alloy quantum barriers in c -plane green light-emitting diode (LED) structures as a means to improve vertical carrier transport and reduce forward voltage ( V F ) . Three-dimensional device simulations that include random alloy fluctuations are used to understand carrier behavior in a disordered potential. The simulated current density–voltage ( J - V ) characteristics and modified electron-hole overlap | F mod | 2 indicate that increasing the indium fraction in the ( In , Ga ) N quantum barriers leads to a reduced polarization discontinuity at the interface between the quantum barrier and quantum well, thereby reducing V F and improving | F mod | 2 . Maps of electron and hole current through the device show a relatively homogenous distribution in the X Y plane for structures using Ga N quantum barriers; in contrast, preferential pathways for vertical transport are identified in structures with ( In , Ga ) N barriers as regions of high and low current. A positive correlation between hole (electron) current in the p -side ( n -side) barrier and indium fraction reveals that preferential pathways exist in regions of high indium content. Furthermore, a negative correlation between the strain ε z z and indium fraction shows that high indium content regions have reduced strain-induced piezoelectric polarization in the Z direction due to the mechanical constraint of the surrounding lower indium content regions. Experimentally, multiple quantum well green LEDs with ( In , Ga ) N quantum barriers exhibit lower V F and blue-shifted wavelengths relative to LEDs with Ga N quantum barriers, consistent with simulation data. These results can be used to inform heterostructure design of low V F , long-wavelength LEDs and provide important insight into the nature of carrier transport in III-nitride alloy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Phase transitions and critical behavior in hadronic transport with a relativistic density functional equation of state

We develop a flexible, relativistically covariant parametrization of the dense nuclear matter equation of state suited for inclusion in computationally demanding hadronic transport simulations. Within an implementation in the hadronic transport code smash, we show that effects due to bulk thermodynamic behavior are reproduced in dynamic hadronic systems, demonstrating that hadronic transport can be used to study critical behavior in dense nuclear matter, both at and away from equilibrium. We also show that two-particle correlations calculated from hadronic transport simulation data follow theoretical expectations based on the second-order cumulant ratio, and constitute a clear signature of the crossover region above the critical point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Lattice model of heavy-light three-body system

We present a study of a $1+1$ dimensional heavy-light three-body system in finite volume. The heavy-light system is simulated by a coupled-channel $\phi^4$ type lattice model, and both ground state and excited states of multiparticle energy spectra are measured on various lattices. The lattice simulation data analysis is performed based on variational approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Neutrino oscillations in supernovae: Angular moments and fast instabilities

Recent theoretical work indicates that the neutrino radiation in core-collapse supernovae may be susceptible to flavor instabilities that set in far behind the shock, grow extremely rapidly, and have the potential to profoundly affect supernova dynamics and composition. Here we analyze the nonlinear collective oscillations that are prefigured by these instabilities. We demonstrate that a zero crossing in $n_{ν_e}$ – $n_{\overline ν_e}$ as a function of propagation angle is not sufficient to generate instability. Our analysis accounts for this fact and allows us to formulate complementary criteria. Using FORNAX simulation data, we show that fast collective oscillations qualitatively depend on how forward peaked the neutrino angular distributions are.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Marginal unbiased score expansion and application to CMB lensing

Here, we present the marginal unbiased score expansion (MUSE) method, an algorithm for generic high-dimensional hierarchical Bayesian inference. MUSE performs approximate marginalization over arbitrary non-Gaussian latent parameter spaces, yielding Gaussianized asymptotically unbiased and near-optimal constraints on global parameters of interest. It is computationally much cheaper than exact alternatives like Hamiltonian Monte Carlo (HMC), excelling on funnel problems which challenge HMC, and does not require any problem-specific user supervision like other approximate methods such as variational inference or many simulation-based inference methods. MUSE makes possible the first joint Bayesian estimation of the delensed Cosmic Microwave Background (CMB) power spectrum and gravitational lensing potential power spectrum, demonstrated here on a simulated data set as large as the upcoming South Pole Telescope 3G 1500 deg 2 survey, corresponding to a latent dimensionality of ~6 million and of order 100 global bandpower parameters. On a subset of the problem where an exact but more expensive HMC solution is feasible, we verify that MUSE yields nearly optimal results. We also demonstrate that existing spectrum-based forecasting tools which ignore pixel-masking underestimate predicted error bars by only ~10%. This method is a promising path forward for fast lensing and delensing analyses which will be necessary for future CMB experiments such as SPT-3G, Simons Observatory, or CMB-S4, and can complement or supersede existing HMC approaches. The success of MUSE on this challenging problem strengthens its case as a generic procedure for a broad class of high-dimensional inference problems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evidence for B + → K + ν ν ¯ decays

We search for the rare decay B + → K + ν ν ¯ in a 362 fb − 1 sample of electron-positron collisions at the ϒ ( 4 S ) resonance collected with the Belle II detector at the SuperKEKB collider. We use the inclusive properties of the accompanying B meson in ϒ ( 4 S ) → B B ¯ events to suppress background from other decays of the signal B candidate and light-quark pair production. We validate the measurement with an auxiliary analysis based on a conventional hadronic reconstruction of the accompanying B meson. For background suppression, we exploit distinct signal features using machine learning methods tuned with simulated data. The signal-reconstruction efficiency and background suppression are validated through various control channels. The branching fraction is extracted in a maximum likelihood fit. Our inclusive and hadronic analyses yield consistent results for the B + → K + ν ν ¯ branching fraction of [ 2.7 ± 0.5 ( stat ) ± 0.5 ( syst ) ] × 10 − 5 and [ 1.1 − 0.8 + 0.9 ( stat ) − 0.5 + 0.8 ( syst ) ] × 10 − 5 , respectively. Combining the results, we determine the branching fraction of the decay B + → K + ν ν ¯ to be [ 2.3 ± 0.5 ( stat ) − 0.4 + 0.5 ( syst ) ] × 10 − 5 , providing the first evidence for this decay at 3.5 standard deviations. The combined result is 2.7 standard deviations above the standard model expectation. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Supernova electron-neutrino interactions with xenon in the nEXO detector

Electron-neutrino charged-current interactions with xenon nuclei were modeled in the nEXO neutrinoless double-𝛽 decay detector (∼5 metric ton, 90% 136 Xe, 10% 134 Xe) to evaluate its sensitivity to supernova neutrinos. Predictions for event rates and detectable signatures were modeled using the Model of Argon Reaction Low Energy Yields (MARLEY) event generator. We find good agreement between MARLEY’s predictions and existing theoretical calculations of the inclusive cross sections at supernova neutrino energies. The interactions modeled by MARLEY were simulated within the nEXO simulation framework and were run through an example reconstruction algorithm to determine the detector’s efficiency for reconstructing these events. The simulated data, incorporating the detector response, were used to study the ability of nEXO to reconstruct the incident electron-neutrino spectrum and these results were extended to a larger xenon detector of the same isotope enrichment. We estimate that nEXO will be able to observe electron-neutrino interactions with xenon from supernovae as far as 5–8 kpc from Earth, while the ability to reconstruct incident electron-neutrino spectrum parameters from observed interactions in nEXO is limited to closer supernovae.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Nuclear structure functions at a future electron-ion collider

The quantitative knowledge of heavy nuclei's partonic structure is currently limited to rather large values of momentum fraction x { robust experimental constraints below x ~ 10 -2 at low resolution scale Q 2 are particularly scarce. This is in sharp contrast to the free proton's structure which has been probed in deep inelastic scattering (DIS) measurements down to x ~ 10 -5 at perturbative resolution scales. The construction of an Electron-Ion Collider (EIC) with a possibility to operate with a wide variety of nuclei, will allow one to explore the low-x region in much greater detail. In the present paper we simulate the extraction of the nuclear structure functions from measurements of inclusive and charm reduced cross sections at an EIC. The potential constraints are studied by analyzing simulated data directly in a next-to-leading order global fit of nuclear parton distribution functions based on the recent EPPS16 analysis. A special emphasis is placed on studying the impact an EIC would have on extracting the nuclear gluon PDF, the partonic component most prone to non-linear e ects at low Q 2 . In comparison to the current knowledge, we find that the gluon PDF can be measured at an EIC with significantly reduced uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Theory of capillary tension and interfacial dynamics of motility-induced phases

The statistical mechanics of equilibrium interfaces has been well-established for over a half century. In the past decade, a wealth of observations have made increasingly clear that a new perspective is required to describe interfaces arbitrarily far from equilibrium. In this work, beginning from microscopic particle dynamics that break time-reversal symmetry, we derive the linear interfacial dynamics of coexisting motility-induced phases. Doing so allows us to identify the athermal energy scale that excites interfacial fluctuations and the nonequilibrium surface tension that resists these excitations. Our theory identifies that, in contrast to equilibrium fluids, this active surface tension contains contributions arising from nonconservative forces which act to suppress interfacial fluctuations and, crucially, is distinct from the mechanical surface tension of Kirkwood and Buff. Here we find that the interfacial stiffness scales linearly with the intrinsic persistence length of the constituent active particle trajectories, in agreement with simulation data. We demonstrate that at wavelengths much larger than the persistence length, the interface obeys surface-area minimizing Boltzmann statistics with our derived nonequilibrium interfacial stiffness playing a role identical to that of equilibrium systems.

36 MATERIALS SCIENCE↗

Life cycle of streaks in the buffer layer of wall-bounded turbulence

In this work, streaks in the buffer layer of wall-bounded turbulence are tracked in time to study their life cycle. Spatially and temporally resolved direct numerical simulation data are used to analyze the strong wall-parallel movements conditioned to low-speed streamwise flow. The analysis of the streaks shows that there is a clear distinction between wall-attached and detached streaks, and that the wall-attached streaks can be further categorized into streaks that are contained in the buffer layer and the ones that reach the outer region. The results reveal that streaks are born in the buffer layer, coalescing with each other to create larger streaks that are still attached to the wall. Once the streak becomes large enough, it starts to meander due to the large streamwise-to-wall-normal aspect ratio, and consequently the elongation in the streamwise direction, which makes it more difficult for the streak to be oriented strictly in the streamwise direction. While the continuous interaction of the streaks allows the superstructure to span extremely long temporal and length scales, individual streak components are relatively small and short-lived. Tall-attached streaks eventually split into wall-attached and wall-detached components. These wall-detached streaks have a strong wall-normal velocity away from the wall, similar to ejections or bursts observed in the literature. Conditionally averaging the flow fields to these split events show that the detached streak has not only a larger wall-normal velocity compared to the wall-attached counterpart, it also has a larger (less negative) streamwise velocity, similar to the velocity field at the tip of a vortex cluster.

42 ENGINEERING↗

Internal energy balance and aerodynamic heating predictions for hypersonic turbulent boundary layers

The elemental equation governing heat transfer in aerodynamic flows is the internal energy equation. For a boundary layer flow, a double integration of the Reynolds-averaged form of this equation provides an expression of the wall heat flux in terms of the integrated effects, over the boundary layer, of various physical processes: turbulent dissipation, mean dissipation, turbulent heat flux, etc. Recently available direct numerical simulation data for a Mach 11 cold-wall turbulent boundary layer allows a comparison of the exact contributions of these terms in the energy equation to the wall heat flux with their counterparts modeled in the Reynolds-averaged Navier-Stokes (RANS) framework. Various approximations involved in RANS, both closure models as well as approximations involved in adapting incompressible RANS models to a compressible form, are assessed through examination of the internal energy balance. There are a number of potentially problematic assumptions and terms identified through this analysis. Here, the effect of compressibility corrections of the dilatational dissipation type is explored, as is the role of the modeled turbulent dissipation, in the context of wall heat flux predictions. The results indicate several potential avenues for RANS model improvement for hypersonic cold-wall boundary-layer flows.

42 ENGINEERING↗

Revival of the Dark Matter Hypothesis for the Galactic Center Gamma-Ray Excess

Statistical evidence has previously implied that the galactic center GeV excess (GCE) originates largely from point sources, and not from annihilating dark matter. Here, we examine the impact of unmodeled source populations on identifying the true origin of the GCE using non-Poissonian template fitting (NPTF) methods. In a proof-of-principle example with simulated data, we discover that unmodeled sources in the Fermi bubbles can lead to a dark matter signal being misattributed to point sources by the NPTF. We discover striking behavior consistent with a mismodeling effect in the real Fermi data, finding that large artificial injected dark matter signals are completely misattributed to point sources. Consequently, we conclude that dark matter may provide a dominant contribution to the GCE after all.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effects of dopants on the glass forming ability in Al-based metallic alloy

The effect of dopants on the metallic glass forming ability is usually considered based on the analysis of changes in the liquid structure or thermodynamics. What is missing in such considerations is an analysis of how a dopant changes the properties of the crystal phases which can form instead of the glass. In order to illuminate this aspect we performed molecular dynamics simulations to study the effects of Mg and Sm dopants on the crystal nucleation in Al. The simulation data were found to be consistent with the experimental observations that addition of Mg to Al does not lead to vitrification but addition of only 8% Sm does. The significant effect of Sm doping was related to the intolerance of Al to this dopant. This leads to increase in the solid-liquid interfacial free energy, and therefore, to increase in the nucleation barrier and to a dramatic decrease in the nucleation rate. Here, the intolerance mechanism also significantly affects the growth kinetics.

36 MATERIALS SCIENCE↗

Universal L - 3 finite-size effects in the viscoelasticity of amorphous systems

We present a theory of viscoelasticity of amorphous media, which takes into account the effects of confinement along one of three spatial dimensions. The framework is based on the nonaffine extension of lattice dynamics to amorphous systems, or nonaffine response theory. The size effects due to the confinement are taken into account via the nonaffine part of the shear storage modulus G'. The nonaffine contribution is written as a sum over modes in k-space. With a rigorous argument based on the analysis of the k-space integral over modes, it is shown that the confinement size L in one spatial dimension, e.g., the z axis, leads to a infrared cutoff for the modes contributing to the nonaffine (softening) correction to the modulus that scales as L -3 . Corrections for finite sample size D in the two perpendicular dimensions scale as ~(L/D) 4 , and are negligible for L«D. For liquids it is predicted that G'~L -3 is in agreement with a previous more approximate analysis, whereas for amorphous materials G'~G' bulk +βL -3 . For the case of liquids, four different experimental systems are shown to be very well described by the L-3 law. The theory can also explain previous simulation data of confined jammed granular packings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optimizers for stabilizing likelihood-free inference

A growing number of applications in particle physics and beyond use neural networks as unbinned likelihood ratio estimators applied to real or simulated data. Precision requirements on the inference tasks demand a high-level of stability from these networks, which are affected by the stochastic nature of training. We show how physics concepts can be used to stabilize network training through a physics-inspired optimizer. In particular, the energy conserving descent (ECD) optimization framework uses classical Hamiltonian dynamics on the space of network parameters to reduce the dependence on the initial conditions while also stabilizing the result near the minimum of the loss function. We develop a version of this optimizer known as , which has few free hyperparameters with limited ranges guided by physical reasoning. We apply to representative likelihood-ratio estimation tasks in particle physics and find on average that it out-performs the widely used Adam optimizer. We expect that ECD will be a useful tool for wide array of data-limited problems, where it is computationally expensive to exhaustively optimize hyperparameters and mitigate fluctuations with ensembling.

Monte Carlo methods↗

Combinatorial appraisal of transition states for in situ pair distribution function analysis

In situ total scattering measurements are increasingly utilized to follow atomic and nanoscale structural details of phase transitions and other transient processes in materials. This contribution presents an automated method and associated tool set to analyze series of diffraction and pair distribution function data with a linear combination of end-member states. It is demonstrated that the combinatorial appraisal of transition states (CATS) software tracks phase changes, relative phase fractions and length scales of interest in experimental data series. Here it is further demonstrated, using a series of local structure data simulations, that the misfit of such a model can reveal details of phase aggregation and growth related to the pair distribution function's sensitivity to interphase correlations. CATS may be applied to quantitative evaluation of many transient processes, including amorphous-to-crystalline phase transitions, the evolution of solid-solution behaviors, the precipitation and growth of aggregates, and other atomic to nanoscale details of crystallization and phase transformation phenomena.

36 MATERIALS SCIENCE↗

Wavelength-resolved Neutron Tomography for Crystalline Materials

Wavelength-resolved (WR) neutron transmission tomography is an emerging technique to characterize engineering materials. While tomographic reconstruction for amorphous samples is straightforward, it is challenging to reconstruct samples with single-crystal domains because the attenuation of the sample varies as a function of its orientation with respect to the incident beam due to Bragg scattering. In this paper, we present an algorithm that can reconstruct samples with single-crystal domains from WR neutron tomographic measurements. In particular, we use a model-based iterative reconstruction (MBIR) technique that reconstructs the volume by identifying and leaving out the regions of the measurement that are affected by Bragg scatter. We combine the output of the MBIR method with an algorithm that matches the reconstruction to the identified Bragg scatter to reconstruct a feature that corresponds to the local crystallography of the sample being measured. Using simulated data, we demonstrate how our algorithm can reconstruct materials with single-crystal domains, thereby adding a powerful new capability for WR neutron imaging instruments.

Singanallur Vaidyanathan, Venkatakrishnan↗