Papers for Friday, May 09 2025

As part of the EU-funded Center of Excellence SPACE (Scalable Parallel Astrophysical Codes for Exascale), seven commonly used astrophysics simulation codes are being optimized to exploit exascale computing platforms. Exascale cosmological simulations will produce large amounts of data (i.e. several petabytes) that will soon be waiting to be analyzed, with enormous potential for scientific breakthroughs. Our tool Spherinator enables the reduction of these complex data sets to a low-dimensional space using Generative Deep Learning to understand the morphological structure of simulated galaxies. A spherical latent space allows the HiPSter module to provide explorative visualization using Hierarchical Progressive Surveys (HiPS) in the Aladin software. Here we present a machine-learning workflow covering all stages, from data collection to preprocessing, training, prediction, and final deployment. This workflow ensures full reproducibility by keeping track of the code, data, and environment. Additionally, the workflow allows for scalability in managing a large amount of data and complex pipelines. We use only open source software and standards that align with the FAIR (Findability, Accessibility, Interoperability and Reproducibility) principles. In this way, we are able to distribute the workflow reliably and enable collaboration by sharing code, data, and results efficiently.

Star formation (SF) in the interstellar medium (ISM) is fundamental to understanding galaxy evolution and planet formation. However, efforts to develop closed-form analytic expressions that link SF with key influencing physical variables, such as gas density and turbulence, remain challenging. In this work, we leverage recent advancements in machine learning (ML) and use symbolic regression (SR) techniques to produce the first data-driven, ML-discovered analytic expressions for SF using the publicly available FIRE-2 simulation suites. Employing a pipeline based on training the genetic algorithm of SR from an open software package called PySR, in tandem with a custom loss function and a model selection technique which compares candidate equations to analytic approaches to describing SF, we produce symbolic representations of a predictive model for the star formation rate surface density ($\Sigma_\mathrm{SFR}$) averaged over both 10 Myr and 100 Myr based on eight extracted variables from FIRE-2 galaxies. The resulting model that PySR finds best describes SF, on both averaging timescales, features equations that incorporates the surface density of gas, $\Sigma_\mathrm{gas}$, the velocity dispersion of gas $\sigma_{\mathrm{gas,~z}}$ and the surface density of stars $\Sigma_\mathrm{*}$. Furthermore, we find that the equations found for the longer SFR timescale all converge to a scaling-relation-like equation, all of which also closely capture the intrinsic physical scatter of the data within the Kennicutt-Schmidt (KS) plane. This observed convergence to physically interpretable scaling relations at longer SFR timescales demonstrates that our method successfully identifies robust physical relationships rather than fitting to stochastic fluctuations.

The ZTF DR2 includes light curves of 3628 Type Ia supernovae (SNe Ia), making it the largest low-redshift SNe Ia sample available. One central question for analyses of SNe Ia is whether the remaining diversity of standardized luminosities arises in part from an intrinsic or extrinsic effect; characterized by the color-independent bimodality in progenitor population or color- and host- dependent diversity in dust extinction, respectively. In the initial analyses of the volume-limited subset (z < 0.06; 945 SNe) of this sample from the ZTF collaboration, the authors reported evidences for the former hypothesis, in contrast to many of the previous evidences for the dust hypothesis found across other largest high- and low- redshift SNe Ia samples. We re-analyze the volume-limited ZTF SNe Ia in the same manner that previous samples were analyzed and report consistency with trends seen in literature samples and in support of the dust hypothesis. We find the following: 1. a color dependency in the canonical `mass-step' size for SNe, with red SNe having a larger host-dependent residual step than blue SNe by 0.18 $\pm$ 0.09 mag; 2. a color-dependent difference in the Hubble residual scatter, with red SNe having a $\sim$ 33% larger scatter than blue SNe at > 3$\sigma$ significance; 3. data's preference of a model that accounts for color-dependency over a simple `step' model at $\sim$ 3$\sigma$; 4. the strongest evidence to date (3.5$\sigma$) that the relationship between SN color, host-galaxy properties, and luminosity is continuous rather than characterized by a discrete step. Accounting for 3 and 4 with our new model, Host2D, yields a 4.0$\sigma$ improvement over the mass-step model. We trace the difference in reported findings to the fitting and analysis methods, in particular the model complexity allowed for the color-luminosity relation, rather than a difference in the sample itself.

The Cosmic Microwave Background (CMB) anisotropies, corrected for foreground effects, form the foundation of cosmology and support the Big Bang model. A previously overlooked foreground component is the formation of massive early-type galaxies (ETGs), which can no longer be ignored, particularly in light of JWST's detection of massive, evolved systems at extreme redshifts (z > 13). The rapid formation of massive ETGs has been advocated in galaxy evolution studies for decades, and recent evidence has compelled even proponents of hierarchical mass assembly to acknowledge the fact that massive ETGs evolve quickly. Constraints from chemical evolution are particularly stringent. Without both intense star formation and a top-heavy galaxy-wide initial mass function of stars (IMF), it is difficult to reconcile stellar population synthesis models with the high metallicity and abundance patterns of alpha elements. We infer from previous studies that the progenitor cloud of each massive ETG must have had a radius of approximately 400 kpc. Comparing this value to the average present-day separation of massive ETGs, their formation may have occurred around 15 < z < 20. We consider this epoch of formation in a flat-LCDM cosmological context, incorporating the known and necessary properties of massive ETGs. Such properties are encapsulated independently by the integrated galaxy-wide IMF (IGIMF) theory. The massive ETG evolution presented in this work is consistent with recent advancements in stellar and galaxy evolution, and is derived entirely without priors or constraints from the CMB. Yet, it emerges as a non-negligible source of CMB foreground contamination. Even in our most conservative estimates, massive ETGs account for 1.4% up to the full present-day CMB energy density.

The Euclid survey aims to measure the spectroscopic redshift of emission-line galaxies by identifying the H$\,{\alpha}$ line in their slitless spectra. This method is sensitive to the signal-to-noise ratio of the line, as noise fluctuations or other strong emission lines can be misidentified as H$\,{\alpha}$, depending on redshift. These effects lead to catastrophic redshift errors and the inclusion of interlopers in the sample. We forecast the impact of such redshift errors on galaxy clustering measurements. In particular, we study the effect of interloper contamination on the two-point correlation function (2PCF), the growth rate of structures, and the Alcock-Paczynski (AP) parameters. We analyze 1000 synthetic spectroscopic catalogues, the EuclidLargeMocks, designed to match the area and selection function of the Data Release 1 (DR1) sample. We estimate the 2PCF of the contaminated catalogues, isolating contributions from correctly identified galaxies and from interlopers. We explore different models with increasing complexity to describe the measured 2PCF at fixed cosmology. Finally, we perform a cosmological inference and evaluate the systematic error on the inferred $f\sigma_8$, $\alpha_{\parallel}$ and $\alpha_{\perp}$ values associated with different models. Our results demonstrate that a minimal modelling approach, which only accounts for an attenuation of the clustering signal regardless of the type of contaminants, is sufficient to recover the correct values of $f\sigma_8$, $\alpha_{\parallel}$, and $\alpha_{\perp}$ at DR1. The accuracy and precision of the estimated AP parameters are largely insensitive to the presence of interlopers. The adoption of a minimal model induces a 1%-3% systematic error on the growth rate of structure estimation, depending on the redshift. However, this error remains smaller than the statistical error expected for the Euclid DR1 analysis.

Young supernova remnants (SNRs) provide crucial insights into explosive nucleosynthesis products and their velocity distribution soon after the explosion. However, these velocities are influenced by the dynamics of the circumstellar medium (CSM), which originates from the progenitor's late-phase mass loss. Cas A, the youngest known Galactic core-collapse SNR, was studied to analyze the spatial distribution of Si and S radial velocities using two high-spectral resolution observations from the XRISM-Resolve imaging this http URL's capabilities enabled the detailed characterization of Si XIII, Si XIV, S XV, and S XVI lines, whose line shapes can be resolved and modeled using Gaussian radial-velocity components. The radial velocities measured generally align with previous CCD-based results, confirming that they were not artifacts caused by blended lines or ionization variations. Modeling line profiles with two-component Gaussians improved fits in some regions, revealing distinct redshifted (backside) and blueshifted (frontside) components only in a few specific areas. In most regions, however, both components were either both redshifted (northwest) or both blueshifted (southeast), consistent with the patchy ejecta shell morphology seen in optically emitting fast-moving knots. The individual line components revealed a line broadening ranging from $\sigma_v \approx 200$ to $\sigma_v \approx 2000$ km/s. Components with $1000 \lesssim \sigma_v \lesssim 2000$km/s are consistent with previously determined reverse shock velocities, suggesting non-equilibrated or partially equilibrated ion temperatures. Narrow components with small radial velocities found near Cas A's projected center likely originate from shocked CSM plasma. But the low radial velocity and small $\sigma_v$ defies identifying these components with either the frontside or backside of the SNR, or both.

After its discovery in 2016, the white dwarf binary AR Scorpii (AR Sco) remained for several years the only white dwarf system to show pulsed radio emission associated with a fast-spinning white dwarf. The evolutionary origin and the emission mechanism for AR Sco are not completely understood, with different models proposed. Testing and improving these models requires observational input. Here we report the results of a targeted search for other binary white dwarf pulsars like AR Sco. Using data from Gaia and WISE, we identified 56 candidate systems with similar properties to AR Sco, of which 26 were previously uncharacterised. These were subject to spectroscopic and photometric follow-up observations. Aside from one new binary white dwarf pulsar found, J191213.72-441045.1, which was reported in a separate work, we find no other systems whose characteristics are akin to AR Sco. The newly characterised systems are primarily young stellar objects (with 10 found) or cataclysmic variables (7 identifications), with the remaining being either blended or non-variable on short timescales.

We present CLUMI+, a self-consistent, multi-probe methodology for reconstructing the mass distribution in and around galaxy clusters by combining gravitational lensing and dynamical observations. Building on the joint-likelihood framework of Umetsu (2013), CLUMI+ integrates weak-lensing shear and magnification data with projected escape velocity measurements in the cluster infall region, yielding tighter constraints on the gravitational potential without relying on equilibrium assumptions. The mass distribution is modeled using a flexible, piecewise-defined convergence profile that characterizes the azimuthally averaged surface mass density within the lensing field, transitioning to a projected power-law form at larger radii where phase-space constraints complement lensing. Additional strong-lensing constraints are incorporated via central aperture-mass measurements, enabling full-scale mass reconstruction from the cluster core to the outskirts. We validate CLUMI+ using synthetic weak-lensing and phase-space data for a massive cluster from the IllustrisTNG simulations, demonstrating unbiased recovery of projected and three-dimensional mass profiles and achieving 10%--30% improvement in precision at large radii. As a case study, we apply CLUMI+ to Abell 2261, combining Subaru and Hubble Space Telescope weak+strong lensing data with spectroscopic measurements from the Hectospec Cluster Survey. This analysis demonstrates the power of multi-probe, equilibrium-free modeling for robust cluster mass reconstruction.

We present analysis of the plateau and late-time phase properties of a sample of 39 Type II supernovae (SNe II) that show narrow, transient, high-ionization emission lines (i.e., "IIn-like") in their early-time spectra from interaction with confined, dense circumstellar material (CSM). Originally presented by Jacobson-Galán et al 2024a, this sample also includes multicolor light curves and spectra extending to late-time phases of 35 SNe with no evidence for IIn-like features at <2 days after first light. We measure photospheric phase light-curve properties for the distance-corrected sample and find that SNe II with IIn-like features have significantly higher luminosities and decline rates at +50 days than the comparison sample, which could be connected to inflated progenitor radii, lower ejecta mass, and/or persistent CSM interaction. However, we find no statistical evidence that the measured plateau durations and $^{56}$Ni masses of SNe II with and without IIn-like features arise from different distributions. We estimate progenitor zero-age main sequence (ZAMS) masses for all SNe with nebular spectroscopy through spectral model comparisons and find that most objects, both with and without IIn-like features, are consistent with progenitor masses <12.5 M$_{\odot}$. Combining progenitor ZAMS masses with CSM densities inferred from early-time spectra suggests multiple channels for enhanced mass loss in the final years before core collapse such as a convection-driven chromosphere or binary interaction. Finally, we find spectroscopic evidence for ongoing ejecta-CSM interaction at radii $>10^{16}$ cm, consistent with substantial progenitor mass-loss rates of $\sim 10^{-4}$--$10^{-5}$ M$_{\odot}$ yr$^{-1}$ ($v_w < 50$ km/s) in the final centuries to millennia before explosion.

Accretion rates from protoplanetary disks onto forming stars are a key ingredient in star formation and protoplanetary disk evolution. Extensive efforts surveying individual star forming regions with spectroscopy and narrow-band photometry have been performed to derive accretion rates on large populations of young stellar objects (YSOs). We use Gaia DR3 XP spectra to perform the first all-sky homogeneous analysis of YSO accretion within 500 pc. We characterise the H$\alpha$ line emission of YSOs by using the H$\alpha$ pseudo-equivalent widths and XP spectra from Gaia DR3. We derive accretion luminosities, mass accretion rates and stellar parameters for 145 975 candidate YSO H$\alpha$ emitters all-sky. We describe filtering strategies to select specific sub-samples of YSOs from this catalogue. We identify a large population of low-accreting YSO candidates untraced by previous surveys. The population of low accreting YSOs is mostly spatially dispersed, away from star forming regions or more clustered environments of star formation. Many YSOs appear disconnected from young populations, reminiscent of 'Peter Pan' YSOs. We find $L_{acc}\propto L_\star^{1.41\pm0.02}$ and $\dot M_{acc}\propto M_\star^{2.4\pm0.1}$ for the purest all-sky sample of YSO candidates. By fitting an exponential to the fraction of accreting stars in clusters of different ages in the Sco-Cen complex, we obtain an accretion timescale of 2.7$\pm$0.4 Myr. The percentage of accretors found by fitting a power-law is 70% at 2 Myr and 2.8% at 10 Myr. With this new catalogue of H$\alpha$ emitters we significantly increase the number of YSO candidates with accretion rate estimations in the local neighbourhood. This allows us to study accretion timescales and the spatial and physical properties of YSO accretion from a large, all-sky, and homogeneous sample for the first time. [abridged]

We present an analysis of the H$\alpha$-emitting ionized gas in the warm phase of the NGC 253 outflow using integral field spectroscopy from the Multi Unit Spectroscopic Explorer (MUSE). In each spaxel, we decompose H$\alpha$, [N II], and [S II] emission lines into a system of up to 3 Gaussian components, accounting for the velocity contributions due to the disk and both intercepted walls of an outflow cone. In the approaching southern lobe of the outflow, we find maximum deprojected outflow velocities down to ~ -500 km/s. Velocity gradients of this outflowing gas range from ~ -350 to -550 km/s/kpc with increasing distance from the nucleus. Additionally, [N II]/H$\alpha$ and [S II]/H$\alpha$ integrated line ratios are suggestive of shocks as the dominant ionization source throughout the wind. Electron densities, inferred from the [S II] doublet, peak at 2100 cm$^{-3}$ near the nucleus and reach $\lesssim 50 $cm$^{-3}$ in the wind. Finally, at an uncertainty of 0.3 dex on the inferred mass of $4\times10^{5}$ M$_{\odot}$, the mass-outflow rate of the H$\alpha$-emitting gas in the southern outflow lobe is ~ 0.4 M$_{\odot}$/year. This yields a mass-loading factor of $\eta$ ~ 0.1 and a ~ 2% starburst energy efficiency.

We revisit the problem of the capture of a primordial black hole (PBH) by a neutron star, accounting for the tidal perturbation from a nearby star or planet. For asteroid-mass PBHs, which could constitute all of the dark matter in the universe, a weakly bound post-capture orbit could be tidally disturbed to the point of preventing the PBH from settling in the neutron star and consuming it within a cosmologically short timescale. We show how this effect depends on environmental parameters and can weaken the proposed constraints based on observations of old neutron stars in high-density dark matter environments for PBH masses $\lesssim 10^{22}\,$g. We also provide approximate analytical formulae for the capture rates.

The cosmological redshift drift promises to be the first observable directly measuring the cosmic expansion rate and should be detectable with upcoming surveys by the Square Kilometre Array and the Extremely Large Telescope. To prepare for these upcoming measurements we study the redshift drift in detail using the relativistic N-body code $\texttt{gevolution}$ focusing on inhomogeneity-induced fluctuations. Using a ray-tracer, we calculate the redshift drift directly from the light cone at two different time steps. To investigate observer-dependent biases we consider 10 different observers. We find that inhomogeneity-induced fluctuations in the redshift drift can in extreme cases be of the same order as the cosmic signal for $z\lesssim0.15$. By comparing our results to first-order perturbation theory, we find that the extreme outliers are due to peculiar motion in over-densities and can be described by first-order perturbation theory to percent precision. We calculate angular power spectra that fit very well with our theoretical predictions based on perturbation theory at linear scales and show a surprisingly large non-linear signal. This shows that redshift drift not only has the power to measure the background expansion, but could also deliver information about the velocity and acceleration fields in clusters.

The long-term dynamical future of the Sun's planets has been simulated and statistically analyzed in great detail, but most prior work considers the solar system as completely isolated, neglecting the potential influence of field star passages. To understand the dynamical significance of field star encounters, we simulate several thousand realizations of the modern solar system in the presence of passing field stars for 5 Gyrs. We find that the impulse gradient of the strongest stellar encounter largely determines the net dynamical effect of field stars. Because the expected strength of such an encounter is uncertain by multiple orders of magnitude, the possible significance of field stars can be large. Our simulations indicate that isolated models of the solar system can underestimate the degree of our giant planets' future secular orbital changes by over an order of magnitude. In addition, our planets and Pluto are significantly less stable than previously thought. Field stars transform Pluto from a completely stable object over 5 Gyrs to one with a ~5% instability probability. Furthermore, field stars increase the odds of Mercury's instability by ~50-80%. We also find a ~0.3% chance that Mars will be lost through collision or ejection and a ~0.2% probability that Earth will be involved in a planetary collision or ejected. Compared to previously studied instabilities in isolated solar systems models, those induced by field stars are much more likely to involve the loss of multiple planets. In addition, they typically happen sooner in our solar system's future, making field star passages the most likely cause of instability for the next 4-4.5 Gyrs.

We report observations of comet C/2014 UN271 (Bernardinelli-Bernstein) carried out on UT 2024 March 8 and 17 at a heliocentric distance ($r_{\mathrm{H}}$) of 16.6 au using the Atacama Large Millimeter/Submillimeter Array (ALMA). The CO ($J$=2-1) line at 230 GHz was detected along with continuum emission from its dust coma and large ($\sim$140 km) nucleus, revealing the nature of the activity drivers and outgassing kinematics of the largest Oort cloud comet discovered to date. This work presents spectrally integrated flux maps, autocorrelation spectra, production rates, and parent scale lengths for CO and a stringent upper limit for the H$_2$CO production rate. CO outgassing displays multiple active jets which evolve from one epoch to the next. The continuum emission is compact and spatially unresolved, and is consistent with thermal emission from the large nucleus and a tentative detection of a dust coma. Complementary optical observations provided activity context for the ALMA epochs, indicating that UN271 underwent an outburst in late February before returning to a quiescent brightness in mid--late March. These results represent the first secure detection of molecular activity reported in the literature for C/2014 UN271 and highlight the dynamic nature of this distantly active small world.

Double Neutron Stars (DNSs) are unique probes to study various aspects of modern astrophysics. Recent discoveries have confirmed direct connections between DNSs and supernova explosions. This provides valuable information about the evolutionary history of these systems, especially regarding whether the second-born Neutron Star (NS) originated from either a Core-Collapse ($CC$) or Electron-Capture Supernovae ($ECSNe$) event. The provided scale diagram illustrates the distribution of different types of DNSs on the basis of their orbital parameters and other factors, including mass loss. As a result, the physical processes in DNSs vary depending on the formation mechanisms of the second-born NS and characteristics of the systems. $ECSNe$ processes are typically associated with merging systems ($e\times{P_{orb}}< 0.05$), while $CC$ processes are more commonly linked to non-merging systems ($e\times{P_{orb}}> 0.05$). Our results suggest a critical mass threshold of 1.30$M_\odot \pm 0.22M_\odot$ (critical value) for the $ECSNe$ process to form an NS, while $CC$ processes might occur at higher masses. Examining the orbital parameters of DNSs in a known gravitational potential can enhance our understanding of the theoretical predictions for DNS progenitor characteristics. It turns out that the $ECSNe$ process predominantly produces DNS systems with short orbital ($P_{orb} \leq 0.25 d$), nearly circular orbits ($e\simeq 0.2$), accompanied by minimal kick velocities imparted on the proto-NS and significant mass loss. In contrast, their orbital dynamics in a known gravitational potential plays a crucial role in enhancing our understanding of the SNe geometry and the formation and evolution processes among different NS samples.

Using the Mikhail Pavlinsky ART-XC onboard the SRG observatory we have detected, for the first time, X-ray pulsations with a period of ~106 s from the poorly-studied high-mass X-ray binary RX J0535.0-6700 located in the Large Magellanic Cloud (LMC), thus proving that the accretor is a neutron star with strong magnetic field. Pulsations with similar period were also found in archival archival data from Chandra and XMM-Newton telescopes. Using photometry from WISE we shown that the source demonstrate significant variability in IR during the last twenty years, which could be caused by a secular evolution of the decretion disk. This discovery makes RX J0535.0-6700 another member of the large family of X-ray pulsars with Be-type companions in the LMC.

Gravitational waves (GWs) from coalescing binary black holes (BBHs) can come from different environments. GWs interact gravitationally with astrophysical objects, which makes it possible to use gravitational lensing by nearby objects (self-lensing) to learn about their environments. We quantify the probability of self-lensing through the optical depth $\tau$ for the main channels of detectable GWs at frequencies $f_{\rm GW}\sim (1-10^3)\,{\rm Hz}$. We then analyze the detectability of the lensing effect (imprint). In star clusters, the probability of self-lensing by stellar-mass black holes (BHs) is low ($\tau\simeq10^{-7}$), even when taking into account nearby BHs in resonant interactions ($\tau\simeq 10^{-5}$). Additionally, the lensing imprint of a stellar-mass lens (diffraction and interference) is too marginal to be detectable by the LIGO-Virgo-KAGRA detectors and most Einstein Telescope signals. For a massive BH lens in the center of a cluster, the probability can reach $\tau\simeq 10^{-4}$ either via von Zeipel-Lidov-Kozai induced mergers of BBHs orbiting a central massive BH or BBHs formed as GW captures in single-single interactions in the Bahcall-Wolf cusp of a nuclear cluster, likely eccentric. For self-lensing by a supermassive BH for BBHs in the migration trap of an AGN (active galactic nucleus) disk, $\tau \simeq 10^{-2}$. The imprint of these massive lenses are two images that are easily detectable already in current detectors. Moreover, AGN disk merger signals have a distinct linear $h_+$ polarization. The probability depends on the extent of the detectability through the threshold impact parameter $y_{\rm max}$, which can increase for future detectors. We conclude that constraining the environment of BBHs is possible by combining self-lensing imprints with other waveform signatures such as eccentricity and polarization.

Current, realistic numerical simulations of the solar atmosphere reproduce observations in a statistical sense; they do not replicate observations such as a movie of solar granulation. Inversions on the other hand reproduce observations by design, but the resulting models are often not physically self-consistent. Physics-informed neural networks (PINNs) offer a new approach to solving the time-dependent radiative hydrodynamics equations and matching observations as boundary conditions. PINNs approximate the solution of the integro-differential equations with a deep neural network. The parameters of this network are determined by minimizing the residuals with respect to the physics equations and the observations. The resulting models are continuous in all dimensions, can zoom into local areas of interest in space and time, and provide information on physical parameters that are not necessarily directly observed such as horizontal velocities. Here we present the first proof of concept of this novel approach, explain the underlying methodology in detail, and provide an outlook to the many applications that PINNs enable.

Given the incomplete sampling of spatial frequencies by radio interferometers, achieving precise restoration of astrophysical information remains challenging. To address this ill-posed problem, compressive sensing(CS) provides a robust framework for stable and unique recovery of sky brightness distributions in noisy environments, contingent upon satisfying specific conditions. We explore the applicability of CS theory and find that for radio interferometric telescopes, the conditions can be simplified to sparse representation. {Building on this insight, we develop a deep dictionary (realized through a convolutional neural network), which is designed to be multi-resolution and overcomplete, to achieve sparse representation and integrate it within the CS framework. The resulting method is a novel, fully interpretable unsupervised learning approach that combines} the mathematical rigor of CS with the expressive power of deep neural networks, effectively bridging the gap between deep learning and classical dictionary methods. {During the deconvolution process, the model image and the deep dictionary are updated alternatively.} This approach enables efficient and accurate recovery of extended sources with complex morphologies from noisy measurements. Comparative analyses with state-of-the-art algorithms demonstrate the outstanding performance of our method, i.e., achieving a dynamic range (DR) nearly 45 to 100 times higher than that of multiscale CLEAN (MS-CLEAN).

Because Cepheid variable stars have long been used as a cosmic benchmark, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. We perform 2D hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations do not model the radial pulsations but focus on the interior structure of Cepheids, which involves an interior convective shell and a convective envelope. We develop a new statistical analysis to examine overshooting in this inner convection zone. Using the extreme value theory, we find that overshooting above the convective shell fills the space between these convectively unstable layers. We develop a new statistical analysis that provides a clearer picture of how overshooting fills this layer, and also allows us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of Gamma distributions. Using our mixture model, we show that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We propose a new form for the diffusion coefficient, which addresses the need for overlapping overshooting layers between convective shells. We introduce the idea of super-mixing layer where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

The gamma-ray burst (GRB) GRB 211211A and GRB 060614, believed to originate from the merger of compact objects, exhibit similarities to the jetted tidal disruption event (TDE) Sw J1644+57, by showing violent variabilities in the light-curve during the decay phase. Previous studies suggest that such fluctuations in TDE may arise from the fallback of tidal disrupted debris. In this paper, we introduce the fluctuations of the mass distribution ${\rm d}M/{\rm d}E$ for the debris ejected during the tidal disruption (with energy $E$) and study their impact on jet power. Turbulence induced by tidal force and the self-gravity of the debris may imprint variabilities in ${\rm d}M/{\rm d}E$ during fallback. We model these fluctuations with a power density spectrum $\propto f_{\rm E}^{\beta}$, where $f_{\rm E} = 1/E$ and $\beta$ is the power-law index. We find that the resulting light curve can preserve the fluctuation characteristics from ${\rm d}M/{\rm d}E$. In addition, the observed fluctuations in the light-curves can be reproduced for a given suitable $\beta$. Based on the observations, we find that the value of $\beta$ should be around $-1$.

In \citet{Guillemot2023} we presented a new method for calibrating pulsar observations conducted with the Nançay decimetric Radio Telescope (NRT), which significantly improved NRT polarimetric measurements and pulsar timing quality for data taken after this method was developed, in November 2019. Results hinted at a dependence of the polarimetric response of the NRT on the observed direction. We investigated this potential dependence, since unaccounted variations of the instrumental response could degrade polarimetric measurements. Additionally, we aimed to develop a method for properly calibrating NRT pulsar observations conducted before November 2019. We conducted three series of observations of bright pulsars over wide declination ranges, in a special observation mode in which the feed horn rotates by $\sim$ 180$^\circ$ degrees across the observation, enabling us to determine the full polarimetric response of the NRT while modeling potential variations of calibration parameters with hour angle and declination. In addition, we used the METM technique to improve the calibration of pre-November 2019 data. From the analysis of the series of observations of bright pulsars with horn rotation, we found that the polarimetric response of the NRT does not appear to vary with hour angle or declination. On the other hand, the new METM-based calibration method appears to significantly improve the calibration of pre-November 2019 data. By analyzing NRT data on a selection of millisecond pulsars we found that the new polarimetric profiles are more homogeneous, they generally have larger signal-to-noise ratios, and found that the TOA data for these MSPs are more accurate and contain lower levels of noise, especially when combining the new calibration method with the \textit{Matrix Template Matching} (MTM) method for extracting TOAs from pulsar observations.

Using data from the SRG/eROSITA all-sky X-ray survey and the GAIA-based catalog of 2,209 members of the Pleiades open star cluster, we found 850 X-ray sources associated with the cluster stars. Over 650 of them were detected in X-rays for the first time. At the distance of the Pleiades, the nominal sensitivity of eROSITA corresponds to a luminosity of $L_X \sim 1.6 \cdot 10^{28}$ erg/s in the 0.3-2.3 keV band. The eROSITA sources associated with Pleiades stars have a total luminosity of $L_{X,tot} \sim 1.3 \cdot 10^{32}$ erg/s , a million times greater than the X-ray luminosity of the quiet Sun. Strong X-ray variability, more than 10 times, was recorded for 27 sources. Most of them are known as eruptive optical variables of the dM class. The value of $R_X=log(L_X/L_{bol})$ increases with decreasing effective temperature of the star from $R_X\approx -5$ to $R_X\approx -2$. The distribution of stars over $R_X$ is bimodal, with the left peak at $R_X\sim-4.3$ being formed by stars of FGK classes, and the right peak at $R_X\sim-3.1$ being mainly populated by M-stars. The relation between $R_X$ and the Rossby number $Ro$ depends on the spectral class. For K- and M- stars, at low Rossby numbers $R_X\sim -3$ and depends weakly on $Ro$. At $Ro \gt 0.25$, a rapid drop in $R_X$ is observed for K stars, while in our sample there are no M stars with large Rossby number. Most of F- and G- stars appear to have smaller $R_X\sim -4.5$, however, our sample size is insufficient for a more detailed characterization of their $R_X-Ro$ dependence.

Metal-poor stars are a rare and ancient type of stars; Carbon-Enhanced Metal-Poor (CEMP) stars are a subset of these celestial bodies that show an enrichment of carbon relative to iron. They are believed to be formed from gas polluted by the first generation of stars after the Big Bang and are important objects for studying the early universe, galaxy evolution, and nucleosynthesis. Due to their rarity, the search for metal-poor stars and CEMP stars is a valuable task. This study investigates the search for CEMP stars based on the low-resolution stellar spectra from LAMOST DR8, and proposes a deep learning scheme. From the LAMOST DR8 spectral library, this work discovered 12,766 CEMP star candidates. For ease of reference and use, we provide the estimated parameters $T_\texttt{eff}$, $\log~g$, [Fe/H] and [C/H] for them.

We present the results of our pipeline for discovering strong gravitational lenses in the ongoing Ultraviolet Near-Infrared Optical Northern Survey (UNIONS). We successfully train the deep residual neural network (ResNet) based on CMU-Deeplens architecture, which is designed to detect strong lenses in ground-based imaging surveys. We train on images of real strong lenses and deploy on a sample of 8 million galaxies in areas with full coverage in the g, r, and i filters, the first multi-band search for strong gravitational lenses in UNIONS. Following human inspection and grading, we report the discovery of a total of 1346 new strong lens candidates of which 146 are grade A, 199 grade B, and 1001 grade C. Of these candidates, 283 have lens-galaxy spectroscopic redshifts from the Sloan Digital Sky Survey (SDSS) and an additional 297 from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1). We find 15 of these systems display evidence of both lens and source galaxy redshifts in spectral superposition. We additionally report the spectroscopic confirmation of seven lensed sources in highquality systems, all with z > 2.1, using the Keck Near-Infrared Echelle Spectrograph (NIRES) and Gemini Near-Infrared Spectrograph (GNIRS).

Swift J1727.8-1613 went into outburst in August 2023 and was one of the brightest black hole X-ray binaries (BHXRBs) in recent years, leading to extensive observing campaigns by NICER and Insight-HXMT. The source exhibited strong X-ray variability and showed type-C quasi-periodic oscillations (QPOs) on a wide range of frequencies. The high data quality over a broad range of X-ray energies (0.5-150 keV) enables us to study the energy-dependence of the QPO waveform and the phase lags at the QPO fundamental and second harmonic frequencies. Using the biphase, we find that the QPO waveform is strongly energy-dependent, with energy bands below and above 15-20 keV showing opposite waveform evolution. We interpret the energy-dependence of the waveform as being due to a pivoting spectral component at the second harmonic frequency, with a pivot energy around 15-20 keV. Using the cross-spectrum, we find that the phase lags between energy bands above 7 keV at the QPO fundamental are small, while those at the harmonic frequency are dominated by a separate lag component that extends over a broader range of frequencies and relates to the broadband noise variability. Comparing the energy-dependent results obtained with the bispectrum and the cross-spectrum, we show that these two Fourier products extract different variability components, e.g. the QPO and the broadband noise, at the same frequencies. Finally, we compare Swift J1727.8-1613 to BHXRB MAXI J1535-571 and find that their spectral-timing properties are similar, indicating that these QPO properties may represent a subset of sources.

We present the discovery and analysis of the sixth microlensing two-planet system, KMT-2022-BLG-1818Lb,c, detected by a follow-up program targeting high-magnification events. Both planets are subject to the well-known ''Close/Wide'' degeneracy, although for the first planet, which has a super-Jovian mass ratio of $q_2 \simeq 5\times 10^{-3}$ in both solutions, the Close topology, with a normalized separation of $s\simeq 0.70$, is clearly preferred by $\Delta\chi^2=26$. However, contrary to all previous two-planet microlensing systems, the mass ratio for the second planet, $q_3$, is substantially (factor of $\sim 10$) different for the Close and Wide topologies of the first planet. While this degeneracy is resolved in the present case due to high-cadence follow-up observations, the appearance of this new degeneracy indicates the need for caution in the analysis of future two-planet systems. A Bayesian analysis suggests that the host is likely a K-dwarf star in the Galactic disk. The first planet is probably a super-Jupiter on a Jupiter-like orbit, while the second planet is a Saturn-class planet on either a Mercury-like or Saturn-like orbit.

EXO 2030+375 is a peculiar high-mass X-ray binary that has been exhibiting Type-I outburst activities consistently over the past decades. The phases of outburst peaks are generally stable and occur near the orbital this http URL, significant orbital phase shifts have occasionally been observed in 1995, 2006, and 2016. In this paper, we report another orbital phase shift triggered by a giant outburst in 2021. During this event, the orbital phase of Type-I outburst peaks changed dramatically, moving from near periastron to 20 days after periastron, followed by a gradual recovery. Additionally, for several cycles after the 2021 giant outburst, some Type-I outbursts were missing, with subsequent outbursts occurring every two orbits. We discuss our observations in the frame of the precession of an eccentric Be disk.

One of the most remarkable results from the \emph{James Webb Space Telescope} has been the discovery of a large population of compact sources exhibiting strong broad H$\alpha$ emission, typically interpreted to be low-luminosity broad-line (Type 1) active galactic nuclei (BLAGN). An important question is whether these observations are in tension with galaxy formation models, and if so how? While comparisons have been made using physical properties (i.e.~black hole mass and accretion rate) inferred from observations, these require the use of SED modelling assumptions, or locally inferred scaling relations, which may be unjustified, at least in the distant high-redshift Universe. In this work we take an alternative approach and forward model predictions from the First Light And Reionisation Epoch Simulations (FLARES) suite of cosmological hydrodynamical zoom simulations to predict the observable properties of BLAGN. We achieve this by first coupling \flares\ with the \qsosed\ model to predict the ionising photon luminosities of high-redshift ($z>5$) AGN. To model the observed broad H$\alpha$ emission we then assume a constant conversion factor and covering fraction, and the fraction of AGN that have observable broad-lines. With a reasonable choice of these parameters, \flares\ is able to reproduce observational constraints on the H$\alpha$ luminosity function and equivalent width distribution at $z=5$.

We dissect the kinematics and morphology of the molecular gas within the near-nuclear region of NGC 1068 to understand the mechanisms in the central AGN that might be fueling it, and the impact of its energy output on the surrounding molecular gas. We present high angular and spectral resolution ALMA observations of the HCO$^+$4->3 and CO 3->2 molecular lines in the near-nuclear region of the prototype Seyfert 2 galaxy NGC 1068. The spatial resolution (1.1~pc) is almost two times better than that of previous works studying the same molecular lines at the same transitions and is the highest resolution achievable with ALMA at these frequencies. Our analysis focuses on moment maps, position-velocity (PV) diagrams, and spectra obtained at the position of the nuclear continuum source, along with a simple kinematic model developed using the 3DBarolo software. Our observations reveal significant asymmetry between the eastern and western sides of the nuclear disc in terms of morphology, velocity, and line intensity. The broad lines seen in the inner 2 pc could be accounted for by either beam smearing or highly turbulent gas in this region. Outside this radius the mean velocities drop to $\pm$30 km/s, which cannot be explained by asymmetric drift. We find low velocity connections extending to 13 pc suggesting interactions with larger scale structures. The CO/HCO$^+$ line ratio at the nucleus reported here are extremely low compared to values in the literature of the same galaxy at lower spatial resolutions.

Published radial velocity (RV) measurements of the K giant star 42 Dra reveal variations consistent with a 3.9 M_Jup mass companion in a 479-d orbit. This exoplanet can be confirmed if these variations are long-lived and coherent. Continued monitoring may also reveal other companions. We have acquired additional RV measurements of 42 Dra spanning fifteen years. Periodogram analyses were used to investigate the stability of the planet RV signal. We also investigated variations in the spectral line shapes using the bisector velocity span as well as infrared photometry from the COBE mission. The new RV measurements do not follow the published planet orbit. An orbital solution using the 2004 - 2011 data yields a period and eccentricity consistent with the published values, but the RV amplitude has decreased by a factor of four from the earlier measurements. Including some additional RV measurements taken between 2014 and 2018 reveal the presence of a second period at 530 d. The beating of this period with the one at 479-d may account for the observed amplitude variations. The planet hypothesis is conclusively ruled out by COBE/DIRBE 1.25 micron photometry that shows variations with the planet orbital period as well as an additional 170 d period. The amplitude variations in the RV as well the COBE/DIRBE photometry firmly establish that there is no giant planet around 42 Dra. The presence of multi-periodic variations suggests that these may be stellar oscillations, most likely oscillatory convection modes. These oscillations may account for some of the long period RV variations attributed to planets around K giant stars. This may skew the statistics of planet occurrence around intermediate mass stars. Long-term monitoring with excellent sampling is required to exclude amplitude variations in the long-periods found in radial velocity of K giant stars.

The ionized interstellar medium disperses pulsar radio signals, resulting in a stochastic time-variable delay known as the dispersion measure (DM) noise. In the wideband paradigm of pulsar timing, we measure a DM together with a time of arrival from a pulsar observation to handle frequency-dependent profile evolution, interstellar scintillation, and radio frequency interference more robustly, and to reduce data volumes. In this paper, we derive a method to incorporate arbitrary models of DM variation, including Gaussian process models, in pulsar timing and noise analysis and pulsar timing array analysis. This generalizes the existing method for handling DM noise in wideband datasets.

We design an uncertainty-aware cost-sensitive neural network (UA-CSNet) to estimate metallicities from dereddened and corrected Gaia BP/RP (XP) spectra for giant stars. This method accounts for both stochastic errors in the input spectra and the imbalanced density distribution in [Fe/H] values. With a specialized architecture and training strategy, the UA-CSNet improves the precision of the predicted metallicities, especially for very metal-poor (VMP; $\rm [Fe/H] \leq -2.0$) stars. With the PASTEL catalog as the training sample, our model can estimate metallicities down to $\rm [Fe/H] \sim -4$. We compare our estimates with a number of external catalogs and conduct tests using star clusters, finding overall good agreement. We also confirm that our estimates for VMP stars are unaffected by carbon enhancement. Applying the UA-CSNet, we obtain reliable and precise metallicity estimates for approximately 20 million giant stars, including 360,000 VMP stars and 50,000 extremely metal-poor (EMP; $\rm [Fe/H] \leq -3.0$) stars. The resulting catalog is publicly available at this https URL. This work highlights the potential of low-resolution spectra for metallicity estimation and provides a valuable dataset for studying the formation and chemo-dynamical evolution of our Galaxy.

The Event Horizon Telescope (EHT) is a network of antennas across the globe currently used to image super-massive black holes (SMBHs) at a frequency of 230 GHz. Since the release of the image of M87$^\ast$ in 2019 and, subsequently, that of Sgr A$^\ast$ in 2022 by the EHT collaboration, the focus has shifted to dynamically imaging SMBHs. This has led to a search for potential sites to extend and fill in the gaps within the EHT network. The Gamsberg Mountain and the H.E.S.S. site are both located within the Khomas highlands and have been identified as potential sites for the Africa Millimetre Telescope (AMT). Precipitable water vapour (PWV) in the atmosphere is the main source of opacity and noise from atmospheric emissions when observing at millimetre to sub-millimetre wavelengths. This study aims to establish the PWV content and the atmospheric transmission at 86, 230, and 345 GHz at the AMT potential sites using Global Navigation Satellite System (GNSS) derived PWV data. Results show both sites have potential for observations at 86 and 230 GHz, with 345 GHz possible at the Gamsberg Mountain during winter. The overall median PWV of 14.27 mm and 9.25 mm was calculated at the H.E.S.S. site and the Gamsberg Mountain, respectively. The EHT window had PWV medians of 16.62 mm and 11.20 mm at the H.E.S.S. site and Gamsberg Mountain, respectively. Among the two sites, the Gamsberg Mountain had the lowest PWV conditions, therefore making it the most suitable site for the AMT.

Little red dots (LRDs) are a newly identified class of broad-line active galactic nuclei (AGN) with a distinctive v-shape spectrum characterized by red optical and blue UV continuum emission. Their high abundance at redshifts of $z\sim6-8$ and decline at lower redshifts suggest a transient origin. We propose that the spectral shape of LRDs originates from compact binary black hole systems, where each black hole is surrounded by a mini-disk and embedded in a larger circum-binary disk. With a binary separation of $\lesssim 10^3$ Schwarzschild radii, the Wien tail of a $T\simeq 5000~{\rm K}$ blackbody spectrum at the inner edge of the circum-binary disk produces the red optical emission, while the mini-disks power the UV continuum. Binary torques carve out a gap between the circum-binary disk and mini-disks, setting the turnover wavelength of the v-shaped spectrum around the Balmer limit. This scenario naturally reproduces LRD spectra requiring only modest dust attenuation ($A_V\lesssim 1$ mag), resolving overestimated luminosities for LRDs in previous studies and alleviating a tension with the so-called Soltan argument. This model predicts a distinct spectral evolution as the binary orbit decays through binary-disk interactions and gravitational waves (GWs), linking early-stage "proto-LRD" binaries to the broader AGN population and late-stage "LRD-descendants" to coalescing binaries detectable in GW experiments.

The High Energy Stereoscopic System (H.E.S.S.) site and the Gamsberg Mountain have been identified as potential sites for the Africa Millimetre Telescope (AMT). The AMT is poised to observe at millimetre and possibly at submillimetre wavelengths. At these wavelengths, precipitable water vapour (PWV) in the atmosphere is the main source of opacity during observations and therefore needs to be accurately assessed at the potential sites for the AMT. In order to investigate the PWV conditions for the AMT, identical Global Navigation Satellite System (GNSS) stations were installed and used to assess the PWV at the two potential sites. In this study, the accuracy of those PWV measurements by the GNSS stations was assessed by comparing the H.E.S.S. installed GNSS station PWV measurements to that from a 210 GHz Water Vapour Radiometer (WVR) also installed at the H.E.S.S. site. A correlation of 98% and an offset of 0.34 mm was found between the GNSS station and the 210 GHz WVR PWV data when on-site pressure and the Nevada Geodetic Laboratory (NGL) weighted mean temperature ($\mathrm{T_m}$) were used calculate the GNSS station PWV data. In comparison, the offset reduces to 0.15 mm when on-site derived $\mathrm{T_m}$ and pressure were used to calculate the GNSS station PWV. The results show that the GNSS station with on-site meteorological data can be used with high accuracy to reliably determine the PWV conditions at the H.E.S.S. site.

The ACT data shows an enhancement in the value of the scalar spectral index, as $n_s = 0.9743 \pm 0.0034$, leading to disfavoring many inflationary models, including the Starobinsky model. To satisfy the constraint made by ACT, we will investigate the Starobinsky potential within the Einstein-Gauss-Bonnet(EGB) gravity theory. EGB gravity is motivated by the higher-dimensional theory, which includes quadratic curvature correction terms and a coupling between the scalar field and the GB term that modifies the dynamical equations. The model is considered by using the slow-roll approximation method and the exact numerical approach for two different coupling functions. The results indicate that the model is in good agreement with the data and the results stand in $1\sigma$ of the ACT $r-n_s$ plane. Considering the running of the scalar spectral index also implies the consistency of the model with data. In addition, the parametric space of the free parameters of the EGB coupling is explored, where we find the acceptable region of the parameters in which the resulting $n_s$ and $r$ stay in $1\sigma$ of the ACT data. Next, the reheating phase is considered. It is determined that the model can simultaneously satisfy the constraint of ACT data and the reheating temperature constraints.

Ultra-long period transients (ULPs) are an elusive class of compact objects uncovered by radio surveys. While magnetars are a leading candidate for isolated ULPs, several observational properties challenge the established evolutionary framework: (i) low quiescent X-ray luminosities, (ii) $\sim$~hour-long rotational periods, and (iii) highly-variable radio flux. It is shown via magnetothermal modelling that, if electric currents thread the fluid core at the time of crust freezing, the neutron star remains multiband silent for an initial period of approximately 0.1 Myr while cooling passively. Once the crust becomes cold enough, the Hall effect begins to dominate the magnetic evolution, triggering crustal failures that inject magnetospheric twist that initiates radio pulsing while depleting rotational kinetic energy from an already-slow star. Depending on where electric currents circulate, such 'late-blooming' magnetars manifesting as Galactic ULPs may thus form a distinct branch from soft gamma repeaters and anomalous X-ray pulsars.

We present the combined 2.5$-$30$\mu$m spectra of four protostars acquired with the infrared camera and the infrared spectrograph on board the AKARI and Spitzer space telescopes, respectively. To analyze the ice absorption features in the 8$-$22$\mu$m, we first performed a continuum determination process on mid-infrared spectra and applied a method to subtract the silicate absorption. We conducted a global fitting process to the absorption features in the combined infrared spectra using the experimental ice absorbance data to identify the intrinsic absorption of each ice component. We first derived the H$_{2}$O ice column densities of both stretch and libration modes at 3.05$\mu$m and 13.6$\mu$m simultaneously. We also identified the absorption features containing NH$_{3}$, CH$_{3}$OH, CO$_{2}$, and CO and decomposed their mixed components and compared their ice abundances at different evolutionary stages of the protostars. We explored possible absorptions of the organic ice species such as HCOOH, CH$_{3}$CHO, and CH$_{3}$CH$_{2}$OH in the mid-infrared ranges. The ice analysis method developed in this study can be applied to the ice spectra obtained by the James Webb Space Telescope.

X-ray emission is from shock-heated plasma in magnetic cataclysmic variables, and is processed by absorption and scattering before reaching the observer. We investigate these effects in the X-ray emission of the Intermediate Polar IGR J17195-4100 by carrying out the X-ray spectral analysis, spin modulation and hardness ratio. We present high-sensitivity broadband X-ray spectral analysis by combining NuSTAR and XMM-Newton observations in the 0.3-78.0 keV energy range. The X-ray spectral analysis is performed using six composite models, including the angle-dependent reflection model (reflect), multi-temperature plasma emission models (CEVMKL, MKCFLOW), and photoionised and/or neutral partially covering absorption models (zxipcf, pcfabs) within XSPEC. We examine also the spin modulation in four energy ranges and the hardness ratio, to determine the absorption and scattering effects. We find that the spectrum is best modelled with a reflection amplitude ($\Omega$) of 0.58 $^{+0.38}_{-0.26}$, an ionisation parameter, log($\xi$), of 1.46$^{+0.44}_{-0.23}$ with an equivalent hydrogen column density of 3.09$^{+2.26}_{-0.68}$ $\times$ 10$^{22}$ cm$^{-2}$, a neutral absorber, and a multi-temperature plasma temperature of 27.14$^{+2.0}_{-2.13}$ keV. In addition, we detect effects of electron scattering in the NuSTAR band, leading to modulation amplitude of about a steady 9$\%$ that increases up to 15$\%$ after 20 keV. We stress that these effects significantly affect the X-ray emission of intermediate polars and should be considered to obtain a good representation of the intrinsic spectrum.

We present new, high frequency radio observations of the merging galaxy clusters PLCK G287.0+32.9, Abell 2744, and Bullet. These clusters are known to host $\sim$Mpc scale sources, known as radio halos, which are formed by the acceleration of cosmic rays by turbulence injected into the intracluster medium during cluster mergers. Our new images reveal previously undetected faint outermost regions of halos, extending to over 2 Mpc. This discovery highlights the presence of radio halos with large extents at high frequencies and suggests that their observable size depends on a combination of the observation sensitivity and uv-coverage, and their radio power. We additionally compare the properties of these three clusters with MACS J0717+3745 and Abell 2142, both of which are known to host prominent large radio halos. Remarkably, all five halos, despite their exceptionally large extents, exhibit properties similar to other classical halos: their radial profiles are described by a single-component exponential fit, they show radial spectral index steepening, and have an average radio emissivity of about $10^{-42}\, \mathrm{erg\,s^{-1}\,cm^{-3}\,Hz^{-1}}$. Our results demonstrate that radio halos can extend to the cluster periphery, without the transition to an observationally distinguishable different halo component in the outermost regions. Our findings highlight that careful subtraction of unrelated sources embedded in the halo is necessary to measure the radio surface brightness accurately, as incomplete subtraction can introduce an apparent secondary component in the peripheral regions.

The number density of UV luminous galaxies discovered by the James Webb Space Telescope at ultra high redshift ($z \gtrsim 10$) is higher, and declines much more slowly with increasing redshift, than expected from extrapolations of lower redshift observations or pre-launch physics-based models. Most of these models assume star formation efficiencies (SFE) of only a few percent, motivated by observations of nearby galaxies. In this work, we incorporate a scaling of SFE with gas surface density (which we refer to as Density Modulated SFE; DMSFE), motivated by cloud-scale simulations and theory, into a semi-analytic cosmological model (SAM) of galaxy formation which is calibrated to match the observed rest-UV sizes of high redshift galaxies. We also model the impact of dust and bursty star formation on the SAM-predicted properties of observed galaxies. We show that with plausible values of the main parameters, such as the fraction of gas in dense clouds $f_{\rm dense}$, our new models easily reproduce or even exceed the observed galaxy number densities at $z\sim 6$-17. While no single value of $f_{\rm dense}$ is able to reproduce the very shallow observed decline of the galaxy number density at $z\gtrsim 12$, it is plausible and even expected for $f_{\rm dense}$ to have some effective dependence on cosmic time, which could bring these models into closer agreement with the data. We show that the combined effects of DMSFE, decreasing dust attenuation, and increasingly bursty star formation at earlier cosmic epochs could conspire to reproduce the observed evolution.

X-Ray Flashes (XRFs) are fast X-ray transients discovered by the BeppoSAX satellite, showing an isotropic sky distribution and a prompt emission duration between 10-1000 seconds. The observed prompt X-ray spectrum is similar to Gamma Ray Bursts (GRBs), but with a softer peak energy of the spectrum. Several pieces of evidence indicate that XRFs are connected to GRBs and likely represent their softer analogues, but their origin is still unclear. Several models have been proposed to explain the observed properties of XRF, mostly in the context of collapsar scenario, similar to GRBs but with different geometrical or physical conditions of the progenitor. These include off-axis GRBs and baryon-loaded explosions, that either produce a low Lorentz factor jet or a spherical mildly (or non-) relativistic ejecta, known as cocoons. In this paper, we present multi-wavelength observations of the afterglow of EP241021a, a soft X-ray transient detected by EP, consistent with being an XRF. We present the results of our multiwavelength campaign from radio (uGMRT, ATCA, e-MERLIN, ALMA), optical (LBT, GTC, CAHA) and X-rays (EP-FXT). EP241021a afterglow is characterized by multiple components. They represent the imprints of the interaction of a jet with the complex environment of the pre-existing progenitor, that is likely shaping its structure from a highly relativistic narrow cone to much broader and mildly relativistic cocoon components.

We re-examine the deep-focus methodology of time-distance helioseismology previously used to estimate the power spectrum of the solar convection at a depth of about 30 Mm, which was found to be significantly weaker than predicted by theory and simulations. The Global Acoustic, Linearized Euler (GALE) and Eulerian Lagrangian (EULAG) codes are used to generate ground-truth simulations through which the accuracy of the convective power spectrum can be evaluated. This validation process shows that the power spectrum diverges significantly from ground truth beyond spatial scales corresponding to the spherical harmonic degree $\ell=15$ - $30$ because of the limited resolution of helioseismic measurements. However, the power estimated at larger spatial scales ($\ell<15$) is sufficiently accurate. We then apply the methodology to solar data and find a spectrum that is substantially stronger than previously reported. We discuss some possible differences in methodology that might have led to the initial under-estimation of solar convective power. The new spectra are in line with recent hydrodynamic and magnetohydrodynamic simulations of solar convection and also consistent with the previous inferences obtained by the ring-diagram local-helioseismology method.

ABRACADABRA-10 cm has had great success as a pathfinder lumped-element axion dark matter experiment, setting limits on axion dark matter at the GUT scale. Now, using the interaction of gravitational waves with electrodynamics and a change in readout strategy, we use the ABRA-10 cm detector for the first search for high-frequency gravitational waves using a modified axion detector. Potential sources at these high frequencies (10 kHz to 5 MHz) include merging primordial black hole binaries or superradiance, among other beyond the standard model phenomena. This paper presents the design, results, and challenges from the ABRA-10 cm high-frequency gravitational wave search, showing it is possible to simultaneously look for axions and high-frequency gravitational waves with both searches matching theoretical expectations for sensitivity. Additionally, we conducted the first time series transient search with data from an axion experiment, achieving sensitivity to $10^{-4}$ in strain. Scaled directly to the next generation axion experiment, DMRadio-GUT, this sensitivity would imply a reach to 0.01 $M_{\odot}$ primordial black hole mergers at distances around 3 pc, with prospects to go significantly further with modifications to the readout.

After 45 years since the discovery of quantum-gravitational birth of the cosmological density perturbations we can try to answer the main question of cosmology what is the origin of the Universe. This has become possible because the observational data are precise enough to build a physical model of the cascade relaxation of the gravitating vacuum as a generator of the evolving Universe. We present a model of the early Universe, which is inspired by the observational data in the spatial wavenumbers $k\in(2\cdot10^{-4}, 10)$ Mpc$^{-1}$ and does not require a solution of the Friedman equations. Based on observations we present a solution in the form of a vacuum attractor of General Relativity, which provides an additional power in the form of a ``bump'' and a blue spectrum of perturbations at $k>10\,$Mpc$^{-1}$. Extending the power spectrum to the scale of hundreds of kiloparsecs and less can solve the problems of the observed cosmology -- the appearance of the early star formation and supermassive black holes at $z>10$, the birth of primordial black holes, $\Lambda$CDM model, and others.

We propose a mechanism for relaxing a gauge theory CP violating phase in discrete steps to very small values. The idea is that the CP violating phase includes the magnetic dual of a $4$-form flux which can discharge by the nucleation of membranes. Inside the bubbles surrounded by the membranes, the total CP violating phase is reduced. When the bubbles are produced rapidly during radiation domination in the early universe, near the chiral symmetry breaking scale, they will collide and percolate, melting away into gauge theory radiation and dramatically relaxing CP violation.

Light primordial black holes (PBHs) may have originated in the early Universe, and could contribute to the dark matter in the Universe. Their Hawking evaporation into particles could eventually lead to the production of antinuclei, which propagate and arrive at Earth as cosmic rays with a flux peaked at GeV energies. We revisit here the antiproton and antideuteron signatures from PBH evaporation, relying on a lognormal PBH mass distribution, state-of-the-art propagation models, and an improved coalescence model for fusion into antideuterons. Our predictions are then compared with AMS-02 data on the antiproton flux. We find that the AMS-02 antiproton data severely constrain the Galactic PBH density, setting bounds that depend significantly on the parameters of the lognormal mass distribution, and that are comparable to or slightly stronger than bounds set from diverse messengers. We also discuss prospects for future detection of antideuterons. Given the bounds from AMS-02 antiproton data, we predict that if antideuterons were to be measured by AMS-02 or GAPS, since the secondary contribution is subdominant, they would clearly be a signal of new physics, only part of which could, however, be explained by PBH evaporation.

We improve limits on $Z^\prime$ extensions of the Standard Model (SM) with light right-handed neutrinos. The presence of shared gauge interactions between the light right-handed neutrinos and other SM fermions allows for production of $\nu_R$ in the early Universe and we use the excess in the effective number of neutrino species $\Delta N_\text{eff}$ to place limits. Our benchmark model is a minimal gauged $U(1)_{B-L}$ that often arises as a building block in other models, and we discuss applicability to more general $U(1)$ extensions. We devise an improved Monte Carlo integration scheme convenient for implementation of generic integrated Boltzmann equations with minimal simplifying assumptions. We sketch our numerical implementation in detail for future reference. Using the new ACT DR6 limit $\Delta N_\text{eff}<0.17$, we improve constraints on the gauge coupling for $1\mathrm{\,GeV} < m_{Z^\prime} < 100\mathrm{\,TeV}$ by orders of magnitude and find the strongest limits thus far, surpassing even current and future colliders, and explore the potential of future CMB experiments to test $U(1)$ extensions up to the GUT scale. We perform a detailed analysis of the robustness of cosmological limits within standard and non-standard thermal histories and find that a strong first order phase transition, early dark energy or early matter domination could dilute $\nu_R$ abundances beyond detection. We investigate the effect of reheating on $\nu_R$-genesis and provide results and prescriptions to apply our bounds to non-standard thermal histories. Limits are generically weakened for reheating $T_\text{reh}\ll m_{Z^\prime}$. Our results suggest that projected limits on $Z^\prime$ with Dirac neutrinos can only be accommodated for in non-standard thermal histories, thus limiting the options to include dark matter candidates or Dirac leptogenesis.

Gravitational wave (GW) signals from extreme mass ratio inspirals (EMRIs) are a key observational target for the Laser Interferometer Space Antenna (LISA). The waveforms may be affected by the astrophysical environment surrounding the central black hole (BH), and in particular by the surrounding dark matter (DM) distribution. In this work, we consider the effect of a rotating DM "spike" around a central Kerr BH, and assess its detectability with LISA. Using a fully relativistic model for the rotating spike, we investigate its effect on the inspiral and hence on the emitted GW signals. We compute dephasings and mismatches to quantify how the spin of the primary BH affects the binary dynamics and the gravitational waveform. We show that the modifications due to the spin of the primary BH improve the detection prospects of DM spikes with LISA, and must be taken into account for future parameter estimation studies. We also estimate within post-Newtonian theory how the environment affects the background metric, and show that this effect is mostly negligible for the systems we consider.

We perform a thorough investigation of (ultra)relativistic freeze-out (UFO) during reheating. While the standard WIMP (non-relativistic freeze-out) and FIMP (freeze-in) paradigms have been explored in detail during the reheating epoch, UFO has not been systematically studied, despite the fact that it is operative in a broad region of parameter space. Although dark matter (DM) is ``hot" at the time of relativistic freeze-out, we show that it can easily undergo enough cooling by the time of structure formation to be compatible with $\Lambda$CDM. Unlike standard WIMP-like freeze-out, there can be significant out-of-equilibrium DM production after UFO, similar to the freeze-in mechanism. However, unlike freeze-in, UFO can accommodate much stronger couplings. The UFO parameter space consistent with $\Omega_{\chi}h^2=0.12$ is quite large, with DM masses spanning about 13 orders of magnitude ($10^{-7} \text{ GeV} \lesssim m_{\chi} \lesssim 10^{6}$ GeV), reheating temperatures spanning 17 orders of magnitude ($10^{-2} \text{ GeV} \lesssim T_{\rm RH} \lesssim 10^{15} \text{ GeV}$) and Beyond the Standard Model (BSM) effective interaction scales spanning 11 orders of magnitude ($10^{3} \text{ GeV} \lesssim \Lambda \lesssim 10^{14}\text{ GeV}$). Interestingly, the most suitable range of couplings for UFO lies precisely between the typical couplings for WIMPs and FIMPs, rendering UFO quite attractive from the standpoint of detection. Particle physics models that are easily amenable to UFO include heavy vector or scalar portal interactions, along with nonrenormalizable effective interactions. Finally, we show there is a distinction between UV UFO and IR UFO, where the relic abundance for the former is sensitive to the freeze-out temperature, while the abundance for the latter is sensitive to the DM mass and the reheating temperature but insensitive to the freeze-out temperature.

The detection of gravitational waves from binary black hole and neutron star mergers by ground-based interferometers, as well as the evidence for a gravitational wave background from pulsar timing array experiments, has marked a new era in astrophysics and cosmology. These experiments also have great potential for discovering new physics through gravitational wave detection. One of the most motivated sources of gravitational waves that can be realized only within the beyond-the-Standard-Model framework is first-order phase transitions. In this work we release PT2GWFinder, a Mathematica package designed to compute phase transition parameters and the gravitational wave power spectrum for an arbitrary scalar theory exhibiting a first-order phase transition in a single direction. PT2GWFinder performs the phase tracing, computes the bounce profile and action using FindBounce, calculates the relevant temperatures and phase transition parameters, and finally evaluates the gravitational wave spectrum. Additionally, it offers a user-friendly interface with DRalgo, which enables the computation of the dimensionally reduced effective potential in the high-temperature regime. This work includes a user manual and two models that demonstrate the capability and performance of PT2GWFinder. As a supplement, for one of these models we obtain the bounce solution and action analytically in the thin-wall approximation and demonstrate excellent agreement with the numerical approach.

The left-right model with the simplest Higgs structure has a $SU(2)_L$ doublet $H_L$, a $SU(2)_R$ doublet $H_R$, and a singlet $\sigma$ that couples to the doublets as $\sigma(|H_L|^2-|H_R|^2)$. The left-right symmetry has $H_L\leftrightarrow H_R$ and $\sigma\leftrightarrow-\sigma$. We study gravitational wave production and baryogenesis in the electroweak phase transitions in the model. For two benchmark points, $\sigma$ first gets a vev, followed by $H_R$ and then $H_L$. An interesting feature is that the vev of $H_L$ is initially much higher than its zero temperature value, leading to a more strongly first-order transition and higher frequency gravitational waves. An unusual benchmark point has the $\sigma$ vev at zero when the electroweak symmetry breaks. This results in both $H_L$ and $H_R$ vevs being equal and in the multi-TeV range. At a lower temperature, the $\sigma$ vev turns on, breaking the left-right symmetry and $H_L$ drops to its standard model value. Since the electroweak symmetry is broken at the multi-TeV scale, the frequency of gravitational waves will be much higher than usual. Baryogenesis is also discussed.

Optical Very Long Baseline Interferometry (VLBI) offers the potential for unprecedented angular resolution in both astronomical imaging and precision measurements. Classical approaches, however, face significant limitations due to photon loss, background noise, and the requirements for dynamical delay lines over large distances. This document surveys recent developments in quantum-enabled VLBI, which aim to address these challenges using entanglement-assisted protocols, quantum memory storage, and nonlocal measurement techniques. While its application to astronomy is well known, we also examine how these techniques may be extended to geodesy -- specifically, the monitoring of Earth's rotation. Particular attention is given to quantum-enhanced telescope architectures, including repeater-based long baseline interferometry and quantum error-corrected encoding schemes, which offer a pathway toward high-fidelity optical VLBI. To aid the discussion, we also compare specifications for key enabling technologies to current state-of-the-art experimental components, including switching rates, gate times, entanglement distribution rates, and memory lifetimes. By integrating quantum technologies, future interferometric networks may achieve diffraction-limited imaging at optical and near-infrared wavelengths, surpassing the constraints of classical techniques and enabling new precision tests of astrophysical and fundamental physics phenomena.

Sparse observations and coarse-resolution climate models limit effective regional decision-making, underscoring the need for robust downscaling. However, existing AI methods struggle with generalization across variables and geographies and are constrained by the quadratic complexity of Vision Transformer (ViT) self-attention. We introduce ORBIT-2, a scalable foundation model for global, hyper-resolution climate downscaling. ORBIT-2 incorporates two key innovations: (1) Residual Slim ViT (Reslim), a lightweight architecture with residual learning and Bayesian regularization for efficient, robust prediction; and (2) TILES, a tile-wise sequence scaling algorithm that reduces self-attention complexity from quadratic to linear, enabling long-sequence processing and massive parallelism. ORBIT-2 scales to 10 billion parameters across 32,768 GPUs, achieving up to 1.8 ExaFLOPS sustained throughput and 92-98% strong scaling efficiency. It supports downscaling to 0.9 km global resolution and processes sequences up to 4.2 billion tokens. On 7 km resolution benchmarks, ORBIT-2 achieves high accuracy with R^2 scores in the range of 0.98 to 0.99 against observation data.

We consider initial conditions leading to Starobinsky inflation in the general quadratic gravity, where the action of the theory contains one more curvature square invariant in addition to $R^2$. We have chosen corresponding coefficients in a way so that the inflationary solution keeps to be stable. Our numerical results show that despite the configuration of initial conditions in the $(H,R)$ plane, leading to Starobinsky inflation can change considerably from the $R^2$ theory, realization of inflation does not need a fine-tuning of the initial conditions.

The equation of state (EoS) for neutron stars is a crucial topic in astrophysics, nuclear physics, and quantum chromodynamics (QCD), influencing their structure, stability, and observable properties. This review classifies EoS models into hadronic matter, hybrid, and quark matter models, analyzing their assumptions, predictions, and constraints. While hadronic models characterize nucleonic matter, potentially including contributions from hyperons or mesons, hybrid models introduce phase transitions to quark matter, and quark models hypothesize the presence of deconfined quark matter cores or entirely quark-composed stars. By synthesizing results from recent theoretical and observational studies, this review aims to offer a comprehensive understanding of the methodologies used in constructing neutron star EoS, their implications, and future directions.

The exterior spacetime geometry surrounding an uncharged, spinning black hole in general relativity depends only upon its mass and spin. However, the exterior geometry surrounding any other rotating compact object, for example a neutron star, will generally depend upon higher moments in its multipole expansion, which will in turn be dependent upon the object's equation of state. Using general relativistic hydrodynamics and electrodynamics simulations, we illustrate that the presence or absence of these higher moments (assuming a physically realistic neutron star equation of state) has a significant qualitative effect near the surface of the compact object on the dynamics of unmagnetized accretion, and a smaller quantitative effect on the electromagnetic field configuration of its magnetosphere. In some places, the discrepancies in energy-momentum density are found to reach or exceed 50%, with electric field strength discrepancies in excess of 10%. We argue that many of these differences are likely to be amplified by the inclusion of more sophisticated plasma physics models, and are therefore likely to be relevant for the dynamics of gravitational collapse, and potentially also for particle acceleration and jet launching. These discrepancies suggest important limitations regarding the use of the Kerr metric when performing numerical simulations around neutron stars.

Quasi-universal relations are known to exist among various neutron star observables that do not depend sensitively on the underlying nuclear matter equations of state. For example, some of these relations imply that the tidally induced multipole moments are approximately characterized by the electric-type quadrupolar tidal deformability. Such relations can be used to reduce the number of independent tidal parameters in gravitational-waveform modeling, thereby allowing us to infer extreme nuclear matter properties more accurately and test General Relativity in an insensitive manner to uncertainties in nuclear physics. We present a comprehensive theoretical investigation into approximate universality of neutron stars. Our approach employs a semi-analytic relativistic stellar interior model, which extends the Tolman VII solution, thereby enabling a refined exploration of the tidal properties of nonrotating stars within a semi-analytic framework. The derived power-law relations among various tidal deformabilities -- referred to as the universal Love relations -- agree well with expressions in previous work found empirically. We elucidate how the equation-of-state dependence is suppressed in a particular combination of macroscopic physical parameters induced by perturbations and demonstrate that the relation between the moment of inertia and electric-type quadrupolar tidal deformability (I-Love relation) rests on the same underlying mechanism. Our findings indicate that the approximate universality of neutron stars can be attributed to low compressibility, consistent with some of the previous studies on the possible origin of the universality.

We study the cosmology of multi-field Dark Energy, using a well-motivated axio-dilaton model that contains the minimal number of fields to have the 2-derivative sigma-model interactions that power-counting arguments show naturally compete with General Relativity at low energies. Our analysis differs from earlier, related, studies by treating the case where the dilaton's couplings to matter are large enough to require screening to avoid unacceptable dilaton-mediated forces in the solar system. We use a recently proposed screening mechanism that exploits the interplay between stronger-than-gravitational axion-matter couplings with the 2-derivative axion-dilaton interactions to suppress the couplings of the dilaton to bulk matter. The required axion-matter couplings also modify cosmology, with the axion's background energy density turning out to resemble early dark energy. We compute the properties of the axion fluid describing the rapid oscillations of the axion field around the time-dependent minimum of its matter-dependent effective potential, extending the usual formalism to include nontrivial kinetic sigma-model interactions. We explore the implications of these models for the Cosmic Microwave Background and the growth of structure and find that for dilaton potentials of the Albrecht-Skordis form (itself well-motivated by UV physics), successful screening can be consistent with the early dark energy temporarily comprising as much as 10% of the total density in the past. We find that increasing the dilaton-matter coupling decreases the growth of structure due to enhanced Hubble friction, an effect that dominates the usual fifth-force effects that amplify structure growth.

We review the role of primordial black holes (PBHs) for illuminating the dark ages of the cosmological evolution and as dark matter candidates. We elucidate the role of phase transitions for primordial black hole formation in the early Universe and focus our attention to the cosmological QCD phase transition within a recent microscopical model. We explore the impact of physics beyond the Standard Model on the cosmic equation of state and the probability distribution for the formation of primordial black holes which serve as dark matter (DM) candidates. We argue that besides primordial black holes also droplet-like quark-gluon plasma inhomogeneities may become gravitationally stabilized for a sufficiently long epoch to distill baryon number and form nuclear matter droplets which upon their evaporation may enrich the cosmos locally with heavy $r$-process elements already in the early Universe.