Papers for Friday, Sep 12 2025

Zodiacal light -- arising from both the thermal emission and scattered sunlight by interplanetary dust -- is the dominant component of the sky background in NIRISS Single Object Slitless Spectroscopy (SOSS) observations. The GR700XD grism disperses the zodiacal background across multiple diffraction orders, producing a characteristic two-dimensional, spectral order-dependent background pattern on the detector that reflects the combined contribution of overlapping orders. Observations reveal significant variability in background intensity driven by JWST's sky pointing and seasonal changes, underscoring the need for precise background subtraction during SOSS data reduction. Current methods rely on a generic background template derived during commissioning, scaled to match individual exposures. However, mis-scaled templates can leave behind structured residuals that may compromise the precision of exoplanet transit depth measurements. To improve background modeling, we conducted two calibration programs (PID 4479 and 6658) using the FULL frame readout mode and a 5-row by 2-column mosaic pattern to sample a range of sky positions. These observations enable empirical reconstruction of the sky background and provide detailed insights into its spatial and spectral characteristics. We present a library of empirically derived background templates and evaluate their performance, alongside the current template, by benchmarking against the PID 2113 dataset, which includes contemporaneous background exposures. This work supports aims to enhance background subtraction for SOSS time-series to achieve higher-precision exoplanet spectroscopy with JWST.

Future lunar missions will depend on an internationally agreed upon timescale that remains accurate under the Moon's unique gravitational environment and its orbital dynamics. This thesis investigates the proposed Lunar Coordinate Time (TCL), derived analogously to Geocentric Coordinate Time (TCG) and thus aligned with current IAU proposals. We first formalise the TCL transformation and quantify its characteristics from solar system simulations. Next, we compute stationary surface-clock drifts caused by gravitational redshift and the Moon's changing orientation parameters, evaluating how accurate atomic clocks deployed on the surface of the Moon (much like for ESA's proposed NovaMoon mission) would have to be to measure these effects. Finally, we simulate relativistic proper time for ESA's Moonlight navigation satellites, identifying average drift and harmonic variations, to better understand the system that will comprise and enable a Lunar PNT (Positioning, Navigation and Timing) architecture. These kinds of investigations are an essential step toward a sustained internationally cooperative operation at the lunar south pole and beyond.

White dwarfs (WDs) receive natal velocity boosts of $\sim 1 \, \mathrm{km \, s^{-1}}$ due to recoil from asymmetric mass loss during the late asymptotic giant branch (AGB) stage. In a wide binary, the acceleration of a proto-WD exerts a torque, modifying the orbital eccentricity. Potential signatures of this effect have been detected in Gaia's astrometric binary sample. At the same time, an AGB star's puffy envelope facilitates strong tidal interactions in binaries with periapsis separations of a few AU, capturing the companion into a tighter orbit and potentially driving the system towards a common-envelope phase. Using an analytical model for wide binary evolution under asymmetric AGB mass loss, we find that recoil can induce tidal interactions in up to $30\%$ of initially wide binaries on the AGB or post-AGB for initial separations of $\sim 100 \mbox{--} 1000$ AU. We relate these interactions to three recent observational puzzles: (i) The dearth of wide WD+MS and WD+WD binaries in Gaia DR3 with eccentricities $\gtrsim 0.9$. (ii) The formation of moderately eccentric WD+MS and WD+WD binaries with orbital periods of $\sim 100 \mbox{--} 1000$ days, which may happen via a high-eccentricity common-envelope phase. (iii) The origin of low-luminosity, long-timescale, dust-obscured transients towards AGB progenitors, such as the ongoing event WNTR23bzdiq in M31. Our findings have potential implications for the survival and dynamical evolution of planetary systems around WD progenitors, to be investigated in future works.

The epoch of reionization (EoR) marks the last phase transition of hydrogen in our Universe, as it evolves from cold and neutral to hot and ionized in the intergalactic medium (IGM). While its endpoint and duration can be estimated from current observations, albeit with large uncertainties, there is no known avenue to constrain its onset. We propose a novel method based on the Pearson cross-correlation coefficient between 21-cm brightness temperature maps and line-intensity maps tracing star-formation (e.g., OIII, CO, CII). This real-space estimator evolves from negative to positive as X-ray heating progresses, and saturates prior to the EoR. We predict a sharp turnover from saturation during the earliest EoR stages, when the IGM ionized fraction reaches $\bar{x}_{\rm HII}\sim 1\%-10\%$. We show that in standard scenarios, where the IGM heating precedes reionization, the turnover is a clear, model-robust signature. Its detection will provide a unique observational anchor for the EoR onset, complementing existing probes and tightening constraints on early galaxy formation models.

Understanding the interaction of dense, cold ejecta clumps with a fast reverse shock, an instance of the "cloud-crushing" problem, is essential to assess whether core-collapse supernovae act as net dust factories or net dust destroyers. This work assesses whether dusty ejecta clumps are destroyed by the reverse shock or instead cool, condense, and grow in mass under realistic supernova-remnant conditions. Cloud-crushing timescales are computed and compared to radiative cooling timescales, including both gas-phase cooling and dust-induced cooling, for a large grid of clump densities, dust-to-gas mass ratios, and shock velocities. When the dust-to-gas mass ratio exceeds $10^{-3}$, gas-grain collisions become efficient enough that the cooling timescale $t_{\rm cool}$ falls below the cloud-crushing timescale $t_{\rm cc}$ over a broad span of clump densities and shock velocities, enabling dusty clumps to survive even fast reverse shocks. For example, at clump densities $\geq 2 \times 10^{4}$ cm$^{-3}$, dust-to-gas mass ratios $\sim 10^{-2}$, and shock velocities up to $2000$ km s$^{-1}$, enhanced gas-grain cooling drives the system into a regime where dusty clumps can gain additional cold mass and increase their dust masses. Strong radiative cooling can shield dust-rich clumps in supernova remnants, enabling a significant fraction of ejecta dust to be injected into the interstellar medium. These results mirror the "growth" regime found in studies of circumgalactic clouds and rapidly cooling shocked stellar winds, implying a larger dust survival in supernova remnants. Indeed, the dusty globules seen in the Crab Nebula occupy the predicted survival regime across a wide range of physical parameters.

We report the discovery of an extremely high-velocity outflow (EHVO) in the most luminous ($L\rm_{Bol}$ $\sim$ 2.29 $\times$ 10$^{48}$ erg/s) QSO, SMSS J2157-3602, at z=4.692. Combined XSHOOTER and NIRES observations reveal that the EHVO reaches a maximum velocity of v$_\mathrm{max} \sim 0.13c$ and persists over rest-frame timescales of a few months to one year. SMSS J2157-3602 also exhibits one of the highest balnicity index discovered in an EHVO so far. In addition, the blueshifted CIV emission traces a high-velocity (v$\rm_{CIV}^{50}\sim$ 4660 km/s) outflow from the broad-line region. Thanks to an XMM-Newton observation, we also discover the X-ray weak nature of this QSO, which likely prevents the overionization of the innermost disk atmosphere and facilitates the efficient launch of the detected EHVO and BLR winds. The extraordinary luminosity of SMSS J2157-3602 and the extreme velocity of the EHVO make it a unique laboratory for testing AGN driven feedback under extreme conditions. Current uncertainties on the outflow's location and column density strengthen the case for dedicated follow-up, which will be essential to assess the full feedback potential of this remarkable quasar.

The Optical Gravitational Lensing Experiment (OGLE) is one of the most productive and influential photometric sky surveys in the history of observational astronomy. Originally designed to detect dark matter through gravitational microlensing events, OGLE has evolved into a cornerstone of time-domain astrophysics, delivering three decades of two-band, high-cadence observations of approximately two billion stars across the Galactic bulge, disk, and Magellanic System. This review summarizes OGLE's key contributions to variable star research, including the discovery, classification and characterization of pulsating stars, eclipsing, ellipsoidal, and rotating variables, or irregular and eruptive stars. Particular emphasis is placed on the OGLE Collection of Variable Stars (OCVS), a publicly available and systematically expanded dataset that has become a fundamental resource for studies of stellar variability and evolution, Milky Way and other galaxies structure, microlensing, compact objects, exoplanets and more. The synergy between OGLE and other major sky surveys, including ASAS, ASAS-SN, ATLAS, Gaia, KMTNet, MACHO, MOA, TESS, PLATO, or ZTF further amplifies its scientific reach.

Tidal disruption events (TDEs) occur when a star is disrupted by the tidal forces of a supermassive black hole (SMBH), producing bright multi-wavelength flares. Among these events, AT2020mot has so far exhibited the highest recorded optical polarisation, with tidal shocks proposed as the primary source of its polarised emission. We present a comprehensive analysis of 13 TDEs with available polarimetric observations, aiming to determine whether the unusually high polarisation of AT2020mot stems from unique physical processes or arises from mechanisms shared by other TDEs. We present new optical polarisation measurements of TDEs obtained from multiple ground-based telescopes, combining them with optical, UV, and X-ray light curves from the Zwicky Transient Facility and the Swift observatory, respectively. We derive intrinsic TDE properties, such as SMBH and stellar masses, using MOSFiT and TDEMass, and compare them with the ones of the sample population. Our population study reveals that AT2020mot agrees with the broader TDE sample in most physical properties, including blackbody temperature, luminosity, and rise timescales. However, its optical polarisation degree is exceptionally high compared to the low or undetected polarisation observed in other events. Additionally, AT2020mot appears to have an elevated column density from our MOSFiT fits, suggesting a more complex environment than is typically assumed. We conclude that although AT2020mot fits well within the general TDE population in terms of global characteristics, its extraordinarily high polarisation and higher column density challenge current models based purely on shock or reprocessing mechanisms. More extensive, time-resolved polarimetric monitoring of newly discovered TDEs will be critical to determine whether AT2020mot represents an outlier or the extreme end of a continuum of TDE properties.

FRBs constitute a unique probe of various astrophysical and cosmological environments via their characteristic dispersion and rotation (RM) measures that encode information about the ionized gas traversed by the FRB sightlines. In this work, we analyse observed RM measured for 14 localized FRBs at $0.05 \lesssim z \lesssim 0.5$, to infer total magnetic fields in various galactic environments. Additionally, we calculate $f_{\rm gas}$ - the average fraction of halo baryons in the ionized CGM. We build a spectroscopic dataset of FRB foreground galaxy halos, acquired with VLT/MUSE and FLIMFLAM survey. We develop a novel Bayesian algorithm and use it to correlate the individual intervening halos with the observed RM. This approach allows us to disentangle the magnetic fields present in various environments traversed by the FRB. Our analysis yields the first direct FRB constraints on the strength of magnetic fields in the ISM and halos of the FRB host galaxies, as well as in halos of foreground galaxies. We find that the average magnetic field in the ISM of FRB hosts is $B_{\rm host}^{\rm local} = 5.44^{+1.13}_{-0.87}\mu{\rm G}$. Additionally, we place upper limits on average magnetic field in FRB host halos, $B_{\rm host}^{\rm halo} < 4.81\mu{\rm G}$, and in foreground intervening halos, $B_{\rm f/g}^{\rm halo} < 4.31\mu{\rm G}$. Moreover, we estimate the average fraction of cosmic baryons inside $10 \lesssim \log_{10} \left( M_{\rm halo} / M_{\odot}\right) \lesssim 13.1$ halos $f_{\rm gas} = 0.45^{+0.21}_{-0.19}$. We find that the magnetic fields inferred in this work are in good agreement with previous measurements. In contrast to previous studies that analysed FRB RMs and have not considered contributions from the halos of the foreground and/or FRB host galaxies, we show that they can contribute a non-negligible amount of RM and must be taken into account when analysing future FRB samples.

Recently discovered massive $\alpha$-enriched (MAE) stars have surface chemistry consistent with members of the older thick disk Milky Way population but high masses ($\sim$ 1.5 - 3 M$_{\odot}$) that suggest these stars are young. The origin of MAE stars is not fully understood although binary interactions are likely an important formation pathway. To better constrain the history of MAE stars, we have measured metallicities, carbon isotope ratios, and CNO abundances in 49 red clump stars and four red giants. Our sample included thin disk, thick disk, and MAE stars to best constrain different formation scenarios. We observed our sample stars using the Tull spectrograph on the McDonald 2.7m telescope and derived abundances using synthetic spectra created by the Turbospectrum radiative transfer code. Overall, we find that 10 of our red clump MAE stars are consistent with the average thick disk carbon isotope ratio of $^{12}$C/$^{13}$C = 8.2 $\pm$ 3.4. We find five MAE stars that have significantly higher carbon isotope ratios ($^{12}$C/$^{13}$C $>$ 15) than stars at similar metallicities. Two of the anomalous stars show abundance patterns consistent with AGB mass transfer while three MAE stars have $^{12}$C/$^{13}$C ratios similar to massive, single RC stars and show no clear signs of binarity from radial velocity monitoring or from the Gaia RUWE measurement. Overall, we find that carbon isotope ratio measurements provide a unique constraint when discerning the possible origins of MAE stars.

The baseline contrast floor from the Roman Coronagraph's High-Order Wavefront Sensing and Control strategy likely degrades over the course of time, requiring periodic recalibration of the dark hole. Here, we propose to consider spatial linear dark field control (sLDFC) on a one-sided deep-contrast region of the focal plane as a potential test. Implementing sLDFC on CGI will likely require some unique data acquisition strategies given the EMCCD's high flux sensitivity in long exposures/high gain: we outline three possible approaches. However, if successful, sLDFC's advances are substantial: (1) enabling us to maintain a fainter, more temporally correlated dark hole which will improve CGI's contrast after post-processing and (2) efficiently providing a reliable signal (bright field) for accurate reconstruction of residual starlight in the dark field, further boosting CGI's detection limit for bright targets.

Radial differential rotation is an important factor in stellar dynamo theory. In the Sun, helioseismology has revealed a near surface shear layer in the upper 5 to 10% of the convection zone. At low to midlatitudes, the rotation velocity gradient decreases sharply near the surface. A depth gradient in rotational velocity was recently detected in the low photosphere using a differential interferometric method on spectroscopic data. Granular structures at different depths in the FeI 630.15 nm line showed a systematic retrograde shift compared to continuum structures, suggesting a height-related decrease in angular velocity, dependent on the assumed granulation coherence time. We use a more direct approach to measure the differential rotational velocity at different photospheric heights. We performed spectroscopic scans of the same granular region in FeI 630.15 nm and CaI 616.2 nm lines and measured displacements of images at different line chords between consecutive scans. These observations require excellent seeing, stable adaptive optics, and scanning times shorter than the granulation lifetime. Adaptive optics stabilizes continuum images but not higher-altitude rotation differences. We used THEMIS and HINODE SOT FeI data to measure formation height differences via perspective shifts observed away from the disk center with the slit radially oriented. Measurements at disk center and $\pm$25° latitude along the central meridian show a parabolic decrease in rotational velocity with height, reaching about 16% slower rotation at 80 km above the continuum. No significant difference is found between the equator and $\pm$25° latitudes. The low photosphere is a transition zone between the convective and radiative layers. Our measurements provide new constraints on its dynamical behavior and valuable boundary conditions for numerical simulations of the Sun s upper convection zone.

Nova explosions occur on accreting white dwarfs. A thermonuclear runaway in the H-rich accreted envelope causes its ejection without destroying the white dwarf, and an increase in the luminosity by several magnitudes. Accretion is re-established some time after the explosion. The explosion of the nova itself is expected to affect the mass-transfer rate from the secondary and the accretion rate, but these effects have been little explored observationally. Most novae are observed only in outburst and the properties of the host systems are unknown. X-ray observations of novae happen mostly during outburst; only a few have been the target of dedicated X-ray observations years or decades after outburst. However, the X-ray emission long after the outburst provides a powerful diagnostics of the accretion rate and the possible magnetic nature of the white dwarf. We have explored the first two years of the SRG/eROSITA All Sky Survey (eRASS) for X-ray sources correlated with Galactic historical novae. We present the first population study of nova hosting systems in X-rays, focus on the evolution of accretion rate as a function of the time since last outburst, and look for new candidates for magnetic systems. In total, 32 X-ray counterparts of novae are found in the western Galactic hemisphere. Combined with 53 nova detections published for the eastern hemisphere, the fraction of X-ray detected novae in quiescence is 18% of the Galactic novae. We have, for the first time, enough statistics to observationally determine the evolution of accretion rate as a function of time since the last nova outburst for a time span of 120 years. The results confirm that magnetic systems remain systematically at higher fluxes and that accretion is enhanced during the first years after the outburst, as predicted theoretically. We also identify new IP candidates in AT Cnc and RR Cha.

Modules for Experiments in Stellar Astrophysics (MESA; (Paxton2011, 2013, 2015, 2018, 2019, Jermyn 2023) is a widely used open-source software for modelling stellar evolution. In many studies, the computational grids covers thousands of models (e.g., Joyce et al 2024), which requires a tremendous amount of time and computational effort to process. To streamline data analysis, the mesalab package was developed. This Python-based pipeline is designed to simplify the post-processing of MESA outputs by automatically identifying various stellar evolutionary phases, with a specific focus on the ``blue loop'' -- a blue-ward excursion in the Hertzsprung-Russell Diagram (HRD) for intermediate-mass stars often associated with peculiar pulsational phenomena like ``strange modes.''

Precessing relativistic jets launched by compact objects are rarely directly measured, and present an invaluable opportunity to better understand many features of astrophysical jets. In this Letter we present MeerKAT radio observations of the neutron star X-ray binary system (NSXB) Circinus X-1 (Cir X-1). We observe a curved S-shaped morphology on $\sim 20''$ $(\sim1\:\text{pc})$ scales in the radio emission around Cir X-1. We identify flux density and position changes in the S-shaped emission on year timescales, robustly showing its association with relativistic jets. The jets of Cir X-1 are still propagating with mildly relativistic velocities $\sim1\:\text{pc}$ from the core, the first time such large scale jets have been seen from a NSXB. The position angle of the jet axis is observed to vary on year timescales, over an extreme range of at least $110°$. The morphology and position angle changes of the jet are best explained by a smoothly changing launch direction, verifying suggestions from previous literature, and indicating that precession of the jets is occurring. Steady precession of the jet is one interpretation of the data, and if occurring, we constrain the precession period and half-opening angle to $>10$ years and $>33°$ respectively, indicating precession in a different parameter space to similar known objects such as SS~433.

The goal of the present paper is to assess the usefulness of quasars as cosmological distance indicators. We calibrate, in a model-independent way, the non-linear relation between X-ray and UV emission to derive quasar luminosity distances. Using this calibration, we construct the Hubble diagram up to redshift $z \sim 7.5$, and test the ability of quasars to constrain the $\Lambda$CDM and $\omega$CDM models at low and high redshifts, in combination with the Type Ia supernova Pantheon Plus sample, as well as the latest results from DESI and Planck Compressed data. We find consistency with previous studies in the values of $\gamma$ and $\beta$ when using a quasar subsample with $z<1.43$. However, when the data is used to constrain the cosmological models, we find that quasars fail as reliable distance indicators, even when combined with other independent cosmological datasets.

Snowball episodes are associated with increases in atmospheric oxygen and the complexity of life on Earth, and they may be essential for the development of complex life on exoplanets. Sustained stable Snowball episodes require a Snowball bifurcation and climate bistability between the globally ice-covered Snowball state and a state with at least some open ocean. We find that climate bistability disappears for an aquaplanet with a slab ocean in the Global Climate Model ExoCAM when the orbital eccentricity is increased to 0.2-0.3. This happens because the Snowball state loses stability as seasonal insolation variations intensify, while the warm state remains stable due to the ocean's large heat capacity. We use a low-order ice-thermodynamic model to show that the Snowball state loses stability as seasonality increases because winter freezing at the ice bottom is reduced relative to summer melting at the ice top due to ice self-insulation. Combined with previous research showing that Snowball climate bistability diminishes for planets orbiting low-mass stars, ones with longer rotation periods, and disappears entirely for tidally locked planets, our work suggests that the Snowball climate bistability may not be as robust to planetary parameters as previously thought, representing one aspect of habitability more consistent with the Rare Earth Hypothesis than the Copernican Principle.

We assess archival high-energy data for key stars on the Habitable Worlds Observatory (HWO) Target Stars and Systems 2025 list, as stellar radiation is critical to shaping and interpreting planetary atmospheres. Using a sample of 98 nearby stars (HWO Tier 1 targets), we compile and evaluate X-ray and ultraviolet (UV) data from archival eROSITA, Chandra, XMM-Newton, ROSAT, EUVE, Swift, FUSE, IUE, GALEX, and HST. We examine spectral and temporal coverage, assess data quality, and identify major gaps. UV data are moderately available, with most coverage coming from near-UV spectra from IUE. Far fewer stars have far-UV spectra, especially from HST. In the X-ray regime, some stars have high-quality spectra, while others are limited to shallow detections or broad-band photometry. A small fraction of the sample has both X-ray and UV spectra of sufficient quality to support full spectral energy distribution modeling. Truly comprehensive coverage across X-ray, extreme-UV, and both UV bands remains extremely rare. Most datasets are single-epoch, limiting assessments of variability and flares - key factors in atmospheric photochemistry and escape. Moreover, the lack of simultaneous or contemporaneous observations across bands adds further uncertainty. Our findings underscore the need for new space-based missions and coordinated multiwavelength campaigns, ideally with overlapping coverage, to improve stellar characterization for HWO. As several key observatories age and face potential decommissioning, there is a narrow window of opportunity to secure these critical data. Investing in this effort now will directly support the science goals of HWO and enhance future studies of planetary habitability.

The typical galaxy in the local universe is expected to be in a self-regulated quasi-equilibrium, displaying a settled morphology that falls within the Hubble Sequence. The Arp and Arp-Madore catalogs are filled with striking examples of galaxies that defy these expectations, making them useful targets for studying the astrophysics that controls dramatic, but short-lived, episodes of disequilibrium that mark galaxies' evolution. In this paper, we greatly expand the available Hubble Space Telescope (HST) imaging of galaxies drawn from the Arp and Arp-Madore catalogs. We present new optical F606W images and point-source photometry for 216 systems, whose sizes are well-matched to the Advanced Camera for Surveys' (ACS) wide field of view. Essentially none of the sample had been previously observed with Hubble. The resulting images display rich morphologies, revealing a variety of massive stars, HII regions, stellar clusters, dust lanes, tidal tails, backlit galaxies, and occasional chance superpositions. We provide a pedagogical guide for interpreting highly-resolved optical galaxy images, which has general application beyond this atlas. The atlas images also provide a superb starting point for more detailed studies with high-resolution imaging in other wavelengths, and spectroscopy to track kinematics and the interstellar medium (ISM). Areas of obvious scientific relevance include feedback and star formation in merging and interacting galaxies, resolved stellar populations at the extremes of stellar density, the properties of young massive stars and stellar clusters, the physics of the cold ISM and dust, and stellar and gas dynamics.

Dark matter mass density profiles and velocity distributions for a set of dwarf spheroidal galaxies (dSphs) have recently been obtained (this http URL, this http URL, this http URL arXiv:2412.20585) by performing a multi-parametric fit to the stellar observations with the help of the GravSphere which solves the Jeans equation. We use these results to calculate the geometrical factors for estimation of the fluxes of cosmic rays expected from decay ($D$-factor) and annihilation ($J_s$-, $J_p$- and $J_d$-factors for $s$-, $p$- and $d$- wave processes) of dark matter particles in galaxies. The general novelty is the account for a possible anisotropy in velocities of dark matter particles. On the basis of this analysis we present empirical scaling approximations to these factors as functions of typical observables: distance to the galaxy $d$, it's half-radius $r_h$ and line-of-sight stellar velocity dispersion $\sigma_{LOS}$. They can be applied to any galaxy, and for $D$- and $J_s$-factors we refine the estimates of (Andrew B. Pace, Louis E. Strigari arXiv:1802.06811): the shifts in the central values remain within 1-2$\sigma$ error bars.

The extended solar corona at 10-30 solar radii is essentially devoid of all waves below 100 kHz other than triggered ion acoustic waves (TIAW), which consist of a low frequency electromagnetic wave at a frequency of a few Hz coupled to one or more electrostatic waves at a few hundred Hz, such that the amplitudes of the higher frequency waves peak at a fixed phase of each low frequency wave period. All the waves in a TIAW event travel at the same phase speed, which is found to be 150 km s (the ion acoustic speed was about 100 km s). It has not been possible to explain the TIAW as a resonant wave-wave interaction, so a non-resonant interaction has been considered in which the loss of energy by the low frequency wave is used to both heat the electrons and grow the higher frequency waves. Evidence in support of this explanation is described and a PIC simulation that discusses this process is summarized. This interplay between a pair of waves, mediated by modifications of plasma parameters and energy conversion, represents a significant nonlinear process in plasma physics, the study of which will deepen the understanding of energy transfer, wave generation, and plasma dynamics in diverse astrophysical environments.

Determining the spatial distribution of Galactic cosmic rays (CRs) is fundamental to understand how these particles propagate in interstellar space and to infer their source spectra. The most sensitive method of studying this problem is observing the gamma--ray and neutrino diffuse fluxes produced in the inelastic interactions of CR protons and nuclei with interstellar gas that encode the energy and spatial distributions of the interacting particles. Most theoretical models assume that the spatial and energy dependencies of the spectra of CR protons and nuclei are factorized, so that the CR energy spectra have the same shape at all points of the Milky Way. However, on the base of the Fermi--LAT observations, some authors have tentatively inferred that particles in the inner regions of the Milky Way have significantly harder spectra. Recently the ground--based telescopes Tibet--AS$\gamma$, LHAASO and HAWC have extended the measurements of the diffuse gamma--ray flux to much higher energies, and the IceCube neutrino telescope has obtained evidence for a $\nu$ flux from the Galactic disk. These new measurements allow for new tests of the factorization hypothesis, and disfavor models in which the spectral shape of CR protons and nuclei is harder in the inner part of the Galaxy. In fact, the LHAASO data above 30~TeV are in better agreement with the opposite hypothesis, that cosmic rays is the inner Galaxy have softer spectra than those in the outer Galaxy. This result, that has currently only modest significance due to large statistical and systematic uncertainties, could be explained by assuming that the CR confinement volume increases with energy. An alternative explanation is that a significant fraction of cosmic rays in the multi--PeV energy range is of extragalactic origin.

Lyman-$\alpha$ (Ly$\alpha$) is the most prominent ultraviolet emission line in low-mass stars, playing a crucial role in exoplanet atmospheric photochemistry, heating, and escape. However, interstellar medium (ISM) absorption typically obscures most of the Ly$\alpha$ profile, requiring reconstructions that introduce systematic uncertainties. We present intrinsic Ly$\alpha$ profiles for 12 high radial velocity K and M dwarfs, where Doppler shifting minimizes ISM contamination, allowing direct measurements of $\sim$50-95\% of the line flux. Our sample spans the K-to-M spectral transition, enabling us to constrain the dependence of self-reversals in Ly$\alpha$ emission profiles on effective temperature ($T_{eff}$). The depth of self-reversal, driven by non-local thermodynamic equilibrium (LTE) effects, decreases with decreasing $T_{eff}$, with M dwarfs exhibiting little to none. Two stars, Ross 1044 and Ross 451, were observed over multiple days, revealing $\sim$20\% Ly$\alpha$ variability confined to the line core - implying that studies relying on reconstructions may underestimate temporal variability. We find strong correlations between Ly$\alpha$ flux, peak-to-trough ratio, and hydrogen departure coefficients with $T_{eff}$, providing empirical constraints for stellar atmosphere models. A comparison of Ly$\alpha$ flux in the habitable zone shows measured values for high radial velocity stars less than the reconstructed values for the rest of the sample, likely due to the older ages of the high-RV stars and/or overestimated reconstructed fluxes due to model deficiency (e.g., neglecting self-reversal). Our results establish an empirical foundation for Ly$\alpha$ emission in K and M dwarfs, reducing uncertainties in reconstructions and improving models of stellar UV emission relevant to exoplanetary studies.

The Hubble tension, a persistent discrepancy between early and late Universe measurements of $H_0$, poses a significant challenge to the standard cosmological model. In this work, we present a new Bayesian hierarchical framework designed to meticulously decompose this observed tension into its constituent parts: standard measurement errors, information loss arising from parameter-space projection, and genuine physical tension. Our approach, employing Fisher matrix analysis with MCMC-estimated loss coefficients and explicitly modeling information loss via variance inflation factors ($\lambda$), is particularly important in high-precision analysis where even seemingly small information losses can impact conclusions. We find that the real tension component ($T_{real}$) has a mean value of 5.94 km/s/Mpc (95\% CI: [3.32, 8.64] km/s/Mpc). Quantitatively, approximately 78\% of the observed tension variance is attributed to real tension, 13\% to measurement error, and 9\% to information loss. Despite this, our decomposition indicates that the observed $\sim$$6.39\sigma$ discrepancy is predominantly a real physical phenomenon, with real tension contributing $\sim$$5.64\sigma$. Our findings strongly suggest that the Hubble tension is robust and probably points toward new physics beyond the $\Lambda$CDM model.

The generation of dense electron-positron pair beams in the laboratory can enable direct tests of theoretical models of $\gamma$-ray bursts and active galactic nuclei. We have successfully achieved this using ultra-relativistic protons accelerated by the Super Proton Synchrotron at CERN. In the first application of this experimental platform, the stability of the pair beam is studied as it propagates through a metre-length plasma, analogous to TeV $\gamma$-ray induced pair cascades in the intergalactic medium. It has been argued that pair beam instabilities disrupt the cascade, thus accounting for the observed lack of reprocessed GeV emission from TeV blazars. If true this would remove the need for a moderate strength intergalactic magnetic field to explain the observations. We find that the pair beam instability is suppressed if the beam is not perfectly collimated or monochromatic, hence the lower limit to the intergalactic magnetic field inferred from $\gamma$-ray observations of blazars is robust.

Recent observations of long period radio transients, such as GLEAM-X J0704-37 and ILTJ1101 + 5521, have revealed a previously unrecognized population of galactic radio transient sources associated with white dwarf - M dwarf binaries. It is an open question how to produce coherent radio emission in these systems, though a model driven by binary interaction seems likely given the nature and correlation of the emission with the binaries' orbital period. Using kinetic plasma simulations, we demonstrate that the relativistic electron cyclotron maser instability (ECMI) is a viable mechanism for generating radio pulses in white dwarf - M dwarf systems, akin to planetary radio emission, such as that from the Jupiter-Io system. We quantify the relativistic ECMI in the nonlinear regime under conditions relevant for white dwarf radio emission for the first time. Our simulations demonstrate that the ECMI can intrinsically produce partially linearly polarized emission relevant to explaining the observed emission spectrum of the two galactic sources, though the precise details will depend on the plasma composition. Our work paves the way for a systematic and fully nonlinear computational modeling of radio emission from interacting white dwarf sources.

The analytical model found that the intrinsic variation in the initial metallicity of the Type Ia supernova (SN Ia) progenitor stars ($Z_{progenitor}$) translates into a 25% variation in the $^{56}$Ni mass synthesized and, therefore, 0.2 mag difference in the observed peak luminosity of SNe Ia. Previous observational studies used the currently-observed global gas-phase metallicity of host galaxies, instead of $Z_{progenitor}$ used in the model, and showed a higher scatter in the $^{56}$Ni mass measurements compared to the model prediction. Here, we use $Z_{progenitor}$ of 34 normal SNe Ia and employ recent SN Ia explosion models with various configurations to cover the observed $^{56}$Ni mass range. Unlike previous studies, our sample covers the $Z_{progenitor}$ range, where most of the $Z_{progenitor}$ effect occurs. Linear regression returns a slope of 0.02+-0.03, which is the opposite trend to the analytical model, but at at low statistical significance level. We find that comparing our sample with SN Ia explosion models on the $Z_{progenitor}$--$^{56}$Ni mass diagram allows us to constrain the progenitor scenarios. We also explore other chemical composition indicators. For $(Fe/H)_{progenitor}$, our sample follows the trend predicted by the analytical models, but at a low significance level. Noticeably, $(\alpha/Fe)_{progenitor}$ shows the opposite trend and a clear gap. When we split the sample at $(\alpha/Fe)_{progenitor}$ = 0.35 $(\alpha/Fe)_{\odot}$, we find a 3$\sigma$ difference in the weighted-means of the $^{56}$Ni mass. Lastly, SNe Ia in different $Z_{progenitor}$ groups show a difference of 0.14+-0.09 mag in the standardized luminosity. The present work highlights a holistic approach (from the progenitor star to the explosion with SN Ia and host galaxy observational data) to understand the underlying physics of SNe Ia for more accurate and precise cosmology.

Within the standard six-parameter Lambda cold dark matter ($\Lambda$CDM) model, a $2$-$3\sigma$ tension persists between baryon acoustic oscillation (BAO) measurements from the Dark Energy Spectroscopic Instrument (DESI) and observations of the cosmic microwave background (CMB). Although this tension has often been interpreted as evidence for dynamical dark energy or a sum of neutrino masses below the established minimum, recent studies suggest it may instead originate from an underestimation of the reionization optical depth, particularly when inferred from large-scale CMB polarization. Jhaveri et al. propose that a suppression of large-scale primordial curvature power could partially cancel the contribution of $\tau$ to the CMB low-$\ell$ polarization power spectrum, leading to a biased low $\tau$ measurement in standard analyses. In this work, we investigate whether punctuated inflation - which generates a suppression of primordial power on large scales through a transient fast-roll phase - can raise the inferred $\tau$ value and thereby reconcile the consistency between CMB and BAO. For simple models with step-like features in the inflaton potential, we find that the constraint on $\tau$ and the CMB-BAO tension remain nearly identical to those in the standard six-parameter $\Lambda$CDM model. We provide a physical explanation for this negative result.

The rapidly expanding catalog of gravitational-wave detections provides a powerful probe of the formation history of compact binaries across cosmic time. In this work, we extend the Binary Compact Object (BCO) phase-space framework to the full set of events in the GWTC-4 catalog to map the observed binary formation scenarios in a data-driven way. Applying this framework, we identify distinct regions of phase-space associated with different channels and discover for the first time a unique mass-cutoff scale in a data-driven way. The mapping of these on different formation channels reveals a population of first-generation (1G) black holes sharply truncated at approximately 45.5 $M_\odot$, consistent with the theoretically predicted pair-instability supernova (PISN) mass gap. These findings demonstrate the capability of the BCO phase-space to disentangle overlapping formation pathways, establish robust connections between gravitational-wave observations and binary evolution, and highlight the potential of upcoming observing runs to reveal rare populations and exotic origins.

With the discovery of pulsations in some of the ultra-luminous X-ray sources (ULXs), it is quite clear that most of the ULXs harbor either a neutron star or a stellar mass black hole as a compact object accreting at super-Eddington rates. In spite of having such a high accretion rate, the reprocessed emission in the ULX sources is quite meagre compared to that observed in Galactic X-ray binaries, except for some absorption lines in the winds. In this work, we investigate the extent of reprocessed emission in ULXs using three diagnostics: (i) searches for Fe $\rm K\alpha$ lines in bright well-known ULXs and Ultra luminous X-ray Pulsars (ULXPs), (ii) evolution of hardness ratio around the eclipse transitions in the eclipsing ULXs, and (iii) the flux ratio between eclipse and out-of-eclipse (OOE) phases in eclipsing ULXs. We placed the most stringent constraints to date on the upper limits on EW of the iron line, 11--20 eV. Furthermore, we have not observed any significant changes in the hardness ratio during the ingress or egress, while in Galactic eclipsing X-ray binaries, an increase in the hardness ratio is observed during the transitions. Finally, the reprocessing efficiency (eclipse to OOE flux ratio) is found to be larger in ULXs compared to Galactic eclipsing X-ray binaries. Based on these results, we discuss the possibility of a metal-poor or highly ionized environment surrounding the ULXs, which suppresses reprocessed emission features.

The LIGO-Virgo-KAGRA Collaboration recently reported GW231123, a black hole merger with total mass of around 190-265 solar mass. This event adds to the growing evidence of "lite" intermediate mass black hole (IMBH) discoveries of post-merger black holes >100 solar mass. GW231123 posed several data analysis challenges owing to waveform-model systematics and presence of noise artifacts called glitches. We present the first comprehensive machine learning analysis to further validate this event, strengthen its astrophysical inference, and characterize instrumental noise in its vicinity. Our approach uses a combination of tools tailored for specific analyses: GW-Whisper, an adaptation of OpenAI's audio transformer, ArchGEM, a Gaussian mixture model-based soft clustering and density approximation software and AWaRe, a convolutional autoencoder. We identify the data segment containing the merger with >70% confidence in both detectors and verify its astrophysical origin. We then characterize the scattered light glitch around the event, providing the first physically interpretable parameters for the glitch. We also reconstruct the real waveforms from the data with slightly better agreement to model-agnostic reconstructions than to quasi-circular models, hinting at possible astrophysics beyond current waveform families (such as non-circular orbits or environmental imprints). Finally, by demonstrating high-fidelity waveform reconstructions for simulated mergers with total masses between 100-1000 solar mass, we show that our method can confidently probe the IMBH regime. Our integrated framework offers a powerful complementary tool to traditional pipelines for rapid, robust analysis of massive, glitch-contaminated events.

We investigate a class of interacting dark energy and dark matter (DM) models, where dark energy is modeled as a $k$-essence scalar field with an inverse-square potential. Two general forms of interaction are considered: one proportional to the Hubble parameter, and another independent of the Hubble parameter, depending instead on combinations of the energy densities and pressures of the dark sectors. The dynamics are analyzed using a dynamical system stability framework by constructing an autonomous system of equations. The models are tested against a wide range of observational datasets, including cosmic chronometers (CC), BAO measurements from DESI DR2, compressed Planck data (PLA), Pantheon+ (PP), DES supernovae, Big Bang Nucleosynthesis (BBN), and strong lensing data from H0LiCOW (HCW). The analysis shows that the models consistently reproduce all major cosmological epochs and yield statistically competitive results compared to the flat $\Lambda$CDM model. The models exhibit stable late-time de-Sitter solutions, ensuring ghost-free evolution, with the Hubble constant in the range $H_0 \sim 67$--$70$ km/s/Mpc.

We propose a new amorphous dust emission model based on the soft-potential (SP) model, applicable in the long-wavelength range from the far-infrared to the microwave. The SP model is widely accepted in material physics to explain amorphous thermal properties and is an extension of the two-level systems (TLS) model, which has been applied in interstellar amorphous dust physics. In the SP model, by assuming that some atoms composing amorphous dust are trapped in a double-well potential (DWP) described by a quartic function, the electric interaction can be solved directly, allowing the absorption cross-section of the amorphous dust to be calculated. We present numerical and analytical solutions for the absorption cross-section of amorphous dust and compare these results, finding good agreement for the DWP with a sufficiently high potential barrier. Our findings show that the SP model can reproduce the observed spectrum of the Perseus molecular cloud with slightly better accuracy than the conventional TLS model. Additionally, the SP model can more effectively explain various long-wavelength dust emission features, such as the spectrum flattening in the submillimeter range and anomalous microwave emission (AME), compared to the TLS model. Further comparison and verification with observational data and laboratory measurements are necessary to refine the model.

Ref. [Symmetry 15 (2023) 709] introduced a Brans-Dicke-like framework wherein the scalar field $\phi$ is composed of both $G$ and $c$ which, for this reason, co-vary according to $c^{3}/G=\text{constant}$. In this paper, we use observational data to constrain the supposed co-varying $G$ and $c$. The datasets include SN Ia, BAO and the value of $\theta$ extracted from CMB data. A proxy function is demanded for the varying $c$ since the framework does not provide a closed set of equations for computing the functional form of either $G$ or $c$ uniquely. Accordingly, we choose three separate parameterizations for $c\left(z\right)$ inspired both by desirable properties of the varying speed of light (VSL) and by successful phenomenological models from the literature -- including the one by Gupta (CCC framework in e.g. Ref. [Mon. Not. R. Astron. Soc., 498 (2020) 4481-4491]. When combined with DESI, Pantheon+ data strongly favor a variable speed of light with more than $3\sigma$ confidence level for all parameterizations considered in this paper, whereas Union2.1 suggests no variation of the speed of light. As we shall demonstrate, this apparent discrepancy is due to a strong correlation that emerges between $H_0$ and VSL.

Long-period radio transients (LPTs) are a new class of coherent radio sources with periods ranging from minutes to hours. Recently, two LPT sources, ILT J1101+5521 and GLEAM-X J0704-37, with periods of 2-3 hours has been confirmed to originate from white dwarf (WD) -- M dwarf (MD) binaries. In this work, we propose that at least some LPTs originate from the magnetic WD -- MD binaries in the pre-polar phase. The asynchronism between the WD's rotation and the binary's orbital motion allows for the unipolar-inductor mechanism or magnetosphere interaction to operate and accelerate radiating particles, with the dominant process depending on the magnetic moment ratio of the two stars. Under asynchronism condition, both the peak flux and the polarization of radio pulses will be modulated by the beat period. The pre-polar phase characterized by an extremely low accretion rate provides the relatively clean magnetospheric environment necessary for a loss-cone-driven maser (LCDM) mechanism to operate, producing the LPT emission. The observed pulse duty cycle of $10^{-3}-10^{-1}$ is attributed to a beaming effect modulated by the binary's magnetic geometry. Furthermore, the magnetized environment of a WD--MD binary is conducive to Faraday conversion with weak coupling, which implies that the polarization state of LPTs should vary significantly at different periods. Finally, we predict that LPTs from WD--MD binaries should exhibit a period distribution following $f_P(P)dP \propto P^{-(1.67-2.33)}dP$ and a luminosity function described by $f_L(L)dL \propto L^{-(1.80-2.67)}dL$, which can be tested by the future large sample.

Deviations from the so-called {\it cosmic distance duality relation} may result from systematic errors in distance measurements or, more interestingly, hint at new physics. Further, it can also be related to the Hubble constant tension between early and local measurements of $H_0$. Based on this, we test validity of this relation through a model-independent parameterization of the Hubble rate via the well-estabilished Bézier polynomials approach. We seek for possible departures from the relation considering three parametrizations, i) a power-law correction, ii) a logarithmic correction and iii) a Padé series $P_{n,m}(z)$ of order (1;2) with $n=1$ being the order of the numerator while $m=2$ is the order of the denominator. Then, assuming a flat scenario, we test them through Monte Carlo -- Markov chain analyses that combine low- and intermediate/high-$z$ data sets, such as observational Hubble data, the Pantheon catalog of type Ia supernovae, galaxy clusters, the second data release from the DESI Collaboration and gamma-ray bursts. In particular, we distinguish between \emph{Analysis A} and \emph{Analysis C}, depending whether the prompt emission $E_{iso}-E_p$ or the prompt-afterglow $L_0-E_p-T$ gamma-ray burst correlations, respectively, is fit together with the other probes previously described. Our results seem to point towards a \emph{no violation} of the cosmic distance duality relation and a preference towards Planck's value of $H_0$.

This research presents unsupervised machine learning and statistical techniques to detect and analyze tidal tails in five open clusters, BH 164, Alessi 2, NGC 2281, NGC 2354, and M67, with ages ranging from 60 Myr to 4 Gyr, using Gaia DR3 data. Color-magnitude diagram (CMD) matching, principal component analysis (PCA), and density-based clustering were used to detect the tails. Tidal tails were found in all clusters, extending 40 to 100 pc and containing 100 to 200 stars. Other statistical methods were used to study the stellar kinematics, photometry, and spatial distribution of the detected features. The tails had luminosity functions similar to their parent clusters, generally lacked massive stars, and showed higher binary fractions. Significant rotation was detected in M67 and NGC 2281 for the first time.

Naked singularities (NkS) are solutions to the Einstein field equations that violate the cosmic censorship conjecture. Recent studies indicate that these objects may serve as compelling mimickers of black hole shadows. In this work, we investigate the accretion dynamics of selected time-like naked singularities using general relativistic magnetohydrodynamic (GRMHD) simulations. Our objective is to determine whether naked singularities exhibit distinct signatures compared to black holes. We find that, unlike black holes, naked singularities exhibit a centrifugal barrier that prevents direct accretion of rotating matter onto the NkS. Despite reduced magnetization in the funnel region, these objects are capable of generating jet powers comparable to those observed in black holes. Additionally, we observe that accreting matter releases gravitational energy as it is driven towards the NkS, powering the strong outflow via local fluid pressure gradient or magnetic pressure forces.

An international team of researchers has come together to undertake an ultra-broadband exploration of the deep radio sky. The superMIGHTEE project combines data from the MIGHTEE project, using the precursor Square Kilometre Array (SKA) MeerKAT telescope in South Africa, with observations from the upgraded Giant Metrewave Radio Telescope (uGMRT) in India to produce deep images at several $\mu$Jy sensitivity over a frequency range of 200 MHz--2.5 GHz, with an angular resolution of a few arcseconds. This paper describes the initial superMIGHTEE uGMRT data release, comprising total intensity continuum images covering a total of 9.9 deg$^2$ at 650 MHz and 6.9 deg$^2$ at 400 MHz in the XMM-LSS, COSMOS, and E-CDFS deep fields. The associated radio source catalogs include 27,101 sources at 650 MHz and 10,946 sources at 400 MHz. The redshift distribution of the sources extends to $z\sim4$ with a median value of $z=1$. An overview of the broadband spectra of the sources, in combination with the MeerKAT MIGHTEE 1280 MHz data, reveals a clear change in spectral properties at the transition from an active galactic nuclei-dominated population to a population dominated by star-forming galaxies at flux densities of a few mJy. At higher frequencies, the star-forming galaxy population exhibits an optically thin synchrotron spectral index indicative of energy injection from supernovae. At lower frequencies, the spectra flatten significantly with decreasing flux density, and the fraction of sources with peaked spectra increases. This is the first superMIGHTEE uGMRT data release. Subsequent releases will include spectropolarimetric and spectral line image cubes, as well as images at lower frequencies. The goal of the superMIGHTEE ultra-wideband dataset is to enhance our understanding of the evolution of active galactic nuclei and star-forming galaxies over cosmic time, (abridged).

It is often argued that gravitational instability of realistic protoplanetary discs is only possible at distances larger than $\sim 50$ au from the central star, requiring high disc masses and accretion rates, and that therefore disc fragmentation results in the production of brown dwarfs rather than gas giant planets. However, the effects of dust growth on opacity can be very significant but have not been taken into account systematically in the models of fragmenting discs. We employ dust opacity that depends on both temperature and maximum grain size to evaluate analytically the properties of a critically fragmenting protoplanetary disc. We find that dust growth may promote disc fragmentation at disc radii as small as $\sim 30$ au. As a result, the critical disc masses and accretion rates are smaller, and the initial fragment masses are in the gas giant planet mass regime. While this suggests that formation of gas giant planets by disc fragmentation may be more likely than usually believed, we caution that numerical models of the process are needed to evaluate the effects not taken into account here, e.g., dust grain mobility and fragment evolution after disc fragmentation.

The disc around the pre-main-sequence triple star system GW Orionis is known from observations to be warped and broken. Theoretical modelling has produced conflicting results regarding the mechanism responsible for breaking the disc. Analytical predictions for the measured parameters of GW Ori suggest the disc is only marginally stable to tearing. We present new high-resolution simulations of GW Ori that replicate the wavelike regime expected in thick, low-turbulence protoplanetary discs for the first time to settle this question. Using the most optimistic values of misalignment and stellar mass ratio allowed by observational constraints, we find that the GW Ori disc can be torn by stellar torques alone, without need for an embedded planet. Even if the disc retains a smooth warp in simulations with similar parameters, it is likely that any small perturbation in the density or temperature structure could cause the disc to break. The new simulations rule out retrograde disc rotation relative to the stellar orbits and tentatively suggest the thicker ($h/r = 0.04$) disc better matches observations. Going forward, we should take care to ensure models of GW Ori and similar systems appropriately represent the propagation of warps. Additionally, analytical predictions are derived from idealized (and often massless) discs and it is useful to assess how each observed disc might deviate from those assumptions, especially in the context of a young and active star-forming neighbourhood.

The origin of ultra-high-energy cosmic rays (UHECRs) is one of the intriguing mysteries in astroparticle physics. In order to identify their sources, we need precise knowledge of the mass composition of UHECRs. The direct detection of UHECRs is not feasible at energies above 0.1 PeV, necessitating the use of mass-sensitive observables of extended air showers induced by UHECRs interacting with the atmosphere. One way to achieve high statistics for these mass-sensitive observables is to use ground-based detector arrays, such as the Surface Detector (SD) of the Pierre Auger Observatory. The SD consists of three sub-arrays of independent detector stations arranged in triangular grids with different spacings. Recently, it has been shown that neural networks (NNs) can extract mass-sensitive observables from data taken by the SD-1500, the largest sub-detector of the SD. In this contribution, we demonstrate the feasibility of using NNs to reconstruct a high-level shower observable, the depth of the shower maximum, from data simulated for and observed by the SD-750, the second-largest detector array nested within the SD-1500. A simulation study shows that the SD-750 NN exhibits behavior similar to that of an SD-1500 NN and outperforms the latter in the energy range [1, 10) EeV. Moreover, we show that, after performing a correction and calibration procedure, the predictions of the SD-750 NN are consistent with the measurement of the depth of the shower maximum obtained by the Fluorescence Detector of the Pierre Auger Observatory.

We present monitoring observations of quasar microlensing in the quadruple quasar HE0435-1223. The microlensing-induced light curves of the quasar images are chromatic, i.e. they depend on the applied filter band. Comparison with microlensing simulations allows us to infer properties of the accretion disk. We determine the R and V band light curves of the four images of HE0435-1223 from 79 and 80 epochs respectively, taken from 2014 to 2024 at the Las Cumbres Observatory using difference imaging analysis. We consider difference light curves to remove the intrinsic quasar variability. This reveals a prominent long-term chromatic microlensing event in image B. We use microlensing light curve simulations with both Gaussian and standard thin accretion disk brightness profiles to analyse this signal. The particularly strong signal observed in image B of HE0435-1223 makes it possible to detect a size ratio of the accretion disk in the R band compared to the V band of $1.24^{+0.08}_{-0.20}$ and $1.42^{+0.11}_{-0.22}$ for the Gaussian and the thin disk model, respectively. These values are in agreement with standard thin disk theory. For the absolute size we find large disk half-light radii of around 0.7 to 1.0 Einstein radii with an uncertainty of about 0.6 dex (depending on the filter bands and the models). Finally, our calculations show that image B undergoes caustic crossings about once per year.

Ultra-compact double white dwarfs (DWDs) represent key targets for multi-messenger astrophysics, as they may be observed both through gravitational waves and the electromagnetic (EM) spectrum. The future Laser Interferometer Space Antenna (LISA) will detect thousands of these systems, and they are predicted to be the most numerous science targets of the mission. We develop a strategy to identify LISA source candidates in multiband photometric surveys. We constructed a synthetic EM catalogue of white dwarf (WD) detections based on a population synthesis code combined with a semi-analytical model of the Milky Way and a consistent cooling model for the evolution. We compared sources in the LISA band with other WD observations in magnitude-colour and colour-colour plots. From a full sky survey with $u \le$24.5, we find that 57$\%$ of the sources in the LISA band occupy a specific region in colour-colour diagrams. Inside this area, we find that $\sim 63\%$ (6.5 $\times 10^4$) of EM observations are LISA candidates, $\sim 31\%$ ($ 3.2 \times 10^4$) are DWDs slightly outside the LISA frequency range, and only a small contamination comes from single WDs and wider binaries. We find that the colour distributions of close DWDs represent a powerful tool to distinguish potential LISA sources from the broader WD population. This is an avenue to select candidates for further follow-up and identification.

The process of galaxy evolution over cosmic time is not yet fully understood, since there is a debate on the impact of galaxy collisions on the star formation and metallicity. The local environment of the galaxy mergers could also have a large impact on the evolution of the galaxies, but it has not yet been possible to examine it in detail. Modern simulations with larger capacity, including the newest physical knowledge and new observations with JWST, help us to answer these questions. Using the IllustrisTNG cosmological simulation, we processed the catalogue data and the merger tree files of the TNG300-1 simulation. We calculated the galaxies average star formation rate (SFR) and mass at redshifts between 0 < z < 15. We investigated the environment of galaxy mergers, with the focus on the local density, and also examined how the SFR changes in merging galaxies. We compared our findings with JWST results and highlighted differences in the star formation rate density (SFRD) history between the models and observations.

We demonstrate a novel beam pattern measurement method for the side lobe characterization of cosmic microwave background telescopes. The method employs a power-variable artificial microwave source under feedback control from the detector under test on the telescope. It enables us to extend the dynamic range of the beam pattern measurement without introducing nonlinearity effects from the detector. We conducted a laboratory-based proof-of-concept experiment, measuring the H-plane beam pattern of a horn antenna coupled to a diode detector at 81 GHz. We achieved a dynamic range of 77.7 dB in the beam pattern measurement, with 60.3 dB attributed to the feedback control. In addition, we verified the measurement by comparing it with other reference measurements obtained using conventional methods. The method is also applicable to general optical measurements requiring a high dynamic range to detect subtle nonidealities in the characteristics of optical devices.

We present an implementation of Quantum Computing for a Markov Chain Monte Carlo method with an application to cosmological functions, to derive posterior distributions from cosmological probes. The algorithm proposes new steps in the parameter space via a quantum circuit whose resulting statevector provides the components of the shift vector. The proposed point is accepted or rejected via the classical Metropolis-Hastings acceptance method. The advantage of this hybrid quantum approach is that the step size and direction change in a way independent of the evolution of the chain, thus ideally avoiding the presence of local minima. The results are consistent with analyses performed with classical methods, both for a test function and real cosmological data. The final goal is to generalize this algorithm to test its application to complex cosmological computations.

The fundamental process of star formation in galaxies involves the interplay between the fueling of star formation via molecular gas and the feedback from recently formed massive stars. This process, by which galaxies evolve, is also closely connected to the intrinsic properties of the interstellar medium (ISM). To study the role that different molecular and atomic phases of the ISM play in star formation, and to characterize their physical conditions, we zoom into our nearest neighboring galaxy, the Large Magellanic Cloud (LMC; 50 kpc). The LMC offers a view of the ISM and star formation conditions in a low metallicity environment similar to, in that regard, the epoch of the peak of star formation in the earlier universe. We present an unprecedentedly detailed analysis of a well-known star-forming regions (SFRs) at a spatial resolution of a few pc. We mapped a 610pcx260pc region in the LMC molecular ridge in [CII] and the [OIII] using the FIFI-LS instrument on the SOFIA telescope. We compare the data with the distribution of the CO (2-1) emission from ALMA, the modeled TIR luminosity as well as Spitzer/MIPS continuum and Halpha. We also provide a detailed description of the observing strategy and the data reduction. We find that [CII] and [OIII] emission is associated with the SFRs in the molecular ridge, but also extends throughout the mapped region, not obviously associated with ongoing star formation. The CO emission is clumpier than the [CII] emission and we find plentiful [CII] present where there is little CO emission, possibly holding important implications for CO-dark gas. We find a clear trend of the [CII]/TIR ratio decreasing with increasing TIR. This suggests a strong link between the [CII]-deficit and the local physical conditions instead of global properties.

Convection in the innermost shells of massive stars plays an important role in initiating core-collapse supernovae. When these convective motions reach the supernova shock, they create extra turbulence, which helps energize the explosion. In our earlier work, we studied the effect of rotation on the hydrodynamic evolution of convective vortices in collapsing stars. This study focuses on how rotation influences the entropy perturbations, which naturally form in turbulent convection. As these perturbations are carried inward with the collapsing star, they generate both vorticity and sound waves. Using linear perturbation theory, we model entropy waves as small disturbances on top of a steady background flow. Our results show that stellar rotation has little effect on the evolution of entropy perturbations during collapse, prior to encountering the supernova shock. This outcome is consistent with our earlier findings on the limited influence of rotation in the accretion of convective eddies.

Dwarf-dwarf galaxy encounters are a key aspect of galaxy evolution as they can ignite or temporarily suppress star formation in dwarfs and can lead to dwarf mergers. However, the frequency and impact of dwarf encounters remain poorly constrained due to limitations of spectroscopic studies, e.g. surface-brightness incompleteness of optical studies and poor spatial resolution of single-dish neutral hydrogen (HI) surveys. We aim to quantify the frequency of isolated gas-rich dwarf galaxy multiples using the untargeted, interferometric Apertif HI survey and study the impact of the interaction on star formation rates of galaxies as a function of the on-sky separation. Our parent dwarf sample consists of 2481 gas-rich galaxies with stellar masses 10^6 < M_* / M_Sun < 5*10^9, for which we identify close companions based on projected separation (r_p) and systemic velocity difference (Del_V_sys). We explore both constant thresholds for r_p and Del_V_sys corresponding to 150 kpc and 150 km/s on all galaxies in our sample, and mass-dependent thresholds based on a stellar-to-halo mass relation. We find the average number of companions per dwarf in our sample to be 13% (20%) when considering mass-dependent (constant) thresholds. In the stellar mass regime of 2*10^8 < M_* / M_Sun < 5*10^9, we find a three times higher frequency (11.6%) of dwarf companions than previously determined from optical spectroscopic studies, highlighting the power of HI for finding dwarf multiples. Furthermore, we find evidence for an increase in star formation rates (SFRs) of close dwarf galaxy pairs of galaxies with similar stellar masses.

We present a new approach to detecting and characterizing a magnetic field in protoplanetary disks through the differential broadening of unpolarized molecular emission from CN. To demonstrate this technique, we apply it to new ALMA observations of the full complement of hyperfine components from the $N=1-0$ transition, achieving a spatial and spectral resolution of ${\approx}\,0.5^{\prime\prime}$ and $80~{\rm m\,s^{-1}}$, respectively. By fitting a model that incorporates the velocity structure of the disk, the potential non-LTE excitation of the molecule, and the Zeeman effect, we recover a radially resolved magnetic field with a strength of ${\sim}10~{\rm mG}$ between 60 and 120~au. The morphology of the field is also inferred through azimuthal variations in the line broadening, revealing a predominantly poloidal field at 60~au, sharply transitioning to one within the disk plane outside of the gap at 82~au. The signal-to-noise ratio of the data meant that the planar component was unable to be decomposed into toroidal and radial components. Lower limits on the local gas density ($n({\rm H_2}) \gtrsim 10^8~{\rm cm^{-3}}$) from the excitation analysis of the CN emission correspond to a lower limit between 0.1 and 0.01 for the plasma $\beta$.

We explore the possibility of axion-photon conversion as a common origin of two low-frequency anomalies: the isotropic radio excess (ARCADE2) and the deep global 21-cm absorption trough (EDGES). From the axion-photon action in an FLRW background with primordial magnetic fields (PMFs), we derive the scale-dependent conversion probability including plasma effects. Resonant conversion, arising when the axion mass matches the plasma-induced photon mass, produces soft photons in the MHz-GHz range. By modeling stochastic PMFs with amplitude $B_0$ and spectral index $n_{\rm B}$, we show that axion-like particles with mass $\sim 10^{-14}$-$10^{-12}\,\mathrm{eV}$ and nanogauss-level nearly scale invariant PMFs can explain both ARCADE2 and EDGES. Heating from PMF dissipation via ambipolar diffusion and turbulent decay reduces the 21-cm trough, shifting the viable parameter space. Our results stem from a consistent theoretical framework developed from first principles and a combined analysis of the radio excess and global 21-cm signal, while remaining consistent with CMB bounds on PMFs and $\Delta N_{\rm eff}$. We conclude that global 21-cm observations may offer potential sensitivity to axions, primordial magnetism, and dark-sector physics.

We present a comprehensive Bayesian spectral analysis of the black hole X-ray binary 4U 1630-47 during its 2022 outburst, using simultaneous \textit{NICER} and \textit{NuSTAR} observations. Using the traditional frequentist approach, we build our model combining reflection spectroscopy with continuum fitting techniques and analyse the data. In the Bayesian framework, we jointly constrain the black hole's spin, mass, inclination, and distance within a unified framework. Employing nested sampling, we capture parameter degeneracies and rigorously propagate both statistical and systematic uncertainties. Our results yield robust and precise spin measurements from both approaches. Our Bayesian analysis fetches spin $a_*= 0.93_{-0.04}^{+0.05}$, mass $M_{\rm BH} = 9.0_{-2.0}^{+2.0} \, M_\odot$, distance $d_{\rm BH} = 10.5_{-1.2}^{+1.3}$~kpc, and inclination angle $i=53.8_{-1.3}^{+1.3}$~deg. It also demonstrates the power of Bayesian inference in fetching valuable insights into the complex physics of black hole accretion and enabling high-confidence measurements of fundamental parameters.

We investigate the relationship between quasars and close companion galaxies using photometry and spectroscopy from the large data set of the SDSS DR16 survey. From the SDSS-QSO catalog of quasars, we selected objects with 0.1 < z < 0.35 and absolute magnitude M(r) <-21.3. For all these targets, we searched for candidate companion galaxies that have a projected distance from the target < 700 kpc and a radial velocity difference from the QSO DeltaV < 1000 km/sec. We find that in 447 QSO at least one companion galaxy is associated. A total of 691 associated galaxies are found in these QSO fields. In the majority of them there is just one associated galaxy but in few cases, many companions are also discovered and in two cases the QSO is found in rich galaxy environments. The possible contamination due to a chance projection of the companion galaxies is found less than 5%. Based on the available data we expect on average to find \sim2 associated companion galaxies for each QSO. A small fraction (13%) of the companion galaxies exhibits [OII] emission lines as signature of recent star formation. However a similar fraction (16%) of unassociated galaxies in the same QSO fields show [OII] emission. This study suggests that there is no significant link between the presence of these close companion galaxies or the signature of recent star formation and the QSO nuclear activity.

Unstable mass transfer in binary systems can lead to transients such as luminous red novae (LRNe). Observations of such transients are valuable for understanding and testing models of mass transfer. For donor stars in the Hertzsprung gap, there can be a long phase of mass-transfer evolution before instability sets in. Only few case studies of such delayed dynamical instability (DDI) mass transfer exist. None consider the full pre-instability evolution and the effects thereof on the observable properties of a binary. We systematically analyse detailed models of stable and unstable mass transfer for Hertzsprung-gap donors. We focus on identifying observable evolutionary features characteristic of ultimately unstable mass transfer and not found in stable mass-transfer binaries. Our binary evolution models, calculated with the MESA code, cover initial donor masses between $2.5 M_{\odot}$ and $10 M_{\odot}$ and initial accretor-to-donor mass ratios between $0.1$ and $1$. We find that the pre-instability evolution is qualitatively the same for all DDI donor stars, consisting of a long slow dimming phase followed by a shorter phase of rapid brightening. The latter phase is powered by recombination of hydrogen and accompanied by a strong increase in effective temperature, unique to unstable mass-transfer binaries. We estimate that a significant fraction of the rapid brighteners should be detectable by Gaia throughout the Galaxy. We model the progenitors of LRNe M31-2015 and V838 Mon, find a higher initial donor mass for M31 2015 than past estimates, and propose a new scenario for V838 Mon in which the known tertiary star dominates the pre-outburst photometry and the outburst results from the DDI of a more massive primary star. This work provides a comprehensive framework linking theory to observations of transients and enables improved classification and prediction of mass-transfer events.

We wish to confirm the nature of YSES 2b, a purportedly faint companion of the young star YSES 2. We used on-sky observations from SPHERE and GRAVITY to measure the astrometric position of 2b with respect to the star YSES 2, and examined the competing hypotheses of (i) a bound substellar companion versus (ii) a distant unrelated background source with a non-zero proper motion. YSES 2b appears to be a late-type M-dwarf star over 2 kiloparsecs behind the star YSES 2. It has a transverse velocity of about 300 km/s and is located within one of the spiral arms of the Galaxy. The main discriminant was multiple epochs of GRAVITY astrometry that identified the sub-milliarcsecond parallactic motion of the star.

We present a simple yet effective method to resolve prior-volume effects, also known as projection effects, in full-shape analyses of the power spectrum multipoles within the Effective Field Theory of Large-Scale Structure (EFTofLSS). By re-defining the EFTofLSS nuisance parameters to incorporate the contribution from the parameters impacting the amplitude of the EFTofLSS modelling components, we substantially mitigate projection effects. With the re-parametrisation the actual posterior maximum values are within the marginalised credible interval, eliminating significant shifts observed in the baseline analysis. We demonstrate the robustness of this method in full-shape $w_0w_a$CDM analyses on synthetic data in BOSS DR12 and DESI DR1 setups. For the evolving dark energy model, we then analyse the BOSS DR12 measurements, in combination with BAO information (from BOSS DR12, 6DF, SDSS DR7 MGS and eBOSS DR16 surveys) and 3$\times$2pt measurements from DES Y3 -- all data combinations are converging into the $w_0-w_a$ parameter region preferred by DESI+CMB+SNIa. From total combination of these large-scale structure probes without additional CMB information we find $w_0=-0.72 \pm 0.21, \, w_a=-0.91^{+0.78}_{-0.64}$. Despite the low significance of deviation from standard cosmology, this result underscores the potential of our re-parametrisation approach in delivering low-redshift cosmological constraints. We argue for the use of this approach in spectroscopic Stage IV surveys, where the potential deviation from standard cosmology can be detected with higher significance.

Astronomical observations at millimeter and submillimeter wavelengths heavily depend on the amount of Precipitable Water Vapor (PWV) in the atmosphere, directly affecting the sky transparency and degrading the quality of the signals received by radio telescopes. Predictions of PWV at different forecasting horizons is crucial to support telescope operations, engineering planning, and observational scheduling and efficiency of radio observatories installed in the Chajnantor area in northern Chile. We developed and validated a Long Short-Term Memory (LSTM) deep learning-based model to predict PWV at forecasting horizons of 12, 24, 36, and 48 hours using historical data from two 183 GHz radiometers and a weather station in the Chajnantor area. We find the LSTM method is able to predict PWV in the 12 and 24 hours forecasting horizons with Mean Absolute Percentage Error (MAPE) of 22% compared to 36% of the traditional Global Forecast System (GFS) method used by Atacama Pathfinder EXperiment (APEX) and the Root Mean Square Error (RMSE) in mm are reduced by 50%. We present a first application of deep learning techniques for preliminary predictions of PWV in the Chajnantor area. The prediction performance shows significant improvements to traditional methods in 12 and 24 hours time windows. We also propose upgrades to improve our method in short (< 1 hour) and long (> 36 hours) forecasting timescales for future work.

We study the build-up and survival of angular momentum in the stellar disc using a statistical suite of cosmological simulations of Milky Way-mass galaxies. Our results show that stellar kinematics at $z=0$ rarely recover the true times of disc spin-up, due to the disruptive impact of massive radial merger events. The proto-disc (i.e. Aurora) and kicked-up disc stars (the Splash) become indistinguishable at low metallicities, and the local fraction of kicked-up disc stars remains $<20$ per cent even in the case of major mergers. In contrast, observations from Gaia and legacy surveys reveal that Galactic $\alpha$-rich populations as old as $\tau=13.5\,\rm{Gyr}$ show significant rotation, with median $\eta > 0.75$. This places strong constraints on the merger ratio between the proto-Milky Way and its last significant merger (Gaia-Sausage Enceladus, GSE), favouring values of $\lesssim 1:7$. We present the age-metallicity relation for the stellar halo and estimate the interaction epoch at $\tau_{\rm{spin\text{-}up}}\simeq\tau_{\rm{GSE}}\sim11\,\rm{Gyr}$. We note an abrupt dearth of halo and Splash stars after a lookback time of $10\,\rm{Gyr}$, marking the end of the merger interaction. Finally, we show that Globular Clusters in the metallicity range $-0.8<\rm{[Fe/H]}<-0.3$ align with a formation time of $\tau_{\rm{starburst}}\sim11\,\rm{Gyr}$, which we interpret as a signature of a starburst triggered by the first pericentric interaction of the GSE. This is remarkable corroboration between our GSE interaction and starburst times of $\tau_{\rm GSE}=\tau_{\rm starburst} \sim 11\,$Gyr.

Globular clusters (GCs) are among the oldest and densest stellar systems in the Universe, yet how they form remains a mystery. Here we present a suite of cosmological simulations in which both dark-matter-free GCs and dark-matter-rich dwarf galaxies naturally emerge in the Standard Cosmology. We show that these objects inhabit distinct locations in the size-luminosity plane and that they have similar ages, age spread, metallicity and metallicity spread to globulars and dwarfs in the nearby Universe. About half of our simulated globulars form by means of regular star formation near the centres of their host dwarf, with the rest forming further out, triggered by mergers. The latter are more tidally isolated and more likely to survive to the present day. Finally, our simulations predict the existence of a new class of object that we call 'globular-cluster-like dwarfs' (GCDs). These form from a single, self-quenching, star-formation event in low-mass dark-matter halos at high redshift and have observational properties intermediate between globulars and dwarfs. We identify several dwarfs in our Galaxy, such as Reticulum II (refs. 2-4), that could be in this new class. If so, they promise unprecedented constraints on dark-matter models and new sites to search for metal-free stars.

Recently, the KM3NeT experiment reported the detection of a neutrino with exceptionally high energy E = 220 PeV, whose origin remains unclear. The corresponding value of the neutrino flux is in tension with the results of other high-energy neutrino experiments. In this study, we discuss the possibility that this neutrino is cosmogenic, i. e., produced by ultra-high energy cosmic rays (UHECR) during their propagation through the intergalactic medium. We adopt the UHECR flux models derived by the Telescope Array experiment, which features a predominantly light mass composition. We show that the predictions of the cosmogenic neutrino flux in these models are consistent with the measurements of the KM3NeT-only and with that of the "global neutrino observatory" at approximately 2${\sigma}$ level. Notably, this result is achieved in a minimal version of the UHECR flux models, that assume one source population with a standard cosmological evolution. We also estimate the corresponding cosmogenic gamma-ray flux and show that it is consistent with Fermi-LAT IGRB measurements and UHE gamma-ray limits; the improvement of the latter can probe these predictions in future.

We report the discovery of a new black-widow millisecond pulsar, PSR J1544-2555, associated with the Fermi-LAT source 4FGL J1544.2-2554. Optical, radio, and gamma-ray observations confirmed its nature as a compact spider binary system. Optical photometry from ULTRACAM revealed a \(\sim\)2.7-hour orbital period, guiding MeerKAT observations that detected \(\sim\)2.4-ms radio pulsations. Subsequent timing campaigns using the Murriyang Parkes Telescope, the Effelsberg 100-m Radio Telescope, and the Nançay Radio Telescope allowed us to obtain a preliminary timing solution, which enabled us to find gamma-ray pulsations. The final timing solution, spanning 16 years of Fermi-LAT gamma-ray data, also displays orbital period variations typical of spider pulsars. X-ray observations from eROSITA indicate non-thermal emission, but the relatively low count rate prohibits the search for X-ray pulsations. Optical light curve modelling using Icarus suggests the asymmetry is best explained by a spot model, where uneven heating creates localised temperature variations on the companion. While the optical spectra we obtained are compatible with the physical properties we infer for the companion star, they were not of sufficient signal-to-noise to allow for radial velocity measurements, thus limiting constraints on the neutron star's mass. The observed bluer colour near the light curve minimum suggests possible non-thermal emission from intra-binary shocks, supported by the presence of an X-ray source. This discovery exemplifies the proven capability of the Fermi-LAT catalogue in identifying millisecond pulsar candidates and highlights the role of optical surveys in detecting variable sources suitable for radio follow-up.

We study the trispectrum in a two-phase USR-SR setup of inflation in which the USR stage is extended in the initial phase of inflation while the second stage of inflation proceeds via a slow-roll phase. A key role is played by the sharpness parameter which controls how quickly the system reaches the final attractor phase after the USR stage. We employ both $\delta N$ and in-in formalisms and calculate trispectrum and the corresponding dimensionless parameters $g_{NL}$ and $\tau_{NL}$. We show that both approaches yield the same results and study the shapes of trispectrum in various configurations. It is shown that the maximum value of trispectrum occurs in the setup with an infinitely sharp transition to the attractor phase while much of trispectrum is washed out in the opposite limit of a mild transition.

Recent experiments have measured the Galactic $\gamma$-ray diffuse emission up to PeV energies, opening a window to study acceleration of Galactic cosmic rays and their propagation up to the cosmic-ray knee. Furthermore, these observations provide a powerful tool to set strong constraints into very-heavy dark matter particles, with masses in the TeV-PeV range. In this paper, we explore the potential of the newest observations of diffuse emissions at the Galactic plane from HAWC and LHAASO to probe this kind of dark matter over a wide mass range. Here, we model secondary emissions (inverse-Compton) from the electrons and positrons produced in the annihilation/decay of dark matter, on top of their prompt $\gamma$-ray emission, including the effects of absorption of high-energy photons via pair production. Furthermore, we show that including the astrophysical backgrounds (namely diffuse emission from cosmic-ray collisions or emission from unresolved sources) can significantly improve these limits. We find that the new measurements provided, specially by LHAASO with the combination of the WCDA and KM2A detectors, allow us to set strong constraints in decaying dark matter, being competitive and even improving the strongest constraints at the moment. We also highlight that these regions lead to constraints that are less affected by uncertainties from the dark matter distribution and discuss how CTA north and SWGO will be able to improve limits in this mass range.

JWST NIRCam and MIRI photometry of NGC 602, a low-metallicity young star cluster in the Small Magellanic Cloud, reveals an extended mid-infrared bright emission feature designated as MZS-1. This feature is prominent between 10 and 25.5 microns, but is extremely faint at 7.7 microns and entirely undetected at shorter wavelengths. MZS-1 exhibits an elliptical morphology with a major axis of approximately 8 arcseconds and a minor axis of about 4 arcseconds. Its elongated shape and multiple emission peaks in the two-dimensional flux map suggest a group of deeply embedded sources with blackbody-like temperatures ranging from 100 K to 140 K. SED fitting using the Robitaille 2017 model grids identifies these sources as Stage I young stellar objects (YSOs) with masses below approximately 3 solar masses and a total stellar mass of the protocluster of about 300 solar masses (based on a Salpeter IMF). The low YSO masses are consistent with their absence in Spitzer-based catalogs due to sensitivity limits. By revealing a deeply embedded, low-mass protocluster invisible in previous surveys, this work highlights JWST's unparalleled resolution and sensitivity in uncovering the earliest stages of low-mass cluster formation in the metal-poor regime.

We present an investigation into the quenching of simulated galaxies across cosmic time, honing in on the role played by both intrinsic and environmental mechanisms at different epochs. In anticipation of VLT-MOONRISE, the first wide-field spectroscopic galaxy survey to target cosmic noon, this work provides clear predictions to compare to the future observations. We investigate the quenching of centrals, high-mass satellites, and low-mass satellites from two cosmological hydrodynamical simulations: IllustrisTNG and EAGLE. Satellites are split according to bespoke mass thresholds, designed to separate environmental and intrinsic quenching mechanisms. To determine the best parameter for predicting quiescence, we apply a Random Forest classification analysis for each galaxy class at each epoch. The Random Forest classification determines supermassive black hole mass as the best predictor of quiescence in centrals and high-mass satellites. Alternatively, the quenching of low-mass satellites is best predicted by group halo mass, at all epochs. Additionally, we investigate the evolution in the dependence of the quenched fraction with various parameters, revealing a more complex picture. There is strong evidence for the rejuvenation of star formation from z = 2 to z = 0 in EAGLE, but not in IllustrisTNG. The starkest discrepancy between simulations rests in the mass threshold analysis. While IllustrisTNG predicts the existence of environmentally quenched satellites visible within the survey limits of MOONRISE, EAGLE does not. Hence, MOONRISE will provide critical data that is needed to evaluate current models, and constrain future models, of quenching processes.

Noise subtraction is a crucial process in gravitational wave (GW) data analysis to improve the sensitivity of interferometric detectors. While linear noise coupling has been extensively studied and successfully mitigated using methods such as Wiener filtering, subtraction of non-linearly coupled and non-stationary noise remains a significant challenge. In this work, we propose a novel independent component analysis (ICA)-based framework designed to address non-linear coupling in noise subtraction. Building upon previous developments, we derive a method to estimate general quadratic noise coupling while maintaining computational transparency compared to machine learning approaches. The proposed method is tested with simulated data and real GW strain data from KAGRA. Our results demonstrate the potential of this framework to effectively mitigate complex noise structures, providing a promising avenue for improving the sensitivity of GW detectors.

Despite their importance, the dynamical properties of massive stars remain poorly understood. Rotation is a key driver of internal mixing and angular momentum transport, significantly influencing massive star evolution. However, constraining rotation from spectroscopy is challenging, as spectral lines often exhibit excess broadening beyond rotation. The origin of this additional broadening, typically attributed to large-scale velocity fields and commonly referred to as macroturbulence, remains uncertain and unconstrained. Here, we present the combined analysis of TESS photometry and rapid time-series spectroscopy using the high-resolution ESPRESSO instrument at the Very Large Telescope of the European Southern Observatory for the rapidly rotating and pulsating massive star Zeta Ophiuchi. Leveraging excellent temporal coverage, our analysis demonstrates that pulsation-induced variability leads to peak-to-peak scatter as large as 88 km/s in the observed macroturbulent velocity time series. We also demonstrate that time-averaged macroturbulent velocities are spectral line specific and can exceed 100 km/s . Furthermore, the macroturbulent velocities from shorter integration times are systematically lower than those derived from stacked spectra mimicking longer exposure times typically needed for fainter stars. These results highlight the role of pulsations in driving variable macroturbulence in massive stars, while also pointing to a potential bias in spectroscopic estimates of macroturbulence for fainter massive stars.

The astrophysical origin of binary black hole mergers discovered by LIGO and Virgo remains uncertain. Efforts to reconstruct the processes that lead to mergers typically rely on either astrophysical models with fixed parameters, or continuous analytical models that can be fit to observations. Given the complexity of astrophysical formation mechanisms, these methods typically cannot fully take into account model uncertainties, nor can they fully capture the underlying processes. Here, we present a merger population analysis that can take a discrete set of simulated model distributions as its input to interpret observations. The analysis can take into account multiple formation scenarios as fractional contributors to the total set of observations, and can naturally account for model uncertainties. We apply this technique to investigate the origin of black hole mergers observed by LIGO Virgo. Specifically, we consider a model of AGN assisted black hole merger distributions, exploring a range of AGN parameters along with several {SEVN} population synthesis models that vary in common envelope efficiency parameter ($\alpha$) and metallicity ($Z$). We estimate the posterior distributions for AGN+SEVN models using $87$ BBH detections from the $O1--O3$ observation runs. The inferred total merger rate is $46.2 {Gpc}^{-3} {yr}^{-1}$, with the AGN sub-population contributing $21.2{Gpc}^{-3}{yr}^{-1}$ and the SEVN sub-population contributing $25.0 {Gpc}^{-3} {yr}^{-1}$.

The Five-hundred-meter Aperture Spherical Telescope (FAST), the world's largest single-dish radio telescope, lists the search for extraterrestrial intelligence (SETI) as one of its key scientific objectives. In this work, we present a targeted SETI observation for 7 nearby active stars utilizing the FAST L-band multibeam receiver, employing a observational strategy that combines position switching with multibeam tracking to balance on-source integration time with the accuracy of the beam response. Using both pulsar and SETI backends, we perform a comprehensive search for narrowband drifting signals with Doppler drift rates within $\pm 4 \ \mathrm{Hz\ s^{-1}}$ and channel-width periodic signal with periods between 0.12 and 100 s and duty cycles between 10\% and 50\%. No credible radio technosignatures were detected from any of the target systems. Based on this null result, we place constraints on the presence of transmitters at a 95\% confidence level, ruling out narrowband transmitters with with EIRP above $3.98\times10^8 \,\mathrm{W}$ and periodic transmitter with EIRP above $1.80\times10^{10} \,\mathrm{W}$,respectively, within the observation band.

One of the most pressing problems in current cosmology is the cause of the Hubble tension. We revisit a two-rung distance ladder, composed only of Cepheid periods and magnitudes, anchor distances in the Milky Way, Large Magellanic Cloud, NGC 4258, and host galaxy redshifts. We adopt the SH0ES data for the most up-to-date and carefully vetted measurements, where the Cepheid hosts were selected to harbour also Type Ia supernovae. We introduce two important improvements: a rigorous selection modelling and a state-of-the-art density and peculiar velocity model using Manticore-Local, based on the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. We infer $H_0 = 71.7 \pm 1.3\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$, assuming the Cepheid host sample was selected by estimated supernova magnitudes. Less plausible selection criteria shift $H_0$ by about one standard deviation. The posterior has a lower central value and a 45 per cent smaller error than a previous study using the same data. The result is also slightly lower than the supernova-based SH0ES inferred value of $H_0 = 73.2 \pm 0.9\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$, and is in $3.3\sigma$ tension with the latest standard cosmological model microwave background results. These results demonstrate that a measurement of $H_0$ of sufficient precision to weigh in on the Hubble tension is achievable using second-rung data alone, underscoring the importance of robust and accurate statistical modelling.

The reionization optical depth $\tau_{\rm reio}$ has interesting connections to existing cosmological anomalies. As first studied in the context of the Hubble tension in our previous paper, a larger $\tau_{\rm reio}$, which could be achieved by removing the Planck low-$\ell$ polarization data, could boost $H_0$ slightly, resulting in a mild reduction of the tension between the early- and late-universe determinations of $H_0$. It has been shown later that a larger $\tau_{\rm reio}$ could also relieve other anomalies including: the tension between BAO and CMB data, the neutrino mass tension, and the latest DESI plus supernovae data's tension with the standard cosmological constant scenario. In this paper, we systematically analyze the correlations between $\tau_{\rm reio}$ and relevant cosmological parameters in the existing cosmic observation anomalies. In addition to Pearson correlation coefficients extracted directly from the covariance matrix, we also study partial correlation coefficients which measure intrinsic relationships between pairs of parameters removing the influence of other parameters. We show that $\tau_{\rm reio}$ has weak intrinsic correlations with the parameters responsible for the tensions and anomalies discussed. The large direct Pearson correlations that allow larger $\tau_{\rm reio}$ inferences to alleviate the cosmological tensions each arise from complicated networks through multiple parameters. As a result, the relationships between $\tau_{\rm reio}$ and each anomaly are not independent of each other. We also employ our method of computing correlations to clarify the impact of large scale polarization data, and comment also on the effects of CMB observations from ACT and SPT.

We study an inflation model driven by a triplet of antisymmetric tensor fields, with minimal and nonminimal couplings to gravity. First, we show that the presence of a triplet of antisymmetric tensor fields can provide inherent background isotropy in the stress-energy tensor contrary to the past studies using an antisymmetric tensor field. Inflation is supported in the presence of non-minimal couplings with gravity. We perform the slow roll analysis and also analyse perturbations to the antisymmetric tensor field as well as the tensor modes of perturbed metric. The speed of gravitational waves manifested from the tensor perturbations is tuned to $c$. We also study the evolution of the gravitational waves, calculate their power spectrum and tensor spectral index.

We investigate the scalar and vector modes arising from cosmological perturbations within the framework of an inflationary scenario driven by an antisymmetric tensor field, minimally coupled to gravity. After eliminating gauge artifacts, there remain four scalar and six vector modes of interest which can be studied separately. We analyze the stability of these modes, while looking for generic instabilities like ghost and gradient instabilities that could potentially plague the theory. Further, we investigate the evolution of these modes across different regimes, particularly subhorizon and superhorizon scales.

Uplink pre-compensation in ground-to-satellite optical links remains a problem, with an appropriate measurement of the wavefront error in the uplink path often not being available. We present a method that uses successive measurements of the readily accessible downlink beam to perform a tomographic reconstruction of the volume of atmosphere common to the downlink beam and the ground-station field-of-view. These measurements are done through the existing downlink wavefront sensor. From here, an estimate of the wavefront error along the uplink path can be obtained. We evaluate this method in simulation over representative atmospheric conditions and find good performance, especially for situations where the satellite point-ahead angle is many times greater than the atmosphere's isoplanatic angle. Compared to pre-compensation directly using the downlink phase, we find that this method estimates the uplink path with a residual mean-square wavefront error that is up to 8.6 times less. The hardware simplicity of this method makes it a promising solution for uplink pre-compensation implementations targeted at optical communications.

Comprehensive Life Cycle Assessment (LCA) as a tool to account for the full range of environmental impacts of resource use in commodities or services is a first step in reducing these impacts. There is an increasing necessity to account for these aspects in the planning, running and end-of-life of scientific experiments and research infrastructure. In the following, the concept for an Open Research Life Cycle Assessment (ORLCA) repository is presented to support this endeavour. It is designed to comply fully with the principles of findability, accessibility, interoperability, and reusability (FAIR).

We place the first constraints on the dark matter fraction contained within dark photon solitons using the absence of their predicted radio-frequency signatures, or radio silence, following mergers around supermassive black holes. In these dense environments, spiky dark matter density profiles can form that enhance the soliton merger rate. We present a novel estimate of this rate by incorporating both the steepened dark matter profile and the soliton velocity dispersion via the Jeans equation. For galaxies with an initial profile $\rho_\mathrm{DM} \propto r^{-1}$, we find the total merger rate across redshifts $0 \leq z \leq 4$ to be $\Gamma_{\text{merg}}^{\text{TOTAL}} \lesssim 10^{-7}f^2_{\text{DM}}\,\text{Mpc}^{-3}\,\text{day}^{-1}$, where $f_\mathrm{DM}$ is the solitonic fraction of dark matter. This enhanced rate leads to more major merger events in which the generated soliton has a mass exceeding a critical threshold, leading to its decay via the parametric resonance phenomenon that produces brief, narrowband, and energetic radio bursts detectable by fast radio burst surveys. Comparing our predictions with the non-observation of such events, we already obtain $f_\mathrm{DM} \lesssim 10^{-1}$ from the first fast radio burst study. This constraint is strengthened to $f_\mathrm{DM} \lesssim 10^{-2}$ from the Parkes HTRU survey, with CHIME projected to tighten this to $f_\mathrm{DM} \lesssim 10^{-3}$. For larger $f_\mathrm{DM}$, we instead constrain the effective coupling strength between the dark and visible sectors to lie outside $10^{-18}\,\mathrm{GeV^{-1}} \lesssim g \lesssim 10^{-8}\,\mathrm{GeV^{-1}}$ for dark photon masses in the range $10^{-6}\,\mathrm{eV} \lesssim m \lesssim 10^{-4}\,\mathrm{eV}$. Our results establish astrophysical transients as powerful probes of dark sectors, opening a window onto the detectability of ultralight vector fields.

Hydrodynamics And Radiation Diffusion} (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion, MPI, and HPX. Node-level parallelism is delegated to Kokkos, providing a single, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression-test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

We construct a higher-order gravity model including all corrections up to mass dimension six. Starting from the Jordan frame, we derive the field equations and specialize to the FLRW background, where the dynamics take the form of a four-dimensional autonomous system. Focusing on the $R+R^{2}+RR_{\mu\nu}R^{\mu\nu}$ case, we obtain linearized equations in the parameter $\gamma_{0}$ and analyze the resulting phase space. The model exhibits the main desirable features of an inflationary regime, with a slow-roll attractor and a stable critical point corresponding to the end of inflation. Analytical expressions for the scalar spectral index $n_{s}$ and the tensor-to-scalar ratio $r$ show that the model is consistent with Planck, BICEP/Keck, and BAO data if $|\gamma_{0}|\lesssim 10^{-3}$. Moreover, negative values of $\gamma_{0}$ restore compatibility with recent ACT, Planck and DESI results, suggesting that higher-order corrections may play a relevant role in refining inflationary cosmology.

We study the cosmological dynamics and predictions in the theory of warped massive gravity. This set-up postulates a five-dimensional ghost-free massive graviton with a brane-localized four-dimensional massive gravity potential, and has the virtue of raising the strong-coupling scale of the 4D theory. We identify two classes of models that lead to decoupled equations for the scale factor on the brane: one characterized by a particular choice of boundary conditions for the Stückelberg fields and one characterized by a special tuning between the coefficients of the 5D and 4D potentials. In the first case, we find interesting solutions including a cosmological bounce without the need of exotic matter. The second case leads to a modified Friedmann equation, and comparison with data shows the potential of the model to alleviate the Hubble tension.

We investigate neutron stars with nonlinear magnetic monopoles in the framework of the Einstein-nonlinear electrodynamics model, specifically within the Bardeen and Hayward models. Solving the modified Tolman-Oppenheimer-Volkoff equations for three different equations of state, we find that upon reaching the critical magnetic charge $q_{c}$, neutron stars enter frozen states characterized by the critical horizon. This extends the concept of frozen states to compact objects composed of ordinary matter (non-field matter), thereby offering a new perspective for related research.

We perform a detailed analysis of a theoretically motivated dark energy quintessence model which interacts with the dark matter sector of the universe. Utilising the current observational datasets from the Cosmic Microwave Background, Baryon Acoustic Oscillations and Type Ia Supernovae, we constrain the parameters that characterise the strength of the time dependent interaction. We also look at the effect of a warm dark matter component in the context of coupled quintessence. Analysis using Deviance Information Criterion indicates strong preference for the quintessence model coupled with warm dark matter. However, Bayesian evidence analysis shows favor in the direction of $\Lambda$CDM model.