Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scalar field 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 127 records · Page 7

The Lorentzian inversion formula and the spectrum of the 3d O(2) CFT

We study the spectrum and OPE coefficients of the three-dimensional critical O(2) model, using four-point functions of the leading scalars with charges 0, 1, and 2 ( s , $\phi$, and t ). We obtain numerical predictions for low-twist OPE data in several charge sectors using the extremal functional method. We compare the results to analytical estimates using the Lorentzian inversion formula and a small amount of numerical input. We find agreement between the analytic and numerical predictions. We also give evidence that certain scalar operators lie on double-twist Regge trajectories and obtain estimates for the leading Regge intercepts of the O(2) model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The cosmological evolution of self-interacting dark matter

We study the evolution of cosmological perturbations in dark-matter models with elastic and velocity-independent self interactions. Such interactions are imprinted in the matter-power spectrum as dark acoustic oscillations, which can be experimentally explored to determine the strength of the self scatterings. Models with self interactions have similarities to warm dark matter, as they lead to suppression of power on small scales when the dark-matter velocity dispersion is sizable. Nonetheless, both the physical origin and the extent of the suppression differ for self-interacting dark matter from conventional warm dark matter, with a dark sound horizon controlling the reduction of power in the former case, and a free-streaming length in the latter. We thoroughly analyze these differences by performing computations of the linear power spectrum using a newly developed Boltzmann code. We find that while current Lyman-α data disfavor conventional warm dark matter with a mass less than 5.3 keV, when self interactions are included at their maximal value consistent with bounds from the Bullet Cluster, the limits are relaxed to 4.4 keV. Lastly, we make use of our analysis to set novel bounds on light scalar singlet dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Light fields during inflation from BOSS and future galaxy surveys

Abstract Primordial non-Gaussianity generated by additional fields present during inflation offers a compelling observational target for galaxy surveys. These fields are of significant theoretical interest since they offer a window into particle physics in the inflaton sector. They also violate the single-field consistency conditions and induce a scale-dependent bias in the galaxy power spectrum. In this paper, we explore this particular signal for light scalar fields and study the prospects for measuring it with galaxy surveys. We find that the sensitivities of current and future surveys are remarkably stable for different configurations, including between spectroscopic and photometric redshift measurements. This is even the case at non-zero masses where the signal is not obviously localized on large scales. For realistic galaxy number densities, we demonstrate that the redshift range and galaxy bias of the sample have the largest impact on the sensitivity in the power spectrum. These results additionally motivated us to explore the potentially enhanced sensitivity of Vera Rubin Observatory's LSST through multi-tracer analyses. Finally, we apply this understanding to current data from the last data release of the Baryon Oscillation Spectroscopic Survey (BOSS DR12) and place new constraints on light fields coupled to the inflaton.

Astronomy & Astrophysics↗

New Machine Learning Techniques for Simulation-Based Inference: InferoStatic Nets, Kernel Score Estimation, and Kernel Likelihood Ratio Estimation

We propose an intuitive, machine-learning approach to multiparameter inference, dubbed the InferoStatic Networks (ISN) method, to model the score and likelihood ratio estimators in cases when the probability density can be sampled but not computed directly. The ISN uses a backend neural network that models a scalar function called the inferostatic potential $\varphi$. In addition, we introduce new strategies, respectively called Kernel Score Estimation (KSE) and Kernel Likelihood Ratio Estimation (KLRE), to learn the score and the likelihood ratio functions from simulated data. We illustrate the new techniques with some toy examples and compare to existing approaches in the literature. We mention en passant some new loss functions that optimally incorporate latent information from simulations into the training procedure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for New Physics through a Long-Lived Di-Muon Resonance in the NuMI Beam with the ICARUS Detector

Many of the unexplained phenomena in particle physics and cosmology today, such as the microphysical nature of dark matter, the strong CP problem, and the origin of the neutrino masses, can be resolved by the existence of a light (~GeV), weakly-coupled hidden sector of new physics. Such hidden sectors often predict the existence of “long-lived” particles (LLPs) that travel a far distance from production before decaying into Standard Model particles. Neutrino oscillation experiments, which combine intense particle beams with precise imaging detectors, are well equipped to probe LLP models with new sensitivity. This thesis details a search for a long-lived particle decaying to two muons with the ICARUS liquid argon time projection chamber (LArTPC) neutrino detector in the Short-Baseline Neutrino program at Fermilab. The calibration of the ICARUS time projection chamber (TPC) which enables the search is also presented. Notably, the calibration measures an angular dependence in electron-ion recombination in argon, a novel effect in the detector physics of LArTPCs. The search is performed using data taken with the Neutrinos at the Main Injector (NuMI) beam, with an exposure of 2.41e20 protons on target. No significant excess over background is observed, and we set world-leading limits on two new physics models that predict this process: the Higgs portal scalar and a heavy axion model. We also present the sensitivity in a model-independent way applicable to any new physics model predicting the process K → π + S(→μμ), for a long-lived particle S.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Geometry-Driven Detection, Tracking and Visual Analysis of Viscous and Gravitational Fingers

Viscous and gravitational flow instabilities cause a displacement front to break up into finger-like fluids. The detection and evolutionary analysis of these fingering instabilities are critical in multiple scientific disciplines such as fluid mechanics and hydrogeology. However, previous detection methods of the viscous and gravitational fingers are based on density thresholding, which provides limited geometric information of the fingers. The geometric structures of fingers and their evolution are important yet little studied in the literature. In this work, we explore the geometric detection and evolution of the fingers in detail to elucidate the dynamics of the instability. Additionally, we propose a ridge voxel detection method to guide the extraction of finger cores from three-dimensional (3D) scalar fields. After skeletonizing finger cores into skeletons, we design a spanning tree based approach to capture how fingers branch spatially from the finger skeletons. Finally, we devise a novel geometric-glyph augmented tracking graph to study how the fingers and their branches grow, merge, and split over time. Feedback from earth scientists demonstrates the usefulness of our approach to performing spatio-temporal geometric analyses of fingers.

97 MATHEMATICS AND COMPUTING↗

Minimal warm inflation

Slow-roll inflation is a successful paradigm. However we find that even a small coupling of the inflaton to other light fields can dramatically alter the dynamics and predictions of inflation. As an example, the inflaton can generically have an axion-like coupling to gauge bosons. Even relatively small couplings will automatically induce a thermal bath during inflation. The thermal friction from this bath can easily be stronger than Hubble friction, significantly altering the usual predictions of any particular inflaton potential. Thermal effects suppress the tensor-to-scalar ratio r significantly, and predict unique non-gaussianities. This axion-like coupling provides a minimal model of warm inflation which avoids the usual problem of thermal backreaction on the inflaton potential. As a specific example, we find that hybrid inflation with this axion-like coupling can easily fit the current cosmological data.

79 ASTRONOMY AND ASTROPHYSICS↗

The second Higgs at the lifetime frontier

We assess the current coverage and the future discovery potential of LHC searches for heavy Higgs bosons decaying into long-lived particles (LLPs), focusing primarily on the production of pairs of LLPs with hadronic final states. These signatures are generic in dark sectors where a heavy scalar decays into pairs of lighter states which subsequently mix with the Standard Model Higgs. We show that a handful of existing analyses provide broad coverage of LLP decay lengths ranging from millimeters to tens of meters, and explore the complementarity between searches for displaced and prompt final states in several simplified models. For both heavy singlet and heavy doublet scalars, LLP searches typically provide the leading sensitivity in current data and exhibit the strongest discovery potential in future LHC runs. We further translate the impact of these searches into the parameter space of various Twin Higgs models, demonstrating that LLP searches are a promising avenue for discovering a Twin Higgs with displaced decays. Finally, we propose a variety of additional search channels that would improve coverage of the second Higgs at the lifetime frontier.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor mixing induced by power-law accelerations in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing turbulent Rayleigh–Taylor mixing driven by a temporal power-law acceleration are derived in the small Atwood number (Boussinesq) limit. The solutions generalize those previously derived for constant acceleration Rayleigh–Taylor mixing for models based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ [O. Schilling, “Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit,” Phys. Fluids 33, 085129 (2021)]. The turbulent fields are expressed in terms of the model coefficients and power-law exponent, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. Mixing layer growth parameters and other physical observables are obtained explicitly as functions of the model coefficients and parameterized by the exponent of the power-law acceleration. Values for physical observables in the constant acceleration case are used to calibrate the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to a canonical (i.e., constant acceleration) Rayleigh–Taylor turbulent flow. The calibrated four-equation model is then used to numerically reconstruct the mean and turbulent fields, and turbulent equation budgets across the mixing layer for several values of the power-law exponent. Finally, the reference solutions derived here can be used to understand the model predictions for strongly accelerated or decelerated Rayleigh–Taylor mixing in the large Reynolds number limit.

42 ENGINEERING↗

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗

On the decaying-sterile-neutrino solution to the electron (anti)neutrino appearance anomalies

We explore the hypothesis that the unexplained data from LSND and MiniBooNE are evidence for a new, heavy neutrino mass-eigenstate that mixes with the muon-type neutrino and decays into an electron-type neutrino and a new, very light scalar particle. We consider two different decay scenarios, one with Majorana neutrinos, one with Dirac neutrinos; both fit the data equally well. We find a reasonable, albeit not excellent, fit to the data of MiniBooNE and LSND. The decaying-sterile-neutrino hypothesis, however, cleanly evades constraints from disappearance searches and precision measurements of leptonic meson decays, as long as $1~{\rm MeV}\gtrsim m_4\gtrsim 10$~keV. The SBN program at Fermilab should be able to definitively test the decaying-sterile-neutrino hypothesis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Models of $J/\varPsi$ photo-production reactions on the nucleon

The J/Ψ photo-production reactions on the nucleon can provide information on the roles of gluons in determining the J/Ψ-nucleon (J/Ψ-N) interactions and the structure of the nucleon. The information on the J/Ψ-N interactions is needed to test lattice QCD (LQCD) calculations and to understand the nucleon resonances such as N*(P c ) recently reported by the LHCb Collaboration. In addition, it is also needed to investigate the production of nuclei with hidden charms and to extract the gluon distributions in nuclei. The main purpose of this article is to review six models of the γ+p→J/Ψ+p reaction which have been and can be applied to analyze the data from Thomas Jefferson National Accelerator Facility (JLab). Here, the formulae for each model are given and used to obtain the results to show the extent to which the available data can be described. The models presented include the Pomeron-exchange model of Donnachie and Landshoff (Pom-DL) and its extensions to include J/Ψ-N potentials extracted from LQCD (Pom-pot) and to also use the constituent quark model (CQM) to account for the quark substructure of J/Ψ (Pom-CQM). The other three models are developed from applying the perturbative QCD approach to calculate the two-gluon exchange using the generalized parton distribution (GPD) of the nucleon (GPD-based), two- and three-gluon exchanges using the parton distribution of the nucleon (2g+3g), and the exchanges of scalar (0 ++ ) and tensor (2 ++ ) glueballs within the holographic formulation (holog). The results of investigating the excitation of the nucleon resonances N*(P c ) in the γ+p→J/Ψ+p reactions are also given. We demonstrate that the differences between these six models can be unambiguously distinguished and the N* can be better studied by using the forthcoming JLab data at large |t| and at energies very near the J/Ψ production threshold. Possible improvements of the considered models are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

58 GEOSCIENCES↗

Spatially Accelerated Winding Numbers for Curved Geometry

The generalized winding number (GWN) is a scalar field that supports robust containment queries on curved geometry, including non-watertight, overlapping, and nested boundary representations. While queries can be easily parallelized over samples, direct evaluation on parametric curves and surfaces remains costly for large and complex models. Fast, state-of-the-art GWN approaches leverage a spatial index to approximate the GWN, typically coupled with a Taylor expansion which approximates the GWN contribution for far clusters of geometric primitives. However, such methods operate only on discrete inputs such as triangle meshes and point clouds, and would introduce containment errors near boundaries if applied to curved input. We extend support for fast GWN evaluation over arbitrary collections of NURBS curves in 2D and trimmed NURBS patches in 3D via a Bounding Volume Hierarchy that stores efficiently precomputed moment data in the hierarchy nodes. When querying the hierarchy, approximations for far clusters are used alongside direct evaluation for nearby NURBS primitives, achieving sub-linear complexity while preserving the geometric features in the vicinity of the query point. Central to our performance improvements is an adaptive subdivision strategy for NURBS primitives during a preprocessing phase, creating better spatial partitions while retaining the same accuracy for containment decisions as a direct evaluation. We demonstrate the performance and accuracy of our approach across a large collection of 2D and 3D datasets.

Computer science↗

Mapping structural heterogeneity at the nanoscale with scanning nano-structure electron microscopy (SNEM)

Here, in this work, we explore the use of scanning electron diffraction (also known as 4D-STEM) coupled with electron atomic pair distribution function analysis (ePDF) to understand the local order (structure and chemistry) as a function of position in a complex multicomponent system, a hot rolled, Ni-encapsulated, Zr 65 Cu 17.5 Ni 10 Al 7.5 bulk metallic glass (BMG), with a spatial resolution of 3 nm. We show that it is possible to gain insight into the chemistry and chemical clustering/ordering tendency in different regions of the sample, including in the vicinity of nano-scale crystallites that are identified from virtual dark field images and in heavily deformed regions at the edge of the BMG. In addition to simpler analysis, unsupervised machine learning was used to extract partial PDFs from the material, modeled as a quasi-binary alloy, and map them in space. These maps allowed key insights not only into the local average composition, as validated by EELS, but also a unique insight into chemical short-range ordering tendencies in different regions of the sample during formation. The experiments are straightforward and rapid and, unlike spectroscopic measurements, don’t require energy filters on the instrument. We spatially map different quantities of interest (QoI’s), defined as scalars that can be computed directly from positions and widths of ePDF peaks or parameters refined from fits to the patterns. We developed a flexible and rapid data reduction and analysis software framework that allows experimenters to rapidly explore images of the sample on the basis of different QoI’s. The power and flexibility of this approach are explored and described in detail. Because of the fact that we are getting spatially resolved images of the nanoscale structure obtained from ePDFs we call this approach scanning nano-structure electron microscopy (SNEM), and we believe that it will be powerful and useful extension of current 4D-STEM methods.

36 MATERIALS SCIENCE↗

Cosmic topology. Part IIa. Eigenmodes, correlation matrices, and detectability of orientable Euclidean manifolds

Abstract If the Universe has non-trivial spatial topology, observables depend on both the parameters of the spatial manifold and the position and orientation of the observer. In infinite Euclidean space, most cosmological observables arise from the amplitudes of Fourier modes of primordial scalar curvature perturbations. Topological boundary conditions replace the full set of Fourier modes with specific linear combinations of selected Fourier modes as the eigenmodes of the scalar Laplacian. We present formulas for eigenmodes in orientable Euclidean manifolds with the topologiesE 1 –E 6 ,E 11 ,E 12 ,E 16 , andE 18 that encompass the full range of manifold parameters and observer positions, generalizing previous treatments. Under the assumption that the amplitudes of primordial scalar curvature eigenmodes are independent random variables, for each topology we obtain the correlation matrices of Fourier-mode amplitudes (of scalar fields linearly related to the scalar curvature) and the correlation matrices of spherical-harmonic coefficients of such fields sampled on a sphere, such as the temperature of the cosmic microwave background (CMB). We evaluate the detectability of these correlations given the cosmic variance of the observed CMB sky. We find that topologies where the distance to our nearest clone is less than about 1.2 times the diameter of the last scattering surface of the CMB give a correlation signal that is larger than cosmic variance noise in the CMB. This implies that if cosmic topology is the explanation of large-angle anomalies in the CMB, then the distance to our nearest clone is not much larger than the diameter of the last scattering surface. We argue that the topological information is likely to be better preserved in three-dimensional data, such as will eventually be available from large-scale structure surveys.

Astronomy & Astrophysics↗

SPT-3G D1: Constraints on inflationary gravitational waves with two years of SPT-3G data

Here, we present a measurement of the 𝐵-mode polarization power spectrum of the cosmic microwave background anisotropies at 32 ≤ ℓ < 502 for three bands centered at 95, 150, and 220 GHz using data from the SPT-3G receiver on the South Pole Telescope. This work uses SPT-3G observations from the 2019 and 2020 winter observing seasons of a ∼1500 deg 2 patch of sky that directly overlaps with fields observed with the BICEP/Keck family of telescopes and covers part of the proposed Simons Observatory and CMB-S4 deep fields. Employing new techniques for mitigating polarized atmospheric noise, the SPT-3G data demonstrates a white noise level of 9.3 (6.7) μ⁢K−arcmin at ℓ ∼500 for the 95 GHz (150 GHz) data, with a 1/ℓ noise knee at ℓ = 128 (182). We fit the observed six auto- and cross-frequency 𝐵-mode power spectra to a model including lensed Λ⁢CDM 𝐵-modes and a combination of Galactic and extragalactic foregrounds. This work characterizes foregrounds in the vicinity of the BICEP/Keck survey area, finding foreground power consistent with that reported by the BICEP/Keck collaboration within the same region and a factor of ∼3 higher power over the full SPT-3G survey area. Using SPT-3G data over the BICEP/Keck survey area, we place a 95% upper limit on the tensor-to-scalar ratio of 𝑟 <0.25 and find the statistical uncertainty on 𝑟 to be 𝜎⁡(𝑟) = 0.067.

Zebrowski, J. A. [University of Chicago; Universit↗

Composition dependence of bulk properties in the Co-intercalated transition metal dichalcogenide $\mathrm{Co_{1/3}TaS_2}$

Spontaneous Hall conductivity has recently been reported in the triangular lattice antiferromagnet Co 1/3 TaS 2 under a zero magnetic field. This phenomenon originates from the distinctive noncoplanar triple-Q magnetic ground state, possessing uniform real-space Berry curvature characterized by scalar spin chirality. We investigated the physical properties of Co 1/3 TaS 2 by judiciously controlling the composition, revealing a drastic change in its bulk properties, even by slight variations in cobalt composition, despite the same crystal structure. For 0.299≤x≤0.325, Co x TaS 2 keeps all the characteristics of the ground state consistent with the previous studies—two antiferromagnetic phase transitions at T N1 and T N2 ( N1 ), a large spontaneous Hall conductivity [σ xy (H=0)], and a weak ferromagnetic moment along the c axis. However, samples with x≥0.330 exhibit distinct bulk properties, including the absence of both σ xy (H=0) and the weak ferromagnetic moment. Our neutron diffraction data reveal that Co x TaS 2 with x≥0.330 develops coplanar helical magnetic order with q m1 =(1/3, 0, 0). This is entirely different from what has been seen in x≤0.325, explaining the observed composition dependence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗