Priority

We present a suite of direct N-body simulations of low mass ($<10^4~M_{\odot}$) globular cluster streams initialized with observationally-motivated binary demographics in order to understand the effect of in-cluster dynamical processing on the stream binary population. The models are initialized with a range of stellar densities and cluster orbits, and Poisson variation in the number of massive and short-lived stars. Wide binaries are disrupted on short timescales by internal tides and on long timescales by two-body encounters. Tides are most important prior to impulsive mass loss-driven cluster expansion. Close binaries ($P_{\rm orb}<10^2~\rm yr$) are most abundant at the stream center due to cluster mass segregation. The wide binary fraction and the degree of binary segregation in the resulting stream are sensitive to the initial cluster density and massive star fraction. In mock radial velocity surveys of the simulated streams, undetectable binaries have velocity amplitudes of $\sim$$0.5$-$1~\rm km~s^{-1}$, adding $\sim0.1~\rm km\ s^{-1}$ of velocity dispersion to the streams, and are dynamically depleted by $\sim10$-$60\%$ compared to the initial binary population. Custom N-body models of Milky Way streams with binaries will allow a holistic understanding of their dynamical structures in advance of upcoming multi-epoch spectroscopic surveys.

Molecular spectral cubes of prestellar cores encode the information on the physical and chemical properties of these objects along the line of sight. To retrieve this information, we need an interpretable model that reproduces the observed spectra. We designed a differentiable 3D geometrical model that produces synthetic observations from the parameterized density and velocity fields, and that can be efficiently optimized to reproduce the real data cubes. The model has been applied to p-NH2D and N2D+ spectral cubes in the prestellar core L1544. The optimized model suggests that to reproduce the observed velocity difference between p-NH2D and N2D+ in L1544, an asymmetric structure in density and velocity is necessary.

Galaxy clusters grow through the accretion of galaxies from groups, filaments, and other clusters. During this process, galaxies may undergo pre-processing in lower-density environments, where galaxy-galaxy mergers and other interactions can significantly alter their properties prior to cluster infall. We investigate the role of galaxy mergers in the pre-processing of galaxies prior to cluster infall by studying the spatial distribution of mergers across the cosmic web. We use a sample of 43,922 galaxies targeted by the 4MOST CHANCES survey in and around 33 low-redshift clusters (z < 0.07). Using Zoobot, a deep-learning framework trained on Galaxy Zoo data, we identify 698 galaxy mergers. We measure their distances to cosmic web filaments and compare them with those of non-merging galaxies. We find that galaxy mergers are significantly closer to filaments than the non-merging galaxy population, with this trend being strongest beyond the cluster virial radius. This suggests that filaments provide conditions conducive to mergers, possibly moderating relative velocities and enhancing gas availability. Our findings support a scenario in which filaments play a key role in transforming galaxies through pre-processing by promoting mergers before they enter cluster cores where star formation quenches.

The relationship between magnetic field strength and gas density is essential to understand the interstellar medium and star formation. Zeeman measurements in dense atomic and molecular gas phases have traditionally been used to directly probe magnetic field strengths in the Milky Way. This allowed derivation of a relationship between magnetic field strength $B$ and gas number density $n$. We recently generalized this relation as a two-part power-law with non-zero slopes and a transition density given as $B/B_0 \propto (n/n_0)^{\alpha_1}$ for $n \le n_0$ and $(n/n_0)^{\alpha_2}$ for $n > n_0$. Here, we extend our previous hierarchical Bayesian framework by incorporating a large body of pulsar observations that probe the diffuse interstellar medium and explicitly modelling density uncertainties through a global log-density correction parameter $R$ applied to all densities. We also account for magnetic field geometry and measurement uncertainties through a magnetic hyperparameter to estimate $B$. This results in a stronger constraint on the diffuse gas part of the $B$–$n$ relation. Our results confirm a non-zero exponent in the diffuse gas and a broad transition density with our best model and data set yielding maximum a posteriori results of $\alpha_1 = 0.18^{+0.02}_{-0.02}$, $\alpha_2 = 0.63^{+0.08}_{-0.05}$, $n_0 = 1630^{+2560}_{-1430}\,\text{cm}^{-3}$, and $B_0 = 7.60^{+2.00}_{-3.47}\,\mu\text{G}$.

We present the results from a 3-year-long Ultraviolet monitoring campaign of the broad line Seyfert 1 galaxy ESO 141-G55 using International Ultraviolet Explorer (IUE). By modelling all individual, extinction-corrected UV spectra in 1150-1978 A and 1850-3348 A wavelength range, we have observed a significant variability in both UV continuum and line fluxes. Variabilities due to ionised UV lines like SiIV, CIV and HeII are delayed with respect to the UV continuum by 2.92$^{+0.54}_{-0.61}$, 4.41$^{+0.44}_{-0.54}$, 4.11$^{+0.35}_{-0.81}$ days, respectively. At a distance of $\sim$0.004c, an outer accretion disc can be a possible site for the origin of UV lines.

Interstellar objects provide the only directly observable samples of icy planetesimals formed around other stars, and can therefore provide insight into the diversity of physical and chemical conditions occurring during exoplanet formation. Here we report isotopic measurements of the interstellar comet 3I/ATLAS, which reveal an elemental composition unlike any Solar System body. The water in 3I/ATLAS is enriched in deuterium, at a level of D/H = (0.95 +- 0.06)%, which is more than an order of magnitude higher than in known comets, while its range of 12C/13C ratios (141-191 for CO2 and 123-172 for CO) exceeds typical values found in the Solar System, as well as nearby interstellar clouds and protoplanetary disks. Such extreme isotopic signatures indicate formation at temperatures $\lesssim30$ K in a relatively metal-poor environment, early in the history of the Galaxy. When interpreted with respect to models for Galactic chemical evolution, the carbon isotopic composition implies that 3I/ATLAS accreted roughly 10-12 billion years ago, following an early period of intense star formation. 3I/ATLAS thus represents a preserved fragment of an ancient planetary system, and provides direct evidence for active ice chemistry and volatile-rich planetesimal formation in the young Milky Way.

Halos surrounding spiral galaxies act as the bridges connecting the galactic disk and the intergalactic medium (IGM). They host a significant fraction of the baryonic mass in the Universe, and feedback from star formation (SF) or active galactic nuclei (AGN) likely plays an important role in regulating this vertical baryonic component. Despite its importance, the contribution of extraplanar molecular gas remains poorly understood. We aim to characterize the vertical extent and the kinematics of molecular gas traced by CO(2-1) emission in the nearby (D = 9.5 Mpc) spiral galaxy NGC 891, one of the best studied edge-on galaxies. We also compare our results with HI, H$\alpha$-traced DIG and dust maps from the literature. Our analysis is based on new CO(2-1) observations of NGC 891 obtained with the IRAM 30m telescope. We mapped two 6 kpc $\times$ 6 kpc regions on the northeastern side and the area surrounding the galactic center. We apply a careful method to estimate and remove the residual contribution of the error beam to the CO cube. The vertical extent of the molecular gas is best described by a two-component Gaussian fit, consisting of a bright thin disk component (deconvolved FWHM $\simeq$ 360 pc) and a fainter thick disk component (deconvolved FWHM $\simeq$ 1.1 kpc). Statistically significant CO(2-1) emission is detected up to 1.3-1.4 kpc above the disk midplane. We estimate that the thick molecular disk component contains up to 27% of the total molecular gas mass of the galaxy. Our results demonstrate that SF-driven feedback in a non-starburst galaxy can lift significant amounts of molecular gas to large vertical distances. We interpret the presence of extraplanar molecular gas in NGC 891 in the framework of a galactic fountain scenario, in which material is expelled from star-forming regions and transported toward the outer halo.

Different models of filament formation predict distinct patterns of angular momentum redistribution toward embedded cores, set by the underlying velocity-field structure, which can set the initial conditions for a preferential orientation between protostellar outflows and filaments. However, the absence of a dominant alignment in observations keeps this connection open to debate. We investigate whether gravity-driven longitudinal flows along filaments can redistribute angular momentum (AM) toward collapse centers and influence outflow-filament alignment. To this end, we analyze the distributions of 3D and 2D-projected angles between sink angular momentum vectors and host filament orientations in an SPH simulation of giant molecular cloud and filament formation. We also characterize the filament velocity field by measuring the angles between SPH particle velocity vectors and filament axes, and the degree of convergent flow toward filament density peaks. No preferred alignment between the sinks' AM and the filament direction is found at early evolutionary stages, neither in 3D nor in 2D. Later, however, a predominantly perpendicular configuration emerges in 3D. Tracking individual sinks indicates that this alignment is not primordial but develops as gravity strengthens. In individual filaments, the onset of perpendicular alignment coincides with the development of convergent longitudinal flows. Finally, we estimate the minimum fraction of perpendicular 3D angles required to reveal a perpendicular 2D alignment for a given sample size. While longitudinal flows develop over extended timescales, once established, they can rapidly reorient the angular momentum vector of the sinks, enabling perpendicular alignments to arise within typical outflow lifetimes.

All detected water reservoirs in the solar system exhibit a deuterium enrichment that links back to the physical environment at the time of stellar birth. Gas-phase and ice-grain D-enrichments occur through chemical processes that operate at low temperatures ($<$ 30 K) pointing towards an origin in the pre-stellar molecular cloud or in the outer parts of the protoplanetary disk. However, not all stars are born in environments similar to our Sun, nor do their subsequent evolutionary histories follow the same path. These environmental differences can be traced by the water D/H ratio. Here we use ALMA observations of the interstellar comet 3I/ATLAS to constrain the water D/H ratio in extrasolar cometary material. With a water D/H value of [D/H]$_{\mathrm{H_2O}} > 6.6\times10^{-3}$, 3I/ATLAS shows a deuterium enrichment exceeding Earth's ocean value by more than a factor of $\gtrsim40$ and typical Solar System cometary values by more than a factor of $\gtrsim30$. The elevated deuterium enrichment points to water that formed under colder, less irradiated conditions and from less thermally processed material, consistent with an origin in a planetary system that formed under different physical and chemical conditions than our own.

We present a new calibration of the broad-line region (BLR) size-luminosity-velocity relation using a sample of 157 AGNs with reliable Hbeta time-delay (\lag) measurements from Wang & Woo 2024. By incorporating the Eddington ratio as a third parameter, we effectively correct the systematic offset of high-Eddington AGNs in the traditional BLR size-luminosity relation. The resulting three-parameter fit defines a fundamental plane in the 3-D space of the \lag, optical luminosity, and Hbeta velocity, with an intrinsic scatter of 0.21 dex. This tight correlation reflects the coupled effects of gas kinematics, photoionization, and BLR geometry. In turn, we develop a new method to infer \lag\ from the combination of optical luminosity and Hbeta velocity, and derive single-epoch black hole mass estimators by adopting either the full-width-at-half-maximum (FWHM) or line dispersion ($\sigma$) of the Hbeta line profile as the velocity indicator. The derived \lag shows a  0.1 dex scatter, depending on the choice of calibrations. We show that the previous mass estimates based on the two-parameter size-luminosity relation with a 0.5 slope can be overestimated by up to 0.5 dex, demonstrating that the new mass estimator substantially changes the cosmic black hole mass density and the growth of black hole seeds in the early universe.

We present 4.7–22.2 GHz Very Long Baseline Interferometry (VLBI) monitoring of the candidate pre-coalescence supermassive black hole binary SDSS J143016.05+230344.4 ($z=0.08105$) from 2022 February to 2024 February, together with quasi-simultaneous 0.7–16.5 GHz connected-array spectra. At all epochs, the radio emission is dominated by a single unresolved milliarcsecond core with $T_{\rm B}\gtrsim10^{7}$ K, confining the variable emission to $\lesssim0.3$ pc. The spectra require two self-absorbed synchrotron components: a persistent low-frequency component with $\nu_{\rm p,steady}\approx0.74$ GHz and $S_{\rm p,steady}\approx1.22$ mJy, and a flare component whose turnover evolves from $(6.35\,{\rm GHz},0.18\,{\rm mJy})$ in 2022 February–May to $(8.61\,{\rm GHz},0.38\,{\rm mJy})$ in 2022 December and then to $(5.83\,{\rm GHz},0.25\,{\rm mJy})$ in 2023 March–April. The 15 GHz flare fraction peaks at $\simeq80\%$ and matches the near-epoch VLBI recovery fraction, showing that the high-frequency brightening arises from a new compact synchrotron component. A second 15.2 GHz VLBI-core brightening is detected from 2023 September to 2024 February while the source remains unresolved. Equipartition scalings imply characteristic radii of $R_{\rm eq}\sim5\times10^{-4}$ pc for the flare and $\sim9\times10^{-3}$ pc for the steady component, and a steep inner circumnuclear density profile, $n\propto R^{-1.7}$. The delayed radio peak is consistent with dissipation of an outflow or jet-base disturbance in a structured circumnuclear medium, while a uniform free–free absorber is disfavored.

We investigate how diversity in the core-halo mass relation affects constraints on the fuzzy dark matter particle mass $m_\psi$ inferred from the internal kinematics of dwarf galaxies. Using stellar line-of-sight velocities and projected positions for eight Milky Way dwarf spheroidal galaxies, we model their dark matter halos as solitonic cores embedded within outer NFW envelopes. We apply both second- and fourth-order Jeans analyses to derive the posterior distribution of $m_\psi$. Our results show that there are two ranges of $m_\psi$ consistent with the observed kinematics: $-20.3 < \log_{10}(m_\psi/\mathrm{eV}) < -19.2$, and a narrower small-mass window $-22.1 < \log_{10}(m_\psi/\mathrm{eV}) < -21.5$, both within the 68\% credible intervals. The latter becomes prominent only if core-halo diversity is taken into account. These constraints pose a challenge to fuzzy dark matter, as the small-$m_\psi$ window is in conflict with Milky Way satellite abundances, and our upper bound largely excludes the parameter space permitted by Lyman-$\alpha$ forest constraints.

Interstellar objects provide a unique opportunity to further our understanding of the planetary formation process by studying in detail material formed around another star. Their ices contain precious clues about the environment and conditions prevailing in their home system. As fractionation processes can be sensitive to the temperature and radiation environment, isotopic ratios are powerful tracers of the origin and evolution of different species. While isotopic ratios have been measured in solar system comets, previously detected interstellar objects have been too faint to measure isotopic ratios. Here we report the measurement of two ratios in 3I/ATLAS from observations of the CN molecule: $^{12}$C/$^{13}$C and $^{14}$N/$^{15}$N. We report $^{12}$C/$^{13}$C=$147^{+87}_{-40}$ and $^{14}$N/$^{15}$N=$343^{+454}_{-124}$. The $^{14}$N/$^{15}$N is higher than the value of $\sim$ 150 usually measured for solar system comets, close to the values measured in the interstellar medium, pre-stellar phases or the outside of protoplanetary discs. The $^{12}$C/$^{13}$C is marginally higher than the values usually measured for solar system comets and in the interstellar medium. These measurements could indicate an origin of 3I in the outer disc around an older low-metallicity star.

The environments around star clusters evolve as stellar feedback reshapes the interstellar medium and dynamical processes reorganize the structure of the surrounding stellar field. As approximately single-age populations, star clusters can serve as clocks to trace these environmental changes. In this exploratory study, we test whether convolutional neural networks (CNNs) can identify age-dependent changes in cluster environments. We take cluster ages as given from basic SED fitting of five-band UV-optical aperture photometry from the PHANGS (Physics at High Angular resolution in Nearby GalaxieS) HST survey. We first show that CNNs can be trained on image cutouts centered on clusters to recover ages directly from imaging. This demonstration provides the foundation for this study, which examines whether the information used by CNNs to predict age is coherent and physically meaningful. We perform controlled image occlusion experiments as an explainable AI method. These show that the CNNs extract age-predictive environmental cues in the absence of cluster light and when information on SED shape is removed by combining the five filters into one image. We find that reliance on environmental information increases at the youngest (<10 Myr) and oldest (>1 Gyr) ages, where clusters can exhibit similarly red colors. Our results are consistent with the long-recognized picture that cluster environments evolve systematically with age. We demonstrate that this information is encoded at a level detectable by machine-learning and recoverable from broadband imaging. This establishes a path for using new techniques to connect image-based age inference to the physical evolution of cluster environments.

We present multi-epoch optical spectroscopy of the interstellar comet 3I/ATLAS obtained from December 2025 to January 2026 (heliocentric distances 1.8 to 3.3 au). From these spectra, we derive post-perihelion production rates and/or mixing ratios for daughter molecules (CN, C$_3$, C$_2$, and CH) and gaseous metals (Fe I and Ni I). We also estimate a lower limit on the CO abundance based on the [O I] $\lambda6300$ line. The resulting outgassing profiles reveal a pronounced perihelion asymmetry, with volatile production declining more gradually outbound than during the inbound phase. In addition, the coma becomes less depleted in C$_2$ following perihelion, and shows a substantial increase, relative to H$_2$O, in the production rates of metals and potentially also CO. These trends may indicate the activation of subsurface material or compositional heterogeneity revealed by seasonal effects. The inferred high CO abundance further suggests possible decoupling of CO from H$_2$O and CO$_2$, similar to the behavior observed in solar system comets. The potential post-perihelion enhancement of both CO and metal production, if confirmed, would also be consistent with metal carbonyls contributing to the release of gaseous metals in cometary comae.

We present a sample of large double radio sources hosted by spiral galaxies (Spiral Double Radio Active Galactic Nuclei, SDRAGNs). Candidates were selected during Radio Galaxy Zoo, and refined using Sloan Digital Sky Survey images. The most promising were targeted in the Zoo Gems Hubble Space Telescope program, yielding images for 36 candidates. We assess the likelihood of each spiral galaxy being the genuine host of the radio emission finding 15 new high-probability SDRAGNs. SDRAGN hosts are seen preferentially close to edge-on. SDRAGNs predominantly show FR II radio structures and optical pseudobulges. Accounting for sample selection effects, the radio-jet axes lie preferentially near the poles of the galaxy disks; we find a constant probability distribution for intrinsic pole-jet angles < 30 degrees, ramping to zero at 60 degrees. We have obtained optical spectra for all these new SDRAGNs. Among both previous and new SDRAGN samples, 8/25 show Seyfert 2 signatures, 6/25 show central star formation, and 5/25 show LINER emission strong enough to indicate AGN or shock ionization, broadly similar to radio galaxies in elliptical hosts with the addition of star formation (diluting or masking weak AGN signatures). SDRAGNs include FR II sources seen at unusually low radio power, and preferentially occur in significant galaxy overdensities on 1-Mpc scales. Our "false alarms" - systems where HST data show the spiral to not be the actual host galaxy - include radio sources seen through large parts of foreground spiral disks, potentially useful for Faraday-rotation studies of disk magnetic fields.

The Galactic Bulge, as the center of the Galaxy, is the closest laboratory for studying galaxy formation and evolution. However, its study faces significant challenges due to heavy dust extinction. This paper is devoted to deriving the average age of the Galactic Bulge and investigating its spatial distribution. We utilize a high-precision PSF-fitting photometric catalogue in the $J$ and $K_{\mathrm{s}}$ bands observed by VISTA to study the average stellar ages within the Bulge. Red giant stars are employed as tracers, with their average distances determined using red clump stars as references. The average ages are fitted with stellar models. Our analysis reveals a systematic age gradient across the Galactic Bulge ($2^{\circ} < |b| < 8^{\circ}$). The mean stellar age increases significantly with galactic latitude, shifting from a younger population ($\sim 4.69^{+0.97}_{-0.81}$ Gyr) prevalent near the plane to a predominantly older population ($\sim 10.48^{+0.93}_{-0.85}$ Gyr) at higher latitudes. We hypothesize that the young stellar population at low latitudes is predominantly composed of a pseudo-bulge formed via disk/bar processes (incorporating contributions from recent star-forming activity in the Galactic center), whereas the older stellar population is associated with spheroidal bulges generated through early-stage collapse or accretion of debris from merged dwarf galaxies.

Low-surface-brightness (LSB) structures provide critical insights into the hierarchical formation of galaxies and galaxy clusters. The KASI Deep Rolling Imaging Fast Telescope (K-DRIFT) is designed to detect such diffuse features through deep, wide-field optical imaging with a surface brightness reaching $\sim$$30~\rm{mag}~\rm{arcsec}^{-2}$. To interpret the observation data expected from K-DRIFT, we have developed the Galaxy Replacement Technique (GRT), an $N$-body simulation framework optimized for tracing the gravitational evolution of stellar components. The GRT works by inserting high-resolution galaxy models, including a dark matter (DM) halo and stellar disk, in place of multiple low-resolution DM halos in the base $N$-body cosmological simulation. It allows us to achieve very high mass ($m_{star}=5.4\times10^4\msun\ h^{-1}$) and spatial resolution (10 $\rm{pc}~h^{-1}$) with shorter computation time compared to full hydrodynamic cosmological simulations. Therefore, this technique is particularly well-suited for studying LSB structures, with a surface brightness reaching $\sim$$31~\rm{mag}~\rm{arcsec}^{-2}$. In this paper, we present the motivation and methodology of the GRT, summarize key results from previous studies, and highlight its synergy with K-DRIFT observations. We further discuss planned science cases using the GRT, aiming to build a theoretical basis for interpreting LSB features in various environments.

We examine to what extent the galaxy spin parameter proxy ($\lambda_R$) is affected by bar formation and how it is related to the strong and weak classifications of stellar bars. By creating mock observations of a simulated galaxy, we show that the emergence of a stellar bar can cause mass-weighted $\lambda_R$ to decrease by up to 16%, depending on the bar's orientation. This decrease can be exaggerated if there is a burst of star formation due to the bar driving gas to the center of the galaxy. We use the SAMI galaxy survey to show that weakly barred galaxies have statistically significant younger average stellar populations, higher galaxy spin proxy and higher specific star formation rates compared to strongly barred galaxies within one effective radius. If we consider galaxies with average light-weighted stellar population age less than 3 Gyr within one effective radius, we still find weakly barred galaxies to have a higher galaxy spin proxy than strongly barred galaxies. Based on these trends found from the SAMI galaxy survey, we suggest weakly barred galaxies are rapidly forming, similar to the bar formation process seen in simulations, while strongly barred galaxies are undergoing slower (secular) evolution.

Red supergiants (RSGs) are crucial for studying the properties and evolution of massive stars. It is representative to conduct a census of RSGs across the Local Group, which spans a broad metallicity range. However, identifying RSGs in distant and metal-poor galaxies remains challenging mainly due to contamination of foreground dwarfs and observational limitations. In this work, we perform PSF photometry on publicly released JWST/NIRCam images of five Local Group galaxies: NGC 6822, Sextans A, NGC 300, WLM, and IC 1613 using the DOLPHOT NIRCam module. We find an optimal color-color diagram (CCD) for metal-poor environments, that is F115W $-$ F200W versus F356W $-$ F444W, which clearly separates RSGs from foreground dwarfs. By using the CCD, we identify 208, 135, and 22 RSG candidates in NGC 6822, Sextans A, and NGC 300, respectively, free from contamination by foreground dwarfs and oxygen-rich asymptotic giant branch stars (O-AGBs). In addition, 40 and 14 RSG candidates are directly selected on the CMD in WLM and IC 1613, respectively. Compared with previous works, the number of RSG candidates within the same luminosity range and sky region increases significantly, demonstrating the advantages of JWST in constructing a more complete RSG sample in the Local Group thanks to its high spatial resolution and photometric quality. In addition, catalogs of O-AGBs and carbon-rich AGBs (C-AGBs) are provided as by-products.

We report the discovery of CGG-z7, the most compact galaxy group at $z\gtrsim7$ identified to the north of the GOODS-North field, observed by the JWST POPPIES program. The system consists of at least six members within a projected size of $7.8\times5.7$ kpc$^2$, four of which are spectroscopically confirmed via [O III] and H$\beta$ emission. The group exhibits a low line-of-sight velocity dispersion ($\approx93.7$ km s$^{-1}$) relative to its substantial stellar mass ($M_* \approx 10^{9.8} M_{\odot}$), yielding a stellar-to-dynamical mass ratio of $M_*/M_{\mathrm{vir}} \approx 0.15$. This ratio, exceeding typical values for virialized halos by a factor of $3$, indicates that the system is highly likely not in dynamical equilibrium. Instead, we interpret CGG-z7 as a pre-virialized structure, likely a major merger caught near apocenter – capturing the rapid, chaotic formation of a massive "Red Nugget". Spectroscopic analysis reveals extreme ionization conditions and low metallicity across the group. In particular, the central galaxy reaches an extraordinarily high [O III]/H$\beta$ ratio of $\sim18$, which is likely indicative of an obscured AGN. CGG-z7 thus serves as a unique laboratory for the physics of pre-virialized galaxy assembly, bridging the gap between turbulent high-$z$ assembly and the quiescent galaxies seen at cosmic noon.

The embedded disks surrounding individual Class 0 protostars are structureless. Disks surrounding Class I stars may be continuous or have a ring-gap substructure, whereas all disks around Class II stars have a ring-gap substructure that gradually disappear as the disks evolve into debris disks. This common sequence in young lone stars requires an explanation. This study aims to show that the physical model Reaction-Diffusion Systems with Moving Reaction Front can be used to describe and classify protostellar disks according to their structure. A comprehensive review of observations made with the ALMA radio telescope shows: first, that the protostar-disk system presents a geometry analogous to that of an reaction-diffusion system with two separate compartments, namely, protostar and disk. Second, that in the protostar, matter is processed at high temperature, resulting in a chemical composition different from that of the disk. Third, that the equatorial outflow emitted by the protostar, rich in highly reactive trihydrogen cation, acts as a moving reaction front, MRF, that triggers the formation of molecules and nuclei in the disk. The time lag of nucleation with respect to the passage of the MRf would be the cause of the formation of the gaps between the rings of particles that form in the disk. The MRF is a transient phenomenon and its passage causes the transformation of a continuous disk, Class 0, into a disk with a ring-gap structure, Class II, whose temporal evolution begins at the interface of the star and moves outwards.

Globular clusters (GCs) exhibit varying degrees of flattening (ellipticity), which may provide insight into their internal dynamics and evolution histories. Commonly used methods to measure ellipticity, such as ellipse fitting of density contours and principal component analysis, often produce biased results, especially for clusters that are nearly round or have few observable stars. Using a combination of ground-based and space-based photometry, we investigate the shapes of 29 Galactic GCs. To that end, we test two commonly used methods: an ellipse fit to a kernel density profile and a principal component analysis. We find that both methods suffer from bias arising when the number of stars is small or the cluster is close to round. To solve this issue, we develop a robust method to measure the ellipticity of GCs, test it extensively on mock data, and apply it to the 29 Milky Way GCs in our sample. Using the $V/\sigma$ diagram used in the isotropic oblate rotator framework, we examine potential causes for the flattening, including rotation and velocity anisotropy. For ten clusters: NGC 104, NGC 1261, NGC 2808, NGC 3201, NGC 5286, NGC 5904, NGC 5986, NGC 6205, NGC 6341, and NGC 7078 we identify a very good agreement between the rotation angle and semi-minor axis of the ellipse, further corroborating the findings that rotation is the main driver of the ellipticity. The $V/\sigma$ diagram reveals that velocity anisotropy or tides could also be important in shaping the GCs. The robust method developed provides reliable measurements of the ellipticity of GCs, emphasising the importance of taking into account the flattening in theoretical models and simulations. It also offers a promising way to investigate the shapes of multiple stellar populations within GCs, where only small samples are usually available.

In recent years, massive star cluster environments have proved to be bright sources of very-high energy gamma-rays, in particular young clusters which are powered by the winds interacting in their cores. In order to understand how these winds can accelerate particles up to very-high energies, it is necessary to model their interactions from small (sub-pc) to large (10s of pc) scales over several millions of years. A key open question concerns the structure and properties of the resulting wind termination shock. By performing 3D magnetohydrodynamic simulations of clustered winds embedded in a superbubble cavity, we demonstrate that the dynamics of stellar wind interactions and the resulting shock structure solely depends on the density and pressure of the cavity. This implies that the initial conditions of the simulation can be tuned in order to simulate star clusters of arbitrary age at a reduced computational cost. This novel method is validated using a toy cluster hosting 30 identical stars. We discuss the properties of the resulting cluster-wind termination shock under various assumptions. In particular, we are able for the first time to obtain a fully decoupled spherical wind termination shock for a 5 Myr old cluster. We further show that radiative cooling increases the sphericity of the shock. In general, the morphology of the outflow depends on the number of dominant stars, on the power of the stars sitting at the edge of the cluster core, and on the compactness of the cluster. We additionally show how a semi-analytical model can be used in order to estimate key morphological properties of the outflow without relying on large-scale simulations.

We introduce a new, discretized grain size evolution model, incorporated into the GIZMO code and coupled with FIRE-3 stellar feedback and ISM physics, to investigate variations in dust abundance, chemical composition, and grain sizes observed in the Local Group. This model tracks the size evolution of specific dust species, and includes stellar production of dust, dust growth through gas-phase metal accretion, dust destruction by sputtering, SNe shocks, and astration, grain-grain collisional shattering and coagulation, and turbulent dust diffusion. Using idealized galaxy simulations, we test the dependence of MW dust properties on variations in each dust process and find that our model uniquely predicts a bimodal grain size distribution. This bimodality is due to our simulation's ability to resolve each dust process and where they occur in the ISM, unlike other works. We find that Local Group dust abundances are determined by dust growth and destruction, with little dependence on coagulation or shattering, explaining why models that do not include these processes can match abundance observations. We also find that variations in Local Group extinction curve slopes are determined by coagulation, with inefficient coagulation leading to steeper slopes. However, inefficient coagulation also results in stronger extinction curve bumps, which are not observed. We also do not predict a population of very small (${<}1$ nm) carbonaceous grains, required for MIR emission features, due to their rapid growth by accretion. These results highlight the possible necessity of “top-down” PAH formation from preexisting grains as a means to inhibit carbonaceous dust growth.

Giant Molecular Clouds (GMCs) are dominated by supersonic turbulence, creating a complex network of shocks and filaments that regulate star formation. While the global inefficiency of star formation is well-observed, predicting exactly which gas parcels within a turbulent cloud will collapse to form stars remains a challenge. In this work, we present a supervised machine learning framework to forecast the Lagrangian history of prestellar cores in magnetohydrodynamic (MHD) turbulence. We utilize Extreme Gradient Boosting (XGBoost) to train a regression model on the trajectories of $\sim 2.1$ million tracer particles evolved within a self-gravitating, turbulent MHD simulation. By mapping the instantaneous phase-space state (position, velocity, and density) of gas parcels to their future coordinates, our model successfully predicts the 3D evolution of star-forming cores over a horizon of $\sim 0.45$ Myr ($0.25~t_{\rm ff}$). We achieve a global coefficient of determination of $R^2 > 0.99$ and demonstrate that the model captures the non-linear convergent flows characteristic of gravitational collapse. Crucially, we show that local phase-space information alone is sufficient to distinguish between transient density fluctuations and bound collapsing cores. This data-driven approach offers a computationally efficient alternative to traditional sink-particle algorithms and provides a pathway for developing high-fidelity subgrid models for galaxy-scale simulations.

The molecular composition inside the dust sublimation zones of protoplanetary disks is mostly unknown but important to understanding terrestrial planet formation. A few molecules have been observed from this region, specifically CO, H2O, OH and SiO. The small surface area makes observing this region difficult, hence modeling is required to disentangle the innermost disk from regions further out. We model a protoplanetary disk around a Herbig-type star including the dust depleted inner region (approx. 0.1-0.3 au) and aim to investigate the chemistry of this region and explain existing and future observations. Methods. We post-process the dust and gas distribution of a magnetohydrostatic model with the radiation thermochemical code ProDiMo to study the chemistry and to produce observables. We find that the dust free inner disk is a molecular rich environment, where besides CO we also find H2, H2O and SiO. The gas temperature profile is complex and fluctuates between 700 and 2000 K, which is warm enough to produce CO overtone line emission. Next to the CO overtone lines we also find strong high J-level fundamental CO lines between 4.3 and 4.6 micron. The elemental enrichment of Si due to dust sublimation leads to 2 orders of magnitude more SiO abundance. The SiO gas has average temperatures of approx. 1000 K resulting in strong SiO overtone emission in the spectral range between 4 and 4.3 micron. We predict that the gas density in the dust depleted inner disk is high enough to allow for H2 formation, resulting in an molecular rich environment. For our representative Herbig model, the dust-depleted inner disk is responsible for at least 90% of the line emission for CO and H2O between 1 and 28 micron. Next to CO overtone lines, SiO overtone lines are expected to be an important tracer of a dust free inner disk.

We present 5 GHz e-MERLIN radio images at 50 mas resolution of the nuclear regions of the Legacy e-MERLIN Multi-band Imaging of Nearby Galaxies survey (LeMMINGs), the deepest statistically complete radio-band survey of the local Universe (<120 Mpc), consisting of 280 galaxies spanning all morphological and nuclear types. We detect nuclear radio emission above a median 5 sigma threshold of 0.33 mJy beam^-1 in 68 of 280 sources (24 percent), with core luminosities in the range 10^35 to 10^41.9 erg s^-1. The radio emission is attributed to active galactic nuclei, circumnuclear star formation, or, in the case of NGC 3690, a tidal disruption event. The brightest radio nuclei, with brightness temperatures >=10^6 K, reside in optically active galaxies such as LINERs and Seyferts. The detection rate for inactive systems (H II and absorption-line galaxies), which may host low-luminosity active galactic nuclei, is 8 percent. Most detections (78 percent) are compact (<10 pc), while the remaining 22 percent show extended jet-like features up to 380 pc. Compared to the 1.5 GHz LeMMINGs data, the 5 GHz observations provide superior resolution and spatial filtering, resolving out large-scale structures and isolating genuine nuclear emission. Our results suggest that low-luminosity active galactic nuclei are the primary manifestation of black hole activity in the local Universe in the form of compact jets and cores, with a preference for early-type hosts. The two LeMMINGs campaigns indicate that up to 30 percent of the local galaxy population hosts a radio-active nucleus, highlighting the necessity of high-resolution, high-sensitivity imaging for uncovering nuclear emission at the lowest luminosities.

Nonequilibrium chemistry is central to many astrophysical environments but remains a major computational bottleneck in simulations because solving the associated stiff ODE systems is expensive. Neural surrogates promise large speedups, yet existing studies rarely provide systematic comparisons of architectures or rigorous optimization toward both accuracy and efficiency. We introduce CODES, a principled framework for optimizing and benchmarking astrochemical surrogate models. Using CODES, we compare four neural surrogate architectures across four KROME-generated datasets spanning primordial and molecular-cloud chemistry with up to 287 reactions across 37 species. Dual-objective optimization reveals pronounced accuracy-efficiency trade-offs across architectures. Fully connected models achieve the highest accuracy and most reliable uncertainty estimates, while latent-evolution models show improved robustness under iterative prediction. Our results highlight the importance of systematic optimization and architectural comparison. The datasets, metrics, and benchmarking procedure are publicly released within CODES to enable reproducible surrogate benchmarking.

Molecular hydrogen (H$_2$) and hydrogen deuteride (HD) are key coolants in primordial gas and regulate the formation of the first stars and proto-galaxies. Recent results from the James Webb Space Telescope provide striking insights into galaxies detected at high redshifts, which are found to be significantly more abundant and luminous than expected from galaxy formation models, thus suggesting a gap in our understanding of the early Universe. Standard pathways for H$_2$ formation in the early Universe proceed through the H$^-$ and H$_2^+$ intermediates, both of which are strongly suppressed at high redshift by the cosmic microwave background. We propose an additional pathway for H2 and HD formation that could be active as early as the end of the epoch of recombination and could enable the formation of the first stars earlier than the current prediction at redshift z   30 - 20. The proposed pathway relies on the manifestation of Jahn-Teller dynamical coupling between electronic states of H$_3^+$. This coupling induces transient three-body recombination in H$^+$, H and H, and charge exchange within the charged atom-dimer complex that directly creates ground-state H$_2$ (and HD), bypassing the fragile intermediates that limit the standard primordial pathways. Our analysis shows that this mechanism could occur under the thermodynamic conditions of the post-recombination epoch, also suggesting that it might be playing a role in the active galactic nuclei feedback processes, regulating the formation rates of the first stars and the accretion rates of the first black holes. Though the global impact on galaxy formation and black-hole growth is not yet determined and will require quantitative assessment in future modeling, the mechanism offers an additional chemical route for H$_2$ and HD formation, with substantial cosmological relevance for primordial chemistry and early structure formation.

Previous studies have established that a weak-to-strong transition occurs in Alfvenic magnetohydrodynamic (MHD) turbulence as energy cascades from large to small scales. However, the spatiotemporal (frequency-wavenumber) properties of compressible MHD turbulence involving all eigenmodes, which encode the strength of nonlinear interactions, remain difficult to characterize observationally. Consequently, whether a similar weak-to-strong transition occurs in compressible turbulence remains elusive. Using a novel multi-spacecraft, polarization-based mode-decomposition technique with measurements from the Cluster spacecraft in Earth's magnetosheath, we obtain spatiotemporal power spectra of all MHD eigenmodes and present the first quantitative assessment of nonlinear frequency broadening. Our results show that slow modes exhibit a weak-to-strong transition, evolving from wave-like peaks to frequency-broadened spectra as nonlinearity increases, whereas fast modes remain weakly turbulent with narrow peaks near their eigenfrequencies. Both Alfvenic and compressible fluctuations contribute significantly to low-frequency, large-scale quasi-two-dimensional structures. These findings provide a comprehensive observational characterization of compressible turbulence across mode composition, spatiotemporal scales, and weak-strong turbulence regimes, with implications for energetic particle transport, turbulent dynamos, plasma heating, and solar wind-magnetosphere coupling.

Observations with the James Webb Space Telescope (JWST) have uncovered a substantial population of high-redshift broad-line active galactic nuclei (BLAGNs) characterized by moderate luminosities, weak X-ray emissions, and faint high-ionization lines. We propose that a subset of these BLAGNs, the so-called "little blue dots" (LBDs) are accreting at super-Eddington rates and use geometrically thick, non-advective disk models to investigate photon scattering and shadowing within the polar funnel. Our models predict extremely blue optical-UV continuum slopes and highly collimated radiation fields where isotropic-equivalent luminosities exceed the Eddington limit in the polar direction, while shadowing suppresses emission at higher inclinations. This "searchlight" configuration naturally generates a stratified ionization structure: coronal and high-excitation narrow lines are produced along the symmetry axis, while the equatorial broad-line region (BLR) remains shielded from the hardest ionizing photons. We show that the anisotropic illumination of the BLR explains the observed faintness of high-ionization lines despite strong Balmer emission. For M_BH=10^{7.5}-10^8 Msun black holes accreting at Eddington ratios  10, standard BLR conditions predict HeII 4686/Hbeta in the range of 0.08-0.28. Notably, because inherently blue disk spectra provide a much higher ratio of ionizing to optical photons than standard quasar composites, the observed large Balmer equivalent widths are matched with typical BLR covering factors without invoking enshrouded geometries. Taken together, these findings support the view that super-Eddington accretion flows, shaped by thick disk geometries, may naturally account for the ionizing SED and emission line diagnostics of high-$z$ LBDs, while offering a plausible pathway to rapid black hole growth at cosmic dawn.

Fuzzy Dark Matter with non-zero quartic self-interaction (SFDM) is shown to be a viable model for simultaneously fitting 17 dark-matter-dominated galaxies from the SPARC database with a single $(m,g)$ point in the space of boson masses and self-coupling constants: $\log_{10}\left(m \,[\mathrm{eV}/c^2] \right) = \log_{10}(1.98)-22^{+0.8}_{-0.6}$ and $\log_{10}\left(g \, [\mathrm{eV \, m}^3/kg] \right) = \log_{10}(9.08)-10^{+0.4}_{-1.2}$. This is based on the combination of an appropriately constructed static super-Gaussian profile for the inner galactic core (`soliton') region, and a Navarro-Frenk-White profile for the surrounding halo region. The explicit identification of a non-zero interaction strength may resolve issues of inconsistent constraints in non-interacting FDM. Our identification of these parameters enables the explicit {\em dynamical} reconstruction of potential host halos for such galaxies through numerical solution of the SFDM equations; we outline a {proof-of-principle procedure via merger simulations} for two galaxies (UGCA444, UGC07866), and show that this yields viable rotation curves over a dynamical period of ${\cal O}(1) \, Gyr$.

Mass loss is a key aspect of stellar evolution, particularly in evolved massive stars, yet episodic mass loss remains poorly understood. To investigate this, we need evolved massive stellar populations across various galactic environments. However, spectral classifications are challenging to obtain in large numbers, especially for distant galaxies. We addressed this by leveraging machine-learning techniques. We combined Spitzer photometry and Pan-STARRS1 optical data to classify point sources in 26 galaxies within 5 Mpc, and a metallicity range 0.07-1.36 Z$_\odot$. Gaia data release 3 (DR3) astrometry was used to remove foreground sources. Classifications are derived using a machine-learning model developed in our previous work. We report classifications for 1,147,650 sources, with 276,657 sources ( 24%) being robust. Among these are 120,479 red supergiants (RSGs;  11%). The classifier performs well even at low metallicities ( 0.1 Z$_\odot$) and distances under 1.5 Mpc, with a slight decrease in accuracy beyond  3 Mpc due to Spitzer's resolution limits. We also identified 21 luminous RSGs (log($L/L_\odot)\ge5.5$), 159 dusty yellow hypergiants in M31 and M33, as well as 6 extreme RSGs (log($L/L_\odot)\ge6$) in M31, challenging observed luminosity limits. Class trends with metallicity align with expectations, although biases exist. This catalog serves as a valuable resource for individual-object studies and James Webb Space Telescope target selection. It enables the follow-up on luminous RSGs and yellow hypergiants to refine our understanding of their evolutionary pathways. Additionally, we provide the largest spectroscopically confirmed catalog of extragalactic massive stars and candidates to date, beyond the Clouds, comprising 5,273 sources (including  330 other objects).

The Star-Forming Main Sequence (SFMS) serves as a critical framework for understanding galaxy evolution, highlighting the relationship between star formation rates (SFR) and stellar masses M_* across cosmic time. Despite its significance, the origin of the 0.3-0.4 dex dispersion in the SFMS remains a key unresolved question. Uncovering the origin of dispersion is crucial for understanding the evolution of galaxies. Using a large sample of approximately 500,000 galaxies, we reveal an unprecedented symmetry in the distribution of key structural properties-effective radius (R_{\rm e}), stellar surface density (M_*/R_{\rm e}^2), and morphology on the SFMS. This symmetry implies that galaxies with high (above SFMS) and low (below SFMS) SFRs share similar fundamental parameters. Moreover, galaxies with smaller R_{\rm e} or higher M_*/R_{\rm e}^2 exhibit greater dispersion in SFR. This dispersion reflects the response to fluctuations in cosmic accretion flows, while the SFR itself represents the time-averaged effect over the gas consumption timescale. Shorter gas consumption timescales, associated with higher M_*/R_{\rm e}^2, lead to greater SFR dispersion. Our results reveal that the variation of SFR originates from the oscillation of accretion flow and is regulated by the stellar surface density.

Using deep JWST/NIRSpec spectra from the Blue Jay survey, we perform the first systematic investigation of neutral gas content in massive galaxies at Cosmic Noon based on the Ca II H, K absorption lines. We analyze a sample of 9 galaxies at 1.8 < z < 2.8 with stellar masses > 10.6, for which we detect neutral gas absorption both in Ca II and in Na I. After removing the stellar continuum using the best-fit model obtained with Prospector, we fit the excess absorption due to neutral gas in the Ca II H, K doublet and in the Na I D doublet, together with nearby emission lines produced by ionized gas. We measure covering fractions between 0.2 and 0.9 from the Ca II H and K lines, which are spectrally well resolved in the NIRSpec R   1000 observations, unlike the absorption lines in the Na I D doublet. We measure the velocity shift, velocity dispersion, and column density separately for Ca II and Na I. About half of the galaxies present blueshifted Ca II, indicative of an outflow of neutral gas, consistent with previous results based on Na I. The velocity shift and the column density measured from Ca II are correlated with those measured from Na I, implying that these absorption lines trace gas in similar physical conditions. However, the column densities are not in a 1:1 relation, meaning that the relative amount of Ca II and Na I atoms along the line of sight varies with the gas column density. After discussing possible reasons for this behavior, we derive an empirical relation between the column density of Ca II and the column density of Na I and, in a more indirect way, of neutral hydrogen H I. This calibration offers a new way to estimate the outflow mass and the mass outflow rate for the neutral phase from current and future JWST observations of massive galaxies at Cosmic Noon and beyond

We investigate the possibility that the recently identified population of high-redshift, obscured quasars - known as "Little Red Dots" (LRDs) - originates from early black hole seed formation driven by ultra-strongly self-interacting dark matter (uSIDM). In this framework, dark matter halos undergo gravothermal core collapse due to large self-interaction cross sections, resulting in the rapid formation of massive black hole (BH) seeds with masses $\gtrsim 10^{5} M_\odot$ at redshifts $z \gtrsim 5$. We develop a semi-analytic model that tracks the evolution of the dark matter halo population, the redshift of collapse $z_{\rm coll}$, and the corresponding BH mass function. Black hole growth is modeled stochastically via a log-normal Eddington ratio distribution and a finite duty cycle. We find that the uSIDM scenario naturally reproduces key observed properties of LRDs, including their abundance, compactness, and characteristic BH masses, while offering a mechanism for early, obscured black hole formation that is difficult to achieve in standard CDM-based models. The predicted SMBH mass function at $z \sim 5$ shows excellent agreement with LRD observational data and SIDM merger-tree simulations, particularly at the high-mass end $(m_{\rm BH} \gtrsim 10^{7} M_\odot)$. These results suggest that LRDs may serve as powerful observational tracers of exotic dark sector physics and that SMBH formation in the early universe could be significantly shaped by non-gravitational dark matter interactions.

Recent JWST observations of massive galaxies at z > 2 have detected blueshifted absorption in Na I D and other resonant absorption lines, indicative of strong gas outflows in the neutral phase. However, the measured mass outflow rates are highly uncertain because JWST observations can only probe the column density of trace elements such as sodium, while most of the gas is in the form of hydrogen. The conversion between the column density of sodium and that of hydrogen is based on observations of gas clouds within the Milky Way, and has not been directly tested for massive galaxies at high redshift. In order to test this conversion, we study a unique system consisting of a massive quiescent galaxy (J1439B) at z = 2.4189 located at a projected distance of 38 physical kpc from the bright background quasar QSO J1439. The neutral outflow from the galaxy is observed as a sub-damped Lyman-alpha absorber in the spectrum of the background quasar, which enables a direct measurement of the hydrogen column density from Lyman transitions. We obtain new near-infrared spectroscopy with Magellan/FIRE and detect Na I D and other resonant absorption lines from Mg II, Mg I, and Fe II. We are thus able to derive new, empirical calibrations between the column density of trace elements and the hydrogen column density, that can be used to estimate the mass and the rate of neutral gas outflows in other massive quiescent galaxies at high redshift. The calibration we derive for Na I is only 30% lower than the local relation that is typically assumed at high redshift, confirming that the neutral outflows observed with JWST at z > 2 are able to remove a large amount of gas and are thus likely to play a key role in galaxy quenching. However, using the local calibration for Mg II yields an order-of-magnitude discrepancy compared to the empirical calibration, possibly because of variations in the dust depletion.

We characterize the JWST spectra of $61$ galaxies at $z=9-14$, including $30$ newly-confirmed galaxies. We directly compare the $z>9$ spectroscopic properties against $401$ galaxies at $6<z<9$, with the goal of identifying evolution in the star formation histories and ISM. We measure rest-UV emission line properties and UV continuum slopes, while also investigating the rest-optical emission lines for the subset of galaxies at $9.0<z<9.6$. With these spectra, we constrain the stellar masses, specific star formation rates, dust attenuation, and the average metallicity and abundance pattern of $z>9$ galaxies. Our dataset indicates that the emission lines undergo a marked change at $z>9$, with extremely large CIII], H$\beta$, and H$\gamma$ EWs becoming $2-3\times$ more common at $z>9$ relative to $6<z<9$. Using the spectra, we infer the distribution of SFRs on short (SFR$_{\rm 3Myr}$) and medium (SFR$_{\rm 3-50Myr}$) timescales, finding that rapid SFR upturns (large SFR$_{\rm 3Myr}$/SFR$_{\rm 3-50Myr}$ ratios) are significantly more likely among $z>9$ galaxies. These results may reflect a larger dispersion in UV luminosity at fixed halo mass and larger baryon accretion rates at $z>9$, although other physical effects may also contribute. We suggest that the shift in star formation conditions explains the prevalence of extreme nebular spectra that have been detected at $z>9$, with hard ionizing sources and nitrogen-enhancements becoming more typical at the highest redshifts. Finally, we identify five $z>9$ spectroscopically confirmed galaxies with red UV colors ($\beta\gtrsim-1.5$), either revealing a small population with moderate dust attenuation ($\tau_V=0.23-0.35$) or very high density nebular-dominated galaxies with hot stellar populations.

The resonant scattering nature of Ly$\alpha$ photons interacting with neutral hydrogen makes Ly$\alpha$ emitters (LAEs) robust tracers of the intergalactic neutral hydrogen fraction, and thus sensitive probes of cosmic reionization. We present an extensive study of the Ly$\alpha$ evolution from galaxies at 4.5 $\leq$ z $\leq$ 11 in the UDS field, observed as part of the CAPERS survey, and complemented with spectra from the DAWN JWST Archive. The combined sample includes 651 spectroscopically confirmed Ly$\alpha$-break galaxies, among which we find 73 S/N>3 LAEs in JWST-NIRSpec PRISM spectra. We trace the redshift evolution of the LAE fraction with EW$_0$ >25 A (X$_{\mathrm{Ly\alpha}}$) between z = 5 and z = 9, extending such an analysis to the UDS field for the first time. At z = 5 and 6, the UDS results agree with the average JWST X$_{\mathrm{Ly\alpha}}$ values from multiple fields. However, JWST measurements are consistently lower than ground-based results. To investigate this, we compare JWST observations to a population of star-forming galaxies at z$\sim$6 observed with VLT-FORS2. We find that a Ly$\alpha$ slit-loss of 35 $\pm$ 10% in JWST spectra accounts for the offset, as the resonant Ly$\alpha$ emission is more spatially extended than the stellar continuum. From z = 6 to 7, the UDS field shows a significant drop in Ly$\alpha$ visibility, from which we infer a neutral hydrogen fraction of X$_{\mathrm{HI}}$ = 0.7–0.9. Finally, we identify two robust ionized bubbles at z = 7.29 and 7.77, with radii of $R_{\mathrm{ion}}$ = 0.6 and 0.5 physical Mpc and photometric overdensities of N/$\langle$N$\rangle$ = 3 and 4, based on candidate counts down to the photometric completeness limit. Compared to the large ionized region at z$\sim$7 in the EGS field, these results indicate significant field-to-field variation, supporting a patchy, inhomogeneous reionization process.

Reconstructing the star formation history (SFH) of disk galaxies is central to understanding their growth and evolution, yet such estimates can be strongly biased by stellar radial migration over cosmic time. Using 186 Milky Way (MW) and Andromeda (M31) analogs from the TNG50 cosmological simulation, we compare star formation rates (SFRs) inferred from present-day stellar positions with those based on stellar birth radii to quantify the magnitude, spatial structure, and temporal evolution of migration-induced biases. We find that radial migration systematically produces artificial star formation in regions that had not yet formed stars. Notably,  80% of galaxies exhibit outer-disk stars older than 10 Gyr, which we find to have formed at radii interior to the outer disk and to have reached their present locations via substantial outward migration. Similar effects appear in  45% of galaxies at intermediate radii during early epochs, and in 30% of quenched inner disks within the past 4 Gyr. Migration also smooths SFHs, washing out localized bursts and suppressions by dispersing stars across neighboring radii. The strength and imprint of these distortions depend sensitively on galactic structure and evolutionary history: strong bars drive mean SFR overestimates of up to 75% in the inner disk and 150% in the outskirts; thinner, dynamically cold disks suffer average outer-disk biases up to 160%; while thick disks exhibit typical inner-disk biases up to 125%. Merger timing further modulates these patterns. Our results demonstrate that failing to account for stellar migration can lead to severe misinterpretations of when and where stars formed, with direct implications for the chemical and evolutionary histories of the MW and external galaxies.

We study the gravitational dynamics of quasi-hierarchical triple systems, where the outer orbital period is significantly longer than the inner one, but the outer orbit is extremely eccentric, rendering the time at pericentre comparable to the inner period. Such systems are not amenable to the standard techniques of perturbation theory and orbit-averaging. Modelling the evolution of these triples as a sequence of impulses at the outer pericentre, we show, by comparing with direct three-body integrations, that such triples lend themselves to a description as an analytical map between subsequent outer pericentre passages. This map exhibits secular oscillations, going beyond the von Zeipel–Lidov–Kozai mechanism. We show that the time to coalescence due to gravitational waves in such systems is modified. We then study the long-term evolution under this map, which lead to a random-walk-like behaviour of the inner eccentricity. While this behaviour is probably absent from isolated triples, it could exist in triples where the outer orbit is weakly coupled to a system with which it can exchange angular momentum, and we describe some properties of this random walk.

We model the absorption of X-rays by gas obscuring the source and photoionized by it. We consider a broad range of column densities, including both Thomson-thin and Thomson-thick media. For the Thomson thin regime, we derive a simple criterion based on the source luminosity and spectrum, as well as the medium radius and column density, that distinguishes between the following cases: (i) The absorption can be modeled well by a neutral medium; (ii) The radiation ionizes its way through the medium, and no absorption is expected; and (iii) A detailed model is required because the column density inferred from modeling the absorption with a neutral gas is much lower than the actual column density, or because the absorption features cannot be fitted by a neutral absorber. We derive the criterion analytically using a toy model of hydrogen and oxygen and calibrate it for realistic compositions with metallicities in the range $Z/Z_{\odot}=0.01-50$, using \textsc{Cloudy}. We generalize the model to the Thomson-thick regime, where we consider, alongside photoabsorption, electron scattering, Compton heating, Comptonization, and photon degradation. In this case, the emergent spectrum depends on the boundary condition experienced by photons scattered back towards the source. We discuss the effect of a reflective boundary and a reprocessing boundary. We provide simple criteria for the expected absorption state and discuss additional effects that alter the spectrum. The main motivation for our modeling is X-ray emission from supernovae interacting with the circumstellar medium; however, we expect it to be useful for many other applications.

In our previous paper, we developed an orbit-superposition method for edge-on barred galaxies and constructed a set of dynamical models based on different mock observations of three galaxies from the Auriga simulations. In this study, we adopted 12 cases with side-on bars (three simulated galaxies, each with four different projections). We decomposed these galaxies into different structures combining the kinematic and morphological properties of stellar orbits. We then compared the model-predicted components to their true counterparts in the simulations. Our models can identify (BP/X-shaped) bars, spheroidal bulges, thin discs, and spatially diffuse stellar halos. The mass fractions of bars and discs are well constrained with absolute biases of $|f_{\rm model}-f_{\rm true}|\le0.15$. We recovered the mass fractions of halos with $|f_{\rm model}-f_{\rm true}|\le0.03$. For the bulge components, 10 out of 12 cases exhibit $|f_{\rm model}-f_{\rm true}|\le0.05$, while the other two cases exhibit $|f_{\rm model}-f_{\rm true}|\le0.10$. Then, by tagging the stellar orbits with ages and metallicities, we derived the chemical properties of each structure. For the stellar ages, our models recovered the negative gradients in the bars and discs, but exhibited relatively larger uncertainties for age gradients in the bulges and halos. The mean stellar ages of all components were constrained with absolute biases $|t_{\rm model}-t_{\rm true}|\rm\lesssim1\,Gyr$. For stellar metallicities, our models reproduced the steep negative gradients of the bars and bulges, as well as all different kinds of metallicity gradients in the discs and halos. Apart from the bulge in the simulated galaxy Au-18, the mean stellar metallicities of all other components were constrained with absolute biases of $|Z_{\rm model}-Z_{\rm true}|\rm\le0.5\,Z_{\odot}$.

RUBIES-UDS-QG-z7 (RQG) is the earliest massive quiescent galaxy identified to date, inferred to have formed its abundant stellar mass in a single burst that ceases rapidly before $z\sim8$. An object of such extreme nature challenges our understanding of galaxy formation, requiring rapid growth and quenching mechanisms only $0.6 \ \rm{Gyr}$ after the Big Bang and implying number densities $2 \ \rm{dex}$ higher than currently predicted by simulations. We use synthetic observables to identify analogous systems within the First Light And Reionisation Epoch Simulations (FLARES) and find two massive galaxies ($M_{\ast}>10^{9} \ \mathrm{M_{\odot}}$) dominated by rapidly quenched bursts. One of these demonstrates excellent agreement with the inferred physical properties of RQG and implies a number density of analogous systems $\log_{10}(\mathrm{N_{Q}} \ / \ \mathrm{Mpc}^{-3}) = -7.92^{\ +0.52}_{\ -0.76}$. Beyond demonstrating that the current FLARES model is capable of producing RQG-like systems, these analogues provide a laboratory within which to study the underlying physics. Their active galactic nuclei (AGN) heat and expel gas, inducing rapid quenching and preventing timely rejuvenation. This causes above-average chemical enrichment at a given stellar mass, with super solar levels predicted for RQG. These metallicities are underestimated by spectral energy distribution fitting and we show that $\alpha$-enhancement cannot be solely responsible. Degeneracies with age and dust attenuation appear the more likely causes. Tensions between observed and simulated number densities can be alleviated in part by considering systematics, but adjustments to AGN feedback, such as allowing super-Eddington accretion rates, may be required for full agreement.

CO gas emission is a fundamental tool for measuring column density, but in cold, dark clouds, much of the CO is locked away in ice. We present JWST results from observations of a star forming filament (G0.342+0.024) that that appears to be associated with the 3 kpc arm. This filament is backlit by the Galactic Center, which has allowed us to construct a high-resolution extinction map (mean separation between stars of  1" outside the filament,  2" in the filament). ALMA Band 3 data reveals embedded star formation within the cloud. Using the CO ice feature covered by the F466N band, we map the CO ice column density of the filament. By combining the extinction map, CO ice column density map, and archival CO observations, we examine the efficacy of standard CO X-factor measurements of mass in star forming gas. We find that 50-88% of the CO is locked away in ice at large column densities ($N_{\rm \rm H_2} \gtrsim 10^{22} \rm ~cm^{-2}, 200 \rm ~M_{\odot} \rm ~pc^{-2}$) in the filament. The primary sources of uncertainty in this estimate are due to uncertainty in the ice composition and lab measurements of ice opacities. This shows that systematic corrections are needed for mass measurements in the Milky Way and nearby galaxies at high column densities.

The correlations between the positions and shapes of galaxies, i.e. intrinsic alignments, have been measured in many observational studies and hydrodynamical simulations. The position-shape correlation measurements of disk galaxies with varying methodologies, samples and hydrodynamical simulations are inconsistent in amplitude and sign. This work compares the correlations of disk and elliptical shapes around all galaxy positions and disk shapes around the positions of ellipticals at $z=0$ and $z=1$ for two different shape definitions in TNG300, Horizon-AGN and EAGLE for multiple morphological definitions in a consistent way. All types of signals are positive and robust in TNG300 and EAGLE and positive or null in Horizon-AGN. The exception to this is the negative correlation of disks around ellipticals in Horizon-AGN. This arises for reduced shapes, which down-weight the outskirts of galaxies, at $z=1$, when disks are identified via a threshold in $|v/\sigma|$, the rotational velocity over the velocity dispersion. A re-weighting of the ellipticals around all galaxies signals in TNG300, according to the underlying stellar mass distributions of the samples, highlights the importance of the influence of (sub-grid) physics at these non-linear scales.

The earliest assembly of the Milky Way remains poorly understood, yet the spatial, chemical, and kinematic properties of its most metal-poor stars provide a unique fossil record of its proto-Galaxy phase. Understanding how this ancient component formed is essential for linking near-field Galactic archaeology to high-redshift galaxy evolution. We construct the currently largest 3D map of inner-Galaxy metal-poor giants by combining several narrow/medium-band photometric surveys, reaching metallicities down to $\mathrm{[Fe/H]}\sim-3.5$. Our final sample contains 5,095,676 giants, including 1,717,610 stars with $\mathrm{[Fe/H]}<-1$. Across $-4\le \mathrm{[Fe/H]}<-1$, the density distribution reveals a centrally concentrated, flattened spheroidal component extending to $r_{\rm gc}\sim15$ kpc, together with a prominent overdensity near $X\sim-5$ kpc that is dominated by metal-poor stars on disklike orbits, with a kinematically hot background also present. The selection-function-corrected metallicity distribution function shows a distinct, very metal-poor component around $\mathrm{[Fe/H]}\sim-2.7$ that becomes most prominent at 1$<r_{\rm gc}<$3 kpc. Stars with $-3.5\lesssim\mathrm{[Fe/H]}\lesssim-1.4$ exhibit weak net rotation and low rotational support within $r_{\rm gc}<15$ kpc. Finally, we briefly note that the centrally enhanced very metal-poor component could be qualitatively consistent with one or more early dissipative build-up episodes (e.g., high-$z$ compaction/"blue-nugget" phases) as one possible interpretation.

We measure and analyze the inflows and outflows of mass, energy, and metals through the interstellar medium (ISM) and circumgalactic medium (CGM) of galaxies in the IllustrisTNG100 simulations. We identify the dominant feedback mechanism in bins of halo virial mass and redshift by computing the integrated energy input from SNe and the “kinetic” and “thermal” mode of AGN feedback. We measure all quantities in a shell at the virial radius (“halo scale”) and one chosen to be approximately at the interface of the CGM and the interstellar medium (ISM; “ISM scale”). We find that galaxies have strong net positive inflows on halo scales, and weaker but still net positive inflows on ISM scales, at $z\gtrsim 2$. At later times, partially due to the onset of kinetic AGN feedback in massive halos, inflows and outflows nearly balance one another, leading to the familiar effects of the slow-down of galaxy growth and the onset of quenching. Halos dominated by SN feedback show only weak evidence of preventative feedback on halo scales, and we see excess ISM scale accretion indicative of rapid gas recycling. Wind mass loadings decrease with increasing halo mass, and with increasing redshift, while energy loadings are nearly independent of both mass and redshift. The detailed catalogs of these mass, metal, and energy inflow and outflow rates on galaxy and halo scales can be used to guide empirical and semi-analytic models, and provide deeper insight into how galaxy growth and quenching is regulated in the IllustrisTNG simulations.

We present a comparative analysis of warm dark matter (WDM) subhalo populations generated by the semi-analytic model {\sc Galacticus} and the COZMIC suite of dark matter-only $N$-body simulations. Using a range of thermal relic WDM particle masses (3–10 keV), we examine key summary statistics – including the subhalo mass function, spatial distribution, maximum circular velocity $V_\text{max}$, and its corresponding radius $ R_\text{max} $ – to evaluate the consistency between these two modeling frameworks. Both models predict a suppression of low-mass subhalos correlated with decreasing WDM particle mass, and that WDM subhalos tend to have lower $V_\text{max} $ and larger $ R_\text{max} $ values than their CDM counterparts at fixed mass. While {\sc Galacticus} provides more statistically precise results due to a larger sample size, the COZMIC simulations display similar qualitative trends. We discuss how differences in halo finder algorithms, simulation resolution, and modeling assumptions affect subhalo statistics. Our findings demonstrate that {\sc Galacticus} can reliably reproduce WDM subhalo distributions seen in $N$-body simulations, offering a computationally efficient tool for exploring the implications of WDM across astrophysical phenomena.

The full extended Gaia mission spans slightly over 10 years, whilst the current data releases represent only a fraction of it, 34 months in Gaia's third data release (DR3). The longer baseline improves the quality of astrometric fits, lowering the noise floor and making consistently bad fits (e.g., due to binarity) more apparent. In this paper, we use simulated binaries from the Gaia Universe Model to examine the long-term astrometric behaviour of single stars and stellar binaries. We calculate nominal upper limits on the spread of goodness of astrometric fits for well-behaved single stars. Specifically, for the RUWE parameter, for upcoming DR4 ($\rm RUWE_{lim}=1.15$) and DR5 ($\rm RUWE_{lim}=1.11$), using the full mission nominal scanning law. These can be used to identify poor astrometric fits and can flag potential binary systems. We show the increase in the number and type of binaries detectable through RUWE. With our updated RUWE thresholds, the number of detectable short-period binaries increases by 5-10% with each subsequent data release, suggesting detections may be possible for orbital periods down to days. The number of detectable long-period systems increases by 10-20%, with periods up to 100 years, causing significant deviations in low and moderate-eccentricity binaries. Very eccentric systems with much longer periods (thousands of years) can still be detected if they pass through periapse during the observing window. Finally, we compare our results to the analytic estimate for the spread in UWE, which we predict from a $\chi$-distribution moderated by the number of observations. These agree with our inferred population limits but suggest that we may be biased by a small number of poorly sampled systems. In regions of the sky that are more frequently observed, lower limits could be employed, potentially bringing even more binaries above the threshold for detectability.

Direct observations of Lyman continuum (LyC) radiation from galaxies during the Epoch of Reionization (EoR) are impeded by absorption in the intergalactic medium, requiring indirect methods to infer the escape fraction of ionizing photons ($f_{\rm esc}^{\rm LyC}$). One approach is to develop and validate such methods on local analogues of the high-redshift galaxies with directly detected LyC leakage. In this study, we constrain $f_{\rm esc}^{\rm LyC}$ using a Bayesian spectral energy distribution (SED) fitting framework built on Prospector, which incorporates a non-parametric star-formation history, a flexible dust attenuation curve, self-consistent nebular emission, and fiber aperture-loss corrections. Our methodology jointly fits broadband photometry and emission line fluxes. We apply six models to the Low-redshift LyC Survey (LzLCS), a sample of local galaxies with physical properties comparable to EoR galaxies, and evaluate them based on their ability to recover the observed flux and their relative Bayesian evidence. The best-performing model is further assessed through a parameter recovery test, demonstrating that $f_{\rm esc}^{\rm LyC}$ can be recovered within uncertainties. Building on these results, we present updated $f_{\rm esc}^{\rm LyC}$ estimates for the LzLCS sample, with a median of 4\%, and values reaching as high as 51\%, with 26 of 64 galaxies having a cosmologically relevant $f_{\rm esc}^{\rm LyC}$($\gtrsim5\%$). Additionally, we present a revised UV $\beta$-slope vs $\log_{10}(f_{\rm esc}^{\rm LyC})$ relation, derived using symbolic regression with PySR trained on a synthetic dataset generated with \texttt{Prospector}:$\log_{10}(f_{\rm esc}^{\rm LyC}) = (-2.30 \pm 1.28)\beta - (6.26 \pm 2.91)$, ($\sigma = 0.43$ dex). The relation successfully reproduces the $f_{\rm esc}^{\rm LyC}$ obtained from full SED fitting of the LzLCS sample within $1\sigma$.

We model the formation of dust in the ejecta of Type Iax supernovae (SNe), which is a low-luminosity subclass of Type Ia SNe. A non-equilibrium chemical kinetic approach is adopted to trace the synthesis of molecules, molecular clusters, and dust grains in the ejecta of thermonuclear SNe. We find that Type Iax SNe provide conditions conducive to the formation of several O-rich dust species in the ejecta. Particularly, iron-rich silicates of chemical type FeSiO3, Fe2SiO4, and MgFeSiO4 are found to form in abundance, suggesting that the ejecta of low-luminosity thermonuclear SNe can be a site where a large fraction of iron is locked up in dust, unlike other stellar sources. The final mass of dust formed in the ejecta ranges between 10^{-5} and 10^{-4} Msun, where most of the dust forms between 1000 and 2000 days post-explosion. Apart from Fe-rich silicates, Mg-silicates, and silicon carbide are also formed in the ejecta of Type Iax SNe. When compared to the dust budget of typical Type Ia SNe, we find that the expected dust-to-ejecta mass ratio is 1 or 2 orders of magnitude larger in Type Iax SNe. We conclude that the ejecta of typical Type Ia SNe form a negligible amount of dust, in agreement with observation, while the low-luminosity subclass Type Iax SNe are potential producers of iron-rich silicates.

Effects of General Relativity are usually neglected in the non-linear evolution of structures, where Newtonian $N$-body simulations are traditionally employed. In the post-Friedmann expansion framework, a weak-field relativistic approximation purpose-built for cosmology, a frame-dragging gravito-magnetic vector potential arises at leading order, sourced by momentum currents. At this order, the vector potential contributes to the metric while leaving the dynamics of the matter fields unaffected, as it does not appear in the Euler equation. It can therefore be extracted a posteriori from standard N-body simulations, where the dynamics is purely Newtonian. Using the Delaunay Tessellation Field Estimator code on the IllustrisTNG simulations, here we extend previous work in order to compute the power spectrum of this vector potential down to galactic scales. The magnitude of the vector potential is two orders of magnitude larger than predicted by perturbation theory, and is a $1\% \sim 0.1\%$ effect compared to the non-linear Newtonian scalar gravitational potential. In the redshift range considered here, the gravito-magnetic effect remains subdominant, without showing any enhancement during a particular phase in the evolution of structures, aside from the continuous growth of non-linearity at low redshift. Although this seems to suggest that, within the $\Lambda$CDM model, no significant gravito-magnetic effects contribute to the non-linear evolution of cosmic structures, i.e. to the dynamics of massive particles, possible observational consequences, e.g. in lensing, deserve further exploration.

The recently reported Cosmic Himalayas (CH) – an extreme quasar overdensity at z 2 – poses an apparent challenge to the Lambda CDM framework, with a reported significance of 16.9-sigma under Gaussian assumptions. Such an event appears improbably rare, with a formal probability of P   10^-68. In this work, we investigate whether CH-like structures can naturally arise in cosmological hydrodynamic simulations. Using the CROCODILE simulation, which self-consistently models galaxy-black hole coevolution, we examine quasar clustering through two complementary approaches: the count-in-cells (CIC) statistic, which probes large-scale overdensities, and the nearest-neighbor distribution (NND), sensitive to small-scale environments. CIC analysis reveals that the underlying distribution is heavy-tailed and non-Gaussian, and that conventional Gaussian-based evaluation substantially overestimates the significance of extreme events. When modeled with an asymmetric generalized normal distribution (AGND), the inferred rarity of the CH is substantially reduced and reconciled with standard Lambda CDM; for instance, regions appearing as 12-sigma outliers under Gaussian assumptions (P   10^-33) are found to occur in the AGND regime with a probability of P   10^-4. NND analysis further demonstrates that extreme overdense regions within the simulation can naturally sustain two-point correlation function values similar to those observed in the CH (r0   30 Mpc/h), suggesting that the strong clustering stems from sample selection biases and local environmental variations. These two analyses conclusively highlight the importance of adopting non-Gaussian statistics when quantifying extreme overdensities of quasars and establish that the CH is not an anomaly, but a natural outcome of structure formation in the Lambda CDM universe.

We present multiband photometric and spectroscopic observations of supernova (SN) 2024abvb, which exhibits early-time prominent photoionized narrow emission lines of C II superposed on a blue continuum. The absence of Balmer features indicates that the SN exploded within hydrogen-poor circumstellar matter (CSM). Together with the lack of explicit evidence of helium signatures, we tentatively identify SN 2024abvb as a Type Icn SN (SN Icn). After correcting for extinction, we estimate an r-band peak absolute magnitude of -19.7, placing SN 2024abvb in the luminous regime of SNe Icn. We adopted a hybrid model that accounts for both the energy released by the ejecta-CSM interaction and the radioactive decay of nickel synthesized in the SN ejecta to fit the light curve of SN 2024abvb. The best-fit model to the multiband light curves within the first   40 days after explosion suggests that the CSM, radioactive nickel, and ejecta masses to be 0.28 Msun, < 3.8 * 10^-2 Msun, and 0.12 Msun, respectively. Such a low ejecta mass indicates that the progenitor star of SN 2024abvb experienced a significant mass-stripping process, consistent with the hydrogen-poor and helium-poor spectral features. SN 2024abvb provides important insights into the physical origins of the rare subclass of SNe Icn.

We investigate how supermassive black hole (SMBH) seeds form in the early Universe at the centers of ultralight dark matter (ULDM) halos. Focusing on the ULDM Jeans scale, we identify the critical conditions under which high-redshift baryonic gas, strongly confined by central solitonic cores of the halos, undergoes direct and monolithic collapse. The solitonic potential naturally drives rapid inflow and shock heating, allowing the gas temperature to exceed the critical atomic-cooling threshold of $\sim 10^4 \,{\rm K}$ required for fragmentation suppression without invoking an external UV background. We derive semi-analytic relations for the halo mass, soliton mass, baryonic core radius, and thermodynamic state of the gas, including the effects of baryonic contraction. These relations simultaneously determine the minimum and maximum SMBH seed masses as functions of redshift. In this framework, pristine gas clouds that satisfy the temperature threshold collapse without fragmentation, forming SMBH seeds with characteristic masses of order $ 10^5M_\odot$, while systems below the threshold are expected to form compact star clusters instead. Our model also implies an upper limit on the attainable SMBH mass, predicting a maximum mass scale of order $10^{10}M_\odot$, consistent with the most massive quasars observed to this http URL ULDM particle mass required to reproduce the inferred seed mass scale, $m \simeq 10^{-22}{\rm eV}$, coincides with the value favored by galactic-scale observations, providing a unified explanation for the characteristic masses of both galactic cores and early SMBH seeds. Our model predicts efficient SMBH seed formation at redshifts $z \gtrsim 10$ and offers a natural interpretation of recently observed little red dots as SMBHs embedded in compact, hot, ionized gas clouds.

Compton-thick active galactic nuclei (CT-AGNs), defined by column density $\mathrm{N_H} \geqslant 1.5 \times 10^{24} \ \mathrm{cm}^{-2}$, are so heavily absorbed that their X-ray emission is often feeble, even undetectable by X-ray instruments. Nevertheless, their radiation is expected to be a substantial contributor to the cosmic X-ray background (CXB), predicting that CT-AGNs should comprise at least $\sim$30% of the total AGN population. In the Cosmological Evolution Survey (COSMOS), the identified CT-AGN fraction falls far below theoretical expectations, indicating that a substantial population of CT-AGNs is hidden due to their low photon counts or their flux below the current flux limits of X-ray instruments. This work focuses on identifying CT-AGNs hidden in mid-infrared (MIR)-selected AGNs. First, we selected a sample of 1,104 MIR-selected AGNs that were covered but individually undetected by X-ray. Next, we reduced the X-ray data in the COSMOS and analyzed multiwavelength data in our sample to derive the key physical parameters required for CT-AGN identification. Using MIR diagnostics, we first find out 7 to 23 CT-AGN candidates. Their subsequent X-ray stacking analysis reveals a clear detection at $>3\sigma$ significance in the soft band and only $>1\sigma$ significance in the hard band. We fit the stacked soft- and hard-band fluxes with a physical model and confirm that these sources are absorbed by Compton-thick material. However, CT-AGNs constitute only 2.1% (23/1104) of our sample, significantly below the fraction predicted by CXB synthesis models, indicating that a considerable population of CT-AGNs remains missed by our selection. A comparison of host-galaxy properties between CT-AGNs and non-CT-AGNs reveals no significant differences.

We present an example cluster follow-up study with Euclid. Our target, a $z\sim 1.74$ candidate cluster nicknamed the `Puddle', was initially discovered by the Massive and Distant Clusters of WISE Survey 2 (MaDCoWS2) as a $z_{phot}\sim 1.65$ candidate cluster. It was also detected independently as a $z_{phot}\sim 1.5$ candidate with both cluster-finding algorithms in Euclid Quick Release 1 (Q1). A Keck MOSFIRE spectrum shows the brightest nucleus is at $z=1.74$ and is AGN-dominated. We focus our analysis on the galaxy population and the Brightest Cluster Galaxy (BCG), using a combination of Euclid and ancillary photometry. Compared to similar fields, we measure an overdensity of $110\pm 14$ galaxies with $H_\mathrm{E}\leq 22.25$ in a 2' radius around the BCG. We estimate that $18\pm 4$% of the completeness-corrected galaxy population is red, which is consistent with some clusters at $z>1.5$ but lower than others. \textit{Euclid} imaging reveals that six or seven galaxies appear to be assembling to form the future BCG. Spectral energy distribution (SED) fitting suggests that the merging BCG has a stellar mass of $5.7\pm 0.3\times 10^{11}\,M_\odot$ and experienced a short burst of star formation about $300\,$Myr ago. Its morphology, stellar mass, and star-formation history suggest that the proto-BCG is a more evolved version of the merging core of SPT2349$-$56. These systems indicate that multiobject mergers might be a common BCG formation process. Assuming a similar density of mergers in the Euclid Wide Survey, we expect that Euclid will discover approximately 400 assembling BCGs by the end of its mission.

Exploring the expansion history of the universe, understanding its evolutionary stages, and predicting its future evolution are important goals in astrophysics. Today, machine learning tools are used to help achieving these goals by analyzing transient sources, which are modeled as uncertain time series. Although black-box methods achieve appreciable performance, existing interpretable time series methods failed to obtain acceptable performance for this type of data. Furthermore, data uncertainty is rarely taken into account in these methods. In this work, we propose an uncertaintyaware subsequence based model which achieves a classification comparable to that of state-of-the-art methods. Unlike conformal learning which estimates model uncertainty on predictions, our method takes data uncertainty as additional input. Moreover, our approach is explainable-by-design, giving domain experts the ability to inspect the model and explain its predictions. The explainability of the proposed method has also the potential to inspire new developments in theoretical astrophysics modeling by suggesting important subsequences which depict details of light curve shapes. The dataset, the source code of our experiment, and the results are made available on a public repository.

One of the most interesting questions that astronomy can hope to answer is: are we alone in our Milky Way galaxy? A detection of an electromagnetic (EM) signal generated by an extraterrestrial technological intelligence (ETI), or the presence in our solar system of an alien probe, would answer this question in the affirmative. Purposeful interstellar communication is a 2-way street - the transmitting and receiving technological intelligence (TI) both need to do its part. As the receiving TI, our EM search programs should incorporate a model of what a transmitting TI is likely to be doing. Published searches for extraterrestrial technological intelligence (SETI) have generally not done so and, thus, have often been sub-optimally designed. We propose an improved search technique that more closely corresponds to astronomical surveys that have been undertaken for reasons that have nothing to do with SETI. Published non-SETI radio and optical surveys are sufficiently extensive that they already supply meaningful constraints on the prevalence of nearby purposely communicative alien civilizations. Purposeful communication can also include the sending of spaceships (probes). The absence of evidence for alien probes in the solar system suggests that no alien civilization has passed within 100 light-years of Earth during the past few billion years.

Gamma-Ray Bursts (GRBs) are the most energetic electromagnetic explosions in the universe. Recently, the Large High Altitude Air Shower Observatory (LHAASO) reported the breakthrough observation of GRB 221009A with gamma-ray energies beyond 13 TeV. This discovery, together with the previous GRB detection well above 100 GeV, confirms the production of very-high-energy (VHE, $\gtrsim 100$ GeV) radiation which might be a common component of all bright GRBs. It is reasonable to expect that bright GRBs are important targets for ground-based gamma-ray experiments. In this work, we estimate the detection rate for current and upcoming ground-based gamma-ray observatories including LHAASO, Large Array of Imaging Atmospheric Cherenkov Telescopes (LACT) and the Southern Wide-field Gamma-ray Observatory (SWGO) under two emission models with GRB 221009A as the template: first, that they all share the same intrinsic VHE spectral shape; second, they have the same environmental parameter and electron spectral index, governing their synchrotron self-Compton (SSC) emission. Using the long GRB luminosity and redshift distribution function obtained from the Fermi-GBM GRB samples, and accounting for the cosmological effects and extra-galactic background light (EBL) absorption, we derive the expected VHE flux at Earth. The sensitivity analysis for LHAASO, the upcoming LACT, and SWGO to evaluate their detection potential across specific redshift and luminosity ranges has been performed. The corresponding 5$\sigma$ detection rates of 221009A-like GRBs for the two emission models are: LHAASO, 0.04-0.05 yr$^{-1}$; LACT, 0.03-0.06 yr$^{-1}$; SWGO, 0.2-0.4 yr$^{-1}$. These rates can vary by up to $\approx 24\%$ due to different EBL models.

The rapid expansion of low-Earth-orbit (LEO) megaconstellations introduces new risks to radio astronomy from unintended electromagnetic radiation (UEMR). In this work, we present an attempt to search for UEMR from Starlink satellites using the 21 Centimeter Array (21CMA). Because the sensitivity of a single pod observation is limited, we focus on developing a robust observing and detection pipeline. Using Two-Line Element (TLE) data, we predict satellite transit times to guide the observations, and we define entry into the field of view (FoV) as an apparent declination greater than $85^{\circ}$ with respect to the 21CMA. We analyze the system equivalent flux density (SEFD) and the resulting single-pod sensitivity limits, which explain the detection of emission originating from the ORBCOMM satellites, rather than any detectable broadband UEMR in our dynamic spectra. To validate the methodology, we developed a Python package, orbdemod, to demodulate ORBCOMM downlink signals in our data. The recovered satellite ID agrees with the satellite predicted by our maximum-declination analysis, thereby validating the accuracy of our transit prediction and identification framework. Furthermore, via modulation power spectrum analysis, we show that the impulsive broadband bursts are produced by power line arcing near the array rather than by satellite UEMR.

The spaceborne laser emission telescope is a core and critical component of the space gravitational wave detection this http URL with ground-based telescopes, the on-orbit space environment is more complex and harsh, presenting higher technical challenges for the design of the optical system and structure - both optical design and structural design face considerable difficulties. To meet the requirements of space gravitational wave detection, this paper designs a laser emission telescope based on an off-axis four-mirror configuration, with a capture field of view of 300{\mu}rad, an optical transmission efficiency of 86.3%, and an optical path stability index of TTL<0.025 nm/{\mu}rad. During the design process, based on existing theories and engineering experience, the primary mirror thickness optimization and lightweight structural design were completed, and a flexible support scheme was adopted to achieve a primary mirror surface figure accuracy of 9.42 nm; the total mass of the entire telescope (excluding mirrors) is only 3.845 kg. Multi-dimensional finite element analysis was conducted on the telescope under actual working conditions: the strength of the telescope's support materials was verified under self-weight and 10G gravity loads; after removing the rigid body displacement of the mirrors using Zernike polynomials, the surface deformation of the primary mirror was controlled within 1/30 wavelength. In the thermal stability analysis, the structural deformation of the telescope under a temperature change of 100 degree celsius was simulated, and key indicators such as eccentricity and tilt between the mirrors all meet the optical design requirements. In the modal analysis, the first-order natural frequency of the telescope reaches 200 Hz under both self-weight and weightless conditions, demonstrating excellent dynamic stability.

We report new high-precision observations of the polarization of light scattered from the atmosphere of Venus, made 100 years after the pioneering studies by Bernard Lyot. The new observations include disk-integrated observations in a range of filters as well as imaging polarimetry. We compare the new results with past observations and models. We have reproduced the 1974 modelling of the Venus polarization by Hansen and Hovenier using modern radiative transfer codes. We show that the new models are in good agreement with the originals, and enable us to calculate the polarization for wavelengths not covered by the original study and to model the polarization distribution across the disk. The new observations are in good agreement with past determinations of the size distribution of the predominant particle mode. They agree with past studies in showing variability of the phase curve between synodic cycles and also polarization variability on short timescales, particularly at higher phase angles (crescent phases). Imaging polarimetry observations show good agreement with models for the redder wavelengths. However, observations in the ultraviolet show very different polarization behavior in the polar regions (within about 30 degrees of the north and south poles). The simplest explanation of this result is that there is a larger Rayleigh scattering component in the polar regions than in the equatorial and mid-latitudes and this could be explained by a lower cloud-top height in agreement with previous spacecraft observations. These ultraviolet polarization observations are inconsistent with horizontally homogeneous atmospheric models.

In ground-based high-contrast instruments, non-common path aberrations (NCPAs) limit detection performance, as they are unseen by the adaptive optics (AO) wavefront sensor but impact the astrophysical image, creating quasi-static speckles. SAXO+, the upgrade of the SAXO (SPHERE AO system) includes a second loop of AO downstream of the SAXO loop that is equipped with a near-infrared pyramid wavefront sensor whose nonlinearities, usually described with modal optical gains, might be challenging for removing quasi-static speckles. We investigated two methods of quasi-static speckle removal : NCPA compensation and a dark hole loop, behind a pyramid AO system, measuring the interest of compensating for the pyramid optical gains. We performed end-to-end numerical simulations under various astrophysical conditions. We offset the pyramid wavefront sensor operating point to apply both the speckle suppression methods, with or without optical gain calibration. We evaluated the performance by measuring the residual starlight in the coronagraph image. A by-product of our study is an on-sky calibration method of measuring the pyramid optical gains. NCPA compensation reduces the residual starlight in the coronagraph image by a factor of 20 for seeing between 0.7" and 1" for a bright star and a factor of 2 at 0.7" for a faint star. Optical gains compensation enhances the performance at poor seeing and small pyramid modulation radius with a bright star, but shows a useless or even negative impact due to estimation inaccuracies at faint targets. On the other hand, the dark hole loop reduces the residual starlight by a factor of 200. The optical gain calibration enhances the dark hole performance behind a single pyramid AO system but is useless behind the SAXO+ system. Our parametric study gives baseline values for the efficient control of the dark hole loop for the SAXO+ system.

We present a GPU-accelerated transient detection pipeline developed for time-domain surveys with the Dark Energy Camera (DECam). It enables real-time-capable image processing, incorporating science-driven candidate filtering to support rapid transient identification in time-critical observing programs. The pipeline serves as the core transient discovery engine for multiple long-term DECam programs, including the GW-MMADS gravitational-wave follow-up campaign and the DESIRT survey for intermediate-redshift transients with DESI synergy. The pipeline ingests calibrated imaging products from the DECam Community Pipeline and performs image differencing using the SFFT algorithm, coupled with CNN-based real-bogus classification, to produce science-ready transient alerts and light curves that are delivered to community brokers. We validate the pipeline using archival DECam data from the DESIRT survey. The real-bogus classifier achieves a completeness of $\sim$ 99\% of real transients while rejecting $\sim$ 96\% of subtraction artifacts, and the workflow typically reduces the candidate load to a manageable level for survey operations. With GPU acceleration, the typical processing time per DECam exposure is $\sim$ 50 s from calibrated image processing to alert generation using a modest allocation of computing resources.

Recent advances in quantum sensing – optical clocks at $5.5\times 10^{-19}$ systematic uncertainty, frequency-dependent squeezing below the standard quantum limit, quantum magnetometers approaching fundamental sensitivity limits – raise a natural question: can these technologies detect gravitational waves directly, or enhance existing detectors beyond current capabilities? We show that the answer is primarily determined by the \emph{coupling mechanism} between the gravitational wave and the sensor. Starting from the tidal Hamiltonian in Fermi normal coordinates, we identify three physically distinct coupling mechanisms and derive their transducer gains within linearized general relativity and non-relativistic quantum mechanics. Internal atomic coupling (tidal distortion of electronic wavefunctions) yields a transducer gain $G_A = 2.4\times 10^{-20}$, with vanishing first-order energy shifts for all $J=0$ clock states – a $\sim\!10^{35}$ deficit relative to laser interferometry that exceeds any projected quantum enhancement. Center-of-mass coupling (Doppler shifts from geodesic motion) reaches strain sensitivities of $\sim\!10^{-18}$, still $10^4$ above LISA requirements. Light propagation coupling (phase accumulation over macroscopic baselines) provides the enormous transducer gain that makes laser interferometry – and atom interferometry – viable. For detectors exploiting this third mechanism, we quantify how much improvement quantum sensors can provide through the detector's noise architecture: LISA's noise budget is $\sim\!91\%$ classical, limiting combined quantum enhancement to $\mathcal{E} \approx 1.04$, while ground-based detectors in the shot-noise-dominated regime achieve $\mathcal{E} = 1.8$–$2.4$. Atom interferometers exploit the same light-propagation mechanism to uniquely access the 0.01–10 Hz band.

One of the main targets for space-borne gravitational wave detectors is the detection of Extreme Mass Ratio Inspirals (EMRIs). The data analysis of EMRIs requires waveform models that are both accurate and fast. The major challenge for the fast generation of such waveforms is the generation of the Teukolsky amplitudes for generic (eccentric and inclined) Kerr orbits. The requirement for the modeling of $\sim10^5$ harmonic modes across a four-dimensional parameter space makes traditional approaches, including direct computation or dense interpolation, computationally prohibitive. To overcome this issue, we introduce a convolutional encoder-decoder architecture for a fast and end-to-end global fitting of the Teukolsky amplitudes. We also adopt a transfer learning strategy to reduce the size of the training dataset, and the model is trained gradually from the simplest Schwarzschild circular orbits to generic Kerr orbits step by step. Within this framework, we obtain a surrogate model based on a semi-analytical Post-Newtonian dataset, and the full harmonic amplitudes can be generated within milliseconds, while the median mode-distribution error for generic orbits is approximately $\sim10^{-3}$. This result indicates that the framework is viable for constructing efficient waveform models for EMRIs.

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

The planned IceCube-Gen2 radio neutrino detector at the South Pole will enhance the detection of cosmic ultra-high-energy neutrinos. It is crucial to utilize the available time until construction to optimize the detector design. A fully differentiable pipeline, from signal generation to detector response, would allow for the application of gradient descent techniques to explore the parameter space of the detector. In our work, we focus on the aspect of signal generation, and propose a modularized deep learning architecture to generate radio signals from in-ice neutrino interactions conditioned on the shower energy and viewing angle. The model is capable of generating differentiable signals with amplitudes spanning multiple orders of magnitude, as well as consistently producing signals corresponding to the same underlying event for different viewing angles. The modularized approach ensures physical consistency of the samples and leads to advantageous computational properties when using the model as part of a bigger optimization pipeline.

Kinetic inductance detectors are widely used in millimeter- and submillimeter-wave astronomy, benefiting from their fast response and relative ease of fabrication. The GroundBIRD telescope employs microwave kinetic inductance detectors at 145 and 220 GHz to observe the cosmic microwave background. As a ground-based telescope, it is subject to inherent environmental systematics, namely atmospheric emission and thermal fluctuations of the focal plane temperature. This study models resonance frequency shifts induced by each source using calibrated on-site measurements of precipitable water vapor and temperature. Comparison with observational data confirms the validity of the models and identifies atmospheric loading as the dominant contributor to frequency variation under typical observation conditions.

We propose a next-generation ground-based gravitational-wave detector, Cryogenic sub-Hz cROss torsion-bar detector with quantum NOn-demolition Speed meter (CHRONOS), optimized for the unexplored $0.1$-$10\,\mathrm{Hz}$ band between the space-based LISA and future ground-based detectors such as Cosmic Explorer and the Einstein Telescope. CHRONOS combines a ring-cavity Sagnac interferometer with torsion-bar test masses to realize the first quantum nondemolition (QND) measurement of angular momentum in a macroscopic system. By implementing a speed-meter readout in the rotational degree of freedom, CHRONOS coherently cancels quantum radiation-pressure noise and enables sub-Hz observations. We calculate, for the first time, that detuned power-recycling and cavity-length optimization can simultaneously relax technical requirements on both torsion bars and speed meters. Assuming a realistic optical design with 1m torsion bar, we estimate strain sensitivities of $h \simeq 5\times10^{-19}\,\mathrm{Hz^{-1/2}}$ at $2\,\mathrm{Hz}$ for detectors with arm lengths of $2.5$ m, $40$ m, and $300$ m. These sensitivities enable (i) direct detection of intermediate-mass black hole binaries up to 340 Mpc with SNR=3, (ii) probing SGWB down to $\Omega_{\mathrm{GW}}\sim\ 3\times10^{-4}$ at 0.2 Hz with 5 year accumulation. Furthermore, CHRONOS enable to prompt detection of gravity-gradient signals from M 5.5 earthquakes even with a $2.5$ m prototype. CHRONOS thus opens new opportunities for quantum-limited geophysical observation and multi-band, multi-messenger gravitational-wave astronomy.

Protoplanetary disks with inner dust cavities (often referred to as "transition disks") are potential signposts of planet formation. We use Gaia astrometry to search for planetary and stellar companions in a sample of 98 transition disks, assessing the occurrence rate of such companions and their potential influence on cavity formation. For the 98 Young Stellar Objects (YSOs), we compute Gaia proper motion anomalies which, together with the RUWE, identify companions with mass ratios $q \gtrsim 0.01$ at $\sim$0.1-30 au. We assess the impact of disk gravity, accretion, disk-scattered light, dippers, starspots, jets, and outflows on the measured proper motion anomalies, concluding that astrometric techniques such as the one of this work can be robustly applied to YSOs. Significant proper motion anomalies are found in 31 transition disks (32% of the sample), indicative of companions. We recover 85% of known companions within our sensitivity range. We model the semi-major axis and mass required for a companion to reproduce the observed astrometric signals. Most inferred companions have $M > 30$ M$\rm{_{J}}$, placing many within or near the stellar mass regime. Seven sources host companions compatible with a planetary mass ($M < 13$ M$\rm{_{J}}$, HD 100453, J04343128+1722201, J16102955-3922144, MHO6, MP Mus, PDS 70, and Sz 76). For the non-detections, we provide the companion masses and semi-major axes that can be excluded in future searches. About half (53%) of detected companions cannot be reconciled with having carved the observed dust cavities. We find that transition disks host as many companions within our sensitivity range as do randomly sampled groups of YSOs and main-sequence stars. If dust cavities are shaped by companions, such companions must reside at larger orbital separations than those of the companions detected here, and we predict them to be of planetary mass. [abridged]

In this paper, we present a description and characterisation of the 2-Channel Focal Reducer Rozhen (FoReRo2), with an emphasis on its polarimetric and spectropolarimetric observation modes. FoReRo2 is a multimode instrument mounted at the Ritchey-Chretien (RC) focus of the 2-meter Ritchey-Chretien-Coude (RCC) telescope at the Bulgarian National Astronomical Observatory (BNAO) - Rozhen. The primary targets of this instrument include novae, Be/X-ray binaries, symbiotic stars, asteroids, and comets. Standard stars with a high degree of polarisation are an essential part of the polarimetric observations. In our sample, HD 204827 shows short-term variability in the position angle (PA), making it unsuitable for calibration. HD 183143 displays intrinsic polarisation variability, but remains stable in terms of the PA. For the first time, we present a statistical study of the K and lambda_max parameters of Serkowski's law, including the mean value, standard deviation, and distribution. To demonstrate FoReRo2's polarimetric capabilities, we present several examples. The recurrent nova RS Oph shows the variability of Serkowski's K parameter, which is due to dust formation within the first few days after the 2021 outburst. These examples also include the imaging polarimetry of comet C/2019 Y4 (ATLAS) and spectropolarimetric observations of the symbiotic star Z And.

Precise sensing and control of spatial mode content is essential for the performance of precision optical systems, particularly interferometric gravitational-wave detectors, where misalignment and mode mismatch can lead to significant optical losses and degraded quantum noise suppression. Conventional approaches, including heterodyne wavefront sensing and phase camera techniques, are effective but can be limited by hardware complexity and systematic uncertainties arising from restricted reference-beam overlap. This paper presents a novel two-step deep learning pipeline for robust beam diagnostics based solely on beam intensity images. In the first stage, a multi-intensity-image convolutional neural network (CNN) performs accurate mode decomposition, recovering the complex modal content of distorted beams. In the second stage, the predicted mode coefficients are fed into a downstream regression network that simultaneously estimates all eight degrees of freedom (DoFs) associated with misalignment and mode mismatch, including beam tilt, lateral offset, and waist size and position mismatches in both transverse directions. The proposed CNN-based framework achieves a mean absolute error (MAE) of 0.0034 in the mode decomposition stage, which propagates to a total MAE of 0.0062 in the recovered beam imperfection parameters at the final stage. This corresponds to an average residual optical loss of 39 ppm per DoF (310 ppm total). This approach relies only on standard CCD imaging and is robust to random intensity noise, eliminating the need for complex interferometric hardware. The results demonstrate that the proposed deep learning pipeline enables real-time, high-accuracy wavefront sensing and mode-mismatch diagnostics, providing a scalable and hardware-efficient tool for improving the stability and sensitivity of precision optical systems.

The search for Galactic binary gravitational waves is a critical challenge for future space-based gravitational wave detectors, such as LISA. We propose an innovative approach to simultaneously explore gravitational waves originating from Galactic binaries by developing a new Local Maxima Particle Swarm Optimization (LMPSO) algorithm. This new approach effectively addresses the inaccuracies often associated with signal subtraction contamination, a challenge for traditional iterative subtraction methods, particularly when dealing with low signal-to-noise ratio (SNR) signals (e.g., SNR $<$ 15). We also account for the effects of overlapping signals and degeneracy noise. To demonstrate the effectiveness of our approach, we use residuals from the LISA mock data challenge (LDC1-4), where 10,982 injected sources with SNR $\ge$ 15 have been removed. For the remaining sources with SNR $<$ 15, our method successfully identifies 6,508 signals, yielding a false alarm rate of $\text{FAS}_{0.8} = 36.8\%$. By focusing on a subset of sources-specifically, those with $f > 3$ mHz and those with $f \le 3$ mHz but SNR $\ge 13$-we identify 3,406 signals, with a reduced false alarm rate of $\text{FAS}_{0.8} = 22.5\%$. We further demonstrate that, within the same detection SNR range, our method achieves a comparable or lower $\text{FAS}$ than other existing methods.

The current KM3NeT/ORCA neutrino telescope, still under construction, has not yet reached its full potential in neutrino reconstruction capability. When training any deep learning model, no explicit information about the physics or the detector is provided, thus they remain unknown to the model. This study leverages the strengths of transformers by incorporating attention masks inspired by the physics and detector design, making the model understand both the telescope design and the neutrino physics measured on it. The study also shows the efficacy of transformers on retaining valuable information between detectors when doing fine-tuning from one configurations to another.