Probing Dark Sectors with Neutron Stars

Tensions in the measurements of neutron and kaon weak decays, such as of the neutron lifetime, may speak to the existence of new particles and dynamics not present in the Standard Model (SM). In scenarios with dark sectors, particles that couple feebly to those of the SM appear. We offer a focused overview of such possibilities and describe how the observations of neutron stars, which probe either their structure or dynamics, limit them. In realizing these constraints, we highlight how the assessment of particle processes within dense baryonic matter impacts the emerging picture -- and we emphasize both the flavor structure of the constraints and their broader connections to cogenesis models of dark matter and baryogenesis.


Introduction
Various astronomical observations speak to the existence of dark matter and dark energy [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16], and the detections of neutrinos from extraterrestrial sources reveal nonzero neutrino masses [17][18][19].Also, the cosmic excess of matter over antimatter [20,21], i.e., the baryon asymmetry of the universe (BAU), cannot be described within the SM, because the electroweak phase transition is not of the first order [22,23] for the empirically determined Higgs mass [24,25].Thus, physics beyond the SM has been established, even if its origins continue to be mysterious.Although these may come from dynamics that is disjoint from the SM, the possibility of identifying its nature mandates the careful consideration of any experimental anomalies that appear, for they could point to new dynamics.In this article, we consider the neutron lifetime anomaly [26,27], which entails the possibility that it could, at least in part, derive from new physics [28][29][30][31], and how that possibility can be constrained by observations of neutron stars [32][33][34][35][36][37].To set the stage, in this introduction we (i) place the neutron lifetime anomaly in the context of other low-energy β decay experiments, (ii) motivate why new physics in neutron decay could be linked to models that explain both dark matter and the BAU, and (iii) discuss how neutron dark decays connect to neutron star structure and dynamics.

Low-Energy Experimental Anomalies in Beta Decay
The expected pattern of quark flavor mixing under charged current interactions in the SM make the Cabibbo angle, i.e., the mixing of the first-and second-generation of quarks, central to a unitarity test of the Cabibbo-Kobayashi-Maskawa (CKM) matrix with elements V ij .This test probes the possibility of additional contributions to quark mixing and/or whether G F , the effective strength of the weak interactions as measured through the muon lifetime, is really universal.The so-called Cabibbo angle anomaly (CAA) has different components, and its resolution gains in importance, because possible non-SM effects need not be very small [38,39].There are disagreements in the assessment of V ud , which is most precisely probed through neutron decays or through superallowed (0 + → 0 + ) transitions in nuclei, as well as in V us , for which K → ℓν ℓ (K ℓ2 ) and K → ℓν ℓ π (K ℓ3 ) decay studies are the most precise.Theoretical corrections and inputs are needed in the analysis of the experiments to determine V ij .The determinations of V us from K µ2 and K ℓ3 decays do not agree [40][41][42], and the extracted values of V ud are also in slight tension with each other [42].The combined results are also at odds with CKM unitarity, with the differences predicated by the computed radiative corrections [43][44][45].In the superallowed case, the effect of the axial vector current enters in the radiative corrections, which can be modified according to nuclear structure as well [46].Recent developments in the theory of the radiative corrections in neutron β decay have led to refined assessments of their value, thus impacting the value of V ud [45]; we refer to [47] in this volume for a summary.It is worth noting that these determinations also differ from the value of |V us /V ud | determined from the ratio of rates in K ℓν and π → ℓν ℓ (π ℓν ) decays [41,42].These anomalies could have various resolutions, thereby arising from either limitations in the experiments themselves (there is only one high quality measurement of K → µ2) or in the SM theoretical inputs; however, particular new physics models, such as those with right-handed currents, could be at work [39,42].
The mixed symmetry nature of neutron β decay implies that at least two different experiments are needed to extract V ud , and our interests make it important to flesh out the details.In the SM, the neutron lifetime τ n is where λ ≡ g A /g V (the weak coupling constants of the nucleon), m e is the electron mass, f is a phase space factor, and η is 4908.6(1.9)s [44], 4903.6(1.0)s [45], or 4905.7(1.5)s [48], depending on the calculation of the radiative correction δ RC (and f ) used.Following common practice, we use λ and g A interchangeably, which is also strongly supported by an analysis of the radiative corrections to g A [49].The uncertainties in δ RC dominate those in η and have been used to determine its reported error.Thus, we see that the right-hand side of Equation ( 1) is strongly constrained, and the parameter λ is determined separately through measurements of the neutron β decay correlations, with the coefficients [50] termed as the "neutron alphabet [51]".Thus, to determine V ud , in order to combine with V us for a test of CKM unitarity, both the lifetime and a measured decay correlation are needed.The most precisely determined correlation coefficient is "A", which is the correlation between the neutron spin and the momentum of the decay electron; thus, that determination has been commonly used.Recently, precise determinations of "a", which is the correlation of the decay electron and antineutrino momenta, have been achieved.Thus, we shall compare the outcomes of combining the measurements of either the A or a correlation coefficients with the neutron lifetime for V ud .
In Figure 1, we show the data set that constitutes the neutron lifetime anomaly-we see that the measured lifetime τ n has been shown to be significantly different if the neutrons decay while stored in a bottle or trap, with the surviving neutrons counted after some storage time, or if the neutrons decay while in a beam within a device that counts the decay protons.The determined neutron lifetimes differ by some 9.5 s, which is a significant variation.One prominent impact this shift could have is on big bang nucleosynthesis (BBN), because the neutron lifetime determines, to a good approximation, the temperature at which the weak interactions freeze out in the early universe-and a shorter lifetime reduces the expected yield of 4 He [52].We refer to [53] for a review of the neutron lifetime's impact on element yields from BBN.Since the bottle lifetime is shorter, the difference between the two determinations can be reasonably attributed to the possibility of neutron dark decays [28] coming from physics beyond the SM, and we refer to the review [54] in this volume for further discussion.could not be explained with known theory [28].[55].In all neutron trap experiments, which use ultracold neutrons (UCNs), the neutron lifetime was determined by counting the neutrons that persist, whereas in the beam experiments, it was determined by counting the decay protons.The central values of the two sets of experiments disagree by 9.5 s, with a significance in excess of 4 σ-this is the neutron lifetime anomaly.The Particle Data Group (PDG) does not include the beam measurements in their fits [56].
Testing the V − A structure of neutron (and nuclear) β decay in the SM has been a focus of intense research for decades [51,[57][58][59].Here, the various decay observables, including the lifetimes and β decay correlations, can constrain the possibility that the neutron lifetime anomaly comes from new physics [36,60,61].That is, the agreement of the value of V ud with CKM unitarity and the value of λ in a neutron lifetime and in a decay correlation limit the possibility of non-V − A contributions, which could include neutron dark decays or the possibility of n − n ′ oscillations, to the lifetime.In that context, the PDG gives an average of |λ| = 1.2754 (13) with a scale factor of 2.7 from the A and a correlation coefficients [56], where the most precise determination, from an A measurement, is from PERKEO III |λ| = 1.27641 (56) [62].The other two most precise determinations of λ, which are also determined from A [63,64], are consistent with that result.Recent and precise a determinations [65][66][67] have emerged, however, that differ significantly from the A results.We note, e.g., |λ| = 1.2677 (28) [65]-and the latest status of the various measurements is shown in Figure 2.This discrepancy cannot yet be challenged by direct calculation in lattice quantum chromodynamics (QCD).As for the lattice QCD results, we note the Flavour Lattice Averaging Group (FLAG) average from simulations with N f = 2 + 1 + 1 flavors, g A = 1.246 (28) [68], as well as the CalLat-18 result of g A = 1.271 (13) [69], which is compatible with the PERKEO III determination, though it has been updated to g A = 1.2642 (93) [70].Placed in the context of a test of CKM unitarity, as shown in Figure 2 using the radiative corrections of [45], we see that the values of |V ud | from the bottle lifetime and the PERKEO III A measurement are entirely consistent with the value determined in superallowed decays.Interestingly, the value of |V ud | from the measured beam lifetime combined with λ from an a measurement is very comparable to its value determined using the bottle lifetime and the λ from A. In the presence of two possible solutions for λ and V ud that are similarly compatible with CKM unitarity, it would seem we cannot (yet) limit a dark decay solution to the neutron lifetime anomaly on the basis of this comparison.This perspective differs from that of earlier work [60,61], and it is driven by the recent a measurements.Ultimately, any solution to the neutron lifetime anomaly using new physics is required to be compatible with the measured V − A structure of the SM currents.1.This figure includes radiative corrections as per [45].Vertical bands illustrate |λ| measurements from aSPECT [65] (light green), aCORN [67] (light blue), and PERKEO III [62] (light purple).The measured value of V ud via superallowed nuclear transitions is shown by a red band [56].Additionally, lattice values of g A are shown comprising the 2021 FLAG average for N f = 2 + 1 + 1 flavors [68] and the 2019 results from the CalLat collaboration [70].

From Anomalies to Dark Matter and the Cosmic Baryon Asymmetry
Although astronomical observations speak to the existence of dark matter that is (i) stable or nearly so on Gyr time scales [13,14], (ii) possessive of no substantial strong or electromagnetic charge [71], and (iii) engenders an "inside-out" formation of large-scale structure [7], many different sorts of dark matter candidates are possible [72,73], and they can impact astrophysical observations at small scales differently [74].It is intriguing that the contributions of baryonic matter and of dark matter to the present cosmic density parameter Ω 0 are not too dissimilar, with the fraction in baryonic matter being some 20% of the dark matter one [13].Perhaps, then, dark matter is asymmetric [75,76]-much as baryonic matter is-thereby suggesting that the two effects could share a common origin.That possibility could be realized through a dark cogenesis model [77][78][79], and interestingly there is also a connection between such models and a dark decay model of the neutron lifetime anomaly [28], thereby yielding the dark decay n → χγ, where χ is a fermion with a baryon number of B = 1.This connection is concrete, even if the ultimate numerical impact of the n → χγ decay rate on the neutron lifetime anomaly finally turns out to be small.Dark cogenesis models can act in counterpoint to other popular mechanisms, such as via weakscale supersymmetry-in the latter, the severity of permanent EDM searches for the electron and neutron limit the possible phase space severely [80].In contrast, the B mesogenesis [78] mechanism and its variants [32,79,[81][82][83][84] do not require sources of CP violation beyond the SM to operate, and they are, moreover, testable mechanisms of baryogenesis [85].The schematic illustration of B mesogenesis is presented in Figure 3.The branching ratio for B → B χ decay cannot be too small [86], and experimental constraints exist from searches using BaBar data, with improvements expected from Belle-2 [87].To generate the noted decays, new scalars ought to exist as well, and they can be probed through collider searches and through constraints from flavor physics [79,88].Here, we discuss what further constraints come from precision studies of neutron stars and their dynamics [36,37].[78]: a late-scale dark cogenesis model.Here, a scalar Φ, produced out of equilibrium, decays and hadronizes to a B 0 q B0 q pair, with q ∈ d, s.The pair mixes in the presence of CP violation, and, finally, a B meson decays preferentially to B → B χ, which is a baryon B-a dark antibaryon χ pair.Since χ carries B = −1, the B χ final states do not break the baryon number, yet this scenario produces a baryon asymmetry in the visible sector.The produced χ(χ) decays to other dark sector states, thus preventing washout of the baryon asymmetry, with the particles of the dark sector final states serving as dark matter candidates.Different dark sector choices, that we note in text, differ in their impacts on the dynamics of a neutron star [37].

Constraints from Neutron Star Structure and Dynamics
Taken at face value, the difference in the beam and bottle lifetimes reflects roughly a 1% change in the neutron decay rate-and if the neutron does decay via n → χγ or n → χϕ, where ϕ is a scalar that could also decay via ϕ → e + e − [28], then the ability to form a protoneutron star is impacted [33,34,89].In particular, its maximum mass is only some 0.7M ⊙ , which is in conflict with observations of neutron stars of much larger masses.This constraint can be sidestepped in different ways, such as if χ self-interactions are included in the model, as also studied in [90] as an example.We note [91] in this volume for a discussion of neutron star constraints on dark scalar boson production in this context.Here, rather, we focus on the new physics channels with dark fermions, particularly with a baryon-number-carrying particle χ and the connections of this possibility to models that describe dark matter and baryon asymmetry.In this context, χ must decay, though not to antibaryons, to avoid a washout of the baryon asymmetry, and the size of the realized baryon asymmetry relies on the rate of B → χB.Here then, the neutron dark decay rate can be smaller and thus be in accord with neutron star observations.It thus would not explain the full size of the neutron lifetime anomaly, but constraints on its size are useful in that they can be interpreted in terms of limits on the flavor structure of the quark-dark baryon couplings.Here we note that although χ must decay to other dark particles, different dark sector choices can be made.The precise energy loss constraints we are able to make from the study of existing neutron stars are sensitive to those choices, and we will explain how in Section 2. There are also dark decay models, such as the n → 3χ model of [92], that we cannot constrain with the methods we review in this paper, though that model is constrained by the would-be agreement between the decay correlations and the bottle lifetime, should the anomaly in the a and A determinations of λ (or g A ) be resolved.
In what follows, we discuss, in turn, the effects of the existence of dark sectors on stars and, particularly, how neutron stars can constrain models of neutron dark decays, how computations of particle processes in dense matter are important to determining those constraints, and, finally, what constraints emerge from neutron stars.We then conclude with a perspective on the neutron lifetime anomaly and a final summary.

Stellar Constraints on Dark Sectors
Neutron stars, which are remnants of stellar evolution, present a unique opportunity to study the intricate interactions of dark matter within extreme astrophysical environments.In this section, we place this in a broader context, thus delving into the multifaceted role of dark matter across different stages of stellar evolution, from the early universe's first stars to the dynamic phenomena observed in supernovae.We explore how dark matter not only influences the life cycle of stars but also plays a crucial role in the formation, structure, and evolution of neutron stars.

Dark Matter's Role in Stellar Evolution
The role of dark matter extends beyond the confines of neutron stars.Within the so-called Λ CDM cosmological framework, cold dark matter (CDM) is a key driver of the universe's "inside-out" structural growth, thereby crucially shaping the formation of galaxies and stars from the earliest density fluctuations [6].Its influence is not limited to gravitational interactions: the dark sector potentially engages in a variety of interactions and phenomena that have a lasting impact on stellar evolution.This section explores these interactions, starting with the pivotal role of dark matter in the life cycle of the first stars.

Impact on First Stars and Population I Stars
In the earliest chapters of cosmic history, population III (Pop III) stars, the universe's very first luminous bodies, emerged from the gravitational collapse of pockets of gas in small halos driven by dark matter [93].These stars have been predicted to form at redshifts of z < 20-30 in halos with masses ranging from 10 6 M ⊙ to 10 8 M ⊙ .Recent reports suggest the possible identification of a Pop III star cluster in the remote fringes of a galaxy, thus highlighting ongoing efforts to deepen our understanding of these primordial objects [94].Beyond its gravitational influence, dark matter, through its annihilation channels, might have significantly impacted the formation of Pop III stars.Such annihilation events could act as an additional heat sources, thereby potentially inhibiting protostar formation by offsetting essential cooling mechanisms [95,96].In scenarios where a high dark matter density exists at the gas cloud's location, if some annihilation products are trapped and thermalized within the cloud, and with this dark matter-induced heating surpassing expected cooling effects, a novel stellar phase, known as the "dark star" phase, could emerge.This phase could manifest itself as a giant hydrogen-helium star (R ≳ 1 AU), which would be powered by dark matter annihilation instead of nuclear fusion [95].Furthermore, with the ongoing capture of dark matter, these objects may evolve into supermassive dark stars, with masses larger than (10 5 -10 7 )M ⊙ and luminosities exceeding (10 9 -10 11 )L ⊙ [97].Such objects could potentially be observable by the James Webb Space Telescope (JWST) [97,98].Among the four metal-poor, high-redshift (10.3 < z < 13.2) sources discovered by the JWST Advanced Deep Extragalactic Survey (JADES) [99,100], three have been identified as potential supermassive dark star candidates [101].
Dark matter can significantly influence the evolution of Pop III stars postformation, especially for those near the Eddington limit.The additional heating from dark matter annihilation could destabilize the most massive Pop III stars, thereby effectively setting an upper mass limit.For instance, the impact of nonthermal, super heavy dark matter particles (mass range 10 8 -10 15 GeV) is explored in [102].It has been found that at a dark matter density of about 10 16 GeV/cm 3 , conforming to the XENON1T one-year detection bounds [103], the maximum mass of a Pop III star could be restricted to a few tens of solar masses.In environments with extreme dark matter densities, such as 10 18 GeV/cm 3 , the growth of Pop III stars might be limited to just a few solar masses due to the Eddington limit.
Pop III stars, which are luminous yet transient features of the early universe, concluded their short lifecycles in supernovae, thereby enriching the cosmos with heavy elements.
However, the influence of dark matter could extend beyond the inception of these stars, thus significantly shaping the endpoint of their evolution.This influence is epitomized by the black hole mass gap (BHMG), which is a range delineated by 45 M ⊙ ≤ M BH ≤ 120 M ⊙ where no black holes emerge from the direct collapse of massive stars due to pair instability.
The upper boundary of the BHMG results from supermassive stars bypassing pairinstability supernovae to form black holes directly, while the lower threshold arises from less massive stars that undergo pulsational pair-instability supernovae, thereby shedding significant mass through pulsations and leading to lower-mass black hole remnants than their progenitors.Modifications to Pop III stellar physics, potentially through new physics such as light particle emission, can shift the BHMG boundaries.Such particles can expedite helium burning and diminish stellar wind mass loss, thereby consequently reducing the available oxygen for burning and suppressing the pulsational pair-instability mechanisms in Pop III progenitors [104].
Transitioning from dark matter considerations in the stellar evolution in the early universe, we now turn our focus to the influence of dark sector physics on more contemporary stars.Luminosity measurements from the Sun and horizontal branch (HB) stars are instrumental in constraining dark sector particles, which are known for their potential to alter stellar cooling rates.The Sun's emission of energy through a "dark channel", such as dark photons, is particularly restricted by solar boron ( 8 B) neutrino flux measurements [105].These observations establish a conservative upper limit on the Sun's dark luminosity, which is mandated to be less than 10% of the Sun's total photon luminosity [106].Furthermore, the emission of light particles such as neutrinos and axions from the cores of red giants can delay helium ignition.This results in the development of more massive cores, thus leading to a brighter phase prior to their evolution into HB stars [107].The tip of the red giant branch (TRGB) is especially sensitive to new energy loss mechanisms, including those prompted by axions [108], dark photons [109], or millicharged particles [110], which can cause noticeable deviations in TRGB brightness from their standard values.Conversely, dark matter annihilation within a star's core might trigger an earlier helium ignition, which could potentially lead to a dimming of the TRGB [111], thereby offering a contrasting perspective on the impact of dark matter within stellar environments.
Significant dark matter accumulation within stars can impact their structures by altering convective and radiative regions due to energy transfer by dark matter particles [112].These structural changes influence stellar pulsations, with implications for asteroseismology [113,114].In asymmetric dark matter (ADM) scenarios, where a primordial asymmetry determines the dark matter relic abundance instead of annihilation processes [115], ADM can accumulate without annihilation counterbalance.This accumulation leads to notable effects on stellar evolution, such as altered thermal transport in main sequence stars [116][117][118].Specifically, ADM interactions, depending on their nature, can modify the temperature gradient in a star's core by evacuating nuclear energy, with spin-dependent interactions playing a significant role during core hydrogen burning [119].

Influence on Supernovae Dynamics
The deaths of stars are also not free from potential dark sector influences.Specifically, the study of core-collapse supernovae, which mark the cataclysmic end of massive stars (with masses exceeding approximately 8 M ⊙ ) offers a unique window for probing possible hidden sector interactions with ordinary matter.The emission of light particles from supernova remnants introduces an additional cooling mechanism, which operates by draining energy from the protoneutron star (PNS).This process is subject to constraints, as it could lead to a reduction in the duration of the neutrino burst observed from SN 1987A, which is a pivotal discriminant among astrophysical observations [120][121][122].To maintain consistency with these observations, a conservative upper limit on the luminosity of dark particles (L d ) has been established, thus ensuring that it remains below the neutrino luminosity observed during the Kelvin-Helmholtz cooling phase, which takes place approximately one second after the core bounce.This constraint results in the formula-tion of a specific condition: L d ≲ 3 × 10 52 erg s −1 , which is imposed one second after the bounce [107].Supernova 1987A has also been instrumental in establishing limits on various extensions of the dark sector.These include constraints on axions [123], dark photons [124], and, notably, on muonic dark forces [125,126].The latter, particularly the gauged L µ − L τ model [127,128], has garnered interest for potentially resolving both the muon (g − 2) µ anomaly [129,130] and the Hubble constant tension [21,131,132].Moreover, SN 1987A has been critical in constraining hyperon decay processes, thereby encompassing both decays into light bosons [133] and those leading to GeV-scale dark sector particles [79] -this is a constraint that we revisit in our study of neutron dark decay constraints in neutron stars.
The exploration of supernovae, which is crucial for understanding the dynamics of dark sectors, culminates in the birth of compact star remnants, which leads us to our next topic.

Neutron Stars as Dark Matter Laboratories
Among these remnants, neutron stars stand out due to their extreme densities and strong gravitational fields, thus making them ideal laboratories for probing dark matter effects in dense matter [134,135].In this discussion, we categorize the effects of the dark sector on neutron stars into two distinct types: static and dynamic.
Static effects pertain to modifications in the equation of state (EoS) that governs the dense matter within neutron stars, with [136] as an explicit example.These alterations reshape the neutron stars' configuration space in hydrostatic equilibrium, as dictated by the EoS.Such effects reach their complete manifestation well within the age (its observation time) of the star, thus acting in a time-independent manner during our observations.Dynamic effects, in contrast, involve changes that occur over time, thereby impacting the neutron star's structure, such as its temperature and spin.These changes evolve with the star, thus reflecting ongoing transformations.
It is important to recognize that a new physics model can induce both static and dynamic effects on neutron stars.There are models that exclusively yield dynamic effects while still preserving the original EoS.In such cases, these models facilitate a quasi-equilibrium transformation of the star, thereby unfolding within a configuration space that is derived from, and remains consistent with, the unaltered EoS [36].

Static Effects of Dark Matter on Neutron Stars
Building on the concept of static dark effects, a notable example is the introduction of a dark baryon.This particle, comparable to a baryon in mass but with markedly weaker interactions, results in a decrease in the pressure within neutron stars.Such a phenomenon effectively "softens" the EoS.This softening consequently reduces the maximum stable mass for neutron stars (M TOV ) as determined by the Tolman-Oppenheimer-Volkoff (TOV) equations for hydrostatic equilibrium [137,138].Such modifications present significant challenges to the simplest models aimed at resolving the neutron lifetime anomaly [54], as we have noted in Section 1.3.They reduce M TOV to about 0.7 M ⊙ [33,34,89], which is markedly below the mass of the most massive observed neutron stars, which exceed 2 M ⊙ [139][140][141].To reconcile these discrepancies, the introduction of repulsive interactions either within the dark sector or between dark particles and baryons becomes necessary [33,34,89,91,142,143].Alternatively, a dark decay channel leading to a significantly lower equilibrium density of dark particles could circumvent this stringent bound [92].In contrast, the condensation of light scalar fields, such as axions, within neutron stars-induced by the star's baryonic content-results in a reduction of the baryon masses.This phenomenon can precipitate the emergence of a new ground state in nuclear matter, thereby potentially leading to an increased maximum mass for neutron stars [144,145].

Dark Dynamics of Neutron Stars
Dynamical effects on neutron stars can manifest themselves in various forms, ranging from abrupt, dramatic events to more gradual changes.A prime example of a sudden effect Universe 2024, 10, 67 9 of 33 is the collapse of a neutron star into a black hole.Such extreme transformations can be triggered by specific interactions with dark matter, as seen in three distinct scenarios.
In one scenario, certain weakly interacting massive particle (WIMP) variants can accumulate within neutron stars.Over time, these WIMPs may gravitate towards the star's center, thus eventually forming a self-gravitating core.This class of possibilities includes massive bosons lacking vectorial repulsive interactions and both bosons and fermions with strong attractive scalar interactions.Such a concentration of WIMPs could precipitate the rapid destruction of neutron stars in a period shorter than their observed lifespans [146].
In a parallel scenario, the accretion of nonannihilating dark matter within a neutron star can result in the formation of a degenerate dark core.If this core expands beyond its Chandrasekhar limit [147], the neutron star hosting it might undergo collapse [148,149].This process potentially leads to the creation of anomalously low-mass black holes within a mass range of approximately (1-2.5)M ⊙ .Such black holes, if they merge, could produce gravitational waves detectable by the LIGO-Virgo-KAGRA network.The absence of detections thus far provides constraints on the interactions of nonannihilating dark matter, as suggested by recent studies [150].However, the possibility of dark matter-induced core collapse can be highly sensitive to dark matter self-interactions [151].Even extremely weak repulsive self-interactions can prevent black hole formation, which is a scenario that is particularly relevant in the context of realistic bosonic dark matter models [152].Furthermore, the coannihilation of dark matter with nucleons, which is a common feature in many asymmetric dark matter models, can also act to inhibit black hole formation [152].
A third avenue explores the impact of dark matter realized as charged massive particles, with masses ranging from 10 2 GeV to 10 16 GeV.These particles can accumulate in collapsing protostellar clouds.After the formation of a neutron star, they fall to its center, thereby potentially forming a black hole on a time scale of the order of a year.This black hole could then grow rapidly through accretion, thus leading to the host neutron star's destruction on a similar time scale.The existence of old neutron stars suggests that these charged massive particles, if present, make up no more than 10 −5 of the dark halo's mass [153].
Gradual dynamical effects can also occur.The thermal history of neutron stars can deviate from the predictions of standard cooling theory [154].This deviation, potentially speaking to new heating mechanisms, has been supported by observations of several old warm neutron stars [155,156].These stars exhibit surface temperatures higher than those predicted by standard cooling models.They could arise from the inclusion of further SM effects [157,158], including nonequilibrium processes [159] such as rotochemical heating, which results from the spin-down behavior of pulsars [160].However, heating could also be incurred through dark matter scattering and capture [158,[161][162][163] or annihilation [164].Alternatively, dark matter effects could alter the structure of dense matter altogether through the appearance of axion quark nuggets [165] or of finite density-seeded vacuum instabilities [166,167].
There are various heating effects.Kinetic heating, resulting from dark matter capture, involves dark matter particles being gravitationally accelerated to velocities exceeding half the speed of light.As these particles scatter within neutron stars, they deposit significant kinetic energy, thus possibly heating the stars to infrared black body temperatures [161].Although the dynamics within neutron star cores, the primary region of such scattering, are not fully understood, the role of the crust has been emphasized in dark matter interactions.Scattering with just the neutron star's low-density crust can still kinetically heat the star to infrared temperatures, which are detectable by forthcoming telescopes, thereby highlighting the crust's importance in dark matter detection and annihilation processes [168].Moreover, inelastic dark matter could significantly influence neutron star heating.Its interactions would be suppressed in standard nuclear recoil direct detection experiments, where small momentum transfers are insufficient relative to the mass splitting between dark matter states.However, neutron stars enable effective inelastic scattering, even with large mass splittings, due to the acceleration of dark matter particles to velocities approaching the speed of light during their infall [162,169].Additionally, the annihilation of WIMPs inside neutron stars can impact their thermal properties.This is particularly evident in neutron stars older than 10 million years that are situated nearer to the Galactic center, where the density of dark matter is expected to be significantly higher than in the vicinity of the Earth.For instance, in regions where the dark matter density is at least 3 GeV/cm 3 , such interactions could lead to an asymptotic internal temperature of around 10 4 K [164].Conversely, in regions with comparatively lower dark matter densities, other heating mechanisms may play a more dominant role in shaping the thermal evolution of neutron stars [157].
Studies of young neutron stars also give constraints.Temperature observations of NS1987A, the remnant of SN1987A [170,171], and of Cassiopeia A (Cas A), the secondyoungest known neutron star with an approximate age of 340 years, were shown to impose constraints on anomalous cooling processes in neutron stars.These observations play a crucial role in constraining the emissions of novel bosonic particles that could serve as additional cooling mechanisms, thereby altering the thermal evolution of stars [172,173].
As a separate tack, binary pulsar systems offer unparalleled insights into both gravitational physics and particle physics, particularly in the context of dark sectors.These systems enable precision tests of fundamental theories, thereby providing constraints on new physics beyond the SM.In gravitational physics, binary pulsars have been instrumental in rigorously testing general relativity and alternative theories of gravity in strong-field regimes [174,175].Notably, the Hulse-Taylor pulsar (PSR B1913+16) [176,177], in a binary system of a pulsar and a neutron star, and the double pulsar system (PSR J0737-3039A/B) [178,179] have played key roles in validating general relativity's predictions [180,181].Furthermore, the sensitivity of binary pulsar orbital periods to energy loss mechanisms [182][183][184] makes them ideal for investigating particle interactions that extend beyond the SM [185].Observations of these periods have led to constraints on ultralight dark matter candidates [186][187][188][189] and limits on dark baryon processes involving neutron-mirror neutron mixing [190], as well as broader constraints on various forms of the quasi-equilibrium baryon number violation (BNV) [36,37].The implications of these constraints, derived from binary pulsar period observations, on quasi-equilibrium BNV and specific dark decays of baryons are further explored in Section 4.
The precision of pulsar timing has become an essential tool in astrophysics, especially through pulsar timing arrays (PTAs), which have been used to find evidence for a stochastic gravitational wave background [191][192][193][194].This technique also plays a pivotal role in exploring ultralight dark matter (ULDM), as ULDM candidates with masses around 10 −22 eV have been predicted to create oscillating gravitational potentials at nanohertz frequencies [195].Such oscillations lead to periodic variations in the arrival times of radio pulses from pulsars, thereby overlapping with the frequency range for gravitational waves detectable by PTAs [196].The European Pulsar Timing Array has effectively used this aspect to constrain scenarios in which dark matter interacts only gravitationally, thereby demonstrating that ultralight particles in the mass range 10 −24.0 eV ≲ m ≲ 10 −23.3 eV cannot fully account for the observed local dark matter density [197].Additionally, pulsar timing exhibits sensitivity to the effects of heavier dark matter candidates, including primordial black holes and compact subhalos, thereby extending its utility in constraining a wide mass range of dark matter candidates [198].In addition, pulsar timing shows potential for probing BNV [199].Anomalies in pulsar spin-down rates and braking indices could be indicative of BNV processes, thereby possibly leading to pulsar spin-up behavior and the reactivation of dormant pulsars.With current timing precision being sensitive to the BNV rates of Γ∼10 −9 yr −1 , PTAs are increasingly becoming a crucial tool for future BNV detection or constraints, particularly as their precision improves in the context of gravitational wave detection.

quasi-equilibrium Dynamics of Baryon Loss in Neutron Stars
In Section 2.2, we explored the ways in which dark sector extensions to the SM can impact neutron stars.Building on this foundation, this section delves further into these phenomena, with a particular focus on the implications of specific dark decays of baryons within neutron stars.
One intriguing example involves neutron dark decays, such as n → χγ.These should occur at a rate of around 10 −5 s −1 to resolve the neutron lifetime anomaly.If the dark decay products, such as χ, remain stable over short timescales, they could swiftly saturate the interior of a neutron star shortly after the collapse of its progenitor star's core.This process is expected to drive dynamical effects whose outcomes would be chemical equilibrium between the decay products and the star's baryons.Modeling the EoS under these circumstances necessitates incorporating the new degrees of freedom, such as χ and other dark sector particles, affecting the star's pressure and energy density.The ground state of this system would be determined by the lowest energy state, subject to charge conservation and an extended baryon number concept, where B(χ) = 1.Should χ maintain stability over long timescales, the neutron star would eventually settle into a static state.This scenario illustrates a static dark sector influence, thereby altering the star's configuration space via its EoS.In this case, static limits such as the maximum mass of neutron stars can be used to constrain the dark decays.
If neutron decay products such as χ further decay or engage in annihilation, leading to particles that escape the star, the situation becomes dynamically complex.It is important to note that these depletion processes need not be direct; they may involve a series of reactions.For instance, if χ decays into leptons, then baryon number conservation would be explicitly violated.This, in turn, would trigger a sequence of baryon-conserving reactions within the neutron star.The resulting depletion of baryons manifests itself not only through neutrino emission but also through contributions to the thermal energy and work done on the stellar fluid.In this case, continuous neutron decays could dynamically influence the neutron star properties over its lifetime.Ongoing processes like these could lead to observable changes in neutron stars, including modifications in thermal evolution, variations in pulsar spins, or alterations in the orbital periods of binary pulsars.In this context, neutron stars offer an observational platform for constraining and potentially revealing active dark decay processes.
Broadening our perspective, the overall influence of baryon dark decays, which should dominate dark baryon processes in neutron stars because of the sheer number of neutrons involved, depends on the production (Γ P ) and depletion (Γ D ) rates of dark particles, as well as the sensitivity (σ obs ) of neutron star observations.While the depletion mechanisms themselves vary-manifesting as either the direct escape of dark particles from the star or an explicit BNV-the crucial factor in modeling these effects is their rate, with one notable exception.The thermal evolution of a neutron star does require a more detailed examination of the depletion mechanism, since it is influenced by the specific sequence of reactions involved in the process.These diverse scenarios necessitate distinct computational approaches, as depicted in Figure 4.
In cases where depletion (Γ D ) exceeds production (Γ P ), the star primarily comprises SM particles.Here, Γ P becomes the effective BNV rate to compare against σ obs .Conversely, if production (Γ P ) surpasses depletion (Γ D ), dark sector particles must be included in the EoS.The depletion rate (Γ D ) then serves as the effective BNV rate, because this determines the time scale over which the star changes.
In scenarios in which the neutron star observations are not sensitive to the BNV rate (Γ BNV < σ obs ), we have the following: 1.
Production dominance (Γ P > Γ D ) leads to static effects via EoS modifications, thus incorporating dark particles.
Conversely, when observations are sufficiently sensitive to possible BNV effects (Γ BNV > σ obs ) and under the assumption that in degenerate cold neutron stars the hydrodynamic response rate exceeds the chemical rate (Γ hyd > Γ chem ), we have the following: I.
The star is in both chemical and hydrostatic quasi-equilibrium if Γ BNV < Γ chem .This state warrants a quasi-equilibrium analysis, as developed in [36].II.The star remains only in hydrostatic quasi-equilibrium if Γ chem < Γ BNV < Γ hyd .A hydrostatic equilibrium description is valid on short timescales (δt < Γ −1 BNV ).However, a complete reaction chain must be considered, as chemical equilibrium is absent in the EoS.III.The absence of both equilibria when Γ BNV > Γ hyd requires comprehensive simulations of the star's hydrodynamic and chemical evolution.A departure from chemical and hydrostatic quasi-equilibrium necessitates numerical simulations for precise stellar modeling, with the expectation that stars transition from non-equilibrium to a quasi-equilibrium state over time.
space via its EoS.In this case, static limits such as the maximum mass of neutron stars can be used to constrain the dark decays.
If neutron decay products such as c further decay or engage in annihilation leading to particles that escape the star, the situation becomes dynamically complex.It's important to note that these depletion processes need not be direct; they may involve a series of reactions.For instance, if c decays into leptons, it would explicitly violate the baryon number conservation.This triggers a sequence of baryon-conserving reactions within the neutron star.The resulting depletion of baryons manifests not only as neutrino emissions but also contributes to the thermal energy and work done on the stellar fluid.In this case, continuous neutron decays could dynamically influence neutron star properties over its lifetime.Ongoing processes like these could lead to observable changes in neutron stars, including modifications in thermal evolution, variations in pulsar spin, or alterations in the orbital periods of binary pulsars.In this context, neutron stars offer an observational platform for constraining and potentially revealing active dark decay processes.
Broadening our perspective, the overall influence of baryon dark decays on neutron The potential pathways and resultant states stemming from neutron decays to dark final states in neutron stars, which are dictated by production (Γ P ) and depletion (Γ D ) rates.Observational sensitivity (σ obs ) and the star's response rates-chemical (Γ chem ) and hydrodynamic (Γ hyd ) ones-further influence these pathways.A departure from chemical and/or hydrostatic quasiequilibrium (denoted by ✗) necessitates numerical simulations for precise stellar modeling, with the expectation that stars transition from a nonequilibrium to a quasi-equilibrium state over time.
Notably, in the quasi-equilibrium scenario (I.), baryon number conservation on short timescales forms a key constraint in the derivation of its EOS.The primary distinction between production-and depletion-dominated states lies in the inclusion or exclusion of dark sector elements and the extended definition of the baryon number.
The quasi-equilibrium framework posits that if dark sector perturbations, such as baryon decays, proceed at a slower rate than the star's responses, the star will evolve within an equilibrium configuration space defined by the EoS.In scenarios where dark particle production dominates, this EoS incorporates the pertinent dark degrees of freedom.Conversely, in cases where depletion prevails, the EoS remains akin to the standard formulation, which is exclusive of dark sector influences.This approach effectively balances the dynamical aspects of dark sector interactions against the inherent stabilizing mecha-nisms of neutron stars.The computational simplicity in depletion-dominated scenarios is significant, as the EoS needs not reflect intricate dark sector details, thus facilitating a model-independent analysis.
Building upon this framework, the configuration space of neutron stars can be effectively parameterized by a set of n independent parameters, which is denoted as {Ξ 1 , Ξ 2 , . . ., Ξ n }.Variations in the neutron star observables, symbolized as O, can then be quantified through the following expression: where each term ∂O/∂Ξ i represents the sensitivity of the observable to changes in the ith parameter, as governed by the EoS.The rates of change for these parameters, dΞ i /dt, are determined by imposing constraints on the system's quasi-equilibrium dynamics.These constraints may include theoretical principles, such as the conservation or specific rate of violation of the baryon number, alongside observational data like the spin-down rates of pulsars.Together, these factors provide a comprehensive framework to model and understand the evolution of neutron star observables within the prescribed parameter space.Before exploring specific examples, we think it crucial to recognize that certain neutron star characteristics, such as its surface temperature, rotational period, or its magnetic field, may not always be central to our framework, particularly when these factors do not significantly influence the observable under study.Neutron stars, which are composed of cold, degenerate matter, experience minimal thermal contributions to their structure.Similarly, rotational effects are often minor for pulsars with longer periods, generally those exceeding a few milliseconds [200].Therefore, dynamical evolution within the configuration space can often be treated as independent from changes in the magnetic field, temperature, or rotation period.However, it is important to note that these elements should be explicitly included in the analysis when they directly impact the observable being studied.For example, when investigating pulsar spin-down behavior, rotation remains a crucial component within the formalism.Similarly, temperature may be a vital factor in processes that are subject to Pauli blocking mechanisms.
We now illustrate the application of the quasi-equilibrium framework with three distinct examples, each of which is succinctly summarized in Table 1.In each case, we note what configuration space is modified through the dynamical change indicated, as well as the constraints that act on those dynamics.Table 1.quasi-equilibrium frameworks for neutron star evolution under various scenarios.Parameters include the central energy densities for the entire star (E c ), visible sector (E V c ), and dark sector (E D c ), as well as the central angular velocity (ω c ).The observed spin-down rate is indicated by Ω, while the rate of baryon number violation or conversion in dark sector processes is denoted by Γ.

Scenario Configuration Space Constraints & Model Inputs
Pulsar Spin-Down Formation of a Dark Core without Depletion: • Configuration Space: This is defined by Ξ = {E V c , E D c }, which represents the central energy densities of the visible (V) and dark (D) sectors at the star's core.The formation rate is generally assumed to be comparable to or slower than the star's age.This rate does not result in thermalization between the two sectors.Consequently, the sectors are characterized by distinct energy densities and pressures, thereby interacting primarily through gravitational influence.
• Constraints: The total baryon number conservation ( ḂD + ḂV = 0) acts as a constraint, thereby ensuring a balanced rate of baryon conversion between the dark and visible sectors.The conversion rate from visible to dark baryons ( ḂD ) is derived from the specifics of the particle physics governing the decay or conversion processes.

3.
Baryon Dark Decays Followed by Efficient Depletion: • Configuration Space: This space, also known as the single-parameter sequence (illustrated in Figure 5), is singularly defined by Ξ = {E c } (if the star is not rapidly rotating).Here, E c denotes the central energy density of the neutron star.This simplified configuration space reflects the scenario where baryon dark decays and subsequent depletion processes dominantly influence the neutron star's structure.

•
Constraints: The evolution is primarily governed by the total baryon loss rate Ḃ.This rate is instrumental in dictating the evolution of the star's core composition and overall structure, as indicated by the green arrow in Figure 5, which represents the quasi-equilibrium evolution along the single-parameter sequence.Transitioning from the comprehensive analysis of neutron stars' quasi-equilibrium dynamics, we now shift our focus to specialized dark sector models, particularly those that expand neutron dark decay scenarios to encompass dark cogenesis.Such models propose a unified origin for the observed baryon asymmetry and a corresponding asymmetry within the dark sector.This approach intriguingly links the physics of neutron stars to wider cosmological narratives.

Illustrative Dark Sector Models
We have considered neutron decays in which a dark baryon χ appears, and the manner in which we can probe this possibility in a neutron star depends on if and how χ decays to other dark sector particles.In this section, we consider just one particular dark sector interaction, but its impact on a neutron star can vary significantly depending on the assumed spectrum of particles in the dark sector.Indeed, with suitable choices, we will see that it can sample all the scenarios outlined in the previous section.
Although we certainly expect χ to decay were it to play a role in mesogenesis, it need not-perhaps χ is stable and a possible dark matter candidate in and of itself.This would require, of course, that the χ is light enough that it cannot decay to light antibaryons, and we presume it does not decay just to mesons and/or leptons alone, so that true baryon number-violating (|∆B| = 1) processes do not appear in this context.Thus, as in [78,79], we introduce two new particles: a scalar baryon ϕ B with B = 1 and a Majorana fermion ξ.It is important that neither ϕ B nor ξ decay, nor that they transform into each other.However, χ may decay via [78] L dark ⊃ y d χϕ B ξ + h.c., where either ϕ B or ξ or both could be dark matter candidates.This interaction can destabilize the proton, produce χ decay, or produce χχ annihilation, with examples as shown in Figure 6a,b.For different choices of the particles' masses, different processes appear.
For example, we have the following: If χ does not decay, the neutron dark decay rate is constrained, or alternatively, χ selfinteractions must be present, as previously noted in Section 1, to align with observed neutron star masses.This scenario corresponds to the production dominance route depicted in Figure 4, thereby leading to static effects such as the modification of the maximum theoretical neutron star mass.However, with a suitable choice of dark sector masses, subsequent to neutron decay to a final state with a χ particle, either ϕ B and/or ξ can appear.These outcomes can act with different effects within a neutron star.Either sort of particle could fall to the core of the star and modify its evolution-and possibly act with a destabilizing effect [149].Alternatively, if ϕ B and ξ are light enough, they can escape the star.Thus, the scenarios described in the previous section can all be sampled through the dark sector interaction of Equation (3).In our own investigations [36,37], we have used a quasi-static approach to determine the evolution of the star in the presence of a slow apparent BNV to dark final states.The consistency of this analysis, that new physics dynamics does not impact the evolution of the star, requires that we have an effective mechanism to remove energy from the star upon decay to a dark final state.Accordingly, our analysis is predicated on the assumptions inherent in the depletion dominance pathway of Figure 4, which culminates in a quasi-equilibrium evolution of the star.As a concrete example of such, we have supposed that M ξ ≫ M ϕ B , M χ , thus realizing χ removal through χχ → ϕ B ϕ B decay, with ϕ B being light enough to escape the star.
Additional physical constraints appear if we work within the context of a dark cogenesis model.For example, if ξ → ϕ B χ could occur, and it would under the conditions we have just outlined for quasi-static evolution of the star, then if χ is heavy enough, χ → pπ − and χ → nπ 0 can also appear.In a broader context, such processes would act to dilute a cosmic baryon asymmetry.Thus, either the flavor structure of the pertinent quark-χ couplings must be sufficiently small, or the mass of χ is further constrained.Following the latter path, we note that successful baryogenesis within the dark sector scenario of Equation (3) (and masses such that χχ → ϕ B ϕ B occurs) can be realized if m χ ≤ 1.07784 GeV.This and nuclear stability considerations gives us the following window on χ [37]: though in what follows, we consider neutron star constraints on χ-baryon mixing over a larger mass window in anticipation of other model building solutions.We now turn to the evaluation of these decays within the dense medium of a neutron star.

Neutron Stars and Particle Decays in Dense Matter
The consideration of dense matter phenomena is crucial for investigating baryon interactions in neutron stars, where intense gravitational forces compress baryonic matter to densities exceeding nuclear saturation.Such extreme conditions push the system into the strongly repulsive domain of nuclear interactions, thus resulting in an elevation of the ground state energies of the baryons.This increase in energy levels potentially allows for decay processes that are otherwise kinematically forbidden in both a vacuum and within nuclei [37].Consequently, neutron stars present themselves as unique astrophysical laboratories, thereby offering the possibility to observe baryon dark decays that bypass the constraints usually encountered in earthbound experiments such as Super-Kamiokande or KamLAND [203,204].To gain a more comprehensive and quantitative insight into this phenomenon, we delve into a particular formalism of dense matter.This exploration aims to uncover the unique role of heavy neutron stars as cosmic laboratories, thus offering new probes of dark sectors.
In this vein, a hadronic relativistic mean field (RMF) theory offers a robust framework for addressing dense matter interactions in such extreme environments [205][206][207][208]. Within RMF theory, baryon interactions are mediated through various Lorentz scalar and vector mesons.The Euler-Lagrange equations that govern the dynamics of these meson fields include baryon currents as source terms, thus taking the form of ψψ for scalar mesons and ψγ µ ψ for vector mesons.As the baryon density increases, these currents intensify, thereby leading to a significant vacuum expectation value (VEV) for the meson fields.Within this framework, the mesonic degrees of freedom are thus effectively represented by their classical VEVs, which is a simplification that sees scalar meson VEVs contributing to the baryon masses and vector mesons to their four-momentums.This theoretical construct provides the foundation for understanding baryon behavior in the dense, locally uniform medium of neutron stars.
In such a medium, the wave function for a baryon with a canonical four-momentum p µ is expressed as ψ(x) = e −ip•x u(p * , λ), (5) where is solved in the Dirac-Pauli representation as where⃗ σ are the Pauli matrices, and η λ is the Pauli spinor with η ↑ = (1, 0) T and η ↓ = (0, 1) T .The baryons' energy spectrum (p 0 ) is delineated by where we note that in the RMF approximation, Σ µ and m * are density-dependent but independent of p µ .The CM energy of the baryons, denoted as E cm (p), is a function of their effective mass m * and momentum ⃗ p in the nuclear matter frame.It is expressed by the following equation: This energy reaches its minimum when the momentum |⃗ p| is zero and reaches its maximum at the Fermi momentum.The numerical values of m * , Σ 0 , and E cm for the baryons within neutron stars under the DS(CMF)-1 EoS [202], which comes from a more sophisticated descendent [208] of the RMF model we have noted, are depicted as functions of density in Figure 7. (We anticipate further evolutionary development of the RMF approach through a determination of its parameters from a description of pure neutron matter within chiral effective theory [209].)As evidenced in Figure 7, a notable reduction in the effective baryon masses (m * ) at higher densities is observed, which is indicative of chiral symmetry restoration in these high-density regimes in this model.Furthermore, the dynamics between m * and Σ 0 (in the nuclear matter frame) are such that while m * decreases, Σ 0 increases, thereby cumulatively leading to an increase in the total energy of the baryons as delineated by Equation ( 8) at higher densities.The core (maximum) density of four pulsars (in binary systems) are denoted by vertical dot-dashed lines for comparison.This sets the stage for understanding how each neutron star, with its unique density profile, represents a distinct astrophysical laboratory.Figure 8 further demonstrates how each neutron star offers a unique probe that is intrinsically linked to its specific density profile.The CM energy spectrum of the baryons (Λ, neutron (n), and proton (p)) in neutron stars in PSR J1614-2230, PSR J0348+0432, and PSR J0737-3039A has been observed to surpass their vacuum rest mass within the stars' cores, as shown in the figure.Notably, heavier neutron stars such as in PSR J1614-2230 and PSR J0348+0432 have been shown to exhibit a higher energy density reach, thereby underscoring the relationship between a star's mass and its potential for exploring baryon dark decays.This characteristic of neutron stars, particularly the more massive ones, highlights their capability to facilitate baryonic decays into states that are kinematically forbidden in vacuum conditions.The unique energy reaches offered by these pulsar systems offer invaluable insights into dark matter models beyond the scope of terrestrial experiments.
In the ensuing section, we demonstrate how observations of binary pulsars, in conjunction with the considerations of dense matter outlined herein, pave the way for novel constraints on the dark decays of baryons.We have chosen to represent only PSR J0737-3039A here and henceforth, as it is nearly identical in mass to PSR J0737-3039B.The plot highlights how the CM energy in high-density regions exceeds the baryons' vacuum rest mass, thus hinting at the possibility of decay into states that are kinematically forbidden in vacuum.Detailed analysis of these phenomena is provided in the main text.

Baryon Loss Constraints from Energy Loss Limits in Pulsar Binaries
Binary pulsar orbital periods, which are sensitive to energy loss mechanisms, present an effective observational tool for probing quasi-equilibrium BNV processes.In scenarios with efficient depletion, such as Example 3 in Section 2.3, the quasi-equilibrium BNV at a rate Ḃ contributes to an energy loss Ṁeff from neutron stars as follows [36]: where ∂ E c represents the partial derivative with respect to the central energy density of stars, M is the mass, I denotes the moment of inertia, and Ω is the angular spin frequency.This energy loss affects the binary pulsar orbital period [182][183][184]: The precision in these terms determines the sensitivity of the BNV limits, thus making binary pulsars with accurate timing ideal for realizing stringent constraints.Moreover, as noted in Section 3, each pulsar, especially the more massive ones, offers a unique opportunity for the study of the BNV.Their larger mass extends the density reach, thereby enabling constraints on higher mass ranges of dark particles such as χ in processes such as n → χγ.Furthermore, selecting binary systems with heavier pulsars is essential for constraining the dark decays of hyperons.In EoSs that incorporate hyperon degrees of freedom, hyperons start populating the medium only in extremely dense environments [210], specifically in the cores of heavier neutron stars (as illustrated in Figure 8).In our pursuit of comprehensive coverage of the model parameter space, a strategic combination of binary pulsars-both heavy and precisely measured-is essential.Ideal candidates are systems with accurately determined masses, orbital periods, and wellquantified contributions to each term in Equation (12).It is critical to select binary systems that are free from mass transfer between their components to avoid complex alterations in the observed orbital period decay rate, as mass transfer can significantly influence Ṗobs b [211].Excluding the heavy candidate PSR J0952-0607, a "black widow" pulsar with substantial mass uncertainty M p = (2.35 ± 0.17) M ⊙ and significant rotational effects on its structure [139], our attention turns to the following candidate binary pulsar systems: . We ultimately selected PSR J0348+0432 and PSR J1614-2230, which offer the advantages of smaller errors in M p and Ṗobs b , thus ensuring more precise constraints.Complementing this selection, we also considered pulsar systems that are renowned for their high-precision measurements.Among these, the double pulsar J0737-3039A/B [179] and the Hulse-Taylor pulsar system [176] stand out.Since both are characterized by similar mass profiles, we focused on the double pulsar J0737-3039A/B, which was chosen for its superior precision.
The constraints derived from our selected pulsar systems were applied to a specific model involving baryon decays B → χγ, which was dictated by an effective mixing parameter ε Bχ .The in-medium Lagrangian governing these decays is represented as: As depicted in Figure 9, the pulsars provided constraints on models with dark baryon χ masses extending several hundred MeV above the mass of the decaying baryon.Specifically, the heaviest pulsar, in PSR J0348+0432, demonstrated the broadest energy reach, while the most precisely measured system, the double pulsar J0737-3039A/B, yielded the most stringent limits for lower χ mass ranges.By varying the choice of EoS, we can explore the robustness of our assessments.We have explored this for the entire family of DS(CMF) EoSs: DS(CMF)-1-8 [202,215] in [37].Interestingly, the combined ε nχ constraint changed negligibly below ∼1 GeV, thus effectively covering our window of interest, Equation (4), whereas variations in ε Λχ of a factor of a few appeared over the entire mass range studiedand only DS(CMF)-1,3,7 supported the appearance of strange baryons for the parameters of our candidate pulsars.We note that the ϵ nχ parameter at 2σ was bounded to be no larger than 2 × 10 −17 for the entire window in the χ mass we have supposed in a dark cogenesis scenario.
We consider the implications of these results in the next section.

Perspectives on the Neutron Lifetime Anomaly and Dark Cogenesis
We now determine the outcomes from our neutron star analysis in regard to the in-vacuum limits on neutron dark decay.We can translate the constraints on ε Bχ that we report in Figure 9 to ones on the in-vacuum decay rate of B → χγ using the following equation: We note that g B signifies the baryon magnetic moment (Landé) g factor, where g n ≃ 3.826 for neutrons, and g Λ ≃ −1.226 for Λ baryons [56].The resulting constraints on the vacuum branching ratios are illustrated in Figure 10.
In Figure 11, we focus on the particular range of masses pertinent to the neutron lifetime anomaly, thus reporting both ϵ nχ and B(n → χγ).It is apparent that terrestrial searches for baryon number-violating processes in underground detectors give the most stringent limits on the possibility of n → χγ, yet the windows on the χ mass open in such searches did not overlap with the possible mass window for the neutron lifetime anomaly-the nuclear stability constraint we have noted yielded the lowest possible χ mass, reported in Equation ( 4), and this precludes that possibility.Nevertheless, we include the KamLAND [204] and Super-Kamiokande [203] limits, as well as the limit on B(Λ → χγ) from BESIII [216] and that from the SN1987A study of [79] for context.Specifically, the KamLAND constraints are relevant to our study of the decay n → χγ; in these scenarios, the prompt photon would not meet the correlation cuts, thus rendering it invisible.We observe that our neutron star limit is about 10 18 times more severe than the accelerator limit and is similarly less sensitive to the BNV than the underground experiments.Additionally, we incorporated constraints arising from the influence of Λ decay on the neutrino signal duration observed from SN1987A [79].These constraints have been derived from a supernova simulation involving an 18.6 M ⊙ progenitor leading to a 1.553 M ⊙ remnant [217].Theoretically, these constraints could extend to lower χ masses as well.Notably, our derived limits significantly surpass the bounds established by the supernova data [79].
The outcome of our analysis appears to exclude a new physics interpretation of the neutron lifetime anomaly, but the context in which we performed this analysis is crucial: different choices within the dark sector could lead to alternate outcomes, as explored in Section 2.3.In our case, we have chosen a dark sector content and interaction as per Equation (3) and the accompanying discussion so that the SM dynamics of the neutron star determine the evolution of the star under a change in its baryon number-with the injected energy from the exotic neutron decay, particularly from χ production, removed through χχ → ϕ B ϕ B decay, with ϕ B subsequently removed from the star.Our approach aligns with the depletion dominance route depicted in Figure 4, thus resulting in a quasi-static evolution of the neutron star.However, in alternative scenarios, the establishing limits would depend on varying pathways of dark decay, as outlined in the same figure.
As for dark cogenesis models, to which our neutron dark decay studies connect, such as B mesogenesis [78], we constrain the flavor structure of the quark-χ couplings through the limits we report [37]-and we show our limits in our favored mass region according to the shaded band in Figure 9.   .Two-σ exclusion limits on the vacuum branching fraction for B → χγ processes are depicted in panels (a) for neutrons and (b) for Λ baryons.The green region represents limits from Super-Kamiokande [203].We highlight that the KamLAND constraint (red) [204] is pertinent to our study of n → χγ decay, as the prompt photon would not pass the correlation cuts and would remain undetected.Panel (b) also presents limits from neutron star studies for the DS(CMF)-7, DS(CMF)-3, and DS(CMF)-1 EoSs in order of increasing severity, along with constraints from the BES-III experiment (blue) [216], noting that the γ in this case would also be undetected, and the SN1987A limits (gold) [79]. .Exclusion limits at 2 σ on (a) ε nχ and (b) the vacuum branching fraction for n → χγ from our quasi-static neutron star analysis in the particular χ mass region pertinent to an explanation of the neutron lifetime anomaly.Additional constraints and expected limits (dashed) as per [218] have been included.This figure has been taken from [37].

Summary
There is an increasing recognition that the solution to the dark matter problem may involve an entire sector of dark particles with dark dynamics.The presence of such dynamics can be probed in terrestrial experiments that are sensitive to the determination of missing momentum [219] in hadron decays, and in this article, we have considered this broader possibility in the context of neutron stars.Particularly, we have noted the consequences from limits on the observation of energy loss in neutron stars, especially those that emerge from the timing of the period of binary pulsars on baryon dark decays, specifically of neutrons and Λ baryons.
Neutron star observables have multiple points of contact with the problem of the neutron lifetime anomaly and with the broader possibility of the simultaneous generation of dark matter and the BAU.We have noted models that would give rise to n → χγ or to all dark decays, as described in [54,91] in this volume, which are limited by the need to form neutron stars of sufficient mass.An unwanted softening of the neutron EoS through such neutron decay effects can be mitigated through the possibility of χ self-interactions or by reducing the neutron dark decay rate.We have studied the possibility of using dark sector interactions to eradicate the χ particles produced via dark decays through χχ annihilation to light, which are dark sector particles that can escape from the star.Here, we have studied the limit in which the dark decays proceed slowly with respect to the usual Urca processes that operate in the neutron star, thereby enforcing chemical and thermal equilibrium.Thus, in this quasi-static framework we are able to study the consequences of the observed energy loss constraints from pulsar timing studies on dark decays in the neutron stars.The constraints we have found thus far are severe, but broader possibilities remain-and they beckon.Certainly, there are models of the neutron lifetime anomaly that are not constrained by our concerns.Here, we note the model that yields n → χχχ [92].That model does not soften the EoS unduely, nor does it give rise to substantial energy loss effects that could be constrained through our pulsar timing studies.Thus, the possibility of a new physics solution to the neutron lifetime anomaly remains, with the future neutron decay correlation and lifetime studies offering the possibility of setting the controversy to rest, though effects in neutron stars, as well as in isolated ones, through observed heating or spin effects may finally serve to identify the effects of dark sectors in a positive way.
Author Contributions: All authors contributed equally to this work.All authors have read and agreed to the published version of the manuscript.
in s.The beam method, ys in flight by counttained a larger value clear Research (CERN).Improved limits on tensor weak forces in helium-6 and lithium-8 decay using the ATLAS Beta Paul Trap have achieved less than uclear and neutron decay.(top) Measurements of the CKM matrix elements V ud and V us es of kaon decays (red, pink, green, and blue, respectively).The yellow ellipse is a global e relationship assuming CKM unitarity.In a world of harmonious experiments described ld all intersect at a single point.Something is amiss.(bottom) Most recent and/or precise time.Measurements using trapped UCNs (red squares and green triangles; red band significantly shorter lifetime than beam-based measurements (blue squares and band).

Figure 1 .
Figure 1.Measured values of the neutron lifetime from beam (blue), material trap (red), and magnetic trap (green) experiments since 1990 from[55].In all neutron trap experiments, which use ultracold neutrons (UCNs), the neutron lifetime was determined by counting the neutrons that persist, whereas in the beam experiments, it was determined by counting the decay protons.The central values of the two sets of experiments disagree by 9.5 s, with a significance in excess of 4 σ-this is the neutron lifetime anomaly.The Particle Data Group (PDG) does not include the beam measurements in their fits[56].

Chemical 7 Hydrostatic 7 3 3 3 Figure 4 .
Figure 4.The potential pathways and resultant states stemming from dark particle decay in neutron stars, dictated by production (G P ) and depletion (G D ) rates.Observational sensitivity (s obs ) and the star's response rates-chemical (G chem ) and hydrodynamic (G hyd )-further influence these pathways.A departure from chemical and hydrostatic quasi-equilibrium necessitates numerical simulations for precise stellar modeling, with the expectation that stars transition from non-equilibrium to a quasi-equilibrium state over time.

Figure 4 .
Figure 4.The potential pathways and resultant states stemming from neutron decays to dark final states in neutron stars, which are dictated by production (Γ P ) and depletion (Γ D ) rates.Observational sensitivity (σ obs ) and the star's response rates-chemical (Γ chem ) and hydrodynamic (Γ hyd ) ones-further influence these pathways.A departure from chemical and/or hydrostatic quasiequilibrium (denoted by ✗) necessitates numerical simulations for precise stellar modeling, with the expectation that stars transition from a nonequilibrium to a quasi-equilibrium state over time.
Total baryon conservation ( ḂD = − ḂV ) Baryon Dark Decay Ξ = {E c } Total baryon loss rate ( Ḃ = −ΓB) 1.Standard Baryon-Conserving Spin-Down Rate of a Pulsar: • Configuration Space: This is defined by Ξ = {E c , ω c }, where E c represents the central energy density, and ω c denotes the central angular velocity of the local inertial frames [200,201].• Constraints: These are governed by the observed spin-down rate ( Ω) while maintaining baryon number conservation ( Ḃ = 0).2.

Figure 5 .
Figure 5. Variations in neutron star properties as functions of central energy density (E c ): mass (M), radius (R), baryon number (B), and moment of inertia (I), all normalized to typical values.Namely, M typical = 1.4 M ⊙ , R typical = 12 km, B typical = 10 57 , and I typical = 70 M ⊙ km 2 .The calculations employ the DS(CMF)-1 EoS [202].Notable features include the quasi-equilibrium baryon loss (indicated by the green arrow) and the onset of hyperon appearance (denoted by cyan vertical dashed lines) at sufficiently large central energy density.
represents the kinetic four-momentum.The vector self-energy (Σ µ ), generated by the vector meson VEVs, is characterized by ⃗ Σ = 0 in the nuclear matter frame.The time component of p * µ is given by E * (p * ) ≡ m * 2 + |⃗ p * | 2 , with m * signifying the baryon's effective mass.The baryon spinor u(p * , λ) satisfying the Dirac equation

Figure 7 .
Figure 7. Variations in baryonic properties under the DS(CMF)-1 EoS depicted as functions of density.(a) Effective masses of various baryons.(b) Vector self-energies of the baryons in the nuclear matter frame.(c) Range of center of mass (CM) frame energies for the baryons.Vertical lines in all panels represent the central number densities (n c ) of the pulsar systems examined in Section 4, while horizontal lines correspond to the vacuum masses of the respective baryons.

Figure 8 .
Figure 8. Energy spectra of baryons (Λ, neutron (n), and proton (p)) in their center of mass (CM) frames within specific neutron stars of (a) PSR J0348+0432, (b) PSR J1614-2230, and (c) PSR J0737-3039A plotted as functions of stellar radius.We have chosen to represent only PSR J0737-3039A here and henceforth, as it is nearly identical in mass to PSR J0737-3039B.The plot highlights how the CM energy in high-density regions exceeds the baryons' vacuum rest mass, thus hinting at the possibility of decay into states that are kinematically forbidden in vacuum.Detailed analysis of these phenomena is provided in the main text.

Figure 10
Figure10.Two-σ exclusion limits on the vacuum branching fraction for B → χγ processes are depicted in panels (a) for neutrons and (b) for Λ baryons.The green region represents limits from Super-Kamiokande[203].We highlight that the KamLAND constraint (red)[204] is pertinent to our study of n → χγ decay, as the prompt photon would not pass the correlation cuts and would remain undetected.Panel (b) also presents limits from neutron star studies for the DS(CMF)-7, DS(CMF)-3, and DS(CMF)-1 EoSs in order of increasing severity, along with constraints from the BES-III experiment (blue)[216], noting that the γ in this case would also be undetected, and the SN1987A limits (gold)[79].

Figure 11
Figure11.Exclusion limits at 2 σ on (a) ε nχ and (b) the vacuum branching fraction for n → χγ from our quasi-static neutron star analysis in the particular χ mass region pertinent to an explanation of the neutron lifetime anomaly.Additional constraints and expected limits (dashed) as per[218] have been included.This figure has been taken from[37].