Priority

Previous Sirocco radiative-transfer models of gas-cocooned AGN predicted lower-column counterparts to little red dots (LRDs): compact, X-ray-weak sources with bluer continuum slopes and Balmer jumps rather than Balmer breaks. The recently identified population of little blue dots (LBDs) closely resembles this predicted phase. Here we explore these lower-column-density cocoons in which nebular recombination emission remains visible while strong Balmer-continuum absorption is avoided. We find that a sequence of increasing column density connects more classical AGN spectra, Balmer-jump LBD-like spectra at $N_{\rm H}\!\sim\!{\rm few}\times10^{24} \mathrm{cm^{-2}}$, and Balmer-break LRD-like spectra at higher columns. In this sequence, electron scattering produces exponential line wings and suppresses X-ray emission before strong Balmer absorption features, characteristic of higher column densities, appear. We therefore propose that LBDs are lower-column analogues of LRDs within a common gas-cocooned AGN sequence. This interpretation predicts that Balmer-jump emission, X-ray weakness, permitted lines with exponential wings, He II $\lambda$4686 emission, smaller H$\alpha$ FWHM values and equivalent widths than in LRDs, and weak or absent absorption features are characteristic of LBDs. We compare to three example LBD spectra and identify Balmer-jump signatures in them.

This chapter presents an overview of the unification models for Active Galactic Nuclei (AGN), focusing on the physical structures, classification schemes, and evolutionary processes that characterize accreting supermassive black holes. We introduce the fundamental components of AGN, including the supermassive black hole, accretion disk, jets, outflows, broad-line and narrow-line regions, polar dust and the dusty anisotropic obscurer. The traditional orientation-based unification model is reviewed, with a focus on the role of the covering factor of the obscuring material in shaping observed properties. We introduce the radiation-regulated unification model, which accounts for the influence of radiative feedback on the nuclear environment of SMBHs. We also examine the evolutionary aspects of AGN unification, including the impact of galaxy mergers, host galaxy properties, and redshift-dependent trends. Finally, we discuss changing-look AGN, which challenge conventional unification frameworks by exhibiting dramatic spectral variability.

The globular cluster $\omega$ Cen has been suggested to have originated in the Gaia Sausage-Enceladus (GSE) merger event, possibly as its nuclear star cluster. However, the present-day orbits of $\omega$ Cen and the GSE debris are very different. We investigate the scenario in which $\omega$ Cen originated in the GSE and migrated to its current position due to perturbations from the Galactic bar. The [$\alpha$/M] distributions of stars located between the GSE debris and $\omega$ Cen in $(L_z,E)$ space tentatively support this scenario, but are not conclusive. We run simulations of the GSE debris and $\omega$ Cen in a realistic Milky Way potential with a decelerating bar at various present-day pattern speeds. We find that $\omega$ Cen can indeed be traced back to the phase space region occupied by the GSE debris. However, this likely requires a pattern speed of $\Omega_\mathrm{b}\lesssim26$ km s$^{-1}$ kpc$^{-1}$, which is much lower than most recent estimates. We conclude that a GSE origin for $\omega$ Cen is dynamically and chemically plausible, but this would require a re-evaluation of the current consensus on the bar's pattern speed.

In the last decades a number of extremely faint and compact Galactic satellites (Ultra Faint Compact Satellites; UFCS) have been discovered by large panoramic surveys. Their nature is uncertain due to their location in the overlapping dwarf galaxy-star cluster region of the $M_V-R_h$ plane and their faintness and distance. Here we show how the deep HST photometry from the Missing Globular Clusters Survey (MGCS), combined with spectroscopic metallicities, provides new insight into the nature of these satellites through accurate distance and age estimates. We consider the case of Koposov 2, currently the most metal-poor bound star cluster known in the entire Milky Way or an extreme case of Ultra Faint Dwarf galaxy. By performing a spectroscopically-informed bayesian isochrone fit on the MGCS data we find $(m-M)_0=16.85\pm 0.06$ ($D=23.4\pm 0.6$ kpc) and age=$13.7^{+0.9}_{-1.3}$ Gyr, showing that, contrary to previous age estimates, Koposov 2 is as old as the oldest Galactic globular clusters. The luminosity function, corrected for incompleteness, is well reproduced by a model with the same age and metallicity and a slope of the mass function $x=-0.35$, suggesting a significant depletion of faint stars. We model the surface stellar density field, deriving new robust estimates of the half-light radius ($R_h=0.39^{+0.06}_{-0.04}$ arcsec, corresponding to $R_h=2.7^{+0.4}_{-0.3}$ pc), of the absolute integrated magnitude ($M_V=-0.95\pm0.22$) and of the stellar mass ($M_{\star}=371.8\pm41.6M_{\odot}$), showing that Koposov 2 is much more compact than dwarf galaxies of similar stellar mass. The new evidence significantly support the hypothesis that Koposov 2 is a star cluster that may have lost a large fraction of its original mass. Finally we show that most UFCS lie in the same locus of the $M_V - R_h$ plane as Galactic open clusters, hinting to a possible additional channel for their formation.

Magnetic fields are a fundamental ingredient of the Universe, influencing the formation and evolution of cosmic structures. While magnetic fields in local galaxy clusters have been studied, their origin, amplification, and strength at high redshift are poorly understood. Proto-clusters represent the early stages of galaxy cluster formation, ideal for investigating the early magnetisation of the intra-cluster medium (ICM). We present a study of CARLA J1510+5958 proto-cluster at z = 1.72, observed with the JVLA in the L-band (1-2 GHz). We aim to investigate the magnetic field strength and structure in the proto-ICM and the role of AGN in magnetising the environment during early cluster formation. We analyse Faraday rotation on the polarised emission from the central radio-loud AGN using the Rotation Measure (RM) synthesis and QU fitting technique. We further interpret the observations with 3D simulations of gas density and turbulent magnetic fields, varying AGN orientation and path length. The two AGN lobes show different polarisation properties. The Western lobe exhibits a uniform RM (average $-115 \pm 32\text{ rad m}^{-2}$, dispersion $36 \pm 11\text{ rad m}^{-2}$), indicating a locally ordered magnetic field likely compressed by the lobe, while the Eastern lobe is depolarised. Although the asymmetry indicates a turbulent, magnetised medium, simulations rule out a purely isotropic random field for the Western lobe RM distribution. The QU fitting further suggests an internal Faraday component, interpreted as magnetised relativistic plasma from the lobe mixed with the surrounding gas, indicating possible magnetisation of the ambient medium by the AGN. From this asymmetry, we constrain the average physical magnetic field in the proto-ICM to a lower limit of 0.4 $\mu$G. These results confirm a magnetised proto-ICM at z = 1.72, proving early field amplification during cluster assembly.

The Coma cluster (Abell 1656; $z=0.023$) is a nearby rich galaxy cluster and a key laboratory for studying cluster assembly in the Cosmic Web. We characterize its projected dark matter distribution and connection to galaxies, the intracluster medium, and reported intracluster filaments (ICFs) with wide-field ($\sim$12-deg$^2$) Subaru/Hyper Suprime-Cam weak-lensing (WL) analysis. We reconstruct the two-dimensional mass distribution, fit Navarro-Frenk-White (NFW) models, derive an aperture mass densitometry profile, and compare the WL signal with optical spectroscopy, eROSITA X-ray observations, radio data, and gas fraction diagnostics. A single-halo NFW fit yields $M_{200\mathrm{c}}=8.2\pm0.7\times10^{14}~M_{\odot}$. The aperture mass profile agrees with the best-fit NFW model and the X-ray hydrostatic mass at $R\gtrsim20'$ ($\sim$560 kpc), suggesting little merger-induced bias in the global WL mass, while the inner region shows substantial hydrostatic bias. A two-halo NFW fit centered on NGC 4874 and NGC 4839 gives masses of $7.8\pm0.6$ and $0.9\pm0.2\times10^{14}~M_{\odot}$, implying a $\sim$1:8 minor merger. The gas mass fraction suggests that the system is returning from first apocenter. We find a positive spatial correlation between the WL signal and X-ray surface brightness, strongest along the ICF directions ($110^{\circ}$ and $340^{\circ}$), where shear-selected subhalos are predominantly detected. The Coma $r$-band mass-to-light ratio is radially constant with $\langle M/L_r\rangle\simeq250\pm66~M_{\odot}/L_{\odot}$ within $R_{200\mathrm{c}}$, whereas the northern and western ICFs show higher values of $\sim1000~M_{\odot}/L_{\odot}$, suggesting stronger dark matter dominance. These results show that joint WL and multiwavelength analyses can effectively probe cluster assembly and the dark matter content of ICFs.

We present a Bayesian analysis of star clusters in the dwarf starburst galaxy NGC 1569 based on high-resolution Hubble Space Telescope imaging combined with integral-field spectroscopy from the Keck Cosmic Web Imager, obtained as part of the DUVET survey. For each cluster identified, we infer posterior probability distributions for mass and age using a forward modelling method that properly accounts for uncertainties due to stochastic sampling of the IMF. We investigate how the inferred properties depend on photometric coverage by repeating the analysis with different filter combinations, including mock extensions to the ultraviolet and near-infrared that emulate the addition of HST UV bands and James Webb Space Telescope imaging. We find that, while inclusion of these wavelength regimes further breaks age and mass degeneracies, the currently available data yields reasonably strong constraints on cluster parameters. We compare inferred cluster properties to the conditions of the local interstellar medium, and find evidence for multiple interesting correlations. The truncation mass of the cluster mass function varies with galactocentric distance, particularly moving off the disk, consistent with a dependence on the density of the interstellar medium. Cluster mass positively correlates with metallicity, suggesting that massive clusters preferentially form in pre-enriched gas, and the ionisation state of the gas, reflecting the increased prevalence of high-mass stars in high-mass clusters. These results demonstrate the power of Bayesian, initial mass function-aware modelling for resolving cluster populations in nearby starburst dwarfs and provide new insight into how cluster formation and feedback respond to local galactic conditions.

We present a spatially resolved chemical abundance analysis of 2291 star-forming galaxies at $z \leq 0.1$, spanning nearly four orders of magnitude in stellar mass ($8 \le \rm log (M_{\star}/M_{\odot}) \le 11.5$), by exploiting the multi-fibre spectra from the Dark Energy Spectroscopic Instrument (DESI) as a sparse integral field spectrograph. In the inner regions ($<2R_e$), the radial gas-phase metallicity profiles show an outward-declining trend for massive galaxies, with the steepest gradient ($\nabla_{log(O/H)}$) $\sim-0.08$ dex/R$_{e}$, whereas low-mass dwarf galaxies exhibit nearly flat profiles ($\nabla_{log(O/H)}\sim-0.02$ dex/R$_{e}$). The large galactocentric radii ($\sim$5 R$_{e}$) probed in this study, reveal flat metallicity profiles near the disk-halo interface. Strikingly, these flat metallicity values are consistent across a wide stellar mass range, likely reflecting the influence of low SFR and metal poor inflows in the outer regions. The metallicity gradient - stellar mass relation exhibits a turnover at $\log(M_\star/M_\odot) \sim 10.5$, beyond which gradients become shallower, possibly driven by the chemical equilibrium in the inner disk of massive galaxies and/or dilution from cosmic gas accretion. At fixed stellar mass, a strong size dependence is observed, where compact galaxies show flatter gradients and higher central enrichment than their extended counterparts. The abundance gradients are further linked with the stellar age distribution within the galactic disk, where galaxies with younger outskirts show steeper gradients than the ones with older outskirts, consistent with ongoing inside-out disc growth sustaining centrally concentrated chemical enrichment. These results underscore the interplay of star formation efficiency, stellar feedback, and metal-poor gas accretion in governing the radial chemical structure in galaxies.

We derive both the mid-plane and off-plane rotation curves, $v_c(R,z)$, and the vertical acceleration, $a_z(R,z)$, of the Milky Way (MW) using \textit{Gaia} DR3 data over the ranges of vertical heights $z \in (-2,2)\,$ kpc and galactocentric distances $R \in (8.5,14)$ kpc where the velocity components are determined with high precision, i.e., with an error $< 5\%$. In contrast, the vertical acceleration $a_z(R,z)$ is dominated by model-dependent systematics, with uncertainties of up to $\sim 20\%$. This level of accuracy allows us to place stringent constraints on the geometry of the MW's dark matter (DM) distribution, as the vertical gradients of the gravitational potential attain their maximum within this range of radial and vertical distances corresponding to the characteristic scales of the disk. We find that models including the observed stellar components together with a spherical DM halo fail to reproduce both the pronounced variation of $v_c(R,z)$ with height and the observed behavior of $a_z(R,z)$. In particular, spherical halos with a scale radius of $r_s \sim 15$ kpc contribute negligibly to the off-plane rotation curve and vertical acceleration in the inner disk, leaving these features primarily determined by the stellar mass distribution. Conversely, models in which DM is confined to a flattened, disk-like configuration predict substantial contributions to both $v_c(R,z)$ and $a_z(R,z)$, resulting in a markedly better agreement with the data. We conclude that disk-like DM distributions are strongly favored over spherical halo models. Forthcoming Gaia data releases will enable even more stringent tests of the geometry and distribution of the MW's DM component.

Modern surveys such as DESI and \textit{Euclid}, which collect data for hundreds of millions of galaxies to map the large-scale structure (LSS) of the Universe, hold the key to determining the cosmological parameters and testing new physics. This ambition, however, is limited by uncertainties in the galaxy-halo connection: the link between observed galaxies and the underlying, unobservable matter field, by accounting for effects such as galaxy bias and assembly bias (AB). These are particularly poorly-understood for modified gravity (MG) models, which are popular alternatives to the cosmological constant to explain accelerated expansion. We approach this problem using mock emission line galaxy (ELG) and luminous red galaxy (LRG) catalogues in $f(R)$ gravity matching the specifications of ongoing Stage-IV galaxy surveys, generated from state-of-the-art MG hydrodynamical simulations. While the interplay between MG – especially the chameleon screening mechanism – and galaxy formation leaves complicated imprints in the galaxy-halo connection, a simple physical picture emerges in which halo and galaxy formation are enhanced for progressively more massive haloes over time. We confirm that the basic galaxy-halo connection model, the halo occupation distribution (HOD), in which galaxy occupation is determined solely by halo mass, underestimates galaxy clustering strength in $\Lambda$CDM by $10$–$20\%$ at $z\lesssim1$ when neglecting AB, and demonstrate that MG introduces further complexity. Extending this model with a suitably-chosen environment density as a secondary HOD variable reduces the AB effect in all models to $2$–$3\%$ for $z\lesssim0.5$. This provides a well-motivated starting point for further works on minimising the impact of AB when testing non-standard cosmological models with LSS.

The inner Milky Way hosts overlapping stellar populations (bar–bulge, inner thin and thick disks, and halo), complicating population assignments along the line of sight. A joint chemical–dynamical approach is required to isolate a clean field bulge sample. Metal-poor bulge stars are valuable as they likely trace the earliest phases of chemical enrichment in the inner Galaxy. We aim to characterize the alpha-element (Si, Mg) and selected Fe-peak abundances of a dynamically defined sample of bulge field stars, and to contrast these trends with the inner-halo tail in the metal-poor regime. We analyze metal-poor candidates from the Pristine Inner Galaxy Survey observed by the bulge Cluster APOGEE Survey. Using APOGEE/ASPCAP abundances (S/N >= 50), we integrate full 6D orbits in a barred Milky Way potential. Bulge membership is defined via orbital confinement using an apocenter cut and the high-density locus in the (E_J, L_z) plane, with uncertainties estimated using Monte Carlo simulations. We identify 98 stars as bulge members. The metallicity distribution spans -2.5 <= [Fe/H] <= -0.4, with a median [Fe/H] = -1.71, offset to higher metallicity than the inner halo (median [Fe/H] = -1.96). The sample follows a high-alpha sequence with slopes d[Si/Fe]/d[Fe/H] = -0.020(+0.012,-0.051) and d[Mg/Fe]/d[Fe/H] = -0.097(+0.062,-0.133). Fe-peak tracers show [Ni/Fe]   0, while [Mn/Fe] declines with [Fe/H] with a mild upturn at higher metallicity. We detect no significant [Si/Fe] gradients with R_apo or Z_max (+0.010(+0.018,0.000) and +0.006(+0.022,-0.011) dex kpc^-1, respectively). Results are insensitive to Omega_p, yielding indistinguishable memberships. The chemo-orbital evidence favors an in situ origin within the inner Galaxy for the metal-poor bulge field, enriched at early times and rearranged by secular bar evolution, with a minor contribution from halo stars at the most metal-poor end.

We constructed the 2D model of the Galaxy which initially includes an analytical bar, bulge, disk and halo. The model disk forms the outer elliptical resonance rings R1 and R2 located near the outer Lindblad resonance of the bar (OLR), as well as the inner resonance ring r located near the corotation radius (CR). As the density of stars in the elliptical rings increased, we introduced additional gravitational perturbations created by the rings. The radial component of gravitational perturbations from the elliptical rings, F_R, at a point with the Galactocentric coordinates (R, theta) was represented as a combination of three polynomials in powers R/Re or Re/R, where Re is the distance to the midline (middle) of the ring at a given angle theta. The azimuthal component of the disturbances, F_T, was calculated using the force F_R. The difference between the values of the force F_R (F_T) calculated using the numerical differentiation of the potential and using the analytical representation does not exceed 5.7% (1.3%) of the maximum value of the force F_R generated by the elliptical rings. In general, the gravity of the elliptical rings has little effect on the process of adjustment of epicyclic motions near the OLR of the bar.

We perform a cosmological zoom-in simulation of a massive galaxy ($M_s\sim10^{10}\rm M_\odot$ at $z\sim10$) using the GIZMO code. By employing $\leq 3\rm pc$ resolution and a $3.4\rm Myr$ supernova feedback delay, we capture the feedback-free starbursts (FFB) in clusters. The simulation reproduces FFB model predictions and super-bright galaxies observed by JWST. At $z\sim10$, cold streams feed a compact galaxy ($R_{\rm e}\sim1\rm kpc$), with stellar and surface densities ($>10^5\rm cm^{-3}$, $>10^5\rm M_\odot pc^{-2}$) exceeding FFB thresholds. The global star-formation efficiency (SFE) is $\varepsilon_s\sim0.2\text{--}0.3$, associated with a fluctuating star-formation history. We identified over $10^5$ star clusters ($M_{\star}>10^{4.5}\rm M_\odot$) with a nearly scale-free mass distribution (${\rm d}N/{{\rm d}\log M}\propto M^{-1.06}$). Approximately 90\% of star formation occurs in clusters, which at a given time constitute $30\text{--}40\%$ of the total stellar mass. The star formation in most of the clusters of masses $<10^7\rm M_\odot$, occurs in bursts of $<3\rm Myr$ and a local SFE $\sim0.5\pm 0.2$. Cluster metallicities ($-2.01<\log (Z/Z_\odot)<-0.45$) indicate rapid baryon recycling. Feedback-driven outflows exhibit typical temperature of $10^7\rm K$ and typical velocities of $\sim 2000\rm km\ s^{-1}$. In the highly dynamic central $1\rm kpc$, clusters undergo rapid orbital decay and merge to assemble the oblate nuclear stellar cluster. Cluster shapes range from oblate to prolate, with a triaxial median. These clusters are consistent with JWST observations, and a fraction of them may survive to yield the globular clusters (GCs) at low redshifts.

Carbon-enhanced metal-poor (CEMP; with $\rm{[Fe/H]} \le -2.0$ and $\rm{[C/Fe]} \ge 0.7$) stars preserve information about early chemical enrichment, low-mass star formation, and the hierarchical assembly of galaxies. In this study, we have compiled an extensive literature sample of 1032 stellar carbon abundances spanning the metal-poor Milky Way halo (437 stars), 21 ultra-faint dwarf galaxies (UFDs; 102 stars), seven classical dwarf spheroidal galaxies (254 stars), three accreted dwarf galaxies (90 stars), the Small Accreted Stellar Systems (SASS; 77 stars), and eleven stellar streams (72 stars). We establish the fractions of CEMP stars for each of these systems and categories. Generally, the low-mass UFDs possess the high fractions at low metallicities, whereas the more massive classical dwarf galaxies have relatively few CEMP stars. This behavior reveals a new low-metallicity Magnitude ($M_{\rm V}$)–CEMP Fraction relation across the dwarf satellite galaxy population. The high CEMP fractions in surviving UFDs suggest their enrichment was dominated by faint supernovae, as higher energy input would likely have quenched star production. The low CEMP fractions in classical dwarfs imply predominantly in situ formation rather than assembly from smaller systems. Using $\rm{[C/H]}$ abundances, we also probe early low-mass star formation. Eight stars lie within or near the theoretical ”forbidden zone”, indicating that dust-induced cooling, alongside fine-structure line cooling, contributed to early star formation. These rare dust-cooled stars may have formed in UFD-like systems that did not survive. Overall, the metal-poor Milky Way halo appears to have assembled from many different dwarf galaxies, with CEMP halo stars being contributed by early UFD-like systems and non-CEMP halo stars by intermediate-sized halos that later formed classical dwarfs.

This study presents a comprehensive X-ray spectroscopic analysis of 642 quasars, obtained by cross-matching the XMM-Newton Serendipitous Source Catalog (DR11) with the Sloan Digital Sky Survey (DR16) quasar catalog. After stringent quality filtering and automated spectral reduction, we derived reliable photon indices ($\Gamma$) and intrinsic 2–10 keV X-ray luminosities. Using multiwavelength data, sources were classified into 561 radio-quiet (RQ) and 81 radio-loud (RL) quasars. We estimate the bolometric luminosity and Eddington ratio ($\lambda_{\mathrm{Edd}}$) from absorption-corrected X-ray measurements and virial black hole masses. Our primary objective is to establish and characterize the fundamental relationship between the photon index and Eddington ratio for RQ population. We find that RQ quasars exhibit systematically higher Eddington ratios, peaking at $\log \lambda_{\mathrm{Edd}} \approx -1.2$, and softer spectra with $\Gamma \approx 2.0$. A statistically significant positive correlation between $\Gamma$ and $\lambda_{\mathrm{Edd}}$ is detected in RQ quasars, supporting disk–corona coupling models. To validate our results within the broader context of AGN evolution, we further examine the dependence of $\lambda_{\mathrm{Edd}}$ on redshift ($z$) and black hole mass ($M_{\mathrm{BH}}$). For RQ quasars, $\lambda_{\mathrm{Edd}}$ increases with redshift and decreases with $M_{\mathrm{BH}}$, in strong agreement with recent results \cite{aggarwal2024evidence}, highlighting the universal nature of these accretion trends. By correlating spectral slope with accretion rate, this work provides new insights into the interplay between accretion physics, jet activity, and the cosmic evolution of quasars.

Using a sample of 151 galaxies with dynamically measured black hole (BH) masses ($M_{\rm BH}$), we investigate the scaling relations between $M_{\rm BH}$ and the stellar velocity dispersion, $\sigma$, and, for a subsample of 30 core-Sérsic galaxies, between $M_{\rm BH}$ and the size of the partially depleted core, $R_{\rm b}$. Core-Sérsic galaxies, identified using high-resolution $Hubble ~ Space ~ Telescope$ imaging and spanning both the normal-core $(R_{\rm b}<0.5$ kpc) and large-core ($R_{\rm b}>0.5$ kpc) regimes, define an updated $M_{\rm BH}$$-$$R_{\rm b}$ relation of the form $M_{\rm BH} \propto R_{\rm b}^{1.16 \pm 0.10}$, with an rms scatter of $\Delta_{\rm rms} \simeq 0.28$ dex in $\log M_{\rm BH}$. We find that Sérsic and normal-core galaxies together follow a common log-linear $M_{\rm BH}$-$\sigma$ relation with a slope of $4.95 \pm 0.29$ and a scatter $\Delta_{\rm rms} \simeq 0.46$ dex. Deviations from this relation arise at the highest BH masses, where large-core galaxies, including six with direct $M_{\rm BH}$ measurements, drive a significant upturn. We find that these galaxies typically host ultramassive black holes whose masses scale more strongly with $R_{\rm b}$ than $\sigma$, and lie $\sim (1-4$) $~\times$ the intrinsic scatter (0.39 dex) above the relation defined by Sérsic and normal-core galaxies. The $M_{\rm BH}$$-$$R_{\rm b}$ relation shows $\sim 30$-$47\%$ less scatter in $\log M_{\rm BH}$ than the corresponding $M_{\rm BH}$$-$$\sigma$ relation for the same sample. We interpret the high-mass upturn in the $M_{\rm BH}$$-$$\sigma$ diagram as a consequence of successive major, dry mergers, a scenario that naturally explains the observed flattening of the $\sigma-L_V$ relation at $M_V < -23.5$ mag.

We investigate how low-ionization gas flows in typical star-forming galaxies at $z\sim3$ depend on galaxy intrinsic properties and viewing angle. For this analysis we use JWST/NIRSpec observations of rest-frame near-UV Fe II and Mg II absorption, and rest-frame optical Na D absorption. This study combines galaxies from the LyC22, EXCELS, and AURORA surveys and contains 176, 197, and 315 galaxies, respectively, with Fe II, Mg II, and Na D coverage. Based on both individual and composite spectra, we find no statistically significant correlations between outflow velocity and galaxy properties. However, galaxies with detected outflows tend towards higher stellar masses, SFR, and $\Sigma_{\rm SFR}$ than those without outflows, suggesting that the two samples are not drawn from the same parent population. Finally, we additionally find that Mg II emission is preferentially detected in galaxies with lower stellar mass and $A_V$, and higher sSFR, consistent with conditions that favor the escape of resonantly scattered line and ionizing continuum radiation. We present the first evidence in $z\sim3$ star-forming galaxies that properties of the absorption lines depend on galaxy inclination, with more face-on systems showing stronger absorption and higher outflow velocities, while inflowing gas is more frequently detected in more highly inclined galaxies. These trends are consistent with observations at $z\lesssim1$ and predictions from cosmological simulations in which galactic winds are launched perpendicular to the galactic disks, while accretion occurs primarily along the disk plane.

The JWST discoveries of supermassive black holes (BHs) at $z \gtrsim 5$ may provide key insights into their seeding origins. Using new $[18{-}72~\rm Mpc]^3$ BRAHMA cosmological simulations, we investigate how variations in heavy-seed prescriptions, coupled with a subgrid dynamical friction model, shape BH populations at $z \sim 5$ and $z \sim 0$. We consider two "lenient” seed models, in which all halos containing sufficient dense & metal-poor gas form $\sim10^4$ and $\sim10^5~M_{\odot}$ seeds, and a "strict” seed model, in which $\sim10^5 M_{\odot}$ seeds form only under additional constraints motivated by direct collapse black hole formation. By $z \sim 5$, all models produce $M_*-M_{\rm BH}$ relations broadly consistent with the observed local Universe for $M_*\gtrsim10^9~M_{\odot}$ galaxies, but only the lenient scenarios generate systems near the upper envelope of the observed local scatter. In galaxies hosting $M_{\rm BH} \sim 10^8$-$10^9~M_{\odot}$ BHs, lenient production of $\sim10^5~M_{\odot}$ seeds also produces multiple overmassive systems with $M_{\rm BH}/M_* \gtrsim 0.01$. Although their growth is dominated by seeding and mergers, these systems reach luminosities of $\sim10^{43}$-$10^{45}\mathrm{erg s^{-1}}$, comparable to those inferred for JWST-detected BHs. As a key observational signature, the lenient seed models yield merger rates of $\gtrsim100\mathrm{yr^{-1}}$ and near-unity local BH occupation fractions even in galaxies with $M_* \lesssim 10^7~M_{\odot}$. In contrast, the strict seed model produces merger rates of only $\sim1\mathrm{yr^{-1}}$ and local occupation fractions of $\lesssim10\%$ for galaxies with $M_* \lesssim 10^8~M_{\odot}$. Future gravitational-wave event rates and measurements of local BH occupation fractions will therefore provide strong constraints on the dominant pathways responsible for high-redshift BH assembly.

A $10^6\,M_\odot$ object has recently been inferred from gravitational imaging of the strong-lensing system JVAS B1938+666, exhibiting an unusually dense inner region embedded within an extended envelope, far exceeding expectations for cold dark matter (CDM) halos. Using gravothermal fluid simulations, we show that such a structure arises naturally in self-interacting dark matter (SIDM) halos evolving into a deep core-collapse phase, where a secondary dense central core forms within an extended profile. The resulting density structure closely matches the inferred properties of the lensing object. We also demonstrate that a similar profile could be reproduced in CDM in the presence of an intermediate-mass black hole, but this requires an early-forming progenitor that subsequently loses $5$ orders of magnitude in mass through tidal stripping by the lens galaxy. Whether such a scenario can be realized in realistic cosmological environments remains an open question.

The first generation of stars formed from pristine, neutral hydrogen gas. The most massive of these exploded as supernovae within a few million years of their birth, producing the first heavier elements and leaving distinct chemical signatures of their origin in the surrounding medium. However, chemical abundance studies have so far mainly relied on emission-line measurements, which are luminosity weighted and hence biased towards the most recently formed stars. Here we analyse near-infrared, medium-resolution spectroscopy from the JWST-SPURS program of three UV-bright galaxies at redshifts 7.8, 8.6, and 9.3, within the first 650 to 520 million years after the Big Bang. The chemical abundance patterns of the metal lines detected in absorption hint at extremely metal-poor gas, substantially lower than inferred from the emission lines tracing the central, star-forming regions. Further, they all exhibit super-solar [C/O] abundances, which is also imprinted in the averaged spectrum of a larger set of galaxies at similar redshifts. These results reveal the distinct chemical signatures of the first Population III supernovae explosions.

This study reports the nature of the companion star to the pulsating ULX NGC 1313 X-2, using long-slit spectroscopic data from Gemini-South observations, based on archival data from 2009. After stacking flux-calibrated spectra from ten nights of observations and fitting the spectra with stellar templates, we find a possible Balmer break in the GMOS-S spectrum below 4000 Angstroms, which is suggestive of an A-type supergiant donor. Using the inferred stellar radii, we report updated constraints on the orbital parameters of the system and on the nature of the binary. We also add some information on the accretion disc size scale by studying the X-ray and optical variability using the lag-frequency spectrum and corroborate on results from earlier studies for the gas bubble expansion rates by modelling the [O III] emission line profiles, allowing constraints on the kinetic power of the wind/jet relative to the accretion power. This study also expands on previous efforts to study the formation history of the binary using multi-wavelength observations.

Optical-to-near infrared interstellar polarisation, induced by aligned dust grains, generally follows a convex wavelength dependence, known as the Serkowski relation. However, observations in the ultraviolet (UV) and at [mid-]infrared wavelengths indicated that some of the spectra do not follow this relation. Specifically, about 25% show an excess in the degree of polarisation at mid-UV wavelengths ($\lambda^{-1} > 3\,\rm \mu m^{-1}$), referred to as the super-Serkowski polarisation. In this study, we re-examine both the Serkowski and super-Serkowski spectra based on the joint effect of paramagnetic relaxation (DG alignment) and radiative torque (RAT) alignment. We used the observational data for HD 30614, HD 204827, HD 37903 and HD 161056 to constrain our modelling. We examined two types of radiation fields: one derived from the scaled interstellar radiation field and another originating from a B-type star. For the super-Serkowski spectra of HD 30614 and HD 204827, our model demonstrates that RAT alignment enhanced by radiation produced from a B-type star below the Lyman limit ($\lambda=912 Å$) can reasonably explain the observations and that a combination with the DG alignments results in a better fit for $\lambda^{-1}\geq 5.5\,\rm \mu m^{-1}$. For the Serkowski spectra in HD 37903 and HD 161056, only RAT alignment by itself under the typical interstellar radiation field above the Lyman limit, within a typical cold neutral medium, can account for the observed spectra, with a combination of a very inefficient DG alignment. The capacity of our model to predict the starlight polarisation spectrum up to from infrared to far-UV is thus a promising tool for interpreting future missions that observe spectrophotometry in the UV bands.

Gas-phase chemical abundances are key observable consequences of galaxy evolution, being intrinsically tied to galaxy formation histories. Gas metallicity rises with increasing stellar mass ($\mathrm{M_*}$), forming the well-known mass-metallicity relation (MZR). MZR residuals have separately been shown to anti-correlate with star-formation rate (the “fundamental” metallicity relation), with gas mass and with optical size, but no single analysis has considered all trends together. We thus perform a combined analysis of all three trends, utilizing optical MaNGA integral field spectroscopy, HI-MaNGA gas masses, and MaNGA DynPop dynamical masses. We estimate inner gas masses for $\sim$1500 star-forming galaxies, finding this to be the most important parameter after $\mathrm{M_*}$ in predicting gas metallicities. We obtain equivalent results for stellar metallicities and gaseous N/O, suggesting that current inner gas masses are intrinsically linked to long-term chemical evolution histories. We show that more compact galaxies have lower dynamical masses, challenging suggestions that deeper gravitational potentials confer higher metallicities. We find a strong correlation between inner gas mass and galaxy size, meaning that short term inflow fluctuations cannot be responsible for the MZR residuals. With chemical evolution models, we show that our results can instead be explained by differences in long-term inflow histories. The earlier inflow histories of compact galaxies lead to lower gas masses and more rapidly declining gas reservoirs at late times, leading to higher metallicities. At fixed stellar mass, galaxy size therefore encodes the link between halo assembly histories, long-term gas inflow histories, current gas reservoirs and metallicity.

Using generalised Laplacian symmetry operators, we construct basis sets or Coulomb resolutions in several separable coordinate systems, including two R-separable systems. This expands the possible geometries in which basis set construction is feasible, a problem which is relevant to both galactic dynamics and computational chemistry. In particular we derive three basis sets (two in prolate spheroidal and one in cylindrical coordinates) which are expressible in closed-form using a single Jacobi polynomial. We also show how any spherical polar or prolate spheroidal basis set may be transformed into a bispherical or toroidal basis set.

We analyse the cluster content of the $\sim 4.5 \text{ deg}^{2}$ XMM-LSS and CDFS VIDEO fields which are expected to be partially covered by the upcoming MOONRISE survey. Using AMICO and WaZP photometric redshift-based cluster finders, we construct a sample of $519$ cluster candidates detected by both finders in the redshift range $z = 0.1-3$, including $74$ detections at $z > 1.5$. For all detections, we identify the Brightest Central Galaxy (BCG) and compute a list of probabilistic cluster memberships. Our photometric redshift measurements of the clusters agree well with spectroscopic redshifts from the literature, when available. From ancillary spectroscopic data, we assign $z_\text{spec}$ measurements to $116$ cluster candidates based on their spectroscopic members and to $204$ based on their likely BCGs. We also show that candidates containing Radio-Loud members are efficiently recovered using the prior-based cluster finder PPM. We perform a preliminary analysis of the galaxy content of these candidates, focusing on the Red-Sequence components of their apparent Colour-Magnitude Diagram. By comparing with models of galaxy evolution, we show that this population is consistent with a model of passive evolution with a formation at high redshift, and is already in place at $z = 1.5-2.0$. Finally, our cluster sample is used to evaluate how these clusters would be detected and characterised, according to various MOONRISE strategies. We show that cluster spectroscopic confirmation and characterisation could be efficiently achieved up to $z\sim1.7$ even with the shallowest survey strategy. This open unprecedented insight into the physical properties of high-redshift galaxy clusters and into galaxy formation in dense environments.

With Hubble Space Telescope (HST) and James Webb Space Telescope (JWST) observations of NGC 628 spanning 0.3–7.7 $\mu$m, we fit the spectral energy distributions (SEDs) of over 12,000 optically-selected star clusters, emerging young star clusters (eYSCs), and MIRI-selected sources with \textsc{cigale} to derive their ages, masses, extinctions, and dust properties. We find that near-infrared selected eYSC-I (compact Pa$\alpha$ and 3.3,$\mu$m PAH emission) and eYSC-II (compact Pa$\alpha$ and diffuse 3.3,$\mu$m PAH emssion) sources peak at $\sim$3–5 Myr, where $\sim 12\%$ of the clusters have an $E(B{-}V)>2$, demonstrating the presence of dust-embedded populations as clusters emerge. Further, the distributions of the fractional polycyclic aromatic hydrocarbon (PAH) abundance ($q_{\mathrm PAH}$) and stellar-to-nebular attenuation ratio ($E(B{-}V)_{\rm \star}/E(B{-}V)_{\rm neb}$) suggest an evolutionary sequence in which sources evolve from eYSC-I to eYSC-II as clusters clear their surrounding dust and gas. The photo-dissociation region (PDR) clearing timescale inferred from the ratio of eYSC-I to optically visible stellar clusters is $\sim$4 Myr. Additionally, we find that star clusters in the spiral arms of NGC 628 are preferentially more massive and more dust-reddened than those in inter-arm regions. Finally, we find that $\sim$65\% of eYSC-I, $\sim$27\% of eYSC-II, and $\sim$40\% of F335M-selected sources coincide with an F770W peak in our MIRI-selected catalog within 4 pixels, confirming that F770W-bright sources preferentially trace the youngest and dustiest regions. Overall, our results highlight the ability of JWST together with \textsc{cigale} model grids to identify and characterize eYSCs during their short-lived embedded phases, and provide constraints on the feedback mechanisms that govern the emergence of stellar clusters.

We investigate the molecular gas content and interstellar medium (ISM) conditions of REBELS-25, a massive, star-forming galaxy at $z=7.31$. Deep VLA Q-band and ALMA Band 3 observations reveal CO(3-2) and CO(7-6) emission (both at $\sim3.5\sigma$), and provide an upper limit on [C I](2-1). From the CMB-corrected CO(3-2) flux-representing the highest-redshift detection of a low-$J$ CO transition to date-we derive a molecular gas mass of $M_{\rm mol}=(1.0\pm0.4)\times10^{11}\,(\alpha_{\rm CO}/(3\,$M$_{\odot}$(K$\,$\kms$\,$pc$^2)^{-1}))\,$M$_{\odot}$, directly confirming the presence of a very massive gas reservoir only $\simeq700\,$Myr after the Big Bang. This implies an extreme gas fraction of $f_{\rm gas}\simeq0.95$, a gas-to-dust ratio of $\delta_{\rm GDR}\simeq6\times10^2$, and a depletion timescale of $\tau_{\rm dep}\simeq1.2\,$Gyr, broadly consistent with extrapolated scaling relations for main-sequence galaxies at lower redshift. Using the radiative transfer code TUNER, we self-consistently model CO and dust continuum emission in the context of the significant CMB background, constraining ISM properties and recovering $M_{\rm mol}= (1.8^{+1.0}_{-0.9})\times10^{11}\,$M$_{\odot}$, independent of assumptions about $r_{31}$ and $\alpha_{\rm CO}$. We further discuss the use of alternative molecular gas tracers at early epochs. Combining CO and [C II] measurements, we infer an empirical [C II]-to-H$_2$ conversion factor of $\alpha_{\rm [C II]}=(60\pm25)\,$M$_{\odot}$/L$_{\odot}$, suggesting [C II] remains a viable molecular gas tracer in the Epoch of Reionization. These results demonstrate the detectability of low-$J$ CO emission even at $z>7$, paving the way for next-generation facilities, and provide critical insights into the rapid mass assembly of galaxies during the first billion years of cosmic history.

Context. Active galactic nuclei (AGN) jets are fundamental drivers of galaxy evolution, injecting kinetic energy into their environments. The large-scale morphology and spectral properties of these radio galaxies are consequences of complex particle acceleration, energy loss, and absorption processes. While the shape of the synchrotron spectrum encodes the plasma's energetic history, understanding the physics of particle acceleration and duty cycles has historically been limited by a lack of well-resolved observations at ultra-low frequencies (< 100 MHz), where the oldest cosmic ray electron populations are traced. Aims. This study aims to perform the first comprehensive multi-frequency analysis of bright extended radio galaxies down to 58 MHz. The goal is to study electron acceleration mechanisms, accurately measure the low-frequency spectral shape, and constrain the injection index for a sample of Fanaroff-Riley (FR) I and II galaxies using spectral ageing models. Methods. Utilising new 58 MHz observations from the LOFAR Low Band Antenna (LBA) combined with LOFAR High Band Antenna (HBA; 144 MHz) and Rapid ASKAP Continuum Survey (RACS, 887 MHz & 943.5 MHz & 1367.5 MHz) data, a sub-sample of 22 extended sources from the 3CRR catalogue was selected, requiring the largest angular size to be at least 2.5'. The analysis involves constructing detailed spectral index maps and utilising radio colour-colour diagrams to interpret spectral shapes and constrain ageing model parameters across the radio lobes. Results. This study presents the ultra-low frequency spectral index maps for this sample. For FR I galaxies, spectral indices range from  0.5 near the core (consistent with first-order Fermi acceleration) to > 1.0 in the lobes. For FR II galaxies, hotspots exhibit steep low-frequency spectra (0.5 - 0.9), suggesting complex acceleration or absorption effects.

Protostellar jets are one of the earliest forms of stellar feedback, but their impact on star formation and cluster assembly in centrally concentrated molecular clouds remains poorly understood. We study how protostellar jets affect the star formation efficiency, the temporal variability of star formation, star cluster structure, and the early dynamical state of centrally concentrated, newly forming star clusters using the Torch star cluster formation framework. We adopt a centrally concentrated initial cloud model with mass M = 2.5 x 10^3 solar masses and compare six pairs of simulations with and without protostellar jets, supplemented by one additional higher resolution pair of simulations. We analyze our simulations using global star formation diagnostics together with structural and dynamical measures of the stellar population. Models with jet feedback achieve star formation efficiencies of 12-16%, while the corresponding models without jets yield higher efficiencies of 19-33%. Jets also cause star formation to occur in discrete bursts rather than continuously, to produce more extended and substructured stellar systems, and to leave behind stellar populations that are less tightly bound and have higher virial parameters. In our centrally concentrated initial conditions, runs with jets form stellar systems that better reproduce the observed range of the projected structural parameter Q_2D in young clusters than runs without jets, indicating that protostellar jets are an important early feedback channel even in centrally concentrated clouds that regulates star formation efficiencies and shapes the emerging cluster structure.

We developed the third-order hydrodynamic equations using relativistic extended thermodynamics of gases with 14 independent fields. The resulting fluid equations are based on the relativity principle, the entropy principle, and the requirement of hyperbolic, and hence finite, propagation of disturbances, which is automatically incorporated. The expressions of entropy, four-current, shear-stress tensor, dynamic pressure, and heat flux are expanded up to third order (cubic). We explicitly present the newly calculated coefficients in the equilibrium properties of an ultra-relativistic gas regime and the non-degenerate relativistic gas. Contrary to the general cases, the non-degenerate regime eliminates fugacity from the coefficients, allowing for the easy normalization of these coefficients, and the ultra-relativistic regime provides us with the upper bounds of these coefficients. We found good agreement on some of the coefficients as compared to calculations from earlier models, specifically in kinetic theory, and other coefficients had slightly different values to those obtained in kinetic theory.

The Large Magellanic Cloud (LMC), the most massive satellite of the Milky Way (MW), can significantly perturb the local dark matter (DM) distribution. We study its impact on directional DM detection using the Auriga cosmological simulations of a MW analogue hosting an LMC analogue. We find that the LMC induces strong anisotropies in directional recoil signals, driven primarily by the non-zero mean azimuthal velocity of the local DM distribution. The characteristic ring-like feature predicted in the Standard Halo Model (SHM) for heavy DM and low recoil energies is strongly distorted, producing an asymmetric recoil pattern concentrated at preferred azimuthal angles. Differences between recoil maps for the MW-LMC analogue and the SHM reach up to $\sim80\%$ near the signal maximum. These distortions significantly enhance directional discovery prospects, reducing the number of events required to reject isotropy by nearly a factor of five for a 100 GeV DM particle in a near-future CYGNUS-like experiment, and by even larger factors for heavier DM. Our results highlight the importance of the LMC for interpreting and optimizing future directional DM searches.

Primordial black holes (PBHs) in the asteroid-mass window ($10^{17}$-$10^{22}$ g) can account for all of the dark matter without violating any observational constraint, yet are routinely dismissed as fine-tuned. I put that dismissal to the test by applying three complementary fine-tuning measures uniformly across a broad landscape: three non-inflationary PBH production mechanisms, six classes of inflationary PBH models, and seven particle dark matter benchmarks, all evaluated against the same observable target. Three distinct naturalness universality classes emerge, determined entirely by the analytic structure of the abundance map rather than by the nature of the dark matter candidate. Biased-domain-wall PBHs are as natural as off-resonance weakly interacting massive particles and freeze-in particles; early-matter-domination and first-order phase transition PBH mechanisms occupy an intermediate tier alongside coannihilating WIMPs, unified by a structural identity in which the fine-tuning measure equals the logarithm of the ratio of the formation scale to the matter-radiation equality scale; and single-field ultra-slow-roll inflationary collapse is severely tuned for a distinct reason: a double exponential in which the power spectrum amplitude is itself exponentially sensitive to the inflaton potential coefficients, on top of the exponential collapse sensitivity of the abundance map. My main conclusion is that {\em the claim that PBH dark matter is generically fine-tuned conflates the worst case with a landscape spanning every naturalness tier}. The three-measure protocol also resolves a tension in the recent literature: the Barbieri-Giudice and Iovino-Riotto fine-tuning measures answer complementary questions and are reconciled within the two-layer decomposition developed here.

The anisotropy of the redshift space bispectrum depends upon the orientation of the triangles formed by three $\mathbf{k}$ modes with respect to the line of sight. For a triangle of fixed size ($k_1$) and shape ($\mu,t$), this orientation dependence can be quantified in terms of angular multipoles $B_\ell^m(k_1,\mu,t)$ which contain a wealth of cosmological information. We propose a fast and efficient FFT-based estimator that computes the bispectrum multipole moments $B_\ell^m$ of a 3D cosmological field for all possible $\ell$ and $m$ (including $m\neq 0$). The time required by the estimator to compute all multipoles from a gridded data cube of volume $N_g^3$ scales as $\sim \mathcal{O}(N_g^4)$ in contrast to the direct computation technique which requires time $\sim \mathcal{O}(N_g^6)$. Here, we demonstrate the formalism and validate the estimator using a simulated non-Gaussian field for which the analytical expressions for all the bispectrum multipoles are known. The estimated results are found to be in good agreement with the analytical predictions for all $16$ non-zero multipoles (up to $\ell= 6, m=6$). We expect the $m \neq 0$ bispectrum multipoles to significantly enhance the information available from galaxy redshift surveys and future redshifted 21-cm observations.

We use the inferred evolution of Lyman-$\alpha$ luminosities of galaxies in the redshift range of $z \sim 9-16$ from the James Webb Space Telescope (JWST) data to predict the power spectrum of Loeb-Rybicki haloes formed by the scattering of Lyman-$\alpha$ photons from neutral hydrogen gas in the intergalactic medium expanding with the Hubble flow, until they Doppler shift out of resonance and escape towards the observer. This leads to the formation of the so-called Loeb-Rybicki intergalactic haloes, which are expected to be prominent even before the epoch of reionization. We find excellent prospects for a statistical detection of the intensity mapping signal from the clustering of these haloes, with current and future experiments such as the SPHEREx and CDIM. We also describe the detectability of the signal in cross-correlation with the 21-cm emission from the neutral hydrogen in the intergalactic medium at these redshifts. We find that the cross-correlation signal should be detectable at a significance of a few to a few tens of standard deviations out to $z \sim 13$ and marginally out to $z \sim 16$, using the above experiments in combination with the Square Kilometre Array (SKA)-LOW and its pathfinder, the Murchison Widefield Array (MWA).

With over two dozen detections in the Milky Way, double neutron stars (DNSs) provide a unique window into massive binary evolution. We use the POSYDON binary population synthesis code to model DNS populations and compare them to the observed Galactic sample. By tracing their origins to underlying single and binary star physics, we place constraints on the detailed evolutionary stages leading to DNS formation. Our study reveals a bifurcation within the well-known common envelope (CE) formation channel for DNSs, which naturally explains an observed split in the orbital periods of the Galactic systems. The two subchannels are defined by whether the donor star has a helium core (Case B mass transfer) or a carbon-oxygen core (Case C) at the onset of the CE, with only the helium core systems eventually merging due to gravitational wave-modulated orbital decay. We find that across different treatments of the CE phase, the formation of DNSs through both subchannels requires either a generous core definition of $\simeq$ 30% H-fraction or a high CE ejection efficiency of $\alpha_{\rm CE}\gtrsim1.2$. By testing different supernova kick velocity models, we find that galactic DNSs are best reproduced using a prescription that favors low velocity kicks ($\lesssim 50 \, \rm km/s$), in agreement with previous studies. Furthermore, our models indicate that merging DNSs are born from a stripped progenitor with a median pre-supernova envelope mass $\sim$ 0.2$M_{\odot}$. Our results highlight the value of detailed evolutionary models for improving our understanding of exotic binary star formation.

We establish a connection between the stochastic self-propagating star-formation model and spatio-temporal point processes by showing that, under suitable discretisation, the SSPSF update law can be represented by a separable spatio-temporal Hawkes process. Building on this connection, we propose a spatio-temporal log-Gaussian Cox-Hawkes process as a continuous point-process model for stochastic star formation. The model represents star-formation events as point patterns driven jointly by deterministic galactic structure, latent spatio-temporal background variation, and dependence on past events. Its key feature is that the deterministic mean field, latent Gaussian random field, and history-dependent interaction field enter through a single log-intensity. This log-scale construction differs from additive Cox-Hawkes formulations and allows the history effect to be signed: past events may either increase or decrease future local intensity while the conditional intensity remains positive. The resulting framework provides an interpretable point-process model for representing latent clustering, self-excitation, local inhibition, and event-driven propagation in stochastic star formation. Beyond linking SSPSF to spatio-temporal point-process theory, it offers a continuous stochastic formulation for analysing the propagation of star formation in galaxies and for interpreting observational surveys of star-forming regions within a unified statistical model.

We present integral field spectroscopic observations of 76 strong gravitational lens candidates identified with a residual neural network in the DESI Legacy Imaging Surveys, obtained with the Multi Unit Spectroscopic Explorer (MUSE) on the ESO's Very Large Telescope. These observations are part of an ongoing effort to build a large, spectroscopically confirmed sample of strong lensing systems for studies on dark matter, galaxy structure, and cosmology. Our MUSE program targets both lens and source redshifts, with particular emphasis on southern hemisphere systems. MUSE's wide spectral coverage and integral field capability allow for efficient identification of multiple sources, lens environments, and weak spectral features. Redshifts for lenses and sources were obtained via manual identification of spectral features in extracted 1D spectra. Our dataset includes systems with complex configurations, such as multiple source planes and group or cluster-scale environments. We extracted and analyzed 223 spectra, successfully determining both the lens and the source redshifts for 55 gravitational lens systems. For an additional 15 targets, we measured the redshifts of the lenses but were unable to determine the redshifts of the background sources. Six targets were confirmed to not be gravitational lenses. The results presented here complement space-based imaging from our HST SNAPshot program and spectroscopic follow-up with DESI and Keck, and have lasting legacy value for identifying interesting high redshift sources and complex lensing configurations.

Galaxy mergers play an important role in galaxy evolution. Therefore, accurate merger identifications are paramount for achieving a complete understanding of how galaxies evolve. As we enter the era of large, deep, high-resolution imaging surveys, we can observe mergers extending to even lower masses and higher redshifts. Despite low-mass galaxies being more common, many previous merger identification methods were calibrated for high-mass galaxies, which are easier to identify. To prepare for upcoming surveys, we train a convolutional neural network (CNN) using mock $\textit{HST}$ CANDELS images at $z\sim1$ created from the IllustrisTNG50 cosmological simulation. We successfully identify galaxy mergers between a wide range of galaxies ($10^8M_\odot < M_\star < 10^{12.5}M_\odot$, and $q\geq1:10$), achieving overall accuracy, purity, and completeness of $\sim65\%$. We show, for the first time, that a CNN trained on this diverse set of galaxies is capable of identifying major mergers, especially at early stages (74% accuracy), similar to that of networks trained at lower redshifts and/or higher masses (with accuracies between $66-80$%). We discuss the inherent limits of galaxy merger identification due to orientation angle, finding 98% of mergers are correctly identified from at least one angle, and 61% from the majority of angles. We additionally explore the confounding variables, such as star formation, to consider when applying to real this http URL network enables the exploration of the impact of previously overlooked mergers of high mass ratio and low stellar masses on galaxy evolution in CANDELS, and can be expanded to surveys from $\textit{JWST}$, Rubin, $\textit{Roman}$, and $\textit{Euclid}$.

Among the emerging excess of massive, bright galaxies at Cosmic Dawn $z \gtrsim 9$ seen by the James Webb Space Telescope, several exhibit spectral features associated with active galactic nuclei (AGN). These AGN candidates suggest that supermassive black holes (SMBHs) grow rapidly in the early Universe. In a series of numerical experiments, we investigate how SMBHs grow within and influence the most massive galaxies at Cosmic Dawn using cosmological hydrodynamic zoom-in simulations run with the adaptive mesh refinement code \textsc{ramses}. Our suite of simulations explore how super-Eddington accretion, seed mass, and the strength of feedback influence SMBH-galaxy co-evolution in the most massive galaxies ($M_\star \gtrsim 10^8 M_\odot$) of the early Universe ($z \sim 15 - 9$). The environment which our numerical experiments reside in is an overdensity that collapses into a $\sim 10^{11} M_\odot$ halo by $z \sim 9$. Within this type of environment we find that SMBH growth is sensitive to stellar feedback which generates a turbulent-multiphase interstellar medium (ISM) that stochastically starves the SMBH. In the absence of AGN feedback, we find that the SMBH is starved $\sim 50\%$ of the time after the onset of star formation in the galaxy. SMBH growth can become self-regulated by AGN feedback if the SMBH becomes massive enough, either by accretion or seeding, for its feedback to dominate the surrounding nuclear region. We find no evidence of galaxy-scale, AGN-driven quenching in the star formation rate (SFR) across all simulations in our suite.

Little Red Dots, discovered by the James Webb Space Telescope, are hypothesized to be active galactic nuclei containing a supermassive black hole, possibly surrounded by a dense stellar cluster, large amounts of gas, and likely by a population of stellar-mass black holes. We develop a simple nuclear star cluster model to evolve the rapid mass growth of black hole seeds into the supermassive regime. The combined processes of tidal disruption events, black hole captures, and gas accretion are accounted for self-consistently in our model. Given the observed number density of Little Red Dots, and under reasonable assumptions, we predict at least a few tens of tidal disruption events and at least a few black hole captures at z=4-6, with a tidal disruption event rate an order of magnitude larger than the black hole capture rate. We also estimate the uncertainties in these estimates. Finally, we comment on the low x-ray luminosity of Little Red Dots.

The dark matter flux in a direct detection experiment depends on its local speed distribution. This distribution has been inferred from simulations of Milky Way-like galaxies, but such models serve only as proxies, given that no simulation directly captures the detailed evolution of our own Galaxy. This motivates alternative approaches that obtain this distribution directly from observations. In this work, we utilize 98 Milky Way analogues from the TNG50 simulation to develop and validate a procedure for inferring the dark matter speed distribution using the kinematics of nearby stars. We find that the dark matter that originated from old mergers, plus that from recent nonluminous accretions, is well described by a Maxwell-Boltzmann speed distribution centered at the local standard-of-rest velocity. Meanwhile, recently accreted dark matter from massive mergers has speeds that can be traced from the associated stellar debris of these events. The stellar populations systematically underestimate the velocity dispersion of their dark matter counterparts, but a simple kinematic boost brings the two into good alignment. Using the TNG50 host galaxies, we demonstrate that combining these two contributions provides an accurate reconstruction of the local dark matter speeds. As an application of the procedure to our own Galaxy, we utilize stellar kinematic data from Gaia to quantify how the dark matter remnants from the Milky Way's last major merger impact its speed distribution in the solar neighborhood.

Due to their different nucleosynthetic origin, a stellar population produces oxygen (O) and iron (Fe) on different timescales and their relative abundance can deviate strongly from solar. Galaxy formation models should treat these elements separately, as they play a distinct role in shaping physical phenomena. For example, oxygen mainly sets the gas cooling rate, while the iron abundance sets stellar atmosphere opacities impacting stellar evolution, spectra and feedback. Observations of star-forming galaxies usually only constrain gas-phase oxygen abundance, vastly limiting our capabilities of separating the cosmic evolution of oxygen and iron. Here, we present an observationally-motivated framework to scale the cosmic evolution of O and Fe abundances. We apply the relation between the alpha-enhancement and galaxies' specific star formation rate ([O/Fe]-sSFR; Chruslinska et al. 2024) to derive the Fe and O-dependent cosmic star formation history (cSFH). We find that star formation with near-solar O/Fe is rare: at least 70% of the integrated cosmic stellar mass forms at non-solar O/Fe. The cosmic average metallicity is generally lower in [Fe/H] than in [O/H] (by up to a factor 3), with the offset increasing from redshifts z=0 to z 3 and then approaching the core-collapse O/Fe ratio. We validate our results against samples that probe the Fe-dependent cSFH in different regimes such as absorption-derived <[Fe/H]> from long gamma-ray bursts. Our results impact the interpretations of stellar and galaxy spectra and the predicted rates of transients, especially those linked to metal-poor progenitors (e.g., black hole mergers).

Interstellar objects (ISOs), particularly those with cometary activity, provide unique insight into the primordial physical and chemical conditions present during the formation of the planetary system in which they originated. Observations in the sub-mm regime allow for direct measurements of several parent molecules released from the comet nucleus into the coma. Here we present observations of the third ISO, 3I/ATLAS, with the `Ū`ū heterodyne receiver on the James Clerk Maxwell Telescope (JCMT), which targeted emission from HCN($J = 3 - 2$) and CO($J = 2 - 1$). Our observations, taken between 16 July 2025 and 21 July 2025 (UT), when 3I/ATLAS was at a heliocentric distance between 4.01 and 3.84 au, provide the earliest sub-mm constraints on its activity. We do not detect HCN or CO in these epochs, with 3$\sigma$ upper-limits on the production rates of $Q(HCN) < 1.7 \times 10^{24}$ s$^{-1}$ at $r_h = 4.01 - 3.97$ au and $Q(CO) < 1.1 \times 10^{27}$ s$^{-1}$ at $r_h = 3.94 - 3.84$ au, respectively. We combine this HCN limit with later JCMT observations of HCN to constrain its temporal evolution. Fitting the HCN detections with a $Q(HCN) \propto r_h^{-n}$ model and accounting for the upper-limits yields $n = 12.7^{+6.9}_{-2.5}$. This slope is steeper than those of typical Solar System comets, but consistent with the production rate slopes measured for other species in the coma of 3I/ATLAS.

We present argon abundances from a sample of 46 star-forming galaxies at $z=2-3.5$ from the Assembly of Ultradeep Rest-Optical Observations Revealing Astrophysics (AURORA) program. Although argon is an $\alpha-$element produced by Core Collapse Supernovae (CCSNe), the latest supernova yield models suggest additional argon production and enrichment by Type Ia supernovae (SNe Ia), unlike other $\alpha-$elements such as oxygen. To study the relationship between argon and oxygen abundances, we construct two median-stacked composite spectra for separate $z=2.0-2.6$ and $z=2.8-3.5$ redshift bins, presenting more representative measurements than previous samples that require individual detection of faint lines. Abundance ratios were determined using an empirical calibration based on the strength of the [ArIII]$\lambda7137$ emission line relative to the [OIII]$\lambda5008$ emission line. With this calibration, we estimate argon abundances (Ar/O) of $0.42^{+0.12}_{-0.10}\mbox{(Ar/O)}_{\odot}$ for the $\langle z \rangle = 2.26$ bin and $0.42^{+0.12}_{-0.11}\mbox{(Ar/O)}_{\odot}$ for the $\langle z \rangle = 3.15$ bin, suggesting minimal SNe Ia and dominant CCSNe enrichment in this sample. Comparison of our abundance measurements of $z\sim 2-3$ AURORA galaxies with chemical evolution modeling of Milky Way stars shows consistency with the Milky Way Bulge component, suggesting a rapid star-formation timescale. However, even larger samples of actively star-forming galaxies with available argon abundances, as well as comparisons between argon abundance and other critical galaxy properties (e.g., sSFR) and models (e.g., one tuned specifically to this redshift range) are needed to draw stronger conclusions on the role of argon in galactic chemical enrichment at Cosmic Noon.

Active galactic nuclei (AGN) hosting water megamaser disks provide precise geometric measurements of black hole mass, distance, maser disk orientation, and dynamical center. In anticipation of space-based very long baseline interferometry, these systems offer a path to black hole shadow (BHS) imaging beyond Sgr A* and M87*. We present new Submillimeter Array continuum observations of water megamaser galaxies, supplemented by archival ALMA and VLA measurements, to assess whether their AGN cores are bright enough for BHS-scale imaging. For a 21-source parent sample, we map the predicted BHS diameters of systems with published SMBH masses to submillimeter/millimeter (submm-mm) baseline requirements, estimate AGN core flux densities at 230 GHz while bounding thermal dust and extended-jet contamination and checking whether variability could affect the continuum estimates, and evaluate the astrometric precision required to detect spin-dependent BHS offsets for NGC 4258. NGC 4258 is the only source resolvable on Earth-L2 baselines; other targets require longer baselines approaching Earth-L4/L5 distances, and only a handful have $S_{230} \gtrsim 10$ mJy beam$^{-1}$. We also find a submillimeter excess in NGC 4258, suggesting that its disk remains geometrically thin to $\lesssim 100$ Schwarzschild radii before transitioning to an advection-dominated flow. Even for maximal spin, the formal 230 GHz BHS centroid precision is not the limiting term: the measurement would require locating the 22 GHz water maser dynamical center and registering it to the 230 GHz BHS image roughly seventy times more precisely than current maser astrometry allows, making the spin-offset measurement infeasible with present data.

Extended warm and cold gas nebulae, with complex morphologies and kinematics, have been observed in the centres of cool-core galaxy clusters. Their origin within the hot intracluster medium (ICM) is still puzzling, and among many mechanisms, positive feedback from the central active galactic nucleus (AGN) has been proposed. In this work, we performed a suite of very high-resolution hydrodynamic simulations of a Perseus-like cool-core galaxy cluster subject to self-regulated AGN jet feedback, which leads to realistic ICM properties. By explicitly following warm ionized, neutral, and molecular gas phases, we studied the complex interplay between AGN activity and the multi-phase ICM. While AGN feedback globally heats the ICM, we find that during the individual AGN jet bursts, hot material is also injected laterally to the jet axis, within the turbulent mixing layer. This material, as it expands, compresses the surrounding hot ICM, reducing the local cooling time, and leads to the formation of cold clumps on a characteristic timescale of $\sim 30$ Myr. By employing tracers, we explicitly track cooling within the affected regions, finding that very hot gas identified in high-compression, low-vorticity zones condenses in situ to form cold clumps. A statistical analysis reveals that the condensation of cold gas is highly promoted once the local turbulent Mach number, $\sigma_{hot}/c_{s,hot}$, in the hot gas component ($T \geq 10^7$ K) takes values around  0.3. The presented process is a further important step in understanding the physical mechanisms that lead to the formation of cold gas in the cluster core. Our measured values of the characteristic turbulent Mach number, together with detailed multi-phase gas kinematics predictions, provide important theoretical tools to interpret future X-ray spectroscopy and deep radio data, ultimately to constrain the origin of cool-core cluster nebulae.

Self-interacting dark matter (SIDM) with a sufficiently large cross section has been shown to naturally produce constant dark matter (DM) cores, as well as core-collapse, at the centers of dwarf halos on cosmic timescales, potentially reducing tensions with observation. Here, we present halos from a new dark matter only (DMO) cosmological (SIDM) simulation: this http URL DMO with a velocity-dependent self-interaction cross section with $\sigma/m_\text{max} = 50$ cm$^2$/g at $v_\text{max} = 35$ km/s. We compare these to the CDM suite of Storm simulations including both DMO and dark matter + hydrodynamics runs, in order to test core-formation (and core-collapse) across different dark matter models. We show that this http URL DMO can reproduce core slopes consistent with observations of isolated dwarf galaxies and more massive ($\text{M}_{vir} \gtrsim 10^{10} M_{\odot}$) CDM dwarf halos that include stellar feedback from the matched CDM run (Storm CDM+baryons). We identify nine this http URL SIDM DMO halos in the core-collapse phase of gravothermal evolution with halo masses below $2\times 10^9 M_{\odot}$. We find that using core slope to measure the core-collapse timescales of this http URL DMO halos agrees well with predicted collapse times estimated with the parametric model for SIDM halos introduced by Yang et al.(2023). Additionally, compared to central density, core slope is less sensitive to both the radius of measurement and halo merger history. These results indicate that the slope of the inner DM density profile more cleanly differentiates core-collapsed versus core-forming halos than central density amplitude.

Finding the first generation of (Population III or Pop III) stars is one of the most ambitious and exciting challenges of astrophysics. JWST opened concrete prospects for their detection during the Epoch of Reionization (EoR), where increasing evidence suggests that residual Pop III formation may persist, even within pristine pockets of high-mass halos, due to inhomogeneous enrichment. However, the identification of Pop III stars within globally enriched environments will be challenging. We investigate the detectability of a subdominant Pop III component in/around massive ($M_\star \gtrsim 10^9 ~\mathrm{M_\odot}$) galaxies at $z \approx 6.5 - 9$ from the dustyGadget cosmological simulation suite, and the confusion arising from second-generation (Pop II) stars in their surroundings. We find that young ($\lesssim 1$ Myr), massive ($M_\mathrm{III} \sim 6 \times 10^5 ~\mathrm{M_\odot}$) Pop III clusters forming within these galaxy environments are responsible for strong HeII1640 line emission ($L_\mathrm{HeII1640} \gtrsim 10^{41} ~\mathrm{erg \, s^{-1}}$), which would be detectable with $\approx 10 (50)$ h of medium-resolution observations with NIRSpec/IFU at $z \approx 6 (10)$. These bright luminosities cannot be produced by standard Pop II populations alone. On the other hand, the dominant Pop II component within massive “hybrid” Pop III hosts powers strong metal line emission ($L_\mathrm{[OIII]5007} \gtrsim 10^{42} ~\mathrm{erg \, s^{-1}}$), indicating that the detection of metal lines alone cannot exclude the presence of Pop IIIs in high-$z$ galaxy environments. We further discuss candidate selection strategies based on Ly$\alpha$, H$\alpha$ and H$\beta$ emission, and how spatially resolved observations may enable the detection of isolated, pristine pockets in the outskirts of massive halos.

We present high-angular-resolution ($\sim0.05^{\prime\prime}$, $\sim 60-250$ au) ALMA Band 6 observations from the HOTDISK project (Hot-Origin Tracer survey of DISKs of massive protostars) aimed at investigating the "hot-disk" chemical pattern traced by vibrationally excited water, NaCl, SiS, and SiO in the innermost regions around massive protostars. Ten targets were selected based on strong CH$_3$CN emission exhibiting clear rotational signatures and centrally concentrated SiO emission from lower-resolution observations. We detect vibrationally excited water emission toward 7 of the 10 sources. In all detections, the blueshifted and redshifted components are compact and located on opposite sides of the 1.3 mm continuum peak, with velocity gradients approximately perpendicular to the outflow axes, consistent with rotation on disk scales. Emission from NaCl and SiS is detected toward 5 of these 7 sources and exhibits similar kinematics, further supporting the presence of compact rotating structures. In contrast, commonly used hot-core tracers (e.g., CH$_3$CN and SO$_2$) primarily probe larger-scale envelope gas. These results demonstrate that vibrationally excited water, NaCl, and SiS are powerful tracers of disk structures on $\sim$100 au scales, when observed at sufficient angular resolution and sensitivity. The high detection rate suggests that hot-disk chemical patterns – and thus compact rotating disks – are common in massive star-forming regions, at least among sources with well-developed rotating envelopes.

Metallicity can play a significant role in massive binary evolution through its impact on the opacity within stellar interiors and wind-driven mass loss. In this work, we investigate how the double neutron star (DNS) delay time distribution (DTD) is shaped by the metallicity-dependent evolution of the helium star$-$NS progenitor system. Drawing from insights rooted in single and binary star physics, we argue that at a given metallicity, the stellar radius during the helium main-sequence sets a lower limit on the size of the DNS orbit at birth. We then perform population synthesis with the detailed binary evolution code POSYDON to illustrate the resulting DTD across a range of metallicities. Our results indicate that, independent of the common envelope efficiency and reasonable natal kicks, the majority of DNS mergers across metallicities occur typically no earlier than $\simeq 40\,\rm{Myr}$ after star formation and peaks strongly between $80-250\,\rm{Myr}$. Roughly $15\%$ of DNSs merge within 80 Myr, which may explain $r$-process enrichment in environments with brief star formation histories, while $\gtrsim 20\%$ merge on delay times $>1$Gyr, providing an explanation for short gamma-ray bursts in old, metal-poor galaxies. The shape of the DTD can be complex, with a metallicity-dependent split in the dominant formation channel imprinting a characteristic double-peaked structure. Although ideally oriented natal kicks can produce very short merging DNS, we find that the required kick magnitudes are inconsistent with observations. Our work has implications for assessing the contribution of DNS mergers to $r$-process enrichment and gamma-ray bursts/kilonovae transients across cosmic time.

Calibration checks for neural posterior estimators in gravitational-wave inference should remain interpretable when observations contain data-quality artifacts. We study marginal interval calibration in a controlled frequency-domain binary-black-hole benchmark with synthetic glitches, frequency masks, and power-spectral-density mismatch. The posterior sampler is a support-aware conditional flow-matching estimator with a circular representation of coalescence phase. We compare raw marginal credible intervals with global rescaling, oracle artifact-stratified rescaling, hard predicted-label rescaling, and soft learned artifact-aware interval rescaling (LAIR). In the 1024-bin evaluation, a single global scale fitted on mixed calibration data transfers poorly to frequency-mask cases, giving MA90CE = 0.1195. Soft LAIR lowers the corresponding error to 0.0672, but it is not uniformly better than the raw FMPE intervals. A 40-seed LAIR audit and a six-checkpoint FMPE training-seed audit show that the frequency-mask behavior is not a single-split artifact. The classifier recognizes frequency masks and PSD mismatch reliably, while glitch recall remains low. Waveform-resolution tests, PyCBC/LAL TaylorF2 backend checks, prior and Gaussian baselines, and controlled-likelihood reference-posterior probes indicate that marginal coverage must be read together with posterior width, geometry, and likelihood-based diagnostics. In this benchmark, LAIR is therefore best viewed as an artifact-structured interval diagnostic rather than as a substitute for posterior validation.

The OSIRIS-REx mission acquired spatially resolved (2-10 m spot sizes) visible-near infrared (VNIR) and thermal infrared (TIR) spectra across four candidate sampling sites on asteroid (101955) Bennu: Nightingale, Osprey, Sandpiper, and Kingfisher. To quantify heterogeneity across a small body (about 500 m radius) like Bennu, we explore remotely observed spectral data to draw conclusions about the mineralogical composition and key physical processes that drive surface variability. We derive diagnostic band parameters from the OSIRIS-REx Visible and Infrared Spectrometer and the OSIRIS-REx Thermal Emission Spectrometer datasets to quantify compositional and physical variability across sites and assess their mineralogical context. The VNIR spectra exhibit similar overall reflectance shapes but systematic differences in spectral slopes and the 2.74 micron OH absorption. TIR emissivity spectra reveal modest but statistically significant shifts in the Christiansen Feature, silicate stretching, and bending band positions, indicating differences in silicate composition, hydration state, and Mg/Fe relative abundance. Principal component analysis separates each site into distinct clusters in multivariate band-parameter space, whereas K-means clustering identifies intra-site spectral sub-populations. Welch's Analysis of Variance and Hotelling's tests confirm that band-parameter variations between sites are significant. These results reveal that Bennu's surface preserves measurable spectral heterogeneity at 2-10 m scales, with site-to-site variations in hydration indicators and silicate band positions. The spectral properties of Nightingale encompass the full range observed across all four sites, establishing a remote sensing baseline for contextualizing laboratory analyses of the returned sample within Bennu's broader composition diversity and alteration history.

This paper describes lessons learned over a long career building astronomical instruments. Although these lessons are based on one person's experiences, they were learned while working on major and minor instrument projects with many scientific and engineering colleagues. When I interact with colleagues entering the field, I am reminded that the "obvious" approaches were in fact learned over many years by observing others develop instruments, by receiving good advice, and by making mistakes. I hope to help others avoid making all of the same mistakes.

The Habitable Worlds Observatory (HWO) aims to enable the detection and characterization of Earth-like planets around Sun-like stars to search for possible signs of life elsewhere in our universe. This requires an incredibly sensitive coronagraph instrument that suppresses the light from the star by a factor of 10 billion, which must contend with error terms that have not previously limited high-contrast instrumentation at lower levels of starlight suppression. Polarization aberrations are one such source of error that is particularly problematic for coronagraphy on a large space telescope. Optical rays in large, compact astronomical observatories can have large changes in angle of incidence over the beam, which induce polarization aberrations that decrease sensitivities to faint signals at small angular separations. Limiting variation in angles of incidence along the optical path could lead to longer, less stable observatories. This could negatively impact the total number of exo-Earths HWO would be able to detect. This study links open-source physical optics modeling tools to an exoplanet yield optimizer to understand how polarization aberrations influence science return for HWO. We also explore how polarization aberrations scale with change in angle of incidence, which could drive the primary-secondary mirror distance and overall observatory stability. In the visible, we find that decreasing the EAC-1 barrel from 16m to 12m results in $\approx 10^{-10}$ contrast at the IWA where we expect exo-Earths to be. In the UV we appear to be less sensitive to polarization because exo-Earths are farther from the IWA. We also find a limited range over which the design reference mission of EAC-1 can be optimized to compensate for polarization aberrations using altruistic yield optimization. We then report on mitigation strategies to minimize the presence of polarization aberrations in HWO.

The standard interpretation of observed Stokes profiles to infer the physical conditions of the solar atmosphere is inherently an ill-defined problem due to observational noise and mathematical degeneracies. Traditional pixel-by-pixel (1D) inversion codes provide point estimates with unreliable uncertainties, at the expense of significant computational time. Recent machine-learning-based Bayesian frameworks are restricted to 1D spatial configurations, ignoring crucial spatial correlations between neighboring pixels. We aim to develop a novel multidimensional inversion framework capable of performing fast and scalable Bayesian inference across an entire 2D field-of-view (FoV). This approach seeks to provide accurate height-dependent atmospheric parameters with reliable posterior distributions while exploiting spatial correlations. We introduce a new generative modeling strategy based on conditional flow matching. The model utilizes multi-scale spatial features extracted from observed Stokes profiles in the Fe I line pair at 630 nm, which then conditions a flow matching generative model to sample from the complex posterior distribution of the atmospheric parameters. The framework is trained using realistic 3D quiet Sun magnetohydrodynamic simulations. Validation on independent synthetic datasets demonstrates that the model accurately captures the true 3D stratification of all thermodynamic and magnetic parameters. Because the code additionally provides a geometrical height scale, it allows for the computation of 3D electric current density maps, Lorentz forces, and Ohmic and ambipolar dissipation maps in the solar photosphere. Application to real Hinode/SP quiet Sun observations yields highly localized electric currents at magnetic boundaries. We also leverage the 3D geometrical information to trace the emergence of small-scale emerging magnetic loops across the solar atmosphere.

Binary asteroids provide key access to fundamental parameters of Solar System remnants and planetary formations. However, the current knowledge of binary asteroids remains strongly biased by observational limitations, and main belt binary systems are still poorly characterised since current techniques preferentially detect either widely separated binaries close and bright systems. In this context, the high-precision astrometry of the Gaia mission has revealed a new population of candidate binaries exhibiting dynamical signatures consistent with unresolved companions. This work is part of the GaiaMoons program, and our aim with it was to characterise a sample of 357 potential binary asteroid targets and confirm or refute their binary nature. The properties of these candidates were derived from the high-precision photometric and astrometric observations provided by Gaia. We adopted stellar occultation as the observational method to study these targets. Between October 2023 and February 2026, we successfully carried out 165 observations for 101 targets. We subsequently analysed these events in the context of the available literature and previously reported observations. Thirty three observation led at least two positives for 24 objects that have undergone unprecedented occultation observation campaigns, with four objects showing indications of binary or contact binary features, namely 1127 Mimi, 35420 1998 AG6, 206 Hersilia, and 36882 2000 SW155. For the vast majority of these objects, the resulting dataset from all reduced observations provides unique physical and astrometric constraints, as they had never been observed through stellar occultations before. GaiaMoons illustrates how stellar occultation campaigns associated with Gaia observations generate a self-improving cycle to find new binary, thereby probing size and shape to constrain future observations.

Recent advances in multimessenger astronomy, particularly through gravitational-wave observations of compact-object mergers, have significantly improved our understanding of dense matter. Nevertheless, the internal composition of compact stars remains uncertain. Depending on the underlying equation of state (EoS), these objects may be neutron stars composed primarily of nucleons, quark stars made of deconfined quark matter, or hybrid stars containing both hadronic and quark phases. More exotic constituents, such as hyperons, meson condensates, or dark matter, have also been proposed. In this work, we investigate whether the internal composition of compact stars can be inferred from observable quantities, including mass, radius, and tidal deformability. To address this problem, we employ machine-learning and deep-learning techniques trained on a larg dataset of EoSs describing both neutron stars and quark stars. From these EoSs, we generate the corresponding mass radius relations spanning a wide range of stellar configurations. The resulting dataset is used to train and evaluate classification models aimed at identifying the nature of compact objects from their macroscopic properties. Our results indicate that suitable combinations of observables can distinguish neutron stars from quark stars with very high accuracy. These findings demonstrate the potential of machine-learning approaches as tools for probing the composition of dense matter. However, further studies incorporating additional scenarios, including hybrid stars and other exotic forms of matter, are required to establish the robustness and general applicability of this methodology.

Context. Non-ideal magnetohydrodynamics (MHD) is a key tool for modeling magnetic flux transport in astrophysical systems such as molecular clouds, protostellar cores, and protoplanetary disks. Conventional explicit methods for non-ideal MHD diffusion are severely limited by timestep constraints, while substepping approaches can be unstable due to truncation errors near boundaries and strong magnetic-field gradients. Aims. Our main goal is to address these limitations by developing robust super-time-stepping methods for Ohmic and ambipolar diffusion. Methods. We present a super-time-stepping method based on the stability of the Gegenbauer polynomials. The method is designed to enhance robustness in the presence of strongly anisotropic resistivity and to reduce sensitivity to truncation errors near boundaries. We implement the scheme in the PLUTO code and assess its performance through dedicated Ohmic and ambipolar diffusion tests. We also compare this novel numerical scheme against two common astrophysical problems, namely magnetic reconnection and the magnetorotational instability. Results. The novel Runge-Kutta-Gegenbauer scheme retains computational efficiency beyond purely explicit schemes while providing excellent stability compared with other traditional substepping methods. It remains stable in the presence of strongly anisotropic diffusion, enabling accurate magnetic-field evolution in regimes characteristic of protoplanetary disks and collapsing dense cores. Benchmark tests, including magnetic reconnection and magnetorotational-instability setups, confirm the method's accuracy, efficiency, and suitability for large-scale non-ideal MHD simulations.

In this paper, we describe Japan's possible contributions for coronagraph of the Habitable Worlds Observatory (HWO) based on our independent study. We are considering to contribute to the HOW coronagraph by science and hardware, based on Japan's experience for the SPICA coronagraph instrument, contributions to the Nancy Grace Roman Space Telescope, and SCExAO for the Subaru telescope. Currently, hardware contributions of various scales, from large-scale to small components, are considered. As an example of the large-scale hardware case, the optical and mechanical layout of the entire infrared coronagraph is presented. Several individual high-contrast technologies are also briefly introduced, for which research is ongoing in Japan. In discussions, it is pointed out that both the inner working angle (IWA) and sensitivity are particularly critical for the NIR coronagraph. In this situation, dedicated observations of a small number of targets close to the solar system can be one of key science program in this situation, and designing consolidating science objectives, requirements, observation targets, and survey plans is important. It is essential to push the development of advanced coronagraphs that provide small IWAs. On the other hand, it is also necessary to prepare solutions that adopt more robust coronagraphs in parallel. How to coexist visible and NIR coronagraphs within constraints of volume, mass, budget etc. is an important issue. The international sharing for the coronagraph development should be carefully decided by international agreement. Although all of our studies may not be realized in contributions to the first generation of HWO instruments, we are considering Japan's multigenerational participation in the HWO to maximize outcomes of the HWO.

The detection of Earth analogs with extreme-precision radial velocities (EPRVs) is limited by spectral variability from stellar activity, telluric absorption, and instrumental systematics. We apply AESTRA, a generative spectrum modeling framework, to NEID Sun-as-a-star observations. AESTRA empirically decomposes the spectra into stellar line-shape variability, micro-telluric absorption, and continuum variability without external atmospheric or stellar templates. After removing the learned telluric and continuum components, we train a low-dimensional representation of the spectrum to infer activity-driven apparent RVs jointly with candidate Doppler signals. We evaluate the method with 500 single-planet injection-recovery tests spanning periods of 2.5 to 400 days and semi-amplitudes of K = 0.1 to 0.7 m s^-1, calibrating the detection criterion to yield zero spurious detections. At this matched confidence level, AESTRA recovers 238 injected planets, including 13 with K < 0.3 m s^-1, whereas traditional CCF-based activity-indicator detrending recovers 9 planets and none below K = 0.5 m s^-1.

Many modern near-infrared instruments employ HAWAII-2RG (H2RG) detectors with integration times that can reach 300-600s. Up-the-ramp (UTR) sampling offers advantages over Fowler sampling, including superior cosmic ray rejection and extended dynamic range, but requires fitting linear ramps from 30-60 reads. Ground-based K-band sky brightness has been reported to vary by 3-10% on timescales of minutes, potentially introducing systematic errors and compromising photometric accuracy. Additionally, UTR data formats involve higher-dimensional FITS files with larger file sizes impacting observatory operations. We present a feasibility study using the GIRMOS Data Simulator with high-fidelity flux budgets and empirical K-band sky variations estimated from Gemini-NIRI at 10-20s cadence. Using a Monte Carlo approach we assess whether linear ramp fitting remains viable under variable sky conditions, quantify SNRs and systematic biases, and report nightly data volume estimates. Our results show that the advantages of the UTR readout hold for read-noise-limited targets placed in the inter-line regions, translating into 4-10% savings in observing time. Over the sky emission lines, UTR fitting remains possible but its performance is compromised, both by a degradation in SNR and by a high rate of pixels falsely flagged by the CR rejection algorithm. Both effects are driven by the higher signal level rather than by sky variability and the latter could be mitigated by adapting CR rejection thresholds to the local signal level. These findings address how ground-based conditions affect UTR implementation in near-infrared spectrographs, with GIRMOS as a concrete case study.

The detection of high-frequency gravitational waves (HFGWs) above 10 kHz provides a crucial probe of exotic astrophysical phenomena and new physics. We report the first search for HFGWs via their conversion to electromagnetic radiation through the inverse Gertsenshtein effect in Earth's magnetic field, utilizing radio telescopes including the Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA). Since no statistically significant signal is observed, we obtain new upper limits on the characteristic strain across the 1 GHz – 1 THz band, with the most stringent constraint reaching $h_c \lesssim 10^{-18}$, improving upon existing bounds by up to three orders of magnitude. These results significantly advance the exploration of uncharted parameter space for exotic gravitational-wave sources, paving the way for future discoveries with next-generation facilities such as the Square Kilometre Array (SKA).

The era of large-scale, high-cadence astronomical surveys demands efficient and robust methods for time-series analysis. ARIMA models provide a versatile parametric description of stochastic variability in this context. However, their practical use is limited by the challenge of selecting optimal model orders while avoiding overfitting. We present a novel solution this problem by combining Autoregressive Integrated Moving Average (ARIMA) models with the Nested Sampling algorithm. Our method yields Bayesian evidences for model comparison and also incorporates an intrinsic Occam's penalty for unnecessary model complexity. Using JAX and Blackjax, a vectorized ARIMA-Nested Sampling framework with GPU-acceleration support is implemented, allowing us to perform model selection across grids of Autoregressive (AR) and Moving Average (MA) orders, with efficient inference of selected model parameters. We validate the approach using simulated time series with known ground-truth parameters and demonstrate accurate recovery of both model order and parameters. We then apply the method to several astronomical datasets, including the historical sunspot number record, stellar light curves of KIC 12008916 and Kepler 17 from the Kepler mission, and quasar light curves of 3C 273 and S4 0954+65 from the TESS mission. For all cases, except Kepler 17, the ARIMA models selected by this method were able to accurately model the stochastic variability in the time series data as well as produce accurate multi-step ahead forecasts for the sunspot number time series. Our results demonstrate that nested sampling offers a rigorous and computationally tractable alternative to autoregressive model selection in astronomical time-series analysis.

Amidst the exoplanet revolution in which multiple techniques have successfully found planets, the Doppler (Radial Velocity, or "RV") technique is unique in its sensitivity to giant planets at long orbital periods around Sun-like stars. The upcoming retirement of Keck-HIRES will incur irreversible changes in the continuation of HIRES's decades-long stable RV baseline and with it, the exoplanet community's ability to detect giant exoplanets with periods longer than Jupiter. With the time elapsed from the last HIRES RV for many stars of interest at  3 years and growing, we tested the impact of a "critical RV", one that would bridge this gap between past HIRES RVs and future stable Keck-KPF RVs, on the recovery of long-period giant exoplanets. We generated 2000 1-planet systems with RVs sampled at a representative timeseries and used the planet-finding code Octofitter to perform injection-recovery experiments including and omitting this critical RV for each system. For the injected long-period super-Jupiters ( 8-55 years, 1-13 $M_J$), including the critical RV induced a $1.5\times$ enhancement in overall planet recovery and a more specific $3.5\times$ enhancement in the recovery of super-Jupiters with Saturn-like periods. These experiments show that gathering a critical RV for stars of interest can help ensure that HIRES's decades-long stable RV baseline in conjunction with future KPF RVs, or indeed that the observationally-gapped RV baselines of any instruments that will undergo an RV zeropoint offset, will continue to be foundational to the discovery of long-period giant exoplanets in years to come.

Large foundation models (FMs) are transforming Earth science by integrating heterogeneous multimodal data, such as multi-platform imagery, gridded reanalysis data, diverse geophysical and geochemical observations, and domain-specific text, to support tasks ranging from basic perception to advanced scientific discovery. This paper provides a unified review of Earth science foundation models (Earth FMs) through two complementary dimensions: depth, which traces the evolution of model capabilities from perception to multimodal reasoning and agentic scientific workflows, and breadth, which summarizes their expanding applications across the atmosphere, hydrosphere, lithosphere, biosphere, anthroposphere, and cryosphere, as well as coupled Earth system processes. Using this framework, we review representative multimodal Earth foundation models and compile more than 200 datasets and benchmarks spanning diverse Earth science tasks and modalities. We further discuss key challenges in multimodal data heterogeneity, scientific reliability and continual updating, scalability and sustainability, and the transition from foundation models to agentic and embodied Earth intelligence, and outline future directions toward more integrated, trustworthy, and actionable AI Earth scientists. Overall, this paper offers a structured roadmap for understanding the development of Earth foundation models from both capability depth and application breadth.

With the next generation of both electromagnetic and gravitational wave observatories beginning to come online, rapid analysis methods for kilonova data are becoming increasingly important in astronomy. Traditional Bayesian parameter estimation using Markov chain Monte Carlo (MCMC) is time-consuming and relies on explicit likelihood approximations that can break down when modeling uncertainties are significant. We develop a simulation-based inference (SBI) framework for kilonova parameter estimation using density-estimation likelihood-free inference. The framework uses a Gaussian process emulator trained on $\sim 1300$ POSSIS simulations. We demonstrate that SBI provides a rapid alternative to MCMC that is robust to likelihood misspecification. The standard Gaussian likelihood approximation fails to capture the non-Gaussian, correlated structure of emulator uncertainty; SBI learns this structure directly from forward simulations. Simulation studies show that the SBI method accurately recovers injected parameters, while the MCMC suffers from systematic bias caused by likelihood misspecification. This problem persists when analyzing AT2017gfo, where a subset of the MCMC posteriors pile up at prior boundaries and the SBI posteriors do not. The SBI framework infers a total ejecta mass of $\sim 0.087 M_{\odot}$ dominated by lanthanide-poor ejecta and excludes toroidal and peanut ejecta geometries at the 99th percentile for both components. The SBI framework generates $\sim 2 \times 10^{4}$ posterior samples in seconds.

I present EXOVEIL, a transit detection system that learns what a star's brightness should look like and flags when reality disagrees. Unlike existing systems that require phase-folded input, EXOVEIL operates on raw flux time series and can detect planets that transit only once.A Transformer world model, trained on 16,499 Kepler light curves with transit-masked self-supervised learning, predicts expected stellar flux. A matched-filter detector with variance weighting extracts transit signals from the prediction residuals. A learned classifier (XGBoost) separates planets from false positives, achieving AUC 0.938 on Kepler DR25. Applied to single-transit injection-recovery, EXOVEIL recovers 32% of transits at 1000 ppm depth a task where all classification-based systems score 0% by construction. A blind search of 3,737 Kepler stars yields 179 new transit-like signals not present in the DR25 TCE catalogue, including 46 monotransit candidates. Applied withoutretraining to 47 confirmed TESS planets in the PLATO LOPS2 field, EXOVEIL achieves 100% recovery, demonstrating zero-shot cross-mission transfer. At PLATO's 25-second cadence, detection reaches 100 ppm – approaching the Earth-analog regime. I provide the first application of conformal prediction to transit detection (95.9% empirical coverage) and release the system as pip install exoveil with pretrained weights and a candidate catalogue.

This paper presents a novel nonlinear backstepping control law for continuous, low-thrust station-keeping in the Earth-Moon system. Quasi-periodic libration point orbits are targeted under a high-fidelity model of the dynamics. Almost global uniform exponential stability guarantees are attained, as shown through Lyapunov's stability theory. Saturation of the actuators is formally included in the controller design, such that these guarantees hold even in the event of saturation. The relationship between saturation threshold, control gains, and deviation is studied and an optimal procedure for gain selection is discussed. The control solution is tested numerically through a Monte Carlo analysis over representative application cases, subject to operational errors, constraints, and external perturbations. Station-keeping under actuation saturation is validated considering a conservative threshold for typical electric propulsion systems.