Vote on papers for Monday, Jul 21 2025

The thermal Sunyaev$-$Zel'dovich (SZ) effect offers a unique probe of the hot and diffuse universe that could help close the missing baryon problem. Traditional extractions of the SZ effect, however, exhibit systematic noise that may lead to unreliable results. In this work, we provide an alternative solution using a three-dimensional Attention Nested U-Net trained end-to-end with supervised learning. Our labeled data consists of simulated SZ signals injected into $\textit{Planck}$ frequency maps, allowing our model to learn how to extract SZ signals in the presence of realistic noise. We implement a curriculum learning scheme that gradually exposed the model to weaker SZ signals. The absence/presence of curriculum learning significantly impacted the amount of bias and variance present in the reconstructed SZ signal. The results from our method were comparable to those from the popular $\textit{needlet internal linear combination}$ (NILC) method when evaluated on simulated data as well as real-world SZ signals. We conclude by discussing future avenues for advancing machine learning extractions of SZ signals.

A crucial ingredient affecting fast neutrino flavor conversion in core-collapse supernovae (SNe) is the shape of the angular distribution of the electron-neutrino lepton number (ELN). The presence of an ELN crossing signals favorable conditions for flavor conversion. However, the dependence of ELN crossings on the SN properties is only partially understood. We investigate a suite of 12 spherically symmetric neutrino-hydrodynamics simulations of the core collapse of a SN with a mass of $18.6 M_\odot$; each model employs different microphysics (i.e., three different nuclear equations of state, with and without muon creation) and includes or not a mixing-length treatment for proto-neutron star convection. We solve the Boltzmann equations to compute the neutrino angular distributions relying on static fluid properties extracted from each of the SN simulations in our suite for six selected post-bounce times. We explore the dependence of the ELN distributions on the SN microphysics and proto-neutron star convection. We find that the latter shifts the proto-neutron star radius outwards, favoring the appearance of ELN crossings at larger radii. On the other hand, muon creation causes proto-neutron star contraction, facilitating the occurrence of ELN crossings at smaller radii. These effects mildly depend on the nuclear equation of state. Our findings highlight the subtle impact of the SN microphysics, proto-neutron star convection, and neutrino transport on the ELN angular distributions.

We analyzed the deepest Hubble Space Telescope (HST) F275W ultraviolet (UV) imaging of M87 to obtain the most robust constraints on its star formation rate (SFR) and star formation history (SFH). After removing the galaxy continuum and globular clusters, we detected an excess of UV point sources near the center. By comparing their colors to young stellar source (YSS) colors generated by stochastically simulated star formation (SF) for various SFRs and SFHs, we ruled out their origin as a UV-upturn population and identified them as YSS. We found an extremely low SFR of $\sim 2\times10^{-5}$ M$_\odot$ yr$^{-1}$ in M87, with evidence of a weak starburst $\sim$125 Myr ago that formed $\sim 1000$ M$_\odot$ of stars. Unlike other cool-core clusters where SF is stronger and directly linked to cooling gas, we found no spatial correlation between YSS and H$\alpha$ filaments. Comparing SF activity with M87's AGN outburst history suggests that recent AGN feedback events ($\lesssim$12 Myr ago) neither triggered nor were associated with any detectable SF, however, earlier outbursts may have triggered weak starbursts. We detected UV filaments co-spatial with H$\alpha$ filaments with similar lengths and widths, though they are obscured by dust near the center. These filaments are likely powered by metal-line emission from collisional ionization, suggesting ongoing low-level precipitation of the intracluster medium. Our results indicate that AGN feedback has quenched SF significantly in M87 for at least 200 Myr, even though some precipitation persists. Additionally, we identified a hotspot created by the counterjet, with the spectral index also constrained.

We present the first independent re-analysis of the galaxy clustering data from DESI Data Release 1, utilizing an effective field theory (EFT)-based full-shape model. We analyze the power spectra and bispectra of the public catalogs using a custom-built pipeline based on window-deconvolved quasi-optimal estimators, carefully accounting for all relevant systematic effects. Compared to the official collaboration analysis, we add the galaxy power spectrum hexadecapole and the bispectrum monopole, and also introduce a novel stochastic estimator for fiber collisions, which facilitates robust bispectrum analyses. As a first application, we perform an EFT-based full-shape analysis of the DESI power spectra and bispectra in the context of the standard cosmological model, $\Lambda$CDM. Using external priors on the physical baryon density and the primordial power spectrum tilt, we constrain the matter density fraction to $\Omega_m=0.284\pm 0.011$, the Hubble constant to $h=0.707\pm 0.011$, and the mass fluctuation amplitude to $\sigma_8=0.811\pm 0.030$. The bispectrum has a noticeable effect on parameter estimation: it sharpens the constraints on $\sigma_8$ and $\Omega_m$ by $\approx 10$\% and shifts $\Omega_m$ by $\approx 1\sigma$ towards the Planck $\Lambda$CDM value. Combining our full-shape likelihood with the official DESI DR2 BAO measurements, cosmological parameters shift further towards the \textit{Planck} values, with $\Omega_m=0.296\pm 0.007$, $h=0.688\pm 0.006$, $\sigma_8=0.818\pm 0.029$ (with tighter constraints obtained in joint analyses). Finally, the galaxy bispectrum data dramatically improves measurements of quadratic bias parameters, which are consistent with predictions from halo occupation distribution models. Our work highlights the importance of higher-order statistics and sets the stage for upcoming full-shape analyses of non-minimal cosmological models.

In this study, we explore the properties of diffuse intergalactic gas residing in the outskirts of the four nearby, low-mass galaxy groups NGC1052, NGC5866, NGC4631, and NGC3992 (all at cz<2000 km/s) beyond their group virial radii. Using archival ultraviolet absorption spectra of bright AGN observed with HST/COS, we search for HI Ly a absorption near the groups' recession velocities along 35 sightlines that pass the outer group medium (OGrM) at normalized impact parameters rho/Rvir=1-3. We derive HI column densities of the absorbers and constrain the physical conditions in the gas by using a hydrostatic toy model of the groups' gas environment and assuming photoionization. HI Ly absorption near the groups' recession velocities is detected along 19 sightlines with HI column densities in the range log N(HI)=12.50-14.34, implying a high OGrM detection rate of more than 50 percent. We transform this value into an incidence rate of OGrM absorbers per unit redshift of dN/dz=this http URL is 25 percent above the value derived for the general population of Ly a absorbers within z = 0 filaments and more than twice the value for the z = 0 Ly a forest. From the modeling, we obtain lower limits for the gas densities from log n_H=-5.00 to -3.72, comparable to densities found in the overall Ly a forest. Our study unveils a large cross section and overdensity of Ly a absorbers in the outskirts of these four nearby groups. Such an overdensity is in line with a previously proposed scenario, in which AGN feedback lifts gaseous material to large distances beyond the virial radius of groups into the OGrM. However, a larger survey of OGrM absorbers and a comparison with hydrodynamical simulations will be necessary to constrain the cosmological mass density of OGrM absorbers and pinpoint their role in cosmological structure formation and galaxy/group evolution.

Sub-Neptune exoplanets are the most abundant type of planet known today. As they do not have a Solar System counterpart, many open questions exist about their composition and formation. Previous spectroscopic studies rule out aerosol-free hydrogen-helium-dominated atmospheres for many characterized sub-Neptunes but are inconclusive about their exact atmospheric compositions. Here we characterize the hot (Teq=1311K) sub-Neptune HD 86226 c, which orbits its G-type host star. Its high equilibrium temperature prohibits methane-based haze formation, increasing the chances for a clear atmosphere on this planet. We use HST data taken with WFC3 and STIS from the Sub-neptune Planetary Atmosphere Characterization Experiment (SPACE) Program to perform near-infrared 1.1-1.7micrometer transmission spectroscopy and UV characterization of the host star. We report a featureless transmission spectrum that is consistent within 0.4 sigma with a constant transit depth of 418+-14ppm. The amplitude of this spectrum is only 0.01 scale heights for a H/He-dominated atmosphere, excluding a cloud-free solar-metallicity atmosphere on HD 86226 c with a confidence of 6.5 sigma. Based on an atmospheric retrieval analysis and forward models of cloud and haze formation, we find that the featureless spectrum could be due to a metal enrichment [M/H] above 2.3 (3 sigma confidence lower limit) of a cloudless atmosphere, or silicate (MgSiO3), iron (Fe), or manganese sulfide (MnS) clouds. For these species, we perform an investigation of cloud formation in high-metallicity, high-temperature atmospheres. Our results highlight that HD 86226c does not follow the aerosol trend of sub-Neptunes found by previous studies. Follow-up observations with the JWST could determine whether this planet aligns with the recent detections of metal-enriched atmospheres or if it harbors a cloud species otherwise atypical for sub-Neptunes.

Numerical simulations have become an indispensable tool in astrophysics. To interpret their results, it is critical to understand their intrinsic variability, that is, how much the results change with numerical noise or inherent stochasticity of the physics model. We present a set of seven realisations of high-resolution cosmological zoom-in simulations of a Milky Way-like galaxy with the Auriga galaxy formation model. All realisations share the same initial conditions and code parameters, but draw different random numbers for the inherently stochastic parts of the model. We show that global galaxy properties at $z=0$, including stellar mass, star formation history, masses of stellar bulge and stellar disc, the radius and height of the stellar disk change by less than $10\%$ between the different realisations, and that magnetic field structures in the disc and the halo are very similar. In contrast, the star formation rate today can vary by a factor of two and the internal morphological structure of the stellar disc can change. The time and orbit of satellite galaxies and their galaxy properties when falling into the main halo are again very similar, but their orbits start to deviate after first pericenter passage. Finally, we show that changing the mass resolution of all matter components in the Auriga model changes galaxy properties significantly more than the intrinsic variability of the model, and that these changes are systematic. This limits detailed comparisons between simulations at different numerical resolutions.

Most general relativistic magnetohydrodynamic simulations of black hole (BH) hot accretion flows are initialized with small rotating tori and produce stable jets with only small fluctuations. However, recent studies using larger scale Bondi-like initial conditions have reported intermittent jet activity and loss of coherent rotation. To investigate the differences, we modify the standard torus setup across four BH spins: $a_*=0$, $0.5$, $0.9$, $-0.9$. First, we increase the torus size significantly (pressure maximum at 500 gravitational radii), allowing long simulations ($2.8\times10^5$ gravitational times) without gas depletion. These runs reproduce the weak variability seen in smaller tori, indicating that a larger dynamic range alone does not cause strong fluctuations. We observe moderate suppression of the accretion rate by factors of $\sim 1.6, ~2.5$ for BH spins $a_*=0.5,~0.9$, respectively, compared to $a_*=0$. Also, the density profile scales as $\rho(r)\propto r^{-1.1}$ for prograde BHs. Next, we considerably strengthen the initial magnetic field in the large torus by setting the plasma-$\beta\approx 1$. This induces strong variability in the evolution. The jet efficiency in the $a_*=0.9$ model, for instance, now varies by over 3 orders of magnitude, and gas rotation reverses directions. Combining these results with prior studies, we propose that a key parameter is the ratio $R$ between the rotational and magnetic energies in the initial state. Strong variability appears later in models with a larger value of $R$. The implication is that all simulations, and by extension all hot accretion flows in Nature, will ultimately develop intermittent jets if evolved long enough.

Recent studies have focused on the low-$z$ Ly$\alpha$ forest as a potential constraint on galactic feedback, as different AGN and stellar feedback models in hydrodynamic simulations produce varying intergalactic medium (IGM) statistics. However, existing low-$z$ Ly$\alpha$ forest data provide insufficient observational constraints for simulations due to their low precision and the lack of observations from $z \sim 0.4$ to $1.8$. The Habitable Worlds Observatory, equipped with an ultraviolet (FUV / NUV) spectrograph, could provide transformative data for this science by increasing both absorbers in the redshift range where current data are lacking and increasing the precision of the Ly$\alpha$ forest observational data at $z \leq 0.4$. This spectrograph should cover the wavelengths $1215-3402$ Å~and have a spectral resolution of $R =$~40,000. Distinguishing between certain active-galactic nuclei feedback models requires high precision ($\sim 5-10$\%) measurements of the Ly$\alpha$ forest 1D transmitted flux power spectrum for $z\leq 1.8$. This requires $\sim 830$ QSO spectra with S/N~$\geq5$ - $25$ depending on redshift. This data would also enable a comprehensive study on the thermal state of the IGM across time, and address the tension between the observed and simulated Ly$\alpha$ forest $b$-value distribution.

In modern-day astronomy, near-infrared, optical, and ultraviolet spectroscopy are indispensable for studying a wide range of phenomena, from measuring black hole masses to analyzing chemical abundances in stellar atmospheres. However, spectroscopic data reduction is often performed using instrument-specific pipelines or legacy software well-established and robust within the community that are often challenging to implement and script in modern astrophysical workflows. In this work, we introduce easyspec, a new Python package designed for long-slit spectroscopy, capable of reducing, extracting, and analyzing spectra from a wide range of instruments--provided they deliver raw FITS files, the standard format for most optical telescopes worldwide. This package is built upon the well-established long-slit spectroscopy routines of the Image Reduction and Analysis Facility (IRAF), integrating modern coding techniques and advanced fitting algorithms based on Markov Chain Monte Carlo (MCMC) simulations. We present a user-friendly open-source Python package that can be easily incorporated into customized pipelines for more complex analyses. To validate its capabilities, we apply easyspec to the active galactic nucleus G4Jy 1709, observed with the DOLORES spectrograph at the Telescopio Nazionale Galileo, measuring its redshift and estimating its supermassive black hole mass. Finally, we compare our results with a previous IRAF-based study.

In spite of their long detection history, the origin of hot Jupiters remains to be resolved. While multiple dynamical evidence suggests high-eccentricity migration is most likely, conflicts remain when considering hot Jupiters as a population in the context of warm and cold Jupiters. Here, we turn to atmospheric signatures as an alternative mean to test the origin theory of hot Jupiters, focusing on population level trends that arise from post-formation pollution, motivated by the upcoming Ariel space mission whose goal is to deliver a uniform sample of exoplanet atmospheric constraints. We experiment with post-formation pollution by planetesimal accretion, pebble accretion, and disk-induced migration and find that an observable signature of post-formation pollution is only possible under pebble accretion in metal-heavy disks. If most hot Jupiters arrive at their present orbit by high-eccentricity migration while warm Jupiters emerge largely in situ, we expect the atmospheric water abundance of hot Jupiters to be significantly elevated compared to warm Jupiters, at a level that is easily detectable by Ariel transit survey. We provide suggestions for future comparative atmospheric characterization between hot Jupiters and wide-orbit directly imaged planets to elucidate the properties of the dust substructures in protoplanetary disks.

The advent of JWST has revolutionized the study of faint satellite galaxies at $z \gtrsim 1$, enabling statistical constraints on galaxy evolution and the galaxy$-$halo connection in a previously unexplored mass and redshift regime. We compare satellite abundances at $1 < z < 3.5$ from recent JWST observations with predictions from cosmological dark matter-only zoom-in simulations. We identify and quantify several sources of biases that can impact theoretical satellite counts, finding that assumptions about subhalo tidal evolution introduce the largest uncertainty in predictions for the satellite mass function. Using a flexible galaxy disruption model, we explore a range of disruption scenarios, spanning hydrodynamically motivated and idealized prescriptions, to bracket plausible physical outcomes. We show that varying galaxy durability can change the predicted satellite mass functions by a factor of $\sim3.5$. The JWST data and our fiducial model are consistent within $1-2\sigma$ across the full redshift ($1 < z < 3.5$) and stellar mass ($M_\star> 10^7~\mathrm{M}_\odot$) range probed. We find evidence that subhalos are at least as long-lived as predicted by hydrodynamic simulations. Our framework will enable robust constraints on the tidal evolution of subhalos with future observations. This work presents the first direct comparison between cosmological models and observations of the high-redshift satellite population in this low-mass regime. These results showcase JWST's emerging power to test structure formation in the first half of the Universe in a new domain and to constrain the physical processes driving the evolution of low-mass galaxies across cosmic time.

Discovery rates of supernovae are expected to surpass one million events annually with the Vera C. Rubin Observatory. With unprecedented sample sizes of both common and rare transient types, photometric classification alone will be insufficient for finding one-in-a-million events and prioritizing the 1% of events for spectroscopic follow-up observations. Here, we present reLAISS, a modified framework for similarity searches of supernovae using extracted features of ZTF light curves and Pan-STARRS host galaxy photometry and built on the original LAISS framework. Unlike its predecessor, reLAISS couples interpretable light curve morphology features with extinction-corrected host-galaxy colors to probe both explosion physics and associated stellar populations simultaneously. The library allows users to customize the number of neighbors retrieved, the weight of host and light curve features, and the use of Monte Carlo simulations to ensure relevant matches when features are poorly constrained. We release reLAISS as a pip-installable package with an accompanying reference set of 20,000 features, and a set of tutorials that demonstrate the code's expanded functionality. All source code can be found at this https URL .

We measure the distribution of matter contained within the emptiest regions of the Universe: cosmic voids. We use the large overlap between the Ultraviolet Near-Infrared Optical Northern Survey (UNIONS) and voids identified in the LOWZ and CMASS catalogues of the Baryon Oscillation Spectroscopic Survey (BOSS) to constrain the excess surface mass density of voids using weak lensing. We present and validate a novel method for computing the Gaussian component of the conventional weak lensing covariance, adapted for use with void studies. We detect the stacked weak lensing void density profile at the $6.2\sigma$ level, the most significant detection of void lensing from spectroscopically-identified voids to date. We find that large and small voids have different matter density profiles, as expected from numerical studies of void profiles. This difference is significant at the $2.3\sigma$ level. Comparing the void profile to a measurement of the void-galaxy cross-correlation to test the linearity of the relationship between mass and light, we find good visual agreement between the two, and a galaxy bias factor of $2.45\pm0.36$, consistent with other works. This work represents a promising detection of the lensing effect from underdensities, with the goal of promoting its development into a competitive cosmological probe.

Context. A comprehensive study of comets over a wide heliocentric distance range helps us understand the physical processes driving their activity and reveals compositional differences across dynamical groups. C/2017 K2 (PANSTARRS) is a Dynamically New Oort Cloud comet (DNC) that showed activity as far as 23.75 au and displayed a CO-rich coma at 6.72 au, making it a key object to investigate pre- and post-perihelion behavior. Aims. We aim to study the long-term activity evolution and chemical composition of C/2017 K2 using photometry and spectroscopy, from October 2017 (r$_h$ = 15.18 au) pre-perihelion to April 2025 (r$_h$ = 8.46 au) post-perihelion. Methods. Broad-band and narrow-band imaging from both TRAPPIST telescopes enabled us to produce an 8-year light curve, color analysis, and derivation of activity slopes. Production rates of OH, NH, CN, C$_3$, and C$_2$ were computed using a Haser model, along with the dust proxy A(0)f$\rho$. High-resolution spectra from CRIRES$^+$ and UVES at three epochs (May - September 2022) provided simultaneous observations of parent and daughter species as the comet crossed the water sublimation zone. Results. The light curve of C/2017 K2 shows a complex evolution with varying slopes and a brightness plateau around perihelion, indicating multiple active species. Coma colors remained constant, suggesting uniform dust properties and similarity to other active long-period comets. Gas production rates indicate a typical C$_2$/CN composition with a high dust-to-gas ratio. Analysis of forbidden oxygen lines shows a transition from CO and CO$_2$-driven activity to water-driven sublimation inside 3 au. Infrared spectra reveal C/2017 K2 as a typical-to-enriched comet, with HCN identified as the main parent of CN, and C$_2$ likely originating from C$_2$H$_2$ rather than C$_2$H$_6$.

Long-period radio transients (LPTs) are an emerging group of radio transients that show periodic polarised radio bursts with periods varying from a few minutes to a few hours. Fewer than a dozen LPTs have been detected so far, and their origin (source and emission mechanism) remains unclear. Here, we report the discovery of a 1.5 hr LPT, ASKAP J144834-685644, adding to the current sample of sources. ASKAP J144834-685644 is one of the very few LPTs that has been detected from X-rays to radio. It shows a steep radio spectrum and polarised radio bursts, which resemble the radio emission in known LPTs. In addition, it also shows highly structured and periodic narrow-band radio emission. Multi-wavelength properties suggest that the spectral energy distribution (SED) peaks at near ultraviolet wavelengths, indicating the presence of a hot magnetic source. Combining multi-wavelength information, we infer that ASKAP J144834-685644 may be a near edge-on magnetic white dwarf binary (MWD), although we can not fully rule out ASKAP J144834-685644 being an isolated white dwarf pulsar or even a transitional millisecond pulsar (despite the lack of radio pulsations). If ASKAP J144834-685644 is an MWD binary, the observed broadband spectral energy distribution can be explained by emission from an accretion disk. This hints that some fraction of optically bright LPTs may be accreting binaries with the radio period being the orbital period. It might further suggest a connection between optically bright non-accreting synchronized LPTs, such as polars, and non-accreting asynchronous WD pulsars, such as AR Sco and J1912-4410.

Extreme emission line galaxies (EELGs) at high redshifts are considered key contributors to cosmic reionization at $z>6$ due to their higher ionization efficiencies. We have identified 119 H$\beta$ + [OIII] emitters at $z\sim7$ selected by a flux excess in the medium-band filter F410M in the public James Webb Space Telescope Cycle-1 fields. Our emitters exhibit a wide range in rest-frame H$\beta$ + [OIII] equivalent width (EWs), 420 $<$ EW$_{0}$/Å $<$ 6850 (with the median value of $\sim1700$ Å). Among them, 19 are EW$_{0}$ $>$ 3000 Å, which represent extreme populations even in the context of recent findings with JWST. They are characterized by (i) low stellar mass ($\sim 3\times10^{7}$ $\mathrm{M_{\odot}}$), (ii) blue colors ($\beta_{\rm UV}\sim -2.2$), and (iii) low dust attenuation ($A_{\mathrm{V}}\sim 0.1$ mag). We discuss the physical mechanisms responsible for the observed high rest-frame H$\beta$ + [OIII] EWs, including (1) photoionization by AGN, (2) stellar photoionization in the vicinity of HII regions, and (3) radiative shocks powered by outflows either from AGN or massive stars. Notably, we find 13 emitters with spatially offset H$\beta$ + [OIII] emission compared to the UV and stellar components. Given the absence of obvious signatures of actively accreting black holes, these emitters are likely under strong feedback-driven winds from massive stars. Lastly, we report a unique overdensity of EELGs in one of the observed fields. The discovery of such a "star-bursting" overdensity supports the idea that large ionizing bubbles formed around some EEGLs in the early Universe.

The Sun continuously expels a fraction of its own mass in the form of a steadily accelerating outflow of ionized gas called the "solar wind." The solar wind is the extension of the Sun's hot (million-degree Kelvin) outer atmosphere that is visible during solar eclipses as the bright and wispy corona. In 1958, Eugene Parker theorized that a hot corona could not exist for very long without beginning to accelerate some of its gas into interplanetary space. After more than half a century, Parker's idea of a gas-pressure-driven solar wind still is largely accepted, although many questions remain unanswered. Specifically, the physical processes that heat the corona have not yet been identified conclusively, and the importance of additional wind acceleration mechanisms continue to be investigated. Variability in the solar wind also gives rise to a number of practical "space weather" effects on human life and technology, and there is still a need for more accurate forecasting. Fortunately, recent improvements in both observations (with telescopes and via direct sampling by space probes) and theory (with the help of ever more sophisticated computers) are leading to new generations of predictive and self-consistent simulations. Attempts to model the origin of the solar wind are also leading to new insights into long-standing mysteries about turbulent flows, magnetic reconnection, and kinetic wave-particle resonances.

We characterize the post-inflationary evolution of the equation of state of the universe from the end of inflation until the onset of radiation domination, when the inflaton is coupled to a daughter field through a trilinear interaction. We consider an inflaton potential that is quadratic near the minimum and flattens in the inflationary regime. By simulating the dynamics in 2+1-dimensional lattices, we have tracked the long-term evolution of the equation of state for about ten e-folds of expansion, for various coupling strengths. The trilinear interaction initially excites daughter field modes through a process of tachyonic resonance immediately after inflation and triggers a temporary deviation of the equation of state from $\bar{w} = 0$ to a maximum value $\bar{w} = \bar{w}_{\rm max} < 1/3$. However, at much later times, the inflaton homogeneous mode once again dominates the energy density, pushing the equation of state towards $\bar{w} = 0$ until the onset of perturbative reheating. By combining the lattice results with a Boltzmann approach, we characterize the entire post-inflationary expansion history, which allows to calculate precise predictions for the inflationary CMB observables. We also accurately compute the redshift of the stochastic gravitational wave background produced during preheating, and show that taking the temporary return of the equation of state towards $\bar{w} = 0$ into account can reduce the amplitude by many orders of magnitude relative to previous estimates.

We present new time-series CCD \emph{VR} photometry of the globular cluster NGC 5634. We aim to use the known RR Lyrae stars, members of the cluster, as indicators of mean metallicity and distance. Accurate coordinates, periods and an identification chart of the variables in the field of our images are provided. A membership analysis was performed, based on $Gaia$-DR3 proper motions, for 3525 point sources within 15 arcmin from the cluster center. The membership status for each known variable was established and it was found that V10, V11 and V16 are most likely field stars. The variability of the RRab star V7, considered non-variable for a number of years, is demonstrated. Via the Fourier decomposition of the light curves of cluster member RR Lyrae, the mean metallicity and distance were calculated independently from RRab and RRc stars to find [Fe/H]$_{\rm ZW}= -1.67 \pm 0.11$ y $D=23.9 \pm 0.8$ kpc, and [Fe/H]$_{\rm ZW}= -1.69 \pm 0.22$ y $D=22.9 \pm 1.0$ kpc respectively.

Waves are thought to play a significant role in the heating of the solar atmosphere and the acceleration of the wind. Among the many types of waves observed in the Sun, the so-called p-modes with a 3 mHz frequency peak dominate the lower atmosphere. In the presence of magnetic fields, these waves can be converted into magnetohydrodynamic modes, which then leak into the corona through magnetic conduits. High-resolution off-limb observations have indeed revealed signatures of ubiquitous and global 3 mHz oscillations in the corona, although limited to low heights and to incompressible modes. We present high-cadence, high-resolution observations of the corona in the range 1.7 - 3.6 solar radii taken in broad-band 580-640 nm visible light by the Metis coronagraph aboard Solar Orbiter. These observations were designed to investigate density fluctuations in the middle corona. The data were acquired over several days in March 2022, October 2022, and for two days in April 2023. We selected representative regions of the corona on three sample dates. Analysis of the data in those regions revealed the presence of periodic density fluctuations. By examining several time-distance diagrams, we determined the main properties (apparent propagation speed, amplitude) of those fluctuations. We also show power spectra in selected locations in order to determine the dominant frequencies. We found wave-like, compressible fluctuations of low amplitude - of the order of 0.1 % of the background - in several large-scale regions in the corona at least up to 2.5 solar radii. We also found that the apparent propagation speeds of these perturbations typically fall in the range 150 - 450 km/s. A power spectrum analysis of time series revealed an excess power in the range 2-7 mHz, often with peaks at 3 or 5 mHz, i.e. in a range consistent with p-mode frequencies of the lower solar atmosphere.

We report the discovery of a possible sub-population of stellar clusters that appear to follow retrograde orbits around the third largest galaxy in the Local Group, M33 (Triangulum). This spiral disk galaxy has apparently had a mostly quiescent existence, although recent discoveries, particularly of a pronounced warp in the gas and stellar disk, suggest that M33's relatively quiet past was interrupted at least once by a dynamical interaction with another galaxy. We suggest that this sub population provides evidence of accretion of one or more dwarf galaxies in M33's history. We estimate a lower limit for the accreted halo virial mass of $M_{\rm vir} \sim (7 \pm 3) \times 10^{10}$\msun, accounting for about 10\% of the virial mass in the halo of M33 that has an accretion origin. We propose one of these accretion events as the source of the observed warp in M33's disk.

JWST observations of quasars in the Epoch of Reionization have revealed that many lie in host galaxies that are severely undermassive relative to the supermassive black holes. It is unclear how these systems will evolve to the tight local relation between stellar mass and black hole mass. We search for companions around the z=7.08 quasar ULAS J1120+0641 using JWST/NIRCam narrow, medium, and wide-band photometry to identify [O III] emitters at the quasar redshift, and explore the potential for growth of the host galaxy through future mergers. We find 22 sources near the quasar's redshift across our two 4.4 arcmin$^2$ fields, indicating that environment of ULAS J1120+0641 is strongly overdense in z~7.1 galaxies relative to the field. We estimate the potential future mass budget of the quasar host galaxy by summing the current stellar and gas masses of the quasar host and surrounding galaxies, correcting for incompleteness and selection effects. With no further black hole growth, ULAS J1120+0641 is unlikely to reach a $M_{\mathrm{BH}}/M_*$ ratio less than ~2.5% at z=0, still much higher than typical for local galaxies. However, such systems -- a quiescent black hole in a low-luminosity galaxy -- may have escaped detection locally if they are sufficiently distant.

Simulation-Based Inference (SBI) offers a principled and flexible framework for conducting Bayesian inference in any situation where forward simulations are feasible. However, validating the accuracy and reliability of the inferred posteriors remains a persistent challenge. In this work, we point out a simple diagnostic approach rooted in ensemble learning methods to assess the internal consistency of SBI outputs that does not require access to the true posterior. By training multiple neural estimators under identical conditions and evaluating their pairwise Kullback-Leibler (KL) divergences, we define a consistency criterion that quantifies agreement across the ensemble. We highlight two core use cases for this framework: a) for generating a robust estimate of the systematic uncertainty in parameter reconstruction associated with the training procedure, and b) for detecting possible model misspecification when using trained estimators on real data. We also demonstrate the relationship between significant KL divergences and issues such as insufficient convergence due to, e.g., too low a simulation budget, or intrinsic variance in the training process. Overall, this ensemble-based diagnostic framework provides a lightweight, scalable, and model-agnostic tool for enhancing the trustworthiness of SBI in scientific applications.

We present early results from a high-resolution analysis ($\sim$100-200pc) of the CO(2-1)/CO(1-0) line ratio in twelve nearby galaxies. We use new ALMA CO(1-0) observations from the Fundamental CO(1-0) Transition Survey (FACTS), and re-imaged CO(2-1) data from PHANGS. We make empirical classifications based on the optical and molecular gas morphologies, which show clear systematic trends in the variation of $R_{21}$ as a function of galactic structure. The sample includes barred and unbarred, and flocculent galaxies. The barred spiral galaxies follow a general trend: $R_{21}$ is high in the center, low along the bar, increases at the bar ends, and then declines in the outer parts of the disk. The structure dependence suggests the importance of galactic dynamics on molecular gas evolution, and consequently on star formation, in galaxies. $R_{21}$ fluctuates in the spiral arms for both barred and unbarred galaxies. HII regions increase $R_{21}$ locally in their surrounding gas and are often associated with galactic structures. Together, $R_{21}$ varies systematically as a function of galactic structure, dynamics, and star formation activity.

I measure and collect timings of phase markers (like eclipse times) for the orbits of 25 X-ray binaries (XRBs) so as to calculate the steady evolutionary period change ($\dot{P}$). I combine these with my observed $\dot{P}$ measures from 52 cataclysmic variables (CVs). Further, I subtract out the contributions from gravitational radiation ($\dot{P}_{\rm GR}$) and mass transfer ($\dot{P}_{\rm mt}$), deriving the period change from the residual unknown angular momentum loss ($\dot{P}_{\rm AML}$=$\dot{P}$-$\dot{P}_{\rm GR}$-$\dot{P}_{\rm mt}$). I have $\dot{P}_{\rm AML}$ measures for 77 XRBs and CVs, with these being direct measures of the driver of binary evolution. The venerable Magnetic Braking Model (MBM) of binary evolution has its most fundamental predictions tested, with most systems having predictions wrong by over one order-of-magnitude. Other proposed mechanisms to explain the AML also fail, so we are left with no known mechanism that dominates the AML. An alternative path to the AML law is empirical, where my $\dot{P}_{\rm AML}$ measures are fitted to a power-law involving the fundamental binary properties. With this, the dominant AML law for systems with orbital periods ($P$) from 0.13--1.0 days is $\dot{P}_{\rm AML} = -1500\times10^{-12} P^{1.29} M_{\rm prim}^{2.75} M_{\rm comp}^{-1.00}\dot{M}^{0.43}_{-8}$, in appropriate units. Similar AML laws for binaries below the Period Gap and for binaries with $P$$>$1.0 day are derived. These three AML laws are of good accuracy and are the best representations of the actual evolution for all 77 XRBs and CVs of all classes, so the three taken together can be called `universal'.

We investigate the feasibility of detecting galactic orbit dark matter passing through Earth by measuring its gravitational coupling with a wire under tension. We do so by exploring the transverse and longitudinal waves induced on the wire to detect a massive particle passing within $\sim 1$ m of the wire. The particle's $r^{-2}$ interaction with the wire provides an initial momentum which develops into a propagating wave carrying a distinctive time dependent displacement. Most interestingly, we find that both transverse and longitudinal waves develop with unique profiles, allowing for a full, three dimensional reconstruction of the particle's trajectory and its mass over velocity ratio. We find that, at interaction distances of 0.1 to 100 mm with a 90 micron diameter copper beryllium wire, Planck scale dark matter with mass $\sim 10^{19}$ GeV/$c^2$ would create immeasurable displacements on the scale of $10^{-24}$ to $10^{-26}$ m. In order to create displacements detectable by modern, commercially available, displacement sensors on the nanometer scale we require dark matter with a particle mass greater than $4 \times 10^7$ kg ($\sim 2 \times 10^{32}$ GeV/$c^2$). This is outside the upper limit of the Planck scale by 13 orders of magnitude and would also have such a low particle flux that a detection event would be implausible. Finally, we perform a similar analysis for a charged wire and an elementary charged particle with their electrostatic interaction, finding that a sufficiently slow charged particle would produce a transverse displacement comparable to the sensitivity of currently available sensors.

Newtonian noise limits the low-frequency sensitivity of ground-based gravitational wave detectors. While seismometers and geophones are commonly employed to monitor ground motion for Newtonian noise cancellation, their limited spatial coverage and high deployment costs hinder scalability. In this study, we demonstrate that distributed acoustic sensing offers a viable and scalable alternative, providing performance comparable to that of conventional seismic instruments. Using data from acoustic sensing and colocated seismometers during both natural and controlled events, we observe a strong correlation, exceeding 0.8, between the two sensor types in the 3 to 20 Hz frequency band relevant for Newtonian noise. Moreover, when distributed acoustic sensing data are used to predict geophone signals, the correlation remains high, above 0.7, indicating that distributed acoustic sensing accurately captures both the spatial and spectral features of ground motion. As a case study, we apply distributed acoustic sensing data to cancel noise recorded by the vertical component of a seismometer and compare the results with those obtained using geophone data for the same task. Both distributed acoustic sensing and geophone-based cancellations yield a residual noise factor of 0.11 at 20 Hz. These findings confirm the feasibility of using distributed acoustic sensing for Newtonian noise mitigation and highlight its potential, in combination with traditional seismic sensors, to improve environmental monitoring and noise suppression in current and next-generation gravitational wave observatories.

We present `James Webb Space Telescope' observations of M4 -- the closest globular cluster -- that probe the lower Main Sequence down to the hydrogen-burning limit. The unveiled stellar sequence reaches much fainter luminosities than previously possible, revealing a few extremely red objects that are consistent with brown dwarfs as cool as T_eff~1000K. However, the lack of a second JWST epoch presently prevents us from verifying the cluster membership of these objects. By cross-matching our data with archival `Hubble Space Telescope' images, we are able to verify cluster membership for a subset of objects down to T_eff~3000K. The observed color distribution indicate that the lower Main Sequence of M4 is likely deficient in oxygen compared to its higher-mass post-Main Sequence members by ~0.5dex. This feature has now been observed in three different globular clusters (M4, NGC6397 and 47Tuc), suggesting a general trend. Finally, we derive the mass function of the Galactic bulge in the background of M4. The mass function was found to have the bottom-heavy slope of alpha=0.88+/-0.36 and appears to terminate at ~0.15 M_Sun, although the latter value may be overestimated due to the limited sample size.

V960~Mon is an FU Orionis object that shows strong evidence of a gravitationally unstable spiral arm that is fragmenting into several dust clumps. We report the discovery of a new substellar companion candidate around this young star, identified in high-contrast $L'$-band imaging with VLT/ERIS. The object is detected at a projected separation of $0.898 \pm 0.01$ arcseconds with a contrast of $(8.39 \pm 0.07) \times 10^{-3}$. The candidate lies close to the clumps previously detected in the sub-mm (at 1.3 mm) and is co-located with extended polarized IR signal from scattered stellar irradiation, suggesting it is deeply embedded. The object is undetected in the SPHERE $H$-band total intensity, placing an upper mass limit of $\sim38~M_\mathrm{Jup}$ from the contrast curve. Using evolutionary models at an assumed age of 1~Myr, we estimate a mass of $\sim660~M_\mathrm{Jup}$ from the L' brightness; however, this value likely includes a significant contribution from a disk around the companion. The discrepancy between near- and mid-infrared results again suggests the source is deeply embedded in dust. This candidate may represent an actively accreting, disk-bearing substellar object in a young, gravitationally unstable environment.

The eight-planet Kepler-90 system exhibits the greatest multiplicity of planets found to date. All eight planets are transiting and were discovered in photometry from the NASA Kepler primary mission. The two outermost planets, g ($P_g$ = 211 d) and h ($P_h$ = 332 d) exhibit significant transit-timing variations (TTVs), but were only observed 6 and 3 times respectively by Kepler. These TTVs allow for the determination of planetary masses through dynamical modeling of the pair's gravitational interactions, but the paucity of transits allows a broad range of solutions for the masses and orbital ephemerides. To determine accurate masses and orbital parameters for planets g and h, we combined 34 radial velocities (RVs) of Kepler-90, collected over a decade, with the Kepler transit data. We jointly modeled the transit times of the outer two planets and the RV time series, then used our two-planet model to predict their future times of transit. These predictions led us to recover a transit of Kepler-90 g with ground-based observatories in May 2024. We then combined the 2024 transit and several previously unpublished transit times of planets g and h with the Kepler photometry and RV data to update the masses and linear ephemerides of the planets, finding masses for g and h of $15.0 \pm 1.3\, M_\oplus$, and $203 \pm 16\, M_\oplus$ respectively from a Markov Chain Monte Carlo analysis. These results enable further insights into the architecturally rich Kepler-90 system and pave the way for atmospheric characterization with space-based facilities.

Gamma-ray bursts (GRBs) are traditionally classified into long (lGRBs) and short (sGRBs) durations based on their $T_{90}$, with lGRBs widely used as tracers of the cosmic star formation rate (SFR) due to their observed association with core-collapse supernovae. However, recent detections of kilonovae accompanying some lGRBs challenge this assumption, suggesting potential contamination from compact binary mergers. Here, we move beyond the conventional $T_{90}$-based classification and focus exclusively on GRBs directly associated with supernovae - the most direct signatures of massive stellar collapse - to reassess their connection to the SFR. Using a sample of SN/GRBs, we construct the luminosity - redshift ($L$-$z$) plane and uncover a significant correlation between these variables. To account for observational biases, we apply the $\tau$ statistic and Lynden-Bell's $C^{-}$ method to derive the intrinsic luminosity function and formation rate. Our analysis reveals that even among this well-defined subsample, the SN/GRB formation rate still exceeds the SFR at low redshifts ($z < 1$). These findings suggest that GRBs at low redshift may not serve as reliable tracers of the SFR, and that larger samples are required to further investigate this discrepancy.

The Fermi and eROSITA bubbles are two pairs of large diffuse multi-band structures detected in our Milky Way galaxy, yet their origins remain uncertain. Using three-dimensional hydrodynamic simulations, we investigate the model where these bubbles were produced by two successive active galactic nucleus jet activities from the Galactic center, with their edges delineated by the correspondingly driven forward shocks. Our simulations successfully reproduce the observed morphology, multi-band X-ray surface brightness distributions between 0.11 and 2.04 keV, and O viii/O vii line ratios, matching ROSAT and eROSITA observations. The results indicate that the first jet pair, launched 15 Myr ago, created the eROSITA bubbles, now extending to ~15 kpc with gas temperatures of ~0.25 keV. The second jet pair, launched 5 Myr ago, produced the Fermi bubbles, currently reaching ~10 kpc with temperatures of 0.3--0.4 keV. The sharp edges of the eROSITA and Fermi bubbles are naturally explained by the forward shocks in this double-episode jet scenario. Our findings support a Galactic-scale origin for these structures and suggest episodic jet activities from the Galactic center.

Observations of accreting black hole (BH) systems, such as microquasars and supermassive black holes, often reveal a precessing jet with changing directions, indicating a misaligned accretion flow relative to the BH spin. The precession is commonly attributed to the Lense-Thirring (LT) effect, which arises from the BH's rotation twisting the surrounding spacetime and accretion flow. In the strongly magnetized regime, which is preferred accretion flow conditions for M~87$^*$ and likely other jet-producing systems, the large-scale magnetic field can significantly influence the flow dynamics. Here, we perform large-scale three-dimensional general relativistic magnetohydrodynamic simulations of tilted accretion onto a rotating BH, and find a never-seen-before new retrograde precession. This precession arises from a magnetic torque on the disk generated by the poloidal magnetic field aligned with the BH's rotation, opposing the LT torque. This finding highlights the unique property of highly magnetized accretion flows around BHs and provides a new interpretation of jet precession observed in many systems.

Compact binary mergers, detected in gravitational waves since 2015, are candidate sources for astrophysical neutrinos in the GeV regime from proton-proton and proton-neutron collisions. This contribution presents the results of the search for such a signal using mergers detected during the fourth observing run of the LIGO, Virgo, and KAGRA interferometers. We use the dense infill array at the center of the IceCube detector, IceCube-DeepCore, to select neutrino candidates in the 0.5-5 GeV energy range. The search for a statistically significant excess associated with an astrophysical signal is performed in a $\pm$ 500 s window around the gravitational wave detection time. We do not observe any statistically significant excess in the neutrino data, and set upper limits on the neutrino emission from these objects. Additionally, we search for subpopulations of neutrino-emitting sources, including merger events detected in previous observing runs; no significant signal has been identified yet.

Context: Helioseismology aims to infer the properties of the solar interior by analyzing observations of acoustic oscillations. The interpretation of the helioseismic data is however complicated by the non-trivial relationship between helioseismic observables and the physical perturbations associated with acoustic modes, as well as by various instrumental effects. Aims: We aim to improve our understanding of the signature of acoustic modes measured in the Helioseismic and Magnetic Imager (HMI) continuum intensity and Doppler velocity observables by accounting for radiative transfer, solar background rotation, and spacecraft velocity. Methods: We start with a background model atmosphere that accurately reproduces solar limb darkening and the Fe I 6173Å spectral line profile. We employ first-order perturbation theory to model the effect of acoustic oscillations on inferred intensity and velocity. By solving the radiative transfer equation in the atmosphere, we synthesize the spectral line, convolve it with the six HMI spectral windows, and deduce continuum intensity (hmi.Ic_45s) and Doppler velocity (hmi.V_45s) according to the HMI algorithm. Results: We analytically derive the relationship between mode displacement in the atmosphere and the HMI observables, and show that both intensity and velocity deviate significantly from simple approximations. Specifically, the continuum intensity does not simply reflect the true continuum value, while the line-of-sight velocity does not correspond to a straightforward projection of the velocity at a fixed height in the atmosphere. Our results indicate that these deviations are substantial, with amplitudes of approximately 10% and phase shifts of around 10 degrees across the detector for both observables. Moreover, these effects are highly dependent on the acoustic mode under consideration and the position on the solar disk.

The Pierre Auger Observatory, located in La Pampa Amarilla, Argentina, has been continuously acquiring data since 2004. It comprises a surface detector array covering 3,000 km$^2$ and 27 fluorescence telescopes, designed to detect extensive air showers initiated by ultra-high-energy cosmic rays. An upgrade to the Observatory was commissioned in 2024, enhancing the existing water-Cherenkov detectors with additional radio antennas, surface scintillator detectors, and a buried scintillator array. This compilation of contributions to the 39th International Cosmic Ray Conference, held in Geneva, Switzerland (July 15-24, 2025), presents recent results from the Pierre Auger Collaboration, addressing a wide range of fundamental questions in astroparticle physics. The included papers cover measurements of the energy spectrum, mass composition, and arrival directions of ultra-high-energy cosmic rays, investigations of hadronic interactions in extensive air showers, and searches for ultra-high-energy photons and neutrinos. Additional topics include radio detection techniques, solar-related phenomena, and atmospheric events such as ELVES and TGFs. The list also contains first results and performance evaluations of the upgraded detectors, AugerPrime, along with reports on outreach and social engagement initiatives conducted by the Collaboration.

Recent observations with the LOw Frequency ARray (LOFAR) have revealed 19 nearby M dwarfs showing bright circularly polarised radio emission. One of the possible sources of such emission is through magnetic star-planet interactions (MSPI) with unseen close-in planets. We present initial results from a spectroscopic survey with the Habitable-zone Planet Finder (HPF) and NEID spectrographs designed to characterize this sample and further investigate the origin of the radio emission. We provide four new insights into the sample. I) We uniformly characterize the stellar properties, constraining their effective temperatures, surface gravities, metallicities, projected rotational velocities, rotation periods, stellar radii, and stellar inclinations where possible. Further, from a homogenous analysis of the HPF spectra, we infer their chromospheric activity and spectroscopic multiplicity states. From this, we identify GJ 625, GJ 1151, and LHS 2395 as single, quiescent stars amenable to precise RV follow-up, making them strong MSPI candidates. II) We show that the distribution of stellar inclinations are compatible with an isotropic distribution, providing no evidence for a preference to pole-on configurations. III) We refine the radial velocity solution for GJ 625 b, the only currently known close-in planet in the sample, reducing the uncertainty in its orbital period by a factor of three, to facilitate future phase-dependent radio analysis. IV) Finally, we identify GJ 3861 as a spectroscopic binary with an orbital period of $P=14.841181_{-0.00010}^{+0.00011}$ d, making it the only confirmed binary with a relatively short orbit in the sample, where we surmise the radio emission is likely related to magnetospheric interactions between the two stars. These results advance our understanding of radio-emitting M dwarfs and establish an observational foundation for identifying MSPI.

Some previous studies have suggested that massive and intermediate-mass primordial black holes (PBHs) could comprise a substantial fraction of dark matter in the Universe. Such black holes, if they existed in the Milky Way halo, would give rise to long-duration microlensing events that may potentially last for years. However, earlier searches were not sufficiently sensitive to detect such events. Here, we present the results of searches for long-timescale gravitational microlensing events toward the Small Magellanic Cloud (SMC) using nearly 20 years of photometric observations collected by the Optical Gravitational Lensing Experiment (OGLE) from 2001 to 2020. We found six events, three of which are new discoveries. We use a sample of five events to measure the microlensing optical depth toward the SMC $\tau = (0.32 \pm 0.18) \times 10^{-7}$ and the event rate $\Gamma = (1.18 \pm 0.57) \times 10^{-7}\,\mathrm{yr}^{-1}\,\mathrm{star}^{-1}$. The properties of the detected events are consistent with lenses originating from known stellar populations within the SMC or in the Milky Way disk. No events with timescales longer than 1 yr were detected, which provides competitive limits on the fraction of massive compact objects, including PBHs, in the Milky Way dark matter halo. Together with the earlier OGLE studies of microlensing events in the direction of the Large Magellanic Cloud, these observations rule out PBHs and other compact objects with masses ranging from $10^{-8}$ to $10^3\,M_{\odot}$ as dominant components of dark matter.

Gamma-ray burst (GRB) afterglows and kilonovae (KNe) are electromagnetic transients that can accompany binary neutron star (BNS) mergers. Therefore, studying their emission processes is of general interest for constraining cosmological parameters or the behavior of ultra-dense matter. One common method to analyze electromagnetic data from BNS mergers is to sample a Bayesian posterior over the parameters of a physical model for the transient. However, sampling the posterior is computationally costly, and because of the many likelihood evaluations needed in this process, detailed models are too expensive to be used directly in inference. In the present article, we address this problem by introducing fiesta, a python package to train machine learning (ML) surrogates for GRB afterglow and kilonova models that can accelerate likelihood evaluations. Specifically, we introduce extensive ML surrogates for the state-of-the-art GRB afterglow models afterglowpy and pyblastafterglow, as well as a new surrogate for the KN emission as modeled by the possis code. Our surrogates enable evaluation of the lightcurve posterior within minutes. We also provide built-in posterior sampling capabilities in fiesta that rely on the flowMC package which scale efficiently to higher dimensions when adding up to tens of nuisance sampling parameters. Because of its use of the JAX framework, fiesta also allows for GPU acceleration during both surrogate training and posterior sampling. We apply our framework to reanalyze AT2017gfo/GRB170817A and GRB211211A using our surrogates, thus employing the new pyblastafterglow model for the first time in Bayesian inference.

Alfvén waves are known to be important carriers of magnetic energy that could play a role in coronal heating and/or solar wind acceleration. As these waves are efficient energy carriers, how they are dissipated still remains one of the key challenges. Using a series of 1.5-D magnetohydrodynamic (MHD) simulations, we explore wave energy trapping associated with field-aligned density enhancements. We examine the parameters which govern the wave reflection and trapping. The goal of our simulations is to find optimal conditions for wave trapping, which would ultimately promote the energisation of the solar atmosphere. In agreement with previous studies, we find that maximum wave reflections happen only for a narrow range of density enhancement widths, namely when it is comparable to the Alfvén wave wavelength. In our paper, we explain this scale-selectivity using a semi-analytical model that demonstrates the importance of wave interference effects. As expected, we find that spatially extended regions of density inhomogeneities favour enhanced wave reflection and trapping. However, wave interference causes saturation of the reflected energy for very extended regions of varying density.

One of the most fundamental relationships in modern cosmology is the cosmic distance duality relation (CDDR), which describes the relationship between the angular diameter distance ($D_{\rm A}$) and the luminosity distance ($D_{\rm L}$), and is expressed as: $\eta(z)=D_{\rm L}(z)(1+z)^{-2}/D_{\rm A}(z)=1$. In this work, we conduct a comprehensive test of the CDDR by combining baryon acoustic oscillation (BAO) data from the SDSS and DESI surveys with type Ia supernova (SN) data from PantheonPlus and DESY5. We utilize an artificial neural network approach to match the SN and BAO data at the same redshift. To explore potential violations of the CDDR, we consider three different parameterizations: (i) $\eta(z)=1+\eta_0z$; (ii) $\eta(z)=1+\eta_0z/(1+z)$; (iii) $\eta(z)=1+\eta_0\ln(1+z)$. Our results indicate that the calibration of the SN absolute magnitude $M_{\rm B}$ plays a crucial role in testing potential deviations from the CDDR, as there exists a significant negative correlation between $\eta_0$ and $M_{\rm B}$. For PantheonPlus analysis, when $M_{\rm B}$ is treated as a free parameter, no evidence of CDDR violation is found. In contrast, fixing $M_{\rm B}$ to the $M_{\rm B}^{\rm D20}$ prior with $-19.230\pm0.040$ mag leads to a deviation at approximately the $2\sigma$ level, while fixing $M_{\rm B}$ to the $M_{\rm B}^{\rm B23}$ prior with $-19.396\pm0.016$ mag remains in agreement with the CDDR. Furthermore, overall analyses based on the SDSS+DESY5 and DESI+DESY5 data consistently show no evidence of the deviation from the CDDR.

We analyze the stellar light, 4000 Angstrom break and emission line profiles of 82 edge-on disk galaxies from the MaNGA survey. We characterize the stellar light profiles perpendicular to the disk plane using two parameters: a) the power law slope of the thick disk component, 2) the transition radius where the profile flattens. The 4000 Angstrom break profiles perpendicular to the plane are characterized by the number of significant changes in slope (breaks) and by the change in D_n(4000) from inner to outer disk. The slope correlates tightly with the stellar mass of the galaxy over the stellar mass range 10^9<log M*<10^10 M_sun. More massive galaxies have more extended thick disks. The slope and transition radius exhibit large scatter for galaxies more massive than 10^{10} M_sun. Half the sample have older stellar populations in their thick disks, a third have flat D$_n$(4000) profiles and 15% have younger thick disks. The D$_n$(4000) profiles exhibit as many as 4 separate breaks. There are more breaks in massive galaxies with bulges and more extended thick disks. This may indicate that the breaks are produced by more frequent accretion events in such systems. The extra-planar Halpha EQW correlates most strongly with the specific star formation rate of the galaxy, and the [OII]/Halpha ratio increases with distance from the disk plane. This increase is most apparent for massive galaxies with extended thick disk components and low SFR/M* These findings support the hypothesis that the larger [OII]/Halpha ratios may be caused by ionization from evolved stars.

Many white dwarfs are observed in compact double white dwarf binaries and, through the emission of gravitational waves, a large fraction are destined to merge. The merger remnants that do not explode in a Type Ia supernova are expected to initially be rapidly rotating and highly magnetized. We here present our discovery of the variable white dwarf ZTF J200832.79+444939.67, hereafter ZTF J2008+4449, as a likely merger remnant showing signs of circumstellar material without a stellar or substellar companion. The nature of ZTF J2008+4449 as a merger remnant is supported by its physical properties: hot ($35,500\pm300$ K) and massive ($1.12\pm0.03$ M$_\odot$), the white dwarf is rapidly rotating with a period of $\approx$ 6.6 minutes and likely possesses exceptionally strong magnetic fields ($\sim$ 400-600 MG) at its surface. Remarkably, we detect a significant period derivative of $(1.80\pm0.09)\times10^{-12}$ s/s, indicating that the white dwarf is spinning down, and a soft X-ray emission that is inconsistent with photospheric emission. As the presence of a mass-transferring stellar or brown dwarf companion is excluded by infrared photometry, the detected spin down and X-ray emission could be tell-tale signs of a magnetically driven wind or of interaction with circumstellar material, possibly originating from the fallback of gravitationally bound merger ejecta or from the tidal disruption of a planetary object. We also detect Balmer emission, which requires the presence of ionized hydrogen in the vicinity of the white dwarf, showing Doppler shifts as high as $\approx$ 2000 km s$^{-1}$. The unusual variability of the Balmer emission on the spin period of the white dwarf is consistent with the trapping of a half ring of ionised gas in the magnetosphere of the white dwarf.

The study of fast variability properties in X-ray binaries advances our understanding of the physical processes and geometric properties of the accretion flow around the compact object. In this work, we study the evolution of the timing properties of Cygnus X-1 with AstroSat/LAXPC, during the transition of the source from the hard to soft state in 2017. We use a novel frequency-segmented technique to fit simultaneously the cross spectra and part of the power spectra and coherence function with a multi-Lorentzian model and predict the phase-lags and the complementary parts of the power spectra and coherence function. We study the evolution of the frequency and power of the main variability components that are present throughout all the states. In particular, we identify previously undetected variability components, one of which manifests as a narrow dip in the coherence function and a broad drop in the phase-lag spectrum at the same frequency. This dip in coherence, which we detected for the first time in Cygnus X-1 at energies above 3 keV, appears in a state in which the source shows high-amplitude radio variability and significant hard X-ray polarization. While the contribution of the compact jet in X-rays is debated in the literature, this study provides a new avenue for investigating jet properties as well as the geometry of the Comptonizing medium.

In this work, we investigate the possibility of generating an inflationary mechanism within the framework of a metric-$f(R)$ modified gravity theory, formulated in the Jordan frame. We explore whether the scalar field, non-minimally coupled to gravity and emerging in the Jordan frame, can play the role of the primordial inflaton. Particular attention is devoted to constructing a dynamical scenario in the Jordan frame that exhibits a slow-rolling phase for the scalar field and admits a quasi-de Sitter solution for cosmic evolution. To ensure consistency with the standard cosmological model, we impose a matching condition with the $\Lambda$CDM model at the end of the inflationary phase. Furthermore, to address the problem of the absence of matter after inflation, we consider a radiation-type particle creation process that maintains an approximately constant energy density. We test our theoretical model against background observational data, specifically Pantheon+ calibrated with SH0ES and DESI calibrated with BBN. We asses its viability and discuss implications for alleviating the Hubble constant tension.

Millisecond pulsars (MSPs) are spun up during their accretion phase in a binary system. The exchange of angular momentum between the accretion disk and the star tends to align the spin and orbital angular momenta on a very short time scale compared to the accretion stage. In this work, we study a subset of $\gamma$-ray MSPs in binaries for which the orbital inclination $i$ has been accurately constrained thanks to the Shapiro delay measurements. Our goal is to constrain the observer viewing angle $\zeta$ and to check whether it agrees with the orbital inclination $i$, thus $\zeta \approx i$. We use a Bayesian inference technique to fit the MSP $\gamma$-ray light-curves based on the third $\gamma$-ray pulsar catalogue (3PC). The emission model relies on the striped wind model deduced from force-free neutron star magnetosphere simulations. We found good agreement between those two angles ($i$ and $\zeta$) for a significant fraction of our sample, between 1/2 and 2/3, confirming the spin-orbit alignment scenario during the accretion stage. However about one third of our sample deviates significantly from this alignment. The reasons are manifold: either the $\gamma$-ray fit is not reliable or some precession and external torque avoid an almost perfect alignment.

Recent JWST observations detected stellar clumps around the z=1.4 gravitationally lensed Sparkler galaxy (of stellar mass $M_*\sim 10^9\,\mathrm{M_\odot}$), with ages and metallicities compatible with globular clusters (GCs). However, most of their masses ($>10^6\,\mathrm{M_\odot}$) and sizes (>30 pc) are about 10 times those of GCs in the local Universe. To assess whether these clumps can evolve into GCs, we performed N-body simulations of their dynamical evolution from z=1.4 to z=0 (9.23 Gyr), under the effect of dynamical friction and tidal stripping. Dynamical friction is studied performing multiple runs of a clump system in a Sparkler-like spherical halo of mass $M_{200}\simeq 5\times 10^{11}\,\mathrm{M_\odot}$ (from the stellar-to-halo mass relation). For the tidal stripping, we simulated resolved clumps, orbiting in an external, static gravitational potential including the same halo as in the dynamical friction simulations and also a Sparkler-like stellar disk. Dynamical friction causes the clumps with mass $>10^7\,\mathrm{M_\odot}$ to sink into the galaxy central regions, possibly contributing to the bulge growth. In absence of tidal stripping, the mass distribution of the surviving clumps (40%) peaks at $5\times 10^6\,\mathrm{M_\odot}$, implying the presence of uncommonly over-massive clumps at z=0. Tidal shocks by the stellar disk strip considerable mass from low-mass clumps, even though their sizes remain larger than those of present-day GCs. When the surviving clumps are corrected for tidal stripping, their mass distribution peak shifts to $2\times 10^6\,\mathrm{M_\odot}$, compatible with massive GCs. Our simulations suggest that a fraction of the Sparkler clumps is expected to fall into the central regions, where they might become bulge fossil fragments or contribute to form a nuclear star cluster. The remaining clumps are too large in size to be progenitors of GCs.

We systematically extend the framework of galaxy bias renormalization to two-loop order. For the minimal complete basis of 29 deterministic bias operators up to fifth order in the density field and at leading order in gradient expansion we explicitly work out one- and two-loop renormalization. The latter is provided in terms of double-hard limits of bias kernels, which we find to depend on only one function of the ratio of the loop momenta. After including stochasticity in terms of composite operator renormalization, we apply the framework to the two-loop power spectrum of biased tracers and provide a simple result suitable for numerical evaluation. In addition, we work out one- and two-loop renormalization group equations (RGE) for deterministic bias coefficients related to bias operators constructed from a smoothed density field, generalizing previous works. We identify a linear combination of bias operators with enhanced UV sensitivity, related to a positive eigenvalue of the RGE. Finally, we present an analogy with the RGE as used in quantum field theory, suggesting that a resummation of large logarithms as employed in the latter may also yield useful applications in the study of large-scale galaxy bias.

We present a study on the detection and characterization of transient quasi-periodic oscillations (QPOs) in the $\gamma$-ray emission of blazars 4C +31.03, MG1 J123931+0443, and PKS 1622$-$253. Using light curves derived from \textit{Fermi} Large Area Telescope data, we investigate oscillatory patterns characterized by periodic multiplicative amplitudes that vary linearly over time. By segmenting the light curves into increasing and decreasing trends, we analyze each segment independently, allowing for precise measurements of both the periodicity and long-term variations. To interpret these QPOs, we explore various theoretical scenarios that could explain their origin and underlying physical mechanisms. The estimated periods for blazars MG1 J123931+0443 and PKS 1622$-$253 are consistent with precessional dynamics expected in binary supermassive black hole systems, although our current results do not allow for a definitive conclusion.

We investigate the James Webb Space Telescope (JWST) MIRI MRS gas molecular content of an externally irradiated Herbig disk, the F-type XUE 10 source, in the context of the eXtreme UV Environments (XUE) program. XUE 10 belongs to the massive star cluster NGC 6357 (1.69 kpc), where it is exposed to an external far-ultraviolet (FUV) radiation $\approx$ 10$^3$ times stronger than in the Solar neighborhood. We modeled the molecular features in the mid-infrared spectrum with Local Thermodynamic Equilibrium (LTE) 0D slab models. We derived basic parameters of the stellar host from a VLT FORS2 optical spectrum using PHOENIX stellar templates. We detect bright CO2 gas with the first simultaneous detection (> 5$\sigma$) of four isotopologues (12CO2, 13CO2, 16O12C18O, 16O12C17O) in a protoplanetary disk. We also detect faint CO emission (2$\sigma$) and the HI Pf$\alpha$ line (8$\sigma$). We also place strict upper limits on the water content, finding a total column density $\lesssim$ 10$^{18}$ cm$^{-2}$. The CO2 species trace low gas temperatures (300-370 K) with a range of column densities of 7.4 $\times$ 10$^{17}$ cm$^{-2}$ (16O12C17O)-1.3 $\times$ 10$^{20}$ cm$^{-2}$ (12CO2) in an equivalent emitting radius of 1.15 au. The emission of 13CO2 is likely affected by line optical depth effects. 16O12C18O and 16O12C17O abundances may be isotopically anomalous compared to the 16O/18O and 16O/17O ratios measured in the interstellar medium and the Solar System. We propose that the mid-infrared spectrum of XUE 10 is explained by H2O removal either via advection or strong photo-dissociation by stellar UV irradiation, and enhanced local CO2 gas-phase production. Outer disk truncation supports the observed CO2-H2O dichotomy. A CO2 vapor enrichment in 18O and 17O can be explained by means of external UV irradiation and early on (10$^{4-5}$ yr) delivery of isotopically anomalous water ice to the inner disk.

Exotic dark matter (EDM) refers to a dark matter species whose equation of state deviates from zero at late times. This behavior enables it to model a variety of non-standard late-time cosmologies, offering alternatives to various dark energy (DE) models, especially when the DE sector violates the null energy condition. In this work, by fitting to CMB, BAO, and Supernovae (SNe) data and comparing models in a Bayesian approach, we show that simple models of exotic dark matter are statistically comparable to the $w_0w_a$CDM DE model in explaining the recent anomaly in the late-time cosmological evolution suggested by DESI and supernova observations, although in both classes of models the evidence against the $\Lambda$CDM model only appears when the DES-Y5 or Union3 SNe dataset is included. The value of $H_0$ remains similar to that in the DE model, except in the no-SNe case, where the DE model predicts lower values than $\Lambda$CDM, thereby worsening the Hubble tension, whereas the EDM models yield values closer to that of $\Lambda$CDM, albeit with larger uncertainty. In addition, the EDM models predict a drastically different energy budget for the present-day universe compared to the standard model, and provide an explanation for a coincident problem in the DE-model explanation of the DESI anomaly.

The fact that strange sea quarks are abundant in the nucleons, but without strangeness, is of great importance for understanding the nature of dense matter condensed by the strong interaction, particularly in the context of the ``gigantic nucleus'' formed by the gravitational collapse of an evolved massive star. We hypothesize that the basic unit of bulk strong matter with the light-flavor symmetry of valence quarks is ``strangeon'', which is the counterpart of the nucleon found in atomic nuclei. In addition to strangeon stars (SnSs) with large baryon number of $A\approx 10^{57}$, strange nuggets (SnNs) with $A\gtrsim 10^{10}$ could also exist in the Universe. Both the SnS and the SnN are explained, with attention to their observational evidence.

We investigate gravitational lensing in the strong-field regime for a black hole in F(R) Euler Heisenberg Gravity with Rainbow gravity modifications. This black hole spacetime is characterized by the Euler Heisenberg parameter, F(R) parameters, Rainbow functions, the black hole charge Q, and mass M . We numerically compute the strong deflection angle and its associated coefficients, and explore their astrophysical implications for supermassive black holes in different galaxies. Our findings show that increasing the Euler Heisenberg parameter enhances key lensing observables such as the photon sphere radius, critical impact parameter, angular position, relative magnification, Einstein ring radius, and time delay for a fixed Q . Conversely, increasing Q decreases these parameters while keeping other quantities fixed. Additionally, a higher Euler Heisenberg parameter reduces the deflection angle and angular separation S , whereas increasing Q causes both to increase. Our study reveals that black holes in this modified gravity framework can act as strong gravitational lenses, producing deflection angles that surpass those of Reissner Nordstrom and Schwarzschild black holes under specific conditions. These results highlight the unique topological features of modified charged black holes and suggest their potential as astrophysical candidates, offering new insights into their observational signatures. Furthermore, we analyze the sensitivity of the lensing predictions with respect to changes in the functional forms of the Rainbow functions. The combined effects of F(R) gravity, Euler Heisenberg electrodynamics, and Rainbow gravity introduce novel qualitative features not present in the individual models.

We present VLT/MUSE observations of comet 67P/Churyumov-Gerasimenko during its 2021 perihelion passage, from which we generated simultaneous maps of dust, [OI], C$_2$, NH$_2$, and CN comae across 12 pre- and post-perihelion epochs. These maps reveal the evolutionary and compositional trends of 67P's coma and further enrich the context and findings of ESA's Rosetta mission. Dust and gas species displayed distinct structures, where NH$_2$ and CN signals were uniquely associated with known dust fans, raising the question of possible correlation to the dust and contributions of extended sources. Localised fitted NH$_2$ scale lengths were 1.5-1.9$\times$ larger than those fitted for the rest of the coma, which is consistent with an extended source component for northern pre-perihelion emissions. In the southern hemisphere, CN was correlated with a prominent and sharp dust structure, potentially revealing an extended source origin via larger dust particles that preserve the CN parent species, as evidenced by higher spectral slopes in the region. Gas maps depicted two distinct evolutionary regimes: (1) evolving H$_2$O ([OI]$^{1}$D) and C$_2$ emissions driven by nucleus sublimation and subsolar insolation, and (2) stable NH$_2$ and CN emissions associated with seasonal dynamics and possible distributed sources. Dust spectral slope maps revealed spectral slope trends consistent with Rosetta findings, while green/red [OI] ratios generally indicate a coma dominated by H$_2$O.

Recent observations of baryon acoustic oscillations (BAO) from the Dark Energy Spectroscopic Instrument (DESI) survey, when combined with measurements of the cosmic microwave background (CMB) and Type Ia supernovae (SNIa), provide compelling evidence for a phantom crossing at late times, along with statistically significant deviations from the standard $\Lambda$CDM model. In this work, we investigate the role of redshift-space galaxy clustering data by employing the pre-reconstruction full-shape (FS) galaxy power spectrum from the Baryon Oscillation Spectroscopic Survey (BOSS) data release 12 (DR12) sample. This dataset is analyzed in combination with BAO measurements from DESI data release 2 (DR2) and various SNIa samples. Our analysis demonstrates that the joint combination of these datasets can yield deviations from $\Lambda$CDM at a significance level of up to $\sim 5\sigma$, suggesting strong indications that the dark energy equation of state parameter $w(z)$ may have crossed the phantom divide ($w = -1$) in the redshift range $z \sim 0.4$--$0.5$. The precise location and strength of this crossing depend on the adopted theoretical parameterizations. Importantly, our results reveal that this trend persists even in the absence of CMB data, underscoring the robustness of the FS power spectrum as a powerful and independent probe for testing dark energy models and for distinguishing between competing cosmological scenarios.

The inflated radii observed in hundreds of hot Jupiters (HJ) represent a long-standing open issue. In this study, we quantitatively investigate this phenomenon within the framework of Ohmic dissipation arising from magnetic induction in the atmosphere, one of the most promising mechanisms for explaining the radius anomaly. We simulate the evolution of irradiated giant planets with {\tt MESA}, spanning the observed range of masses and equilibrium temperatures, incorporating an internal source of Ohmic dissipation that extends to deep layers of the envelope. We infer average atmospheric wind intensities, averaged in the region $p\lesssim 10$ bar, in the range 0.01-1 km/s in order to reproduce the range of observed radii, decreasing roughly linearly with planetary mass, and much more steeply with equilibrium temperature. This is consistent with the expected effects of magnetic drag from the induced field, which is higher for more intense irradiation, via conductivity, and for larger masses, which have higher dynamo fields. Due to the evolution of the dynamo field and the proportionality of the induced currents on it, the Ohmic efficiency typically decreases by at least one order of magnitude from 0.1 to 10 Gyr, at contrast with the common assumption of a constant-in-time value. Notably, the extent of the main convective region, and the associated heat flux supporting the dynamo, is reduced in the presence of strong Ohmic dissipation, which in turn depends on the dynamo field strength, generating a non-trivial coupling of the latter with the atmospheric induction, potentially leading to an oscillatory behaviour of the field strength. These findings remain generally valid even when accounting for a long-term increase in the main-sequence host star luminosity, although this case can more readily lead to HJ re-inflation, consistent with previous studies.

Magnetic reconnection in current layers that form intermittently in radiatively inefficient accretion flows onto black holes is a promising mechanism for particle acceleration and high-energy emission. It has been recently proposed that such layers, arising during flux eruption events, can power the rapid TeV flares observed from the core of M87. In this scenario, inverse Compton scattering of soft radiation from the accretion flow by energetic electron-positron pairs produced near the reconnection layer was suggested as the primary emission mechanism. However, detailed calculations show that radiation from pairs alone cannot account for the GeV emission detected by the Fermi observatory. In this work, we combine analytic estimates with 3D radiative particle-in-cell simulations of pair-proton plasmas to show that the GeV emission can be naturally explained by synchrotron radiation from protons accelerated in the current sheet. Although the exact proton content of the layer is uncertain, our model remains robust across a broad range of proton-to-pair number density ratios. While protons are subdominant in number compared to pairs, our simulations demonstrate that they can be accelerated more efficiently, leading to a self-regulated steady state in which protons dominate the energy budget. Ultimately, proton synchrotron emission accounts for approximately 5% to 20% of the total dissipation power. The majority is radiated as MeV photons via pair synchrotron emission, with a smaller fraction emitted as TeV photons through inverse Compton scattering.

Gamma-Ray Bursts (GRBs) are known to be unpredictable in time and position. A few (observationally) exceptional events have been observed, as GRB221009A that stands out for its fluence and peak flux, being orders of magnitude higher than what measured so far. Analyzing the observed fluence, peak flux or duration distributions typically requires one to assume some scenarios, and the consistency of the observed data with the predictions turns out to be an important model diagnostic. However, it is also of interest to model these distributions using general statistical properties that do not rely on specific model assumptions, allowing one to derive inferences only based on the consistency of the observed distributions with the hypothesis of one single population of events that generate them. We obtained fluences, peak fluxes and durations from the catalogues of GRBs observed by the CGRO-BATSE and Fermi-GBM instruments. We selected the extreme values in slots of equal duration and modelled their distributions by the generalized extreme value (GEV) formalism. The GEV distribution is a limit distribution naturally arising when the number of observations is large and is essentially independent of the phenomena producing the observed data. The distributions of extreme values for fluences, peak fluxes and durations are consistent with being extracted from a single population of events but the fluence and peak flux recorded for GRB221009A constitutes a striking exception. The probability to observe such an event, assuming it is a cosmological GRB, is low, with a median value of about one event per millennium for the fluence and about one event per century for the peak flux.

We analyze two different scenarios for the late Universe dynamics, resulting into Hubble parameters deviating from the $\Lambda$CDM, mainly for the presence of an additional free parameter, which is the dark energy parameter. The first model consists of a pure evolutionary dark energy paradigm, as result of its creation by the gravitational field of the expanding Universe. The second model also considers an interaction of the evolutionary dark energy with the matter component, postulated via the conservation of the sum of their ideal energy-momentum tensors. These two models are then compared \textit{via} the diagnostic tool of the effective running Hubble constant, with the binned data of the so-called ``Master sample'' for the Type Ia Supernovae. The comparison procedures, based on a standard MCMC analysis, led to a clear preference of data for the dark energy - matter interaction model, which is associated to a phantom matter equation of state parameter (very close to $-1$) when, being left free by data (it has a flat posterior), it is fixed in order to reproduce the decreasing power-law behavior of the effective running Hubble constant, already discussed in literature.

Classifying stars, galaxies, and quasars is essential for understanding cosmic structure and evolution, but the vast data from modern surveys make manual classification impractical, while supervised learning methods remain constrained by the scarcity of labeled spectroscopic data. We aim to develop a scalable, label-efficient method for astronomical classification by leveraging semi-supervised learning (SSL) to overcome the limitations of fully supervised approaches. We propose a novel SSL framework combining K-means clustering with Random Forest classification. Our method partitions unlabeled data into 50 clusters, propagates labels from spectroscopically confirmed centroids to 95\% of cluster members, and trains a Random Forest on the expanded pseudo-labeled dataset. We apply this to the CPz catalog, containing multi-survey photometric and spectroscopic data, and compare performance with a fully supervised Random Forest. Our SSL approach achieves F1-scores of 98.8\%, 98.9\%, and 92.0\% for stars, galaxies, and quasars, respectively, closely matching the supervised method with F1-scores of 99.1\%, 99.1\%, and 93.1\% while outperforming traditional color-cut techniques. The method demonstrates robustness in high-dimensional feature spaces and superior label efficiency compared to prior work. This work highlights SSL as a scalable solution for astronomical classification when labeled data is limited, though performance may degrade in lower-dimensional settings.

Our understanding of the chemical evolution of galaxies has advanced through measurements from both distant galaxies across redshift, and our own Milky Way (MW). To form a comprehensive picture, it is essential to unify these constraints, placing them on a common scale and parlance and to understand their systematic differences. In this study, we homogenize oxygen and iron measurements from star-forming galaxies at Cosmic Noon ($z{\sim}2-3$) with resolved stellar abundances from the Local Group. The MW is divided into four components, assuming the outer halo is dominated by debris from the Gaia-Sausage-Enceladus (GSE) progenitor. After converting all abundances to a common Solar scale, we identify clear $\alpha$- and iron-enhancement trends with mass in the $z{\sim}2-3$ galaxies and find good agreement between these galaxies and the MW high-$\alpha$ disc in [O/Fe] vs. [Fe/H]. We also find excellent agreement between the [O/Fe] trends seen in the MW high- and low-$\alpha$ discs with O-abundances seen in old and young planetary nebulae in M~31 respectively, supporting the existence of $\alpha$-bimodality in the inner regions of M~31. Finally, we use globular cluster ages to project the MW and GSE back in time to $z{\sim}3$ and find that their estimated mass, oxygen and iron abundances are strikingly consistent with the mass-metallicity relation of star-forming galaxies at $z{\sim}3$. In the future, increased transparency around the choice of Solar scale and abundance methodology will make combining chemical abundances easier -- contributing to a complete picture of the chemical evolution of all galaxies.

Cosmic voids have emerged as powerful probes for cosmology, providing complementary information on the large-scale structure of the universe. We present the first application of a hydrodynamical framework to model the evolution of cosmic voids. This approach offers a physically intuitive characterization of void dynamics and can naturally be applied to non-standard cosmologies. We derive the cosmology-dependent mapping that relates the linear (Lagrangian) and fully non-linear (Eulerian) evolution of the matter density contrast, a central component for accurate theoretical modeling of void statistics. Furthermore, we present a new method for determining the shell-crossing epoch across arbitrary cosmological backgrounds, thereby extending previous treatments restricted to the Einstein-de Sitter universe. Motivated by recent DESI results hinting at dynamical dark energy, we investigate void evolution in $ w_0w_a$CDM cosmologies by varying $ w_0$ and $w_a$. We also consider the impact of varying the matter density parameter, $ \Omega_{\mathrm{m},0}$. We find that the evolution of isolated, spherically symmetric cosmic voids is most sensitive to $ \Omega_{\mathrm{m},0} $ and $ w_0 $, which can alter the non-linear density contrast by up to 20-30%. Variations in $w_a$ have a smaller impact, but may still lead to measurable effects. We also show that the cosmology-dependent mapping between linear and non-linear density contrasts may provide a sensitive probe of dynamical dark energy in precision void analyses.

Double white dwarf binaries are an important remnant of binary evolution as they are possible type Ia supernova progenitors and strong sources of gravitational waves in the low-frequency regime. The double-lined double white dwarf (DBL) survey searches for compact double white dwarfs where both stars are spectrally disentangleable. Candidates are identified by being overluminous compared to the cooling sequence of a typical mass, single white dwarf. In this second DBL survey instalment, we present full orbital solutions of 15 double white dwarf binaries from our ongoing campaign to accurately measure a magnitude-limited mass-period distribution. 12 of these systems are fully solved for the first time. A long-standing bias in the full population has been evident, favouring systems with orbital periods up to a few hours, with little exploration of the majority of the compact double white dwarf population, whose orbital period distribution centres at approximately 20hr. The 15 systems in this study span the orbital period range 5-75hr, significantly augmenting the number of well-characterised systems over these periods, and in general have two similar mass stars combining to approximately 1.0 solar masses. We witness that the orbitally derived mass ratios generally show an excellent agreement with those deduced from atmospheric fits to double-lined spectra in previous work, emphasising the power of wide-scale spectroscopic surveys to efficiently locate the highest mass, double-lined double white dwarfs in the local Galaxy.

Massive young clusters with rich populations of high-mass stars are ideal laboratories to explore their evolutionary paths. Despite being the most prominent cluster in the Perseus-arm Cas OB8 association, NGC 663 remains comparatively little studied. We present a comprehensive investigation of its properties, integrating astrometric, photometric and spectroscopic data for the cluster and its surroundings, including accurate spectral classification for over 150 members. Gaia astrometry indicates over 300 B-type members, possibly rendering NGC 663 the most massive cluster in the Perseus arm, with initial mass likely exceeding 10000 M_\sun . This large population makes NGC 663 an excellent laboratory for studying massive star evolution. Spectral analysis of the earliest members reveals approximately solar metallicity and a turn-off mass of approximate 8.5 M_\sun, consistent with the photometric age of 23 Ma. We identify five spectroscopic blue stragglers, including the Be/X-ray binary RX J0146.9$+$6121. We outline its evolutionary history and compare its properties with other Be stars. Although the cluster contains many Be stars, their relative fraction is not particularly high. Intriguingly, four of the six blue supergiant members appear to have significantly higher masses than the brightest giants near the Hertzsprung gap. These observations suggest that most mid-B supergiants may form via mergers, unless stars of 10-12 M_\sun born as primaries in binaries rarely undergo supernova explosions. Similarly, if Be stars form through the binary channel, then either most are produced through case A evolution or supernovae are uncommon among primaries in this mass range.

We present an updated atlas of spectra from the Low-Resolution Spectrometer (LRS) on the Infrared Astronomical Satellite (IRAS), which took spectra from 7.67 to 22.73 um with a spectral resolving power (lambda / Delta lambda) of 20-60. The updated atlas includes 11,238 spectra, including 5425 spectra published in the original LRS Atlas, 5796 spectra published in three later papers, and 17 spectra previously available online but not published. The updated atlas has significantly more sources close to the Galactic plane than the original atlas. We have applied an improved spectral correction to remove an artifact at 8 um in the original database. While the IRAS mission flew over 40 yr ago, the extended LRS atlas remains the single most complete database of mid-infrared spectra of nearby and bright objects in the Galaxy.

We present new constraints on the halo masses and matter density profiles of DESI galaxy groups by cross-correlating samples of Luminous Red Galaxies (LRGs) and Bright Galaxy Survey (BGS) galaxies with the publicly available CMB lensing convergence map from ACT DR6. This provides an independent, lensing-based calibration of halo masses, complementary to methods relying on clustering or dynamics. We derive constraints on the mean halo mass for three DESI-selected samples, finding $\log(M_{\rm halo}/(M_\odot/h)) \approx 13.18$, 13.03 and 13.02 for the Main LRG, Extended LRG, and BGS samples, respectively. Using a halo model approach, we also compare the projected galaxy-matter density profiles with previously reported gas profiles inferred from measurements of the kinematic Sunyaev-Zel'dovich (kSZ) effect. This work addresses one of the key uncertainties in interpreting kSZ signals -- the unknown host halo mass distribution -- by providing an independent and consistent mass calibration. The agreement between the gas and total mass profiles at large aperture suggests that sufficiently far from the group center (2--3 virial radii), we recover all the baryons, offering a resolution to the 'missing baryon' problem. We further study the cumulative gas fractions for all galaxies as well as for the most massive galaxy groups in the sample ($\log(M_{\rm halo}/(M_\odot/h)) \approx 13.5$), finding values that are physically sensible and in agreement with previous findings using kSZ and X-ray data: compared to the TNG300 simulation, the observed gas fractions are systematically lower at fixed radius by $\gtrsim$4$\sigma$, providing compelling, independent evidence for stronger baryonic feedback in the real Universe. These findings highlight the power of combining CMB lensing with galaxy surveys to probe the interplay between baryons and dark matter in group-sized halos.

The physics responsible for the production of observed high-energy neutrinos have not been established so far, neither for the diffuse astrophysical ones nor for those detected from single blazars. We recently proposed that both could be explained by deep inelastic scatterings between sub-GeV dark matter (DM) around blazars and protons within their jets. Here, we compute the proton-recoil signals at the neutrino detectors Super-Kamiokande, KamLAND, Borexino, JUNO, Hyper-Kamiokande and DUNE induced by DM that is itself boosted by the scatterings with protons in blazar jets. We do it for the four cases of vector, axial, scalar and pseudoscalar mediators of DM-quark interactions. We perform the analysis for the single blazar TXS 0506+056 and for a sample of more than 300 stacked blazars. We find that searches for such blazar-boosted DM leave room for a variety of DM models to explain observations of high-energy neutrinos. We check that the depletion of the DM spike induced by DM-proton and DM-DM interactions does not compromise the DM interpretation for high-energy neutrinos, but challenges other blazar-DM signals.

The decay of light bosonic dark matter particles can produce a bright electron/positron ($e^+e^-$) flux that can be strongly constrained by local Voyager observations of the direct $e^+e^-$ flux, as well as 511 keV Line and X-ray continuum observations of $e^+e^-$ emission. We carefully analyze the $e^+e^-$ yield and resulting cosmic-ray and X-ray spectra from theoretically well-motivated light dark matter models, including: (a) electrophilic axion-like particles, (b) dark photons, (c) scalars, and (d) $B-L$ and $L_{i}-L_{j}$ vector bosons. We use the morphology and spectrum of the INTEGRAL 511 keV line data, the eROSITA X-ray continuum spectrum and the Voyager $e^+e^-$ spectrum to constrain the decay lifetime and coupling of each dark matter model. We find that 511 keV observations typically set world-leading limits on bosonic dark matter decay below masses of $\sim$1 GeV, while eROSITA observations provide the strongest constraints in the range from 1--10 GeV. Finally, we forecast future limits from 21 cm line searches with next-generation HERA data.

Recent study[1] has suggested that warm inflation may be realized with Standard Model gauge interaction alone. Motivated by this framework, we investigate the gravitational wave spectrum and graviton-portal dark matter production through the freeze-in process generated during warm inflation scenarios. We perform a comparative analysis for different dissipation terms, focusing on their distinct gravitational wave signatures in the high-frequency regime. Our findings reveal qualitative and quantitative differences in the spectral behavior, offering a preliminary pathway for discriminating among various inflationary and dark matter models through high-frequency gravitational wave signals.

The production site of heavy r-process elements, such as Gold and Uranium, is uncertain. Neutron star mergers are the only astrophysical phenomenon in which we have witnessed their formation. However, the amount of heavy elements resulting from the merger remains poorly constrained, mainly due to uncertainties on the mass and angular momentum of the disk formed in the merger remnant. Matter accretion from the disk is also thought to power gamma ray-bursts. We discover from numerical relativity simulations that the accretion disk influences the ringdown gravitational-wave signal produced by binaries that promptly collapse to black-hole at merger. We propose a method to \emph{directly} measure the mass of the accretion disk left during black hole formation in binary mergers using observatories such as the Einstein Telescope or Cosmic Explorer with a relative error of 10\% for binaries at a distance of up to 30~Mpc, corresponding to an event rate of 0.001 to 0.25 events per year.