Results for “Sigma Models” — Page 6

DOE OSTI2025

Proton spin from small-x with constraints from the valence quark model

We apply the valence quark model [1] to constrain the non-perturbative initial condition for the small-x helicity evolution. The remaining free parameters are constrained by performing a global analysis akin to [2] to the available polarized small-x deep inelastic scattering data. A good description of the world data is obtained with only 8 free parameters. The model parameters are tightly constrained by the data, allowing us to predict the proton polarized structure function $g^p_1$ to be negative at small x. Furthermore, we obtain the small-x quark and gluon spins to give a contribution $\int^{0.1}_{10^{-5}} dx (\frac{1}{2}\Delta\Sigma + \Delta{G}) = 0.63 ± 0.10$ or $1.35 ± 0.16$ (or to the proton spin, depending on the applied running coupling prescription.

Color-class condensate↗

DOE OSTI2026

Improving MicroBooNE's Inclusive Single-Photon Search with Low-Energy Hadronic Identification

This analysis aims to further investigate MicroBooNE’s inclusive single-photon search results, which reported a 2.2$\sigma$ excess below 600 MeV in shower energy for events with no reconstructed protons using roughly half of MicroBooNE's dataset. Taken together with MiniBooNE’s long-standing low-energy excess and MicroBooNE’s recent electron-like search results showing no observable excess with respect to Standard Model predictions, this result provides strong impetus for expanded exploration of the single-photon channel in Fermilab’s short-baseline liquid-argon time projection chamber (LArTPC) experiments. By identifying and classifying isolated MeV-scale energy depositions, or blips, in the vicinity of single-photon events selected by the Wire-Cell reconstruction framework, we establish a more comprehensive labeling scheme for nearby hadronic content. In particular, blips found backwards along the shower axis indicate the presence of previously-unidentified final-state protons, while elevated blip counts at wide angles signal the presence of final-state neutrons. By applying this new technique to its full dataset, MicroBooNE will perform a purer and higher-statistics test of the truly isolated nature of its modest photon-like excess, furthering its hunt for the presence of unexpected new physics beyond the Standard Model.

Andrade Aldana, Diego [Los Alamos; IIT, Chicago (m↗

DOE OSTI2024

Cosmological Constraints from Combining Photometric Galaxy Surveys and Gravitational Wave Observatories

Spatial variations in survey properties due to selection effects generate substantial systematic errors in large-scale structure measurements in optical galaxy surveys on very large scales. On such scales, the statistical sensitivity of optical surveys is also limited by their finite sky coverage. By contrast, gravitational wave (GW) sources appear to be relatively free of these issues, provided the angular sensitivity of GW experiments can be accurately characterized. We quantify the expected cosmological information gain from combining the forecast LSST 3$\times$2pt analysis (combination of three 2-point correlations of galaxy density and weak lensing shear fields) with the large-scale auto-correlation of GW sources from proposed next-generation GW experiments. We find that in $\Lambda$CDM and $w$CDM models, there is no significant improvement in cosmological constraints from combining GW with LSST 3$\times$2pt over LSST alone, due to the large shot noise for the former; however, this combination does enable a $\sim6\%$ constraint on the linear galaxy bias of GW sources. More interestingly, the optical-GW data combination provides tight constraints on models with primordial non-Gaussianity (PNG), due to the predicted scale-dependent bias in PNG models on large scales. Assuming that the largest angular scales that LSST will probe are comparable to those in Stage III surveys ($\ell_{\rm min}\sim50$), the inclusion of next-generation GW measurements could improve constraints on the PNG parameter $f_{\rm NL}$ by up to a factor of $\simeq6.6$ compared to LSST alone, yielding $\sigma(f_{\rm NL})=8.5$. These results assume the expected capability of a network of Einstein Telescope-like GW observatories, with a detection rate of $10^6$ events/year. We investigate the sensitivity of our results to different assumptions about future GW detectors as well as different LSST analysis choices.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and an efficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad midwall temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. There are many development shortcomings in DASSH detailed in this document, but this version is functional for modern analysis needs. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology for pin lattice models is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. The user can alternatively provide their own geometry and power profile instead of inheriting it from DIF3D. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and a sufficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Observation of $\gamma\gamma\to\tau\tau$ in proton-proton collisions and limits on the anomalous electromagnetic moments of the $\tau$ lepton

The production of a pair of τ leptons via photon–photon fusion, ${{\unicode{x03B3}\unicode{x03B3}\to\unicode{x03C4}\unicode{x03C4}}}$, is observed for the first time in proton–proton collisions, with a significance of 5.3 standard deviations. This observation is based on a data set recorded with the CMS detector at the LHC at a center-of-mass energy of 13 TeV and corresponding to an integrated luminosity of 138 fb$^{−1}$. Events with a pair of τ leptons produced via photon–photon fusion are selected by requiring them to be back-to-back in the azimuthal direction and to have a minimum number of charged hadrons associated with their production vertex. The τ leptons are reconstructed in their leptonic and hadronic decay modes. The measured fiducial cross section of ${{\unicode{x03B3}\unicode{x03B3}\to\unicode{x03C4}\unicode{x03C4}}}$ is $\sigma^\text{fid}_\text{obs} = 12.4 ^{+3.8}_{-3.1}\,\textrm{fb}$. Constraints are set on the contributions to the anomalous magnetic moment (${a_{\unicode{x03C4}}}$) and electric dipole moments (${d_{\unicode{x03C4}}}$) of the τ lepton originating from potential effects of new physics on the $\unicode{x03B3}\unicode{x03C4}\unicode{x03C4}$ vertex: ${a_{\unicode{x03C4}}} = 0.0009_{-0.0031}^{+0.0032}$ and $|{d_{\unicode{x03C4}}}| \lt 2.9\times 10^{-17}\,{e}\,\textrm{cm}$ (95% confidence level), consistent with the standard model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Measurement of the t-channel single top quark cross section in proton-proton collisions at $\sqrt{s}=5.02$ TeV

The single top quark t-channel production cross section is measured in proton-proton collisions at the CERN LHC at $\sqrt{s}=5.02$ TeV, using data recorded with the CMS detector in 2017, corresponding to an integrated luminosity of 302 pb −1 , and resulting in the first CMS measurement of the process at that energy. Events with one electron or muon and two or more jets, among which at least one is identified as originating from a b quark fragmentation, are analyzed. The combined cross section of single top quark (tq) and single top antiquark ($\overline{\mathrm{t}}\mathrm{q}$) production is $\sigma \left(\mathrm{tq}+\overline{\mathrm{t}}\mathrm{q}\right)={25.5}_{-3.5}^{+3.6}{\left(\mathrm{stat}\right)}_{-3.9}^{+4.2}\left(\mathrm{syst}\right)\pm 0.5\left(\mathrm{lumi}\right)$ pb. The individual cross sections are measured to be $\sigma \left(\mathrm{tq}\right)={17.7}_{-2.7}^{+2.9}{\left(\mathrm{stat}\right)}_{-2.4}^{+2.6}\left(\mathrm{syst}\right)\pm 0.3\left(\mathrm{lumi}\right)$ pb and $\sigma \left(\overline{\mathrm{t}}\mathrm{q}\right)={6.7}_{-1.6}^{+2.4}{\left(\mathrm{stat}\right)}_{-2.5}^{+2.1}\left(\mathrm{syst}\right)\pm 0.1\left(\mathrm{lumi}\right)$ pb. Their ratio is measured to be ${\mathcal{R}}_{\mathrm{t}-\mathrm{ch}}={2.7}_{-0.8}^{+1.5}{\left(\mathrm{stat}\right)}_{-0.3}^{+1.2}\left(\mathrm{syst}\right)$. The absolute value of the Cabibbo-Kobayashi-Maskawa matrix element is found to be |f LV V tb | = 0.92 ± 0.09(exp) ± 0.01(theo). The measurements are in good agreement with the standard model predictions at next-to-next-to-leading order accuracy in quantum chromodynamics.

Hadron-Hadron Scattering↗

DOE OSTI2025

Clustering of DESI galaxies split by thermal Sunyaev-Zeldovich effect

The thermal Sunyaev-Zeldovich (tSZ) effect is associated with galaxy clusters - extremely large and dense structures tracing the dark matter with a higher bias than isolated galaxies. We propose to use the tSZ data to separate galaxies from redshift surveys into distinct subpopulations corresponding to different densities and biases independently of the redshift survey systematics. Leveraging the information from different environments, as in density-split and density-marked clustering, is known to tighten the constraints on cosmological parameters, like $\Omega_m$, $\sigma_8$ and neutrino mass. We use data from the Dark Energy Spectroscopic Instrument (DESI) and the Atacama Cosmology Telescope (ACT) in their region of overlap to demonstrate informative tSZ splitting of Luminous Red Galaxies (LRGs). We discover a significant increase in the large-scale clustering of DESI LRGs corresponding to detections starting from 1-2 sigma in the ACT DR6 + Planck tSZ Compton-$y$ map, below the cluster candidate threshold (4 sigma). We also find that such galaxies have higher line-of-sight coordinate (and velocity) dispersions and a higher number of close neighbors than both the full sample and near-zero tSZ regions. We produce simple simulations of tSZ maps that are intrinsically consistent with galaxy catalogs and do not include systematic effects, and find a similar pattern of large-scale clustering enhancement with tSZ effect significance. Moreover, we observe that this relative bias pattern remains largely unchanged with variations in the galaxy-halo connection model in our simulations. This is promising for future cosmological inference from tSZ-split clustering with semi-analytical models. Thus, we demonstrate that valuable cosmological information is present in the lower signal-to-noise regions of the thermal Sunyaev-Zeldovich map, extending far beyond the individual cluster candidates.

Astronomy data analysis↗

DOE OSTI2024

Suppression of diffraction in deep-inelastic scattering on nuclei and dynamical mechanism of leading twist nuclear shadowing

Abstract Using the leading twist approach (LTA) to nuclear shadowing, we calculate the ratios of diffractive and usual parton distributions for a heavy nucleus (Pb) and the proton,$$ {R}_{A/p}=\left({f}_{i/A}^{D(3)}/{f}_{i/A}\right)/\left({f}_{i/p}^{D(3)}/{f}_{i/p}\right) $$ R A / p = f i / A D 3 / f i / A / f i / p D 3 / f i / p , for coherent and summed (coherent plus quasi-elastic) nuclear deep-inelastic scattering. We find thatR A/p ≈ 0.5 − 1 for quarks as well as for the ratio of the diffractive and total cross sections$$ {\left[\left({d\sigma}_{\textrm{diff}}/{d M}_X^2\right)/{\sigma}_{\textrm{tot}}\right]}_{eA}/{\left[\left({d\sigma}_{\textrm{diff}}/{d M}_X^2\right)/{\sigma}_{\textrm{tot}}\right]}_{ep} $$ dσ diff / dM X 2 / σ tot eA / dσ diff / dM X 2 / σ tot ep andR A/p ≈ 0.5 − 1.3 for gluons in a broad range ofx, including the kinematics of the Electron-Ion Collider, which reaffirms the difference from the nuclear enhancement ofR A/p predicted in the gluon saturation framework. We demonstrate that the magnitude ofR A/p is controlled by the cross section of the interaction of hadronic fluctuations of the virtual photon with target nucleons, which explains an enhancement ofR A/p in the color dipole model and its suppression in LTA. We argue that the black disk limit corresponds toR A/p = 1 and$$ {R}_{A/p}^{\textrm{coh}} $$ R A / p coh = 0.86 for the summed and coherent scattering, respectively. Relying on an intuitive definition of the saturation scale, we show that the ratio of the saturation scales of a heavy nucleus and proton$$ {Q}_{sA}^2(b)/{Q}_{sp}^2(b)\approx 1 $$ Q sA 2 b / Q sp 2 b ≈ 1 at small impact parametersbdue to the strong leading twist nuclear shadowing and diluteness of the nuclear density.

Physics↗

DOE OSTI2025

Constraints on Inflationary Gravitational Waves with Two Years of SPT-3G Data

We present a measurement of the $B$-mode polarization power spectrum of the cosmic microwave background anisotropies at 32 $\le$$\ell$$<$ 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 $\sim$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) $\mu$K-arcmin at $\ell \sim 500$ for the 95 GHz (150 GHz) data, with a $1/\ell$ noise knee at $\ell$=128 (182). We fit the observed six auto- and cross-frequency $B$-mode power spectra to a model including lensed $\Lambda$CDM $B$-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 $\sim$ 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 $r < 0.25$ and find the statistical uncertainty on $r$ to be $\sigma(r) = 0.067$.

Zebrowski, J. A. [Chicago U., KICP; Chicago U., As↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

In-Situ Magnetic Field Reconstruction in the MAGIS-100 Experiment

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.; Ferm↗

DOE OSTI2025

CMB low multipole alignments across WMAP and Planck data releases

ABSTRACT The first observations of the cosmic microwave background (CMB) from NASA's Wilkinson Microwave Anisotropy Probe (WMAP) led to finding ‘alignment’ anomalies not expected from fluctuations in the isotropic cosmological model. We study the data of all 8 full-sky public releases since then to test for anomalous alignments and shapes of the first 60 multipoles, i.e. over the range $2\le l \le 61$. We use rotationally invariant and covariant statistics to test isotropy of all subsequent WMAP data releases, along with those from the ESA’s Planck mission. Anomalous alignments among the multipoles $l=1, 2, 3$ are very consistent and robust. More alignments are detected, some of them new, while significance is diluted by the large range of the search. Power entropy, a measure of the randomness of the multipoles, is consistently anomalous at about $2\sigma$ level or better across all data releases. It appears that the CMB is not as random as the cosmological principle predicts on large angular scales.

Patel, Sanjeet Kumar↗

DOE OSTI2026

pop-cosmos : redshifts and physical properties of KiDS-1000 galaxies

ABSTRACT Principled Bayesian inference of galaxy properties has not previously been performed for wide-area weak-lensing surveys with millions of sources. We address this gap by applying the pop-cosmos generative model to perform spectral energy distribution (SED) fitting for 4 million KiDS (Kilo-Degree Survey)-1000 galaxies. Calibrated on deep COSMOS2020 photometric data, pop-cosmos specifies a physically motivated prior over the galaxy population up to $z \simeq 6$ in stellar population synthesis (SPS) parameter space. Using the Speculator SPS emulator with GPU (graphics processing unit)-accelerated Markov Chain Monte Carlo sampling, we perform full posterior inference at 8.2 GPU seconds per galaxy, obtaining joint constraints on galaxy redshifts and physical properties. We validate photometric redshifts against $\sim \!185\,\!000$ KiDS galaxies cross-matched to Dark Energy Spectroscopic Instrument Data Release 1 spectroscopic samples, achieving low bias ($2\times 10^{-3}$), scatter ($\sigma _{\mathrm{MAD}}=0.03$), and outlier fraction (3.2 per cent) for the Bright Galaxy Survey, with comparable performance (bias $3\times 10^{-2}$, $\sigma _{\mathrm{MAD}}=0.05$, 1.0 per cent outliers) for luminous red galaxies (LRGs). Within the LRG sample, we identify massive, dusty, star-forming contaminants at $z \simeq 0.4$ satisfying standard colour selections for quenched populations. We infer trends in stellar mass, star formation, metallicity, and dust across five tomographic redshift bins consistent with established scaling relations. Using specific star formation rate constraints, we identify $\sim$7 per cent of KiDS-1000 galaxies as quenched, versus 37 per cent implied by conservative colour cuts. This enables the construction of weak-lensing samples defined by physical properties while mitigating intrinsic alignment systematics and preserving statistical power. Our analysis validates pop-cosmos out of sample, establishing it as a scalable approach for galaxy evolution and cosmological analyses with photometric surveys.

Halder, Anik [Institute of Astronomy and Kavli Ins↗

DOE OSTI2025

Search for dark matter from the center of the Earth with 10 years of IceCube data

The nature of dark matter remains unresolved in fundamental physics. Weakly Interacting Massive Particles (WIMPs), which could explain the nature of dark matter, can be captured by celestial bodies like the Sun or Earth, leading to enhanced self-annihilation into Standard Model particles including neutrinos detectable by neutrino telescopes such as the IceCube Neutrino Observatory. This article presents a search for muon neutrinos from the center of the Earth performed with 10 years of IceCube data using a track-like event selection. We considered a number of WIMP annihilation channels (χχ→τ+τ-$$\chi \chi \rightarrow \tau ^+\tau ^-$$/W+W-$$W^+W^-$$/bb¯$$b\bar{b}$$) and masses ranging from 10 GeV to 10 TeV. No significant excess over background due to a dark matter signal was found while the most significant result corresponds to the annihilation channel χχ→bb¯$$\chi \chi \rightarrow b\bar{b}$$ for the mass mχ=250$$m_{\chi }=250$$ GeV with a post-trial significance of 1.06σ$$1.06\sigma $$. Our results are competitive with previous such searches and direct detection experiments. Our upper limits on the spin-independent WIMP scattering are world-leading among neutrino telescopes for WIMP masses mχ>100$$m_{\chi }>100$$ GeV.

Abbasi, R↗

DOE OSTI2025

Leveraging intermediate resonances to probe CP violation at colliders

We explore the phenomenological impact of interference in tree-level contributions to three-body final states in $2\rightarrow 3$ scattering processes. This work introduces a novel search strategy leveraging asymmetries to enable sensitivity to CP-violating effects in less well-explored regions of phase space. Analytically, we demonstrate the effectiveness of this observable in probing interference between Standard Model charged-current decays and effective left-handed vector interactions, illustrated in a toy model featuring a scalar leptoquark, $S_1 \sim (3, 1, -\,1/3)$. Numerically, we apply this framework to studying the process $pp\rightarrow b \tau \nu $; unlike traditional high-$p_T$ searches or “bump hunts”, this approach utilizes an intermediate energy regime – where new physics is neither light enough to be produced on shell or heavy enough to justify an effective field theory treatment. A proof-of-principle analysis at parton level demonstrates a percent-level asymmetry, with sensitivity also to BSM weak-CP phase. While the specific phase sensitivity is diminished at particle level due to showering and detector effects, a machine learning classifier can recover sensitively to the presence of SM-BSM interference, significantly outperforming standard analysis methods. Notably discrimination between BSM signal and SM background could be achieved at the 2$\sigma $ level for the current LHC dataset and 8$\sigma $ at the High-Luminosity LHC. Moreover, this asymmetry observable as defined can also be more broadly applied to other searches for CP-violation in $2\rightarrow 3$ processes in present and future collider environments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗