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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effect of stress on Laves phase precipitation in creep ruptured Grade 92 ferritic martensitic steel characterized by a novel accessible method

At high temperature conditions relevant to fossil and nuclear energy plants, Laves phase (Fe 2 X, X=Mo, W) precipitation is observed in common ferritic martensitic (FM) structural steels, with various reported effects on creep behavior. Despite being valuable metrics to correlate with mechanical properties and other precipitate phases, the volume fraction and number density of Laves phase precipitates has been difficult to quantify accurately using common techniques such as transmission electron microscopy (TEM) due to the relatively large size (∼0.25 μm) and low number density (∼10 11 cm −3 ) of Laves precipitates. Here, to address this characterization challenge, we developed and demonstrated a high-throughput and widely accessible method to quantify the volume fraction and number density of the Laves phase based on scanning electron microscope (SEM) images with a backscattered electron signal and the information depth (ID) of backscattered electrons. We applied this new technique in creep ruptured Grade 92 FM steel to study the effect of Laves phase on creep properties and determine the influence of stress on Laves phase precipitation. The quantitative accuracy of the SEM-based volume fraction and number density values was verified using synchrotron high energy X-ray diffraction and serial sectioning tomography. Stress did not significantly affect the Laves phase size or volume fraction during creep testing at 550 – 650°C and stress levels of 90 – 260 MPa (vs. unstressed conditions). Conversely, a moderate but statistically significant stress-enhanced increase in Laves phase number density, corresponding to an increase in nucleation rate, occurred during creep exposure above 110 MPa.

Creep↗
DOE OSTI2024

X-ray imaging and electron temperature evolution in laser-driven magnetic reconnection experiments at the national ignition facility

We present results from x-ray imaging of high-aspect-ratio magnetic reconnection experiments driven at the National Ignition Facility. Two parallel, self-magnetized, elongated laser-driven plumes are produced by tiling 40 laser beams. A magnetic reconnection layer is formed by the collision of the plumes. A gated x-ray framing pinhole camera with micro-channel plate detector produces multiple images through various filters of the formation and evolution of both the plumes and current sheet. As the diagnostic integrates plasma self-emission along the line of sight, two-dimensional electron temperature maps ⟨Te⟩Y are constructed by taking the ratio of intensity of these images obtained with different filters. The plumes have a characteristic temperature ⟨Te⟩Y=240 ± 20 eV at 2 ns after the initial laser irradiation and exhibit a slow cooling up to 4 ns. The reconnection layer forms at 3 ns with a temperature ⟨Te⟩Y=280 ± 50 eV as the result of the collision of the plumes. The error bars of the plumes and current sheet temperatures separate at 4 ns, showing the heating of the current sheet from colder inflows. Using a semi-analytical model, we survey various heating mechanisms in the current sheet. We find that reconnection energy conversion would dominate at low density (ne≲7×1018 cm−3) and electron-ion collisional drag at high-density (≳1019 cm−3).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Modeling noise in global Mølmer-Sørensen interactions applied to quantum approximate optimization

Many-qubit Mølmer-Sørensen (MS) interactions applied to trapped ions offer unique capabilities for quantum information processing, with applications including quantum simulation and the quantum approximate optimization algorithm (QAOA). Here, we develop a physical model to describe many-qubit MS interactions under four sources of experimental noise: vibrational mode frequency fluctuations, laser power fluctuations, thermal initial vibrational states, and state preparation and measurement errors. The model parametrizes these errors from simple experimental measurements, without free parameters. Herein, we validate the model in comparison with experiments that implement sequences of MS interactions on two 171 Yb + ions. The model shows reasonable agreement after several MS interactions as quantified by the reduced chi-squared statistic $χ^{2}_{red}$≈2. As an application we examine MaxCut QAOA experiments on three and six ions. The experimental performance is quantified by approximation ratios that are 91% and 83% of the optimal theoretical values. Our model predicts $0.93^{+0.03}_{-0.02}$ and $0.95^{+0.04}_{-0.03}$, respectively, with disagreement in the latter value attributable to secondary noise sources beyond those considered in our analysis. With realistic experimental improvements to reduce measurement error and radial trap frequency variations, the model achieves approximation ratios that are 99% of the optimal. Incorporating these improvements into future experiments is expected to reveal new aspects of noise for future modeling and experimental improvements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The completed SDSS-IV extended baryon oscillation spectroscopic survey: pairwise-inverse probability and angular correction for fibre collisions in clustering measurements

ABSTRACT The completed extended Baryon Oscillation Spectroscopic Survey (eBOSS) catalogues contain redshifts of 344 080 quasars at 0.8 < z < 2.2, 174 816 luminous red galaxies between 0.6 < z < 1.0, and 173 736 emission-line galaxies over 0.6 < z < 1.1 in order to constrain the expansion history of the Universe and the growth rate of structure through clustering measurements. Mechanical limitations of the fibre-fed spectrograph on the Sloan telescope prevent two fibres being placed closer than 62 arcsec in a single pass of the instrument. These ‘fibre collisions’ strongly correlate with the intrinsic clustering of targets and can bias measurements of the two-point correlation function resulting in a systematic error on the inferred values of the cosmological parameters. We combine the new techniques of pairwise-inverse probability and the angular upweighting (PIP+ANG) to correct the clustering measurements for the effect of fibre collisions. Using mock catalogues, we show that our corrections provide unbiased measurements, within data precision, of both the projected $\rm {\mathit{ w}_p}\left(\mathit{ r}_p\right)$ and the redshift-space multipole ξ(ℓ = 0, 2, 4)(s) correlation functions down to $0.1\, h^{-1}{\rm Mpc}$, regardless of the tracer type. We apply the corrections to the eBOSS DR16 catalogues. We find that, on scales $s\gtrsim 20\, h^{-1}{\rm Mpc}$ for ξℓ, as used to make baryon acoustic oscillation and large-scale redshift-space distortion measurements, approximate methods such as nearest-neighbour upweighting are sufficiently accurate given the statistical errors of the data. Using the PIP method, for the first time for a spectroscopic program of the Sloan Digital Sky Survey, we are able to successfully access the one-halo term in the clustering measurements down to $\sim 0.1\, h^{-1}{\rm Mpc}$ scales. Our results will therefore allow studies that use the small-scale clustering to strengthen the constraints on both cosmological parameters and the halo occupation distribution models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Persistent Protein Motions in a Rugged Energy Landscape Revealed by Normal Mode Ensemble Analysis

Proteins are allosteric machines that couple motions at distinct, often distant, sites to control biological function. Low-frequency structural vibrations are a mechanism of this long-distance connection and are often used computationally to predict correlations, but experimentally identifying the vibrations associated with specific motions has proved challenging. Spectroscopy is an ideal tool to explore these excitations, but measurements have been largely unable to identify important frequency bands. The result is at odds with some previous calculations and raises the question what methods could successfully characterize protein structural vibrations. Here we show the lack of spectral structure arises in part from the variations in protein structure as the protein samples the energy landscape. However, by averaging over the energy landscape as sampled using an aggregate 18.5 μs of all-atom molecular dynamics simulation of hen egg white lysozyme and normal-mode analyses, we find vibrations with large overlap with functional displacements are surprisingly concentrated in narrow frequency bands. These bands are not apparent in either the ensemble averaged vibrational density of states or isotropic absorption. However, in the case of the ensemble averaged anisotropic absorption, there is persistent spectral structure and overlap between this structure and the functional displacement frequency bands. We systematically lay out heuristics for calculating the spectra robustly, including the need for statistical sampling of the protein and inclusion of adequate water in the spectral calculation. The results show the congested spectrum of these complex molecules obscures important frequency bands associated with function and reveal a method to overcome this congestion by combining structurally sensitive spectroscopy with robust normal mode ensemble analysis.

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

BPS chaos

Black holes are chaotic quantum systems that are expected to exhibit random matrix statistics in their finite energy spectrum. Lin, Maldacena, Rozenberg and Shan (LMRS) have proposed a related characterization of chaos for the ground states of BPS black holes with finite area horizons. On a separate front, the “fuzzball program” has uncovered large families of horizon-free geometries that account for the entropy of holographic BPS systems, but only in situations with sufficient supersymmetry to exclude finite area horizons. The highly structured, non-random nature of these solutions seems in tension with strong chaos. We verify this intuition by performing analytic and numerical calculations of the LMRS diagnostic in the corresponding boundary quantum system. In particular we examine the 1/2 and 1/4-BPS sectors of \mathcal{N}=4 𝒩 = 4 SYM, and the two charge sector of the D1-D5 CFT. We find evidence that these systems are only weakly chaotic, with a Thouless time determining the onset of chaos that grows as a power of N N . In contrast, finite horizon area BPS black holes should be strongly chaotic, with a Thouless time of order one. In this case, finite energy chaotic states become BPS as N N is decreased through the recently discovered “fortuity” mechanism. Hence they can plausibly retain their strongly chaotic character.

Chen, Yiming (ORCID:0000000218613230)↗
DOE OSTI2020

Stochastic nonlinear analysis of unidirectional fiber composites using image-based microstructural uncertainty quantification

Here, we present a data-driven nonlinear uncertainty quantification and propagation framework to study the microstructure-induced stochastic performance of unidirectional (UD) carbon fiber reinforced polymer (CFRP) composites. The proposed approach integrates (1) microscopic image characterization, (2) stochastic microstructure reconstruction, and (3) efficient multiscale finite element simulations enabled by self-consistent clustering (SCA) analysis. To model the complex microstructural variability, the proposed UQ methods take the non-Gaussian uncertainty sources into account through a distribution-free sampling approach leveraging nonparametric and asymptotic statistical tools. A hierarchical conditional sampling strategy enables the simultaneous sampling of multiple sources of uncertainties. Our approach provides insights into the impact of microstructural variabilities, which are shown to have an increasing impact on the nonlinear responses of UD CFRP parts under progressive compression loading and ultimately on the failure rate over time. We discover that before CFRP parts start to fail, a characteristic time period emerges with distinctive uncertainty distributions specific to the microstructure variability and the probability of failure. Identifying the failure time period is crucial to the reliability prediction, which is an essential component of CFRP design.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine Learning-Enabled Wearable Piezoelectric Acoustic Sensor for Real-Time Breast Abnormality Detection

In contemporary society, breast health has become a significant public health concern, particularly among women. According to statistics from the World Health Organization, both the incidence and mortality rates of breast tumors have steadily increased in recent years. Therefore, effective early-stage screening and postoperative monitoring are essential for maintaining breast health. However, conventional clinical diagnostic modalities are typically bulky, operationally complex, and unsuitable for continuous real-time monitoring, which limits their use in portable and everyday health management applications. To address these limitations, this study proposes a machine learning-integrated wearable piezoelectric sensing platform as an auxiliary tool for breast health assessment. The device consists of a PDMS matching layer embedded with flexible silver nanowires, a P(VDF-TrFE) piezoelectric layer, and a multi-channel low-noise signal acquisition circuit. It is capable of acquiring acoustic echo signals from tissue-mimicking environments and automatically evaluating signal validity using a convolutional neural network (CNN). By integrating piezoelectric sensing with deep learning-based signal analysis, the proposed system achieves a signal-to-noise ratio exceeding 70 dB and a real-time classification accuracy above 96% under controlled conditions. These results demonstrate that the platform provides a compact, portable, and intelligent approach for wearable sensing of mechanical heterogeneity and highlight its potential for future development in continuous biomedical monitoring technologies.

He, Shuaitong↗
DOE OSTI2023

Subseasonal Clustering of Atmospheric Rivers Over the Western United States

Abstract The serial occurrence of atmospheric rivers (ARs) along the US West Coast can lead to prolonged and exacerbated hydrologic impacts, threatening flood‐control and water‐supply infrastructure due to soil saturation and diminished recovery time between storms. Here a statistical approach for quantifying subseasonal temporal clustering among extreme events is applied to a 41‐year (1979–2019) wintertime AR catalog across the western United States (US). Observed AR occurrence, compared against a randomly distributed AR timeseries with the same average event density, reveals temporal clustering at a greater‐than‐random rate across the western US with a distinct geographical pattern. Compared to the Pacific Northwest, significant AR clusters over the northern Coastal Range of California and Sierra Nevada are more frequent and occur over longer time periods. Clusters along the California Coastal Range typically persist for 2 weeks, are composed of 4–5 ARs per cluster, and account for over 85% of total AR occurrence. Across the northwest Coast‐Cascade Ranges, clusters account for ∼50% of total AR occurrence, typically last 8–10 days, and contain 3–4 individual AR events. Based on precipitation data from a high‐resolution dynamical downscaling of reanalysis, the fractions of total and extreme hourly precipitation attributable to AR clusters are largest along the northern California coast and in the Sierra Nevada. Interannual variability among clusters highlights their importance for determining whether a particular water year is anomalously wet or dry. The mechanisms behind this unusual clustering are unclear and require further research.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Consistency testing for robust phase estimation

In this work, we present an extension to the robust phase estimation protocol, which can identify incorrect results that would otherwise lie outside the expected statistical range. Robust phase estimation is increasingly a method of choice for applications such as estimating the effective process parameters of noisy hardware, but its robustness is dependent on the noise satisfying certain threshold assumptions. We provide consistency checks that can indicate when those thresholds have been violated, which can be difficult or impossible to test directly. We test these consistency checks for several common noise models, and identify two possible checks with high accuracy in locating the point in a robust phase estimation run at which further estimates should not be trusted. One of these checks may be chosen based on resource availability, or they can be used together in order to provide additional verification.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

What drives the scatter of local star-forming galaxies in the BPT diagrams? A Machine Learning based analysis

ABSTRACT We investigate which physical properties are most predictive of the position of local star forming galaxies on the BPT diagrams, by means of different Machine Learning (ML) algorithms. Exploiting the large statistics from the Sloan Digital Sky Survey (SDSS), we define a framework in which the deviation of star-forming galaxies from their median sequence can be described in terms of the relative variations in a variety of observational parameters. We train artificial neural networks (ANN) and random forest (RF) trees to predict whether galaxies are offset above or below the sequence (via classification), and to estimate the exact magnitude of the offset itself (via regression). We find, with high significance, that parameters primarily associated to variations in the nitrogen-over-oxygen abundance ratio (N/O) are the most predictive for the [N ii]-BPT diagram, whereas properties related to star formation (like variations in SFR or EW(H α)) perform better in the [S ii]-BPT diagram. We interpret the former as a reflection of the N/O–O/H relationship for local galaxies, while the latter as primarily tracing the variation in the effective size of the S+ emitting region, which directly impacts the [S ii] emission lines. This analysis paves the way to assess to what extent the physics shaping local BPT diagrams is also responsible for the offsets seen in high redshift galaxies or, instead, whether a different framework or even different mechanisms need to be invoked.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗
DOE OSTI2023

Irradiation creep measurement and microstructural analysis of chromium nitride–coated zirconium alloy using pressurized tubes

Environmental barrier coatings for Zr-based materials are currently under development to reduce oxidation and embrittlement in light-water reactors. Chromium nitride is one such candidate for this application, particularly as accident-tolerant fuel cladding. However, quantifying the impact of coatings on the irradiation-induced creep of zircaloy (Zry) is critical as this mechanism often exceeds thermal creep rates under light-water reactor operating conditions and can be a limiting design characteristic. Additionally, examining irradiation effects in the microstructure at the coating interface is key to understanding the compatibility of the material system. Here, to accelerate the experimental measurement of irradiation creep and microstructure evolution in CrN-Zry, compact, pressurized creep tubes were fabricated and irradiated in the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory. Miniature, thin-walled rodlets fabricated from annealed Zr-Sn barstock were coated with CrN using physical vapor deposition (PVD) to nominal thicknesses of 4 and 8 μm. Coated and uncoated rodlet specimens were internally pressurized and welded, generating nominal circumferential hoop stresses of 0, 90, or 180 MPa under 300°C irradiation conditions. Twelve specimens were measured diametrically prior to irradiation using a low-cost, automated, contactless laser profilometer developed for this work. Specimens were irradiated in sealed capsules for one 25-day HFIR cycle, accumulating approximately 1.8 $\times$ 10 21 n/cm 2 fast fluence (E n > 10 MeV). The irradiated samples were retrieved and remeasured using the same profilometry system in a shielded hot cell facility. Irradiation creep between specimens was compared using standard statistical tests and showed that both thicknesses of CrN coating had a negligible effect on the irradiation creep strain of the Zry material. Microstructure characterization of pre- and post-irradiated CrN-Zry specimens showed minimal changes due to irradiation but did show a substantial O-rich region at the Zry-CrN interface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modular Processing of Flare Gas for Carbon Nanoproducts

This project demonstrated the technical viability and economic promise of a modular system for converting flared natural gas into valuable carbon nanoproducts (CNPs) through catalytic chemical vapor deposition (CVD). All major project milestones were successfully completed, including reactor design and commissioning, catalyst development, process optimization, technoeconomic analysis, and application testing in concrete systems. The overarching goal was to create a scalable, field-deployable process that valorizes stranded methane by producing high-value carbon materials for use in cementitious composites. At the lab scale, the team designed and built a fluidized bed reactor optimized for use with silica fume-supported nickel catalysts synthesized via atomic layer deposition (ALD). A statistically designed sintering study enabled precise tuning of nickel nanoparticle size, identifying the influence of oxygen partial pressure, time, and temperature on catalyst morphology and performance. These insights allowed the team to target catalyst conditions that maximize carbon nanofilament growth. Subsequent CVD experiments achieved up to 31.8 wt% carbon deposition under optimized conditions, with TEM confirming the presence of nanofilament structures and sustained hydrogen evolution during reaction. Reactor upgrades and empirical fluidization studies supported the development of reliable, repeatable experimental protocols. The modular pilot-scale skid reactor was fully constructed, instrumented, and commissioned. Capable of operating at 675–800°C and pressures up to 290 psig, the system was designed for continuous operation at a carbon production rate of 1 kg/hr. Initial demonstration runs confirmed solids handling, thermal control, and system leak-tightness, although a critical reactor component (the downfeed tube) was inadvertently omitted during final assembly. This omission limited gas–solid contact and prevented meaningful carbon deposition during pilot-scale CVD runs. Nonetheless, the system operated safely under design conditions, and the root cause of performance limitations was clearly identified. Complementary work on UHPC formulations demonstrated that small additions of carbon nanoproducts, including those derived from flare gas, can significantly enhance mechanical performance while preserving workability. A comprehensive study of CNF dispersion techniques and mix design optimization led to a clear protocol for integrating these nanomaterials into concrete. Incorporation of CNPs improved flexural toughness and reduced porosity, supporting their use in high-performance infrastructure applications. A technoeconomic analysis (TEA) confirmed that this process can produce CNP-loaded catalyst material at a levelized cost below $\$$7/kg across a range of catalyst loadings and reaction yields. With estimated market values for the carbon composite product ranging from $\$$14 to over $\$$60/kg, and the ability to blend CNPs into concrete at sub-percent levels with less than 10% added cost, the system presents a compelling economic case. While additional engineering work is needed to optimize fluidization and heat transfer at scale, this project establishes a strong foundation for commercial development. The process is not only technically sound but also economically promising, representing a viable pathway for flare gas mitigation through modular carbon nanomaterial production.

03 NATURAL GAS↗
DOE OSTI2023

Impacts of divalent cations (Mg 2+ and Ca 2+ ) on PFAS bioaccumulation in freshwater macroinvertebrates representing different foraging modes

Per- and polyfluoroalkyl substances (PFAS) have extensively contaminated freshwater aquatic ecosystems where they can be transported in water and partition to sediment and biota. In this paper, three freshwater benthic macroinvertebrates with different foraging modes were exposed to environmentally relevant concentrations of eight perfluoroalkyl carboxylates (PFCA), three perfluoroalkyl sulfonates (PFSA), and three fluorotelomer sulfonates (FTS) at varying divalent cation concentrations of magnesium (Mg 2+ ) and calcium (Ca 2+ ). Divalent cations can impact PFAS partitioning to solids, especially to sediments, at higher concentrations. Sediment dwelling worms (Lumbriculus variegatus), epibenthic grazing snails (Physella acuta), and sediment-dwelling filter-feeding bivalves (Elliptio complanata) were selected due to their unique foraging modes. Microcosms were composed of synthetic sediment, culture water, macroinvertebrates, and PFAS and consisted of a 28-day exposure period. L. variegatus had significantly higher PFAS bioaccumulation than P. acuta and E. complanata, likely due to higher levels of interactions with and ingestion of the contaminated sediment. “High Mg 2+ ” (7.5 mM Mg 2+ ) and “High Ca 2+ ” (7.5 mM Ca 2+ ) conditions generally had statistically higher bioaccumulation factors (BAF) than the “Reference Condition” (0.2 mM Ca 2+ and 0.2 mM Mg 2+ ) for PFAS with perfluorinated chain lengths greater than eight carbons. Long-chain PFAS dominated the PFAS profiles of the macroinvertebrates for all groups of compounds studied (PFCA, PFSA, and FTS). Furthermore, these results indicate that the study organism has the greatest impact on bioaccumulation, although divalent cation concentration had observable impacts between organisms depending on the environmental conditions. Elevated cation concentrations in the microcosms led to significantly greater bioaccumulation in the test organisms compared to the experimental reference conditions for long-chain PFAS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Fast & accurate emulation of two-body scattering observables without wave functions

We combine Newton's variational method with ideas from eigenvector continuation to construct a fast & accurate emulator for two-body scattering observables. The emulator will facilitate the application of rigorous statistical methods for interactions that depend smoothly on a set of free parameters. Our approach begins with a trial K or T matrix constructed from a small number of exact solutions to the Lippmann–Schwinger equation. Subsequent emulation only requires operations on small matrices. We provide several applications to short-range potentials with and without the Coulomb interaction and partial-wave coupling. It is shown that the emulator can accurately extrapolate far from the support of the training data. When used to emulate the neutron-proton cross section with a modern chiral interaction as a function of 26 free parameters, it reproduces the exact calculation with negligible error and provides an over 300x improvement in CPU time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Broadband quadrature-squeezed vacuum and nonclassical photon number correlations from a nanophotonic device

We report demonstrations of both quadrature-squeezed vacuum and photon number difference squeezing generated in an integrated nanophotonic device. Squeezed light is generated via strongly driven spontaneous four-wave mixing below threshold in silicon nitride microring resonators. The generated light is characterized with both homodyne detection and direct measurements of photon statistics using photon number–resolving transition-edge sensors. We measure 1.0(1) decibels of broadband quadrature squeezing (~4 decibels inferred on-chip) and 1.5(3) decibels of photon number difference squeezing (~7 decibels inferred on-chip). Nearly single temporal mode operation is achieved, with measured raw unheralded second-order correlations g (2) as high as 1.95(1). Multiphoton events of over 10 photons are directly detected with rates exceeding any previous quantum optical demonstration using integrated nanophotonics. These results will have an enabling impact on scaling continuous variable quantum technology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

The equilibrium configuration of a plasma in an axially symmetric reactor is described mathematically by a free boundary problem associated with the celebrated Grad-Shafranov equation. The presence of uncertainty in the model parameters introduces the need to quantify the variability in the predictions. This is often done by computing a large number of model solutions on a computational grid for an ensemble of parameter values and then obtaining estimates for the statistical properties of solutions. In this study, we explore the savings that can be obtained using multilevel Monte Carlo methods, which reduce costs by performing the bulk of the computations on a sequence of spatial grids that are coarser than the one that would typically be used for a simple Monte Carlo simulation. We examine this approach using both a set of uniformly refined grids and a set of adaptively refined grids guided by a discrete error estimator. Numerical experiments show that multilevel methods dramatically reduce the cost of simulation, with cost reductions typically on the order of 60 or more and possibly as large as 200. Furthermore, adaptive griding results in more accurate computation of geometric quantities such as x-points associated with the model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗