Verifying the fully kinematic nature of the cosmic microwave background (CMB) dipole is of fundamental importance in cosmology. In the standard cosmological model with the Friedman-Lemaitre-Robertson-Walker (FLRW) metric from the inflationary expansion the CMB dipole should be entirely kinematic. Any non-kinematic CMB dipole component would thus reflect the preinflationary structure of spacetime probing the extent of the FLRW applicability. Cosmic backgrounds from galaxies after the matter-radiation decoupling, should have kinematic dipole component identical in velocity with the CMB kinematic dipole. Comparing the two can lead to isolating the CMB non-kinematic dipole. It was recently proposed that such measurement can be done using the near-IR cosmic infrared background (CIB) measured with the currently operating Euclid telescope, and later with Roman. The proposed method reconstructs the resolved CIB, the Integrated Galaxy Light (IGL), from Euclid's Wide Survey and probes its dipole, with a kinematic component amplified over that of the CMB by the Compton-Getting effect. The amplification coupled with the extensive galaxy samples forming the IGL would determine the CIB dipole with an overwhelming signal/noise, isolating its direction to sub-degree accuracy. We develop details of the method for Euclid's Wide Survey in 4 bands spanning 0.6 to 2 mic. We isolate the systematic and other uncertainties and present methodologies to minimize them, after confining the sample to the magnitude range with negligible IGL/CIB dipole from galaxy clustering. These include the required star-galaxy separation, accounting for the extinction correction dipole using the method newly developed here achieving total separation, accounting for the Earth's orbital motion and other systematic effects. (Abridged)

Massive black holes (MBHs) are typically hosted in the centres of massive galaxies but they appear to become rarer in lower mass galaxies, where nuclear star clusters (NSCs) frequently appear instead. The transition region, where both an MBH and NSC can co-exist, has been poorly studied to date and only a few dozen galaxies are known to host them. One avenue for detecting new galaxies with both an MBH and NSC is to look for accretion signatures of MBHs. Here, we use new SRG/eROSITA all-sky survey eRASS:4 data to search for X-ray signatures of accreting MBHs in NSCs, while also investigating their combined occupation fraction. We find significant detections for 18 galaxies (~8.3%), including one ultra-luminous X-ray source; however, only three galaxies (NGC2903, 4212, and 4639) have X-ray luminosities that are higher than the expected value from X-ray binaries, indicative of the presence of an MBH. In addition, the X-ray luminosity of six galaxies (NGC2903, 3384, 4321, 4365, 4639, and 4701) differs from previous studies and could indicate the presence of a variable active galactic nucleus. The combined occupation fraction of accreting MBHs and NSCs becomes non-zero for galaxy masses above ~10^7.5 M_sun and this result is slightly elevated as compared to the literature data. Our data extend, for the first time, towards the dwarf elliptical galaxy regime and identify promising MBH candidates for higher resolution follow-up observations. At most galaxy masses (and with the exception of three cases), the X-ray constraints are consistent with the expected emission from binary systems or an Eddington fraction of at most 0.01%, assuming a black holes mass of 10^6.5 M_sun. This work confirms the known complexities in similar-type of studies, while providing the appealing alternative of using X-ray survey data of in-depth observations of individual targets with higher resolution instruments.

The CGM around unobscured AGN has received much attention in recent years. Comparatively, nebulae associated with obscured AGN are less studied. Here, we simulate the Ly$\alpha$, H$\alpha$, and HeII nebulae around the two types of AGN at $z=2-3$ with ten massive systems from the FIRE simulations based on the unified model to show their differences and to test if they can be used to constrain the AGN model. We post-process the data with the CLOUDY and the Ly$\alpha$ radiative transfer code, RASCAS. Overall, we find that the Ly$\alpha$ nebulae around the unobscured AGN (type-I nebulae) and obscured AGN (type-II nebulae) do not exhibit significant differences in the luminosity, area, and HeII/Ly$\alpha$ when the simulated cutout is set to the halo virial radius. Whereas, the type-II nebulae exhibit less symmetric morphologies, flatter surface brightness profiles, and larger emission line widths (at $R\geq 10$ kpc) than those of the type-I nebulae. These nebulae properties exhibit complicated correlations with the AGN, indicating that nebulae observations can be applied to constrain the AGN engine. However, independent observations on nebulae in the mentioned emissions are insufficient to test the unified model as a priori in observations is not possible to know the direction and opening angle of the ionization cone. We prompt that the joint observations of Ly$\alpha$ nebulae and radio jets can help to reveal the ionization cone to probe the unified model. Our calculations suggest that this method requires $\geq 75$ type-II Ly$\alpha$ nebulae with current instruments to reach a confidence level of $\geq 95\%$.

We create the first large-scale mock spectroscopic survey of gas absorption sightlines traversing the interstellar medium (ISM), circumgalactic medium (CGM), and intergalactic medium (IGM) surrounding galaxies of virtual Universes. That is, we create mock, or synthetic, absorption spectra by drawing lines-of-sight through cosmological hydrodynamical simulations, using a new mesh-free Voronoi ray-tracing algorithm. The result is the Synthetic Absorption Line Spectral Almanac (SALSA), which is publicly released on a feature-rich online science platform (this http URL). It spans a range of ions, transitions, instruments, observational characteristics, assumptions, redshifts, and simulations. These include, but are not limited to: (ions) HI, OI, CI, MgI, MgII, FeII, SiII, CaII, ZnII, SiIII, SiIV, NV, CII, CIV, OVI; (instruments) SDSS-BOSS, KECK-HIRES, UVES, COS, DESI, 4MOST, WEAVE, XSHOOTER; (model choices) with/without dust depletion, noise, quasar continua, foregrounds; (redshift) from z=0 to z~6; (ancillary data) integrated equivalent widths, column densities, distances and properties of nearby galaxies; (simulations) IllustrisTNG including TNG50, TNG-Cluster, EAGLE, and SIMBA. This scope is not fixed, and will grow and evolve with community interest and requests over time -- suggestions are welcome. The resulting dataset is generic and broadly applicable, enabling diverse science goals such as: (i) studies of the underlying physical gas structures giving rise to particular absorption signatures, (ii) galaxy-absorber and halo-absorber correlations, (iii) virtual surveys and survey strategy optimization, (iv) stacking experiments and the identification of faint absorption features, (v) assessment of data reduction methods and completeness calculations, (vi) inference of physical properties from observables, and (vii) apples-to-apples comparisons between simulations and data.

The population of the Little Red Dots (LRDs) may represent a key phase of supermassive black hole (SMBH) growth. A cocoon of dense excited gas is emerging as key component to explain the most striking properties of LRDs, such as strong Balmer breaks and Balmer absorption, as well as the weak IR emission. To dissect the structure of LRDs, we analyze new deep JWST/NIRSpec PRISM and G395H spectra of FRESCO-GN-9771, one of the most luminous known LRDs at $z=5.5$. These reveal a strong Balmer break, broad Balmer lines and very narrow [O III] emission. We unveil a forest of optical [Fe II] lines, which we argue is emerging from a dense ($n_{\rm H}=10^{9-10}$ cm$^{-3}$) warm layer with electron temperature $T_{\rm e}\approx7000$ K. The broad wings of H$\alpha$ and H$\beta$ have an exponential profile due to electron scattering in this same layer. The high $\rm H\alpha:H\beta:H\gamma$ flux ratio of $\approx10.4:1:0.14$ is an indicator of collisional excitation and resonant scattering dominating the Balmer line emission. A narrow H$\gamma$ component, unseen in the other two Balmer lines due to outshining by the broad components, could trace the ISM of a normal host galaxy with a star formation rate $\sim5$ M$_{\odot}$ yr$^{-1}$. The warm layer is mostly opaque to Balmer transitions, producing a characteristic P-Cygni profile in the line centers suggesting outflowing motions. This same layer is responsible for shaping the Balmer break. The broad-band spectrum can be reasonably matched by a simple photoionized slab model that dominates the $\lambda>1500$ Å continuum and a low mass ($\sim10^8$ M$_{\odot}$) galaxy that could explain the narrow [O III], with only subdominant contribution to the UV continuum. Our findings indicate that Balmer lines are not directly tracing gas kinematics near the SMBH and that the BH mass scale is likely much lower than virial indicators suggest.

Major mergers of galaxies are likely to trigger bursty star formation activities. The accumulation of dense gas and the boost of star formation efficiency (SFE) are considered to be the two main drivers of the starbursts. However, it is still unclear how each process operates on the scale of individual star-forming clouds. Here, we present a high-resolution (2 Msun) RHD simulation of a gas-rich dwarf galaxy merger using the RIGEL model to investigate how mergers affect the properties of the structure of dense star-forming gas and the cloud-scale SFE. We track the evolution of sub-virial dense clouds in the simulation by mapping them across successive snapshots taken at intervals of 0.2 Myr. We found that the merger triggers a 130 times higher SFR and shortens the galaxy-wide gas-depletion time by two orders of magnitude compared to those of two isolated galaxies. However, the depletion time of individual clouds and their lifetime distribution remained unchanged over the simulation period. The cloud life cycles and cloud-scale SFE are determined by the local stellar feedback rather than environmental factors regardless of the merger process, and the integrated SFE ($\epsilon_{\rm int}$) of clouds in complex environments is still well described by an $\epsilon_{\rm int}-\Sigma_{\rm tot}$ relation found in idealized isolated-cloud experiments. During the peak of the starburst, the median integrated SFE only changed by 0.17-0.33 dex lower compared to the value when the galaxies are not interacting. The merger boosts the SFR primarily through the accumulation and compression of dense gas fueling star formation. Strong tidal torques assemble $>10^ 5$ Msun clouds that seed massive star clusters. The average separation between star-forming clouds decreases during the merger, which in turn decreases the cloud--cluster spatial decorrelation from >1 kpc to 0.1 kpc depicted by tuning fork diagrams.

While the James Webb Space Telescope (JWST) now allows identifying quiescent galaxies (QGs) out to early epochs, the photometric selection of quiescent galaxy candidates (QGCs) and the derivation of key physical quantities are highly sensitive to the assumed star-formation histories (SFHs). We aim to quantify how the inclusion of JWST/MIRI data and different SFH models impacts the selection and characterisation of QGCs. We test the robustness of the physical properties inferred from the spectral energy distribution (SED) fitting, such as M*, age, star formation rate (SFR), and AV, and study how they impact the quiescence criteria of the galaxies across cosmic time. We perform SED fitting for ~13000 galaxies at z<6 from the CEERS/MIRI fields with up to 20 optical-mid infrared (MIR) broadband coverage. We implement three SFH prescriptions: flexible delayed, NonParametric, and extended Regulator. For each model, we compare results obtained with and without MIRI photometry and dust emission models. We evaluate the impact of these configurations on the number of candidate QGCs, selected based on rest UVJ colours, sSFR and main-sequence offset, and on their key physical properties such as M*, AV, and stellar ages. The number of QGCs selected varies significantly with the choice of SFH from 171 to 224 out of 13000 galaxies, depending on the model. This number increases to 222-327 when MIRI data are used (up to ~45% more QGCs). This enhancement is driven by improved constraints on dust attenuation and M*. We find a strong correlation between AV and M*, with massive galaxies (M*~10^11 M\odot) being 1.5-4.2 times more attenuated in magnitude than low-mass systems (M*~10^9 M\odot), depending on SFH. Regardless of the SFH assumption, ~13% of QGCs exhibit significant attenuation (AV > 0.5) in support of recent JWST studies challenging the notion that quiescent galaxies are uniformly dust-free.

It is often understood that supernova (SN) feedback in galaxies is responsible for regulating star formation and generating gaseous outflows. However, a detailed look at their effect on the local interstellar medium (ISM) on small mass scales in simulations shows that these processes proceed in clearly distinct channels. We demonstrate this finding in two independent simulations with solar-mass resolution, LYRA and RIGEL, of an isolated dwarf galaxy. Focusing on the immediate environment surrounding SNe, our findings suggest that the large-scale effect of a given SN on the galaxy is best predicted by its immediate local density. Outflows are driven by SNe in diffuse regions expanding to their cooling radii on large (\sim kpc) scales, while dense star-forming regions are disrupted in a localized (\sim pc) manner. However, these separate feedback channels are only distinguishable at very high numerical resolutions capable of following scales <<10^3 \msun. On larger scales, ISM densities are greatly mis-estimated, and differences between local environments of SNe become severely washed out. We demonstrate the practical implications of this effect by comparing with a mid-resolution simulation (Mptcl. \sim 200 \msun) of the same dwarf using the SMUGGLE model. The coarse-resolution simulation cannot self-consistently determine whether a given SN is responsible for generating outflows or suppressing star formation, suggesting that emergent galaxy physics such as star formation regulation through hot-phase outflows is fundamentally unresolvable by subgrid stellar feedback models, without appealing directly to simulations with highly resolved ISM.

We present a rigorous identification of candidates for dormant black holes (BHs) and neutron stars (NSs) in binaries using summary statistics from Gaia DR3, rather than full orbital solutions. Although Gaia astrometric orbits have already revealed a small sample of compact object binaries, many systems remain undetected due to stringent quality cuts imposed on the published orbits. Using a forward-modelling framework that simulates Gaia observables, in particular the renormalised unit weight error (ruwe) and radial velocity (RV) scatter, we infer posterior distributions for companion mass and orbital period via MCMC sampling, marginalising over nuisance orbital parameters. We validate our approach by comparing the predicted masses and periods against full orbit solutions from DR3, and by successfully recovering known compact object binaries as promising candidates. The method is best suited for systems with red giant primaries, which have more reliable Gaia RV scatter and a light centroid more likely dominated by one component, compared to main-sequence stars. And they are less likely to be triples with short-period inner binaries, which produce confounding signatures. We apply the method to three million giants and identify 556 systems with best-fit companion masses $\gtrsim 3\,M_\odot$. Recovery simulations suggest our selection method is substantially more sensitive than the DR3 non-single-star catalogue, particularly for binaries with periods below 1 year and above $\sim 6$ years. These candidates represent promising targets for spectroscopic follow-up and Gaia DR4 analysis to confirm the presence of compact objects. Candidate main-sequence stars with massive companions face a larger set of confounding effects. Therefore, we present an analogous catalogue of 279 additional `main sequence' candidates only as an appendix.

CARMENES is a dual-channel high-resolution spectrograph at the 3.5 m Calar Alto telescope designed to detect low-mass planets around late-type dwarfs by measuring their radial velocities (RVs). High thermal stability in both the visible (VIS) and near infrared channels is essential to achieve the precision required for these measurements. In particular, stabilising the NIR channel to the millikelvin level, which operates at cryogenic temperatures (140 K), poses significant engineering this http URL CARMENES-PLUS project was initiated to improve the instruments intrinsic RV precision. In this article, we focus on the thermal stability improvements made to the NIR channels cooling system. The NIR cooling system was originally conceived to operate with a discontinuous flow of cryogenic nitrogen gas. As part of CARMENES-PLUS, this was upgraded to a continuous flow configuration. Additional changes included the installation of an automatic vacuum system, a proportional control valve, and a pressure regulation system. These upgrades were designed to reduce thermal fluctuations and enhance long-term stability. The implemented upgrades significantly improved the intrinsic RV precision of the NIR channel. We quantified this improvement using Fabry Perot calibration spectra, obtaining an intrinsic RV precision of 0.67 ms after the interventions, an improvement of nearly 2 ms . We also assessed the stability of the nightly zero points, finding a reduced scatter of 3.9 ms post upgrade, compared to 6.1 ms before. For a sample of slowly rotating stars (vsin i below 2 kms), the median scatter decreased from 8.8 ms to 6.7 ms after the upgrades. These results demonstrate that the thermal control upgrades introduced in CARMENES PLUS have enhanced the NIR channels RV performance, bringing it closer to the VIS channels stability and reinforcing CARMENES capabilities for exoplanet detection around M dwarfs.

Exoplanetary systems with multiple giant planets present an opportunity to understand planet formation, migration processes, and long-term system-wide dynamical interactions. In particular, they provide constraints to distinguish between smooth disk-driven migration or more dynamically excited system evolution pathways. We report the discovery and characterization of a unique multi-planet system hosting three gas giant planets orbiting the post-main sequence star TOI-375. The innermost planet, TOI-375 b, was initially detected by the TESS mission and then confirmed with photometric follow-up observations conducted using MEarth and LCOGT, and radial velocity measurements obtained with FEROS and CHIRON. The radial velocity data revealed the presence of two additional planetary candidates, TOI-375 c and TOI-375 d. We find that TOI-375 b is a hot Jupiter with an orbital period of $9.45469 \pm 0.00002$ days, mass $0.745 \pm 0.053,M_\mathrm{J}$, radius $0.961 \pm 0.043, R_\mathrm{J}$, and eccentricity $0.087 \pm 0.042$. The outer two planets, TOI-375 c and TOI-375 d, are warm-cold and cold Jupiters with orbital periods of $115.5^{+2.0}_{-1.6}$ days and $297.9^{+28.9}_{-18.6}$ days, and minimum masses of $2.11 \pm 0.22, M_\mathrm{J}$ and $1.40 \pm 0.28, M_\mathrm{J}$, respectively. In terms of formation and overall system architecture, the physical properties of TOI-375 b are consistent with the core accretion scenario, while the current configuration of the system could be explained by both disk-driven and high-eccentricity migration scenarios. The discovery of TOI-375 as the first known system hosting three or more fully evolved gas giants, with at least one transiting planet, makes it an excellent candidate for testing formation and migration theories.

We review the state of the evidence for the existence and observational appearance of supermassive black hole binaries. Such objects are expected from standard hierarchical galaxy evolution to form after two galaxies, each containing a supermassive black hole, have merged, in the centre of the merger remnant. A complex interaction is predicted to take place with stars and gas in the host galaxy, leading to observable signatures in weakly as well as actively accreting phases. Direct observational evidence is available and shows examples of dual active galactic nuclei from kpc scales down to parsec scales. Signatures of possibly closer supermassive black hole binaries may be seen in jetted black holes. The interaction with stars and gas in a galaxy significantly affects the hardening of the binary and hence contributes to uncertainties of the expected gravitational wave signal. The Laser Interferometer Space Antenna (LISA) should in the future detect actual mergers. Before the launch of LISA, pulsar timing arrays may have the best chance to detect a gravitational wave signal from supermassive black hole binaries. The first signs of the combined background of inspiralling objects might have been seen already.

The 21cm line from the spin-flip transition of neutral hydrogen (HI) provides a unique window into the Epoch of Reionization (EoR), the final phase transition of our Universe. The Square Kilometre Array (SKA) enables precise measurements of 21cm fluctuations that trace ionization, temperature, and density fluctuations of the intergalactic medium (IGM). Nevertheless, a direct reconstruction of the timeline of the EoR in terms of the progress of ionization remains an ongoing challenge due to the highly non-Gaussian nature and thus intractable likelihood of the 21cm signal. Here, we present EoRFlow, a simulation-based inference (SBI) framework for reconstructing the global neutral hydrogen fraction $x_{\mathrm{HI}}(z)$ directly from 2D cylindrically averaged power spectra (2DPS) of the 21cm signal. We validate our method on realistic mock datasets for SKA-Low. Bypassing the need for explicit likelihood formulations, our approach enables fast, unbiased posterior estimation of the $x_{\mathrm{HI}}$ evolution in narrow redshift slices, allowing for piecewise reconstruction of the global reionization history. By directly inferring the reionization history from 21cm power spectra, our framework provides a scalable and robust path forward for 21cm cosmology in the SKA era.

We report the discovery and characterization of three transiting giant planets in the TIC118798035 system. The three planets were identified as transiting candidates from data of the TESS mission, and confirmed with ground-based photometric transit observations along with radial velocity variations obtained with FEROS, HARPS and ESPRESSO. The three planets present transit timing variations (TTVs). We performed a N-body orbital fitting to the TTVs and radial velocities finding that TIC118798035 b is as warm low-density Neptune with a mass of 0.0250$\pm$0.0023 $M_J$, a radius of 0.655$\pm$0.018 $R_J$, and an orbital period of 11.507 d; TIC118798035 c is a warm Saturn with a mass of 0.403$\pm$0.024 $M_J$, a radius of 0.973$\pm$0.023 $R_J$, and an orbital period of 22.564 d; and TIC118798035 d is a warm Jupiter with a mass of 0.773$\pm$0.052 $M_J$, a radius of 0.923$\pm$0.044 $R_J$, and an orbital period of 48.925 d. The bulk metallicities of the three planets don't fully follow the mass-metallicity correlation found for the giant planets of the solar system, which hints at a somewhat different formation history for the planets of the TIC118798035 system. TIC118798035 is the only system having more than two transiting planets larger than 0.5 $R_J$ with a precise orbital and physical characterization, amenable for future atmospheric studies.

We explore constraints on gauge bosons of a weakly coupled $U(1)_{B-L}$, $U(1)_{L_\mu-L_e}$, $U(1)_{L_e-L_\tau}$ and $U(1)_{L_\mu-L_\tau}$. To do so we apply the full constraining power of experimental bounds derived for a hidden photon of a secluded $U(1)_{X}$ and translate them to the considered gauge groups. In contrast to the secluded hidden photon that acquires universal couplings to charged Standard Model particles through kinetic mixing with the photon, for these gauge groups the couplings to the different Standard Model particles can vary widely. We take finite, computable loop-induced kinetic mixing effects into account, which provide additional sensitivity in a range of experiments. In addition, we collect and extend limits from neutrino experiments as well as astrophysical and cosmological observations and include new constraints from white dwarf cooling. We discuss the reach of future experiments in searching for these gauge bosons.

Little Red Dots (LRDs) are compact sources at $z>5$ discovered through JWST spectroscopy. Their spectra exhibit broad Balmer emission lines ($\gtrsim1000\rm~km~s^{-1}$), alongside absorption features and a pronounced Balmer break -- evidence for a dense, neutral hydrogen medium with the $n=2$ state. When interpreted as arising from AGN broad-line regions, inferred black hole masses from local scaling relations exceed expectations given their stellar masses, challenging models of early black hole--galaxy co-evolution. However, radiative transfer effects in dense media may also impact the formation of hydrogen emission lines. We model three scattering processes shaping hydrogen line profiles: resonance scattering by hydrogen in the $n=2$ state, Raman scattering of UV radiation by ground-state hydrogen, and Thomson scattering by free electrons. Using 3D Monte Carlo radiative transfer simulations with multi-branching resonance transitions, we examine their imprint on line shapes and ratios. Resonance scattering produces strong deviations from Case B flux ratios, clear differences between H$\alpha$ and H$\beta$, and encodes gas kinematics in line profiles but cannot broaden H$\beta$ due to conversion to Pa$\alpha$. While Raman scattering can yield broad wings, scattering of UV continuum is disfavored given the absence of strong FWHM variations across transitions. Raman scattering of higher Lyman-series emission can produce H$\alpha$/H$\beta$ wing width ratios of $\gtrsim1.28$, agreeing with observations. Thomson scattering can reproduce the observed $\gtrsim1000~\rm km\, s^{-1}$ wings under plausible conditions, e.g., $T_{\rm e} \sim 10^4\rm \, K$ and $N_{\rm e}\sim10^{24}\rm~cm^{-2}$ -- and lead to black hole mass overestimates by factors $\gtrsim10$. Our results provide a framework for interpreting hydrogen lines in LRDs and similar systems.