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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Search for production of four top quarks in final states with same-sign or multiple leptons in proton–proton collisions at $\sqrt{s}=13$ ${TeV}$

The standard model (SM) production of four top quarks ($\text {t} {}{\overline{\text {t}}} \text {t} {}{\overline{\text {t}}} $) in proton–proton collisions is studied by the CMS Collaboration. The data sample, collected during the 2016–2018 data taking of the LHC, corresponds to an integrated luminosity of 137$\,\text {fb}^{-1}$ at a center-of-mass energy of 13$\,\text {TeV}$. The events are required to contain two same-sign charged leptons (electrons or muons) or at least three leptons, and jets. The observed and expected significances for the $\text {t} {}{\overline{\text {t}}} \text {t} {}{\overline{\text {t}}} $ signal are respectively 2.6 and 2.7 standard deviations, and the $\text {t} {}{\overline{\text {t}}} \text {t} {}{\overline{\text {t}}} $ cross section is measured to be $12.6^{+5.8}_{-5.2}\,\text {fb} $. The results are used to constrain the Yukawa coupling of the top quark to the Higgs boson, $y_{\text {t}}$, yielding a limit of $|y_{\text {t}}/y_{\text {t}}^{\mathrm {SM}} | < 1.7$ at $95\%$ confidence level, where $y_{\text {t}}^{\mathrm {SM}}$ is the SM value of $y_{\text {t}}$. They are also used to constrain the oblique parameter of the Higgs boson in an effective field theory framework, $\hat{H}<0.12$. Limits are set on the production of a heavy scalar or pseudoscalar boson in Type-II two-Higgs-doublet and simplified dark matter models, with exclusion limits reaching 350–470$\,\text {GeV}$ and 350–550$\,\text {GeV}$ for scalar and pseudoscalar bosons, respectively. Upper bounds are also set on couplings of the top quark to new light particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Nucleon isovector axial form factors

We present results for the isovector axial vector form factors obtained using thirteen 2 + 1 + 1 -flavor highly improved staggered quark (HISQ) ensembles generated by the MILC collaboration. The calculation of nucleon two- and three-point correlation functions has been done using Wilson-clover fermions. In the analysis of these data, we quantify the sensitivity of the results to strategies used for removing excited state contamination and invoke the partially conserved axial current relation between the form factors to choose between them. Our data driven analysis includes removing contributions from multihadron N π states that make significant contributions. Our final results are g A = 1.292 ( 53 ) stat ( 24 ) sys for the axial charge; g S = 1.085 ( 50 ) stat ( 103 ) sys and g T = 0.991 ( 21 ) stat ( 10 ) sys for the scalar and tensor charges; ⟨ r A 2 ⟩ = 0.439 ( 56 ) stat ( 34 ) sys fm 2 for the mean squared axial charge radius, g P * = 9.03 ( 47 ) stat ( 42 ) sys for the induced pseudoscalar charge; and g π N N = 14.14 ( 81 ) stat ( 85 ) sys for the pion-nucleon coupling. We also provide a parametrization of the axial form factor G A ( Q 2 ) over the range 0 ≤ Q 2 ≤ 1 GeV 2 for use in phenomenology and a comparison with other lattice determinations. We find that the various lattice data agree within 10% but are significantly different from the extraction of G A ( Q 2 ) from the ν -deuterium scattering data. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Probing New Hadronic Forces with Heavy Exotic Atoms

We explore the potential of precision spectroscopy of heavy exotic atoms where electrons are substituted by negative hadrons to detect new force carriers with hadronic couplings. The selected transitions are unaffected by nuclear contact terms, thus enabling highly accurate calculations using bound-state QED, provided that the nuclear polarization is under control. Alternatively, we demonstrate that the dipole polarizability, a fundamental property of nuclei, can be extracted from the spectroscopy of exotic atoms in a novel way by combining two transitions while maintaining high sensitivity to new physics. Based on existing data, we extracted world-leading bounds on mediator masses ranging from 0.1 to 10 MeV for two benchmark models and show that forthcoming experiments could enhance the sensitivity to new physics by 2 orders of magnitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Differential equations for cosmological correlators

Cosmological fluctuations retain a memory of the physics that generated them in their spatial correlations. The strength of correlations varies smoothly as a function of external kinematics, which is encoded in differential equations satisfied by cosmological correlation functions. In this work, we provide a broader perspective on the origin and structure of these differential equations. As a concrete example, we study conformally coupled scalar fields in a power-law cosmology. The wavefunction coefficients in this model have integral representations, with the integrands being the product of the corresponding flat-space results and “twist factors” that depend on the cosmological evolution. Similar twisted integrals arise for loop amplitudes in dimensional regularization, and their recent study has led to the discovery of rich mathematical structures and powerful new tools for computing multi-loop Feynman integrals in quantum field theory. The integrals of interest in cosmology are also part of a finite-dimensional basis of master integrals, which satisfy a system of first-order differential equations. We develop a formalism to derive these differential equations for arbitrary tree graphs. The results can be represented in graphical form by associating the singularities of the differential equations with a set of graph tubings. Upon differentiation, these tubings grow in a local and predictive fashion. In fact, a few remarkably simple rules allow us to predict — by hand — the equations for all tree graphs. While the rules of this “kinematic flow” are defined purely in terms of data on the boundary of the spacetime, they reflect the physics of bulk time evolution. We also study the analogous structures in tr ϕ 3 theory, and see some glimpses of hidden structure in the sum over planar graphs. This suggests that there is an autonomous combinatorial or geometric construction from which cosmological correlations, and the associated spacetime, emerge.

Cosmological models↗
DOE OSTI2026

Mesogenesis through the ephemeral dark decay of Beauty

Mesogenesis provides a path for generating the baryon asymmetry of the Universe, using only the CP violation furnished by the Standard Model in the decay of B mesons. While this is an intriguing possibility, it is largely constrained by the data on B meson branching fractions into baryons and missing energy carried into the dark sector. We point out that it is possible to make this branching fraction dominant only in the early Universe, through an ultralight scalar coupled to the dark sector and the Standard Model leptons. A scenario is examined where the thermal density of muons in the early Universe temporarily lowers the mass of a dark fermion, allowing for efficient B meson decays. This `dark’ decay channel is shut off later when the muon number density falls, making the scenario compatible with flavor data. Our model can be consistent with the LHC constraints on color-charged heavy bosons required to implement Mesogenesis; such states may be discovered in the future runs as their masses cannot be far above the current bounds. We also outline other possible signals, which can arise in future displaced vertex searches, long range force searches, and observations of neutron star binary mergers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Demonstration of <2 ns timing resolution for neutrino interaction in the MicroBooNE detector

The MicroBooNE detector, located in the Booster Neutrino Beamline (BNB) at Fermilab, collected data from 2015 to 2021. MicroBooNE’s Liquid Argon Time Projection Chamber is accompanied by a Photon Detection System consisting of 32 PMTs used to measure the argon scintillation light and determine the timing of the neutrino interactions. This work employs improved analysis techniques combining light signals and reconstructed tracks achieving neutrino interaction timing resolution <2 ns. The result obtained allows MicroBooNE to access the 2 ns neutrino pulse structure of the BNB for the first time. The timing resolution achieved enables significant enhancement of cosmic background rejection for all neutrino analyses. Furthermore, the ns timing resolution opens the door for searching new long-lived-particles (e.g. Heavy Neutral Leptons, Higgs Portal Scalars) in MicroBooNE, as well in future large LArTPC experiments, namely SBN and DUNE.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

tapir: A tool for topologies, amplitudes, partial fraction decomposition and input for reductions

The demand for precision predictions in the field of high energy physics has dramatically increased over recent years. Experiments conducted at the LHC, as well as precision measurements at the intensity frontier such as Belle II require equally precise theoretical predictions to make full use of the acquired data. To match the experimental precision, second-, third- and, for certain quantities, even higher-order calculations in perturbative quantum field theory are required. To facilitate such calculations, computer software automating as many steps as possible is required. Yet, each calculation poses different challenges and thus, a high level of configurability is required. In this context we present tapir: a tool for identification, manipulation and minimization of Feynman integral families. It is designed to integrate in toolchains based on the computer algebra system FORM, the use of which is common practice in the field. tapir can be used to reduce the complexity of multi-loop problems with cut-filters, topology mapping, partial fraction decomposition and alike. Program Title:tapir CPC Library link to program files:https://doi.org/10.17632/ptc9t46xyn.1 Developer's repository link:https://gitlab.com/tapir-devs/tapir Licensing provisions: GPLv3 Programming language:python 3, C++ Nature of problem: Multi-loop computations require the automatization of a large number of different tasks related to Feynman integral topologies. Among them are the identification and minimization of integral topologies, partial fraction decomposition of topologies in the case of linearly dependent propagators as well as mapping scalar products of loop momenta to scalar functions. Solution method: The minimization of topologies is performed by comparison of their respective Nickel indices [1], even further minimization utilizes Pak's algorithm [2]. To efficiently map scalar products of loop momenta to scalar functions FORM [3] code is generated. Additional comments including restrictions and unusual features: Minimization based on Pak's algorithm slows down for many lines and scales. A coarser minimization using the Nickel indices, however, is still possible. [1]B. Nickel, D. Meiron, G.A.J. Baker, Compilation of 2-pt and 4-pt graphs for continuous spin model, Report, University of Guelph, 1977.[2]A. Pak, J. Phys. Conf. Ser. 368 (2012) 012049, https://doi.org/10.1088/1742-6596/368/1/012049, arXiv:1111.0868.[3]B. Ruijl, T. Ueda, J. Vermaseren, FORM version 4.2, arXiv:1707.06453, 7 2017.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Probing parity violation with the four-point correlation function of BOSS galaxies

Parity-violating physics in the early universe can leave detectable traces in late-time observables. While vector- and tensor-type parity violation can be observed in the B -modes of the cosmic microwave background, scalar-type signatures are visible only in the four-point correlation function (4PCF) and beyond. This work presents a blind test for parity violation in the 4PCF of the baryon oscillation spectroscopic survey (BOSS) CMASS sample, considering galaxy separations in the range [20; 160]h -1 Mpc The parity-odd 4PCF contains no contributions from standard Λ CDM physics and can be efficiently measured using recently developed estimators. Data are analyzed using both a nonparametric rank test (comparing the BOSS 4PCFs to those of realistic simulations) and a compressed χ 2 analysis, with the former avoiding the assumption of a Gaussian likelihood. These find similar results, with the rank test giving a detection probability of 99.6% (2.9 σ). This provides significant evidence for parity violation, from either cosmological sources or systematics. We perform a number of systematic tests: although these do not reveal any observational artifacts, we cannot exclude the possibility that our detection is caused by the simulations not faithfully representing the statistical properties of the BOSS data. Our measurements can be used to constrain physical models of parity violation. As an example, we consider a coupling between the inflaton and a U (1) gauge field and place bounds on the latter’s energy density, which are several orders of magnitude stronger than those previously reported. Upcoming probes such as DESI and Euclid will reveal whether our detection of parity violation is due to new physics, and strengthen the bounds on a variety of models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Leptonic cascade decays of a heavy Higgs boson through vectorlike leptons at the LHC

We demonstrate the potential of fully leptonic cascade decays of a heavy neutral Higgs boson through vectorlike leptons as a simultaneous probe for extended Higgs sectors and extra matter particles at the LHC. The processes we explore are unique in that their event topologies lead to di-boson-like leptonic final states with a lepton pair which does not reconstruct the mass of a gauge boson. By recasting existing $2\ell$ + $E^{\text{miss}}_{R}$ and $3/4\ell$ searches channels using run2 data from the LHC we obtain model independent bounds on the masses of heavy scalars and vectorlike leptons and use these results to explore future prospects at the HL-LHC. Our results can be directly applied to any kind of new physics scenarios sharing the final states and the event topology. For concreteness, we apply our results to a benchmark scenario: a two Higgs doublet model type-II augmented with vectorlike leptons. Remarkably, even with current data the sensitivity of our analysis shows a reach for masses of a heavy neutral Higgs and vectorlike leptons up to 2 TeV and 1.5 TeV, respectively. Even for low tan β ≳ 1, the analysis retains sensitivity to heavy Higgs masses slightly above 1 TeV. The future sensitivities at the HL-LHC extend the reach for heavy Higgses and new leptons to 2.7 TeV and 2 TeV, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

MOOSE ProbML: Parallelizable Probabilistic Machine Learning and Uncertainty Quantification Capabilities

The Multiphysics Object Oriented Simulation Environment (MOOSE) is a widely used open- source finite element software for performing multiphysics multiscale simulations in a massively parallel fashion. Recently, the computational team at Idaho National Laboratory (INL) has implemented Probabilistic Machine Learning (ProbML) capabilities in MOOSE—in a parallelized fashion—and enable active learning with large-scale computational models for tasks such as surrogate model development, scale bridging, forward/inverse uncertainty quantification (UQ), Bayesian optimization, etc. This presentation summarizes these developments in MOOSE along with demonstrations on several real applications relevant to nuclear energy. At the fundamental level, samplers like Monte Carlo/Latin Hypercube, variance reduction, parallelized Markov Chain Monte Carlo (MCMC) support uncertainty propagation in both forward and inverse settings. These samplers can be integrated with the Gaussian processes (GP) suite in MOOSE, which offer several variants like scalar GPs, multi-output GPs, and deep GPs, to enable active learning. These GPs can be tuned using gradient-based optimization methods like Adam and its variants or gradient-free methods like the elliptical slice sampler (a variant of MCMC adept under Gaussian settings) for more complex covariance kernels or likelihoods whose gradient computations can be cumbersome. A variety of batch acquisition functions permit parallelized evaluation of the computational model and support different learning objectives with high efficiency like Bayesian inference, global surrogate development, optimization, etc. Furthermore, libtorch integration supports training, evaluation, and re-training of neural networks and other complex machine learning models in active learning settings. The impacts of these developments are shown on several real applications: (1) nuclear fuel inverse UQ and model inadequacy assessment using the Kennedy O’Hagan framework; (2) uncertainty aware surrogate modeling for additive manufacturing to predict field quantities; (3) nuclear reactor rare events analysis; and (4) complex fluid flow prediction using a global surrogate with quantified prediction uncertainty. Finally, the outlook of MOOSE ProbML is discussed for both outer-loop and inner-loop computations in the broad view to accelerate fuels and materials qualification, address gaps in knowledge and data, and assess new reactor/fuel systems.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Investigation of CMB constraints for dark matter-helium scattering

In this work, we study dark matter-helium scattering in the early Universe and its impact on constraints from cosmic microwave background (CMB) anisotropy measurements. We describe possible theoretical frameworks for dark matter-nucleon interactions via a scalar, pseudoscalar, or vector mediator; such interactions give rise to hydrogen and helium scattering, with cross sections that have a power-law dependence on relative velocity. Within these frameworks, we consider three scenarios: dark matter coupling to only neutrons, to only protons, and to neutrons and protons with equal strength. For these various cases, we use Planck 2018 temperature, polarization, and lensing anisotropy data to place constraints on dark matter scattering with hydrogen and/or helium for dark matter masses between 10 keV and 1 TeV. For any model that permits both helium and hydrogen scattering with a non-negative power-law velocity dependence, we find that helium scattering dominates the constraint for dark matter masses well above the proton mass. Furthermore, we place the first CMB constraints on dark matter that scatters dominantly/exclusively with helium in the early Universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Bringing the Peccei-Quinn mechanism down to Earth

It is conventionally assumed that the physics underlying the Peccei-Quinn (PQ) mechanism for addressing the strong C P problem is at very high energies, orders of magnitude above the weak scale. However, this may not be the case in general and the associated PQ boson ϕ , besides the signature state, i.e., the ultralight axion a , may emerge well below the weak scale. We consider this possibility and examine some of the conditions for its viability. The example model proposed here may also provide the requisite Standard Model Higgs mass parameter, without invoking new scalars above the giga-electron-volt (GeV) scale. The corresponding parameter space can maintain against quantum corrections. This scenario, depending on the choice of parameters, can potentially be constrained by flavor data. We point out that the current mild excess in B + → K + ν ν ¯ , reported by the Belle II experiment, could in principle be explained in this setup as B + → K + ϕ and B + → K + a , with both ϕ and a escaping the detector as missing energy. For a sufficiently heavy PQ boson, in the GeV regime (as the preferred explanation of the Belle II excess), one can separate these two contributions, due to the difference in K + momenta. In this case, the axion may also affect lighter meson, e.g., kaon, decays while ϕ would not be a kinematically allowed final state. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Local bilinear computation of Jacobi sets

Abstract We propose a novel method for the computation of Jacobi sets in 2D domains. The Jacobi set is a topological descriptor based on Morse theory that captures gradient alignments among multiple scalar fields, which is useful for multi-field visualization. Previous Jacobi set computations use piecewise linear approximations on triangulations that result in discretization artifacts like zig-zag patterns. In this paper, we utilize a local bilinear method to obtain a more precise approximation of Jacobi sets by preserving the topology and improving the geometry. Consequently, zig-zag patterns on edges are avoided, resulting in a smoother Jacobi set representation. Our experiments show a better convergence with increasing resolution compared to the piecewise linear method. We utilize this advantage with an efficient local subdivision scheme. Finally, our approach is evaluated qualitatively and quantitatively in comparison with previous methods for different mesh resolutions and across a number of synthetic and real-world examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Kinematic Flow and the Emergence of Time

Perhaps the most basic question we can ask about cosmological correlations is how their strength changes as we smoothly vary kinematic parameters. The answer is encoded in differential equations that govern this evolution in kinematic space. In this Letter, we introduce a new perspective on these differential equations. We show that, in the simplified setting of conformally coupled scalars in power-law Friedmann-Robertson-Walker spacetimes, the equations for arbitrary tree-level processes can be obtained from a small number of simple combinatorial rules. While this “kinematic flow” is defined purely in terms of boundary data, it reflects the physics of bulk time evolution. The unexpected regularity of the equations suggests the existence of an autonomously defined mathematical structure from which cosmological correlations and the time evolution of the associated spacetime emerge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dark sink enhances the direct detection of freeze-in dark matter

We describe a simple dark sector structure which, if present, has implications for the direct detection of dark matter (DM); the dark sink. A dark sink transports energy density from the DM into light dark-sector states that do not appreciably contribute to the DM density. As an example, we consider a light, neutral fermion ψ which interacts solely with DM Χ via the exchange of a heavy scalar Φ. We illustrate the impact of a dark sink by adding one to a DM freeze-in model in which Χ couples to a light dark photon γ' which kinetically mixes with the Standard Model (SM) photon. This freeze-in model (absent the sink) is itself a benchmark for ongoing experiments. In some cases, the literature for this benchmark has contained errors; we correct the predictions and provide them as a public code. We then analyze how the dark sink modifies this benchmark, solving coupled Boltzmann equations for the dark-sector energy density and DM yield. We check the contribution of the dark sink ψ’s to dark radiation; consistency with existing data limits the maximum attainable cross section. For DM with a mass between MeV –Ο⁡(10 GeV), adding the dark sink can increase predictions for the direct detection cross section all the way up to the current limits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Dark Matter Reconstruction in LBAI Experiments with Imperfect Data

Long-baseline atom interferometer (LBAI) experiments offer unprecedented sensitivity to ultralight scalar dark matter (DM) [1], however reconstruction of a putative DM signal with traditional frequency-domain analysis requires ``perfect data (i.e., regularly-sampled with no missing samples). In a real LBAI experiment, there will undoubtedly be imperfections in the data leading to downtime. This downtime can arise from operational considerations (e.g., maintenance), the operational environment (motion of people and animals [2] or elevators), and robustness of the experimental apparatus (e.g., bad atom launches). In this work, we investigate the impact of various downtime models on the overall DM sensitivity of such an experiment. We compare the sensitivity for each downtime model as determined by a ``compound FFT analysis to a baseline no-downtime case. We also show how much sensitivity can be regained by moving to a Lomb-Scargle frequency analysis, as in [2]. Furthermore, we demonstrate reconstruction of the DM wave s phase as well as its frequency. [1] D. Antypas, et al, ``New Horizons: Scalar and Vector Ultralight Dark Matter (2022). arXiv:2203.14915 [2] J. Carlton and C. McCabe, ``From RATs to riches: mitigating anthropogenic and synanthropic noise in atom interferometer searches for ultra-light dark matter (2023). arXiv:2308.101731

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

QES-Plume v1.0: a Lagrangian dispersion model

Low-cost simulations providing accurate predictions of transport of airborne material in urban areas, vegetative canopies, and complex terrain are demanding because of the small-scale heterogeneity of the features influencing the mean flow and turbulence fields. Common models used to predict turbulent transport of passive scalars are based on the Lagrangian stochastic dispersion model. The Quick Environmental Simulation (QES) tool is a low-computational-cost framework developed to provide high-resolution wind and concentration fields in a variety of complex atmospheric-boundary-layer environments. Part of the framework, QES-Plume, is a Lagrangian dispersion code that uses a time-implicit integration scheme to solve the generalized Langevin equations which require mean flow and turbulence fields. Here, QES-Plume is driven by QES-Winds, a 3D fast-response model that computes mass-consistent wind fields around buildings, vegetation, and hills using empirical parameterizations, and QES-Turb, a local-mixing-length turbulence model. In this paper, the particle dispersion model is presented and validated against analytical solutions to examine QES-Plume’s performance under idealized conditions. In particular, QES-Plume is evaluated against a classical Gaussian plume model for an elevated continuous point-source release in uniform flow, the Lagrangian scaling of dispersion in isotropic turbulence, and a non-Gaussian plume model for an elevated continuous point-source release in a power-law boundary-layer flow. In these cases, QES-Plume yields a maximum relative error below 6 % when compared with analytical solutions. In addition, the model is tested against wind-tunnel data for a uniform array of cubical buildings. QES-Plume exhibits good agreement with the experiment with 99 % of matched zeros and 59 % of the predicted concentrations falling within a factor of 2 of the experimental concentrations. Furthermore, results also emphasize the importance of using high-quality turbulence models for particle dispersion in complex environments. Finally, QES-Plume demonstrates excellent computational performance.

58 GEOSCIENCES↗
DOE OSTI2021

Sneutrino tribrid inflation, metastable cosmic strings and gravitational waves

Here, we present a successful realization of sneutrino tribrid inflation model based on a gauged U(1)_B-L extension of Minimal Supersymmetric Standard Model (MSSM). A single interaction term involving the B-L Higgs field and the right-handed neutrinos serves multiple purposes. These include the generation of heavy Majorana masses for the right-handed neutrinos to provide an explanation for the tiny neutrino masses via the seesaw mechanism, a realistic scenario for reheating and non-thermal leptogenesis with a reheat temperature as low as 10 6 GeV, and a successful realization of inflation with right-handed sneutrino as the inflaton. The matter parity which helps avoid rapid proton decay survives as a Z 2 subgroup of a U(1) R-symmetry. Depending on the choice of model parameters yields the following predicted range of the tensor to scalar ratio, 3 x 10 -11 ≲ r≲ 7x 10 -4 (6 x 10 -7 ≲ r ≲ 0.01), and the running of the scalar spectral index, -0.00022 ≲ dn_s/dln k ≲ -0.0026 (-0.00014 ≲ dn_s/dln k ≲ 0.005), along with the B-L breaking scale, 3 x 10 14 ≲ M/ GeV≲ 5 x 10 15 (6 x 10 15 ≲ M/ GeV≲ 2 x 10 16 ), calculated at the central value of the scalar spectral index, n_s =0.966, reported by Planck 2018. The possibility of realizing metastable cosmic strings in a grand unified theory setup is briefly discussed. The metastable cosmic string network admits string tension values in the range 10 -8 ≲ Gμ s ≲ 10 -6 , and predicts a stochastic gravitational wave background lying within the 2-σ bounds of the recent NANOGrav 12.5-yr data.

supersymmetry and cosmology↗