Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Systematic Error”

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 919 records · Page 51

Correlations of π N partial waves for multireaction analyses

In the search for missing baryonic resonances, many analyses include data from a variety of pion- and photon-induced reactions. For elastic πN scattering, however, usually the partial waves of the SAID (Scattering Analysis Interactive Database) or other groups are fitted, instead of data. We provide the partial-wave covariance matrices needed to perform correlated χ 2 fits, in which the obtained χ 2 equals the actual χ 2 up to nonlinear and normalization corrections. For any analysis relying on partial waves extracted from elastic pion scattering, this is a prerequisite to assess the significance of resonance signals and to assign any uncertainty on results. Lastly, the influence of systematic errors is also considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Precision mass measurements of magnesium isotopes and implications for the validity of the isobaric mass multiplet equation

If the mass excess of neutron-deficient nuclei and their neutron-rich mirror partners are both known, it can be shown that deviations of the isobaric mass multiplet equation (IMME) in the form of a cubic term can be probed. Such a cubic term was probed by using the atomic mass of neutron-rich magnesium isotopes measured using the TITAN Penning trap and the recently measured proton-separation energies of 29 Cl and 30 Ar. The atomic mass of 27 Mg was found to be within 1.6σ of the value stated in the Atomic Mass Evaluation. The atomic masses of 28,29 Mg were measured to be both within 1σ, while being 7 and 33 times more precise, respectively. Using the 29 Mg mass excess and previous measurements of 29 Cl, we uncovered a cubic coefficient of d = 28(7)keV, which is the largest known cubic coefficient of the IMME. This departure, however, could also be caused by experimental data with unknown systematic errors. Hence there is a need to confirm the mass excess of 28 S and the one-neutron separation energy of 29 Cl, which have both come from a single measurement. Lastly, our results were compared with ab initio calculations from the valence-space in-medium similarity renormalization group, resulting in a good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

α -decay energies of superheavy nuclei: Systematic trends

Background: New superheavy nuclei are often identified through their characteristic α -decay energies, which requires accurate calculations of Q α values. While many Q α predictions are available, little is known about their uncertainties, and this makes it difficult to carry out extrapolations to as-yet-unknown systems. Purpose: This work aims to analyze several models, compare their predictions to available experimental data, and study their performance for the unobserved α-decay chains of 296 120 and 298 120, which are of current experimental interest. Our quantified results will also serve as a benchmark for future, more sophisticated statistical studies. Methods: We use nuclear superfluid density functional theory (DFT) with several Skyrme energy density functionals (EDFs). To estimate systematic model uncertainties, we employ uniform model averaging. Results: We investigated the Q α values for even-even nuclei from Fm to Z = 120 . For well-deformed nuclei between Fm and Ds, we find excellent consistency between different model predictions, and a good agreement with experimental results. For transitional nuclei beyond Ds, intermodel differences grow, resulting in an appreciable systematic error. In particular, our models underestimate Q α for the heaviest nucleus 294 Og. Conclusions The robustness of DFT predictions for well-deformed superheavy nuclei supports the idea of using experimental Q α values, together with theoretical predictions, as reasonable ( Z, A ) indicators. Unfortunately, this identification method is not projected to work well in the region of deformed-to-spherical shape transition as one approaches N = 184 . The use of Q α values in the identification of new superheavy nuclei will benefit greatly from progress in developing both new spectroscopic-quality EDFs and more sophisticated statistical techniques of uncertainty quantification.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Prompt fission product yields in the U 238 ( n , f ) reaction

Background: Significant yield discrepancies (500-600%) were reported recently between experi- mental results and predictions (from the GEF model) and evaluations (from the JEFF-3.1.1 and ENDF/B-VII.1 libraries) for Mo and Sn fission-fragment yields in fast-neutron induced reactions on 238 U using γ-γ-γ coincidence spectroscopy. The model/evaluations also predict Mo and Sn fragments that are on average ~1 to 2 neutrons richer than the experimental results. Purpose: γ-γ-γ coincidence spectroscopy favors detection of higher-multiplicity γ-ray cascades. An alternative approach is determining the fragment yields using single-γ-ray spectroscopy, as it was attempted here for selected cases where it was feasible. Advantages/drawbacks in both approaches need to be understood and potential systematic errors in the experimental results should be addressed using theoretical models. Methods: Fast neutrons from the LANSCE WNR facility were used to induce fission on 238 U. The emitted γ rays were measured with the GEANIE spectrometer. Results: The yield of selected even-even fission fragments was determined. The selection was based on the ability to reliably determine excitation functions for the detected γ rays. Conclusions: Our single-γ-ray results provide better agreement between experiment and predictions/evaluations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Machine-learning prediction for quasiparton distribution function matrix elements

There have been rapid developments in the direct calculation in lattice QCD (LQCD) of the Bjorken-x dependence of hadron structure through large-momentum effective theory (LaMET). LaMET overcomes the previous limitation of LQCD to moments (that is, integrals over Bjorken x) of hadron structure, allowing LQCD to directly provide the kinematic regions where the experimental values are least known. LaMET requires large-momentum hadron states to minimize its systematics and allow us to reach small- x reliably. This means that very fine lattice spacing to minimize lattice artifacts at order $(P_za)^n$ will become crucial for next-generation LaMET-like structure calculations. Furthermore, such calculations require operators with long Wilson-link displacements, especially in finer lattice units, increasing the communication costs relative to that of the propagator inversion. In this work, we explore whether machine-learning algorithms can make predictions of correlators to reduce the computational cost of these LQCD calculations. We consider two algorithms, gradient-boosting decision tree and linear models, applied to LaMET data, the matrix elements needed to determine the kaon and $η_s$ unpolarized parton distribution functions (PDFs), meson distribution amplitude (DA), and the nucleon gluon PDF. We find that both algorithms can reliably predict the target observables with different prediction accuracy and systematic errors. The predictions from smaller displacement z to larger ones work better than those for momentum p due to the higher correlation among the data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dark Energy Survey Year 1 Results: Cosmological constraints from cluster abundances and weak lensing

In this work we perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =\sigma_8(\Omega_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $\Omega_{\rm m}=0.179^{+0.031}_{-0.038}$, with $\sigma_8-\Omega_{\rm m}$ posteriors in $2.4\sigma$ tension with the DES Y1 3x2pt results, and in $5.6\sigma$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($\lambda \ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Nucleon mass with highly improved staggered quarks

We present the first computation in a program of lattice-QCD baryon physics using staggered fermions for sea and valence quarks. For this initial study, we present a calculation of the nucleon mass, obtaining 964 ±16 MeV with all sources of statistical and systematic errors controlled and accounted for. This result is the most precise determination to date of the nucleon mass from first principles. We use the highly improved staggered quark action, which is computationally efficient. Three gluon ensembles are employed, which have approximate lattice spacings 𝑎 ≈ 0.09, 0.12, and 0.15 fm, each with equal-mass 𝑢/𝑑, 𝑠, and 𝑐 quarks in the sea. Further, all ensembles have the light valence and sea 𝑢/𝑑 quarks tuned to reproduce the physical pion mass, avoiding complications from chiral extrapolations. Our work opens a new avenue for precise calculations of baryon properties, which are both feasible and relevant to experiments in particle and nuclear physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Bayesian model averaging for analysis of lattice field theory results

Statistical modeling is a key component in the extraction of physical results from lattice field theory calculations. Although the general models used are often strongly motivated by physics, many model variations can frequently be considered for the same lattice data. Model averaging, which amounts to a probability-weighted average over all model variations, can incorporate systematic errors associated with model choice without being overly conservative. We discuss the framework of model averaging from the perspective of Bayesian statistics, and give useful formulae and approximations for the particular case of least-squares fitting, commonly used in modeling lattice results. In addition, we frame the common problem of data subset selection (e.g. choice of minimum and maximum time separation for fitting a two-point correlation function) as a model selection problem and study model averaging as a straightforward alternative to manual selection of fit ranges. Numerical examples involving both mock and real lattice data are given.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Lattice determination of 𝐼=0 and 2 𝜋⁢𝜋 scattering phase shifts with a physical pion mass

Phase shifts for 𝑠-wave 𝜋⁢𝜋 scattering in both the 𝐼 =0 and 𝐼 =2 channels are determined from a lattice QCD calculation performed on 741 gauge configurations obeying G-parity boundary conditions with a physical pion mass and lattice size of 32 3 ×64. These results support our recent 2021 study of direct 𝐶⁢𝑃 violation in 𝐾 → 𝜋⁢𝜋 decay, improving our earlier 2015 calculation. The phase shifts are determined for both stationary and moving 𝜋⁢𝜋 systems, at three (𝐼 = 0) and four (𝐼 = 2) different total momenta. We implement several 𝜋⁢𝜋 interpolating operators including a scalar bilinear “𝜎” operator and paired single-pion bilinear operators with the constituent pions carrying various relative momenta. Several techniques, including correlated fitting and a bootstrap determination of p-values have been used to refine the results and a comparison with the generalized eigenvalue problem method is given. A detailed systematic error analysis is performed which allows phase shift results to be presented at a fixed energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comprehensive symmetric-hybrid ring design for a proton EDM experiment at below 10 –29 e · cm

A concise demonstrative summary of the Symmetric Hybrid ring design for the storage ring proton electric dipole moment experiment is presented. In this paper, critical issues such as lattice design, background electrical fields, geometrical phase, general relativity, spin coherence time and polarimeter systematics are presented. Overall, we find that with the currently proposed design iteration, systematic error sources are reduced by orders of magnitude and that the ring alignment requirements are within the currently available technology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

BOSS DR12 full-shape cosmology: Λ CDM constraints from the large-scale galaxy power spectrum and bispectrum monopole

Here, we present a full Λ CDM analysis of the BOSS DR12 dataset, including information from the power spectrum multipoles, the real-space power spectrum, the reconstructed power spectrum and the bispectrum monopole. This is the first analysis to feature a complete treatment of the galaxy bispectrum, including a consistent theoretical model and without large-scale cuts. Unlike previous works, the statistics are measured using window-free estimators: this greatly reduces computational costs by removing the need to window-convolve the theory model. Our pipeline is tested using a suite of high-resolution mocks and shown to be robust and precise, with systematic errors far below the statistical thresholds. Inclusion of the bispectrum yields consistent parameter constraints and shrinks the σ 8 posterior by 13% to reach < 5 % precision; less conservative analysis choices would reduce the error bars further. Our constraints are broadly consistent with Planck: in particular, we find H 0 = 69.6$^{+1.1}_{–1.3}$ km s –1 Mpc –1 , σ 8 = 0.692$^{+0.035}_{–0.041}$ and n s = 0.870$^{+0.067}_{–0.064}$, including a BBN prior on the baryon density. When n s is set by Planck, we find H 0 = 68.31$^{+0.83}_{–0.86}$ km s –1 Mpc –1 and σ 8 = 0.722$^{+0.032}_{–0.036}$. Our S 8 posterior, 0.751 ± 0.039, is consistent with weak lensing studies, but lower than Planck. Constraints on the higher-order bias parameters are significantly strengthened from the inclusion of the bispectrum, and we find no evidence for deviation from the dark matter halo bias relations. These results represent the most complete full-shape analysis of BOSS DR12 to-date, and the corresponding spectra will enable a variety of beyond-Λ CDM analyses, probing phenomena such as the neutrino mass and primordial non-Gaussianity.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of e + e - → η Φ via initial state radiation at Belle

Using 980 fb -1 of data collected on and around the (n = 1, 2, 3, 4, 5) resonances with the Belle detector at the KEKB collider, we measure the cross section of e + ⁢e - → η⁢Φ from threshold to 3.95 GeV via initial state radiation. There are clear Φ⁡(1680) and J/ψ signals but no significant Φ⁡(2170) signal in the ηΦ final state. The branching fraction $\mathscr{B}$⁡[J/ψ →η⁢Φ] is measured to be (7.2 ±0.8 ±0.5) ×10 -4 . The resonant parameters of Φ⁡(1680) are determined to be m Φ(1680) = (1696±8±10) MeV/c 2 (statistical and systematic errors, respectively), Γ Φ⁡(1680) =(175±13±16) MeV and, depending on the possible presence of predominantly constructive or destructive interference between Φ⁡(1680) and continuum production, Γ$^{e^+e^-}_{Φ⁡(1680)}$·$\mathscr{B}$⁡[Φ⁡(1680)→ηΦ] and $\mathscr{B}$⁡[Φ(1680)→η⁢Φ] are determined to be (75 ±10 ±11) eV and (25±12±2)% or (207±16±20) eV and (23 ±10 ±2)%, respectively. The upper limit for Γ$^{e^+e^-}_{Φ⁡(2170)}$·$\mathscr{B}$⁡[Φ⁡(2170) → ηΦ] is determined to be either 0.17 or 18.6 eV at the 90% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dark Energy Survey Year 3 results: Constraints on extensions to Λ CDM with weak lensing and galaxy clustering

We constrain six possible extensions to the Λ cold dark matter (CDM) model using measurements from the Dark Energy Survey’s first three years of observations, alone and in combination with external cosmological probes. The DES data are the two-point correlation functions of weak gravitational lensing, galaxy clustering, and their cross-correlation. We use simulated data vectors and blind analyses of real data to validate the robustness of our results to astrophysical and modeling systematic errors. In many cases, constraining power is limited by the absence of theoretical predictions beyond the linear regime that are reliable at our required precision. The Λ CDM extensions are dark energy with a time-dependent equation of state, nonzero spatial curvature, additional relativistic degrees of freedom, sterile neutrinos with eV-scale mass, modifications of gravitational physics, and a binned σ 8 ( z ) model which serves as a phenomenological probe of structure growth. For the time-varying dark energy equation of state evaluated at the pivot redshift we find ( w p , w a ) = ( - 0.9 9 - 0.17 + 0.28 , - 0.9 ± 1.2 ) at 68% confidence with z p = 0.24 from the DES measurements alone, and ( w p , w a ) = ( - 1.0 3 - 0.03 + 0.04 , - 0. 4 - 0.3 + 0.4 ) with z p = 0.21 for the combination of all data considered. Curvature constraints of Ω k = 0.0009 ± 0.0017 and effective relativistic species N eff = 3.1 0 - 0.16 + 0.15 are dominated by external data, though adding DES information to external low-redshift probes tightens the Ω k constraints that can be made without cosmic microwave background observables by 20%. For massive sterile neutrinos, DES combined with external data improves the upper bound on the mass m eff by a factor of 3 compared to previous analyses, giving 95% limits of ( Δ N eff , m eff ) ≤ ( 0.28 , 0.20 eV ) when using priors matching a comparable Planck analysis. For modified gravity, we constrain changes to the lensing and Poisson equations controlled by functions Σ ( k , z ) = Σ 0 Ω Λ ( z ) / Ω Λ , 0 and μ ( k , z ) = μ 0 Ω Λ ( z ) / Ω Λ , 0 , respectively, to Σ 0 = 0. 6 - 0.5 + 0.4 from DES alone and ( Σ 0 , μ 0 ) = ( 0.04 ± 0.05 , 0.0 8 - 0.19 + 0.21 ) for the combination of all data, both at 68% confidence. Overall, we find no significant evidence for physics beyond Λ CDM .

79 ASTRONOMY AND ASTROPHYSICS↗

Extraction of the neutron F 2 structure function from inclusive proton and deuteron deep-inelastic scattering data

The available world deep-inelastic scattering (DIS) data on proton and deuteron structure functions F 2 p , F 2 d , and their ratios are leveraged to extract the free neutron F 2 n structure function, the F 2 n / F 2 p ratio, and associated uncertainties using the latest nuclear effect calculations in the deuteron. Special attention is devoted to the normalization of the proton and deuteron experimental datasets and to the treatment of correlated systematic errors, as well as the quantification of procedural and theoretical uncertainties. The extracted F 2 n dataset is utilized to evaluate the Q 2 dependence of the Gottfried sum rule and the nonsinglet F 2 p − F 2 n moments. To facilitate replication of our study, as well as for general applications, we provide a comprehensive DIS database including all recent Jefferson Lab 6 GeV measurements, the extracted F n 2 , a modified CTEQ-JLab global parton distribution function fit named CJ15nlo_mod, and grids with calculated proton, neutron, and deuteron DIS structure functions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Constraining gravity with a new precision 𝐸 𝐺 estimator using Planck + SDSS BOSS data

The 𝐸 𝐺 statistic is a discriminating probe of gravity developed to test the prediction of general relativity (GR) for the relation between gravitational potential and clustering on the largest scales in the observable Universe. We present a novel high-precision estimator for the 𝐸 𝐺 statistic using CMB lensing and galaxy clustering correlations that carefully matches the effective redshifts across the different measurement components to minimize corrections. A suite of detailed tests is performed to characterize the estimator’s accuracy, its sensitivity to assumptions and analysis choices, and the non-Gaussianity of the estimator’s uncertainty is characterized. After finalization of the estimator, it is applied to Planck CMB lensing and SDSS CMASS and LOWZ galaxy data. We report the first harmonic space measurement of 𝐸 𝐺 using the LOWZ sample and CMB lensing and also updated constraints using the final CMASS sample and the latest Planck CMB lensing map. We find $\hat{𝐸}$$^{Planck+CMASS}_{𝐺}$ = 0.3⁢6$^{+0.06}_{−0.05}$⁢(68.27%) and $\hat{𝐸}$$^{Planck+LOWZ}_{𝐺}$ = 0.4⁢0$^{+0.11}_{−0.09}$⁢(68.27%), with additional subdominant systematic error budget estimates of 2% and 3%, respectively. Using Ω m,0 constraints from Planck and SDSS BAO observations, Λ⁢CDM-GR predicts 𝐸$^{GR}_ {𝐺}$⁡(𝑧 =0.555) = 0.401 ± 0.005 and 𝐸$^{GR}_{𝐺}$⁡(𝑧 =0.316) = 0.452 ± 0.005 at the effective redshifts of the CMASS and LOWZ based measurements. We report the measurement to be in good statistical agreement with the Λ⁢CDM-GR prediction and report that the measurement is also consistent with the more general GR prediction of scale independence for 𝐸 𝐺 . Furthermore, this work provides a carefully constructed and calibrated statistic with which 𝐸 𝐺 measurements can be confidently and accurately obtained with upcoming survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

Trace anomaly form factors from lattice QCD

The hadron mass can be obtained through the calculation of the trace of the energy-momentum tensor in the hadron which includes the trace anomaly and sigma terms. The anomaly due to conformal symmetry breaking is believed to be an important ingredient for hadron mass generation and confinement. In this work, we will present the calculation of the glue part of the trace anomaly form factors of the pion up to Q 2 ∼ 4.3 GeV 2 and the nucleon up to Q 2 ∼ 1 GeV 2 . The calculations are performed on a domain wall fermion ensemble with overlap valence quarks at seven valence pion masses varying from ∼ 250 to ∼ 540 MeV , including the unitary point ∼ 340 MeV . We calculate the radius of the glue trace anomaly for the pion and the nucleon from the z expansion. By performing a two-dimensional Fourier transform on the glue trace anomaly form factors in the infinite momentum frame with no energy transfer, we also obtain their spatial distributions for several valence quark masses. The results are qualitatively extrapolated to the physical valence pion mass with systematic errors from the unphysical sea quark mass, discretization effects in the renormalization sum rule, and finite-volume effects to be addressed in the future. We find the pion’s form factor changes sign, as does its spatial distribution, for light quark masses. This explains how the trace anomaly contribution to the pion mass approaches zero toward the chiral limit. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Kaon mixing beyond the standard model with physical masses

We present nonperturbative results for beyond the standard model kaon mixing matrix elements in the isospin symmetric limit ( m u = m d ) of QCD, including a complete estimate of all dominant sources of systematic error. Our results are obtained from numerical simulations of lattice QCD with N f = 2 + 1 flavors of dynamical domain wall fermions. For the first time, these quantities are simulated directly at the physical pion mass m π ∼ 139 MeV for two different lattice spacings. We include data at three lattice spacings in the range a = 0.11 – 0.07 fm and with pion masses ranging from the physical value up to 450 MeV. Compared to our earlier work, we have added both direct calculations at physical quark masses and a third lattice spacing making the removal of discretization effects significantly more precise and eliminating the need for any significant mass extrapolation beyond the range of simulated data. We renormalize the lattice operators nonperturbatively using RI-SMOM off-shell schemes. These schemes eliminate the need to model and subtract nonperturbative pion poles that arises in the RI-MOM scheme and, since the calculations are performed with domain wall fermions, the unphysical mixing between chirality sectors is suppressed. Our results for the bag parameters in the MS ¯ scheme at 3 GeV are B K ≡ B 1 = 0.5240 ( 17 ) ( 54 ) , B 2 = 0.4794 ( 25 ) ( 35 ) , B 3 = 0.746 ( 13 ) ( 17 ) , B 4 = 0.897 ( 02 ) ( 10 ) and B 5 = 0.6882 ( 78 ) ( 94 ) , where the first error is from lattice uncertainties and the second is the uncertainty due to the perturbative matching to MS ¯ . Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Exploration of the polarization angle variability of the Crab Nebula with POLARBEAR and its application to the search for axionlike particles

The Crab Nebula, also known as Tau A, is a polarized astronomical source at millimeter wavelengths. It has been used as a stable light source for polarization angle calibration in millimeter-wave astronomy. However, it is known that its intensity and polarization vary as a function of time at a variety of wavelengths. Thus, it is of interest to verify the stability of the millimeter-wave polarization. If detected, polarization variability may be used to better understand the dynamics of Tau A, and for understanding the validity of Tau A as a calibrator. One intriguing application of such observation is to use it for the search of axionlike particles (ALPs). Ultralight ALPs couple to photons through a Chern-Simons term, and induce a temporal oscillation in the polarization angle of linearly polarized sources. After assessing a number of systematic errors and testing for internal consistency, we evaluate the variability of the polarization angle of the Crab Nebula using 2015 and 2016 observations with the 150 GHz P instrument. We place a median 95% upper bound of polarization oscillation amplitude A < 0.06 5 ° over the oscillation frequencies from 0.75 year − 1 to 0.66 hour − 1 . Assuming that no sources other than ALP are causing Tau A’s polarization angle variation, that the ALP constitutes all the dark matter, and that the ALP field is a stochastic Gaussian field, this bound translates into a median 95% upper bound of ALP-photon coupling g a γ γ < 2.16 × 10 − 12 GeV − 1 × ( m a / 10 − 21 eV ) in the mass range from 9.9 × 10 − 23 eV to 7.7 × 10 − 19 eV . This demonstrates that this type of analysis using bright polarized sources is as competitive as those using the polarization of cosmic microwave background in constraining ALPs. Published by the American Physical Society 2024

Adachi, Shunsuke (ORCID:0000000204007555)↗