Papers for Monday, Sep 02 2024

Discoveries of close-in young puffy (R$_p \gtrsim$ 6 R$_\oplus$) planets raise the question of whether they are bona fide hot Jupiters or puffed-up Neptunes, potentially placing constraints on the formation location and timescale of hot Jupiters. Obtaining mass measurements for these planets is challenging due to stellar activity and noisy spectra. Therefore, we aim to provide independent theoretical constraints on the masses of these young planets based on their radii, incident fluxes, and ages, benchmarking to the planets of age $<$1 Gyr detected by Kepler, K2 and TESS. Through a combination of interior structure models, considerations of photoevaporative mass loss, and empirical mass-metallicity trends, we present the range of possible masses for 24 planets of age $\sim$10-900 Myr and radii $\sim$6-16 R$_\oplus$. We generally find that our mass estimates are in agreement with the measured masses and upper limits where applicable. There exist some outliers including super-puffs Kepler-51 b, c and V1298 Tau d, b, e, for which we outline their likely formation conditions. Our analyses demonstrate that most of the youngest planets ($\lesssim$ 100 Myr) tend to be puffed-up, Neptune-mass planets, while the true hot Jupiters are typically found around stars aged at least a few hundred Myr, suggesting the dominant origin of hot Jupiters to be late-stage high eccentricity migration.

Transiting planets around young stars are hard to find due to the enhanced stellar activity. Only a few transiting planets have been detected around stars younger than 100 Myr. We initially detected a transit-like signal in the K2 light curve of a very cool M dwarf star (EPIC 211101996) in the Pleiades open cluster, with an estimated age of about 100 Myr. Our detailed analysis of the per-pixel light curves, detrending with the Wōtan software and transit search with the Transit Least Squares algorithm showed that the source of the signal is a contaminant source (Gaia DR3 66767847894609792) 20" west of the target. The V-like shape of its phase-folded light curve and eclipse depth of ~15% suggest that it is a grazing eclipsing binary. The contaminant has hitherto been listed as a single star, which we now identify as an eclipsing stellar binary with a period of about 6 days.

The Virgo Cluster is the nearest substantial cluster of galaxies to the Milky Way and a cornerstone of the extragalactic distance scale. Here, we present JWST/NIRCam observations that simultaneously cover the cores and halos of ten galaxies in and around the Virgo Cluster and are designed to perform simultaneous measurements of the tip of the red giant branch (TRGB) and surface brightness fluctuations (SBF). Seven of the targets are within the Virgo Cluster and where we are able to resolve some of the cluster's substructure, while an additional three provide important constraints on Virgo infall. The seven galaxies within Virgo itself all have SBF measurements from the Advanced Camera for Surveys Virgo Cluster Survey (ACSVCS). After adjusting the ACSVCS measurements for the offset from our new JWST TRGB measurements, we determine a distance to the Virgo Cluster of d $=$ 16.17 $\pm$ 0.25 (stat) $\pm$ 0.47 (sys) Mpc. The work presented here is part of a larger program to develop a Population II distance scale through the TRGB and SBF that is completely independent of the prominent Cepheid + Type Ia supernova ladder. The TRGB distances to the galaxies presented here, when combined with future SBF measurements, will provide a crucial step forward for determining whether or not systematic errors can explain the well-known "Hubble tension'' or if there is significant evidence for cracks in the $\Lambda$CDM model.

The MIRI Exoplanets Orbiting White dwarfs (MEOW) Survey is a cycle 2 JWST program to search for exoplanets around dozens of nearby white dwarfs via infrared excess and direct imaging. In this paper, we present the detection of mid-infrared excess at 18 and 21 microns towards the bright (V = 11.4) metal-polluted white dwarf WD 0310-688. The source of the IR excess is almost certainly within the system; the probability of background contamination is $<0.1\%$. While the IR excess could be due to an unprecedentedly small and cold debris disk, it is best explained by a $3.0^{+5.5}_{-1.9}$ M$_{\rm Jup}$ cold (248$^{+84}_{-61}$ K) giant planet orbiting the white dwarf within the forbidden zone (the region where planets are expected to be destroyed during the star's red giant phase). We constrain the source of the IR excess to an orbital separation of 0.1-2 AU, marking the first discovery of a white dwarf planet candidate within this range of separations. WD 0310-688 is a young remnant of an A or late B-type star, and at just 10.4 pc it is now the closest white dwarf with a known planet candidate. Future JWST observations could distinguish the two scenarios by either detecting or ruling out spectral features indicative of a planet atmosphere.

The low-frequency component of the upcoming Square Kilometre Array Observatory (SKA-Low) will be sensitive enough to construct 3D tomographic images of the 21-cm signal distribution during reionization. However, foreground contamination poses challenges for detecting this signal, and image recovery will heavily rely on effective mitigation methods. We introduce \texttt{SERENEt}, a deep-learning framework designed to recover the 21-cm signal from SKA-Low's foreground-contaminated observations, enabling the detection of ionized (HII) and neutral (HI) regions during reionization. \texttt{SERENEt} can recover the signal distribution with an average accuracy of 75 per cent at the early stages ($\overline{x}_\mathrm{HI}\simeq0.9$) and up to 90 per cent at the late stages of reionization ($\overline{x}_\mathrm{HI}\simeq0.1$). Conversely, HI region detection starts at 92 per cent accuracy, decreasing to 73 per cent as reionization progresses. Beyond improving image recovery, \texttt{SERENEt} provides cylindrical power spectra with an average accuracy exceeding 93 per cent throughout the reionization period. We tested \texttt{SERENEt} on a 10-degree field-of-view simulation, consistently achieving better and more stable results when prior maps were provided. Notably, including prior information about HII region locations improved 21-cm signal recovery by approximately 10 per cent. This capability was demonstrated by supplying \texttt{SERENEt} with ionizing source distribution measurements, showing that high-redshift galaxy surveys of similar observation fields can optimize foreground mitigation and enhance 21-cm image construction.

We present a novel method to constrain the formation time of the Milky Way disc using the chrono-kinematic signatures of the inner Galaxy. We construct an O-rich Mira variable sample from the Gaia Long-period Variable catalogue to study the kinematic behaviour of stars with different ages in the inner Galaxy. From the Auriga suite of cosmological zoom-in simulations, we find that the age of the oldest stellar population with imprints of the bar in density and kinematics matches the disc spin-up epoch. This is because stars born before the spin-up show insufficient rotation and are not kinematically cold enough to be efficiently trapped by the bar. We find that the bar kinematic signature disappears for Mira variables with a period shorter than 190 days. Using the period-age relation of Mira variables, we constrain the spin-up epoch of the Milky Way to be younger than $\sim11-12$~Gyr (redshift $\sim3$). We also discuss and compare our method and result to other evidence of the Milky Way spin-up epoch under the context of a realistic age uncertainty. Age uncertainty leads to an overestimation of the disc formation time when performing backward modelling. Our constrain of the spin-up epoch is independent from previous studies because it relies on the kinematics of the inner Galaxy instead of the solar vicinity.

Resolved rotation curves (RCs) are our best probe of the dark matter distribution around individual galaxies. However their acquisition is resource-intensive, rendering them impractical for large-scale surveys and studies at higher redshift. Spatially integrated HI flux profiles on the other hand are observationally abundant and also probe dynamics across the whole HI disc. Despite this, they are typically only studied using the highly compressed linewidth summary statistic, discarding much of the available information. Here we construct a Bayesian model to infer halo properties from the full shape of the spatially integrated 21-cm line profile of a galaxy, utilising all the available information. We validate our model by assessing the consistency of halo parameters obtained from the flux profile with those obtained from RC fits for galaxies where both are available, finding good agreement provided the profile is not strongly asymmetric. We study the relative constraining power (quantified using the Kullback--Leibler divergence of the posterior from the prior), finding the flux profile inference recovers posteriors on generalised Navarro--Frenk--White halo parameters on average three times tighter than those from the linewidth, and in some cases as tight as those from resolved RCs. Finally we introduce and validate a probabilistic empirical model for the spatial distribution of HI, enabling our model to be applied to datasets for which no spatially resolved HI information is available. As the next-generation of HI observatories comes online, our framework will enable mass modelling in new regimes, with particular utility for constraining the dark matter content of galaxies across cosmic time.

We use the inferred evolution of Lyman-$\alpha$ luminosities of galaxies in the redshift range of $z \sim 9-16$ from the $James \ Webb \ Space \ Telescope$ (JWST) data to predict the power spectrum of Loeb-Rybicki haloes formed by scattered Lyman-$\alpha$ radiation in the diffuse intergalactic hydrogen before reionization. We find excellent prospects for a statistical detection of the intensity mapping signal with current and future experiments such as the JWST, SPHEREx and CDIM. We also describe the detectability of the signal in cross-correlation with the 21-cm emission from atomic hydrogen in the intergalactic medium at these redshifts. We find that the cross-correlation signal should be detectable at a significance of a few tens of standard deviations out to $z \sim 16$ using the above experiments in combination with the Square Kilometre Array (SKA)-LOW and its pathfinder, the Murchison Widefield Array (MWA).

The energy released by AGN has the potential to heat or remove the gas of the ISM, thus likely impacting the cold molecular gas reservoir of host galaxies at first, with star formation following on longer timescales. Previous works on high-z galaxies have yielded conflicting results, possibly due to selection biases and other systematics. To provide a reliable benchmark for galaxy evolution models at cosmic noon (z=1-3), two surveys were conceived: SUPER and KASHz, both targeting unbiased X-ray-selected AGN at z>1 that span a wide bolometric luminosity range. In this paper, we assess the effects of AGN feedback on the molecular gas content of host galaxies in a statistically robust, uniformly selected, coherently analyzed sample of AGN at z=1-2.6, drawn from the KASHz and SUPER surveys. By using ALMA data in combination with dedicated SED modeling, we retrieve CO and FIR luminosity as well as $M_*$ of SUPER and KASHz AGN. We selected non-active galaxies from PHIBBS, ASPECS and multiple ALMA/NOEMA surveys of sub-mm galaxies. By matching the samples in z, $M_*$ and $L_{FIR}$, we compared the properties of AGN and non-active galaxies within a Bayesian framework. We find that AGN hosts at given $L_{FIR}$ are on average CO depleted compared to non-active galaxies, confirming what was previously found in the SUPER survey. Moreover, the molecular gas fraction distributions of AGN and non-active galaxies are statistically different, with that of of AGN being skewed to lower values. Our results indicate that AGN can indeed reduce the total cold molecular gas reservoir of their host galaxies. Lastly, by comparing our results with predictions from three cosmological simulations (TNG, Eagle and Simba) filtered to match the observed properties, we confirm already known discrepancies and highlight new ones between observations and simulations.[Abridged]

We present the updated design of HighSpec, a high-resolution $\mathcal{R} \sim 20,000$ spectrograph designed for the Multi Aperture Spectroscopic Telescope (MAST). HighSpec offers three observing modes centered at the Ca II H&K, Mg b triplet, and H$\alpha$ lines. Each mode is supported by a highly optimized ion-etched grating, contributing to an exceptional instrument peak efficiency of $\gtrsim85\%$ for the two latter bands ($\gtrsim55\%$ for the Ca II H&K band). Optimizing throughput over wavelength coverage ($\Delta \lambda=10-17$ nm), HighSpec enables the precise measurement of spectral lines from faint targets. This approach is especially relevant for stellar object studies, specifically of WDs, which are intrinsically faint and have few spectroscopic lines. Each observing mode was tailored to target spectral features essential for WD research. Its integration with MAST, an array of 20 custom-designed telescopes that can function as a single large telescope (equivalent to a $2.7$ m telescope in collecting area) or multiplexing over the entire sky, provides unique adaptability for extensive and effective spectroscopic campaigns. Currently in its final assembly and testing stages, HighSpec's on-sky commissioning is scheduled for 2025.

Pair-instability supernovae (PISNe) have crucial implications for many astrophysical topics, including the search for very massive stars, the black hole mass spectrum, and galaxy chemical enrichment. To this end, we need to understand where PISNe are across cosmic time, and what are their favourable galactic environments. We present a new determination of the PISN rate as a function of redshift, obtained by combining up-to-date stellar evolution tracks from the PARSEC and FRANEC codes, with an up-to-date semi-empirical determination of the star formation rate and metallicity evolution of star-forming galaxies throughout cosmic history. We find the PISN rate to exhibit a huge dependence on the model assumptions, including the criterion to identify stars unstable to pair production, and the upper limit of the stellar initial mass function. Remarkably, the interplay between the maximum metallicity at which stars explode as PISNe, and the dispersion of the galaxy metallicity distribution, dominates the uncertainties, causing a $\sim$ seven-orders-of-magnitude PISN rate range. Furthermore, we show a comparison with the core-collapse supernova rate, and study the properties of the favourable PISN host galaxies. According to our results, the main contribution to the PISN rate comes from metallicities between $\sim 10^{-3}$ and $10^{-2}$, against the common assumption that views very-low-metallicity, Population III stars as exclusive or dominant PISN progenitors. The strong dependencies we find offer the opportunity to constrain stellar and galaxy evolution models based on possible future (or the lack of) PISN observations.

Extended Ly$\alpha$ emission is routinely found around single quasars (QSO) across cosmic time. However, few studies have investigated how such emission changes in fields with physically associated QSO pairs, which should reside in dense environments and are predicted to be linked through intergalactic filaments. We present VLT/MUSE snapshot observations (45 min./source) to unveil extended Ly$\alpha$ emission on scales of the circumgalactic medium (CGM) around the largest sample of physically associated QSO pairs to date, encompassing 8 pairs (14 observed QSOs) at $z$~3 with $i$-band magnitude between 18 and 22.75. The pairs are either at close (~50-100 kpc, 5 pairs) or wide (~450-500 kpc, 3 pairs) separation with velocity differences of $\Delta$v < 2000 km s$^{-1}$. We detect extended emission around 12 of the 14 targeted QSOs and investigate the luminosity, size, kinematics and morphology of these Ly$\alpha$ nebulae. On average, they span 90 kpc and are 2.8 $\times 10^{43}$ erg s$^{-1}$ bright. Irrespective of the QSOs' projected distance, the nebulae often (~45 %) extend toward the other QSO in the pair, leading to asymmetric emission whose flux-weighted centroid is at an offset position from any QSO location. We show that large nebulae are preferentially aligned with the large-scale structure as traced by the two QSOs, and conclude that the cool gas (10$^4$ K) in the CGM traces well the direction of cosmic web filaments. Additionally, the radial profile of the Ly$\alpha$ surface brightness around QSO pairs can be described by a power law with a shallower slope (~$-1.6$) with respect to single QSOs (~$-2$), indicative of increased CGM densities out to large radii and/or enhanced contribution from the intergalactic medium (IGM). The sample presented in this study contains excellent targets for ultra-deep observations to directly study filamentary IGM structures in emission.

Characterization of search selection effects comprises a core element of gravitational-wave data analysis. Knowledge of selection effects is needed to predict observational prospects for future surveys and is essential in the statistical inference of astrophysical source populations from observed catalogs of compact binary mergers. Although gravitational-wave selection functions can be directly measured via injection campaigns -- the insertion and attempted recovery of simulated signals added to real instrumental data -- such efforts are computationally expensive. Moreover, the inability to interpolate between discrete injections limits the ability to which we can study narrow or discontinuous features in the compact binary population. For this reason, there is a growing need for alternative representations of gravitational-wave selection functions that are computationally cheap to evaluate and can be computed across a continuous range of compact binary parameters. In this paper, we describe one such representation. Using pipeline injections performed during Advanced LIGO & Advanced Virgo's third observing run (O3), we train a neural network emulator for $P(\mathrm{det}|\theta)$, the probability that given a compact binary with parameters is successfully detected, averaged over the course of O3. The emulator captures the dependence of $P(\mathrm{det}|\theta)$ on binary masses, spins, distance, sky position, and orbital orientation, and it is valid for compact binaries with components masses between $1$--$100\,M_\odot$. We test the emulator's ability to produce accurate distributions of detectable events, and demonstrate its use in hierarchical inference of the binary black hole population.

A common setting in astronomy is the availability of a small number of high-quality observations, and larger amounts of either lower-quality observations or synthetic data from simplified models. Time-domain astrophysics is a canonical example of this imbalance, with the number of supernovae observed photometrically outpacing the number observed spectroscopically by multiple orders of magnitude. At the same time, no data-driven models exist to understand these photometric and spectroscopic observables in a common context. Contrastive learning objectives, which have grown in popularity for aligning distinct data modalities in a shared embedding space, provide a potential solution to extract information from these modalities. We present Maven, the first foundation model for supernova science. To construct Maven, we first pre-train our model to align photometry and spectroscopy from 0.5M synthetic supernovae using a constrastive objective. We then fine-tune the model on 4,702 observed supernovae from the Zwicky Transient Facility. Maven reaches state-of-the-art performance on both classification and redshift estimation, despite the embeddings not being explicitly optimized for these tasks. Through ablation studies, we show that pre-training with synthetic data improves overall performance. In the upcoming era of the Vera C. Rubin Observatory, Maven serves as a Rosetta Stone for leveraging large, unlabeled and multimodal time-domain datasets.

Radial velocity (RV) follow-up is a critical complement of transiting exoplanet surveys like the Transiting Exoplanet Survey Satellite (TESS ), both for validating discoveries of exoplanets and measuring their masses. Stellar activity introduces challenges to interpreting these measurements because the noise from the host star, which is often correlated in time, can result in high RV uncertainty. A robust understanding of stellar activity and how its timescales interact with the observing cadence can optimize limited RV resources. For this reason, in the era of over-subscribed, high-precision RV measurements, folding stellar activity timescales into the scheduling of observation campaigns is ideal. We present gaspery, an open-source code implementation to enable the optimization of RV observing strategies. Gaspery employs a generalized formulation of the Fisher Information for RV time series that also incorporates information about stellar correlated noise. We show that the information contained in an observing strategy can be significantly affected by beat frequencies between the orbital period of the planet, the stellar rotation period, and the observation epochs. We investigate how the follow-up observing strategy will affect the resulting radial velocity uncertainty, as a function of stellar properties such as the spot decay timescale and rotation period. We then describe two example use cases for gaspery: 1) calculating the minimum number of observations to reach an uncertainty tolerance in a correlated noise regime and 2) finding an optimal strategy given a fixed observing budget. Finally, we outline a prescription for selecting an observing strategy that is generalizable to different targets.

Here we critically examine the geophysical and geochemical properties of the Moon in order to identify the extent to which dynamical scenarios satisfy these observations. New joint inversions of existing lunar geophysical data (mean mass, moment of inertia, and tidal response) assuming a laterally- and vertically homogeneous lunar mantle show that, in all cases, a core with a radius of 300$\pm$20 km ($\sim$0.8 to 1.5 % the mass of the Moon) is required. However, an Earth-like Mg# (0.89) in the lunar mantle results in core densities (7800$\pm$100 kg/m$^3$) consistent with that of Fe-Ni alloy, whereas FeO-rich compositions (Mg# = 0.80--0.84) require lower densities (6100$\pm$800 kg/m$^3$). Geochemically, we use new data on mare basalts to reassess the bulk composition of the Moon for 70 elements, and show that the lunar core likely formed near 5 GPa, 2100 K and $\sim$1 log unit below the iron-wüstite buffer. Moreover, the Moon is depleted relative to the Earth's mantle in elements with volatilities higher than that of Li, with this volatile loss likely having occurred at low temperatures (1400$\pm$100 K), consistent with mass-dependent stable isotope fractionation of moderately volatile elements (e.g., Zn, K, Rb). The identical nucleosynthetic (O, Cr, Ti) and radiogenic (W) isotope compositions of the lunar and terrestrial mantles, strongly suggest the two bodies were made from the same material, rather than from an Earth-like impactor. Rb-Sr in FANs and Lu-Hf and Pb-Pb zircon ages point Moon formation close to $\sim$4500 Ma. Taken together, there is no unambiguous geochemical or isotopic evidence for the role of an impactor in the formation of the Moon, implying perfect equilibration between the proto-Earth and Moon-forming material or alternative scenarios for its genesis.

Mid-infrared imaging of exoplanets and disks is now possible with the coronagraphs of the MIRI on the JWST. This wavelength range unveils new features of young directly imaged systems and allows us to obtain new constraints for characterizing the atmosphere of young giant exoplanets and associated disks. These observations aim to characterize the atmosphere of the planet HD 95086 b by adding mid-infrared information so that the various hypotheses about its atmospheric parameters values can be unraveled. Improved images of circumstellar disks are provided. We present the MIRI coronagraphic imaging of the system HD 95086 obtained with the F1065C, F1140, and F2300C filters at central wavelengths of 10.575, 11.3, and 23 microns, respectively. We explored the method for subtracting the stellar diffraction pattern in the particular case when bright dust emitting at short separation is present. Furthermore, we compared different methods for extracting the photometry of the planet. Using the atmospheric models Exo-REM and ATMO, we measured the atmospheric parameters of HD 95086 b. The planet HD 95086 b and the contribution from the inner disk are detected at the two shortest MIRI wavelengths F1065C and F1140C. The outer colder belt is imaged at 23 microns. The mid-infrared photometry provides better constraints on the atmospheric parameters. We evaluate a temperature of 850-1020 K, consistent with one previous hypothesis that only used NIR data. The radius measurement of 1.0-1.13 RJup is better aligned with evolutionary models, but still smaller than predicted. These observations allow us to refute the hypothesis of a warm circumplanetary disk. HD 95086 is one of the first exoplanetary systems to be revealed at mid-infrared wavelengths. This highlights the interests and challenges of observations at these wavelengths.

One of the main obstacles for the signal extraction of the three point correlation function using photometric surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST), will be the prohibitive computation time required for dealing with a vast quantity of sources. Brute force algorithms, which naively scales as $\mathcal{O}(N^3)$ with the number of objects, can be further improved with tree methods but not enough to deal with large scale correlations of Rubin's data. However, a harmonic basis decomposition of these higher order statistics reduces the time dramatically, to scale as a two-point correlation function with the number of objects, so that the signal can be extracted in a reasonable amount of time. In this work, we aim to develop the framework to use these expansions within the Limber approximation for scalar (or spin-0) fields, such as galaxy counts, weak lensing convergence or aperture masses. We develop an estimator to extract the signal from catalogs and different phenomenological and theoretical models for its description. The latter includes halo model and standard perturbation theory, to which we add a simple effective field theory prescription based on the short range of non-locality of cosmic fields, significantly improving the agreement with simulated data. In parallel to the modeling of the signal, we develop a code that can efficiently calculate three points correlations of more than 200 million data points (a full sky simulation with Nside=4096) in $\sim$40 minutes on a single high-performance computing node, enabling a feasible analysis for the upcoming LSST data.

We present photometric redshifts for 1 341 559 galaxies from the Physics of the Accelerating Universe Survey (PAUS) over 50.38 ${\rm deg}^{2}$ of sky to $i_{\rm AB}=23$. Redshift estimation is performed using DEEPz, a deep-learning photometric redshift code. We analyse the photometric redshift precision when varying the photometric and spectroscopic samples. Furthermore, we examine observational and instrumental effects on the precision of the photometric redshifts, and we compare photometric redshift measurements with those obtained using a template method-fitting BCNz2. Finally, we examine the use of photometric redshifts in the identification of close galaxy pairs. We find that the combination of samples from W1+W3 in the training of DEEPz significantly enhances the precision of photometric redshifts. This also occurs when we recover narrow band fluxes using broad bands measurements. We show that DEEPz determines the redshifts of galaxies in the prevailing spectroscopic catalogue used in the training of DEEPz with greater precision. For the faintest galaxies ($i_{\rm AB}=21-23$), we find that DEEPz improves over BCNz2 both in terms of the precision (20-50 per cent smaller scatter) and in returning a smaller outlier fraction in two of the wide fields. The catalogues were tested for the identification of close galaxy pairs, showing that DEEPz is effective for the identification of close galaxy pairs for samples with $i_{\rm AB} < 22.5$ and redshift $0.2 < z < 0.6$. In addition, identifying close galaxy pairs common between DEEPz and BCNz2 is a promising approach to improving the purity of the catalogues of these systems.

We present an analysis of the cold gas phase in a low metallicity starburst generated in a high-resolution hydrodynamical simulation of a gas-rich dwarf galaxy merger as part of the GRIFFIN project. The simulations resolve (4 M$_\odot$ gas phase mass resolution, $\sim$ 0.1 pc spatial resolution) the multi-phase interstellar medium with a non-equilibrium chemical heating/cooling network at temperatures below $10^4$ K. Massive stars are sampled individually and interact with the ISM through the formation of HII regions and supernova explosions. In the extended starburst phase, the ISM is dominated by cold ($T_\mathrm{gas} < 300$ K) filamentary clouds with self-similar internal structures. The clouds have masses of $10^{2.6}$ - $10^{5.6}$ M$_\odot$ with a power law mass function, $dN/dM \propto M^\alpha$ with $\alpha = -1.78 (\pm 0.08)$. They also follow the Larson relations, in good agreement with observations. We trace the lifecycle of the cold clouds and find that they follow an exponential lifetime distribution and an e-folding time of $\sim$ 3.5 Myr. Clouds with peak masses below $10^4$ M$_\odot$ follow a power law relation with their average lifetime $\tau_\mathrm{life} \propto M^{0.3}_\mathrm{max}$ which flattens out for higher cloud masses at $ < 10$ Myr. A similar relation exists between cloud size at peak mass and lifetime. This simulation of the evolution of a realistic galactic cold cloud population supports the rapid formation and disruption of star-forming clouds by stellar radiation and supernovae on a timescale less than 10 Myr.

The entire classical cosmological history between two extreme de Sitter vacuum solutions is discussed based on Einstein's equations and non-equilibrium thermodynamics. The initial non-singular de Sitter state is characterised by a very high energy scale which is equal or smaller than the reduced Planck mass. It is structurally unstable and all the continuous created matter, energy and entropy of the material component comes from the irreversible flow powered by the primeval vacuum energy density. The analytical expression describing the running vacuum is obtained from the thermal approach. It opens a new perspective to solve the old puzzles and current observational challenges plaguing the cosmic concordance model driven by a rigid vacuum. Such a scenario is also modelled through a non-canonical scalar field. It is demonstrated that the resulting scalar field model is shown to be step by step a faithful analytical representation of the thermal running vacuum cosmology.

Blue large-amplitude pulsators (BLAPs) are a recently discovered group of hot pulsating stars whose evolutionary status remains uncertain. We study the pulsation characteristics of BLAPs for the two main mass scenarios, $ 0.3 - 0.4 \: \rm M_{\odot} $ and $ 0.7 - 1.1 \: \rm M_{\odot} $, and compare them with observations to find evidence for either scenario. We compute about a half million linear BLAP models using MESA-RSP and compare the linear pulsation periods with the observed ranges of BLAP stars. For the low-mass scenario, BLAPs are in a region of the HR diagram where the growth rates are positive for the fundamental mode, and the model periods correspond well to the observed ones, assuming the reported luminosities of $ \sim 200 \: \rm L_{\odot} $. For the high-mass scenario, pulsations in the first overtone dominate. Assuming a larger range of luminosities, high-mass models could also explain all BLAPs, including the high-gravity BLAPs. Furthermore, we provide the first seismically constrained mass estimate for the first double-mode BLAP star, OGLE-BLAP-030. We find that a linear model with a mass of $ 0.62 \: \rm M_{\odot} $ matches the reported parameters of this star exactly, placing it in between the two mass scenarios. We also derive new period relations based on our models and all available observed BLAPs, and we find that the derived relations also support the low-mass scenario.

We present high-resolution observations of a C4.1-class solar flare (SOL2023-05-03T20:53) in AR 13293 from the ViSP and VBI instruments at the DKIST. The fast cadence, good resolution, and high polarimetric sensitivity of ViSP data provide a unique opportunity to explore the photospheric magnetic fields before and during the flare. We infer the magnetic field vector in the photosphere from the Fe I 6302 line using Milne-Eddington inversions. Combined analysis of the inverted data and VBI images reveals the presence of two oppositely-polarity pores exhibiting rotational motion both prior to and throughout the flare event. Data-driven simulations further reveal a complex magnetic field topology above the rotating pores, including a null-point-like configuration. We observed a 30% relative change in the horizontal component ($\delta F_h$) of Lorentz force at the flare peak time and roughly no change in the radial component. We find that the changes in $\delta F_h$ are the most likely driver of the observed pore rotation. These findings collectively suggest that the back-reaction of magnetic field line reconfiguration in the corona may influence the magnetic morphology and rotation of pores in the photosphere on a significantly smaller scale.

The cosmic variance on the star formation history may lead to bias to the merger rate density estimation of binary neutron star (BNS) mergers by the compact binary population synthesis. In this paper, we take the advantage of the large boxsize of the Millennium Simulation combined with the semi-analytic galaxy formation model GABE, and the parameterized population binary star evolution (BSE) model to examine how much effect will the cosmic variance introduce on the estimation of merger rate density of BNS mergers. We find that for sub-box size of $100\rm Mpc$ and $200\rm Mpc$, the variance of merger rate density $\sigma_{\rm R}/\rm R$ at different redshift is about $23\%-35\%$ and $13\%-20\%$ respectively. On one hand, as for the variance of the detection rate on BNS mergers with current LIGO-Virgo-KAGRA (LVK) detector network, this value is very small $\lesssim 10\%$, which indicates ignoring the cosmic variance is reasonable for estimating the merger rate density from current LVK observation. On the other hand, with next-generation gravitational wave detectors, it is possible to localize BNS mergers within sub-boxes possessing length of $\rm 40 Mpc$ for source redshift $z_{s}<0.2$. In such a small box, the cosmic variance of the merger rate density is significant, i.e., the value of $\sigma_{\rm R}/\rm R$ is about $\sim 55\%$. This hints that estimating the merger rate density of BNS in different sky areas may provide useful information on the cosmic variance.

The periodic coronal rain and in-phase radiative intensity pulsations have been observed in multiple wavelengths in recent years. However, due to the lack of three-dimensional coronal magnetic fields and thermodynamic data in observations, it remains challenging to quantify the coronal heating rate that drives the mass cycles. In this work, based on the MURaM code, we conduct a three-dimensional radiative magnetohydrodynamic simulation spanning from the convective zone to the corona, where the solar atmosphere is heated self-consistently through dissipation resulting from magneto-convection. For the first time, we model the periodic coronal rain in an active region. With a high spatial resolution, the simulation well resembles the observational features across different extreme ultraviolet wavelengths. These include the realistic interweaving coronal loops, periodic coronal rain and periodic intensity pulsations, with two periods of 3.0~h and 3.7~h identified within one loop system. Moreover, the simulation allows for a detailed three-dimensional depiction of coronal rain on small scales, revealing adjacent shower-like rain clumps $\sim500$~km in width and showcasing their multi-thermal internal structures. We further reveal that these periodic variations essentially reflect the cyclic energy evolution of the coronal loop under thermal non-equilibrium state. Importantly, as the driver of the mass circulation, the self-consistent coronal heating rate is considerably complex in time and space, with hour-level variations in one order of magnitude, minute-level bursts, and varying asymmetry reaching ten times between footpoints. This provides an instructive template for the ad hoc heating function, and further enhances our understanding of the coronal heating process.

New observations are presented of millimetre line emissions of the circumstellar envelope (CSE) of AGB star $\chi$ Cygni using the recently upgraded NOEMA array. $\chi$ Cygni is an S-type Mira variable, at the border between oxygen-rich and carbon-rich, that has been observed for over 40 years to display features giving evidence for the strong role played by pulsation-associated shock waves in the generation of its wind. The new observations give evidence for a bright H$^{12}$CN(3-2) line emission confined to the very close neighbourhood of the star, however significantly more extended in 2024 than in 2023. Interpretation of such variability in terms of maser emission is considered but found to raise significant unanswered questions. Moreover, other unexpected features are observed in the very close neighbourhood of the star, including low Si$^{16}$O(6-5)/Si$^{17}$O(6-5), $^{28}$SiO(5-4)/$^{29}$SiO(5-4) and $^{12}$CO(2-1)/$^{13}$CO(2-1) line emission ratios. Possibly confirming the important role played by shocks, we discuss a measurement of the SiO(5-4)/SiO(6-5) emission ratio, evidence for a recent mass ejection particularly enhanced in the north-western red-shifted octant of the circumstellar envelope, leaving a depression of emission in its wake, and patterns of enhanced CO(2-1) line-emission suggesting an interpretation in terms of episodic outflows enhanced over solid angles associated with the surface of convective cells and on a time scale of a few decades. Unravelling the mechanisms underlying such newly observed features is very challenging and confirmation of the reported observations with improved sensitivity and angular resolution would be highly welcome. Observation of SiO maser emission in the ($\nu$=1,$J$=6-5) transition is reported for the first time.

How and where the first generation of inner solar system planetesimals formed remains poorly understood. Potential formation regions are the silicate condensation line and water-snowline of the solar protoplanetary disk. Whether the chemical compositions of these planetesimals align with accretion at the silicate condensation line (water-free and reduced) or water-snowline (water-bearing and oxidized) is, however, unknown. Here we use Fe/Ni and Fe/Co ratios of magmatic iron meteorites to quantify the oxidation states of the earliest planetesimals associated with non-carbonaceous (NC) and carbonaceous (CC) reservoirs, representing the inner and outer solar system, respectively. Our results show that the earliest NC planetesimals contained substantial amounts of oxidized Fe in their mantles (3-19 wt% FeO). In turn, we argue that this required the accretion of water-bearing materials into these NC planetesimals. The presence of substantial quantities of moderately and highly volatile elements in their parent cores is also inconsistent with their accretion at the silicate condensation line and favors instead their formation at or beyond the water-snowline. Similar oxidation states in the early-formed parent bodies of NC iron meteorites and those of NC achondrites and chondrites with diverse accretion ages suggests that the formation of oxidized planetesimals from water-bearing materials was widespread in the early history of the inner solar system.

We present the results of bicoherence analysis for Swift J1727.8-1613 during its 2023 outburst, using data from Insight-HXMT. Our analysis focused on observations with quasi-periodic oscillations (QPOs) of frequencies greater than 1 Hz, revealing that all of them belong to type C QPOs. We found a strong correlation between the QPO frequency and the hardness ratio, as well as a linear relationship between the QPO RMS and the hardness ratio. The bicoherence analysis revealed a transition from a "web" pattern to a "hypotenuse" pattern in the LE and HE energy bands. In the bicoherence patterns, there are correlations between horizontal and vertical bicoherence at $f_1=f_2=f_{\rm QPO}$ with count rates. The diagonal structure at $f_1+f_2=f_{\rm QPO}$ becomes more prominent with increasing energy. Additionally, we discovered a new bicoherence pattern in the medium energy band from 10 -- 20 keV, the diagonal structure at $f_1+f_2=f_{\rm har}$ is prominent only in this energy band, which we refer to as the "parallel" pattern. The bicoherence analysis indicates that the source is likely a low-inclination source.

As lensing of coherent astrophysical sources e.g. pulsars, fast radio bursts, and gravitational waves becomes observationally relevant, the mathematical framework of Picard-Lefschetz theory has recently been introduced to fully account for wave optics effects. Accordingly, the concept of lensing images has been generalized to include complex solutions of the lens equation referred to as "imaginary images", and more radically, to the Lefschetz thimbles which are a sum of steepest descent contours connecting the real and imaginary images in the complex domain. In this wave-optics-based theoretical framework of lensing, we study the "Stokes phenomena" as the change of the topology of the Lefschetz thimbles. Similar to the well-known caustics at which the number of geometric images changes abruptly, the corresponding Stokes lines are the boundaries in the parameter space where the number of effective imaginary images changes. We map the Stokes lines for a few lens models. The resulting Stokes line-caustics network represents a unique feature of the lens models. The observable signature of the Stokes phenomena is the change of interference behavior, in particular the onset of frequency oscillation for some Stokes lines. We also demonstrate high-order Stokes phenomena where the system has a continuous number of effective images but with an abrupt change in the way they are connected to each other by the Lefschetz thimbles. Their full characterization calls for an analogy of the catastrophe theory for caustics.

Most exoplanets with a few Earth radii are more inflated than bare-rock planets with the same mass, indicating a substantial volatile amount. Neither the origin of the volatiles nor the planet's bulk composition can be constrained from the mass-radius relation alone, and the spectral characterization of their atmospheres is needed to solve this degeneracy. Previous studies showed that chemical interaction between accreted volatile and possible molten rocky surface (i.e., magma ocean) can greatly affects the atmospheric composition. However, a variety in the atmospheric compositions of such planets with different properties remains elusive. In this work, we examine the dependence of atmospheric H, O, and C on planetary parameters (atmospheric thickness, planetary mass, equilibrium temperature, and magma properties such as redox state) assuming nebula gas accretion on an Earth-like core, using an atmosphere-magma chemical equilibrium model. Consistent with previous work, we show that atmospheric $\rm H_{2}O$ fraction on a fully molten rocky interior with an Earth-like redox state is on the order of $10^{-2}$-$10^{-1}$ regardless of other planetary parameters. Despite the solubility difference between H- and C-bearing species, C/H increases only a few times above the nebula value except for atmospheric pressure $\lesssim$1000 bar and $\rm H_{2}O$ fraction $\gtrsim$10\%. This results in a negative O/H-C/O trend and depleted C/O below one-tenth of the nebula gas value under an oxidized atmosphere, which could provide a piece of evidence of rocky interior and endogenic water. We also highlight the importance of constraints on the high-pressure material properties for interpreting the magma-atmospheric interaction of inflated planets.

We conducted the first light curve study of the binary star EL Tuc within the Binary Systems of South and North (BSN) Project framework. The photometric observations were made using standard multiband BVRI filters at an observatory in Argentina. We presented a new ephemeris for EL Tuc and a linear fit to the O-C diagram, utilizing our extracted times of minima and additional literature. We employed the PHysics Of Eclipsing BinariEs (PHOEBE) Python code and the Markov chain Monte Carlo (MCMC) approach for the system's light curve analysis. The target system's light curve solution required a cold starspot on the hotter component. We conclude that EL Tuc is a total contact binary system with a low mass ratio of q=0.172+_0.002, an orbital inclination of i=83.74+_0.40 degree, and a fillout factor of f=53.7+_1.6%. We used the P-a relationship and the Gaia Data Release 3 (DR3) parallax method to determine the absolute parameters of EL Tuc to compare the precision of our results. This system was classified as W-type based on the mass and effective temperature of the companion stars. The positions of the systems were depicted on the M-L, M-R, T-M, and q-L_ratio diagrams. The relationship between the spectroscopic and photometric mass ratios of binaries was discussed.

(abridged) Together with gravity, turbulence, and stellar feedback, magnetic fields (B-fields) are thought to play a critical role in the evolution of molecular clouds and star formation processes. We aim to map the morphology and measure the strength of B-fields of the nearby molecular cloud, rho Ophiuchus-A ($\rho$ Oph-A), and then to understand the role of B-fields in regulating star formation and shaping the cloud. We have analyzed the far-infrared (FIR) polarization of thermal dust emission observed by SOFIA/HAWC$+$ at 89 and 154 $\mu$m toward the densest part of $\rho$ Oph-A, which is irradiated by the nearby B3/4 star, Oph-S1. The cloud exhibits well-ordered B-fields with magnetic orientations mainly perpendicular to the ridge of the cloud toward the densest region and B-field strengths are in the range of 0.2-2.5 mG, using the Davis-Chandrasekhar-Fermi method. The B-fields are strongest at the densest part of the cloud, which is associated with the starless core SM1, and decreases toward the outskirts of the cloud. By calculating the map of the mass-to-flux ratio, Alfvén Mach number, and plasma $\beta$ parameter in $\rho$ Oph-A, we find that the cloud is predominantly magnetically sub-critical, sub-Alfvénic, which implies that the cloud is supported by strong B-fields that dominate over gravity, turbulence, and thermal gas energy. Measured B-field strengths at two densest subregions using other methods that account for the compressible mode are relatively lower than that measured with the DCF method but do not significantly change our conclusions on the roles of B-fields relative to gravity and turbulence on star formation. A virial analysis suggests that the cloud is gravitationally unbound. We find that B-fields are sufficiently strong to support the cloud against radiative feedback and to regulate the shape of the cloud.

We present the Python pipeline that was developed to simulate the high contrast mode images of the MICADO instrument. This mode will comprise three classical Lyot coronagraphs with different occulting spot sizes, one vector apodized phase plate, and two sparse aperture masks. One critical aspect of these modes lies in their sensitivity to aberrations, requiring careful integration of non common path aberrations, as well as turbulent adaptive optics (AO) residuals. Besides, they will be operated following a specific observing strategy based on pupil tracking mode. For these reasons, we have developed the MICADO SimulaTor for HIgh Contrast (MISTHIC) pipeline to simulate realistic image cubes and derive the expected performance of these modes. Several sources of aberrations can be included: turbulent AO residuals, but also static or rotating aberrations, amplitude aberrations, simple Zernike screens... Including MICADO's specific features, such as the absence of atmospheric dispersion prior to the coronagraphic focal plane mask, our intent is to make this pipeline available to the community. This tool can be used to prepare scientific observations with the high contrast mode of MICADO, by predicting the performance to be expected within the most realistic assumptions.

Star formation in the central molecular zone (CMZ) is suppressed with respect to that of the Galactic disk, and this is likely related to its high turbulent environment. This turbulence impedes the potential detection of prestellar cores. We present the temperature, density, and spatial structure of the CMZ molecular cloud G+0.693-0.027, which has been proposed to host a prestellar cluster in the Sgr B2 region. We analysed multiple HC$_{3}$N rotational transitions that were observed with the IRAM 30m, APEX, Yebes 40m, and GBT radio telescopes, together with SMA+APEX spatially resolved maps. The spectral shape of HC$_{3}$N lines shows three velocity components: a broad component with a line width of 23 km s$^{-1}$ (C1), and two narrow components with line widths of 7.2 and 8.8 km s$^{-1}$ (C2 and C3). This suggests that a fraction of the molecular gas in this cloud is undergoing turbulence dissipation. From a non-local thermodynamic equilibrium analysis we have found H$_{2}$ densities of 2$\times$10$^{4}$ cm$^{-3}$, 5$\times$10$^{4}$ cm$^{-3}$, and 4$\times$10$^{5}$ cm$^{-3}$ and kinetic temperatures of 140 K, 30 K, and 80 K for C1, C2, and C3, respectively. The spatially resolved maps confirm that the colder and high-density condensations C2 and C3, which peak in the 70-85 km s$^{-1}$ velocity range, are embedded in a more diffuse and warmer gas (C1). The larger-scale structure of the Sgr B2 region shows a hole at 40-50 km s$^{-1}$ that is likely due to a small cloud that shocked the Sgr B2 region and is spatially related with a massive cloud at 60-80 km s$^{-1}$. We propose that the impacting small cloud sequentially triggered the formation of Sgr B2(M), (N), and (S) and the condensations in G+0.693-0.027 during its passage. Based on the analysis of the masses of the two condensations and on the virial parameters, C2 might expand, while C3 might further fragment or collapse.

MICADO, the European Extremely Large Telescope first light imager will feature a dedicated high contrast imaging mode specifically designed for observing and characterizing exoplanets and circumstellar disks. Two sparse aperture masks (SAM) will be included, consisting of opaque masks with a set of holes arranged in a non-redundant configuration. Pupil masking transforms a monolithic telescope into an interferometer, with the aim of recovering spatial information down to the diffraction limit of the telescope and below, even in the presence of residual aberrations, such as turbulent AO residuals, and non common path aberrations. On the ELT, SAM will enable the detection of features down to 3.3mas in the J band, 12mas in the K band. Two designs have been chosen and will be reviewed, with complementarity in terms of sensitivity and spatial frequency coverage for image reconstruction. In this contribution, the technical choices will be detailed, such as the hole diameter and arrangement, given the technical constraints such as spectral filter bandwidths, detector dimensions, sampling, read-out-noise and frame rate. We will also report on simulations performed to assess the expected capabilities of this mode, with application examples of close companion detection and contrast curves.

Solar active regions (ARs) are often surrounded by dark large areas of reduced emission compared to the quiet Sun, observed at various wavelengths corresponding to chromosphere, transition region (TR) and corona, and known as Dark Halos (DHs). DHs have been insufficiently studied, and the mechanisms behind their darker emission remain unclear. This study aims to investigate for the first time the fine structure of a DH observed by the EUV High Resolution Imager (HRI$_{EUV}$) onboard the ESA's Solar Orbiter (SO) mission and its appearance in the TR. We utilized the extensive 1-hour dataset from SO on 19 March 2022, which includes high-resolution observations of NOAA 12967 and part of the surrounding DH. We analyzed the dynamics of the HRI$_{EUV}$ DH fine structure and its appearance in the HRI$_{Ly\alpha}$ image and the Spectral Imaging of the Coronal Environment (SPICE) Ly$\beta$, C III, N VI, O VI and Ne VIII lines, which sample the TR in the logT (K) $\sim$ 4.0 - 5.8 range. This analysis was complemented with a simultaneous B$_{LOS}$ magnetogram taken by the High Resolution Telescope (HRT). We report the presence of a peculiar fine structure which is not observed in the quiet Sun, characterized by combined bright EUV bundles and dark regions, arranged and interconnected in such a way that they cannot be clearly separated. They form a spatial continuum extending approximately radially from the AR core, suggesting a deep connection between the DH and the AR. Additionally, we find that the bright EUV bundles are observed in all the SPICE TR lines and the HRI$_{Ly\alpha}$ band and present photospheric B$_{LOS}$ footprints in the HRT magnetogram. This spatial correlation indicates that the origin of the 174 Å DH may lie in the low atmosphere, i.e. photosphere/chromosphere.

We report the discovery and characterization of a transiting sub-Saturn exoplanet TOI-6651b using PARAS-2 spectroscopic observations. The host, TOI-6651 ($m_{V}\approx 10.2$), is a sub-giant, metal-rich G-type star with $[{\rm Fe/H}] = 0.225^{+0.044}_{-0.045}$, $T_{\rm eff} = 5940\pm110\ \mathrm{K}$, and $\log{g} = 4.087^{+0.035}_{-0.032}$. Joint fitting of the radial velocities from PARAS-2 spectrograph and transit photometric data from Transiting Exoplanet Survey Satellite (TESS) reveals a planetary mass of $61.0^{+7.6}_{-7.9}\ M_\oplus$ and radius of $5.09^{+0.27}_{-0.26}\ R_\oplus$, in a $5.056973^{+0.000016}_{-0.000018}$ day orbit with an eccentricity of $0.091^{+0.096}_{-0.062}$. TOI-6651b has a bulk density of $2.52^{+0.52}_{-0.44}\ \mathrm{g\ cm^{-3}}$, positioning it among the select few known dense sub-Saturns and making it notably the densest detected with TESS. TOI-6651b is consistent with the positive correlation between planet mass and the host star's metallicity. We find that a considerable portion $\approx$ 87% of the planet's mass consists of dense materials such as rock and iron in the core, while the remaining mass comprises a low-density envelope of H/He. TOI-6651b lies at the edge of the Neptunian desert, which will be crucial for understanding the factors shaping the desert boundaries. The existence of TOI-6651b challenges conventional planet formation theories and could be a result of merging events or significant atmospheric mass loss through tidal heating, highlighting the complex interplay of dynamical processes and atmospheric evolution in the formation of massive dense sub-Saturns.

We investigated the effect of time-dependent ice growth on dust grains on the opacity and hence on the dust temperature in a collapsing molecular cloud core, with the aim of quantifying the effect of the dust temperature variations on ice abundances as well as the evolution of the collapse. We employed a one-dimensional collapse model that self-consistently and time-dependently combines hydrodynamics with chemical and radiative transfer simulations. The dust opacity was updated on-the-fly based on the ice growth as a function of location in the core. The results of the fully dynamical model were compared against simulations assuming fixed ice thickness. We found that the ice thickness increases fast and reaches a saturation value of approximately 90 monolayers in the central core (volume density $\sim$$10^4\,\rm cm^{-3}$), and several tens of monolayers at a volume density of $\sim$$10^3\,\rm cm^{-3}$, after only a few $10^5\,\rm yr$ of evolution. The results thus exclude the adoption of thin ($\sim$10 monolayer) ices in molecular cloud simulations except at very short timescales. The differences in abundances and dust temperature between the fully dynamic simulation and those with fixed dust opacity are small; abundances change between the solutions generally within a factor of two. The assumptions on the dust opacity do have an effect on the collapse dynamics through the influence of the photoelectric effect on the gas temperature, and the simulations take a different time to reach a common central density. In conclusion, carrying out chemical simulations using a dust temperature corresponding to a fixed opacity seems to be a good approximation. Still, although at least in the present case its effect on the overall results is limited - as long as the grains are monodisperse - ice growth should be considered to obtain the most accurate representation of the collapse dynamics.

Galaxy cluster counts have historically been important for the measurement of cosmological parameters and upcoming surveys will greatly reduce the statistical errors. To exploit the potential of current and future cluster surveys, theoretical uncertainties on the predicted abundance must be smaller than the statistical errors. Models used to predict cluster counts typically combine a model for the dark matter only (DMO) halo mass function (HMF) with an observable - mass relation that is assumed to be a power-law with lognormal scatter. We use the FLAMINGO suite of cosmological hydrodynamical simulations to quantify the biases in the cluster counts and cosmological parameters resulting from the different ingredients of conventional models. For the observable mass proxy we focus on the Compton-Y parameter quantifying the thermal Sunyaev-Zel'dovich effect, which is expected to result in cluster samples that are relatively close to mass-selected samples. We construct three mock surveys based on existing (Planck and SPT) and upcoming (Simons Observatory) surveys. We ignore measurement uncertainties and compare the biases in the counts and inferred cosmological parameters to each survey's Poisson errors. We find that widely used models for the DMO HMF differ significantly from each other and from the DMO version of FLAMINGO, leading to significant biases for all three surveys. For upcoming surveys, dramatic improvements are needed for all additional model ingredients, i.e. the functional forms of the fits to the observable-mass scaling relation and the associated scatter, the priors on the scaling relation and the prior on baryonic effects associated with feedback processes on the HMF.

Large-scale cosmic-ray detectors like the Giant Radio Array for Neutrino Detection (GRAND) are pushing the boundaries of our ability to identify air shower events. Existing trigger schemes rely solely on the timing of signals detected by individual antennas, which brings many challenges in distinguishing true air shower signals from background. This work explores novel event-level radio trigger methods specifically designed for GRAND, but also applicable to other systems, such as the Radio Detector (RD) of the Pierre Auger Observatory. In addition to an upgraded plane wave front reconstruction technique, we introduce orthogonal and complementary approaches that analyze the radio-emission footprint, the spatial distribution of signal strength across triggered antennas, to refine event selection. We test our methods on mock data sets constructed with simulated showers and real background noise measured with the GRAND prototype, to assess the performance potential in terms of sensitivity and background rejection in GRAND. Our preliminary results are a first step to identifying the most discriminating radio signal features at event-level, and optimizing the techniques for future implementation on experimental data.

Context: In the context of Galactic archaeology, the outer halo remains relatively unexplored with respect to its metallicity distribution, merger debris, and the abundance of known very/extremely metal-poor ([Fe/H]<-2.5) stars. Aim: We utilize the Pristine survey's publicly available, Pristine data release 1 (PDR1) and Pristine-Gaia synthetic (PGS) catalogues of photometric metallicities, to select Red Giant Branch (RGB) stars in the outer halo. Methods: The RGB selection pipeline selects giants based on the absence of a well-measured parallax in the brightness range where dwarfs have reasonable parallax estimate from Gaia DR3 data. The photometric distances are calculated using a BaSTI-isochrone fitting code and the photometric metallicities. Results: Photometric distances derived from PDR1- and PGS-giants show typical uncertainties of 12% and a scatter of up to 20% and 40% respectively, when validated against inverted-parallax and Starhorse-code distances. The PDR1-giants catalogue provides a low-to-no bias view of the metallicity structure versus distance compared to the PGS-giants catalogue (with a distance-metallicity selection bias), while the PGS-giants catalogue provides an all-sky view of the outer halo, especially in the VMP end. The PDR1-giants catalogue is used to study the metallicity distribution function of the halo out to ~100 kpc, using Gaussian Mixture Model decomposition. Additionally, we use the PDR1-giants with radial velocities from spectroscopic surveys to map the metallicity view of the integrals-of-motion space where accreted dwarf galaxy debris conserve their orbital parameters for a long time. Using the PGS-giants catalogue, we associate 41 stars tentatively to the stellar counterpart of the Magellanic stream in the VMP end, out to 70 kpc. Conclusions: We publish two RGB catalogues between -4.0<[Fe/H]<+0.1 with reliable photometric distances and metallicities.

Solar jets are collimated plasma ejections along magnetic field lines observed in hot (EUV jets) and cool (chromospheric surges) temperature diagnostics. Their trigger mechanisms and the relationship between hot and cool jets are still not completely understood. We aim to investigate the generation of a sequence of active region solar jets and their evolution from the photospheric to the coronal heights. Using the synergy of high spatial and temporal resolution observations by the SST, along with the SDO, we analyze a sequence of solar jets originating in a mixed polarity region between the leading and following sunspots of an active region. We use a NFFF extrapolation technique for deriving the magnetic field topology of the active region. A mixed polarity region is formed over a long period (24 hours) with persistent magnetic flux emergence. This region has been observed as an arch filament system (AFS) in chromospheric SST observations. In this region, negative polarities surrounded by positive polarities create a fan-surface with a null point at a height of 6 Mm detected in the NFFF extrapolation. SST observations in H-beta spectral line reveal a large flux rope over the AFS and moving from the North to South, causing successive EUV and cool jets to move in the East-West direction and later towards the South along the long open loops. The high resolution SST observations (0.038 arcsec per pixel) resolve the dark area observed at the jet base and reveal the existence of an AFS with an extended cool jet which may be the result of a peeling-like mechanism of the AFS. Based on the combined analysis of SST and AIA observations along with extrapolated magnetic topology, it is suggested that the magnetic reconnection site may move southward by approximately 20 Mm until it reaches a region where the open magnetic field lines are oriented North-South.

Gravitational waves (GWs) can undoubtedly serve as a messenger from the early Universe acting as well as a novel probe of the underlying gravity theory. In this work, motivated by one-loop vacuum-polarization effects on curved spacetime, we investigate a gravitational theory with non-minimal curvature-electromagnetic coupling terms of the form $\xi R F_{\mu\nu}F^{\mu\nu}$, where $R$ is the scalar curvature and $F_{\mu\nu}$ the Faraday tensor, which can be responsible for the generation of primordial electromagnetic fields. We study then the GW signatures of such coupling terms deriving in particular for the first time to the best of our knowledge the modified tensor modes equation of motion. Remarkably, we find a universal infrared (IR) frequency scaling $f^5$ of the electromagnetically induced GW (EMIGW) signal, which, depending on the energy scale of inflation, the duration of inflation and reheating as well as the dynamical behaviour of the gauge coupling function $\xi$, can be well within the detection sensitivity bands of GW experiments such as SKA, LISA, ET and BBO, being thus potentially detectable in the future by GW observatories.

In the era of massive photometric surveys, we explore several approaches to estimate the masses of supermassive black holes (SMBHs) in active galactic nuclei (AGN) from optical ground-based imaging, in each case comparing to the independent SMBH mass measurement obtained from spectroscopic data. We select a case-study sample of 28 type 1 AGN hosted by nearby galaxies from the Galaxy And Mass Assembly (GAMA) survey. We perform multi-component spectral decomposition, extract the AGN component and calculate the SMBH mass from the broad H$\alpha$ emission line width and luminosity. The photometric $g$ and $i$ band data is decomposed into AGN+spheroid(+disc)(+bar) components with careful surface brightness fitting. From these, the SMBH mass is estimated using its relation with the spheroid Sérsic index or effective radius (both used for the first time on ground-based optical imaging of AGN); and the more widely used scaling relations based on bulge or galaxy stellar mass. We find no correlation between the H$\alpha$-derived SMBH masses and those based on the spheroid Sérsic index or effective radius, despite these being the most direct methods involving only one scaling relation. The bulge or galaxy stellar mass based methods both yield significant correlations, although with considerable scatter and, in the latter case, a systematic offset. We provide possible explanations for this and discuss the requirements, advantages and drawbacks of each method. These considerations will be useful to optimise stategies for upcoming high quality ground-based and space-borne sky surveys to estimate SMBH masses in large numbers of AGN.

A significant barrier to achieving the radial velocity (RV) measurement accuracy and precision required to characterize terrestrial mass exoplanets is the existence of time-variable features in the measured spectra, from both telluric absorption and stellar variability, which affect measured line shapes and can cause apparent RV shifts. Reaching the desired accuracy using traditional techniques often requires avoiding lines contaminated by stellar variability and/or changing tellurics, and thus discarding a large fraction of the spectrum, lowering precision. New data-driven methods can help achieve extremely precise and accurate RVs by enabling the use of a larger fraction of the available data. While there exist methods for modeling telluric features or the stellar variability individually, there is a need for additional tools that are capable of modeling them simultaneously at the spectral level. Here we present StellarSpectraObservationFitting.jl (SSOF), a Julia package for measuring Doppler shifts and creating data-driven models (with fast, physically-motivated Gaussian Process regularization) for the time-variable spectral features for both the telluric transmission and stellar spectrum, while accounting for the wavelength-dependent instrumental line-spread function. We demonstrate SSOF's state-of-the-art performance on data from the NEID RV spectrograph on the WIYN 3.5m Telescope for multiple stars. We show SSOF's, ability to accurately identify and characterize spectral variability and provide $\sim$2-6x smaller photon-limited errors over the NEID CCF-based pipeline and match the performance of SERVAL, a leading template-based pipeline, using only observed EPRV spectra.

We implement a standard thin disk model with the outer disk radius ($R_{\rm out}$) as a free parameter, integrating it into standard X-ray fitting package to enable self-consistent and simultaneous fitting of X-ray spectra and UV/optical/NIR photometry. We apply the model to the late-time data ($\Delta t \approx 350-1300$ days) of the tidal disruption event (TDE) ASASSN-14li. We show that at these late-times the multi-wavelength emission of the source can be fully described by a bare compact accretion disk. We obtain a black hole mass ($M_{\rm BH}$) of $7^{+3}_{-2}\times10^{6} M_{\odot}$, consistent with host-galaxy scaling relations; and an $R_{\rm out}$ of $45 \pm 13 \, R_{\rm g}$, consistent with the circularization radius, with possible expansion at the latest epoch. We discuss how simplistic models, such as a single-temperature blackbody fitted to either X-ray spectra or UV/optical photometry, lead to erroneous interpretations on the scale/energetics of TDE emission. We also apply the model to the soft/high state of the intermediate-mass black hole (IMBH) candidate HLX-1. The model fits the full spectral energy distribution (from X-rays to NIR) without needing an additional stellar population component. We investigate how relativistic effects improve our results by implementing a version of the model with full ray tracing calculations in the Kerr metric. For HLX-1, we find $M_{\rm BH} = 4^{+3}_{-1} \times 10^{4} M_{\odot}$ and $R_{\rm out} \approx {\rm few} \times 10^{3} \, R_{\rm g}$, in agreement with previous findings. The relativistic model can constrain the inclination ($i$) of HLX-1 to be $20^o \leq i \leq 70^o$.

New constraints on the expansion rate of the Universe seem to favor evolving dark energy in the form of thawing quintessence models, i.e., models for which a canonical, minimally coupled scalar field has, at late times, begun to evolve away from potential energy domination. We scrutinize the evidence for thawing quintessence by exploring what it predicts for the equation of state. We show that, in terms of the usual Chevalier-Polarski-Linder parameters, ($w_0$, $w_a$), thawing quintessence is, in fact, only marginally consistent with a compilation of the current data. Despite this, we embrace the possibility that thawing quintessence is dark energy and find constraints on the microphysics of this scenario. We do so in terms of the effective mass $m^2$ and energy scale $V_0$ of the scalar field potential. We are particularly careful to enforce un-informative, flat priors on these parameters so as to minimize their effect on the final posteriors. While the current data favors a large and negative value of $m^2$, when we compare these models to the standard $\Lambda$CDM model we find that there is scant evidence for thawing quintessence.

Contamination from singly ionized iron emission is one of the greatest obstacles to determining the intensity of emission lines in the UV and optical wavelength ranges. This study presents a comprehensive analysis of the FeII emission in the bright quasar RM 102, based on the most recent version of the CLOUDY software, with the goal of simultaneously reproducing UV and optical FeII emission. We employ a constant pressure model for the emitting clouds, instead of the customary constant density assumption. The allowed parameter range is broad, with metallicity up to 50 times the solar value and turbulent velocity up to 100 km s$^{-1}$ for a subset of models. We also consider geometrical effects that could enhance the visibility of the non-illuminated faces of the clouds, as well as additional mechanical heating. Our investigation reveals that the broad line region of RM 102 is characterized by highly metallic gas. The observed FeII features provide strong evidence for an inflow pattern geometry that favours the dark sides of clouds over isotropic emission. This study confirms the presence of chemically enriched gas in the broad line region of bright quasars, represented by RM 102, which is necessary to explain the strong FeII emission and its characteristic features. Additionally, we report that CLOUDY currently still lacks certain transitions in its atomic databases which prevents it from fully reproducing some observed FeII features in quasar spectra.

The Giant Radio Array for Neutrino Detection (GRAND) is a distributed, sparse ground antenna array designed to detect the radio emission from highly inclined extensive air showers induced by ultra-high-energy particles in the atmosphere. We use CoREAS air-shower simulations to adapt an existing signal model of the radio emission of inclined showers to the 50-200 MHz frequency band GRAND is sensitive to. The model uses a parameterisation of the charge excess emission to isolate the geomagnetic component. By fitting a one-dimensional lateral distribution function to the geomagnetic energy fluence of a shower, we reconstruct its radiation energy. This work details the signal model and the intrinsic event reconstruction of our method, as well as the adaptations to the new frequency band. This work is part of the NUTRIG project.

Galaxy evolution depends on the environment in which galaxies are located. The various physical processes (ram-pressure stripping, tidal interactions, etc.) that can affect the gas content in galaxies have different efficiencies in different environments. In this work, we examine the gas (atomic \ce{HI} and molecular \ce{H2}) content of local galaxies inside and outside clusters, groups, and filaments as well as in isolation using observational and simulated data. We exploited a catalog of galaxies in the Virgo cluster(including the surrounding filaments) and compared the data against the predictions of the Galaxy Evolution and Assembly(GAEA) semi-analytic model, which has explicit prescriptions for partitioning the cold gas content in its atomic and molecular phases. We extracted from the model a mock catalog that mimics the observational biases and one not tailored to observations to study the impact of observational limits on the results and predict trends in regimes not covered by the current observations. The observations and simulated data show that galaxies within filaments exhibit intermediate cold gas content between galaxies in clusters and isolation. The amount of \ce{HI} is typically more sensitive to the environment than \ce{H2} and low-mass galaxies ($\log_{10} [{\rm M}_{\star} / \rm{M}_{\sun} ] < 10$) are typically more affected than their massive ($\log_{10} [{\rm M}_{\star} / \rm{M}_{\sun} ] > 10$) counterparts. Considering only model data, we identified two distinct populations among filament galaxies present in similar proportions: those simultaneously lying in groups and isolated galaxies. The former has properties more similar to cluster and group galaxies, and the latter is more similar to those of field galaxies. We therefore did not detect filaments' strong effects on galaxies' gas content, and we ascribe the results to the presence of groups in filaments.

Using the light curve time series data for more than 11.7 million variable sources published in the Gaia Data Release 3, the average magnitudes, colors, and variability parameters have been computed for 0.836 million Gaia CRF objects, which are mostly quasars and active galactic nuclei (AGNs). To mitigate the effects of occasional flukes in the data, robust statistical measures have been employed, namely, the median, median absolute deviation, and Spearman correlation. We find that the majority of the CRF sources have moderate amplitudes of variability in the Gaia $G$ band just below 0.1 mag. The heavy-tailed distribution of variability amplitudes (quantified as robust standard deviations) does not find a single analytical form, but is closer to Maxwell distribution with a scale of 0.078 mag. The majority of CRF sources have positive correlations between $G$ magnitude and $G_{\rm BP}-G_{\rm RP}$ colors, meaning that these quasars and AGNs become bluer when they are brighter. The variations in the BP and RP bands are also mostly positively correlated. Dependencies of all variability parameters with cosmological redshift are fairly flat for the more accurate estimates above redshift 0.7, while the median color shows strong systematic variations with redshift. Using a robust normalized score of magnitude deviations, a sample of 5000 most variable quasars is selected and published. The intersection of this sample with the ICRF3 catalog shows a much higher rate of strongly variable quasars (mostly, blazars) in ICRF3.

The GAia Map of the Brightness Of the Natural Sky (GAMBONS) is a model to map the natural night brightness of the sky in cloudless and moonless nights. It computes the star radiance from the photometric data in Gaia and Hipparcos catalogues, adding the contributions of the diffuse galactic and extragalactic light, zodiacal light and airglow, and taking into account the effects of atmospheric attenuation and scattering. The model allows computing the natural sky brightness in any given photometric band for a ground-based observer, if appropriate transformations from the Gaia bands are available. In this work we present the most recent improvements of the model, including the use of Gaia EDR3 data, the inclusion of a wide set of photometric bands and the derivation of additional sky brightness indicators, as the horizontal irradiance and the average hemispheric radiance

Context. Recent evidences from spectroscopic surveys point towards the presence of a metal-poor, young stellar population in the chemical thin disk. In this context, the investigation of the spatial distribution and time evolution of precise, unbiased abundances is fundamental to disentangle the scenarios of evolution of the Galaxy. Aims. We study the evolution of abundance gradients in the Milky Way by taking advantage of a large sample of open clusters, which are among the best tracers for this purpose. In particular, we use data from the last release of the Gaia-ESO survey. Methods. We perform careful selection of open cluster stars excluding those that may be affected by biases in spectral analysis. The cleaned open clusters sample is compared with detailed chemical evolution models for the Milky Way, using well tested stellar yields and prescription for radial migration. Different scenarios of Galaxy evolution are tested to explain the data, i.e. the two-infall and the three-infall frameworks, suggesting that the chemical thin disk is formed by one or two subsequent gas accretion episodes, respectively. Results. With the performed selection in cluster stars, we still find a metallicity decrease between intermediate age (1<Age/Gyr<3) and young (Age<1Gyr) open clusters. This decrease cannot be explained in the context of the two-infall scenario, even by accounting for the effect of migration and yield prescriptions. The three-infall model, with its late gas accretion in the last 3 Gyr, can explain the low metallic content in young clusters. However, we invoke a milder metal dilution relative to previous findings. Conclusions. To explain the observed low metallic content in young clusters, we propose that a late gas accretion episode triggering metal dilution should have taken place, extending the framework of the three-infall model for the first time to the entire Galactic disk.

Over the past three decades, multiple reverberation mapping (RM) campaigns conducted for the quasar PG 2130+099 have exhibited inconsistent findings with time delays ranging from $\sim$10 to $\sim$200 days. To achieve a comprehensive understanding of the geometry and dynamics of the broad-line region (BLR) in PG 2130+099, we continued an ongoing high-cadence RM monitoring campaign using the Calar Alto Observatory 2.2m optical telescope for an extra four years from 2019 to 2022. We measured the time lags of several broad emission lines (including He II, He I, H$\beta$, and Fe II) with respect to the 5100 Å continuum, and their time lags continuously vary through the years. Especially, the H$\beta$ time lags exhibited approximately a factor of two increase in the last two years. Additionally, the velocity-resolved time delays of the broad H$\beta$ emission line reveal a back-and-forth change between signs of virial motion and inflow in the BLR. The combination of negligible ($\sim$10%) continuum change and substantial time-lag variation (over two times) results in significant scatter in the intrinsic $R_{\rm H\beta}-L_{\rm 5100}$ relationship for PG 2130+099. Taking into account the consistent changes in the continuum variability time scale and the size of the BLR, we tentatively propose that the changes in the measurement of the BLR size may be affected by 'geometric dilution'.

We describe a novel proposal for inflationary magnetogenesis by identifying the non-Abelian sector of Spectator Chromo Natural Inflation (SCNI) with the $\rm{SU(2)}_{\rm L}$ sector of the Standard Model. This mechanism relies on the recently discovered attractor of SCNI in the strong backreaction regime, where the gauge fields do not decay on super-horizon scales and their backreaction leads to a stable new trajectory for the rolling axion field. The large super-horizon gauge fields are partly transformed after the electroweak phase transition into electromagnetic fields. The strength and correlation length of the resulting helical magnetic fields depend on the inflationary Hubble scale and the details of the SCNI sector. For suitable parameter choices we show that the strength of the resulting magnetic fields having correlation lengths around $1\, {\rm {Mpc}}$ are consistent with the required intergalactic magnetic fields for explaining the spectra of high energy $\gamma$ rays from distant blazars.

Context. To date, three nuclear transients have been associated with high-energy neutrino events. These transients are generally thought to be powered by tidal disruptions of stars (TDEs) by massive black holes. However, AT2019aalc, hosted in a Seyfert-1 galaxy, was not yet classified due to a lack of multiwavelength observations. Interestingly, the source has re-brightened 4 years after its discovery. Aims. We aim to classify the transient and explain the mechanism responsible for its second optical flare. Methods. We conducted a multi-wavelength monitoring program (from radio to X-rays) of AT2019aalc during its re-brightening in 2023. Results. The observations revealed a uniquely bright UV counterpart and multiple X-ray flares during the second optical flaring episode of the transient. The second flare, similarly to the first one, is also accompanied by IR dust echo emission. A long-term radio flare is found with an inverted spectrum. Optical spectroscopic observations reveal the presence of Bowen Fluorescence lines and strong high-ionization coronal lines indicating an extreme level of ionization in the system. Conclusions. The results suggest that the transient can be classified as a Bowen Fluorescence Flare (BFF), a relatively new sub-class of flaring active galactic nuclei (AGN). AT2019aalc can be also classified as an extreme coronal line emitter (ECLE). We found that, in addition to AT2019aalc, another BFF AT2021loi is spatially coincident with a high-energy neutrino event. The multi-wavelength properties of these transients suggest a possible connection between ECLEs, BFFs and TDEs in AGN.

Understanding space weather is vital for the protection of our terrestrial and space infrastructure. In order to predict space weather accurately, large amounts of data are required, particularly in the extreme ultraviolet (EUV) spectrum. An exquisite source of information for such data is provided by the Solar Dynamic Observatory (SDO), which has been gathering solar measurements for the past 13 years. However, after a malfunction in 2014 affecting the onboard Multiple EUV Grating Spectrograph A (MEGS-A) instrument, the scientific output in terms of EUV measurements has been significantly degraded. Building upon existing research, we propose to utilize deep learning for the virtualization of the defective instrument. Our architecture features a linear component and a convolutional neural network (CNN) -- with EfficientNet as a backbone. The architecture utilizes as input grayscale images of the Sun at multiple frequencies -- provided by the Atmospheric Imaging Assembly (AIA) -- as well as solar magnetograms produced by the Helioseismic and Magnetic Imager (HMI). Our findings highlight how AIA data are all that is needed for accurate predictions of solar irradiance. Additionally, our model constitutes an improvement with respect to the state-of-the-art in the field, further promoting the idea of deep learning as a viable option for the virtualization of scientific instruments.

High-redshift ($z\sim2-3$) galaxies accrete circumgalactic gas through cold streams. Recent high-resolution MHD simulations of these streams showed a significant amplification of the intergalactic magnetic field in the shear layer around them. In this work we estimate the magnetisation of high-redshift galaxies that would result purely due to the accretion of already magnetised gas from cold streams. We use the mass inflow rates and saturated magnetic field values from cold stream simulations as input to a simple analytic model that calculates the galactic magnetic field purely from mass accretion. Our model predicts average magnetic field strengths that exceed $\rm\mu G$ values at $z\sim 2-3$ for inflow rates above $0.1 \, \rm{M_{\odot} yr^{-1}}$. For high inflow rates, our model results are consistent with the recent detection of a strong magnetic field in a $z=2.6$ galaxy by Geach et al. (2023). Within the assumptions of our simple model, magnetised cold streams emerge as a viable mechanism for seeding a dynamically important galactic magnetic field.

We introduce picasso, a model designed to predict thermodynamic properties of the intracluster medium based on the properties of halos in gravity-only simulations. The predictions result from the combination of an analytical gas model, mapping gas properties to the gravitational potential, and of a machine learning model to predict the model parameters for individual halos based on their scalar properties, such as mass and concentration. Once trained, the model can be applied to make predictions for arbitrary potential distributions, allowing its use with flexible inputs such as N-body particle distributions or radial profiles. We present the model, and train it using pairs of gravity-only and hydrodynamic simulations. We show that when trained on non-radiative hydrodynamic simulations, picasso can make remarkably accurate and precise predictions of intracluster gas thermodynamics. Training the model on full-physics simulations yields robust predictions as well, albeit with slightly degraded performance. We further show that the model can be trained to make accurate predictions from very minimal information, at the cost of modestly reduced precision. picasso is made publicly available as a Python package, which includes trained models that can be used to make predictions easily and efficiently, in a fully auto-differentiable and hardware-accelerated framework.