Main

Recent LIGO-Virgo-KAGRA (LVK) analyses have revealed complex structure in the binary black hole (BBH) population, including distinct features in the primary mass spectrum and nontrivial spin-mass correlations. However, the phenomenological models used to capture these features often lack analytic transparency, making it difficult to isolate robust physical laws from modeling artifacts. To address this, symbolic regression is applied to the posterior inference products of the GWTC-4 catalog, discovering compact, closed-form analytic expressions for four key population relationships: (i) the merger-rate evolution with redshift; (ii) the mass-ratio dependence of the effective-spin distribution; (iii) the redshift evolution of the effective-spin distribution; and (iv) the conditional mass-ratio distributions associated with the 10 solar mass and 35 solar mass primary mass peaks. This framework successfully compresses both rigid and highly flexible models into differentiable phenomenological laws, dynamically recovering a consistent low-redshift merger-rate slope without assuming an a priori power-law form. The exact analytic derivatives provided by symbolic regression show that the mass ratio–effective spin and redshift–effective spin correlations are robustly driven by broadening of the posterior widths rather than shifts in the mean. Furthermore, qualitatively distinct functional forms for the mass-ratio distributions conditioned on the 10 solar mass and 35 solar mass primary mass peaks are identified. These closed-form expressions enable exact analytic gradient diagnostics and compact surrogate summaries, particularly for flexible numerical posteriors that are not otherwise available in low-dimensional analytic form. They also facilitate rapid downstream calculations for rate forecasting, formation channel comparison, and stochastic background estimation.

By monitoring the times of arrival of radio pulses from millisecond pulsars, Pulsar Timing Arrays (PTAs) serve as unique gravitational wave (GW) laboratories in the nanohertz band. To date, the primary astrophysical sources of GWs targeted in this frequency range have been inspiraling supermassive black hole binaries (SMBHBs) on circular and eccentric orbits. In this work, we demonstrate that, thanks to the so-called pulsar term in the timing residual waveform of GW signals, PTAs can probe individual SMBHBs that merged before timing observations began. We refer to the latter as \emph{zombie binaries}. Using SMBHB population models consistent with current PTA constraints, we find that while the probability of detecting such systems in existing PTA datasets remains low, the Square Kilometer Array observatory is expected to achieve sufficient sensitivity to have a few zombie binaries with a signal-to-noise ratio exceeding 3 in its data. Although their confident identification might be challenging, this new class of PTA sources opens a novel window for studying the most massive SMBHBs in our local universe.

Relativistic effects are sensitive to subtle changes in dark energy. These effects grow on very large scales and at high redshifts, which will be the reach of upcoming surveys. We investigate these effects in both the linear and the angular galaxy power spectra in a late-time universe dominated by cold dark matter and k-essence, focusing on three core models (dilaton, tachyon, and DBI scalar fields) and contrasting their predictions with those of the concordance model. By enforcing identical present-day cosmological parameters, we isolate the imprints of k-essence dynamics and perturbations on very large scales. We found that relativistic corrections dominate on very large scales and grow with redshift, but are largely insensitive to k-essence microphysics in Fourier space, leading to strong degeneracies among the models. However, in the angular power spectrum, where line-of-sight integrals are naturally included, relativistic effects are significantly amplified, yielding better sensitivity to clustering k-essence. In particular, the tachyon exhibits clear deviations across multipoles and redshifts, with distinct imprints in the Doppler and the combined (velocity and gravitational) potentials contributions. Furthermore, our results show that neglecting relativistic corrections can lead to systematic misestimation of deviations of k-essence from the cosmological constant. Our results show the relativistic angular galaxy power spectrum as a more consistent and robust probe of ultra-large-scale physics. These findings underscore the need for full relativistic modelling in next-generation surveys that are targeting horizon-scale modes, where the imprint of non-standard dark energy is most pronounced.

Pulsar timing arrays record gravitational waves from supermassive black hole binaries at two spacetime points: an Earth term, measured when the wave passes the Earth, and a pulsar term, measured when the wave passed each pulsar at an earlier epoch. We show that a future $\mu$Hz-band detection of a nearby massive binary by a mission such as $\mu$Ares would turn PTA pulsar terms into targeted probes of binary evolution. In analogy with supernova light echoes, each pulsar term acts as a gravity echo: a dated snapshot of the binary at an earlier stage of its inspiral. Together, the $\mu$Hz Earth-term measurement and the nHz pulsar-term echoes provide a temporal baseline that neither detector could access alone. For a fiducial equal-mass binary with total mass $10^9\,M_\odot$ at 80 Mpc, we find a combined pulsar timing array echo signal-to-noise ratio of 33, with up to 24 pulsars individually resolving the signal among pulsars with 50-year baselines. The angular dependence of the single-pulsar echo sensitivity alone enables independent sky localization of the source to $\sim$10–100 deg$^2$, and the resolved pulsar-term frequencies directly measure the binary inspiral rate hundreds to thousands of years ago. With sufficient pulsar distance precision, a small set of anchor pulsars could additionally phase-connect the array and trace the post-Newtonian evolution coherently over kpc baselines. The source population required for gravity echoes is drawn from the same massive-end census responsible for the observed nanoHertz stochastic background.

While massive black hole binaries (MBHBs) merge at gravitational-wave frequencies above the pulsar timing array (PTA) sensitivity band, we show that they leave orphaned low-frequency contributions in the PTA pulsar term. Due to the light-propagation time between each pulsar in the array and Earth, the pulsar term acts as a time-delayed probe of a chirping merger with a specific frequency response determined by the direction of origin and intrinsic properties of the MBHB. We provide a detailed consideration of how such a multiband signal would manifest in a full PTA, demonstrate an approach to stack these orphaned pulsar terms across the array, and discuss prospects for an archival, multiband search in conjunction with MBHB mergers observed in astrometric data or spaceborne interferometers like LISA.

We compute the abundance of cosmological gravitational waves produced during the evolution of an ultralight vector (spin-1) dark matter field. A homogeneous background vector field breaks spatial isotropy, requiring a Bianchi I geometry and inducing a mixing between the scalar, vector, and tensor perturbation sectors. We derive the perturbation equations in this background and show that, as a consequence of this mixing, scalar perturbations act as a source of tensor modes, generating a stochastic GW background. The production and cosmological evolution of these gravitational waves are implemented in a modified version of CLASS, from which we obtain the present-day spectrum.

In core-collapse supernovae (CCSNe), deleptonization normally favors $\nu_e$ over $\bar{\nu}_e$ emission. However, lepton-number emission self-sustained asymmetry (LESA) can make the energy-integrated emitted lepton-number flux negative along some directions. We derive a simple diagnostic for this transition and test it in 33 independent 3D CCSN simulations: 25 Princeton/Fornax models ($8.1$–$100\,M_\odot$) and 8 Garching models, including non-, slow-, and fast-rotating $15\,M_\odot$ cases. Of 23 non-black-hole-forming Princeton models, 22 cross the threshold, with median onset $t_c=225\,\mathrm{ms}$, IQR $162$–$264\,\mathrm{ms}$, and cross-model scatter $\mathrm{CV}=18.6\%$. Full-sky flux-sign searches show that the threshold identifies the anti-LESA-pole transition, distinguishing the global LESA-driven crossing from early localized turbulent crossings. The fast-rotating Garching $15\,M_\odot$ model, where rapid rotation suppresses the LESA dipole, is correctly classified as a non-crosser without using any rotation parameter. Both black-hole-forming Princeton models cross near $250\,\mathrm{ms}$ post-bounce and remain above threshold for $1807$ and $2463\,\mathrm{ms}$ before collapse. Thus, in the next nearby CCSN, the emitted $\bar{\nu}_e$ energy flux may exceed the $\nu_e$ flux along some lines of sight. Such directions may also correlate with sustained fast flavor instability, although testing this requires local phase-space distributions or dedicated linear stability analysis. The relevant quantity here is the energy-integrated emitted flux field, i.e. a luminosity difference per steradian, not a neutrino number flux.

We present SDO2IRIS$^2$: a visual transformer model that translates a combination of images of the chromosphere and transition region (TR), observed by AIA, and a line-of-sight magnetogram, provided by HMI, into temperature, line-of-sight velocity (v$_{los}$), velocity of the turbulent motions (v$_{turb}$), and electron density (n$_{e}$) in the chromosphere. Using the thermodynamic variables obtained from the inversion of the chromospheric lines Mg II h&k, observed by IRIS, as the target of the model, and the intensity images in the chromosphere and TR, and the photospheric magnetogram as the input, the predicted T and n$_{e}$ show a strong correlation ($\approx 0.80$) for $\approx$80% of the test inverted data, a moderate-to-strong correlation ($\approx0.63$) for 70% of the v$_{turb}$ of the target test inverted data, while for the $v_{los}$, the correlation is weak. Therefore, the predicted values by SDO2IRIS$^2$ may be used as an estimation of the thermodynamics in the chromosphere, either as a stand-alone result or as complementary information to other chromospheric data observed simultaneously. The execution time employed by SDO2IRIS$^2$ to obtain the thermodynamic values in the chromosphere is of the order of a few minutes, being $\le10$ minutes when using a CPU, and $\le5$ minutes when using a GPU. SDO2IRIS$^2$ opens a new avenue for the use of SDO data thanks to the inversions provided by IRIS observables.

We investigate dense-matter equations of state (EOSs) within a Bayesian framework, with particular emphasis on whether recent small-radius compact-star candidates can be accommodated in a twin-star scenario. For the hadronic sector, we adopt a meta-modeling EOS constrained by the NICER mass–radius measurements of PSR J0030$+$0451, PSR J0437$-$4715, PSR J0614$-$3329, and the massive pulsar PSR J0740$+$6620. The hadronic inference indicates that PSR J0614$-$3329 favors a somewhat softer EOS than the other two \(\sim1.4\,M_\odot\) pulsars, while the \(\sim2\,M_\odot\) constraint prevents the EOS from becoming too soft. We then introduce a strong first-order phase transition through a constant-speed-of-sound quark-matter segment. Using HESS J1731$-$347 and XTE J1814$-$338 to constrain the phase-transition parameters, we find a preferred transition density of \(n_\mathrm{t}\sim2.7\text{--}2.8\,n_0\), a sizable energy-density jump of \(600\text{--}700\) MeV, and a relatively large post-transition sound speed of \(c_s^2/c^2\sim0.85\). Such a phase transition generates a disconnected hybrid branch with radii of about \(6\text{--}7\) km at masses around \(1.2\text{--}1.4\,M_\odot\), and strongly suppresses the dimensionless tidal deformability relative to the purely hadronic branch. This pronounced change in tidal deformability is a characteristic signature of the twin-star mechanism and may provide an important observational tool for identifying phase transitions in neutron-star matter in future multimessenger measurements. These results show that small-radius compact stars can provide direct constraints on both the strength of a first-order phase transition and the stiffness of the post-transition phase in dense matter.

Black-hole accretion systems exhibit a characteristic coexistence of activities: broad-band X-ray variability, hot coronae, wide-angle winds, and both steady and discrete jets. This coexistence suggests a persistently time-dependent magnetic background in which noisy fluctuations and explosive release are both essential. In this paper, we connect them all to intermittent magnetic reconnection and propose a Synchronized Spin Model (SSM) in which multiple local dynamos in a rotating accretion flow are represented as interacting macro-spins. Their synchronization, partial synchronization, excursion, and reversal define a compact set of collective variables that organize both timing statistics and large-scale morphology. In this picture, multiscale magnetic reconnection sustains coronal heating, flares, intermittent outflows, and discrete jet activity, while the same synchronization dynamics produce amplitude modulation and demodulation, providing a route to $1/f$-like variability, rms–flux/Taylor-like scaling, and approximately log-normal statistics of the demodulated envelope. We further argue that, although the continuous flux distribution in black-hole systems is more naturally discussed in multiplicative or log-normal terms, broader event-catalog statistics remain useful for describing suitably defined burst hierarchies, particularly by analogy with solar and stellar flare systems. The hard/soft cycle of X-ray binaries is then interpreted as motion through magnetic state space.

Be stars are rapid rotators generally produced by binary interactions. The single Be stars in the observations pose challenges to the Be star formation theory. In this paper, we propose a new pathway for the formation of single Be stars, in which the Be star is taken as the ejected companion star from a Type Ia supernova (SN Ia) explosion. Recent numerical simulations suggest that explosive oxygen burning, initiated via the convective Urca process in certain helium (He) stars near the Chandrasekhar mass limit, can set off a SN Ia. Based on this proposition, we further demonstrate that about $0.4\%$ of He star + main-sequence (MS) star binaries may evolve into single Be stars, where the MS star is spun up due to the mass accretion from the He star, and then the He star explodes as a SN Ia. We employ detailed binary evolutionary simulations and find the parameter space that would produce single Be stars via the SN Ia channel. Around $22\%$ of Be stars from the SN Ia progenitor channel exhibit peculiar tangential velocities exceeding $24\ \rm km/s$, classifying them as runaway stars. This suggests that the SN Ia channel plays a meaningful role in forming single Be stars, particularly within the runaway star population.

We investigate electron–positron pair creation through the Schwinger process in newborn magnetars with millisecond spin periods and surface dipole fields close to or above the QED critical field, $B_{\rm Q} = 4.414\times10^{13}\,\mathrm{G}$. In the unscreened field scenario, we derive the analytical global pair creation flux and recast it into a compact form with accurate analytic approximations. For a fiducial model with $B_{\rm p} = 10^{14}\,\mathrm{G}$ and $P_0 = 1\,\mathrm{ms}$, the Schwinger channel exceeds the classical Goldreich–Julian particle supply by many orders of magnitude and becomes the dominant source of charges at the earliest stage of the magnetar. The associated discharge removes about $90\%$ of the initial rotational energy within 30 ms, suppresses the gravitational-wave loss channel, and implies that the observable millisecond phase is extremely short in this unscreened scenario. The rapid energy release over such a short timescale may also provide a viable power source for astrophysical transients. Extending the same fiducial model to $10^4\,\mathrm{yr}$ gives spin periods of order seconds, linking newborn millisecond magnetars to the mature magnetar population.

White-light flares (WLFs) are usually tend to be those very large flares. Nevertheless, several small and compact WLFs have been reported and thought to be produced by low-height magnetic reconnection. However, whether low-height magnetic reconnection can efficiently produce WLFs remains unclear. For the first time, we conduct a statistical study of the WL emission in compact flares to address this question. Using over a decade observations from the \textit{Solar Dynamics Observatory} (SDO), we identify 28 compact flares, including 19 C-class and 9 B-class flares. We find these compact flares can be classified into three types based on the magnetic configuration of the flare, corresponding to the U-shape loop (type I), the flux emergence near sunspot (type II), and the fan-spine like structure (type III). For each type, the flares numbers are 9 (7 C-calss and 2 B-class), 9 (3 C-calss and 6 B-calss) and 10 (9 C-calss and 1 B-calss), respectively. We find the occurrence rate of WLFs in compact flares is $\sim60.7\%$ (17/28), and for the C-class the rate can be up to $\sim89.5\%$ (17/19). No WLF was found in B-class compact flares. The occurrence rates for three types are $\sim77.8\%$ (7/9), $\sim11.1\%$ (1/9) and 90\% (9/10), respectively. And for the C-class flares, the occurrence rates for three types are 100\% (7/7), $\sim33.3\%$ (1/3) and 100\% (9/9), respectively. Our results suggest type-I and type-III compact flares are more likely to produce WL emissions.

We present optical B and V band photometry of GP Com - an ultracompact binary consisting of an accreting white dwarf and helium secondary component. Our data set contains 7.7 hours observations in V band with the 2.0m telescope and 2.9 hours simultaneous observations in B and V bands with the 1.5m telescope of the Rozhen National Astronomical Obsevatory, Bulgaria. The observations cover of about 13 orbital periods. We find an orbital modulation with amplitude 0.04-0.05 mag in B and V bands. Adopting that it is due to a bright spot, we estimate its temperature 19700 \pm 3000 K. We estimate mass accretion rate onto the white dwarf of about 2.10^{-12} M_sun/yr, consistent with the predicted rate for a cool donor. The data are available on Zenodo: this http URL .

Context. For stellar evolution models we rely on mass-loss rate prescriptions that show features that lack empirical confirmation, such as the bi-stability jump. This jump is an increase in mass loss in the decreasing temperature regime Teff 28-21 kK. Although papers compared empirical results to prescriptions,a large observational sample of B stars for which the wind has been scrutinised over different metallicities is still lacking. Aims. By modelling of both ultraviolet (ULLYSES) and optical (XShootU) spectra, we determined the stellar and wind parameters, of 24 SMC B stars ranging in Teff from 13 to 29 kK. By combining this sample with LMC studies, we compared the wind behaviour of B stars in two different metallicity regimes. We compared our results to common mass-loss rate prescriptions to test features present in these models and their metallicity dependence. Methods. We have used the model atmosphere code fastwind and the genetic algorithm code Kiwi-GA to fit the UV and optical spectra simultaneously. This allows us to determine wind properties including clumping behaviour. Results. The metallicity trends present in the mass-loss prescriptions (Z^(0.41-1.4)) explored here overestimate the empirical metallicity dependence in the B-star regime, which appears very weak. We do not find an increase in mass-loss rate at approximately spectral type B1. We show that on average 40% of the wind mass is located in the wind medium between the clumps. We compiled a sample of 80+ O and B stars in the SMC and LMC. From a comparison we find a clear difference in O- and B-type metallicity dependence. Conclusions. The lack of a bi-stability jump in the B-star regime and a weak metallicity dependence for the same stars offers new empirical constraints to models of line-driven winds. As differences between these models are large (reaching orders of magnitude) such constraints are much needed.

Non-conservative mass transfer plays a central role in close-binary evolution, yet its effects on mass-transfer stability are uncertain. One widely adopted prescription, isotropic re-emission, is often assumed to promote stability compared to conservative mass transfer. We investigate the impact of isotropic re-emission on the stability of mass transfer in binaries with radiative envelopes that undergo delayed dynamical instability (DDI). We assess whether simplified criteria used in binary population synthesis codes accurately capture stability boundaries under varying mass-transfer efficiencies. We perform detailed stellar evolution calculations for a set of representative binaries undergoing DDI. Varying the mass-transfer efficiency beta, we track the onset of instability and quantify the corresponding critical mass ratio. We compare our results with predictions from the commonly used zeta-method, which is based on mass-radius exponents. We find that a lower mass-transfer efficiency destabilizes mass transfer in DDI systems, whereas the zeta-method predicts that isotropic re-emission should stabilize it. The discrepancy arises because the zeta-method fails to capture the full evolution of the orbit and mass ratio during pre-instability mass transfer. In some cases, the critical mass ratio is underestimated by nearly a factor of two. Our findings show that isotropic re-emission can reduce, rather than enhance, DDI stability, underscoring the limitations of using fixed critical mass ratios and zeta-based criteria. This highlights the need for calibrated prescriptions that capture the time-dependent evolution of mass ratio and orbital separation, with direct implications for modelling X-ray binaries, symbiotic stars, and double white dwarfs, including their transient rates and delay-time distributions.

We investigate the viability of Tsallis holographic dark energy (THDE) models, focusing on the role of the infrared (IR) cutoff in the growth of cosmic structures. Considering two commonly used choices of the cutoff, the particle horizon and the future event horizon, we analyze the evolution of linear matter perturbations and compute the growth factor, growth rate, and the observable $f\sigma_8(z)$. These predictions are compared with observational data from redshift-space distortion measurements. We find that the growth history is highly sensitive to the choice of IR cutoff. Models based on the future event horizon are consistent with observational data and can provide a fit comparable to, or slightly better than, the $\Lambda$CDM model for suitable values of the Tsallis parameter $\delta$. In contrast, models constructed using the particle horizon generally fail to reproduce the observed growth of structure. These results demonstrate that the viability of THDE models depends crucially on the choice of IR cutoff and highlight the importance of structure formation as a stringent test of generalized holographic dark energy scenarios.

We present the first high-resolution X-ray spectral analysis of Cygnus X-1 using XRISM. The observation was carried out from April 7 to 10, 2024, covering the orbital phase range 0.65–0.17 during its low/hard state. Taking advantage of the exceptional energy resolution of the Resolve instrument, we examined highly ionized iron absorption lines and characterized the ionization states, column densities, and line-of-sight velocities of the absorbing plasma. Spectral analysis revealed an ionization parameter of approximately 3, column densities of a few times 10^21 cm^-2, and a blueshifted velocity of approximately 100 km s^-1. The observation was divided into two phases: before and after orbital phase phi_orb = 0.9, corresponding to non-dipping and dipping intervals. While only weak absorption features were present before phi_orb = 0.9, strong absorption by He-like and H-like Fe appeared during the dipping phase. We measured equivalent widths of 2.3 eV, 0.4 eV, and 1.2 eV for He-like Fe K-alpha, and H-like Ly-alpha1 and Ly-alpha2, respectively, demonstrating the capability of XRISM Resolve to securely detect narrow absorption features of only a few eV. These measurements trace the motion of the absorbing material and offer insight into the kinematics and spatial distribution of the wind in the vicinity of the black hole. These findings enhance our understanding of wind-fed accretion in Cygnus X-1 and highlight the importance of continued high-resolution X-ray observations to further constrain the physical properties of winds and accretion flows in high-mass X-ray binaries.

We demonstrate that discrete self-similarity (DSS), originally discovered in the collapse of a massless scalar field in an asymptotically flat system, survives in primordial black hole (PBH) formation within an expanding cosmological background. Using fully relativistic numerical simulations of massless scalar-field collapse in an Friedmann-Lemaître-Robertson-Walker universe, we resolve the critical regime down to $|p-p_c|\sim 10^{-8}$, where $p$ and $p_c$ respectively are a parameter of the family of initial data and its threshold value, and find clear log-periodic oscillations in the PBH mass scaling relation. The detailed structure of these oscillations differs from that previously reported in the asymptotically flat case, exhibiting a more pronounced asymmetry between peaks and troughs. Analyzing two distinct families of initial data (Gaussian and piecewise rational curvature profiles), we find critical exponents and DSS periods that differ slightly but are broadly consistent within uncertainties. The presence of DSS implies characteristic log-periodic modulations in the PBH mass spectrum, with potential consequences for PBH abundances and the spectrum of induced gravitational waves.

We present the first high-resolution spectroscopy of the black hole high-mass X-ray binary Cygnus X-1 with XRISM, including orbital-phase-resolved analyses and tentative evidence of short-term variability in the Fe-K band on second timescales. Using data from the Performance Verification phase in April 2024, we analyzed spectral variability across orbital phases with the Resolve microcalorimeter and the Xtend CCD imager. The unprecedented resolution of Resolve reveals variability in highly ionized Fe absorption lines. The absorption features show orbital-phase-dependent variability in column density, ionization state, and blueshifted velocity, suggesting structural variations in the focused stellar wind along the line of sight. We also find indications of subtle broadening of the neutral Fe emission profile. In addition, intensity-sorted spectroscopy during dip phases suggests possible variability on timescales of a few seconds in the absorption features, consistent with cooler, denser and lower-ionized gas clumps. Although the statistical significance is limited, these results hint that the stellar wind and the X-rays from the accretion disk around the black hole may interact on timescales as short as a few seconds. These XRISM results constrain wind-fed accretion in Cyg X-1 and highlight Resolve's capability to probe plasma environments in high-mass X-ray binaries.

The mass-loss process of red supergiant (RSG) and asymptotic giant branch (AGB) stars and its relation to variability is poorly constrained. We study two evolved stars, the Mira-type AGB star R Car and the extreme RSG VX Sgr. Our sample comprises 54 VLTI-GRAVITY snapshots taken over 7 years, being the largest VLTI time-series dataset to date. We determine the angular diameter as a function of time. The radii of the photosphere ($R_{\star}$) and atomic atmospheric layers are variable and relate to the light curve with phase shifts, showing a maximum radius near visual brightness minima. The more extended CO layers show longer, irregular periods and maximum extensions of $\sim 1.3-1.7 \: R_{\star}$ for R Car, and of $\sim 1.5-2.2 \: R_{\star}$ for VX Sgr. Comparison with CO5BOLD simulations revealed a similar behavior. Furthermore, during 2020-2021, VX Sgr exhibited an extreme mass-loss event similar to that of Betelgeuse, preceded by two strong shocks and culminating with the extreme expansion of H$_2$O and CO layers, both up to $\sim 2.2 \: R_{\star}$. During this event, we detected Brackett $\gamma$ and Balmer emission lines, both of which are signatures of a shock propagating through the atmosphere. The Mira R Car showed a photospheric radius $R_{\star} = 280 \pm 25 \: \rm R_\odot$, with a fundamental mode (FM) pulsation amplitude $\sim13 \%$ of $R_{\star} $. During its active cycle, the RSG VX Sgr showed $R_{\star} = 1556 \pm 110 \: \rm R_\odot$ with FM amplitude $ \sim13 \%$ of $R_{\star} $, the same as R Car. During its quiescent cycle, it showed $R_{\star}= 1456 \pm 108 \: \rm R_\odot$ and low-amplitude pulsations near the first overtone, only $\sim4 \%$ of $R_{\star} $. This supports a steady mass loss for Miras related to stable, large-amplitude FM pulsation, whereas the mass-loss process for RSGs may be dominated by extreme events connected to changes in the pulsation mode.

We present multi-wavelength study of the $\gamma$/X-ray transient EP250416a (also designated GRB 250416C), triggered by the Einstein Probe (EP) Wide-field X-ray Telescope and also by SVOM and Konus-Wind. Observations spanning the gamma-ray, X-ray, and optical bands facilitated detailed analysis of the burst's prompt emission, afterglow evolution, and physical origin. EP250416a exhibits a burst duration of 30 s in X-ray and 17.7 s in gamma-rays, with joint spectral fitting of 0.5-5000 keV data gives $E\rm_{peak}=342_{-232}^{+90}$ keV. Optical spectroscopy of the afterglow, acquired with the Gemini Multi-Object Spectrograph (GMOS) on Gemini South, yielded a redshift of $z=0.963$. Accounting for the measured redshift, the isotropic energies are $E\rm_{X,iso}=2.7_{-0.5}^{+0.9}\times10^{50}$ erg and $E\rm_{\gamma,iso}=7.34_{-2.1}^{+5.1}\times10^{51}$ erg, aligning with the Amati relation for long GRBs. The fluence ratio $\rm S(25-50~keV)/S(50-100~keV)=0.78_{-0.15}^{+0.1}$ classifies EP250416a as an X-ray rich (XRR) GRB. The X-ray afterglow shows an initial shallow decay ($\alpha \approx -0.5$) transitioning to a canonical decay phase ($\alpha \approx -1$), with a very late jet break at $t\sim 1.5\times 10^6$ s, corresponding to a jet half-opening angle of $\theta _j=10.6_{-1.8}^{+1.9}$ degrees. EP250416a is optically dark, as it shows only a faint $r$-band detection ($r=24.16$ mag) from Gemini South-GMOS and a low optical-to-X-ray spectral index $\beta_{\rm OX} = 0.3$. This may be attributed to significant host-galaxy extinction, with a required $A_V^{\text{host}}=5.5\ \text{mag}$ derived from the extinction curve model.

The distribution of galaxies, halo abundance, and peculiar velocities are influenced by non-linear gravitational interactions, making the study of non-linear evolution crucial for accurate cosmological predictions. We explore these aspects using N-body simulations. Theoretical models of the halo mass function (HMF) can be formulated without referencing a cosmological model or input power spectrum. HMF obtained from N-body simulations show systematic deviations of 5-20\% from theoretical predictions. The physical origin of deviations may result from cosmology, the power spectrum, or both. We examine HMF deviations from universality for scale-free power spectra with an Einstein-de Sitter cosmology. We demonstrate that the mass function exhibits an explicit dependence on the slope of the input power spectrum. We find that an effective index of the $\Lambda$CDM model can correspond to the HMF from scale-free cosmologies as a first approximation. Furthermore, structure formation has led to deviations from homogeneity and isotropy on scales up to at least $100$ Mpc/h, expected to affect measurements of $H_0$. We revisit this issue of the concordance model. We find a correlation between errors in $H_0$ estimates and the density around the observer. Further, our mock observations reveal that deviations of up to 5\% can occur in Milky Way-sized halos. While this finding alone does not fully resolve the Hubble tension, it may account for part of it. It is essential to understand the limitations of N-body simulations to avoid misinterpreting data. We show that the missing power at small scales introduces errors in the root-mean-square fluctuations and in the simulated mass function. Our analytical calculation indicates that mode coupling between small and large scales depends on resolving collapsed halos. Therefore, accurate mode coupling estimates require sufficient halos in the simulation.

We apply the slip-back mapping method of Titov et al. 2009 and Lionello et al. 2020 to a thermodynamic MHD simulation to track topological changes in the magnetic field at a range of temporal cadences. The method constitutes the logical successor to a simple open-field map for a steady-state model, as it tracks changes in the open and closed fields for a time-dependent model by tracking individual magnetic elements as they advect across the map, rather than simply tracing field line connectivity from each cell. Through careful categorization of the slip-back mapping values and analysis of the flux changes, we not only effectively track the open flux but can recover the flux processed through interchange reconnection as well. The field lines involved in these processes are shown to follow lines of high squashing factor, as proposed by interchange reconnection-driven slow solar wind theory. The time-dependent model, which is scaled to solar minimum-like activity, projects that a median value of 3.5% of the total open flux in any given 24-hour interval has been processed through interchange reconnection. This corresponds to a relatively high proportion of the total open flux changes over time in the heliosphere. Our results show that not only is this method a useful tool for accurately tracking topological change in time-dependent simulations, but that its inherent complexity can be visually reduced into an intuitive 2D plot that simply and effectively communicates temporal changes.

Using the ssrAp star 33 Lib (HD137949) as an example, we show that magnetic braking caused by the interaction of a strong magnetic field with a stellar wind can play a key role in slowing the rotation of ssrAp stars. Numerical modeling of stellar rotation spin-down in the MESA package, taking into account the evolution of magnetic fields and the stellar wind, shows that stars with rotation periods of up to 80 years and longer can form. Moreover, braking by a magnetized wind makes it possible to estimate the mass-loss rate for stars of moderate mass ($1.25 M_\odot < M < 2 M_\odot$), which is difficult to do by other methods. We introduce the dimensionless parameter $\Xi$, which reflects the spin-down time and the stellar lifetime. Thus, when $\Xi \gg 1$, braking is important, whereas when $\Xi \ll 1$, it is negligible.

The 21-cm signal, one of the most promising probes of the high-redshift Universe, has traditionally been modelled without accounting for the effects of active galactic nuclei (AGN) in the pre-JWST era, primarily due to the lack of observational evidence for AGNs at z > 6. However, following the discovery of several AGNs at redshifts as high as z   10 by JWST, it has become imperative to incorporate the impact of these early AGNs when predicting the 21-cm signal. Supposing that these AGNs are seeded by primordial black holes (PBHs), we study their impact with a semi-numerical model setup. Specifically, we extended the explicitly photon-conserving reionization framework, SCRIPT, including essential cosmic dawn physics and PBH contributions. This enables us to compute both the global signal and the power spectrum of the 21-cm line over the redshift range z   30 - 5 within a self-consistent framework. Building on this setup, we then investigate the impact of different PBH mass functions (obeying current observational constraints) on the resulting signal. The X-ray heating from PBHs can substantially make the depth of the global 21-cm signal shallower and suppress the expected power amplitude during cosmic dawn. We also find that the choice of mass function plays a crucial role in shaping the 21-cm signal, and can, in fact, lead to significantly different predictions.

We introduce a family of phenomenological cosmological models featuring an interacting dark sector modulated by a sparseness scale parameter, in order to describe the late-time accelerated expansion of the universe. The sparseness scale, inspired by well-established saturation mechanisms in ecology and biology, is introduced in the interaction as a half-saturation constant that bounds the energy exchange between dark matter and dark energy, controls the dynamical behaviour of the physical variables and can prevent the phantom crossing. We consider three nonlinear interacting models, where two of them recover the linear interacting scenarios when the sparsity parameter vanishes. We examine the phase-space of the cosmological field equations by using the Hubble normalization approach. We determine the stationary points and their stability properties in order to reconstruct the asymptotics behaviour of the field equations. Such an analysis allows us to demonstrate the effects of the sparseness scale on the background dynamics. We test the interacting models with observational data. Specifically, we employ Supernovae catalogues, cosmic chronometers, Baryon Acoustic Oscillation measurements from DESI DR2, and redshift-space distortion measurements of the growth of large-scale structure through the $f$ and $f\sigma_8$ observables. The Bayesian analysis suggests that, for two of the three models, a vanishing sparsity parameter is disfavoured at more than the 95\% confidence interval, providing observational support for a nonzero sparseness scale in the dark sector interaction.

The Seyfert 1 galaxy NGC 3783 hosts a multiphase warm absorber (WA) that has been extensively studied in the X-ray band. High-resolution spectra from 2000-2001 revealed a complex outflow with multiple ionization and velocity components. Two decades later, new XMM-Newton and XRISM observations allow us to investigate the long-term evolution of these outflows. We perform joint spectral modeling of the XMM-Newton/RGS and XRISM/Resolve time-averaged spectra using the pion photoionization code within SPEX. We derive the ionization parameter, column density, turbulent velocity, and outflow velocity for each absorption component, and investigate their thermal stability and Absorption Measure Distribution (AMD) to characterize the physical and dynamical properties of the WA in NGC 3783 in 2024. We compare these results with the 2000-2001 epoch to assess long-term variability, stability, and possible changes in the absorber population. We identify eight WA components spanning log $\xi =$ 1.08-3.38 and outflow velocities of 480-1230 km s$^{-1}$. The ranges of column densities and turbulent velocities remain broadly consistent with the WAs from 2000-2001, but the earlier data contained more low-ionization, high-velocity components. The total column density in 2024 is 1.5 times larger than in 2000-2001, requiring replenishment by fresh material. The dominant Unresolved Transition Array (UTA) absorber (Comp. B3) has increased its column density by a factor of three while maintaining a similar ionization parameter. The WAs in NGC 3783 have undergone significant structural and dynamical evolution over the past 24 years.

High radio frequencies observations with the Italian network of large single-dish radio telescopes resulted in  450 solar images between 2018 and 2023 in K-band frequency range (18-26 GHz). Solar radio mapping at these frequencies allows the probing of the Active Regions (ARs) chromospheric magnetic field close to the Transition Region, where strong flares and coronal mass ejection events occur. Enhanced magnetic fields up to 1500-2000 G determine anomalous spectra in the ARs brightness compared to pure free-free emission, due to the addition of a steeper gyro-resonance component also associated with circular polarisation up to  40%. When a significant AR spectral flattening is detected, the probability of a strong flare occurrence within  30 hours is high ( 89% in terms of statistical precision). Despite an approximate weekly cadence of our observations, only  12% of strong flares are missed/unpredicted within this time interval. Through a correlation analysis, we assess the trade-off on the sensitivity and the robustness of this physics-based flare forecast method.

We model the multi-band lightcurves of 80 SNe Ic-BL, including 11 associated with lGRBs, using a magnetar engine model with $^{56}$Ni decay. We find that the data are all consistent with a magnetar central engine, and such a model yields high-quality fits across the sample. The medians with $1\sigma$ regions of the key parameters are $P_{\rm{i}}\sim2.04^{+1.84}_{-0.96}\,{\rm{ms}}$, $B_{\rm{p}}\sim3.96^{+3.28}_{-1.40}\times10^{15}\,{\rm{G}}$, $M_{\rm{ej}}\sim2.30^{+1.48}_{-1.02}\,M_\odot$, and $M_{\rm{Ni}}\sim0.18^{+0.14}_{-0.09}\,M_\odot$, with strong and statistically significant correlations observed for both $M_{\rm{ej}}-P_{\rm{i}}$ (anti-correlation) and $M_{\rm{Ni}}-M_{\rm{ej}}$ (correlation). Comparing the SN Ic-BL samples with and without lGRB association using fitting parameters, we find no significant difference between them, although the GRB-associated sample is slightly brighter, possibly due to an observational bias. Relative to ordinary SNe Ic, SNe Ic-BL have similar $^{56}$Ni and ejecta masses, suggesting comparable pre-SN progenitor properties, with differences possibly arising from the presence of a magnetar engine. In comparison with other possible magnetar-powered SESNe, including SLSNe Ic and FBOTs, we confirm a strong universal $M_{\rm{ej}}-P_{\rm{i}}$ correlation, indicating a common origin. SNe Ic-BL and SLSNe Ic have similar ejecta mass distributions, typically $M_{\rm ej}\gtrsim0.5\,M_\odot$, while FBOTs mostly lie below this value. Differences between SNe Ic-BL and SLSNe Ic may arise from magnetar properties, with SN Ic-BL magnetars rotating faster and having stronger fields. Moreover, the $P_{\rm{i}}-B_{\rm{p}}$ distribution of lGRB magnetars largely overlaps with that of SN Ic-BL magnetars. In connection with binary simulation results, we propose a unified physical classification and progenitor framework for magnetar-powered and ordinary SESNe.

Constraining primordial non-Gaussianity via its scale-dependent imprint on galaxy clustering requires knowledge of the bias parameter $b_{\phi}$, which is exactly degenerate with $f^{\rm{loc}}_{\rm{NL}}$ at leading order. To break this degeneracy, current analyses adopt the relation $\left(b_{\phi} = 2\delta_c\left(b_1 - 1\right)\right)$ based on the assumption of a universal mass function. This relation is known to break down for physically motivated galaxy selections, introducing systematic errors in the inferred $f^{\rm{loc}}_{\rm{NL}}$ that scale directly with the assumed $b_{\phi}$ prior. We present a framework to construct physically motivated, observation-conditioned priors on $b_{\phi}$ by marginalizing over galaxy formation uncertainties. We use the CAMELS-SAM simulation suite, augmented by separate Universe simulations, to measure galaxy formation observables, like the stellar mass function (SMF) and the stellar-to-halo mass relationship (SHMR), and $b_{\phi}$ across a range of galaxy formation parameters. From these measurements, we construct a distribution of $b_{\phi}$ conditioned on observations, and we select our galaxy sample to resemble the DESI Emission Line Galaxy (ELG) sample. Conditioning on the SMF or SHMR decreases $\sigma_{b_{\phi}}$ from $0.69$ to $0.08$ and $0.02$ respectively – reductions of $88\%$ and $97\%$ – with consistent results when conditioning on the observed data directly. Despite substantial shifts in the galaxy formation posteriors driven by known SC-SAM discrepancies at high halo masses, the resulting $b_{\phi}$ distributions remain mutually consistent across all observables. The SMF and SHMR are found to carry sufficient constraining power to reduce the galaxy formation uncertainty in $b_{\phi}$ relevant for $f^{\rm{loc}}_{\rm{NL}}$ inference with next-generation spectroscopic surveys

We have searched for the 3.55 keV line from sterile neutrino decay using 3.1 megaseconds of RXTE cluster data. A 2.5$\sigma$ excess of emission over a thermal model is found over the energy span of the 3.55 keV line in the combined spectra of the eight clusters that individually have an excess. The residuals are added to increase the signal to noise ratio of the excess, which is then modeled with a Gaussian to simulate the instrumental spectral response. We find a significant correlation (r = 0.76) for a line centered at 3.6 keV with a model flux of 3.07 x 10$^{-5}$ ph cm$^{-2}$ s$^{-1}$. Mixing angle for detected clusters ranges from 0.35 to 6.2 x 10$^{-10}$. The decay rate inferred from the line flux is strongly correlated (r = 0.87) with cluster temperature, which is due to hotter, more massive clusters having a larger amount of dark matter. Approximately half of the decay line total flux comes from the Coma cluster. We fit the Coma cluster spectrum with two different three-component models. The first includes a Gaussian fixed at 3.55 keV to model soft emission. The second three-component model uses a second thermal component to model soft emission. The model fit parameters indicate that the second thermal component is modeling high-energy residuals rather than low ones, where the Gaussian is. Though our line fluxes exceed most reported detections and upper limits, they do not overproduce the dark matter. We conclude that some fraction of the marginally detected excess could be attributed to the decay line since low-temperature thermal emission and systematics fail to model it completely.

The James Webb Space Telescope and the future Habitable Worlds Observatory aim to discover exoplanet atmospheric spectra that detect life. Currently, most existing spectral "retrieval" algorithms focus on inferring the abundances of biogenic gases from these spectra. However, abundances are hard to interpret as signatures of life because they are modified by photochemistry, climate, and atmospheric escape. To address this problem, we develop a method for inferring the fluxes of gases at a planetary surface by inverting a coupled photochemical-climate model. As a proof-of-concept, we apply the approach to a synthetic 10-transit JWST NIRSpec Prism spectrum of TRAPPIST-1 e assuming it hosts a biosphere similar to the Archean Earth's. The retrieval confidently detects CO$_2$ and CH$_4$ and can constrain the flux of CH$_4$ into the atmosphere to within approximately 1.5 orders of magnitude (68$\%$ credible interval) provided that TRAPPIST-1's near-UV spectrum is accurately known. We demonstrate how inferred surface gas fluxes naturally fold into a probabilistic assessment of life, finding that   80$\%$ of the surface gas flux posterior is consistent with a CH$_4$-producing metabolism for our nominal test case. As with any inverse problem, these results are conditional on a number of assumptions in our forward model. Overall, we argue that increasing the robustness of life detection on exoplanets requires moving beyond atmospheric abundances toward inference of the surface fluxes that sustain them.

We use GRAVITY near-infrared (NIR) flare astrometry to test whether Sagittarius A* could be a solitonic boson star. We consider five spherically symmetric solitonic boson-star models with different effective radii, together with the Schwarzschild black hole. Treating the flares as hot spots on circular equatorial orbits, we analyze their centroid motions and images in these spacetimes and use them for parameter fitting. We perform the fitting using both $\chi^2$ analysis and Markov Chain Monte Carlo (MCMC) methods, which yield consistent results: the inferred masses of boson-star models are systematically larger than the established value of $4.3\times10^6M_\odot$. Notably, more diffusive boson stars exhibit imaging properties closer to those of a black hole, leading to mass estimates that are correspondingly closer to the established value. Overall, our results place stringent constraints on solitonic boson star interpretations of Sagittarius A*, although do not completely rule them out.

Massive stars with their strong ionizing radiation and strong stellar winds are the key feedback agents of the universe. Stellar winds of massive stars are often measured by fitting resonance lines in the UV using non-LTE stellar atmosphere models. So far, the line formation regions of these lines have not been measured empirically, preventing a comparison to the model's structures. We aim to conduct the first measurement of the resonance line formation regions in an early-type eclipsing binary in the SMC, namely AzV 75. We employ TESS and ASAS-SN photometry in combination with radial velocity measurements from multi-epoch HST UV spectra to derive the ephemeris. We examine the intensity changes in the C IV and N V resonance lines in the UV and combine them with a light-curve analysis to estimate the region in the wind where these lines are formed. AzV 75 has an orbital period P=165.66d, eccentricity e=0.42, mass ratio q=0.72, and inclination i=85.77°. With this orbital configuration, no secondary eclipse is expected. We report that the optically thick UV resonance lines exhibit flattening and shortening of the absorption trough, and weakening of their emission features, as they approach the phase of the expected secondary eclipse, while the continuum UV flux appears to remain unaffected. We illustrate that this can be explained by the primary's optically thick wind eclipsing the secondary star. The C IV and N V resonance line formation regions in the primary star extend up to 316 Rsol. The measured extend of the formation regions of resonance lines in a stellar wind are important benchmarks for 1D as well as 3D non-LTE stellar atmosphere models. A first comparison to 1D-stellar atmosphere models indicates that a classical beta-law with an exponent of beta=0.5 instead of beta=0.8 might be favoured for the primary star's velocity field.

The spectral kernel field equation R[k] = T[k] lacks a conservation-law analog. We prove (i) the fixed-point flow is strictly volume-expanding (tr DF > 0), precluding automatic conservation, and (ii) the conservation deficit per mode equals the Hessian stability margin exactly: D_m = -Delta'. Closing the deficit requires a scene-side compensating contribution, which we formalise as the distinction dynamics equation dc/dt = G[c, h_t], with MaxCal-optimal realisation G_opt. On fixed-topology 3D surface graphs we derive a conditional topology-preserving compression theorem: retaining k >= beta_0 + beta_1 modes (under a spectral-ordering assumption) preserves all Betti-number charges; we include a worked short-cycle counterexample (figure-eight) calibrating when the assumption fails. A triple necessary spectral diagnostic – Fiedler-mode concentration, elevated curl energy, anomalous beta_1 – is derived for planetary drainage networks at O(N) cost. Two internal real-data sequences serve as preliminary consistency checks; full benchmarks and adaptive-topology extensions are deferred.

We present a covariant and gauge-invariant formulation of the theory of radial adiabatic linear perturbations of self-gravitating, non-dissipative imperfect fluids within the theory of general relativity. By codifying the thermodynamical properties of the source into an equation of state and an ansatz on anisotropic pressure that involves both matter and kinematic variables, we obtain a set of equations that is directly applicable to a wide variety of thermodynamic theories for matter fields. As examples, we evaluate and compare the predictions of the Eckart theory, the Bemfica-Disconzi-Noronha-Kovtun theory, and the Truncated Israel-Stewart theory on the properties and evolution of radial adiabatic perturbations of stellar compact objects modeled by classical equilibrium solutions. Introducing a new solution of the Einstein field equations, and imposing causality, we propose an upper bound for the maximum compactness of dynamically stable stars with non-trivial radial and tangential pressures.

We use the 3D fully kinetic simulation to study different turbulence modes and turbulence anisotropy of relativistic turbulence in magnetically dominated collisionless plasmas. We extend the method developed by Cho & Lazarian (2002) for decomposing non-relativistic magnetohydrodynamic (MHD) turbulence into Alfvén, fast, and slow modes to the regime of collisionless plasmas. We find that Alfvén and slow modes are anisotropic, following the Goldreich & Sridhar (1995) scaling, while fast modes are isotropic. We observe a larger kinetic energy fraction of fast modes compared to that in the non-relativistic MHD turbulence, suggesting a stronger coupling of Alfvén and fast modes in relativistic magnetized turbulence in collisionless plasmas. We further examine the dynamic alignment and find a weaker scale dependence of the alignment angle than previously proposed. The dominant thermal fluctuations in the kinetic range can cause flattening of the turbulent velocity structure function and weakening of the turbulence anisotropy and dynamic alignment near the kinetic scales.

Astrophysics plays a pivotal role in the quest for axions and axion-like particles, offering guidance to experimental efforts and enabling the investigation of axion properties that cannot be probed otherwise. In this context, the extreme conditions in the interiors of compact stellar objects – such as core-collapse supernovae, neutron stars, and binary neutron star mergers – significantly enhance axion production, providing unparalleled sensitivity to extremely feeble couplings to Standard Model particles. In this context, the techniques of multimessenger astrophysics deepens the understanding of powerful transient events, maximizing the capabilities of current instruments to identify possible signatures of axion emission.

Ultra high energy cosmic rays have been see coming from the direction of the local cosmic void. We use this fact to argue that at least some of these these cosmic rays are relatively light magnetic monopoles and that their relative fraction above 1020 eV can be found from full sky observations.

We present a dynamical solution to the dark matter-baryon coincidence problem based on the neutron portal operator connecting the visible and dark sector asymmetries. This framework is motivated by the possibility that a strongly supercooled dark confinement phase transition accounts for the nano-Hz stochastic gravitational wave signal observed by pulsar timing arrays, while also generating the dark matter and baryon asymmetry in the Universe. We show that the GeV-scale mass of asymmetric dark matter can be naturally correlated with the (multi-)TeV scale cut-off for the neutron portal through its ultraviolet completion. The dark sector is governed by an approximate fixed point and confines once the heavy portal states are integrated out, dynamically generating a scale of $\mathcal{O} ({\rm GeV})$. We analyze both tree and loop-level ultraviolet completions and demonstrate how the resulting confinement scale is linked to the effective neutron portal scale. We also discuss cosmological constraints and experimental prospects in beam dump searches and colliders for probing the neutron portal.

We present a comprehensive statistical study of the radial evolution of solar wind turbulence near Mercury's orbit using long-term magnetic field measurements from the MESSENGER mission. Owing to Mercury's highly elliptical orbit and the spacecraft's repeated, extended residence in the upstream solar wind, the data set provides more than 17,000 hours of observations, enabling robust statistics across well-defined heliocentric distance intervals (0.31-0.47 au). We find that inertial-range spectral slopes remain close to -3/2 throughout Mercury's orbit, showing no significant radial evolution. Combined with low magnetic compressibility, this result indicates a stable, predominantly Alfvenic inertial-range cascade already established here. In contrast, kinetic-range spectral slopes exhibit clear radial evolution, becoming progressively shallower with increasing heliocentric distance, highlighting the greater sensitivity of kinetic-scale turbulence to heliocentric conditions. The ion-scale spectral break frequency decreases with distance in the spacecraft frame, while its normalized form increases relative to the local proton cyclotron frequency, demonstrating that the break is not tied to a single ion scale but reflects evolving local plasma conditions. Magnetic compressibility shows a similar frequency dependence at all distances, with a subtle radial enhancement of compressive fluctuations at kinetic scales. Autocorrelation analysis reveals strong anisotropy, with the correlation times of field-aligned magnetic fluctuations increasing with heliocentric distance, while those of perpendicular fluctuations remain shorter and nearly invariant. Together, these results demonstrate a clear scale-dependent radial evolution of solar wind turbulence near Mercury's orbit, providing new constraints on the development of kinetic processes in the inner heliosphere.

Superradiant instability of ultralight bosons can produce clouds around rotating black holes, whose continuous gravitational wave (GW) emission is a promising observational target. Precise predictions of the signal frequency and its evolution are essential for detecting such continuous GWs. For axions, self-interactions can populate multiple superradiant modes via nonlinear couplings, and GW emission can occur through various channels. To calculate the frequency shifts of GWs emitted through these channels, we employ relativistic perturbation theory based on a bilinear form. We apply this framework to self-interaction effects for the first time, and also revisit the treatment of the self-gravity contribution. Our results provide a simple and unified framework for calculating frequency shifts, including cases in which multiple modes are excited, and are relevant for next-generation GW observations.

There has been a recent revival in understanding the spontaneous scalarization phenomenon in scalar-tensor gravity as a phase transition. Using the tools of the Landau theory, we now know that first-order transitions where scalarization occurs in a discontinuous manner is more prominent than what had been considered in the literature, and this might lead to novel observation channels. However, the examples so far have been restricted to specific quadratic scalar coupling terms and spherically symmetric stars. Here we explore the phase transition structure of scalarization for more general couplings, considering linear as well as quadratic terms in the conformal scaling factor of the theory. Moreover, we also investigate the effect of rotation on the scalarization phase transition. Both of these considerations are natural choices since the coupling in a scalar-tensor theory can appear at all orders, and astrophysical neutron stars commonly have angular momentum. The introduction of linear coupling leads to a complex solution space which is harder to explore. However, we demonstrate that the Landau model of scalarization enables us to systematically find the branches of scalarized solutions that are commonly overlooked in numerical searches, providing a novel tool. On the other hand, the main effect of stellar rotation is shifting the stellar masses at which the phase transition occurs to higher values, but the qualitative picture remains similar to what happens under spherical symmetry.

Since the dawn of the space age, observations of energetic particles in planetary radiation belts have been interpreted within a diffusive transport framework, even though the processes that populate and deplete these belts produce highly structured and spatially localized distributions. This exposes a fundamental problem: how can coherent phase-space structures evolving under collisionless dynamics give rise to observational signatures that appear consistent with diffusion-based transport? Here we show that diffusion-like behaviour can arise from an observational phase-mixing effect, independent of stochastic wave-particle transport. As spacecraft sample neighbouring drift shells while particles undergo electromagnetic drifts, spatially localized drift-phase structures are converted into rapidly decorrelating temporal signals, making them observationally indistinguishable from stochastic processes. We show that the effective lifetime of these structures is only a few drift periods, preventing the resolution of fine-scale structure. These results demonstrate that collisionless dynamics can mimic diffusive transport on short timescales, limiting the inference of particle acceleration processes and biasing transport estimates. This calls for a reassessment of diffusion-based interpretations of radiation belts at Earth, across the solar system, and in the radiation belts of ultra-cool brown dwarfs.

The recent detection of a stochastic gravitational wave background by pulsar timing arrays has opened a new window in understanding supermassive black hole binaries and in probing the universe at the early time. Recently, pulsar timing array (PTA) collaborations have been further paving the way to probe anisotropies in the stochastic gravitational wave background. This study investigates dipole-type statistical anisotropy in the primordial power spectrum within a phenomenological framework. We demonstrate that the primordial dipole induces both dipolar and quadrupolar anisotropies in the energy density spectrum of scalar-induced gravitational waves (SIGWs), without generating extra polarization modes. Based on this anisotropic spectrum, we derive the corresponding PTA overlap reduction functions (ORFs), which exhibit frequency dependence, with the anisotropies enhanced on small scales. Furthermore, owing to the non-uniform distribution of millisecond pulsars over the sky in current PTA dataset, the ORFs exhibit a morphology that explicitly depends on the preferred direction of the anisotropy. However, our bayesian analysis of the NANOGrav 15-year dataset still yields no significant evidence for a preferred direction and a weak upper limit on anisotropy amplitude $(g\lesssim0.5)$. This result arises because the observational frequency band lies below the spectral peak, where our models predict suppressed anisotropic contributions. This limitation highlights the potential of future PTA observations. Specifically, datasets with broader frequency coverage are expected to tighten constraints on dipole-type anisotropy.

In this work, we systematically investigate the quantum-information diagnostics of cosmological perturbations with a nontrivial sound speed, utilizing a normalized open two-mode squeezed-state framework. Rather than introducing new observables, our analysis focuses on how a modified sound speed dynamically reshapes the Schrödinger evolution of the squeezing parameters ($r_k$ and $\phi_k$). We demonstrate how these dynamical changes are inherited by the reduced density matrix of the observable sector. By employing a sound-speed-resonance parametrization, we derive and evaluate the purity, von Neumann entropy, Rényi entropies, and logarithmic negativity. To overcome the intrinsic multiscale stiffness of the post-inflationary equations, we introduce a bounded variable $x = \tanh r_k$ as a partial regularization, which enables reliable numerical simulations exclusively within the inflationary regime. Our numerical results reveal that a nontrivial sound speed significantly suppresses the purity of the reduced state, indicating enhanced effective mixedness. Simultaneously, it strongly amplifies and modulates both the entropic and entanglement diagnostics. More precisely, a nontrivial sound speed postpones the onset of classicality by modulating the decoherence process. Ultimately, we show that a nontrivial sound speed leaves distinct and identifiable quantum-information signatures within the entanglement structure of the early universe.

This article opens new window to obtain accelerating scaling attractors without any need of dark energy. We study cosmological dynamics in a two-fluid system where pressureless dark matter (DM) undergoes adiabatic particle creation and exchanges energy with a barotropic fluid. Considering six widely used interaction prescriptions, we formulate the corresponding autonomous systems in a compact phase space and perform a unified dynamical analysis. We find that accelerating scaling attractors, namely late-time states where both fluids coexist with fixed energy fractions, arise only when the interaction is controlled by the DM density and energy flows from DM to the second fluid. Such attractors appear in the global and local DM-based interactions, and in the global mixed case, but are entirely absent when the interaction depends on the second fluid or on local mixed terms, which instead drive the universe to a DM-dominated accelerating phase. These results clarify the unique conditions under which matter creation can mimic dark-energy-like behaviour without introducing a dark-energy component.

Scalar-induced gravitational waves (SIGWs) are ubiquitous in many early-Universe processes accompanied by non-Gaussianity; while Gaussian perturbation can generate significant SIGWs, computations of SIGWs can be significantly affected and enhanced if the scalar perturbations have some degree of non-Gaussianity; hence, precise calculations of these kinds of SIGWs involve a full understanding of non-Gaussianity. In this Letter, we propose to use the lattice simulations to directly calculate the energy density spectra of SIGWs with non-Gaussianity up to all orders. Our proposal has been first verified to match the existing semi-analytical results with non-Gaussianity, and then applied to more general cases, including high-order primordial non-Gaussianities, the logarithmic dependence in curvature perturbations, the curvaton model, and the ultra slow-roll model. We find that even a modest non-Gaussianity can significantly alter ultraviolet behaviors in SIGW spectra, necessitating special cautions in future detections as well as mutual constraints on/from primordial black holes.

Though Rossby waves have been observed on the Sun, their radial eigenfunctions remain a mystery. The prior theoretical work either considers quasi-2D systems, which do not apply to the solar interior, or only considers fully radiative or fully convective atmospheres. This project calculates the radial eigenfunctions for Rossby waves in a deep atmosphere for a general stratification. Here, we use the $\beta$-plane approximation to derive a vertical equation in terms of the Lagrangian pressure fluctuation $\delta P$, and we then calculate radial eigenfunctions for Rossby waves in a standard solar model, Model S. We find that working in the Lagrangian pressure fluctuation results in cleaner wave equations that lack internal singularities that have been encountered in prior work. The resulting radial wave equation makes it abundantly clear that there are two wave cavities in the solar interior, one in the radiative interior and another in the convection zone. Surprisingly, our calculated radial vorticity eigenfunctions for the radiative interior modes are nearly constant throughout the convection zone, raising the possibility that they may be observable at the solar surface.

We present a joint cosmological analysis combining data from the Planck satellite, the Atacama Cosmology Telescope, and the South Pole Telescope. We construct a unified likelihood that reproduces the measured temperature and polarisation power spectra by jointly modelling the cosmic microwave background (CMB) signal, Galactic and extragalactic foregrounds, and instrumental systematics across all datasets. We reduce reliance by combining datasets and improve the robustness of parameter estimation by marginalising over the choice of foreground templates. Within this joint analysis, $\Lambda$CDM parameters exhibit remarkable stability with respect to variations in foreground modelling. Parameters for cosmological extensions are more sensitive to these assumptions, with uncertainties increasing by up to 35% in the neutrino sector after marginalising over foreground models. In contrast, the determination of foreground parameters depends more strongly on the assumptions made about the underlying foreground models. Overall, this work demonstrates the feasibility and reliability of a fully joint analysis of current CMB experiments and emphasizes the importance of consistent and accurate foreground modelling for the scientific goals of next-generation, high-sensitivity CMB surveys.

Observing [O I] emission to calculate an "oxygen line ratio" has been proposed as a potential proxy for direct CO$_2$ measurement in comets. However, the photochemistry governing [O I] release into the coma is not well understood, and using theoretical release rates often yields different results than using empirical release rates determined in conjunction with direct space-based measurements of CO$_2$. We hypothesize that the accuracy of the release rates could depend on the level of solar activity at the time the comet is observed, which will be influenced by the solar cycle. We present observations and analysis of [O I] emission in one comet observed near solar maximum, C/2014 Q2 (Lovejoy), and one near solar minimum, C/2007 N3 (Lulin). Our [O I] measurements were obtained using two high spectral resolution optical spectrographs: the Tull Coudé spectrometer at McDonald Observatory and the ARCES spectrometer at Apache Point Observatory. We use empirical and theoretical models for [O I] emission from the literature to derive multiple sets of inferred CO$_2$ abundances for these comets and compare to contemporaneous space-based measurements of CO$_2$. We find that the empirical model, which was developed based on comet observations obtained near solar maximum, reproduces the directly measured CO$_2$ abundances better for Lovejoy. Neither model accurately reproduces the direct measurement for Lulin. We discuss the implications of our findings for the accuracy and dependencies of the oxygen line ratio method for inferring CO$_2$ abundances in cometary comae.

Several observatories designed to detect ultrahigh-energy neutrinos are planned for the next decade. The most imminent of these is the Payload for Ultrahigh Energy Observations (PUEO), a long-duration balloon-based experiment that will provide unprecedented sensitivity to neutrinos with energies in the range of   1 - 1000 EeV. In this work, we assess the scientific reach of PUEO. In particular, we evaluate the sensitivity of this observatory to cosmogenic neutrinos and, in turn, to the proton fraction of the ultrahigh-energy cosmic-ray spectrum. We also consider the potential of PUEO to probe scenarios in which neutrinos are produced through the decays of ultraheavy dark matter particles or are radiated from cosmic strings. We find that PUEO will be able to constrain the proton composition of ultrahigh-energy cosmic rays in scenarios that feature very strong source evolution and in which protons are accelerated to extremely high energies. Although gamma-ray observations are generally more sensitive to decaying particles than neutrino observations, PUEO is expected to set the strongest neutrino-detector constraints above 10^19 eV. PUEO will also provide the strongest constraints on some models of cosmic strings.

Our understanding of the $\gamma$-ray sky has greatly advanced, yet studying the unresolved $\gamma$-ray background (UGRB) can unveil the nature of the faintest $\gamma$-ray source populations in the Universe. Statistical cross-correlations between the UGRB and tracers of large-scale cosmic structure allow us to infer which sources contribute the most to this emission. In this work, we examine the angular correlation between the UGRB and the matter distribution traced by galaxies, using twelve years of Fermi Large Area Telescope (LAT) observations along with three years of Dark Energy Survey (DES) data. We detect a correlation with a signal-to-noise ratio of 7.85, primarily driven by large angular scales. We then perform a multi-tracer analysis that combines this measurement with the cross-correlation between $\gamma$ rays and DES weak lensing. The two single-tracer results are mutually consistent, and their combination yields a total significance of 10.31, firmly establishing the extragalactic origin of the UGRB. Intriguingly, the properties inferred for the sources contributing to the UGRB show departures from those of the resolved {\gamma}-ray population, suggesting that the faint end of the $\gamma$-ray sky is not a simple extrapolation of currently resolved sources.

Supernovae characterized by enduring narrow optical hydrogen emission lines (SNe IIn) are believed to result primarily from the core-collapse of massive stars undergoing sustained interaction with a dense circumstellar medium (CSM). While the properties of SN IIn progenitors have relatively few direct constraints, the ongoing ejecta-CSM interaction provides unique information about late-stage stellar mass-loss preceding core-collapse. We present late-time X-ray and radio observations of four $\geq$3000 day-old SNe IIn: SN 2013L, SN 2014ab, SN 2015da, and KISS15s. The radio and X-ray emission from KISS15s indicate a mass-loss rate of \eq{\dot M\sim4\times 10^{-3} \rm{M_{\odot} yr^{-1}}} at $\sim$450 years pre-supernova – 2 orders of magnitude below earlier optical estimates (which probed the mass loss immediately preceding the supernova). We find hints of a spectral inversion in the radio SED of KISS15s; a possible signature of a secondary shock due to a binary system or the emergence of a pulsar wind. For SN 2013L, we obtain a mass-loss rate of \eq{\dot M\sim2 \times 10^{-3} \rm{M_{\odot} yr^{-1}}} at $\sim$400 years pre-explosion based on the X-ray detection. For SN 2014ab and SN 2015da, we find a upper limits on the mass-loss rates of \eq{\dot M<2\times10^{-3} M_{\sun} yr^{-1}} explosion at $\sim$300 and 250 years pre-explosion, respectively. All four objects display mass-loss rates lower than estimates from earlier optical analyses by at least 1-2 orders of magnitude, necessitating a rapidly evolving progenitor process over the last centuries pre-explosion. Our analysis reveals how X-ray and radio observations can elucidate progenitor evolution when these objects have faded at optical wavelengths.

The systematic variation of solar active region (AR) properties with their magnetic flux has been the subject of numerous studies but the proposed scaling laws still vary rather widely. A correct representation of these laws and the deviations from them is important for modelling the source term in surface flux transport and dynamo models of space climate variation, and it may also help constrain the subsurface origin of active regions. Here we determine active region scaling laws based on the recently constructed ARISE active region data base listing bipolar ARs for cycle 23, 24 and 25. For the area $A$, pole separation $d$ and tilt angle $\gamma$ we find scalings against magnetic flux $\Phi$ and heliographic latitude $\lambda$. Residuals from these relations are also modelled. These scaling relations are recommended for use in space climate research for the modelling of future data or missing past data, as well as for the identification of candidate rogue ARs.

In this work, we present Fisher matrix forecast of the tomographic weak lensing power spectrum and bispectrum for three physically distinct types of models of beyond-$\Lambda$CDM: the CPL parametrisation of dynamical dark energy, interacting dark energy (IDE) with a dark sector energy-momentum exchange, and Hu-Sawicki models of $f(R)$ gravity. We find that for all three models, including the bispectrum significantly tightens the Fisher constraints: the bispectrum reduces the marginalised $1\sigma$ error on the CPL equation of state parameter from $\sigma(w_0) = 0.2511$ (power spectrum only) to $\sigma(w_0) = 0.1557$, on the IDE coupling from $\sigma(\alpha) = 2.6895$ to $\sigma(\alpha) = 0.2944$, and on the scalaron amplitude from $\sigma(\ln|f_{R0}|) = 2.236$ to $\sigma(\ln|f_{R0}|) = 2.237$ after full marginalisation over nuisance parameters e.g., photo-z error $\sigma_z$ and intrinsic alignment amplitude $\mathcal{A}_{\rm IA}$. We find that $f(R)$ models are the most sensitive to systematics and especially in bispectrum. The results also demonstrates the importance of higher order weak lensing statistics as a practical necessity to maximise the scientific return of Stage IV surveys.

Rotating radio transients (RRATs) represent a significant but poorly understood component of the Galactic neutron star population, characterized by sporadic emission first detectable only through single-pulse searches. We present the RRATalog, an up-to-date catalogue of 335 RRATs, and utilize a uniform sample of RRATs discovered in four Parkes telescope surveys to model their Galactic population. Accounting in detail for observational selection effects, we find a radial density profile similar to pulsars, but identify a significantly steeper luminosity function (power-law index $\alpha \simeq -1.3$) than previously assumed. For sources beaming towards Earth, we estimate $34000 \pm 1600$ potentially observable RRATs above a peak luminosity of 30 mJy kpc$^2$. At these high luminosities, the RRAT population is comparable in size to that of canonical pulsars. Consistent with the observed distribution, the underlying period distribution is significantly shifted toward longer periods compared to canonical pulsars, suggesting RRATs represent a more evolved population. We find evidence for a turnover in the luminosity function below 30 mJy kpc$^2$, and predict that the total number of potentially observable RRATs is $\lesssim 70,000$. Applying the Tauris & Manchester beaming model, we estimate the total Galactic RRAT population to be $\lesssim 400,000$. The implied birth rate of $\lesssim 1.4$ RRATs per century is consistent with the Galactic core-collapse supernova rate, suggesting RRATs can be reconciled with known progenitor rates without requiring a separate evolutionary origin. We provide predictions for RRAT discoveries in ongoing and future surveys.

The upper atmospheres of temperate sub-Neptunes are strongly influenced by atmospheric dynamics due to their cool equilibrium temperature and thereby longer chemical timescales than the atmospheric dynamical timescales. In this study, we used a three-dimensional (3D) general circulation model to investigate the transport-induced disequilibrium chemistry and vertical mixing on temperate gas-rich mini-Neptunes, using K2-18b as an example. We model K2-18b assuming 180 times solar metallicity and consider it as either a synchronous or an asynchronous rotator, exploring spin-orbit resonances of 2:1, 6:1, and 10:1. We find that the vertical transport affects the chemical structure significantly, making CO$_2$ and CO more abundant ($\sim$10$^{-3}$) in the upper atmosphere compared to the chemical equilibrium abundance (<10$^{-15}$), and horizontal winds further homogenize the chemical composition zonally in this region. Molecular abundances in the photosphere generally agree across different rotation periods. We employ a passive tracer in the model to estimate the one-dimensional (1D) equivalent eddy-diffusion coefficient ($K_{zz}$) of K2-18b, providing a parameter useful for future 1D atmospheric models. Additionally, synthetic transmission spectra generated from our model are compared with the JWST observations, and we find that our model can provide a comparable fit to the observations. This work offers a 3D perspective on transport-induced chemistry on a temperate sub-Neptune and derives vertical mixing parameters to support 1D modelling.

We present a set of numerical-relativity binary neutron star merger simulations incorporating muons and muonic reactions for two baseline baryonic equations-of-state. In order to investigate the possible impact of muons and muonic weak reactions, we treat neutrinos with a gray (energy-independent) truncated moments scheme and an implicit-explicit time integrator. Newly computed neutrino rates are employed within the full kinematics approach for a set of relevant reactions, and pair-processes are modeled via opacities computed using reaction kernels, that allow a consistent treatment of neutrino interaction rates. We find that equilibration between matter and radiation is successfully captured by a novel two timescales approach. Of astrophysical interest is the general agreement between our muonic and non-muonic results regarding the remnant evolution, disk and outflow properties. Average electron fractions, asymptotic velocities and temperatures are different by less than $\sim 6\%$, while the main impact of muons is a reduction in ejecta masses by at most $\sim 17\%$. Therefore, based on our findings, accounting for the presence of muons and muonic reactions might result much less severe consequences regarding nucleosynthetic yields and electromagnetic counterparts than previously reported in the literature.

Spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI) and Euclid are mapping the spatial distribution of millions of galaxies, with Emission Line Galaxies (ELGs) serving as the dominant tracer in the redshift range $0.8<z<1.6$. Standard approaches for extracting cosmological information from galaxy clustering, however, typically discard highly constraining measurements from the nonlinear regime. We apply SHAMe-SF - a modification of Subhalo Abundance Matching tailored for star-forming galaxy samples - to analyse the three-dimensional clustering of DESI ELGs from the One-Percent data release, extending their cosmological analysis deep into the nonlinear regime. We validate our pipeline using two mock ELG samples drawn from the state-of-the-art cosmological hydrodynamical simulation MillenniumTNG, demonstrating that our model yields unbiased constraints on $\sigma_8$ and $\Omega_{\rm m}h^2$ down to scales of $0.3~h^{-1}$Mpc on both samples. We find that including scales below $0.8~h^{-1}$Mpc is critical for mitigating projection effects and obtaining unbiased constraints on $\sigma_8$. Applied to the DESI One-Percent measurements, our analysis yields $\sim6$% constraints on $\sigma_8 = 0.81^{+0.05}_{-0.06}$ and $\Omega_{\rm m}h^2=0.146^{+0.009}_{-0.009}$. Remarkably, the accuracy of these constraints is similar to that obtained from the combined full-shape analysis of all DESI DR1 tracers, yet using only 1% of the survey volume. A naive extrapolation of our results from the One-Percent to the full survey area suggests that the complete survey could deliver roughly an order-of-magnitude improvement in precision - a prospect that, while subject to significant practical challenges, illustrates the cosmological potential encoded in the nonlinear regime.

We present high-sensitivity Very Long Baseline Interferometry (VLBI) observations of four ultraluminous X-ray sources (ULXs): Holmberg II X-1, IC 342 X-1, NGC 6946 X-1, and NGC 925 X-1. No compact emission was detected on milliarcsecond scales, with rms noise levels reaching approximately 5–20 $\mu$Jy. The corresponding $5\sigma$ flux density upper limits reach $\sim 26\,\mu\mathrm{Jy}$, implying radio luminosity limits $L_{\rm R} \lesssim 2 \times 10^{33}\,\mathrm{erg\,s^{-1}}$. This disfavors any persistently bright hard-state-like compact core at our sensitivity level. The previously reported VLBI core in Holmberg II X-1 exhibits significant long-term variability, broadly consistent with an overall decline over the past decades. This behavior is consistent with emission from optically-thin ejecta undergoing adiabatic expansion. The VLBI non-detections may reflect intrinsically weak/intermittent compact emission, and/or low–surface–brightness structure that is resolved out by VLBI, and/or absorption/propagation effects such as free–free absorption in dense, ionized winds.

Axions and axion-like particles can be probed through gravitational waves indirectly, often referred to as "audible axions". The usual concept of audible axion relies on the coupling between the axions and the gauge fields. Here we consider an axion-like mechanism with coupling to the Nieh-Yan term. This interaction leads to the direct and efficient production of gravitational waves during the radiation-dominated era, originating from the tachyonic instability of the gravitational perturbations with the Nieh-Yan term. We calculate the energy spectral density of the chiral gravitational wave background and the comoving energy density of axion-like fields. Based on the numerical results, we explore the parameter space of axion masses and decay constants for detectable gravitational wave signals, either in pulsar timing arrays or space-based gravitational wave detections.

A large primordial density perturbation of the Hubble scale will gravitationally collapse, generating an outgoing sound shell, whether or not a primordial black hole (PBH) is formed. In this Letter, we report a hybrid numerical analysis of the stochastic gravitational wave background induced by the collision of sound shells in the early Universe. The peak frequency and amplitude in the GW spectrum depend on the Hubble horizon and the abundance of sound shells. Abundant density perturbations would lead to GW backgrounds potentially detectable for future pulsar timing arrays and ground-based/space-borne detectors. For those perturbations that collapse into PBHs, future null detection of the corresponding high-frequency GW background could put new observational constraints on those PBHs that have already evaporated.

We examine the behaviour of the gauge invariant scalar field perturbations in an analytic inflationary model that transitions from slow-roll to an ultra-slow-roll (USR) phase. We find that the numerical solution of the Mukhanov-Sasaki equation is well described by Hamilton-Jacobi (HJ) theory, as long as the appropriate branches of the Hamilton-Jacobi solutions are invoked: Modes that exit the horizon during the slow-roll phase evolve into the USR as described by the first HJ branch, up to a subdominant $\mathcal{O}(k^2/H^2)$ correction to the Hamilton-Jacobi prediction for their final amplitude that we compute, indicating the influence of neglected gradient terms. Modes that exit during the USR phase are described by a separate HJ branch once they become sufficiently superhorizon, obtained by the shift $\left(\epsilon_1,\epsilon_2\right) \simeq \left(0,-6+\Delta \right) \rightarrow \left(\tilde{\epsilon}_1,\tilde{\epsilon}_2\right)\simeq (0,-\Delta)$ and corresponding to a slow-roll solution (very close to de Sitter) supported by the same potential. This transition is similar to the conveyor belt concept put forward in our previous work [1] and suggests that the limit $\epsilon_2\rightarrow -6$ is unphysical as an asymptotic value for the background/long wavelength solution. We further discuss implications for the validity of the stochastic equations arising from the Hamilton-Jacobi formulation. Our work suggests that if Hamilton-Jacobi attractors are appropriately used, they can successfully describe the dynamics of long wavelength inflationary inhomogeneities for potentials with USR regions.

Gravitational wave signals from binary neutron star (BNS) mergers and binary low-mass black hole (BLMBH) mergers are highly similar in the early inspiral phase. Consequently, the astrophysical origin of recently detected low-mass compact binary coalescences has remained ambiguous, particularly in the absence of electromagnetic counterparts. In this work, we demonstrate that proposed detectors with increased high-frequency sensitivity $-$ including NEMO, Cosmic Explorer, and the Einstein Telescope $-$ will reliably distinguish these two source classes in the late inspiral and postmerger regimes. We further show how these detections can be used to disentangle the individual contributions of BNS and BLMBH systems to the compact binary merger rate, while accounting for misclassification probabilities. Finally, we show this can lead to constraints on the interaction of heavy, non-annihilating dark matter with nucleons. This is achieved by noting that capture of such dark matter particles into neutron stars would lead to transmuted black holes (TBHs), formed via neutron star collapse, which would contribute to the BLMBH rate.

The recently proposed first-order viscous relativistic hydrodynamics formulation by Bemfica, Disconzi, Noronha, and Kovtun (commonly known as the BDNK formulation) has been shown to be causal, stable, strongly hyperbolic, and thus locally well-posed. It is now a viable new option for modelling out-of-equilibrium effects in fluids, and has attracted wide attention in its potential applications to astrophysical systems. In this work, we present the first non-linear numerical simulation of spherically symmetric neutron stars using the BDNK formulation under the Cowling approximation. Using a simplified equation of state, we show that stable evolutions can be constructed within a restricted parameter space up to the simulation time we explored. From these simulations, we analyse the frequency content of the quasi-normal modes and the decay rate of the fundamental mode. This analysis serves as a first step towards constructing a fully consistent model of neutron stars using the BDNK formulation.

The Parity-Doublet Model (PDM) is a chirally invariant effective theory for strong-interaction matter involving nucleons and their opposite-parity partners in a parity-doubling framework. We introduce a multiplicatively renormalizable mean-field approach to include the baryonic vacuum contributions to the resulting grand-canonical potential in an explicitly renormalization-group invariant form. As an application, we evaluate the pertinent thermodynamics of two-flavor symmetric and asymmetric nuclear matter, focusing on the restoration of spontaneously broken chiral symmetry at baryon densities and temperatures relevant for the astrophysics of neutron stars. Special attention is paid to the effect of the baryonic vacuum fluctuations on the evolution of chiral condensate with baryon density and temperature for specific choices of the chirally invariant baryon mass $m_0$ to demonstrate the importance of consistently including these vacuum fluctuations in the PDM.

The final equilibrium stage of stellar evolution can result in either a black hole or a compact object such as a white dwarf or neutron star. In general relativity, both stationary black holes and stationary stellar configurations are known to be axisymmetric, and black hole rigidity has been extended to several higher curvature modifications of gravity. In contrast, no comparable result had previously been established for stationary stars beyond general relativity. In this work we extend the stellar axisymmetry theorem to a broad class of diffeomorphism invariant metric theories. Assuming asymptotic flatness and standard smoothness requirements, we show that the Killing symmetry implied by thermodynamic equilibrium inside the star uniquely extends to the exterior region, thereby enforcing rotational invariance. This demonstrates that axisymmetry of stationary stellar configurations is not a feature peculiar to Einstein gravity but a universal property of generally covariant gravitational theories, persisting even in the presence of higher curvature corrections.

In previous work, we introduced the $\ell$-boson stars, a generalization of standard boson stars, which are parameterized by an angular momentum number $\ell$, while still preserving the spacetime's spherical symmetry. In this article, we present and study the properties of $\ell$-boson stars in spacetimes with a negative cosmological constant, such that they are asymptotically anti-de Sitter.

We formulate within Palatini Supergravity a model of induced-gravity inflation excellently consistent with ACT DR6. The inflaton belongs in the decomposition of a conjugate pair of Higgs superfields which lead to the spontaneous breaking of a U(1)B-L symmetry at a scale close to the range (0.102-5.85)x10^16 GeV. The inflaton field is canonically normalized thanks to a real and shift-symmetric contribution into the Kaehler potential. It also includes two separate holomorphic and antiholomorphic logarithmic terms, the argument of which can be interpreted as the coupling of the inflaton to the Ricci scalar. The attainment of inflation allows for subplanckian inflaton values and energy scales below the cut-off scale of the corresponding effective theory. Embedding the model in a B-L extension of the MSSM we show how the mu parameter can be generated and non-thermal leptogenesis can be successfully realized. An outcome of our scheme is split SUSY with gravitino mass in the range (40-60) PeV, which is consistent with the results of LHC on the Higgs boson mass.