Priority

Determining molecular abundances in astrophysical environments is crucial for interpreting observational data and constraining physical conditions in these regions. Chemical modelling tools are essential for simulating the complex processes that govern molecular evolution. We present SIMBA, a new Python-based single-point astrochemical modelling package designed to solve chemical reaction networks across diverse astrophysical environments. The software follows standardised rate equation approaches to evolve molecular abundances under specified physical conditions, incorporating gas-phase chemistry, grain-surface processes, and photochemistry. While leveraging Python for accessibility, performance-critical routines utilise just-in-time compilation to achieve computational efficiency suitable for research applications. A key feature of SIMBA is its graphical interface, which enables rapid investigation of chemical evolution under varying physical conditions. This makes it particularly valuable for exploring parameter dependencies and complementing more computationally intensive multi-dimensional models. We demonstrate the package's capabilities by modelling chemical evolution in a photoevaporative flow driven by external FUV irradiation. Using simplified gas dynamics, we chain multiple SIMBA instances to create a dynamic 1D model where gas evolves both chemically and dynamically. Comparing this approach to typical static models - where chemistry in each grid cell evolves independently - reveals that molecular ices, especially those with relatively high binding energies like H2O, can survive much farther into the flow than static models predict. This example case highlights how SIMBA can be extended to higher dimensions for investigating complex chemical processes. The package is open-source and includes comprehensive documentation.

Recent studies show that light seeds of black holes, which grow into massive black holes (MBHs) over time, often struggle to remain at the centers of their birthplaces in high-redshift galaxies, limiting their ability to accrete gas and merge with other black holes. In this work, we investigate how off-center formation of the first seeds affects the evolution of the MBH and massive black hole binary (MBHB) populations over cosmic history. To this end, we use the $\texttt{L-Galaxies}{\it BH}$ semi-analytical model, which includes multiple seed formation mechanisms, with light Population III remnants being the most significant contributors. To incorporate off-center formation, we modify the model to track the initial seed location, the sinking timescales toward the galactic center, and any growth during this phase. The results indicate that seed formation occurring away from the galactic center has a negligible impact on the MBH population at $z<1$, but causes significant differences at higher redshifts. Particularly, the abundance of $10^5 M_{\odot}$ MBHs at $z>4$ can be up to 2-10 times smaller compared to a nuclear seed formation model. Quasar luminosity functions with $\rm L_{bol}>10^{44} \rm erg/s$ are similarly affected, although they still align with observational constraints. The off-centre formation also alters the galaxy-MBH mass relation. At $z>5$, the amplitude of the relation can be up to 2 dex smaller than in nuclear seed models. These differences fade by $z \sim 2$ for galaxies $>10^{11} M_{\odot}$, and by $z=0$ for smaller galaxies. Notably, the overmassive MBH population recently unveiled by JWST is still present in the model, suggesting they can form independently of the seed dynamics. Finally, the merging rate of MBHs within LISA sensitivity band is strongly impacted. Specifically, there is a suppression of events at high-$z$ and an enhancement at low-$z$.

The role of interactions and mergers in the rapid quenching of massive galaxies in the early Universe remains uncertain, largely due to the difficulty of directly linking mergers to quenching. Collisional ring galaxies provide a unique opportunity, as their morphology allows precise dating of the interaction, which can then be compared to quenching timescales inferred from star formation histories. We study a gravitationally bound system at $z=1.61$ in the UDS field, composed of a Host galaxy ($M_\star = 10^{11.4} M_\odot$) with a collisional ring and an X-ray AGN, and the Bullet galaxy ($M_\star = 10^{11.2} M_\odot$), located at a projected distance of $\sim 8$ kpc. Combining JWST and HST imaging with Keck/MOSFIRE spectroscopy, we find compelling evidence for an ongoing starburst in the Host concurrent with rapid quenching in the Bullet. The ring, $\sim 20$ kpc in diameter, is expanding at $127^{+72}_{-29}$ km s$^{-1}$, implying the galaxies first collided 47–96 Myr ago. This timeline is consistent with the Host's current starburst and the Bullet's sudden quenching, strongly suggesting both phenomena were triggered by the interaction. Crucially, the Bullet shows no evidence of a preceding starburst, ruling out rapid gas consumption as the primary quenching channel. Instead, we suggest that merger-driven processes – such as enhanced turbulence and disk instabilities – may have suppressed star formation. An additional possibility, which we term the “Dragon Effect,” is that AGN-driven outflows from the Host disrupted the Bullet's low-density molecular gas, thereby preventing efficient star formation and accelerating quenching.

We investigate the X-shaped radio galaxies (XRGs) with optical double-peaked narrow emission (DPNEL) as potential hosts of dual or binary supermassive black holes (SMBHs). Using a sample of 187 XRGs selected from SDSS and DESI optical spectroscopic surveys, we check the AGN nature of both emission components using the BPT diagnostics of multiple emission lines, namely {[O III]$\lambda\lambda$4959,5007}, H$\alpha$, {[N II]$\lambda\lambda$6548,6584}, {[S II]$\lambda\lambda$6716,6731}, and H$\beta$, and mid-infrared colors. We find that the detection rate of [O III] DPNEL features in XRGs is 30% compared to just 1% in the general galaxy population (mostly radio quiet). The dual AGN fraction in DPNEL galaxies is found to depend strongly on the radio luminosity, increasing from $\sim$25% for radio-undetected to $\sim$58% in the radio-detected sample of general DPNEL galaxies. In contrast, the DPNEL XRGs and FR-II radio galaxies having higher radio power show a $\sim$95% likelihood of hosting a dual AGN. Secondly, the detection of companion galaxies in more than 30% of DPNEL XRGs suggests a vital role of mergers in the XRG formation. We also investigate the parsec-scale radio structure of the nuclei of several XRGs using Very Long Baseline Array (VLBA) maps at 1.4 GHz, 4.3 GHz or 7.6 GHz and find a resolved core for only one of the XRGs. However, the flat spectral indices of the VLBA cores along with the DPNEL components exhibiting AGN characteristics, together with the detection of radio-optical offsets between the VLBA and Gaia position, are strongly indicative of XRGs being likely candidates for hosting dual/binary AGNs.

The properties of Milky Way satellite galaxies have important implications for galaxy formation, reionization, and the fundamental physics of dark matter. However, the population of Milky Way satellites includes the faintest known galaxies, and current observations are incomplete. To understand the impact of observational selection effects on the known satellite population, we perform rigorous, quantitative estimates of the Milky Way satellite galaxy detection efficiency in three wide-field survey datasets: the Dark Energy Survey Year 6, the DECam Local Volume Exploration Data Release 3, and the Pan-STARRS1 Data Release 1. Together, these surveys cover $\sim$13,600 deg$^2$ to $g \sim 24.0$ and $\sim$27,700 deg$^2$ to $g \sim 22.5$, spanning $\sim$91% of the high-Galactic-latitude sky ($|b| \geq 15^\circ$). We apply multiple detection algorithms over the combined footprint and recover 49 known satellites above a strict census detection threshold. To characterize the sensitivity of our census, we run our detection algorithms on a large set of simulated galaxies injected into the survey data, which allows us to develop models that predict the detectability of satellites as a function of their properties. We then fit an empirical model to our data and infer the luminosity function, radial distribution, and size-luminosity relation of Milky Way satellite galaxies. Our empirical model predicts a total of $265^{+79}_{-47}$ satellite galaxies with $-20 \leq M_V \leq 0$, half-light radii of $15 \leq r_{1/2} (\rm pc) \leq 3000$, and galactocentric distances of $10 \leq D_{\rm GC} (\rm kpc) \leq 300$. We also identify a mild anisotropy in the angular distribution of the observed galaxies, at a significance of $\sim$$2\sigma$, which can be attributed to the clustering of satellites associated with the LMC.

We present 32 new spectroscopic-binary orbits from our extended radial-velocity (RV) survey of the old ($6.4 \pm 0.2$ Gyr) open cluster NGC 188. Using data from the WIYN Open Cluster Study (WOCS) and APOGEE-2, this work nearly doubles the temporal baseline of the previous RV study of NGC 188. We obtain orbital solutions within a stellar sample that spans a magnitude range of $10.8 \leq \mathrm{G} \leq 16.5 \; (0.9\mbox{-}1.2 \; {M_\odot})$. With revised membership determinations using Gaia DR3 proper-motions and parallaxes, we reassess the cluster binary frequency and period-eccentricity distribution, finding an incompleteness-corrected binary frequency of $33.1 \% \pm 3.8\%$ for periods less than $10^4$ days. We also find a tidal-circularization period of $14.4^{+0.14}_{-0.11}$ days. We find evidence that giants are deficient in short-period orbits, and suggest that the missing giants may have undergone mass transfer and in part formed the population of blue straggler stars and blue lurkers. Among the binaries of note, we highlight WOCS 3953 and WOCS 4945, long-period candidate mass-transfer systems that exhibit possible UV excess.

Spatially resolved stellar populations of massive, quiescent galaxies at cosmic noon provide powerful insights into star-formation quenching and stellar mass assembly mechanisms. Previous photometric work has revealed that the cores of these galaxies are redder than their outskirts. However, spectroscopy is needed to break the age-metallicity degeneracy and uncover the driver of these colour gradients. Here, we derive age and elemental abundance gradients for 8 distant ($1.2 \lesssim z \lesssim 2.2$), massive ($10.3\lesssim\log({\rm M}_*/{\rm M}_\odot)\lesssim 11.1$), quiescent galaxies, by fitting full-spectrum models to ultra-deep NIRSpec-MSA spectroscopy from the JWST-SUSPENSE survey. We find that these galaxies have negative age, positive [Mg/H] and [Mg/Fe], and flat [Fe/H] gradients, implying that galaxy cores are older and Mg-deficient compared to galaxy outskirts. The age gradients indicate inside-out quenching, while the Mg-deficient cores suggest rapid gas expulsion as the central quenching mechanism. Thus, galaxy cores formed faster and quenched more efficiently than their outskirts. In this scenario, however, our [Fe/H] and [Mg/Fe] gradients are still puzzling. Our results contrast lower-redshift studies, which find flat age and [Mg/Fe] gradients and negative metallicity gradients. Additionally, we find a positive trend between age gradients and rotational support, and marginal trends between gradients and galaxy velocity dispersions and ages. We discuss our findings in the context of galaxy growth scenarios, including minor mergers and progenitor bias, and the possible occurrence of different quenching mechanisms across redshift. With this work, we present the first stellar population gradients from NIRSpec-MSA spectroscopy, in the largest current sample of distant, quiescent galaxies.

Close encounters between stars and black holes can trigger micro-tidal disruption events (micro-TDEs) in dense young star clusters (YSCs). Using direct N-body simulations with PETAR, we found that most micro-TDEs arise from few-body multiple encounters. The inferred rate is approximately 350-450 Gpc$^{-3}$ yr$^{-1}$. Micro-TDEs could be detected both by upcoming surveys such as LSST, expected to observe roughly 10-100 events per year, and by their gravitational-wave (GW) signals peaking in the deci-Hertz band, detectable with future instruments such as LGWA and DECIGO.

MACS J0138-2155 is the only known cluster to strongly lens two supernovae (SNe), Requiem and Encore, from the same host galaxy at z=1.949. We present seven independent mass models of the galaxy cluster built using six software packages. By conducting a blind analysis (no exchanges of results between modeling teams), we quantified uncertainties due to modeling and software. Through HST, JWST and MUSE observations, we assembled high-quality data products, including eight "gold" lensed image systems consisting of 23 images with secure spectroscopic redshifts, and one "silver" system with a likely redshift value. Restricting to the gold images, we obtain overall consistent model predictions of the positions, magnifications and time delays of SN Encore and SN Requiem images, especially for models with $\chi^2 \leq 25$. We predict the appearance of the next images of SNe Encore and Requiem with a time delay of > 3000 days and of  3700 to 4000 days, respectively, based on a fiducial cosmological model of $H_0 = 70 {\rm\ km\ s^{-1}\ Mpc^{-1}}$ and $\Omega_{\rm m} = 0.3$. We obtain relations between $H_0$ and the time delays of SNe Encore and Requiem. In particular, for $H_0 = 73 {\rm\ km\ s^{-1}\ Mpc^{-1}}$, the four lowest $\chi^2$ models predict SN Requiem to reappear in  Apr-Dec 2026; for $H_0 = 67 {\rm\ km\ s^{-1}\ Mpc^{-1}}$, in  Mar-Nov 2027. Using the newly measured time delay between the two detected images of SN Encore by Pierel et al. (submitted) and our mass models, we jointly infer $H_0 = {\rm 66.9^{+11.2}_{-8.1}\ km\ s^{-1}\ Mpc^{-1}}$, where the uncertainty is dominated by that of the time delay. The long delays of the next-appearing SN Requiem and SN Encore images provide excellent opportunities to measure $H_0$ with an uncertainty of 2-3%. Our mass models form the basis for cosmological inference from this unique lens cluster with two strongly lensed SNe. (Abridged)

Recent Pulsar Timing Arrays (PTAs) results provided strong evidence for a stochastic gravitational wave background (sGWB), consistent with a population of merging massive black holes (MBHs) at $z<1$. Meanwhile, JWST observations at $z>5$ suggest a higher number density of accreting MBHs than previously estimated. Together with constraints from local MBHs and high-$z$ quasars, these findings offer a unique opportunity to test MBH seeding and early growth models. We explore this using ${\tt L-Galaxies}\textit{BH}$, a new extension of the galaxy formation model ${\tt L-Galaxies}$, developed to explicitly model all stages of MBH evolution, including seeding, accretion, and binary dynamics. To take advantage of both the high resolution of the ${\tt MillenniumII}$ and the large volume of the ${\tt Millennium}$ simulations, we run ${\tt L-Galaxies}\textit{BH}$ on the former and use its outputs as initial conditions for the latter, via our $\textit{grafting}$ method. We find that reproducing the number density of high-$z$ active MBHs observed by JWST requires either a heavy seed formation rate significantly higher than that predicted by current models ($\gtrsim 0.01 Mpc^{-3}$ at $z \sim 10$), or widespread formation of light seeds undergoing multiple phases of super-Eddington accretion. Furthermore, matching the amplitude of the PTA sGWB signal requires nearly all galaxies with stellar masses $M_{*}> 10^9 M_\odot$ to host central MBHs by $z\sim0$. Given the extreme heavy seed densities required to satisfy both PTA and JWST constraints, our results favor a scenario in which MBHs originate from light seeds that grow rapidly and efficiently in the early universe. This work demonstrates the power of combining multi-messenger data with physical models to probe the origins and evolution of MBHs across cosmic time.

We present Py2DJPAS, a Python-based tool to automate the analysis of spatially resolved galaxies in the \textbf{miniJPAS} survey, a 1 deg$^2$ precursor of the J-PAS survey, using the same filter system, telescope, and Pathfinder camera. Py2DJPAS streamlines the entire workflow: downloading scientific images and catalogs, performing PSF homogenization, masking, aperture definition, SED fitting, and estimating optical emission line equivalent widths via an artificial neural network. We validate Py2DJPAS on a sample of resolved miniJPAS galaxies, recovering magnitudes in all bands consistent with the catalog ($\sim 10$ \% precision using SExtractor). Local background estimation improves results for faint galaxies and apertures. PSF homogenization enables consistent multi-band photometry in inner apertures, allowing pseudo-spectra generation without artifacts. SED fitting across annular apertures yields residuals $<10$ \%, with no significant wavelength-dependent bias for regions with $S/N>5$. We demonstrate the IFU-like capability of J-PAS by analyzing the spatially resolved properties of galaxy 2470-10239 at $z = 0.078$, comparing them to MaNGA data within 1 half-light radius (HLR). We find excellent agreement in photometric vs. spectroscopic measurements and stellar mass surface density profiles. Our analysis extends to 4 HLR (S/N $\sim$ 5), showing that J-PAS can probe galaxy outskirts, enabling the study of evolutionary processes at large galactocentric distances.

The formation mechanism of Brown Dwarfs (BDs), whether akin to stars or ejected planetary-mass objects, remains debated. We present the first 3D radiation-MHD simulations of magnetized, turbulent, gravitationally unstable low-mass cores ($0.05-0.1\ \mathrm{M_{\odot}}$) collapsing into proto-BDs. Using the {\ttfamily RAMSES} code with adaptive mesh refinement, we model the full dynamical range ($10^{5}~-10^{22}\ \mathrm{cm^{-3}}$), including radiative transfer (flux limited diffusion) and non-ideal MHD (ambipolar diffusion). Our simulations self-consistently follow the isothermal collapse, first hydrostatic core formation, H$_{2}$ dissociation, and BD birth. The resulting BDs have initial radii $\approx 0.75\ \mathrm{R_{\odot}}$ and masses $\approx 0.8\ \mathrm{M_{Jup}}$, growing via accretion as we follow the early evolution of the object. Crucially, we find that BDs may form similarly to low-mass stars but with a prolonged first-core phase, supporting a star-like formation scenario.

Stellar collisions in dense galactic nuclei might play an important role in fueling supermassive black holes (SMBHs) and shaping their environments. The gas released during these collisions can contribute to SMBH accretion, influencing phenomena such as active galactic nuclei and tidal disruption events of the remnants. We address the challenge of rapidly and accurately predicting the outcomes of stellar collisionsincluding remnant masses and unbound gasacross a broad parameter space of initial conditions. Existing smoothed-particle-hydrodynamic (SPH) simulation techniques, while detailed, are too resource-intensive for exploratory studies or real-time applications. We develop a machine learning framework trained on a dataset of $\sim 16,000$ SPH simulations of main-sequence star collisions. By extracting physically meaningful parameters (e.g., masses, radii, impact parameters, and virial ratios) and employing gradient-boosted regression trees with Huber loss, we create a model that balances accuracy and computational efficiency. The method includes logarithmic transforms to handle dynamic ranges and regularization to ensure physical plausibility. The model achieves predictions of collision outcomes (remnant masses, and unbound mass) with very low mean absolute errors respect to the typical mass scale. It operates in fractions of a second, enabling large-scale parameter studies and real-time applications. Parameter importance analysis reveals that the impact parameter and the relative velocity dominate outcomes, aligning with theoretical expectations. Our approach provides a scalable tool for studying stellar collisions in galactic nuclei. The rapid predictions facilitate investigations into gas supply for SMBH accretion and the cumulative effects of collisions over cosmic time, particularly relevant to address the growth of SMBHs.

At high metallicity, a majority of massive stars have at least one close stellar companion. The evolution of such binaries is subject to strong interaction processes, heavily impacting the characteristics of their life-ending supernova and compact remnants. For the low-metallicity environments of high-redshift galaxies constraints on the multiplicity properties of massive stars over the separation range leading to binary interaction are crucially missing. Here we show that the presence of massive stars in close binaries is ubiquitous, even at low metallicity. Using the Very Large Telescope, we obtained multi-epoch radial velocity measurements of a representative sample of 139 massive O-type stars across the Small Magellanic Cloud, which has a metal content of about one fifth of the solar value. We find that 45% of them show radial velocity variations which demonstrate that they are members of close binary systems, and predominantly have orbital periods shorter than one year. Correcting for observational biases indicates that at least 70[+11:-6]% of the O stars in our sample are in close binaries, and that at least 68[+7:-8]% of all O stars interact with a companion star during their lifetime. We found no evidence supporting a statistically significant trend of the multiplicity properties with metallicity. Our results indicate that multiplicity and binary interactions govern the evolution of massive stars and determine their cosmic feedback and explosive fates.

Using the PMO 13.7m telescope, we present large-field and high-sensitivity CO(1-0) line observations toward the Crab Nebula, in order to better understand the interstellar gas environment of this well-known historical supernova remnant. The CO observations show molecular clouds toward the Crab Nebula at a velocity range from about 0 to 16 km/s. After checking the CO spectra, we find shocked signatures in the clouds extending at a velocity of roughly [5, 11] km/s. These shocked molecular clouds, with an angular distance of about 0.4-0.5 degree toward the Crab Nebula, are located at the shell of a bubble discovered in the GALFA-HI (and HI4PI) images at the same velocity range. The dimension of the bubble is roughly 2.3$\times$2.6 degree and the expansion velocity is about 5 km/s. The kinetic energy referred from the shocked molecular clouds (roughly 3.5$\times$10$^{51}$ erg), together with the HI bubble, support the picture that the Crab Nebula belongs to a typical core-collapse supernova remnant. Nevertheless, due to the large uncertainty in the distance measurement, further observations are needed to verify the physical association between the shocked molecular clouds and the Crab Nebula.

The thermodynamic properties of the intracluster medium (ICM) at the outskirts of galaxy clusters provide valuable insights into the growth of the dark matter halo and the heating of the ICM. Considering the results of the soft X-ray background study of non-cluster Suzaku fields, we revisit 65 Suzaku pointing observations of the Perseus cluster in eight azimuthal directions beyond 1 Mpc (0.8 $r_{500}$). A possible foreground component, whose spectrum is modeled as a 1 keV collisional ionization equilibrium plasma, significantly affects the temperature and density measurements of the ICM in cluster outskirts. The emission measures in the six arms are similar, showing that the radial slopes of temperature and density follow $r^{-0.67\pm0.25}$ and $r^{-2.21\pm 0.06}$, respectively. The radial pressure profile is close to the average profile measured by the Planck satellite. The resulting entropy slope is $\propto r^{0.81\pm 0.25}$, consistent with the theoretical slope of 1.1. The integrated gas fraction, the ratio of the integrated gas mass to the hydrostatic mass, is estimated to be 0.13$\pm$0.01 and 0.18$\pm$0.02 at $r_{500}$ and $r_{200}$, respectively, consistent with the cosmic baryon fraction. These results suggest that the ICM at the cluster outskirts is quite regular and close to hydrostatic equilibrium. The remaining two arms show that the emission measure is higher by a factor of 1.5-2, possibly due to accretion from filaments from the large-scale structure. A sudden drop in the emission measure also occurs in a direction toward one of the filaments.

Context. The radiation field consisting of hydrogen recombination lines and continuum emission might significantly affect the hydrogen-level populations in ultra- and hypercompact (U/HC) H II regions. The escape probability approximation was used to estimate the effect of the radiation field in previous models for calculating hydrogen-level populations. The reliability of this approximation has not been systematically studied, however. Aims. We investigate the appropriate ranges of previous models with the escape probability approximation and without the effects of the radiation field. We create a new model for simulating the integrated characteristics and the spatially resolved diagnostics of the hydrogen recombination lines throughout H II regions. Methods. We developed a new nl model with a full radiative transfer treatment of the radiation field causd by hydrogen recombination lines and continuum emission to calculate the hydrogen-level populations and hydrogen recombination lines. We then compared the level populations and the corresponding hydrogen recombination line intensities simulated by the new model and previous models. Results. We studied the applicability and the valid parameter ranges of previous models. Radiation fields exhibit negligible effects on the level populations in classical and UC H II regions. With the modified escape probability, the model with the escape probability approximation is suitable for most HC H II regions. The improved new model performs better in the HC H II region with an extremely high emission measure. To address the high computational costs inherent in numerical models, we trained a precise machine-learning model to enable a rapid estimation of hydrogen-level populations and the associated hydrogen recombination lines.

We determined HI parameters for eleven nearby late-type dwarf galaxies using FASHI data cubes, despite the fact that the first version of the FASHI catalog does not list any radio sources that could correspond to these galaxies. Four of them are probable peripheral satellites of the bright spiral galaxies: NGC 3556, NGC 4258, NGC 4274, NGC 4490, while others are isolated objects. The considered sample has the following median parameters: a heliocentric velocity of $V_\mathrm{h} = 542 \ km/s$, an HI-line width of $W_\mathrm{50} = 28 \ km/s$, an hydrogen mass of $\log (M_{HI} / M_\odot) = 6.83$, a stellar mass of $\log (M_\star / M_\odot) = 7.19$, and a specific star formation rate of $\mathrm{sSFR} = -10.17 \ yr^{-1}$.

We present a new determination of the evolving far-infrared galaxy luminosity function (FIR LF) and the resulting inferred evolution of dust-obscured star-formation rate density (SFRD) out to redshift z 6. To establish the evolving co-moving number density of FIR-bright objects, we make use of the high-resolution ALMA follow-up study (AS2UDS), of the JCMT SCUBA-2 Cosmology Legacy Survey (S2CLS) sub-mm imaging in the UKIDSS UDS survey field. In order to estimate the contributions of faint/low-mass sources we implement a method in which the faint-end of the IR LF is inferred by stacking (in stellar mass and redshift bins) the optical/near-infrared samples of star-forming galaxies into the appropriate FIR Herschel and sub-mm JCMT maps. Using this information we determine the faint-end slope of the FIR LF in two intermediate redshift bins (where it can be robustly established) and then adopt this result at all other redshifts. The evolution of the characteristic luminosity of the galaxy FIR LF, L*, is found to be increase monotonically with redshift, evolving as z^1.38+-0.07, while the characteristic number density is well fitted by double power-law function, constant at z<2.24 and declining as z^-4.95+-0.73 at higher redshifts. The evolution of the corresponding dust-obscured star-formation rate density was then calculated and is here compared with the results from a number of recent studies in the literature. Our analysis confirms that dust-obscured star-formation activity dominates SFRD at cosmic noon, but then becomes progressively less important with increasing redshift: while dusty star-forming galaxies are still found out to the highest redshifts explored here, UV-visible star formation dominates at z>4, and dust-obscured activity contributes <25% of SFRD by z 6.

We investigate how stellar disks sustain their ultrathin structure throughout their evolution. We follow the evolution of ultrathin stellar disks with varying dark matter (DM) halo concentration ($c$) using collisionless $N$-body simulations with \texttt{AREPO}. We test models embedded in steep ($c = 12$), shallow ($c = 2$), and intermediate ($c = 6$) DM concentrations. Our models match the observed structural properties of the stellar disk in the low surface brightness (LSB) ultrathin galaxy FGC 2366, specifically its surface brightness, disk scalelength, and vertical thinness ($h_{z}/R_{D} = 0.1$), while excluding gas, allowing us to isolate the effects of DM. The internal disk heating mechanism driven by bars is suppressed in the LSB ultrathin stellar disks regardless of the DM concentration. The ratio of disk thickness ($h_z$) to scalelength ($R_D$) remains constant at $\leq 0.1$ throughout their evolution. To clearly establish that the LSB nature of stellar disks is the key to preventing disk thickening, we construct the initial conditions by increasing the stellar mass fraction from $f_{s} \sim 0.01$ to $0.02$ and $0.04$, respectively, while keeping the total mass equal to $10^{11} M_\odot$ and $h_z/R_D \leq 0.1$ unchanged. We find that models with a higher stellar mass fraction embedded in a shallow DM potential ($c = 2$) form bars and undergo significant disk thickening ($h_{z}/R_{D} \gg 0.1$) concurrent with the bar growth. We conclude that if the LSB disks are thin to begin with, they remain so throughout their evolution in isolation, regardless of the concentration of the DM halo.

In cold, dense astrophysical environments dust grains are mixed with molecular ices. Chemistry in those dust/ice mixtures is determined by diffusion and reaction of molecules and radicals. However, investigations of diffusion of astrophysically relevant radicals and molecules across the surface and through the pores of cosmic dust grains and of surface reactions consequent to such diffusion is largely uncharted territory. This paper presents results of a study of a solid-state reaction of two molecular species, CO2 and NH3, separated by a layer of porous silicate grain aggregates, analogues of cosmic dust. The experiments demonstrate that the presence of the dust layer was necessary for a pure thermal CO2 + 2NH3 reaction to proceed, leading to the formation of ammonium carbamate (NH4+NH2COO-), an ionic solid containing a complex organic moiety of prebiotic interest recently detected in a protoplanetary disk. This result speaks for: (i) efficient diffusion of molecules on/within cosmic dust, (ii) an underestimated role for surface catalysis in the astrochemistry of cosmic dust, and (iii) potentially efficient dust-promoted chemistry in warm cosmic environments, such as protostellar envelopes and protoplanetary disks.

The nature and evolution of hydrocarbonaceous grains within interstellar and circumstellar media is still far from resolved, perhaps owing to the rather complex nature of their seemingly simple binary atomic compositions. This work explores the fine details of amorphous hydrocarbon nanoparticle, a-C(:H), composition and the evolution of the inherent sub-structures under extreme conditions, focusing on the characteristic CH$_n$ bands in the 3-4 micron wavelength region. Particular attention is paid to the role of dehydrogenation and its effects on the sp^3 and sp^2 hybridisations, leading to an extensive conjugated domain functionalisation of the contiguous structural network within a-C(:H) nanoparticles. Qualitatively this approach is able to explain the origin and evolution, including the appearance and disappearance, of emission bands observed in the 3-4 micron wavelength regime without a significant aromatic moiety content within the structures. A diatomic a-C(:H) phase is likely at the heart of the observed dust evolution in the interstellar medium, and circumstellar and photodissociation regions, as observed at short wavelengths. It appears that we have some way to go in fully understanding these complex materials. Much laboratory work will be required in order to elucidate their chemical and structural evolution at nanoparticle sizes under extreme conditions.

Intensity interferometry (II) offers a powerful means to observe stellar objects with a high resolution. In this work, we demonstrate that II can also probe internal stellar kinematics by revealing a time-asymmetric Hanbury Brown and Twiss (HBT) effect, causing a measurable shift in the temporal correlation peak away from zero delay. We develop numerical models to simulate this effect for two distinct astrophysical scenarios: an emission-line circumstellar disk and an absorption-line binary system. Our simulations reveal a clear sensitivity of this temporal asymmetry to the system's inclination angle, velocity symmetry, and internal dynamics. This suggests that, with sufficiently high time resolution, II can be used to extract quantitative information about internal kinematics, offering a new observational window on stellar dynamics.

Local Universe dwarf galaxies are both cosmological and mass assembly probes. Deep surveys have enabled the study of these objects down to the low surface brightness (LSB) regime. In this paper, we estimate Euclid's dwarf detection capabilities as well as limits of its MERge processing function (MER pipeline), responsible for producing the stacked mosaics and final catalogues. To do this, we inject mock dwarf galaxies in a real Euclid Wide Survey (EWS) field in the VIS band and compare the input catalogue to the final MER catalogue. The mock dwarf galaxies are generated with simple Sérsic models and structural parameters extracted from observed dwarf galaxy property catalogues. To characterize the detected dwarfs, we use the mean surface brightness inside the effective radius SBe (in mag arcsec-2). The final MER catalogues achieve completenesses of 91 % for SBe in [21, 24], and 54 % for SBe in [24, 28]. These numbers do not take into account possible contaminants, including confusion with background galaxies at the location of the dwarfs. After taking into account those effects, they become respectively 86 % and 38 %. The MER pipeline performs a final local background subtraction with small mesh size, leading to a flux loss for galaxies with Re > 10". By using the final MER mosaics and reinjecting this local background, we obtain an image in which we recover reliable photometric properties for objects under the arcminute scale. This background-reinjected product is thus suitable for the study of Local Universe dwarf galaxies. Euclid's data reduction pipeline serves as a test bed for other deep surveys, particularly regarding background subtraction methods, a key issue in LSB science.

Boxy/peanut and X-shaped (BP/X) bulges are prominent features in edge-on disk galaxies and are believed to be vertically thickened bars. Despite their relevance in bar evolution, a statistically robust census of these structures in large surveys has been lacking. We aim to provide the largest catalog of BP/X structures in edge-on galaxies to date, and to investigate their properties and role in shaping galaxy scaling relations. We selected a sample of 6684 edge-on galaxies from SDSS DR8 using Galaxy Zoo classifications, requiring a high edge-on probability ($> 0.9$) and a minimum of 10 independent votes. Two-dimensional image decomposition is performed using GALFIT to obtain structural parameters. Residual images are visually inspected to classify BP/X features into four categories: strong both-sided, both-sided, one-sided, and control (no BP/X). We also estimated stellar mass, distance, and physical size for each galaxy. Out of 6653 classified galaxies, we identified 1675 ($\sim$25%) with both-sided BP/X features-545 ($\sim$8%) strong and 1130 ($\sim$17%) faint-as well as 1108 ($\sim$17%) one-sided structures, making up a total of 2783 BP/X-hosting galaxies ($\sim$42%). One-sided structures, likely signatures of ongoing buckling, are more frequent than strong both-sided bulges across all stellar masses. The fraction of BP/X bulges increases with stellar surface mass density, indicating a connection with bar formation in dense disks. We also find that galaxies with strong BP/X bulges contribute to increased scatter in the stellar mass-size and stellar mass-surface density relations, particularly at higher masses.

The Vera C. Rubin Observatory LSST is expected to discover tens of millions of new Active Galactic Nuclei (AGNs). The survey's exceptional cadence and sensitivity will enable UV/optical/NIR monitoring of a significant fraction of these objects. The unprecedented number of sources makes spectroscopic follow-up for the vast majority of them unfeasible in the near future, so most studies will have to rely on photometric redshifts estimates which are traditionally much less reliable for AGN than for inactive galaxies. This work presents a novel methodology to constrain the photometric redshift of AGNs that leverages the effects of cosmological time dilation, and of the luminosity and wavelength dependence of AGN variability. Specifically, we assume that the variability can be modeled as a damped random walk (DRW) process, and adopt a parametric model to characterize the DRW timescale ($\tau$) and asymptotic amplitude of the variability (SF$_\infty$) based on the redshift, the rest-frame wavelength, and the AGN luminosity. We construct variability-based photo-$z$ priors by modeling the observed variability using the expected DRW parameters at a given redshift. These variability-based photometric redshift (VAR-PZ) priors are then combined with traditional SED fitting to improve the redshift estimates from SED fitting. Validation is performed using observational data from the SDSS, demonstrating significant reduction in catastrophic outliers by more than 10% in comparison with SED fitting techniques and improvements in redshift precision. The simulated light curves with both SDSS and LSST-like cadences and baselines confirm that, VAR-PZ will be able to constrain the photometric redshifts of SDSS-like AGNs by bringing the outlier fractions down to below 7% from 32% (SED-alone) at the end of the survey.

Using NIRSpec on JWST, we studied a sample of 15 intermediate-mass (1.8-4.1 Msun) young stellar objects (YSOs) previously identified with MIRI photometry in the low-metallicity NGC 346 star-forming cluster in the Small Magellanic Cloud (SMC). All objects, observed in the 1.7-5.3 micron range, show strong hydrogen recombination lines in the Paschen, Brackett, Pfund, and Humphreys series, confirming their very young ages. The spectra of 11 YSOs show prominent absorption bands from the three most important ice species (H2O, CO2, CO), marking the first detection of these ices in intermediate-mass YSOs beyond our Galaxy. In three YSOs, water ice appears to be in crystalline form. In some objects, we also detect 13CO2 and OCS ices – never before observed beyond the Milky Way (MW) – and methanol ice in at least one star. We compared the column densities of H2O, CO2, and CO ices with those measured in more and less massive protostars in the MW and Large Magellanic Cloud, finding that in NGC 346 ice column densities reach values nearly an order of magnitude lower than in more massive objects ( 1x10^{17} cm-2 for water and  1x10^{16} cm-2 for CO2 and CO). However, the relative proportions of the ice species abundances do not differ from those in massive MW YSOs. This suggests that metallicity may not significantly affect ice chemistry in protoplanetary discs and that, shielded by the protostellar envelope or deep in the midplane, circumstellar material is likely impervious to the radiation environment.

Mass segregation is seen in many star clusters, but whether massive stars form in the center of a cluster or migrate there dynamically is still debated. N-body simulations have shown that early dynamical mass segregation is possible when sub-clusters merge to form a dense core with a small crossing time. However, the effect of gas dynamics on both the formation and dynamics of the stars could inhibit the formation of the dense core. We aim to study the dynamical mass segregation of star cluster models that include gas dynamics and self-consistently form stars from the dense substructure in the gas. Our models use the Torch framework, which is based on AMUSE and includes stellar and magnetized gas dynamics, as well as stellar evolution and feedback from radiation, stellar winds, and supernovae. Our models consist of three star clusters forming from initial turbulent spherical clouds of mass $10^{4,5,6}\rm~M_\odot$ and radius $11.7\rm~pc$ that have final stellar masses of $3.6\times10^3\rm~M_\odot$, $6.5\times10^4\rm~M_\odot$, and $8.9\times10^5\rm~M_\odot$, respectively. There is no primordial mass segregation in the model by construction. All three clusters become dynamically mass segregated at early times via collapse confirming that this mechanism occurs within sub-clusters forming directly out of the dense substructure in the gas. The dynamics of the embedded gas and stellar feedback do not inhibit the collapse of the cluster. We find that each model cluster becomes mass segregated within $2~$Myr of the onset of star formation, reaching the levels observed in young clusters in the Milky Way. However, we note that the exact values are highly time-variable during these early phases of evolution. Massive stars that segregate to the center during core collapse are likely to be dynamically ejected, a process that can decrease the overall level of mass segregation again.

Recent observations suggest that the extended stellar halos of low-redshift massive galaxies are tightly connected to the assembly of their dark matter halos. In this paper, we use the Illustris, IllustrisTNG100, and IllustrisTNG300 simulations to compare how different stellar aperture masses trace halo mass. For massive central galaxies ($M_\star\geq 10^{11.2}M_\odot$), we find that a 2D outskirt stellar mass measured between 50 to 100 kpc ($M_{\star,[50,100]}$) consistently outperforms other aperture-based stellar masses. We further show that $M_{\star,[50,100]}$ correlates better with halo mass than the total amount of accreted stars (the ex situ mass), which suggests that not all accreted stars connect to halo assembly equally. While the galaxy formation recipes are different between Illustris and IllustrisTNG100, the two simulations yield consistent ex situ outskirt fractions for massive galaxies (about 70% in $M_{\star,[50,100]}$). These results demonstrate the potential of using the outskirt stellar mass to deepen our understanding of galaxy-halo connection in massive dark matter halos and trace dark matter halos better.

Synergies between JWST and ALMA are unveiling a population of bright, super-early ($z>10$) galaxies, including systems like GS-z14-0 ($z=14.2$) and GHZ2 ($z=12.3$) with extreme FIR line ratios ([OIII] 88$\,$um / [CII] 158$\,$um $>3$) that challenge galaxy formation models. To address this, we identify a synthetic analog of these sources, "Amaryllis", within the SERRA zoom-in simulations, and track its evolution from $z=16$ to $z=7$. During this period, Amaryllis grows from $\log(M_\star/M_{\odot}) \sim 7.4$ to $10.3$, linking super-early progenitors to the massive galaxy population at the end of reionization. At $z \sim 11.5$, Amaryllis closely matches the observed properties of GS-z14-0, including $M_\star$, SFR, and the luminosity of FIR ([OIII] 88$\,$um) and UV (e.g. CIII]$\,1908$) lines. We find that high [OIII]/[CII] ratios appear during short, merger-driven starburst episodes, when low metallicity ($Z \sim 0.02\,Z_{\odot}$) and high ionization conditions ($U_{\mathrm{ion}} \sim 0.1$) push the ISM far from equilibrium. These extreme FIR line ratios are thus transient and linked to major mergers that ignite strong ionized gas outflows. Strikingly, despite this dynamical activity, Amaryllis develops a rotation-supported gaseous disk ($V/\sigma \sim 4$-6) by $z \sim 11$, while stars remain dispersion-dominated. This coexistence of ordered gas rotation and merger-driven disturbances occurs within a massive yet typical $\Lambda$CDM halo, enabling disk formation even at cosmic dawn.

We investigate how AGN disk turbulence affects the orbital dynamics of a stellar-mass black hole (BH) initially located at a migration trap, focusing on the long-term behavior of eccentricity and inclination in the quasi-embedded regime. We develop a semi-analytical framework in which turbulence is modeled as a stochastic velocity field acting through a modified drag force. We integrate the resulting stochastic differential equations both in Cartesian coordinates and in orbital elements using a linearized perturbative approach, and compare these results with full numerical simulations. Eccentricity and inclination evolve toward steady-state Rayleigh distributions, with variances determined by the local disk properties and the ratio of the gas damping rate to the orbital frequency. The analytical predictions agree well with the numerical simulations. We provide closed-form expressions for the variances in both the fast and slow damping regimes. These results are directly applicable to Monte Carlo population models and can serve as physically motivated initial conditions for hydrodynamical simulations. Turbulent forcing prevents full circularization and alignment of BH orbits in AGN disks, even in the presence of strong gas drag. This has important implications for BH merger and binary formation rates, which are sensitive to the residual eccentricity and inclination. Our results highlight the need to account for turbulence-induced stochastic heating when modeling the dynamical evolution of compact objects in AGN environments.

In the era of JWST, with its unprecedented sensitivity and spectral resolution, infrared spectral surveys have revealed a rich inventory of ices, including complex organic molecules (COMs), in young stellar objects (YSOs). However, robust methods to decompose and quantify these absorption features particularly across broad spectral ranges, are still under investigation. We present INDRA (Ice-fitting with NNLS-based Decomposition and Retrieval Algorithm), a fully Python-based tool that performs continuum and silicate removal, global ice fitting using Weighted Non-Negative Least Squares (NNLS), and estimates column densities and statistical significance. We apply INDRA to NGC 1333 IRAS 2A, a target from the JWST Observations of Young protoStars (JOYS+) program previously studied using local fitting. We derive optical depths via polynomial continuum subtraction and remove silicate absorption using a synthetic model, isolating ice features for global MIRI fitting. Our results are consistent with previous local fits, confirming simple species and COMs, and expand the inventory by identifying additional absorption features from CO2 and NH4+. We also propose the presence of organic refractories contributing up to 9.6% in the spectral region of 5-8 microns among the various ice components, whose inclusion significantly improves the global spectral fitting. These broad absorption features, extending across 5.5-11 microns, are likely produced by large, complex molecules containing carbonyl (C=O), hydroxyl (O-H), amine (N-H), and C-H bending modes. Our expanded inventory, now incorporating these organic residues, offers new insights into the chemical evolution of ices in star-forming regions and highlights the importance of global spectral fitting in constraining ice compositions.

The quenching of star formation in galaxies is an important aspect of galaxy evolution, but the physical mechanisms that drive it are still not understood. Measuring the spatial distribution of quenching can help determine these mechanisms. We present the star-formation histories (SFHs) and stellar metallicity evolution of rapidly quenched regions in 86 local post-starburst (PSB) galaxies from the MaNGA integral field survey, obtained through Bayesian full spectral fitting of their rest-frame optical spectra. We found that regardless of spatial location, the PSB regions have similar past SFHs and chemical evolution, once radial metallicity gradients are accounted for. This suggests that all PSB regions are regulated by a common set of local scale processes in the interstellar medium, regardless of the broader triggering mechanism. We show that the centres of galaxies with outer PSB regions are also quenching. The central specific star-formation rate (sSFR) has declined by $\sim1.0\;$dex on average during the last 2 Gyr, a significantly steeper decline than main sequence galaxies over the same period ($\approx0.2\;$dex). This central quenching can be either synchronous, outside-in or inside-out, and slower or as fast as the outer regions, highlighting the diversity of quenching pathways for local galaxies. Our results imply a primary quenching mechanism that is both catastrophic and global in rapidly halting star formation in local galaxies. We suggest the predominant cause is galaxy mergers or interactions, with large scale feedback from a starburst or a central supermassive black hole playing a lesser role.

It is commonly supposed that quenched field dwarfs near the edge of the Local Group (LG) are backsplash galaxies, having previously orbited within the Milky Way (MW) or M31's virial radius, whereas galaxies on first infall should still have gas and star formation. We measured proper motions (PMs) for six dwarf galaxies located 400-1000 kpc from the MW using the Hubble Space Telescope. For four galaxies (Aquarius, Cetus, Pisces, Tucana), we report the first PMs. For the remaining two (Leo T and Pegasus), we measure PMs with order-of-magnitude improvement. We compute orbital histories to assess whether any of the six are backsplash galaxies. While some have non-zero likelihoods of past interaction with the MW or M31, these are weak and typically occur at large distances (e.g., $>$ 2Rvir). The properties of Aquarius, Leo T, Pisces, and Pegasus are consistent with first passage through a massive halo. Cetus, which shows a low probability ( 4-6%) of interacting with the MW or M31 in the last 6 Gyr, is more likely a backsplash galaxy resulting from an interaction with M31 over 6 Gyr ago, in the same regime where rigid orbital models become less reliable. Tucana has been thought to be a backsplash galaxy, but our orbits indicate it cannot have interacted with a massive LG host. Our results highlight the diversity of evolutionary pathways for isolated, intermediate-mass dwarfs ($M_* \approx 10^5-10^7 M_{\odot}$) and the need to reassess quenching mechanisms beyond environmental interactions with massive hosts.

Dark matter (DM) remains one of the most compelling unresolved problems in fundamental physics, motivating the search for new detection approaches. We propose a network-based quantum sensor architecture to enhance sensitivity to ultralight DM fields. Each node in the network is a superconducting qubit, interconnected via controlled-Z gates in symmetric topologies such as line, ring, star, and fully connected graphs. We investigate four- and nine-qubit systems, optimizing both state preparation and measurement using a variational quantum metrology framework. This approach minimizes the quantum and classical Cramér-Rao bounds to identify optimal configurations. Bayesian inference is employed to extract the DM-induced phase shift from measurement outcomes. Our results show that optimized network configurations significantly outperform conventional GHZ-based protocols while maintaining shallow circuit depths compatible with noisy intermediate-scale quantum hardware. Sensitivity remains robust under local dephasing noise. These findings highlight the importance of network structure in quantum sensing and point toward scalable strategies for quantum-enhanced DM detection.

We present \textbf{VADER} (Variational Autoencoder for Disks Embedded with Rings), for inferring both planet mass and global disk properties from high-resolution ALMA dust continuum images of protoplanetary disks (PPDs). VADER, a probabilistic deep learning model, enables uncertainty-aware inference of planet masses, $\alpha$-viscosity, dust-to-gas ratio, Stokes number, flaring index, and the number of planets directly from protoplanetary disk images. VADER is trained on over 100{,}000 synthetic images of PPDs generated from \texttt{FARGO3D} simulations post-processed with \texttt{RADMC3D}. Our trained model predicts physical planet and disk parameters with $R^2 > 0.9$ from dust continuum images of PPDs. Applied to 23 real disks, VADER's mass estimates are consistent with literature values and reveal latent correlations that reflect known disk physics. Our results establish VAE-based generative models as robust tools for probabilistic astrophysical inference, with direct applications to interpreting protoplanetary disk substructures in the era of large interferometric surveys.

Near-Earth asteroid 2024 YR4 was discovered on 2024-12-27 and its probability of Earth impact in December 2032 peaked at about 3% on 2025-02-18. Additional observations ruled out Earth impact by 2025-02-23. However, the probability of lunar impact in December 2032 then rose, reaching about 4% by the end of the apparition in May 2025. James Webb Space Telescope (JWST) observations on 2025-03-26 estimated the asteroid's diameter at 60 +/- 7 m. Studies of 2024 YR4's potential lunar impact effects suggest lunar ejecta could increase micrometeoroid debris flux in low Earth orbit up to 1000 times above background levels over just a few days, possibly threatening astronauts and spacecraft. In this work, we present options for space missions to 2024 YR4 that could be utilized if lunar impact is confirmed. We cover flyby & rendezvous reconnaissance, deflection, and robust disruption of the asteroid. We examine both rapid-response and delayed launch options through 2032. We evaluate chemical and solar electric propulsion, various launch vehicles, optimized deep space maneuvers, and gravity assists. Re-tasking extant spacecraft and using built spacecraft not yet launched are also considered. The best reconnaissance mission options launch in late 2028, leaving only approximately three years for development at the time of this writing in August 2025. Deflection missions were assessed and appear impractical. However, kinetic robust disruption missions are available with launches between April 2030 and April 2032. Nuclear robust disruption missions are also available with launches between late 2029 and late 2031. Finally, even if lunar impact is ruled out there is significant potential utility in deploying a reconnaissance mission to characterize the asteroid.

Mitigation of the threat from airbursting asteroids requires an understanding of the potential risk they pose for the ground. How asteroids release their kinetic energy in the atmosphere is not well understood due to the rarity of significant impacts. Ordinary chondrites, in particular L chondrites, represent a frequent type of Earth-impacting asteroids. Here, we present the first comprehensive, space-to-lab characterization of an L chondrite impact. Small asteroid 2023 CX1 was detected in space and predicted to impact over Normandy, France, on 13 February 2023. Observations from multiple independent sensors and reduction techniques revealed an unusual but potentially high-risk fragmentation behavior. The nearly spherical 650 $\pm$ 160 kg (72 $\pm$ 6 cm diameter) asteroid catastrophically fragmented around 28 km altitude, releasing 98% of its total energy in a concentrated region of the atmosphere. The resulting shockwave was spherical, not cylindrical, and released more energy closer to the ground. This type of fragmentation increases the risk of significant damage at ground level. These results warrant consideration for a planetary defense strategy for cases where a >3-4 MPa dynamic pressure is expected, including planning for evacuation of areas beneath anticipated disruption locations.

High-resolution spectroscopic measurements of OB stars are important for understanding processes like stellar evolution, but require labor-intensive observations. In contrast, photometric missions like the Transiting Exoplanet Survey Satellite (TESS) can monitor hundreds of thousands of stars with a range of temporal resolutions, but do not provide such detailed measurements. With surveys like the Legacy Survey of Space and Time promising unprecedented photometric coverage over the next ten years, it is increasingly important to develop methods that connect large-scale time-series photometry with the detailed stellar parameter measurements typically derived from spectroscopy. In this paper, we test whether machine learning can recover such parameters by combining TESS light curves with spectroscopic measurements from the IACOB project, using a sample of 285 light curves from 106 unique O stars. Using both multilayer perceptrons and convolutional neural networks, we demonstrate that (1) O star light curves contain sufficient information to meaningfully infer stellar parameters and (2) periodograms derived from light curves capture substantially more information than previously identified correlation parameters. Our best model achieves moderate success in predicting both spectroscopic luminosity ($R^2 = 0.641_{-0.167}^{+0.107}$) and effective temperature ($R^2 = 0.443_{-0.234}^{+0.056}$), key stellar parameters for determining positions of stars on the spectroscopic Hertzsprung-Russell diagram, despite the small dataset size. Further progress will require expanded datasets of matched photometric and spectroscopic observations.

Large-scale research infrastructures (LSRIs) are central to contemporary science policy, combining massive capital investments with international access regimes. Yet whether open access to these infrastructures translates into more equitable scientific authority remains contested. Astronomy provides a critical case: world-leading observatories are globally shared but embedded in specific national contexts. We compile a novel country–year dataset (1955–2025) linking the location of astronomical facilities with publication usage and authorship roles. This enables us to distinguish between hosting, using, and leading in telescope-based research. Our analysis reveals: (i) usage and impact are heavily concentrated in a small number of facility hubs; (ii) scientific leadership is even more unequal than access or usage (Gini coefficient 0.91 for first/corresponding authorship versus 0.85 for facilities and usage); (iii) hosting and leadership often decouple–countries such as Chile and South Africa mediate large publication volumes without commensurate gains in leading roles; and (iv) global leadership has shifted from U.S. dominance to a multi-hub system centered in the United States, Western Europe, China, Japan, and Australia. These findings challenge the assumption that international access alone democratizes science. We argue that converting participation into leadership requires domestic PI programs, investments in instrumentation and data pipelines, and governance models that distribute credit more equitably. The study highlights how the governance of LSRIs shapes global scientific hierarchies and offers design principles for infrastructures that seek not only to share data but also to broaden scientific authority.

Periodic orbits (POs) play a central role in the circular restricted three-body problem (CRTBP). This paper introduces a method to search for POs by identifying single- and multiple-revolution fixed points in chosen Poincare maps that describe the CRTBP dynamics, with a theoretical capability to detect all fixed points across arbitrary revolution counts this http URL , high-order transfer maps (HOTMs), represented as polynomials, are constructed within the differential algebra (DA) framework for both planar and spatial CRTBP to map states between successive Poincare section crossings, with the Jacobi constant used to reduce the number of independent variables. Next, an automatic domain splitting (ADS) strategy is employed to generate subdomains, preserving HOTM accuracy, with an integrated feasibility estimation to reduce ADS's computation this http URL , a two-stage HOTM-based polynomial optimization framework is introduced, first identifying combinable subdomain sequences and then refining the fixed point solutions. Finally, the method is applied to the Earth-Moon CRTBP, identifying POs up to nine revolutions in the planar case and four in the spatial case. Known families such as distant retrograde orbits (DROs) and Lyapunov orbits are recovered, along with a previously undocumented family that exhibits a hybrid character between DROs and Lyapunov orbits.

Our ability to observe, detect, and characterize exoplanetary atmospheres has grown by leaps and bounds over the last 20 years, aided largely by developments in astronomical instrumentation; improvements in data analysis techniques; and an increase in the sophistication and availability of spectroscopic models. Over this time, detections have been made for a number of important molecular species across a range of wavelengths and spectral resolutions. Ground-based observations at high resolution are particularly valuable due to the high contrast achievable between the stellar spectral continuum and the cores of resolved exoplanet absorption features. However, the model-independent retrieval of such features remains a major hurdle in data analysis, with traditional methods being limited by both the choice of algorithm used to remove the non-exoplanetary components of the signal, as well as the accuracy of model template spectra used for cross-correlation. Here we present a new algorithm TSD (Transmission Spectroscopy Decomposition) formulated as an inverse problem in order to minimize the number of assumptions and theoretically modelled components included in the retrieval. Instead of cross-correlation with pre-computed template exoplanet spectra, we rely on high spectral resolution and instrument stability to distinguish between the stellar, exoplanetary, and telluric components and velocity frames in the sequence of absorption spectra taken during multiple transits. We demonstrate the performance of our new method using both simulated and real K band observations from ESO's VLT/CRIRES+ instrument, and present results obtained from two transits of the highly-inflated super-Neptune WASP-107 b which orbits a nearby K7V star.

Resolving high-contrast targets is a fundamental yet highly challenging task in astronomy. Using quantum estimation theory, we demonstrate that the ultimate limit for estimating the separation between two unequal-brightness thermal sources via interferometry remains constant, enabling the potential for superresolution. We give a comparative analysis of two primary stellar interferometric schemes: amplitude interferometry and intensity interferometry. Notably, the nulling strategy employed in amplitude interferometry, a configuration specifically proposed for exoplanet detection by leveraging destructive interference to suppress the brighter source, is quantum optimal for separation estimation. While intensity interferometry is less effective than amplitude interferometry in lossless scenarios and fails to achieve superresolution, it becomes competitive when optical loss in large-scale interferometry is considered. By applying these methodologies to modern stellar interferometry, we highlight the promise of large-scale interferometry for advancing high-resolution astronomical observation.

The performance of future observatories such as the Extremely Large Telescope is mainly limited by atmospheric turbulence and structural vibrations of the optical assembly. To further enhance the mitigation performance of adaptive optics, real-time information about the disturbances acting on the control loop is needed. Current systems therefore employ a combination of wavefront sensor- and accelerometer-based filters. In this work, methods using only data from natural- and laser guide star (NGS, LGS) measurements are presented, as telescopes like the Very Large Telescope already have multiple fast and high-resolution wavefront sensors installed. This approach also avoids the costly installation and operation of additional accelerometers on the optical elements. We introduce two innovative disturbance observer schemes to sense both turbulence and vibration information. A multi-rate estimator for atmospheric influences is based on Kalman filter theory and can incorporate NGS and LGS signals at different loop rates. The estimator for structural perturbations uses Gaussian process regression and can be implemented in an offline and online configuration. We validate the filter designs with data from a realistic end-to-end adaptive optics model with randomly generated turbulence and vibrations. The simulation is fed with on-sky data from the Adaptive Optics Facility of the Very Large Telescope. The presented disturbance observer schemes demonstrate promising results and may be considered as potential alternatives or extensions to existing techniques such as linear-quadratic controllers with Kalman filtering (LQG).

Recent advances in ground-based astronomy have made it possible to create optical telescopes with primary mirrors up to 40 m in size. With growing mirror diameter, the suppression of non-atmospheric disturbances becomes increasingly important. Precise knowledge of the movement of telescope mirrors is essential for understanding and compensating for vibration-based perturbations. A model from VLT accelerometer data for each individual mirror is developed, while the influence of wind buffeting is accounted for by a von Karman wind model. To describe the relevant rigid body motion, we consider the piston, tip and tilt modes of the mirrors. The identification is validated by comparing the power spectral density of the measured and identified modes. Additionally, we assess the robustness of the approach by calculating the identification error over different sections of the data. The study indicates that the employed methods are adequate for the identification of modal telescope vibrations. It is anticipated that said findings will serve as a significant foundation for the development of advanced model-based AO controllers for large telescopes, such as linear quadratic Gaussian control.

The goal of this study is to analyze the photometric properties of Deimos using Mars Express (MEX) observations, to improve the photometric properties and provide new insights into the texture and composition of the surface of Deimos, in preparation for the MMX mission. We analyzed the data obtained by the HRSC and the SRC cameras onboard MEX. The HRSC data, obtained through the use of four filters (blue, green, red, IR), provides 390 to 800 m/px resolution, while the SRC data reach 85 to 300 m/px and cover a wide phase angle range (0.06-138°). We performed the disk-integrated and disk-resolved photometric analysis using the Hapke model. The Deimos surface is dark and predominantly backscattering, with a single-scattering albedo (SSA) value (6.8%-7.5%) comparable to Phobos. The Deimos phase curve shows a strong opposition effect due to shadow-hiding, with negligible coherent backscattering. The amplitude and the half-width of the shadow-hiding opposition surge were found to be 2.14 +/- 0.14 and 0.065 +/- 0.004, respectively. We found a high porosity of 86% at the top-layer surface, consistent with complex-shaped grains or fractal aggregates, suggesting a thick dust layer. We did not observe significant variations of the opposition surge across the surface. A blue unit on Deimos, located on streamers of the equatorial ridge, shows reflectance increases up to 58%, and a spectral slope decrease of 50% in comparison with the average surface. This blue unit may be due to a different texture of the surface between the two units, with finer grain and/or higher porosity. Deimos photometric properties, including SSA, opposition surge, and phase integral, are very similar to Phobos. The presence of a blue unit on Deimos reinforces the idea that the Martian moons have a common origin, making the capture of two different bodies with such similar properties unlikely.

An international collaboration composed of Italian, Japanese, Spanish and Swiss institutes, is developing the advanced camera (AdvCam), the next-generation camera for Imaging Atmospheric Cherenkov Telescopes, designed specifically for the Large-Sized Telescopes (LST) of the Cherenkov Telescope Array Observatory. AdvCam incorporates cutting-edge Silicon Photomultipliers (SiPMs) and a fully digital readout system, setting new standards for performance and efficiency. The upgraded camera will feature four times more pixels for the same field of view as the existing PMT-based camera, enabling finer image resolution and significantly improving angular precision and background noise rejection. To cope with the increase in number of pixels, many technological challenges are being tackled, from low power and high speed integrated chip design to real-time data processing on hardware accelerators. This technological leap will lower the energy threshold by allowing operation at lower observation threshold and providing brighter images. The increase in effective area, angular and energy resolution will enhance the sensitivity, unlocking new potential for gamma-ray astronomy. In this work, we present the performance of the AdvCam's core building blocks and its innovative architecture capable of enabling unprecedented triggering capabilities. We also showcase the latest performance results based on Monte-Carlo data that has been tuned to reflect the latest stages of the on-going technological developments, highlighting the transformative capabilities of this next-generation instrument.

Current models of binary systems often depend on simplified approach of the radiation field, which are unlikely to accurately capture the complexities of asymmetric environments. We investigate the dynamical and chemical implications of a 3D asymmetric radiation field that accounts for the optical properties of sub-structures present in a protoplanetary disk, as well as the inclusion of a secondary radiation source in binary systems. We conducted a series of 3D-SPH hydrodynamical simulations using PHANTOM, coupled with the 3D Monte Carlo radiative transfer code MCFOST, to compute disc temperatures on-the-fly. We explored different binary-disk orientations (0$^o$ and 30$^o$) for an eccentric binary, along with a constant dust-to-gas ratio and dust as a mixture prescription. We also simulated an outburst event as an example of a drastic increase in luminosity. Heating from the secondary star inflates the outer disk, increasing the aspect ratio facing the companion by about 25% in inclined cases compared to 10% in coplanar ones. Dust settling in the mid-plane enhances extinction along the disk plane, making the coplanar case cooler than the inclined one on the side of the disk facing the companion. Besides, heating causes a shift in the snow line for species with freeze-out temperatures below 50 K, depending on the disk-binary inclination and binary phase. During outbursts, the aspect ratio doubles on the star-facing side and increases by 50% on the opposite side in inclined cases. The snow line shift would impact all the species considered in the outburst case. Disk heating in binaries depends on stellar properties, orbital phase, and disk local and global characteristics. This results in temperature asymmetries, especially during secondary star outbursts, leading to variations in aspect ratio and snow lines that can affect chemistry and planet formation.

To achieve the sensitivity required to detect signals from neutral hydrogen from the Cosmic Dawn and Epoch of Reionisation it is critical to have a well-calibrated instrument which has a stable calibration over the course of the observation. Previous calibration methods do not explicitly use the time information available and make assumptions on the impedance matching of the reference sources. Here we present a new calibration method based on noise wave parameters which fits a calibration solution over time and frequency to the data, interpolating the solutions to the times at which the antenna is being measured. To test this method we simulate a dataset using measurements of the REACH receiver, modelling a low noise amplifier which is drifting over time. Fitting a polynomial surface in frequency and time to the simulated data demonstrates that we can remove the drift in the calibrated solution over time but leaves a chromatic residual. We further show that we can remove assumptions on the reflection coefficients of the reference noise source and the cold load, reducing degeneracies in the parameter fits. Applying this new calibration equation and surface fitting method to the simulated data removes the chromatic residual in the calibrated spectrum and recovers the parameters to within 0.06% of the truth and a 97% reduction in the RMSE of the spectrum of the validation source compared with previous calibration methods. For two parameters we report up to six times smaller fit error after the degeneracies are removed from the time-based calibration.

The surface array of the IceCube Neutrino Observatory, IceTop, measures cosmic rays in the PeV-EeV primary energy range. Stations comprising radio antennas and scintillation detectors will be added to enhance the existing surface detectors. A prototype station, consisting of eight scintillation detectors and three radio antennas, has been in operation in with the instrumentation in final design since the beginning of 2023. Radio signals from air showers are measured by antennas that are read-out when the trigger condition from the scintillation detectors is met. This contribution reports on air-shower coincidence measurements of these radio antennas with IceTop. Geometric shower parameters reconstructed from the radio antennas are compared with those from IceTop to determine the angular resolution. We also present details on the two new stations that were tested, deployed and commissioned with their respective data acquisition systems during the latest field season at the South Pole.

The CUbesat Solar Polarimeter (CUSP) project is a CubeSat mission planned for a launch in low-Earth orbit and aimed to measure the linear polarization of solar flares in the hard X-ray band by means of a Compton scattering polarimeter. CUSP will allow us to study the magnetic reconnection and particle acceleration in the flaring magnetic structures of our star. CUSP is a project in the framework of the Alcor Program of the Italian Space Agency aimed at developing new CubeSat missions. It is undergoing a 12-month Phase B that started in December 2024. The Compton polarimeter on board CUSP is composed of two acquisition chains based on plastic scintillators read out by Multi-Anode PhotoMultiplier Tubes for the scatterer part and GAGG crystals coupled to Avalanche PhotoDiodes for the absorbers. An event coincident between the two readout schemes will lead to a measurement of the incoming X-ray's azimuthal scattering angle, linked to the polarization of the solar flare in a statistical manner. The current status of the CUSP mission design, mission analysis, and payload scientific performance will be reported. The latter will be discussed based on preliminary laboratory results obtained in parallel with Geant4 simulations.

While astronomical twilight closes the observing window for optical astronomers, the infrared sky remains dark even through sunrise, allowing IR astronomers to observe through twilight. The Slicer Combined with an Array of Lenslets for Exoplanet Spectroscopy (SCALES) instrument is a 2-5 micron coronagraphic integral field spectrograph scheduled to arrive at Keck in early 2026. SCALES has the potential to execute exciting science and support the astronomical community and upcoming NASA missions through a dedicated cadenced twilight observing program. We estimate that the current twilight observing program on Keck conducts 18+-1 hours per year of science observations; a facilitized twilight observing program that is prioritized by the observatory could yield 151+-2 hours of science time per year. This work presents the scientific motivation and high-level feasibility of two primary SCALES twilight science cases, monitoring of Solar System objects and a high-contrast imaging search for exoplanets around bright nearby stars, taking lessons from the existing NIRC2 and OSIRIS Twilight Zone program and considering increases in program scope. We also consider technical and operational challenges to overcome before the SCALES instrument begins its twilight observing program.

The beam-crossing is a novel technique aimed at reducing residual tropospheric Doppler noise for microwave radiometer calibrations. In this work, we report the findings of the first test of this technique using ESA's Tropospheric Delay Calibration System (TDCS) at the complex in Malargue. The data consists in 14 tracking passes of the BepiColombo spacecraft collected between October 2023 and March 2024 during two separate test campaigns. We analyzed the performance of the beam-crossing technique and compared it with the nominal radiometer pointing through the analysis of the Doppler residuals extracted from the orbit determination process. Results show that the beam-crossing performed similarly to the standard pointing, with modest noise reductions and improved stability only at time scales between 100 s and 300 s. Key factors affecting the results include the antenna elevation and the boundary layer height, indicating the need to revisit initial test assumptions, which comprised a fixed boundary layer. Furthermore, comparing the beam-crossing test results with those obtained during the first two BepiColombo superior solar conjunction experiments highlights a potential application of this technique during periods of solar conjunction. However, technical challenges, adverse weather, and limited Ka-band transponder use, reduced the number of analyzed tracking passes. Future studies should therefore expand the dataset to consolidate the results. Furthermore new theoretical studies and test campaigns should elaborate on the selection process for the optimal crossing height.

In this work, we explore a possible application of a machine learning classifier for candidate events in a template-based search for gravitational-wave (GW) signals from various compact system sources. We analyze data from the O3a and O3b data acquisition campaign, during which the sensitivity of ground-based detectors is limited by real non-Gaussian noise transient. The state-of-the-art searches for such signals tipically rely on the signal-to-noise ratio (SNR) and a chi-square test to assess the consistency of the signal with an inspiral template. In addition, a combination of these and other statistical properties are used to build a 're-weighted SNR' statistics. We evaluate a Random Forest classifiers on a set of double-coincidence events identified using the MBTA pipeline. The new classifier achieves a modest but consistent increase in event detection at low false positive rates relative to the standard search. Using the output statistics from the Random Forest classifier, we compute the probability of astrophysical origin for each event, denoted as $p_\mathrm{astro}$. This is then evaluated for the events listed in existing catalogs, with results consistent with those from the standard search. Finally, we search for new possible candidates using this new statistics, with $p_\mathrm{astro} > 0.5$, obtaining a new subthreshold candidate (IFAR =0.05) event at $gps: 1240423628$ .

The field of gravitational wave (GW) detection is progressing rapidly, with several next-generation observatories on the horizon, including LISA. GW data is challenging to analyze due to highly variable signals shaped by source properties and the presence of complex noise. These factors emphasize the need for robust, advanced analysis tools. In this context, we have initiated the development of a low-latency GW detection pipeline based on quantum neural networks (QNNs). Previously, we demonstrated that QNNs can recognize GWs simulated using post-Newtonian approximations in the Newtonian limit. We then extended this work using data from the LISA Consortium, training QNNs to distinguish between noisy GW signals and pure noise. Currently, we are evaluating performance on the Sangria LISA Data Challenge dataset and comparing it against classical methods. Our results show that QNNs can reliably distinguish GW signals embedded in noise, achieving classification accuracies above 98\%. Notably, our QNN identified 5 out of 6 mergers in the Sangria blind dataset. The remaining merger, characterized by the lowest amplitude, highlights an area for future improvement in model sensitivity. This can potentially be addressed using additional mock training datasets, which we are preparing, and by testing different QNN architectures and ansatzes.

Space radiation is one of the major obstacles to space exploration. If not mitigated, radiation can interact both with biological and electronic systems, inducing damage and posing significant risk to space missions. Countermeasures can only be studied effectively with ground-based accelerators that act as a proxy for space radiation, typically with a harsher radiation field that worsen the effects of space radiation. Following an in-silico design and optimization process we have developed a galactic cosmic ray (GCR) simulator using a hybrid active-passive methodology. In this approach, the primary beam energy is actively switched and the beam interacts with specifically designed passive modulators. In this paper, we present the implementation of such a GCR simulator and its experimental microdosimetric characterization. Measuring the GCR field is of paramount importance, both before providing it to the user as a validated radiation field and for achieving the best possible radiation description. The issue is addressed in this paper by using a tissue-equivalent proportional counter to measure radiation quality and by comparing experimental measurements with Monte Carlo simulations. In conclusion, we will demonstrate the GCR simulator's capability to reproduce a GCR field.

We report optical evidence of cesium (Cs) evaporation from a bialkali (SbKCs) photo- cathode during controlled heating of a photomultiplier tube (PMT). A DFB laser scanned across the 852.113 nm Cs D2 line reveals absorption features only above 60 degrees Celsius, indicating thermal desorption. The absorption correlates with temperature and offers a non-invasive method to monitor photocathode degradation in sealed detectors.

Gravitational-wave astronomy has developed enormously over the last decade with the first detections and continuous development across broad frequency bands. However, the decihertz range has largely been left out of this development. Gravitational waves in this band are emitted by some of the most enigmatic sources, including intermediate-mass binary black hole mergers, early inspiraling compact binaries$\unicode{x2014}$whose late evolution and merger are seen by Earth-based detectors$\unicode{x2014}$, and possibly primordial gravitational waves. To tap this exciting band, we propose the construction of a detector based on pulsar timing principles, the Artificial Precision Timing Array (APTA). We envision APTA as a solar system array of artificial “pulsars”$\unicode{x2014}$precision-time-reference-carrying satellites that emit periodic electromagnetic signals towards Earth or other satellite constellation centrum. In this fundamental study, we estimate the clock precision needed for APTA to be able to detect gravitational waves. Our results suggest that 6 satellites and a clock relative uncertainty of $10^{-18}$ at 1 s of averaging, which is currently attainable with atomic clocks, would be sufficient for APTA to reach pristine sensitivity in the decihertz band and be sensitive to $10^3\unicode{x2013}10^4$ $\mathrm{M}_\odot$ black hole mergers and the early inspiral of heavy LIGO-Virgo-KAGRA sources. Future time reference, oscillator, and clock technologies realistically expected in the next decade(s) would enable the detection of an increasingly diverse set of sources and allow APTA to reach a better sensitivity than other detector concepts proposed for the decihertz band. This work opens up a new area of research into designing and constructing gravitational-wave detectors relying on principles used successfully in pulsar timing.

Separating a stochastic gravitational wave background (SGWB) from noise is a challenging statistical task. One approach to establishing a detection criterion for the SGWB is using Bayesian evidence. If the evidence ratio (Bayes factor) between models with and without the signal exceeds a certain threshold, the signal is considered detected. We present a formalism to compute the averaged Bayes factor, incorporating instrumental-noise and SGWB uncertainties. As an example, we consider the case of power-law-shaped SGWB in LISA and generate the corresponding Bayesian sensitivity curve. Unlike existing methods in the literature, which typically neglect uncertainties in both the signal and noise, our approach provides a reliable and realistic alternative. This flexible framework opens avenues for more robust stochastic gravitational wave background detection across gravitational-wave experiments.

Imminent impactors may be detected only a few hours before their impact with Earth, providing a brief opportunity to characterize them before impact. We describe the characterization of imminent impactor 2024 RW$_1$, which was discovered by the Catalina Sky Survey on 2024 September 4 at 05:43 UTC, before it entered the atmosphere near the northern Philippines at 16:39 UTC. We observed 2024 RW$_1$ with the Astrophysical Research Consortium Telescope Imaging Camera on the Apache Point Astrophysical Research Consortium's 3.5-m telescope on 2024 September 4 10:16 UTC. We obtained g, r, i, and z photometry of 2024 RW$_1$, yielding color indices of g-r = 0.47$\pm$0.04, r-i = 0.13$\pm$0.04, i-z = -0.11$\pm$0.07, and g-i = 0.60$\pm$0.04, corresponding to a spectral slope of 0.67$\pm$0.40 $\%$/100 nm. The closest match to an asteroid spectral type is with B-type asteroids from the C-complex. We detect variations in the time series photometry of the asteroid with an amplitude of $\sim$0.75, and a double-peaked rotation period of $\sim$1900 s. Assuming a visible albedo of 0.07$\pm$0.03, a density of $\sim$1500 kg/m$^3$, and a calculated absolute magnitude of 30.92$\pm$0.05, we estimate that the asteroid has a diameter of 3.3$\pm$0.7 m and a total mass of $\sim$28,000 kg. Comparing our astrometric orbital solutions to NEOMOD3, the most likely source of 2024 RW$_1$ is the 3:1 main belt mean motion resonance (77\% probability) followed by the $\nu_6$ resonance (13\% probability), consistent with its organic B-type nature.

Counting the number and brightness of ionizing radiation sources out to a redshift of z   1.2 will revolutionize our understanding of how the ionizing background is created and sustained by the embedded growth of meta-galactic structures. The sheer number of sparsely separated targets required to efficiently construct redshift binned luminosity functions is industrial in scale, driving the need for low spectral resolution multi-object spectroscopy (MOS) with a short wavelength cut-off   1000 Å, a sensitivity in the far-UV to better than 30 abmag, and an instantaneous field-of-view   (2')$^2$. A MOS on Habitable Worlds Observatory is the only instrument that could conceivably carry out such an ambitious observing program. This program will quantify how much of the ionizing radiation produced by galaxies is attenuated by intervening neutral H, He and dust, and how much escapes to maintain the universe in a mostly ionized state.

As Radial velocity (RV) spectrographs reach unprecedented precision and stability below 1 m/s, the challenge of granulation in the context of exoplanet detection has intensified. Despite promising advancements in post-processing tools, granulation remains a significant concern for the EPRV community. We present a pilot study to detect and characterise granulation using the High-Accuracy Radial-velocity Planet Searcher for the Northern hemisphere (HARPS-N) spectrograph. We observed HD166620, a K2 star in the Maunder Minimum phase, intensely for two successive nights, expecting granulation to be the dominant nightly noise source in the absence of strong magnetic activity. Following the correction for a newly identified instrumental signature arising from illumination variations across the CCD, we detected the granulation signal using structure functions and a one-component Gaussian Process (GP) model. The granulation signal exhibits a characteristic timescale of 43.65$\pm$15.8 minutes, within one $\sigma$, and a standard deviation of 22.9$\pm$0.77 cm/s, with in three $\sigma$ of the predicted value. By examining spectra and RVs as a function of line formation temperature , we investigated the sensitivity of granulation-induced RV variations across different photospheric layers. We extracted RVs from various photospheric depths using both the line-by-line (LBL) and cross-correlation function (CCF) methods to mitigate any extraction method biases. Our findings indicate that granulation variability is detectable in both temperature bins, with the cooler bins, corresponding to the shallower layers of the photosphere, aligning more closely with predicted values.

The Cherenkov Telescope Array Observatory (CTAO) will greatly improve upon sensitivities in the field of very-high-energy gamma-ray astrophysics. The CTAO northern site (CTAO-North, La Palma, Spain) currently hosts LST-1 with the remaining three large-sized telescopes (LSTs) expected in mid-2026 and one medium-sized telescope (MST) expected in mid-2027. The CTAO southern site (CTAO-South, Paranal, Chile) expects the delivery of five small-sized telescopes (SSTs) and two MSTs in early 2026 with on-site construction beginning in mid-2026. The dual-mirrored Schwarzschild-Couder Telescope (SCT) is a candidate MST for CTAO-South and is capable of observations in the energy range of 100 GeV to 10 TeV, the core of CTAO's 20 GeV to 300 TeV energy range. Inaugurated in January 2019, the prototype SCT (pSCT) located at the Fred Lawrence Whipple Observatory in southern Arizona observed gamma-ray emission from the Crab Nebula at a significance of 8.6 sigma in 2020. The pSCT utilizes a novel dual-mirror optics design and a densely packed focal plane of silicon photomultipliers (SiPMs). An upgrade of the pSCT camera is underway to fully instrument the camera with 11,328 pixels and an 8-degree diameter FoV. In addition, upgraded electronics will lower the front-end electronics noise, allowing for a lower trigger threshold and improved event reconstruction and background rejection. This work will present the status of the upgrade of the pSCT and discuss the future of the SCT.

Strong gravitational lensing can produce copies of gravitational-wave signals from the same source with the same waveform morphologies but different amplitudes and arrival times. Some of these strongly-lensed gravitational-wave signals can be demagnified and become subthreshold. We present TESLA-X, an enhanced approach to the original GstLAL-based TargetEd Subthreshold Lensing seArch (TESLA) method, for improving the detection efficiency of these potential subthreshold lensed signals. TESLA-X utilizes lensed injections to generate a targeted population model and a targeted template bank. We compare the performance of a full template bank search, TESLA, and TESLA-X methods via a simulation campaign, and demonstrate the performance of TESLA-X in recovering lensed injections, particularly targeting a mock event. Our results show that the TESLA-X method achieves a maximum of $\sim 10\%$ higher search sensitivity compared to the TESLA method within the subthreshold regime, presenting a step towards detecting the first lensed gravitational wave. TESLA-X will be employed for the LIGO-Virgo-KAGRA's collaboration-wide analysis to search for lensing signatures in the fourth observing run.

Technical and environmental noise in ground-based laser interferometers designed for gravitational-wave observations like Advanced LIGO, Advanced Virgo and KAGRA, can manifest as narrow (<1Hz) or broadband ($10'$s or even $100'$s of Hz) spectral lines and features in the instruments' strain amplitude spectral density. When the sources of this noise cannot be identified or removed, in cases where there are witness sensors sensitive to this noise source, denoising of the gravitational-wave strain channel can be performed in software, enabling recovery of instrument sensitivity over affected frequency bands. This noise hunting and removal process can be particularly challenging due to the wealth of auxiliary channels monitoring the interferometry and the environment and the non-linear couplings that may be present. In this work, we present a comprehensive analysis approach and corresponding cyberinfrastructure to promptly identify and remove noise in software using machine learning techniques. The approach builds on earlier work (referred to as DeepClean) in using machine learning methods for linear and non-linear regression of noise. We demonstrate how this procedure can be operated and optimized in a tandem fashion close to online data taking; it starts off with a coherence monitoring analysis that first singles out and prioritizes witness channels that can then be used by DeepClean. The resulting denoised strain by DeepClean reflects a 1.4\% improvement in the binary neutron star range, which can translate into a 4.3\% increase in the sensitive volume. This cyber infrastructure we refer to as Coherence DeepClean, or CDC, is a significant step toward autonomous operations of noise subtraction for ground-based interferometers.

We propose an O(100)m Atom Interferometer (AI) experiment - AICE - to be installed against a wall of the PX46 access shaft to the LHC. This experiment would probe unexplored ranges of the possible couplings of bosonic ultralight dark matter (ULDM) to atomic constituents and undertake a pioneering search for gravitational waves (GWs) at frequencies intermediate between those to which existing and planned experiments are sensitive, among other fundamental physics studies. A conceptual feasibility study showed that this AI experiment could be isolated from the LHC by installing a shielding wall in the TX46 gallery, and surveyed issues related to the proximity of the LHC machine, finding no technical obstacles. A detailed technical implementation study has shown that the preparatory civil-engineering work, installation of bespoke radiation shielding, deployment of access-control systems and safety alarms, and installation of an elevator platform could be carried out during LS3, allowing installation and operation of the AICE detector to proceed during Run 4 without impacting HL-LHC operation. These studies have established that PX46 is a uniquely promising location for an AI experiment. We foresee that, if the CERN management encourages this Letter of Intent, a significant fraction of the Terrestrial Very Long Baseline Atom Interferometer (TVLBAI) Proto-Collaboration may wish to contribute to AICE.