A tight correlation between the baryonic and observed acceleration of galaxies has been reported over a wide range of mass (10^8 < M_{\rm bar}/{\rm M}_\odot < 10^{11}) - the Radial Acceleration Relation (RAR). This has been interpreted as evidence that dark matter is actually a manifestation of some modified weak-field gravity theory. In this paper, we study the radially resolved RAR of 12 nearby dwarf galaxies, with baryonic masses in the range 10^4 < M_{\rm bar}/{\rm M}_\odot < 10^{7.5}, using a combination of literature data and data from the MUSE-Faint survey. We use stellar line-of-sight velocities and the Jeans modelling code GravSphere to infer the mass distributions of these galaxies, allowing us to compute the RAR. We compare the results with the EDGE simulations of isolated dwarf galaxies with similar stellar masses in a $\Lambda$CDM cosmology. We find that most of the observed dwarf galaxies lie systematically above the low-mass extrapolation of the RAR. Each galaxy traces a locus in the RAR space that can have a multi-valued observed acceleration for a given baryonic acceleration, while there is significant scatter from galaxy to galaxy. Our results indicate that the RAR does not apply to low-mass dwarf galaxies and that the inferred baryonic acceleration of these dwarfs does not contain enough information, on its own, to derive the observed acceleration. The simulated EDGE dwarfs behave similarly to the real data, lying systematically above the extrapolated RAR. We show that, in the context of modified weak-field gravity theories, these results cannot be explained by differential tidal forces from the Milky Way, nor by the galaxies being far from dynamical equilibrium, since none of the galaxies in our sample seems to experience strong tides. As such, our results provide further evidence for the need for invisible dark matter in the smallest dwarf galaxies.

In this study, we report conclusive evidence for an ancient star cluster that has been accreted by the Large Magellanic Cloud (LMC). By leveraging observations from the Hubble Space Telescope (HST), we investigate the chrono-dynamical structure of a sample of seven old star clusters within the LMC in a self-consistent way. The multi-epoch nature of the dataset allowed the determination of high-precision proper motions for the clusters. Employing an isochrone-fitting methodology, we additionally infer from the deep high-resolution HST data homogeneous and robust estimates for their distances, ages and metallicities. Supplementing these data with literature line-of-sight velocities, we investigate the full 3-dimensional dynamics of the clusters within the frame of the LMC. With respect to the other clusters in our sample, NGC 1841 depicts a peculiar case. Its position in the age-metallicity plane, that makes it about 1 Gyr younger than the other metal-poor LMC clusters, but also its dynamical properties with a radial orbit almost perpendicular to the LMC disc plane, clearly advocates for a different origin. We thus conclude that NGC 1841 has likely been accreted by the LMC from a smaller galaxy. The other clusters in our sample show disc-like kinematics, with the case of NGC 2210 being peculiar, based on its inclined orbit. Their coherent age-metallicity relation closely resembles that of Gaia-Sausage-Enceladus globular clusters, thus suggesting a similar early evolution for the two dwarf galaxies. We do not find clear-cut chrono-kinematic evidence that NGC 2005 has been accreted by the LMC as suggested by a previous study based on its chemical abundance pattern. Regardless of its nature, its very old age illustrates that peculiar chemical evolutions already emerge at very early times.

Weakened magnetic braking (WMB) was originally proposed in 2016 to explain anomalously rapid rotation in old field stars observed by the Kepler mission. The proximate cause was suggested to be a transition in magnetic morphology from larger to smaller spatial scales. In a series of papers over the past five years, we have collected spectropolarimetric measurements to constrain the large-scale magnetic fields for a sample of stars spanning this transition, including a range of spectral types from late F to early K. During this time, we gradually improved our methods for estimating the wind braking torque in each of our targets, and for evaluating the associated uncertainties. Here, we reanalyze the entire sample with a focus on uniformity for the relevant observational inputs. We supplement the sample with two additional active stars to provide more context for the evolution of wind braking torque with stellar Rossby number (Ro). The results demonstrate unambiguously that standard spin-down models can reproduce the evolution of wind braking torque for active stars, but WMB is required to explain the subsequent abrupt decrease in torque as Ro approaches a critical value for dynamo excitation. This transition is seen in both the large-scale magnetic field and the X-ray luminosity, indicating weakened coronal heating. We interpret these transitions as evidence of a rotational threshold for the influence of Coriolis forces on global convective patterns and the resulting inefficiency of the global stellar dynamo.

Between 40$\%$ and 80$\%$ of cataclysmic variables (CVs) are expected to have evolved past the period-minimum and contain a degenerate donor. However, observational surveys for CVs have only been able to detect a few of these highly evolved "period-bouncers", most likely due to the intrinsic faintness associated with their predicted low mass accretion rates. We have produced an initial selection of 137 high-likelihood period-bounce candidates from WD catalog based on our multiwavelength period-bouncer scorecard and selection cuts including X-ray data from the extended ROentgen Survey with an Imaging Telescope Array (eROSITA) on board the Spektrum-Roentgen-Gamma spacecraft (SRG). We have laid out three main requirements (classification as a CV, determination of an orbital period and detection of a very late-type donor) that should result in the confirmation of several of these candidates. Our path for confirming these candidates has already produced its first successful result with the confirmation of GALEX J125751.4-283015 as a new period-bouncer. Several other candidates have already fulfilled at least one of our three requirements making their future confirmation likely. Our search for period-bouncers using the X-ray eROSITA emission of objects in optical WD catalogs has led to the confirmation of six new period-bouncers identified from the Gaia DR3 WD catalog (five previously known CVs and one WD candidate), a 18$\%$ increase that brings the present population to 39 systems. Both the selection method for period-bounce candidates and the confirmation path that we outlined will aid in future searches for new period-bounce candidates, contributing to the goal of resolving the discrepancy between the predicted high number of period-bouncers and the low number of these systems successfully observed to date.

While cosmic reionization has been broadly constrained by global observables, the interplay between internal sources (Milky Way, M31, and their satellites) and external ionization fronts remains poorly understood in a realistic Local Group (LG) context. To address this issue, we perform radiative transfer post-processing on the original HESTIA LG constrained simulation. We calibrate our source models using a uniform 1024^3 particle, dark-matter only, HESTIA simulation coupled with a subgrid collapse-fraction model to match the global reionization observables. These source models are then applied to the HESTIA zoom-in simulations, which consist of a 4096^3 particle effective resolution in the zoom region centered on the Milky Way (MW) and M31 haloes, which resolves haloes down to 10^8 solar masses. We find that in all scenarios, reionization within the LG proceeds in an inside-out manner with the progenitors of the MW and M31 having 50 percent of their material ionized by z ~ 9-8.6, significantly earlier than the global midpoint at z ~ 7-7.7, noting that external fronts from large-scale structure play a negligible role, even under the most permissive feedback model. We further show that present-day satellite galaxies exhibit only a weak correlation between their reionization redshift and their present-day radial distance from their host halo, with somewhat tighter trends around M31 than the MW. Finally, we find that satellites which assembled before reionization are systematically more massive today, suggesting that the oldest stellar populations preferentially reside in the most massive subhaloes.

We perform a frequentist analysis using the standard profile likelihood method for clustering measurements from Data Release 1 of the Dark Energy Spectroscopic Instrument (DESI). While Bayesian inferences for Effective Field Theory models of galaxy clustering can be highly sensitive to the choice of priors for extended cosmological models, frequentist inferences are not susceptible to such effects. We compare Bayesian and frequentist constraints for the parameter set $\{\sigma_8, H_0, \Omega_{\rm{m}}, w_0, w_a\}$ when fitting to the full-shape of the power spectrum multipoles, the post-reconstruction Baryon Acoustic Oscillation (BAO) measurements, as well as external datasets from the CMB and type Ia supernovae measurements. Bayesian prior effects are very significant for the $w_0w_a$CDM model; while the $1 \sigma$ frequentist confidence intervals encompass the maximum a posteriori (MAP), the Bayesian credible intervals almost always exclude the maximum likelihood estimate (MLE) and the MAP - indicating strong prior volume projection effects - unless supernovae data are included. We observe limited prior effects for the $\Lambda$CDM model, due to the reduced number of parameters. When DESI full-shape and BAO data are jointly fit, we obtain the following $1\sigma$ frequentist confidence intervals for $\Lambda$CDM ($w_0w_a$CDM): $\sigma_8 = 0.867^{+0.048}_{-0.041} , \ H_0 = 68.91^{+0.80}_{-0.79} \ \rm{km \ s^{-1}Mpc^{-1}} , \ \Omega_{\rm{m}} = 0.3038\pm0.0110$ ($\sigma_8 = 0.793^{+0.069}_{-0.048} , \ H_0 = 64.9^{+4.8}_{-2.8} \ \rm{km \ s^{-1}Mpc^{-1}} , \ \Omega_{\rm{m}} = 0.369^{+0.029}_{-0.059}$ , $w_0 = -0.24^{+0.17}_{-0.64}$ , $w_a = -2.5^{+1.9}_{}$), corresponding to 0.7$\sigma$, 0.3$\sigma$, 0.7$\sigma$ (1.9$\sigma$, 3.4$\sigma$, 5.6$\sigma$, 5.5$\sigma$, 5.6$\sigma$) shifts between the MLE relative to the Bayesian posterior mean for $\Lambda$CDM ($w_0w_a$CDM) respectively.

We introduce the THESAN-ZOOM project, a comprehensive suite of high-resolution zoom-in simulations of $14$ high-redshift ($z>3$) galaxies selected from the THESAN simulation volume. This sample encompasses a diverse range of halo masses, with $M_\mathrm{halo} \approx 10^8 - 10^{13}~\mathrm{M}_\odot$ at $z=3$. At the highest-resolution, the simulations achieve a baryonic mass of $142~\mathrm{M}_\odot$ and a gravitational softening length of $17~\mathrm{cpc}$. We employ a state-of-the-art multi-phase interstellar medium (ISM) model that self-consistently includes stellar feedback, radiation fields, dust physics, and low-temperature cooling through a non-equilibrium thermochemical network. Our unique framework incorporates the impact of patchy reionization by adopting the large-scale radiation field topology from the parent THESAN simulation box rather than assuming a spatially uniform UV background. In total, THESAN-ZOOM comprises $60$ simulations, including both fiducial runs and complementary variations designed to investigate the impact of numerical and physical parameters on galaxy properties. The fiducial simulation set reproduces a wealth of high-redshift observational data such as the stellar-to-halo-mass relation, the star-forming main sequence, the Kennicutt-Schmidt relation, and the mass-metallicity relation. While our simulations slightly overestimate the abundance of low-mass and low-luminosity galaxies they agree well with observed stellar and UV luminosity functions at the higher mass end. Moreover, the star-formation rate density closely matches the observational estimates from $z=3-14$. These results indicate that the simulations effectively reproduce many of the essential characteristics of high-redshift galaxies, providing a realistic framework to interpret the exciting new observations from JWST.

Ginnungagap is a fully parallel (MPI+OpenMP) code designed to generate cosmological initial conditions for simulations involving very large numbers of particles. It operates in several modes, including the creation of initial conditions with either uniform or spatially varying resolution (for "zoom-in" simulations). The initial conditions can be fully random or derived by extending the resolution of existing ones while preserving the large-scale structures. Ginnungagap is open source and modular, consisting of a collection of independent tools that can be used for a variety of tasks. In this paper, we describe the main features of Ginnungagap and present test results for different types of simulations prepared with it.

Solar eclipses offer unparalleled opportunities for public engagement in astronomy. Large groups of people often gather to view eclipses, and these events require affordable and easy to use tools to safely observe the Sun. One unique way to observe a solar eclipse is by using a disco ball. Here, we present an analysis of the experiences of educators who used a disco ball as a solar projector during various public outreach events. Through a survey conducted shortly after the April 2024 total solar eclipse and the March 2025 partial solar eclipse, we collected data on the use, engagement, and perceived educational value of a disco ball projector from 31 individual events. The results suggest that disco balls were not only affordable and safe, but also popular and educational.

The third gaia data release (DR3) provides a wealth of new data products. The early part of the release, Gaia EDR3, already provided the astrometric and photometric data for nearly two billion sources. The full release now adds improved parameters compared to Gaia DR2 for radial velocities, astrophysical parameters, variability information, light curves, and orbits for Solar System objects. The improvements are in terms of the number of sources, the variety of parameter information, precision, and accuracy. For the first time, Gaia DR3 also provides a sample of spectrophotometry and spectra obtained with the Radial Velocity Spectrometer, binary star solutions, and a characterisation of extragalactic object candidates. Before the publication of the catalogue, these data have undergone a dedicated transversal validation process. The aim of this paper is to highlight limitations of the data that were found during this process and to provide recommendations for the usage of the catalogue. The validation was obtained through a statistical analysis of the data, a confirmation of the internal consistency of different products, and a comparison of the values to external data or models. Gaia DR3 is a new major step forward in terms of the number, diversity, precision, and accuracy of the Gaia products. As always in such a large and complex catalogue, however, issues and limitations have also been found. Detailed examples of the scientific quality of the Gaia DR3 release can be found in the accompanying data-processing papers as well as in the performance verification papers. Here we focus only on the caveats that the user should be aware of to scientifically exploit the data.

Metal absorbers represent a powerful probe of galaxy feedback and reionization, as highlighted by both observational and theoretical results showing an increased abundance of low-ionised metal species at higher redshifts. The origin of such absorbers is currently largely unknown because of the low number of galaxy counterparts detected, suggesting that they might be surrounded by low-mass faint sources that fall below the detection threshold of current instruments. We use the THESAN-ZOOM radiation hydrodynamic simulations to investigate the connection between properties of neutral oxygen (OI) absorbers and galaxies within the redshift range $z = 5 - 8$. We find that the circumgalactic medium of galaxies becomes progressively ionised with cosmic time, leading to a decrease of $\approx 0.2$ in the covering fraction of neutral oxygen, while the total oxygen covering fraction remains constant. The observable absorbers ($N_{\rm OI} \gtrsim 10^{13}\,\text{cm}^{-2}$) are not confined to haloes: at $z \geq 5$ the majority ($\gtrsim 60\%$) arise beyond $R_{\rm{vir}}$, and including these systems is essential to reproduce the observed increase in absorber incidence with redshift. The simulated absorbers preferentially reside in overdensities rich in low-mass galaxies ($M_\star \leq 10^8\,\rm{M}_\odot$), explaining the scarcity of detected counterparts, while not excluding the possibility of nearby star-forming sources ($\geq 5\,\text{M}_\odot\,\text{yr}^{-1}$) similar to those suggested by the latest ALMA observations and, at larger distances, by the JWST. These results establish OI absorbers as sensitive tracers of the evolving ionisation structure around faint galaxies to be probed by forthcoming deep spectroscopic surveys.

[abridged] Accurately accounting for the AGN phase in galaxy evolution requires a large, clean AGN sample. This is now possible with SRG/eROSITA. The public Data Release 1 (DR1, Jan 31, 2024) includes 930,203 sources from the Western Galactic Hemisphere. The data enable the selection of a large AGN sample and the discovery of rare sources. However, scientific return depends on accurate characterisation of the X-ray emitters, requiring high-quality multiwavelength data. This paper presents the identification and classification of optical and infrared counterparts to eRASS1 sources using Gaia DR3, CatWISE2020, and Legacy Survey DR10 (LS10) with the Bayesian NWAY algorithm and trained priors. Sources were classified as Galactic or extragalactic via a Machine Learning model combining optical/IR and X-ray properties, trained on a reference sample. For extragalactic LS10 sources, photometric redshifts were computed using Circlez. Within the LS10 footprint, all 656,614 eROSITA/DR1 sources have at least one possible optical counterpart; about 570,000 are extragalactic and likely AGN. Half are new detections compared to AllWISE, Gaia, and Quaia AGN catalogues. Gaia and CatWISE2020 counterparts are less reliable, due to the surveys shallowness and the limited amount of features available to assess the probability of being an X-ray emitter. In the Galactic Plane, where the overdensity of stellar sources also increases the chance of associations, using conservative reliability cuts, we identify approximately 18,000 Gaia and 55,000 CatWISE2020 extragalactic sources. We release three high-quality counterpart catalogues, plus the training and validation sets, as a benchmark for the field. These datasets have many applications, but in particular empower researchers to build AGN samples tailored for completeness and purity, accelerating the hunt for the Universe most energetic engines.

High-redshift active galactic nuclei (AGN) have long been recognized as key probes of early black hole growth and galaxy evolution. However, modeling this population remains difficult due to the wide range of luminosities and black hole masses involved, and the high computational costs of capturing the hydrodynamic response of gas and evolving radiation fields on-the-fly. In this study, we present a new suite of simulations based on the IllustrisTNG galaxy formation framework, enhanced with on-the-fly radiative transfer, to examine AGN at high redshift (z > 5) in a protocluster environment extracted from the MillenniumTNG simulation. We focus on the co-evolution of black holes and their host galaxies, as well as the radiative impact on surrounding intergalactic gas. The model predicts that black holes form in overdense regions and lie below the local black hole-stellar mass relation, with stellar mass assembly preceding significant black hole accretion. Ionizing photons are primarily produced by stars, which shape the morphology of ionized regions and drive reionization. Given the restrictive black hole growth in the original IllustrisTNG model, we reduce the radiative efficiency from 0.2 to 0.1, resulting in higher accretion rates for massive black holes, more bursty growth, and earlier AGN-driven quenching. However, the resulting AGN remain significantly fainter than observed high-redshift quasars. As such, to incorporate this missing population, we introduce a quasar boosted model, in which we artificially boost the AGN luminosity. This results in strong effects on the surrounding gas, most notably a proximity effect, and large contributions to He ionization.

We have built and characterized the first six-telescope discrete beam combiner (DBC) for stellar interferometry in the astronomical J-band. It is the DBC with the largest number of beam combinations and was manufactured using ultrafast laser inscription (ULI) in borosilicate glass, with a throughput of approximately 56%. For calibration of the visibility-to-pixel matrix (V2PM), we use a two-input Michelson interferometer and extract the complex visibility. A visibility amplitude of 1.05 and relative precision of 2.9% and 3.8% are extracted for 1328 nm and 1380 nm, respectively. Broadband (< 40 nm) characterization is affected by dispersion, but shows similar performance.

This study explores the dynamical impact of cosmic rays (CRs) in Milky Way-like galaxies using the Rhea simulation suite. Cosmic rays, with their substantial energy density, influence the interstellar medium (ISM) by supporting galactic winds, modulating star formation, and shaping ISM energetics. The simulations incorporate a multi-phase ISM, self-consistent CR transport in the advection-diffusion approximation, and interactions with magnetic fields to study their effect on galaxy evolution. Key findings reveal that CRs reduce star formation rates, and drive weak but sustained outflows with mass loading factors of $\sim0.2$, transporting a substantial fraction (20%-60%) of the injected CR energy. These CR-driven outflows are launched not just from the galactic center but across the entire disk, illustrating their pervasive dynamical influence. Galactic disks supported by CRs exhibit broader vertical structures compared to magnetic-field-dominated setups, though the scale heights are similar. CR feedback enhances magnetic flux transport to the circumgalactic medium (CGM), leading to a magnetically enriched CGM with field strengths of $\sim0.5\mu\mathrm{G}$ while reducing gas temperatures to $\lesssim10^5\,\mathrm{K}$. The CR energy is relatively smoothly distributed in the disk, with gradient lengths exceeding the typical size of molecular clouds, indicating that the CR behavior is not adiabatic.

We present a new integral-field spectroscopic dataset of the central part of the Orion Nebula (M 42), observed with the MUSE instrument at the ESO VLT. We reduced the data with the public MUSE pipeline. The output products are two FITS cubes with a spatial size of ~5.9'x4.9' (corresponding to ~0.76 pc x 0.63 pc) and a contiguous wavelength coverage of 4595...9366 Angstrom, spatially sampled at 0.2". We provide two versions with a sampling of 1.25 Angstrom and 0.85 Angstrom in dispersion direction. Together with variance cubes these files have a size of 75 and 110 GiB on disk. They represent one of the largest integral field mosaics to date in terms of information content. We make them available for use in the community. To validate this dataset, we compare world coordinates, reconstructed magnitudes, velocities, and absolute and relative emission line fluxes to the literature and find excellent agreement. We derive a two-dimensional map of extinction and present de-reddened flux maps of several individual emission lines and of diagnostic line ratios. We estimate physical properties of the Orion Nebula, using the emission line ratios [N II] and [S III] (for the electron temperature $T_e$) and [S II] and [Cl III] (for the electron density $N_e$), and show two-dimensional images of the velocity measured from several bright emission lines.

(abridged) We have obtained precise radial velocities for a sample of 373 G and K type giants at Lick Observatory regularly over more than 12 years. Planets have been identified around 15 giant stars; an additional 20 giant stars host planet candidates. We investigate the occurrence rate of substellar companions around giant stars as a function of stellar mass and metallicity. We probe the stellar mass range from about 1 to beyond 3 M_Sun, which is not being explored by main-sequence samples. We fit the giant planet occurrence rate as a function of stellar mass and metallicity with a Gaussian and an exponential distribution, respectively. We find strong evidence for a planet-metallicity correlation among the secure planet hosts of our giant star sample, in agreement with the one for main-sequence stars. However, the planet-metallicity correlation is absent for our sample of planet candidates, raising the suspicion that a good fraction of them might indeed not be planets. Consistent with the results obtained by Johnson for subgiants, the giant planet occurrence rate increases in the stellar mass interval from 1 to 1.9 M_Sun. However, there is a maximum at a stellar mass of 1.9 +0.1/-0.5 M_Sun, and the occurrence rate drops rapidly for masses larger than 2.5-3.0 M_Sun. We do not find any planets around stars more massive than 2.7 M_Sun, although there are 113 stars with masses between 2.7 and 5 M_Sun in our sample (corresponding to a giant planet occurrence rate < 1.6% at 68.3% confidence in that stellar mass bin). We also show that this result is not a selection effect related to the planet detectability being a function of the stellar mass. We conclude that giant planet formation or inward migration is suppressed around higher mass stars, possibly because of faster disk depletion coupled with a longer migration timescale.

We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.

Characterizing the evolution of the star-forming main sequence (SFMS) at high redshift is crucial to contextualize the observed extreme properties of galaxies in the early Universe. We present an analysis of the SFMS and its scatter in the THESAN-ZOOM simulations, where we find a redshift evolution of the SFMS normalization scaling as $\propto (1+z)^{2.64\pm0.03}$, significantly stronger than is typically inferred from observations. We can reproduce the flatter observed evolution by filtering out weakly star-forming galaxies, implying that current observational fits are biased due to a missing population of lulling galaxies or overestimated star-formation rates. We also explore star-formation variability using the scatter of galaxies around the SFMS ($\sigma_{\mathrm{MS}}$). At the population level, the scatter around the SFMS increases with cosmic time, driven by the increased importance of long-term environmental effects in regulating star formation at later times. To study short-term star-formation variability, or ''burstiness'', we isolate the scatter on timescales shorter than 50 Myr. The short-term scatter is larger at higher redshift, indicating that star formation is indeed more bursty in the early Universe. We identify two starburst modes: (i) externally driven, where rapid large-scale inflows trigger and fuel prolonged, extreme star formation episodes, and (ii) internally driven, where cyclical ejection and re-accretion of the interstellar medium in low-mass galaxies drive bursts, even under relatively steady large-scale inflow. Both modes occur at all redshifts, but the increased burstiness of galaxies at higher redshift is due to the increasing prevalence of the more extreme external mode of star formation.