Papers for Thursday, Apr 10 2025

Galaxy surveys are crucial for studying large-scale structure (LSS) and cosmology, yet they face limitations--imaging surveys provide extensive sky coverage but suffer from photo-$z$ uncertainties, while spectroscopic surveys yield precise redshifts but are sample-limited. To take advantage of both photo-$z$ and spec-$z$ data while eliminating photo-$z$ errors, we propose a deep learning framework based on a dual UNet architecture that integrates these two datasets at the field level to reconstruct the 3D photo-$z$ density field. We train the network on mock samples representative of stage-IV spectroscopic surveys, utilizing CosmicGrowth simulations with a $z=0.59$ snapshot containing $2048^3$ particles in a $(1200~h^{-1}\rm Mpc)^3$ volume. Several metrics, including correlation coefficient, MAE, MSE, PSNR, and SSIM, validate the model's accuracy. Moreover, the reconstructed power spectrum closely matches the ground truth at small scales ($k \gtrsim 0.06~h/\rm Mpc$) within the $1\sigma$ confidence level, while the UNet model significantly improves the estimation of photo-$z$ power spectrum multipoles. This study demonstrates the potential of deep learning to enhance LSS reconstruction by using both spectroscopic and photometric data.

We present a new, Bayesian analysis of the highest-resolution optical spectrum of the supermassive black hole (SMBH) binary candidate PG 1302-102, obtained with ESPRESSO@VLT (R \simeq 138, 000). Our methodology, based on robust Bayesian model selection, reveals the presence of multiple narrow emission lines at the expected redshift of the source and confirms (for H\{beta}) and detects (for H{\gamma}) the presence of redshifted broad components. Additionally, we have discovered a very broad and, if it is associated with the H\{beta}, very redshifted component at {\lambda} \simeq 5000Å. We evaluate two scenarios for explaining the observed broad emission line (BEL) in PG 1302-102. In the case in which the redshifted BEL asymmetry arises from the orbital motion of a putative binary, our measurements coupled with simple estimates of the broad-line region (BLR) sizes suggest that the individual black hole BLRs are either settled in a single BLR or in the process of merging and, therefore truncated and highly disturbed. Alternatively, in the scenario of a single SMBH, we explain the distorted emission of the BELs with a nonsymmetric distribution of the BLR clouds; namely, a thin disk with a spiral perturbation. This BLR configuration is statistically preferred over any empirical multi-Gaussian fit and simultaneously explains the asymmetric emission of the H{\beta} and H{\gamma} close to the bulk of the line and any additional excess (or the lack of it, in the case of the H{\gamma}) at much longer wavelengths. The physical origins of the perturbation are unclear, and a connection with the possible presence of a black hole binary cannot be ruled out. Given the growing evidence from theoretical and observational works demonstrating the common presence of disturbed BLRs in active galactic nuclei, we argue that an origin related to self-gravitating instabilities may be more plausible.

The origin of planetesimals ($\sim$100 km planet building blocks) has confounded astronomers for decades, as numerous growth barriers appear to impede their formation. In a recent paper we proposed a novel interaction where the streaming instability (SI) and dust coagulation work in tandem, with each one changing the environment in a way that benefits the other. This mechanism proved effective at forming planetesimals in the fragmentation-limited inner disk, but much less effective in the drift-limited outer disk, concluding that dust traps may be key to forming planets at wide orbital separations. Here we explore a different hypothesis: That vortices host a feedback loop in which a vortex traps dust, boosting dust coagulation, which in turn boosts vortex trapping. We combine an analytic model of vortex trapping with an analytic model of fragmentation limited grain growth that accounts for how dust concentration dampens gas turbulence. We find a powerful synergy between vortex trapping and dust growth. For $\alpha \le 10^{-3}$ and solar-like metallicity this feedback loop consistently takes the grain size and dust density into the planetesimal formation region of the streaming instability (SI). Only in the regime of strong turbulence ($\alpha \ge 3\times 10^{-3}$) does the system often converge to a steady state below the SI criterion. The combination of vortex trapping with dust coagulation is an even more powerful mechanism than the one involving the SI. It is effective at lower metallicity and across the whole disk -- anywhere that vortices form.

The most metal-poor stars found in the Galaxy and in nearby galaxies are witnesses of the early evolution of the Universe. In a general picture in which we expect the metallicity to increase monotonically with time, as a result of the metal production in stars, we also expect the most metal-poor stars to be the most primitive objects accessible to our observations. The abundance ratios in these stars provide us important information on the first generations of stars that synthesised the nuclei that we observe in these stars. Because they are so primitive the modelling of their chemical inventory can be often satisfactorily achieved by assuming that all the metals were produced in a single supernova, or just a few. This is simpler than modelling the full chemical evolution, using different sources, that is necessary at higher metallicity. The price to pay for this relative ease of interpretation is that these stars are extremely rare and require specifically tailored observational strategies in order to assemble statistically significant samples of stars. In this review we try to summarise the main observational results that have been obtained in the last ten years.

Active galactic nuclei (AGN) can accelerate protons to energies of $\sim$ 10-100 TeV, with secondary production of high-energy neutrinos. If the acceleration is driven by magnetized turbulence, the main properties of the resulting proton and neutrino spectra can be deduced based on insights from particle-in-cell simulations of magnetized turbulence. We have previously shown that these properties are consistent with the TeV~neutrino signal observed from the nearby active galaxy NGC 1068. In this work, we extend this result to a population study. We show that the produced neutrino flux depends mainly on the energetics of the corona -- the relative fraction of X-ray, magnetic, and non-thermal proton energy -- and on the spectral energy distribution of the AGN. We find that coronae with similar properties can explain neutrinos from the candidate AGN for which IceCube has reported an excess, albeit less significant than NGC 1068. Building on this framework, we show how the neutrino signal evolves with the AGN luminosity, and use this AGN sequence to predict the diffuse neutrino flux from the extragalactic population, showing that it can account for the diffuse neutrino signal observed by IceCube in the $\sim$1-100 TeV energy range.

Protoplanetary discs with cavities, also known as ``transition discs'', constitute ~10% of protoplanetary discs at sub-mm wavelengths. Among several explanations, one hypothesis suggests these cavities are carved by undetected stellar or planetary companions. We present a novel approach to quantify the likelihood that a cavity-carving companion goes undetected because it is either too close to the star (i.e., has a small projected separation) or too faint to be resolved. We generate two independent samples of stellar and planetary companions with random sky orientations, assuming distributions in eccentricity, mass ratio, and time-weighted orbital phases, to study the statistical properties of the cavities they produce. We calculate the likelihood that a companion appears with a projected separation $d$ relative to its semi-major axis $a_{bin}$ ($d/a_{bin}$). Then, using a disc truncation model, we calculate the likelihood that companions carve a cavity with size $a_{cav}$ relative to its semi-major axis $a_{bin}$ and projected separation $d$, deriving distributions of $a_{bin}/a_{cav}$ and $d/a_{cav}$. We find that stellar companions carve cavities with median sizes ~3 times larger than their projected separation $d$ ($a_{cav}\sim3d$, $a_{cav}\sim1.7 d$ for planets), but with a statistically significant tail towards larger values ($a_{cav}\gg 3d$). We estimate the likelihood that cavity-carving companions go undetected due to projection effects when the system is observed with spatial resolution $R$, $P(d< R)$. Considering observational constraints on companion masses, we apply this framework to 13 well-known transition discs. We find that undetected stellar companions are unlikely in 8 out of 13 systems we considered, with 5 notable exceptions: ABAur, MWC758, HD135344B, CQTau and HD169142. The presence of undetected planets cannot be excluded in any of the transition discs considered.

The chemical evolution of galaxies is shaped by their star formation histories and the exchange of gas with their environments. Metallicity provides key insights into these processes, reflecting the interplay between star formation and gas flows. A fundamental aspect of this evolution is the mass-metallicity relation, which captures the strong correlation between a galaxy stellar mass ($M_\star$) and its gas-phase oxygen abundance. In this study, we use MUSE observations to analyze star-forming disc galaxies in 12 clusters within the redshift range $0.3 < z < 0.5$. Galaxies were classified into three groups: ram-pressure stripping (RPS), control cluster, and control field. For the first time, we investigate the impact of RPS on gas-phase metallicities across a wide mass range of galaxies at intermediate redshift, comparing RPS galaxies to counterparts in both cluster and field environments. By analyzing the integrated flux within galactic disks, our results reveal that, on average, RPS induces a metallicity enhancement of 0.2 dex over non-stripped galaxies. Contrary to the prevailing view that cluster membership alone drives metallicity enrichment, we find that control cluster galaxies exhibit metallicities comparable to field galaxies at a given $M_\star$, with only RPS galaxies displaying significantly higher metal content, highlighting the unique role of RPS in shaping the chemical properties of galaxies. These differences become more pronounced at lower $M_\star$, indicating that environmental influences play a more critical role in shaping the chemical evolution of lower-mass galaxies. Our findings suggest that both enhanced star formation rates and suppressed gas inflows -- consequences of ram pressure stripping -- drive the elevated metallicity observed in RPS galaxies.

Gas within the influence sphere of accreting massive black holes is responsible for the emission of the broad lines observed in optical-UV spectra of unobscured active galactic nuclei. Since the region contributing the most to the broad emission lines (i.e. the broad line region) depends on the active galactic nucleus luminosity, the study of broad line reverberation to a varying continuum can map the morphology and kinematics of gas at sub-pc scales. In this study, we modify a preexisting model for disc-like broad line regions, including non-axisymmetric structures, by adopting an emissivity profile that mimics the observed luminosity-radius relation. This makes our implementation particularly well suited for the analysis of multi-epoch spectroscopic campaigns. After validating the model, we use it to check if strongly non-axisymmetric, single broad line regions could mimic the short time-scale evolution expected from massive black hole binaries. We explore different orientations and anisotropy degrees of the broad line region, as well as different light curve patterns of the continuum to which the broad line region responds. Our analysis confirms that recently proposed algorithms designed to search for massive black hole binaries in large multi-epoch spectroscopic data are not contaminated by false positives ascribed to anisotropic broad line regions around single MBHs.

A plethora of evidence suggests that $\omega$ Centauri ($\omega$ Cen) is the nuclear star cluster of a galaxy that merged with the Milky Way in early times. We use APOGEE, Gaia, MUSE, and HST data supplemented by galaxy chemical evolution models to place constraints on the assembly and chemical enrichment history of $\omega$ Cen. The APOGEE data reveal three stellar populations occupying separate loci on canonical chemical planes. One population resembles metal-poor halo field stars (P1), a second shows light-element abundance anti-correlations typical of metal-poor globular clusters (IM), and a third population (P2) is characterised by an extreme "second-generation" abundance pattern. Both P1 and P2 populations cover a broad range of metallicity, consistent with extended histories of bursty star formation (SF), which is also evident from their light- and $\alpha$-element abundance patterns. Conversely, the IM stars exhibit a narrow metallicity spread, combined with the Al-Mg, Na-O, and C-N anti-correlations common to metal-poor Galactic globular clusters. Moreover, these three populations alone seem to account for the distribution of $\omega$ Cen stars in the chromosome map. We discuss these findings in context of a scenario according to which $\omega$ Cen formed by a combination of in situ SF within the host galaxy (P1), followed by the spiralling in of gas-rich globular clusters (IM), leading to another burst of SF (P2). We perform a robust comparison of the chemical composition of $\omega$ Cen with those of halo substructures well represented in APOGEE DR17, finding no chemical associations to a high confidence level.

The properties of active region nests -- locations on the Sun with recurring flux emergence -- are poorly constrained by observations from Earth alone. ESA's Solar Orbiter now monitors the Sun's far-side for extended periods of time; facilitating observations of the entire Sun. We combined observations from near-Earth satellites and Solar Orbiter to evaluate the contribution of a long-lived active region nest to the Sun's global flaring activity. We identified a location in Carrington coordinates with episodic bursts of flux emergence throughout 2022. The combined observations allowed near-continuous monitoring of this region from April to October; during its most active period. GOES and SolO/STIX were used to compare its flaring activity to that of the entire Sun. The region's morphology was extracted from SDO/AIA and SolO/EUI extreme ultraviolet images and combined with magnetic field measurements from SDO/HMI and SolO/PHI to assess its unsigned magnetic flux. The active region nest grew in complexity from January to May due to repeated flux emergence events; with a peak unsigned magnetic flux of $5\times10^{22}$Mx. The region was responsible for 40% of the observed solar flares in 2022, including five months where it produced 50-70% of all flares over the entire Sun (during the near-continuous monitoring window). Of the 17 complex flaring NOAA active regions in 2022 this region contained 10 but occupied less than 20% of the area in the active latitudes. Active region nests can maintain a high rate of flaring activity for several solar rotations and are more likely to produce complex active regions that can trigger X-class solar flares. Improving the identification and monitoring of long-lived active region nests would benefit short to medium term space weather forecasts.

We present high resolution follow-up data of the planetary microlensing event MOA-2010-BLG-328, using Keck and the Hubble. Keck data, taken 8 years after the event, reveal a strong lens detection enabling a direct measurement of lens flux and source-lens relative proper motion. We find the relative source-lens proper motion to be $\mu_{\rm rel, Hel} = 4.07 \pm 0.34\ \rm mas\ yr^{-1}$, with the lens being $\sim10$ times fainter than the source. The lens was very faint in the Hubble passbands, and the small lens-source separation of $\sim$35 mas made its detection difficult. However, we obtained estimates of the lens magnitudes in Hubble bands by constraining its location to match the Keck K-band detection. The original analysis by \citet{Furusawa2013} reports a degenerate light curve, with several viable models depending on higher-order effects. We attempt to break the degeneracy by remodeling the event using constraints from follow-up data. Our best fit model includes parallax, orbital motion, xallarap and the magnification of a source companion. Models omitting any of these are excluded. However, even with a lens detection the solution remains unclear, as the degeneracy between a nearby late M dwarf and a distant early M dwarf in the disk persists, and cannot be broken with NIR data alone. We conclude the lens is either a $\sim0.2\ M_{\odot}$ star at $2-3$kpc, or a $\sim0.5\ M_{\odot}$ star at $4-5$kpc. This study highlights the importance of multi-band data and comprehensive modeling to resolve microlensing degeneracies.

Emission from two massive black holes (MBHs) bound in a close binary is expected to be modulated by different processes, such as Doppler boost due to the orbital motion, accretion rate variability generated by the interaction with a circumbinary disc, and binary gravitational self-lensing. When the binary is compact enough, the two black holes are thought to be surrounded by a common broad line region that reprocesses the impinging periodically varying ionising flux, creating broad emission lines with variable lineshapes. Therefore, the study of broad emission line variability through multi-epoch spectroscopic campaigns is of paramount importance for the unambiguous identification of a binary. In this work, we study the response of a disc-like broad line region to the Doppler-boosted ionising flux emitted by $\mathrm{sub-milli-pc}$ MBH binaries on a circular orbit and compare it with the response of a broad line region illuminated by a single MBH with a periodically but isotropically varying intrinsic luminosity. We show that in the binary case, the time lags of the blue and red wings of the broad emission lines, arising from diametrically opposite sides of the circumbinary disk, are out of phase by half of the binary's orbital period, as they each respond to the periodic "light-house" modulation from the binary's continuum emission. This asymmetric time-lag represents a new binary signature that can not be mimicked by a single MBH.

The presence of a 'butterfly' or a quadrupole structure in the stellar mean radial velocity field of the Milky Way is well known from the Gaia and the APOGEE surveys. Past studies indicated that a stellar bar can excite such a quadrupole feature in the $<V_R>$ distribution. However, a systematic study investigating the co-evolution of bar and quadrupole structure is largely missing. Furthermore, whether this quadrupole structure in $<V_R>$ can be used as a robust kinematic diagnostic to constrain bar properties, particularly for the Milky Way, is still beyond our grasp. Here, we investigate the bar-induced quadrupole feature using a suite of isolated N-body models forming prominent bars and a sample of Milky Way-like barred galaxies from the TNG50 cosmological simulation. We demonstrate that the properties of the quadrupole (strength, length, and orientation) are strongly correlated with the bar properties, regardless of the choice of the stellar tracer population; thereby making the quadrupole feature an excellent kinematic diagnostic for constraining the bar properties. In presence of spirals, the estimator which takes into account the phase-angle of m=4 Fourier moment, serves as a more appropriate estimator for measuring the length of the quadrupole. Further, we constructed a novel Gaia-like mock dataset from a simulated bar model while incorporating the dust extinction and the broad trends of observational errors of the Gaia survey. The quadrupole properties (strength and length) estimated from those Gaia-like mock data are larger (~35-45 percent) when compared with their true values. We determined that the majority of this effect is due to the uncertainty in parallax measurement. This demonstrates the potential caveat of inferring Milky Way's bar properties by using the stellar kinematic information from the Gaia DR3 without properly accounting for the observational uncertainties.

Galactic open clusters (OCs) are subject to internal and external destructive effects that gradually deplete their stellar content, leaving imprints on their structure. To investigate their dynamical state from an observational perspective, we employed Gaia DR3 data to perform a comprehensive analysis of 174 OCs (~10% of Dias et al.'s 2021 catalogue). We employed radial density profiles and astrometrically decontaminated colour-magnitude diagrams to derive structural parameters, distance, mass and time-related quantities. We explored the parameters space and searched for connections relating the clusters' structure with the internal evolutionary state and the external Galactic tidal field. Correlations were verified after segregating the sample according to the Galactocentric distance and half-light to Jacobi radius ratio ($r_h/R_J$). This tidal filling ratio decreases with both the cluster mass and dynamical age. At a given evolutionary stage, OCs with larger $r_h/R_J$ tend to present larger fractions of mass loss due to dynamical effects. Regarding the impact of the external conditions, we identified different evaporation regimes: for ambient densities ($\rho_{\rm{amb}}$) larger than ~0.1M$_{\odot}$/pc$^3$, clusters tend to be more tidally filled as they are subject to weaker tidal stresses. For $\rho_{\rm{amb}}$ $\lesssim$ 0.1M$_{\odot}$/pc$^3$, the opposite occurs: $R_J$ increases for smaller $\rho_{\rm{amb}}$, causing $r_h/R_J$ to decrease. In turn, two-body relaxation tends to compact the cluster core, which is less sensitive to variations of the external potential. The higher the degree of central concentration, the larger the number of relaxation times a cluster takes until its dissolution.

This study sets new constraints on Cold+Warm Dark Matter (CWDM) models by leveraging the small-scale suppression of structure formation imprinted in the Lyman-$\alpha$ forest. Using the Sherwood-Relics suite, we extract high-fidelity flux power spectra from simulated Lyman-$\alpha$ forest data, spanning a broad range of cosmologies and thermal histories. This enables precise constraints on the warm dark matter (WDM) fraction, $f_{\mathrm{WDM}}$, and the mass of the WDM particle, $m_{\mathrm{WDM}}$. A key advancement of our analysis is the integration of a neural network emulator directly at the likelihood level, significantly accelerating Bayesian parameter inference. With new observations of high-redshift ($z$ = 4.2$-$5.0) quasar spectra from UVES and HIRES, we establish stringent upper limits: for $m_{\mathrm{WDM}}$ = 1 keV, we find $f_{\mathrm{WDM}} < 0.16$ (2$\sigma$), with constraints loosening to 35\%, 50\%, and 67\% for $m_{\mathrm{WDM}}$ = 2, 3, and 4 keV, respectively. Our results for pure WDM reaffirm the lower bounds of previous work. Crucially, we account for the fixed resolution of simulations and the impact of patchy reionization, demonstrating their minimal influence on mixed dark matter constraints. This robustness paves the way for tighter bounds with improved statistical samples in the future. Our findings suggest that CWDM models can naturally accommodate mild suppression of matter clustering in the high redshift Lyman-$\alpha$ forest 1D flux power, potentially offering a resolution to some of the ongoing cosmological tensions at low redshifts, namely the $S_{8}$ tension.

Recent high-cadence transient surveys have uncovered a rare subclass of Type Ibc supernovae (SNe) that exhibit an early, blue peak lasting a few days before the main, radioactively powered peak. Since progenitors of Type Ibc SNe are typically compact and lack an extended envelope, this early peak is commonly attributed to shock cooling emission from circumstellar matter (CSM) surrounding the progenitor star. As such, these SNe provide a unique opportunity to constrain the pre-explosion activity of Type Ibc SN progenitors. We present the first systematic study of this Type Ibc SN population that uses hydrodynamic modelling. We simulated Type Ibc SNe exploding in a CSM using the multi-group radiation-hydrodynamics code STELLA, exploring a range of SN and CSM properties. By comparing the resulting theoretical multi-band light curves to a sample of seven Type Ibc SNe with early peaks, we constrained their CSM properties. Assuming a wind-like density distribution for the CSM, we found CSM masses of $10^{-2} - 10^{-1} \ \mathrm{M}_\odot$ and CSM radii of $(1 - 5) \times 10^3 \ \mathrm{R}_\odot$. While the masses were roughly consistent with a previous estimate obtained using an analytical model, the radii were significantly different, likely due to our assumption of spatially spread out CSM. We infer that the progenitors could have created CSM via late-time binary mass transfer or pulsational pair instability. We also estimate that, in the planned ULTRASAT high-cadence survey, $\sim 30$ shock cooling peaks from Type Ibc SNe will be observed.

The dynamical friction force acting on a spatially extended probe (globular clusters and dwarf galaxies) moving in the environment of ultralight bosonic dark matter in the state of the Bose-Einstein condensate is determined. Modeling the probe as a Plumer sphere of radius $l_p$, the radial and tangential components of the dynamic friction force are found in an analytic form, which reduce in the limit $l_p \to 0$ to the corresponding analytic expressions obtained in the literature in the case of a point probe. The dependence of dynamical friction force on boson particle mass $m$ was analyzed and found to be non-monotonous in the interval $10^{-23} - 10^{-21}$eV.

Recovering high-fidelity images of the night sky from blurred observations is a fundamental problem in astronomy, where traditional methods typically fall short. In ground-based astronomy, combining multiple exposures to enhance signal-to-noise ratios is further complicated by variations in the point-spread function caused by atmospheric turbulence. In this work, we present a self-supervised multi-frame method, based on deep image priors, for denoising, deblurring, and coadding ground-based exposures. Central to our approach is a carefully designed convolutional neural network that integrates information across multiple observations and enforces physically motivated constraints. We demonstrate the method's potential by processing Hyper Suprime-Cam exposures, yielding promising preliminary results with sharper restored images.

Spiral galaxies are thin and susceptible to being disrupted vertically. The largest star clusters, and nuclear starbursts, generate enough energy from winds and supernovae to send disk material to the halo. % METHODS Observations of edge-on galaxies allow for the clearest view of vertical disruptions. We present new observations of the nearby, edge-on galaxy NGC 5775 with SOFIA [C II] 157.7 micron and archival images from Hubble in Halpha to search for extraplanar gas. The extraplanar [C II] extends 2 kpc from the midplane over much of the star-forming disk. The extraplanar [C II] at 2 kpc from the midplane approximately follows the rotation of the disk, with a lag of approximately 40 km/s; this lag is similar to what has been previously reported in Halpha. Significant vertical extensions (to 3 kpc) are seen on the northeast side of the galaxy, potentially due to super star clusters in the NGC 5775 disk combined with gravitational interaction with the companion galaxy NGC 5774. The Halpha narrow-band image reveals a narrow plume that extends 7 kpc from the nucleus and is almost exactly perpendicular to the disk. The plume shape is similar to that seen from the comparable galaxy NGC 3628 and may arise from the nuclear starburst. Alternatively, the Halpha plume could be a relic of past activity.

Super-sensitive observations of bright pulsars by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) have revealed weak radio emission continuously emerged in the rotation phases between the main pulse and interpulse of an rotating neutron star. We develop a model for the polarized radio emission radiated from different heights in the pulsar magnetosphere and examine emission intensity distribution over the whole rotation phases of pulsars seen from all directions by the line of sight. We find that for pulsars with small periods and the magnetosphere filled with much more relativistic particles, the polarized radio emission can be generated in all rotation phases for both the aligned and perpendicular rotating neutron stars. When the line of sight cuts the pulsar emission beam between the rotation and magnetic axes, the polarization angles have the same sense of variation gradient for the ``main'' pulse and ``interpulse''. If the line of sight cuts the beams between the inclined magnetic axis and the equator, the opposite senses can be found for the main pulse and interpulse. In addition to the pulsed emission, we find persistent radio emission generated in the pulsar magnetosphere. The model can naturally explain the emission across the entire rotation phases.

Mergers are fundamental to our understanding of the processes driving the evolution of the structure and morphology of galaxies, star formation, AGN activity, and the redistribution of stellar mass in the Universe. Determining the fraction and properties of mergers across cosmic time is critical to understanding the formation of the Universe we observe today. This fraction and its evolution also provide inputs and constraints for cosmological simulations, crucial for theoretical models of galaxy evolution. We present robust estimates of major close-pair fractions and merger rates at $0.2 < z < 0.9$ in the Deep Extragalactic VIsible Legacy Survey (DEVILS). We identify major mergers by selecting close-pairs with a projected spatial separation $r_{\mathrm{sep}} < 20$ h$^{-1}$ kpc and a radial velocity separation $v_{\mathrm{sep}} < 500$ km s$^{-1}$. For galaxies with stellar masses of log$_{10}$($M_\star$/$M_\odot$) = 10.66 $\pm$ 0.25 dex, we find a major close-pair fraction of $\approx 0.021$ at $0.2 < z < 0.34$ using a highly complete, unbiased spectroscopic sample. We extend these estimates to $0.2 < z < 0.9$ by combining the full probability distribution of redshifts for galaxies with high-quality spectroscopic, photometric, or grism measurements. Fitting a power-law $\gamma_{m} = A(1 + z)^m$, we find $A = 0.024 \pm 0.001$ and $m = 0.55 \pm 0.22$. Consistent with previous results, the shallow slope suggests weak redshift evolution in the merger fraction. When comparing with large hydrodynamical simulations, we also find consistent results. We convert close-pair fractions to merger rates using several literature prescriptions for merger timescales and provide all measurements for future studies.

As an interesting subclass of gamma-ray burst (GRB), Type IL GRB (such as GRB 211211A and GRB 230307A) features a long-duration prompt emission but originating from compact binary merger. The "long duration" emisison of Type IL GRB are dominately composed of the main burst, rather than the extended emission, differentiating them from the traditional "long-short" GRB (e.g., GRB 060614). Previous study has reported several Type IL GRBs by visual inspection of their light curves. In this work, we established a detailed criterion to identify Type IL GRBs by light curve, and then systematically searched the archival \textit{Fermi}/GBM data with this criterion, resulting in a sample of 5 type IL GRBs from January 1, 2014 to January 1, 2024, i.e. GRB 230307A, GRB 211211A, GRB 200914A, GRB 200311A and GRB 170228A. Apart from the light curve pattern, we find that the temporal and spectral properties of these 5 GRBs also support this classification. Interestingly, we find that the energy ratio between extended emission and main emission is almost constant ($\sim0.7$, with small scattering) for these GRBs, which have strong implication on the mechanism of Type IL burst. We discuss theoretical models to interpret the progenitor, central engine, and extended emission of these Type IL bursts.

This study presents calculations of rate coefficients, resonance strengths, and cross sections for the dielectronic recombination (DR) of Y^+, Sr^+, Te^2+, and Ce^2+--low-charge ions relevant to kilonovae and non-local thermodynamic equilibrium (non-LTE) plasmas. Using relativistic atomic structure methods, we computed DR rate coefficients under conditions typical of these environments. Our results highlight the critical role of low-lying DR resonances in shaping rate coefficients at kilonova temperatures (~ 10^4 K) and regulating charge-state distributions. Pronounced near-threshold DR resonances significantly influence the evolving ionization states and opacity of neutron star merger ejecta. Comparisons with previous studies emphasize the necessity of including high-n Rydberg states for accurate DR rate coefficients, especially for complex heavy ions with dense energy levels. Discrepancies with existing datasets underscore the need for refined computational techniques to minimize uncertainties. These results provide essential input for interpreting spectroscopic observations of neutron star mergers, including James Webb Space Telescope data. We also put forward suitable candidates for experimental studies, recognizing the challenges involved in such measurements. The data presented here have potential to refine models of heavy-element nucleosynthesis, enhance plasma simulation accuracy, and improve non-LTE plasma modeling in astrophysical and laboratory settings.

Fast radio bursts (FRBs) are bright millisecond radio events of unknown extra-galactic origin. Magnetars are one of the main contenders. Some sources, the repeaters, produce multiple events but so far generally without the characteristic periodicity that one could associate with the spin of a neutron star. We fit a geometrical model to the two main repeaters of the CHIME/FRB catalogue, namely FRB 20180814A and FRB 20180916B. Assuming the bursts originate from a magnetar's magnetosphere, we constrain the spin and magnetic parameters of the star which are encoded into burst spectro-temporal morphologies. We estimate that a very strong toroidal magnetic component together with spin periods of respectively $2.3_{-0.5}^{+0.5} ~ \rm s$ and $0.8_{-0.2}^{+0.1} ~ \rm s$ best explain the data. We argue that this points towards young magnetars with super-twisted magnetospheres.

The RS CVn type binary star GT Mus was observed during its quiescence using the Resolve X-ray microcalorimeter spectrometer onboard XRISM. The main and satellite lines of the Fe XXIV--XXVI K-shell transitions were resolved for the first time from stellar sources. We conducted line ratio analysis to investigate any deviations from collisional onization equilibrium (CIE) and Maxwell electron energy distribution with a single-temperature. By using five combinations of direct excitation lines and dielectronic recombination satellite lines in three line complexes (Fe He$\alpha$, Ly$\alpha$, and He$\beta$), we found that the plasma is well characterized by two-temperature thermal plasmas with temperatures of 1.7 and 4.3 keV, which is consistent with a thermal broadening of Fe XXV and the broadband fitting results in the 1.7--10 keV band. Other forms of deviation from a single-temperature plasma, such as different ionization and electron temperatures or the $\kappa$ distribution for the electron energy distributions, are not favored, which is reasonable for stellar coronae at quiescence. This study demonstrates the utility of the Fe K-shell line ratio diagnostics to probe plasma conditions using X-ray microcalorimeters.

Multiscale studies of the morphology and strength of the magnetic field are crucial to properly unveil its role and relative importance in high-mass star and cluster formation. G31.41+0.31 (G31) is a hub-filament system that hosts a high-mass protocluster embedded in a hot molecular core (HMC). G31 is one of the few sources showing a clear hourglass morphology of the magnetic field on scales between 1000 au and a few 100 au in previous interferometric observations. This strongly suggests a field-regulated collapse. To complete the study of the magnetic field properties in this high-mass star-forming region, we carried out observations with the James Clerk Maxwell Telescope $850 \mu$m of the polarized dust emission. These observations had a spatial resolution of $\sim$0.2 pc at 3.75 kpc. The aim was to study the magnetic field in the whole cloud and to compare the magnetic field orientation toward the HMC from $\sim$50,000 au to $\sim$260 au scales. The large-scale ($\sim$5 pc) orientation of the magnetic field toward the position of the HMC is consistent with that observed at the core ($\sim$4,000 au) and circumstellar ($\sim$260 au) scales. The self-similarity of the magnetic field orientation at these different scales might arise from the brightest sources in the protocluster, whose collapse is dragging the magnetic field. These sources dominate the gravitational potential and the collapse in the HMC. The cloud-scale magnetic field strength of the G31 hub-filament system, which we estimated using the Davis-Chandrasekhar-Fermi method, is in the range 0.04--0.09 mG. The magnetic field orientation in the star-forming region shows a bimodal distribution, and it changes from an NW--SE direction in the north to an E--W direction in the south [abridged abstract].

We analyse the properties of the Comptonizing medium in the black-hole X-ray binary Swift J1727.8$-$1613 using the time-dependent Comptonization model vkompth, using NICER observations of type-C QPOs in the hard and hard-intermediate states. During the 2023 outburst of the source, we measure the rms and phase lags of the QPO across 45 observations as the QPO frequency, $\nu_{\rm QPO}$, evolves from $\sim 0.3$ Hz to $\sim 7$ Hz. By simultaneously fitting the time-averaged spectrum of the source and the rms and lag spectra of the QPO, we derive the evolution of the disk and corona parameters. At $\nu_{\rm QPO} = 0.34$ Hz, the QPO phase lags are hard, with 10 keV photons lagging 0.5 keV photons by $\sim 0.5$ rad. As $\nu_{\rm QPO}$ increases, the lags for the same energy bands decrease, reaching near zero at $\nu_{\rm QPO} \sim 1.2$ Hz, and then reverse to soft lags of $\sim -1.1$ rad at $\nu_{\rm QPO} \sim 7$ Hz. Initially, the inner radius of the accretion disk is truncated at $\sim 30-40 R_g$ (assuming a 10 solar-mass black hole) and, as the QPO frequency increases, the truncation radius decreases down to $\sim 10 R_g$. Initially, two coronas of sizes of $\sim 6.5 \times 10^3$ km and $\sim 2 \times 10^3$ km, extend over the disk and are illuminated by different regions of the disk. As the QPO frequency increases, both the coronas shrink to $\sim 2 \times 10^3$ km at $\nu_{\rm QPO} = 2.5$ Hz. Following a data gap, one corona expands again, peaking at a size of $\sim 2 \times 10^4$ km. We interpret the evolution of the coronal size in the context of accompanying radio observations, discussing its implications for the interplay between the corona and the jet.

We propose a new method to tomographically probe cosmic birefringence using radio galaxies. We show that the redshift evolution of the cosmic birefringence angle induced by a slow-rolling pseudoscalar field, which is a candidate for dynamical dark energy, is independent of the detailed model of the pseudoscalor field. This universal profile evolves predominantly at $z\lesssim10$. In contrast, if the origin is a dark matter-like pseudscalor field, the resulting birefringence angle tends to be negligible in the low-redshift regime. This new insight provides a strong motivation to independently test the cosmic birefringence using polarized astrophysical sources such as radio galaxies. We find that a sample size of $\order{10^5-10^6}$ is required to distinguish the profiles, which is achievable with ongoing and upcoming radio surveys such as ASKAP or SKA.

RR Lyrae stars have long been considered reliable tracers of old, metal-poor populations, primarily due to their prevalence in globular clusters and the Galactic halo. However, the discovery of a metal-rich subpopulation in the Galactic disc, kinematically colder and more rotationally supported, challenges this classical view. Understanding the age of these metal-rich RR Lyrae stars is crucial for constraining their formation pathways and assessing what Galactic populations they are tracing. In this work, we leverage the unprecedented astrometric precision of Gaia DR3 to infer the age distribution of metal-rich RR Lyrae stars through a kinematic comparison with O-rich Mira variables. Mira variables, with their well-established period-age relation, serve as a natural clock, allowing us to transfer age information to RR Lyrae stars via their phase-space properties. By applying this approach across different metallicity bins, we find that the most metal-rich RR Lyrae stars ($[\rm Fe/H] > -0.5$) exhibit kinematics consistent with a population significantly younger ($\approx 6-7$ Gyr) than typically assumed for RR Lyrae stars. In contrast, those with $-1 < [\rm Fe/H] < -0.5$ show properties more aligned with older ($\approx 9-11$ Gyr) populations. Interestingly, we also find evidence of a possible double age populations for the most metal-rich RR Lyrae, one younger with ages between 3 and 6 Gyr, and another one older ranging from 8 to 11 Gyr. These results provide strong evidence that metal-rich RR Lyrae stars in the Galactic field do not exclusively trace ancient populations. This finding challenges the current model of RR Lyrae formation and supports alternative formation scenarios, such as binary evolution.

Modern detector manufacturing allows spectral and polarimetric filters to be directly integrated on top of separate detector pixels. This enables the creation of CubeSat-sized spectro-polarimetric instruments that are not much larger than the detector and a lens. Redundancy inherent to the observed scene, offers the opportunity for sparse sampling in the form of not scanning all filters at every location. However, when there are fewer pushbroom steps than filters, data are missing in the resulting data cube. The missing, largely redundant data can be filled in with interpolation methods, often called demosaicers. The choice of filters and their precise layout influences the performance of the instrument after the demosaicing process. In these proceedings we describe a part of a design toolbox for both the filter layout and the optimum parameters for the reconstruction to a full spectro-polarimetric data cube. The design tool is based on training a (neural) network and jointly updating the values of the filters and demosaicer. We optimized a filter layout by training on spectro-polarimetric remote observations of the Earth acquired by SPEX airborne. This optimised filter layout could reconstruct a validation scene from five overlapping snapshots (pushbroom steps), which would take 109 pushbroom steps when measuring with a classical layout and no reconstruction.

Nuclear star clusters (NSCs) are the densest stellar systems in the Universe. They can be found at the center of all galaxy types, but tend to favor galaxies of intermediate stellar mass around 10$^9\,$M$_{\odot}$[1, 2]. Currently, two main processes are under debate to explain their formation: in-situ star-formation from gas infall[3] and migration and merging of globular clusters (GCs) caused by dynamical friction[4]. Studies[5-9] of NSC stellar populations suggest that the former predominates in massive galaxies, the latter prevails in dwarf galaxies, and both contribute equally at intermediate mass. However, up to now, no ongoing merger of GCs has yet been observed to confirm this scenario. Here we report the serendipitous discovery of five dwarf galaxies with complex nuclear regions, characterized by multiple nuclei and tidal tails, using high resolution images from the Hubble Space Telescope. These structures have been reproduced in complementary N-body simulations, supporting the interpretation that they result from migrating and merging of star clusters. The small detection rate and short simulated timescales (below 100 Myr) of this process may explain why this has not been observed previously. This study highlights the need of large surveys with high resolution to fully map the migration scenario steps.

One of the most important effects shaping small-scale galaxy clustering is galaxy assembly bias, which refers to the dependence of galaxy clustering on halo properties. We investigate this effect using galaxy samples selected according to stellar mass, r-band magnitude, and broad-band colors from the largest hydrodynamical simulation of the IllustrisTNG suite. We find that galaxy assembly bias depends strongly upon the selection criteria, number density, and redshift of the sample, increasing or decreasing the clustering by as much as 25%. Interestingly, no single secondary halo property fully captures the strength of this effect for any galaxy population. Therefore, empirical approaches modeling galaxy assembly bias as a function of a single halo property cannot reproduce predictions from hydrodynamical simulations. We then study how galaxy assembly bias emerges from the interplay of halo assembly bias -- the dependence of halo clustering on properties other than mass -- and occupancy variation -- the correlation between galaxy occupation and secondary halo properties -- and provide an analytical expression that predicts the amount of galaxy assembly bias caused by any halo property. This expression facilitates understanding the dependence of galaxy assembly bias on halo properties and enables the straightforward incorporation of this effect into halo model approaches.

One of the primary objectives in modern astronomy is to discover and study planets with characteristics similar to Earth. This pursuit involves analyzing the spectra of exoplanets and searching for biosignatures. Contamination of spectra by nearby objects (e.g., other planets and moons in the same system) is a significant concern and must be addressed for future exo-Earth searching missions. The aim is to estimate, for habitable planets, the probability of spectral contamination by other planets within the same star system. This investigation focuses on the Large Interferometer for Exoplanets (LIFE). Since the Rayleigh criterion is inapplicable to interferometers such as those proposed for LIFE, we present new criteria based on the principle of parsimony, which take into account two types of issues: contamination or blending of point sources, and cancellation of point sources due to destructive interference. We define a new spatial resolution metric associated with contamination or cancellation that generalizes to a broader family of observing instruments. In the current baseline design, LIFE is an X-array architecture nulling interferometer. Our investigation reveals that its transmission map introduces the potential for two point sources to appear as one, even if they do not appear in close proximity. We find that LIFE has a spatial resolution comparable to that of a traditional telescope with a diameter of $D = 600\,\text{m}$, observing at $\lambda = 4 \,\mu\text{m}$. Our survey of a star system population shows that, out of 73.4 expected habitable planets detected, 71.3 are not contaminated on average.

Future telescopes will survey temperate, terrestrial exoplanets to estimate the frequency of habitable ($\eta_{\text{Hab}}$) or inhabited ($\eta_{\text{Life}}$) planets. This study aims to determine the minimum number of planets ($N$) required to draw statistically significant conclusions, particularly in the case of a null result (i.e., no detections). Using a Bayesian framework, we analyzed surveys of up to $N=100$ planets to infer the frequency of a binary observable feature ($\eta_{\text{obs}}$) after null results. Posterior best fits and upper limits were derived for various survey sizes and compared with predicted yields from missions like the Large Interferometer for Exoplanets (LIFE) and the Habitable Worlds Observatory (HWO). Our findings indicate that $N=20-50$ ``perfect'' observations (100\% confidence in detecting or excluding the feature) yield conclusions relatively independent of priors. To achieve 99.9\% upper limits of $\eta_{\text{obs}} \leq 0.2/0.1$, approximately $N \simeq 40/80$ observations are needed. For ``imperfect'' observations, uncertainties in interpretation and sample biases become limiting factors. We show that LIFE and HWO aim for sufficiently large survey sizes to provide statistically meaningful estimates of habitable environments and life prevalence under these assumptions. However, robust conclusions require careful sample selection and high-confidence detection or exclusion of features in each observation.

We investigate the realization of quintom scenario for dynamical dark energy within modified gravity theories that can efficiently fit the recent observational datasets. Starting from a general effective field theory formulation of dark energy in metric-affine geometry, we derive the background action in unitary gauge and we demonstrate how both $f(T)$ and $f(Q)$ gravity can naturally realize quintom behavior through appropriate forms and parameter choices. Additionally, using the Gaussian process reconstruction of the latest DESI DR2 BAO data combined with SNe and CMB observations, we extract the reconstructed dark-energy equation-of-state parameter, showing that it exhibits quintom-type evolution, crossing the phantom divide from below. Moreover, through detailed parameter estimations and application of information criteria, we compare the model with the quadratic one. Our results show that, due to its rich structure, modified gravity stands as one of the main candidates for the realization of the data-favoured dynamical dark energy.

Using cosmological hydrodynamic simulations with radiative transfer, we investigate star formation and overdensity ($\delta$) in Coma-type cluster progenitors from $z=14$ to 6. Our simulations reproduce observed $M_{\rm star}$-SFR relations and $\delta$ at these redshifts. We find: (1) protocluster (PC) and mean-density field (MF) galaxies show similar $M_{\rm star}$-SFR relations, with PC galaxies extending to higher $M_{\rm star}$ and SFR. (2) UV-bright PC galaxies ($M_{\rm UV}\lesssim -20$~mag) have $>2$ mag higher UV attenuation and shallower UV slopes than MF galaxies. (3) $\delta$ increases with redshift, depending on observational parameters (e.g., $\delta\sim50$ at $z=14$ to $\delta\sim3$ at $z=6$ for a search volume of $\sim3000$~cMpc$^3$ and a limiting magnitude of $M_{\rm UV}=-17$~mag). These results indicate that enhanced star formation in PCs is driven by massive galaxy overdensity, not anomalously high specific SFR. While simulated $\delta$ agrees with observed PC candidates (potential Coma progenitors), some MF galaxies show comparable $\delta$. We propose a robust PC identification method using both $\delta$ and $M_{\rm star}$ of the most massive member. Critical $M_{\rm star}$ thresholds for Coma progenitors are estimated ($10^{7.1}$ to $10^{10.2}$ M$_\odot$ from $z=14$ to 6). Comparison with JWST observations suggests GS-z14-0 and GS-z14-1, the current highest redshift holders, are likely progenitors of Coma-type clusters.

Hot cores represent critical astrophysical environments for high-mass star formation, distinguished by their rich spectra of organic molecular emission lines. We aim to utilize high-angular resolution molecular line data from ALMA to identify hot cores, with a particular focus on weak-emission candidates, and to provide one of the largest samples of hot core candidates. We propose to use spectral stacking and imaging techniques of complex organic molecules (COMs) in the ALMA-ATOMS survey, including line identification & weights, segmentation of line datacubes, resampling, stacking and normalization, moment 0 maps, and data analysis, to search for hot core candidates. We classify cores with dense emission of CH3OH and at least one molecule from the other six molecules as hot core candidates. In addition to the existing sample of 60 strong hot cores from the ALMA-ATOMS survey, we have detected 40 new weak candidates through stacking. All hot core candidates display compact emission from at least one of the other six COM species. For the strong sample, the stacking method provides molecular column density estimates that are consistent with previous fitting results. For the newly identified weak candidates, all species except CH3CHO show compact emission in the stacked image, which cannot be fully resolved spatially. These weak candidates exhibit column densities of COMs that are approximately one order of magnitude lower than those of the strong sample. The entire hot core sample, including the weak candidates, reveals tight correlations between the compact emission of CH3OH and other COM species, suggesting they may share a similar chemical environment for COMs, with CH3OH potentially acting as a precursor for other COMs. The molecular line stacking technique is used to identify hot core candidates in this work, leading to the identification of 40 new hot core candidates.

We present the results of our study of the long-period eclipsing binary star \Aur. The results are based on spectroscopic data obtained with the UFES échelle spectrograph and photometric observations from TESS. The derived radial velocity curve is based on 17 spectra obtained between 2021 and 2023, covering all orbital phases of this binary system. The orbital period determined from TESS data, $P = 27.019803 \pm 0.000003$ days, agrees within uncertainties with the period established in previous studies. The model constructed for the TESS photometric light curve achieves a precision of 0.01\%. The effective temperatures of both components, as well as the system metallicity, were directly derived from the spectra and are $T_\mathrm{eff, A} = 6250 \pm 50$\,K, $T_\mathrm{eff, B} = 5855 \pm 50$\,K, and $\mathrm{[Fe/H]} = -0.10 \pm 0.08$, respectively. Our analysis of the photometric and spectroscopic data allowed us to directly compute the luminosities of the components, $L_A = 1.82\,L_\odot$ and $L_B = 1.07\,L_\odot$, their radii, $R_A = 1.15\,R_\odot$ and $R_B = 1.00\,R_\odot$, and their masses, $M_A = 1.137\,M_\odot$ and $M_B = 1.023\,M_\odot$, with uncertainties below 1\%. Comparison with evolutionary tracks indicates that the system's age is $1.18 \pm 0.10$\,Gyr, and both components are still on the main sequence. The \Aur\ system is particularly interesting due to the partial eclipse of the primary component, which results in the ``inversion'' of the primary and secondary minima in the photometric light curve.

One of the observational challenges in the standard cosmological model is known as the Hubble tension. This $\sim$ 5$\sigma$ discrepancy between early and late measurements of the Hubble Constant arises from observations that rely on cosmological distance estimates, either explicitly or implicitly. In this study, we relax the assumption of the Friedmann-Lemaître-Robertson-Walker (FLRW) distance-redshift relation and explore the influence of small-scale inhomogeneities on the propagation of light from distant sources, using the Zeldovich-Kantowski-Dyer-Roeder (ZKDR) approximation as an alternative approach to address this tension. We employ the ZKDR equation along with a modified version to test our hypothesis using recent Type Ia supernovae data from the Pantheon+ compilation and the SH0ES collaboration and six gravitational lens systems from the H0LiCOW collaboration. Our findings indicate that a background model characterized by the ZKDR approximation and its modifications does not solve or alleviate the Hubble tension.

Fast Radio Bursts (FRBs) have emerged as powerful probes in cosmology. An optimized method for extracting the cosmic baryon density from localized FRBs, based on maximizing the joint likelihood function of the extragalactic dispersion measure ($\mathrm{DM}_{\mathrm{ext}}$), was proposed by Macquart et al. [Nature 581, 391 (2020)]. In this Letter, we identify a crucial term that was omitted in their derivation of the probability density function (PDF) for $\mathrm{DM}_{\mathrm{ext}}$. Using simulated FRB data, we demonstrate that neglecting this term leads to a systematic bias in the inferred cosmic baryon density, with deviations exceeding the $1\sigma$ confidence level. This highlights the importance of the missing term for the reliable cosmological application of FRBs. Furthermore, employing a sample of 88 real localized FRBs, we find that the baryon density derived using the original PDF by Macquart et al. is inconsistent with the Planck 2018 CMB data, while our corrected PDF yields a result in excellent agreement. We conclude that the omitted term is essential and must be included in order to obtain accurate cosmological constraints from FRB observations.

CONTEXT: Combining high-contrast imaging with high-resolution spectroscopy offers a powerful way to detect and characterize exoplanets around nearby stars, despite challenges linked to their faintness. Instruments like VLT/SPHERE are state of the art in high-contrast imaging, but their spectral resolution (R=50) limits them to basic characterization of close companions. These systems can detect planets down to 5-10 Mjup at 10 AU from their stars. Detection limits are mainly constrained by speckle noise, which dominates over photon and detector noise at short separations, even with advanced differential imaging. Space-based high-contrast imaging is also limited by image stability. Speckle noise can, however, be mitigated through molecular mapping, a technique that leverages high-resolution spectroscopic data. AIMS: We aim to predict detection limits in spectro-imaging after molecular mapping, analyzing how photon and detector noise propagate and comparing predictions with real data to assess performance losses from instrumental effects. We also propose mitigation strategies and validate our model using observations. METHODS: We analyzed JWST/MIRI/MRS data with FastCurves, an numerical tool, and compared results to outputs from the MIRI simulator. We also applied principal component analysis (PCA) to identify and isolate systematic effects, with and without molecular mapping. RESULTS: We studied various systematic effects and their impacts on signal and noise. PCA helped highlight and reduce straylight, fringes, and aliasing. We further compared observed and modeled companion spectra. CONCLUSIONS: FastCurves was improved to account for systematics and validated with real data. In high-flux regimes, systematics impose contrast limits even with molecular mapping. Our approach could benefit other instruments and inform the planning of future facilities like ELT/ANDES and ELT/PCS.

The Sun is the most studied of all stars, and thus constitutes a benchmark for stellar models. However, our vision of the Sun is still incomplete, as illustrated by the current debate on its chemical composition. The problem reaches far beyond chemical abundances and is intimately linked to microscopic and macroscopic physical ingredients of solar models such as radiative opacity, for which experimental results have been recently measured that still await theoretical explanations. We present opacity profiles derived from helioseismic inferences and compare them with detailed theoretical computations of individual element contributions using three different opacity computation codes, in a complementary way to experimental results. We find that our seismic opacity is about 10% higher than theoretical values used in current solar models around 2 million degrees, but lower by 35% than some recent available theoretical values. Using the Sun as a laboratory of fundamental physics, we show that quantitative comparisons between various opacity tables are required to understand the origin of the discrepancies between reported helioseismic, theoretical and experimental opacity values.

It is shown for the first time that the stripped helium stars with masses 2 to 7 solar mass which are formed in close binary systems in the so-called case B of mass-exchange and retained low-mass hydrogen-helium envelopes, experience nonlinear radial pulsations. Pulsations are excited by the kappa-mechanism due to helium ionization. The region of pulsation instability extends over Hertzsprung-Russel diagram from the red giants branch to the region of effective temperatures from about 30,000 K to about 50,000 K. Variations of stellar luminosity should be observed mostly in the ultraviolet. The amplitudes of pulsations of the studied models reach 0.8 stellar magnitude and increase, as the stellar radii decrease. Pulsation periods of stars with effective temperatures exceeding 10,000 K range from 0.17 to 8.5 day and decrease with decreasing radii. The stars have substantially larger effective temperatures than their companions, which could be Be-stars. They are components of relatively wide binaries with orbital periods up to several years. The number of pulsating moderate-mass stripped helium stars in the Galaxy is about 1000.

Detection of spectral line in gamma-ray bursts (GRBs) is importance for studying GRB physics, as it provides insights into the composition and physical conditions of the GRB environment. However, progress in detecting X-ray or gamma-ray emission and absorption lines in GRB spectra has been relatively slow, only the narrow emission line feature of about 10 MeV found in GRB 221009A has exhibited a significance exceeding $5 \sigma$. Here, we report the probable evidence of a narrow emission feature at about 2.1 mega-electron volts (MeV) in the spectrum of GRB 221023A. The highest statistical significance of this feature is observed in the time interval between 8 and 30 seconds after Fermi Gamma-Ray Burst Monitor trigger, with the chance probability value $<2.56 \times 10^{-5}$ (after accounting for the look-elsewhere effect), corresponding to a Gaussian-equivalent significance $> 4.20 \sigma$. We interpret this feature as being generated through the de-excitation of excited electrons in the relativistic hydrogen-like high-atomic-number ions entrained in the GRB jet.

We identified a new post-interaction binary, HIP 15429, consisting of a stripped star and a recently formed, rapidly rotating Be star companion ($v \sin i \approx 270$ km/s) sharing many similarities with recently identified bloated stripped stars. From orbital fitting of multi-epoch radial velocities we find a 221-day period. We also find an eccentricity of $e=0.52$, which is unexpectedly high as tides are expected to have circularised the orbit efficiently during the presumed recent mass transfer. The formation of a circumbinary disk during the mass transfer phase or the presence of an unseen tertiary companion might explain the orbit's high eccentricity. We determined physical parameters for both stars by fitting the spectra of the disentangled binary components and multi-band photometry. The stripped nature of the donor star is affirmed by its high luminosity at a low inferred mass ($\lesssim 1 \mathrm{M}_\odot$) and imprints of CNO-processed material in the surface abundances. The donor's relatively large radius and cool temperature ($T_{\mathrm{eff}} = 13.5 \pm 0.5$ kK) suggest that it has only recently ceased mass transfer. Evolutionary models assuming a 5-6 $\mathrm{M}_\odot$ progenitor can reproduce these parameters and imply that the binary is currently evolving towards a stage where the donor becomes a subdwarf orbiting a Be star. The remarkably high eccentricity of HIP 15429 challenges standard tidal evolution models, suggesting either inefficient tidal dissipation or external influences, such as a tertiary companion or circumbinary disk. This underscores the need to identify and characterise more post-mass transfer binaries to benchmark and refine theoretical models of binary evolution.

The understanding of galaxy properties and evolution is contingent on knowing the initial mass function (IMF), and yet to date, the IMF is constrained only to local galaxies. Individual stars are now becoming routinely detected at cosmological distances, where luminous stars such as supergiants in background galaxies critically lensed by galaxy clusters are temporarily further magnified by huge factors up to $10^{4}$ by intra-cluster stars, thus being detected as transients. The detection rate of these events depends on the abundance of luminous stars in the background galaxy and is thus sensitive to the IMF and the star formation history (SFH), especially for the blue supergiants detected as transients in the rest-frame UV/optical filters. As a proof of concept, we use simple SFH and IMF models constrained by spectral energy distribution (SED) to see how well we can predict the {\it HST} and {\it JWST} transient detection rate in a lensed arc dubbed ``Spock'' ($z = 1.0054$). We find that demanding a simultaneously fit of SED and rest-frame UV/optical transient detection rate places constraints on the IMF, independent of the assumed simple SFH model. We conclude our Bayesian likelihood analysis indicates that the data definitively prefer the ``Spock'' galaxy to have a Salpeter IMF ($\alpha = 2.35$) rather than a Top-heavy IMF ($\alpha = 1$) -- what is thought to be the case in the early universe -- given our methodology and assumptions with no clear excess of supergiants above the standard IMF.

Highly magnified stars ($\mu$ $>$ 100) are now outinely identified as transient events at cosmological distances thanks to microlensing by intra-cluster stars near the critical curves of galaxy clusters. Using the {\it James Webb} Space Telescope (JWST) in combination with the {\it Hubble} Space Telescope (HST), we outline here an analytical framework that is applied to the Warhol arc (at $z=0.94$) in the MACS 0416 galaxy cluster (at $z=0.396)$ where over a dozen microlensed stars have been detected to date. This method is general and can be applied to other lensed arcs. Within this lensed galaxy we fit the spatially resolved SED spanned by eight JWST-NIRCam filters combined with three ACS filters, for accurate lensed star predictions in 2D. With this tool we can generate 2D maps of microlensed stars for well resolved arcs in general, including dependence on wavelength and limiting apparent magnitude, for comparison with with planned cadenced campaigns for JWST and Hubble, for constraining directly the IMF and the level of dark matter substructure.

X-ray synchrotron radiation is expected to be highly polarized. Thanks to the Imaging X-ray Polarimetry Explorer (IXPE), it is now possible to evaluate the degree of X-ray polarization in sources such as supernova remnants (SNRs). Jointly using IXPE data and high-resolution Chandra observations, we perform a spatially resolved spectropolarimetric analysis of SNR Cassiopeia A (Cas A). We focus in the 3-6 keV energy band on regions near the shell dominated by nonthermal synchrotron emission. By combining IXPE's polarization sensitivity with Chandra's higher spatial and spectral resolution, we constrain the local polarization degree (PD) and polarization angle (PA) across the remnant. Our analysis reveals PD values ranging locally from 10% to 26%, showing significant regional variations that underscore the complex magnetic field morphology of Cas A. The polarization vectors indicate a predominantly radial magnetic field, consistent with previous studies. Thanks to the improved modeling of thermal contamination using Chandra data, we retrieve higher PD values compared to earlier IXPE analysis and more significant detections with respect to the standard IXPEOBSSIM analysis. Finally, we also estimate the degree of magnetic turbulence {\eta} from the measured photon index and PD, under the assumption of an isotropic fluctuating field across the shell of Cas A.

Understanding the mechanisms driving the escape of ionizing or Lyman continuum (LyC) emission from the interstellar medium of galaxies is necessary to constrain the evolution of Reionization, and the sources responsible for it. While progress has been made into identifying the global galaxy properties linked to the escape fraction of ionizing radiation, f$_{esc}^{LyC}$, little is currently known about how spatially resolved galaxy properties impact this parameter. We present Hubble Space Telescope (HST) imaging data obtained as part of the Lyman $\alpha$ and Continuum Origins Survey (LaCOS). LaCOS consists of HST imaging in 5 filters covering rest-frame optical and UV bands for a subsample of 42 galaxies in the Low redshift Lyman Continuum Survey, 22 being Lyman continuum emitters ($f_{esc}^{LyC}=0.01-0.49$). These data allow for investigations of the connection between sub-kpc stellar and nebular properties, including Ly$\alpha$ emission, and $f_{esc}^{LyC}$. Here, we describe the sample selection, observations and data reduction methods. Additionally, we present results on the link between global and resolved Ly$\alpha$ photometry and $f_{esc}^{LyC}$. We find similar trends between global photometric observables ($L_{Ly\alpha}$, $EW_{Ly\alpha}$, $f_{esc}^{Ly\alpha}$, $r_{50}$, $\Sigma_{SFR}$) and $f_{esc}^{LyC}$ as previously found with spectroscopy, but the correlations generally show a slightly smaller degree of correlations. However, we do find strong correlations between Ly$\alpha$ observables ($L_{Ly\alpha}$,$EW_{Ly\alpha}$) and $f_{esc}^{LyC}$ when measured in a small aperture around the brightest UV source in each galaxy. We interpret these results as evidence that LyC photons escaping on the line-of-sight are contributed by a small number of UV-bright compact regions in most galaxies in our sample.

We investigate the non-Gaussianity of second-order matter density perturbations induced by primordial gravitational waves (GWs). These tensor-induced scalar modes arise from local fluctuations in the GWs energy density, which is quadratic in tensor perturbations. The resulting second-order density contrast follows a chi-squared distribution, naturally exhibiting significant non-Gaussianity. We compute the bispectrum of these tensor-induced scalar modes and analyze its dependence on various primordial GWs power spectra, including scale-invariant, blue-tilted, Gaussian-bump, and monochromatic sources. We find that the bispectrum shape is inherently sensitive to the underlying GWs spectrum by construction. In particular, Gaussian-bump and monochromatic sources produce a strong signal peaking in the equilateral configuration, similar to the effect of scalar-induced tensor modes. Our findings reveal a new way to probe primordial GWs via galaxy surveys and highlight a unique feature of tensor-induced density perturbations, otherwise mimicking linear ones on sub-horizon scales.

We present the results of our study of ASASSN-14dx, a previously known but poorly characterised cataclysmic variable (CV). The source was observed as part of an ongoing high-time-resolution photometric survey of CVs, which revealed that, in addition to the known 82.8min orbital period, it also exhibits other transient periods, the strongest of which around 4 and 14 min. Here, we report our findings resulting from a multifaceted follow-up programme consisting of optical spectroscopy, spectropolarimetry, imaging polarimetry, and multicolour fast photometry. We find that the source displays complex optical variability, which is best explained by the presence of a massive white dwarf exhibiting non-radial pulsations. An intermediate polar-like scenario involving a spinning magnetic white dwarf can be ruled out based on the detected changes in the observed periods. Based on our optical spectroscopy, we can constrain the mass and effective temperature of the white dwarf to be ~1.1 Msol and 16 100 K, respectively. The overall intrinsic flux level of the source is unusually high, suggesting that there remains significant residual emission from the accretion disc and/or the white dwarf even ten years after the 2014 outburst. Finally, we cannot detect any spectroscopic signatures from the donor star, making ASASSN-14dx a possible period bouncer system evolving towards a longer orbital period.

The winds of massive stars remove a significant fraction of their mass, strongly impacting their evolution. As a star evolves, the rate at which it loses mass changes. In stellar evolution codes, different mass-loss recipes are employed for different evolutionary stages. The choice of the recipes is user-dependent and the conditions for switching between them are poorly defined. Focusing on hot stars, we aim to produce a physically motivated, empirically calibrated mass-loss recipe suitable for a wide range of metallicities. We want to provide a ready-to-use universal recipe that eliminates the need for switching between recipes for hot stars during stellar evolution calculations. We compile a sample of hot stars with reliable stellar and wind parameters in the Galaxy and the Magellanic Clouds. The sample is used to determine the dependence of the mass-loss rate on the basic stellar parameters. We find that independent of evolutionary stage and temperature, the wind mass-loss rate is a function of the electron-scattering Eddington parameter ($\Gamma_e$) and metallicity (Z), being in line with expectations of radiation-driven wind theory. Our derived scaling relation provides an adequate ($\Delta$log($\dot{M}$/(M$_\odot$/yr)) = 0.43) and broadly applicable mass-loss recipe for hot stars. The newly derived mass-loss recipe covers nearly the entire parameter space of hot stars with UV radiation-driven winds and eliminates the need for interpolation between mass-loss formulae at different evolutionary stages when applied in stellar evolution models. Examples of stellar evolution calculations using our new recipe reveal that the predictions on the ionizing fluxes and final fates of massive stars, especially at low metallicity, differ significantly from models that use the standard mass-loss rates, impacting our understanding of stellar populations at low metallicity and in the young Universe.

One of the current challenges in galaxy evolution studies is to establish the mechanisms that govern the escape of ionizing radiation from galaxies. In this work, we investigate the connection between Lyman Continuum (LyC) escape and the conditions of the Circumgalactic Medium (CGM), as probed by Ly$\alpha$ halos (LAHs) in emission. We use Ly$\alpha$ and UV continuum imaging data from the Lyman alpha and Continuum Origins Survey (LaCOS), targeting 42 nearby ($z \simeq 0.3$), star-forming galaxies with LyC observations (escape fractions of $f_{\rm esc}^{\rm LyC} \simeq 0.01-0.49$). LaCOS galaxies show extended Ly$\alpha$ emission ubiquitously, with LyC emitters (LCEs) having more compact Ly$\alpha$ morphologies relative to the UV size than non-LCEs, and Ly$\alpha$ spatial offsets that do not exceed the extent of the UV continuum. We model the diffuse LAHs using a combined Sersic plus exponential 2D profile, and find that the characteristic scale length of the Ly$\alpha$ is ten times the scale length of the UV, on average. We unveil a tight anti-correlation between $f_{\rm esc}^{\rm LyC}$ and the Ly$\alpha$ Halo Fraction (HF, or contribution of the halo to the total Ly$\alpha$ luminosity), that we propose as a new LyC indicator. Our observations also show that the HF scales positively with the neutral gas in the ISM, revealing a picture in which Ly$\alpha$ and LyC photons in LCEs emerge through clear sight-lines directly from the central starbursts and, in the case of Ly$\alpha$, minimizing the number of scattering interactions in the CGM. The properties of LAHs in LaCOS resemble those of LAHs at $z \geq 3$, suggesting a lack of evolution in the $f_{\rm esc}^{\rm LyC}$ predictors that rely on the spatial properties of Ly$\alpha$, and ensuring the applicability of these indicators to observations of high-redshift galaxies.

We present a differentiable, end-to-end Bayesian forward modeling framework for line intensity mapping cosmology experiments, with a specific focus on low-frequency radio telescopes targeting the redshifted 21 cm line from neutral hydrogen as a cosmological probe. Our framework is capable of posterior density estimation of the cosmological signal jointly with foreground and telescope parameters at the field level. Our key aim is to be able to optimize the model's high-dimensional, non-linear, and ill-conditioned parameter space, while also sampling from it to perform robust uncertainty quantification within a Bayesian framework. We show how a differentiable programming paradigm, accelerated by recent advances in machine learning software and hardware, can make this computationally-demanding, end-to-end Bayesian approach feasible. We demonstrate a proof-of-concept on a simplified signal recovery problem for the Hydrogen Epoch of Reionization Array experiment, highlighting the framework's ability to build confidence in early 21 cm signal detections even in the presence of poorly understood foregrounds and instrumental systematics. We use a Hessian-preconditioned Hamiltonian Monte Carlo algorithm to efficiently sample our parameter space with a dimensionality approaching $N\sim10^5$, which enables joint, end-to-end nuisance parameter marginalization over foreground and instrumental terms. Lastly, we introduce a new spherical harmonic formalism that is a complete and orthogonal basis on the cut sky relevant to drift-scan radio surveys, which we call the spherical stripe harmonic formalism, and it's associated three-dimensional basis, the spherical stripe Fourier-Bessel formalism.

High-precision corner-cube retroreflectors (CCRs) are critical for advanced lunar laser ranging (LLR) because they enable sub-millimeter-scale measurements of the Earth-Moon distance -- a level of precision essential for rigorous tests of relativistic gravitation and for advancing our understanding of lunar geophysics. In this work, we develop a comprehensive two-dimensional Fourier-optics model for single CCRs with apertures ranging from 80-110 mm. Our model incorporates realistic thermal-mechanical wavefront errors, detailed diffraction effects, and velocity aberration offsets. Our analysis reveals a strong coupling between aperture size and aberration angular offset: while larger CCRs deliver high on-axis flux under near-ideal conditions, their narrow diffraction lobes suffer significant flux loss at moderate aberration offsets, thereby favoring smaller apertures with broader main lobes. Furthermore, comparisons between solid fused-silica and hollow silicon-carbide (SiC) CCRs show that hollow designs not only achieve competitive or superior photon return -- particularly at 1064 nm, where phase errors are relatively reduced -- but also offer nearly an order-of-magnitude mass reduction for the same aperture sizes. These results establish a robust quantitative framework for optimizing CCR designs to perform at the sub-millimeter level under realistic lunar conditions and underscore the advantages of precision hollow SiC CCRs for next-generation LLR operations.

One of the oldest problems in physics is that of calculating the motion of $N$ particles under a specified mutual force: the $N$-body problem. Much is known about this problem if the specified force is non-relativistic gravity, and considerable progress has been made by considering the problem in one spatial dimension. Here, I review what is known about the relativistic gravitational $N$-body problem. Reduction to one spatial dimension has the feature of the absence of gravitational radiation, thereby allowing for a clear comparison between the physics of one-dimensional relativistic and non-relativistic self-gravitating systems. After describing how to obtain a relativistic theory of gravity coupled to $N$ point particles, I discuss in turn the two-body, three-body, four-body, and $N$-body problems. Quite general exact solutions can be obtained for the two-body problem, unlike the situation in general relativity in three spatial dimensions for which only highly specified solutions exist. The three-body problem exhibits mild forms of chaos, and provides one of the first theoretical settings in which relativistic chaos can be studied. For $N\geq 4$, other interesting features emerge. Relativistic self-gravitating systems have a number of interesting problems awaiting further investigation, providing us with a new frontier for exploring relativistic many-body systems.

We study the dynamics of gravitons in a squeezed vacuum state in a thermal radiation background. Unlike traditional treatments that rely on the Boltzmann equation, we employ the Heisenberg equation and average it over general quantum states. In contrast to the usual Boltzmann-based descriptions, our approach captures the subtleties arising from quantum coherence in different number eigenstates, which is essential for soft graviton modes in the squeezed vacuum state. Our new method successfully reproduces the previous one-loop results within the in-in formalism when the expansion parameter is small and deviates significantly as the parameter increases, indicating that our results extend beyond the one-loop in-in formalism. We examine the implications of graviton emission effects stimulated by quantum coherence in both flat and expanding backgrounds. In the flat background, it is found that backreaction of radiation on the spacetime dynamics is crucial for significant stimulated emission. In the expanding background, to avoid the subtleties associated with superhorizon modes, we investigate the effect of emission within the horizon immediately after reheating and find a significant effect. We examined the IR graviton evolution from a symmetry perspective and propose a regularization prescription to eliminate the secular growth problem.

We show that Earth's natural environment can serve as a powerful probe for ultralight axion dark matter. In the presence of global geomagnetic fields, the axions with masses ranging from $10^{-15}\,{\rm eV}-10^{-13}\,{\rm eV}$ induce electromagnetic waves in the (sub-) extremely low-frequency band ($0.3-30\,{\rm Hz}$) through the axion-photon coupling. We predict the amplitude of induced magnetic fields in the Earth-ionosphere cavity, taking the finite conductivity of the atmosphere into account. This allows us to constrain the axion-photon coupling parameter, $g_{\rm a\gamma}$, from the long-term monitoring data of the low-frequency magnetic fields, resulting in a significant improvement from the previous constraints down to $g_{\rm a\gamma} \lesssim 4\times10^{-13}\,{\rm GeV}^{-1}$ for axion mass $\sim 3 \times 10^{-14}\,{\rm eV}$.

Next-generation gravitational-wave detectors are expected to constrain the properties of extreme density matter via observations of static and dynamical tides in binary neutron star inspirals. The required modelling is straightforward in Newtonian gravity -- where the tide can be represented in terms of a sum involving the star's oscillation modes -- but not yet fully developed in general relativity -- where the mode-sum approach is problematic. As a step towards more realistic models, we are motivated to explore the post-Newtonian (pN) approach to the problem (noting that the modes should still provide an adequate basis for a tidal expansion up to 2pN order). Specifically, in this paper we develop the pN framework for neutron star oscillations and explore to what extent the results remain robust for stars in the strong-field regime. Our numerical results show that the model is accurate for low-mass stars ($\lesssim 0.8M_{\odot}$), but becomes problematic for more massive stars. However, we demonstrate that the main issues can be resolved (at the cost of abandoning the consistency of the pN expansion) allowing us to extend the calculation into the neutron star regime. For canonical neutron stars ($\approx 1.4M_\odot$) our adjusted formulation provides the fundamental mode of the star with an accuracy comparable to that of the relativistic Cowling approximation. For lower mass stars, our approach is significantly better.

A universal experimental challenge when studying radiation effects on cryogenic devices is to precisely and accurately characterize the position-dependent device response very near the energy detection threshold. We have developed a compact cryogenic optical beam steering system that can be used to generate O({\mu}s) pulses of small numbers of photons over the energy range of 1.2 - 4.5eV at room temperature, and deliver those photons via fiber optic to any specified location on the surface of a detector operating at cryogenic temperatures. This new system will allow for robust calibration of any photon-sensitive detector, including supercondcting devices. The system can be used efficiently to explore the physics of target materials, quantify the position sensitivity of different sensor designs, measure phonon transport, and study the effects of quasiparticle poisoning on detector operation. We describe the design of this pulsed calibration method and present first results obtained with a second-generation system operated at room temperature and sub-Kelvin temperatures.