Complex organic molecules (COMs) are thought to be the precursors of pre-biotic molecules and are observed in many protostellar sources. For this paper we studied the formation of COMs during star formation and their evolution in the midplane of the circumstellar disk up to the end of the Class I stage. We used the Analytical Protostellar Environment (APE) code to perform analytical simulations of star formation and the Nautilus code to model the chemical evolution. Most COMs mainly form during the collapse or in the disk, except the lightest (CH3CCH, C3H6, CH3OH, CH3CHO, CH3OCH3, C2H5OH, CH3CN, CH3NC, C2H3CN, and CH3SH), which are significantly inherited by the disk from the prestellar phase. Over the first 150 kyr of the disk, the abundances of several COMs in the midplane vary negligibly (e.g., CH3CCH, CH3OH, and CH3CN), while others experience a variation of one order of magnitude (e.g., C2H3CHO HOCH2CHO, and CH3COCH2OH). Changing physical conditions also have an impact on the abundance profiles of COMs in the disk, and their inheritance. For example, increasing the temperature of the molecular cloud from 10 K to 15 K significantly promotes the formation of COMs in the prestellar phase, notably c-C2H4O and N-bearing species. Conversely, increasing the cloud mass from 2 Msol to 5 Msol only has a minor effect on the disk abundances in the early stages.

This is the second paper in the HOWLS (higher-order weak lensing statistics) series exploring the usage of non-Gaussian statistics for cosmology inference within \textit{Euclid}. With respect to our first paper, we develop a full tomographic analysis based on realistic photometric redshifts which allows us to derive Fisher forecasts in the ($\sigma_8$, $w_0$) plane for a \textit{Euclid}-like data release 1 (DR1) setup. We find that the 5 higher-order statistics (HOSs) that satisfy the Gaussian likelihood assumption of the Fisher formalism (1-point probability distribution function, $\ell$1-norm, peak counts, Minkowski functionals, and Betti numbers) each outperform the shear 2-point correlation functions by a factor $2.5$ on the $w_0$ forecasts, with only marginal improvement when used in combination with 2-point estimators, suggesting that every HOS is able to retrieve both the non-Gaussian and Gaussian information of the matter density field. The similar performance of the different estimators\inlinecomment{, with a slight preference for Minkowski functionals and 1-point probability distribution function,} is explained by a homogeneous use of multi-scale and tomographic information, optimized to lower computational costs. These results hold for the $3$ mass mapping techniques of the \textit{Euclid} pipeline: aperture mass, Kaiser--Squires, and Kaiser--Squires plus, and are unaffected by the application of realistic star masks. Finally, we explore the use of HOSs with the Bernardeau--Nishimichi--Taruya (BNT) nulling scheme approach, finding promising results towards applying physical scale cuts to HOSs.

Recent JWST observations of the temperate sub-Neptune K2-18 b with NIRISS SOSS/NIRSpec G395H and MIRI LRS have yielded apparently inconsistent results: the MIRI spectra exhibit spectral features nearly twice as large as those seen at shorter wavelengths, challenging the high-metallicity, CH4-rich non-equilibrium model that fits the NIRISS/NIRSpec data. We perform a suite of atmospheric retrievals on both datasets, including free-chemistry, non-equilibrium, and aerosol models, using laboratory-derived complex refractive indices for a variety of photochemical haze analogues. Free retrievals systematically return lower metallicities than inferred by self-consistent chemical disequilibrium models, and the inclusion of absorbing aerosols, especially CH4-dominated, nitrogen-poor tholins, can further reduce the inferred metallicity by over an order of magnitude. These hazes reproduce the observed NIRISS slope through scattering and match MIRI features via C-H bending absorption near 7 um, while yielding particle properties consistent with photochemical production in H2-rich atmospheres. Although their inclusion improves the joint fit and reduces tension between datasets, it also significantly lowers the retrieved CH4 abundance, highlighting degeneracies between metallicity, composition, and aerosol properties. Our results underscore the importance of aerosol absorption in interpreting temperate sub-Neptune spectra, and motivate future JWST observations and laboratory work to break these degeneracies.

The different spatial distributions of N-bearing and O-bearing species, as is well known towards Orion~KL, is one of the long-lasting mysteries. We conducted a survey observation and chemical modeling study to investigate if the different distributions of O- and N-bearing species are widely recognized in general star-forming regions. First, we report our observational results of complex organic molecules (COMs) with the 45~m radio telescope at the Nobeyama Radio Observatory towards eight star-forming regions. Through our spectral survey ranging from 80 to 108~GHz, we detected CH$_3$OH, HCOOCH$_3$, CH$_3$OCH$_3$, (CH$_3$)$_2$CO, CH$_3$CHO, CH$_3$CH$_2$CN, CH$_2$CHCN, and NH$_2$CHO. Their molecular abundances were derived via the rotation diagram and the least squares methods. We found that N-bearing molecules, tend to show stronger correlations with other N-bearing molecules rather than O-bearing molecules. While G10.47+0.03 showed high fractional abundances of N-bearing species, those in NGC6334F were not so rich, being less than 0.01 compared to CH$_3$OH. Then, the molecular abundances towards these sources were evaluated by chemical modeling with NAUTILUS three-phase gas-grain chemical code. Through the simulations of time evolutions for the abundances of COMs, we suggest that observed correlations of fractional abundances between COMs can be explained by the combination of the different temperature structures inside the hot cores and the different evolutionary phase. Since our modeling could not fully explain the observed excitation temperatures, it is important to investigate the efficiency of grain surface reactions and their activation barriers, and the binding energy of COMs to further promote our understanding.

We produce a clean and well-characterised catalogue of objects within 100\,pc of the Sun from the \G\ Early Data Release 3. We characterise the catalogue through comparisons to the full data release, external catalogues, and simulations. We carry out a first analysis of the science that is possible with this sample to demonstrate its potential and best practices for its use. The selection of objects within 100\,pc from the full catalogue used selected training sets, machine-learning procedures, astrometric quantities, and solution quality indicators to determine a probability that the astrometric solution is reliable. The training set construction exploited the astrometric data, quality flags, and external photometry. For all candidates we calculated distance posterior probability densities using Bayesian procedures and mock catalogues to define priors. Any object with reliable astrometry and a non-zero probability of being within 100\,pc is included in the catalogue. We have produced a catalogue of \NFINAL\ objects that we estimate contains at least 92\% of stars of stellar type M9 within 100\,pc of the Sun. We estimate that 9\% of the stars in this catalogue probably lie outside 100\,pc, but when the distance probability function is used, a correct treatment of this contamination is possible. We produced luminosity functions with a high signal-to-noise ratio for the main-sequence stars, giants, and white dwarfs. We examined in detail the Hyades cluster, the white dwarf population, and wide-binary systems and produced candidate lists for all three samples. We detected local manifestations of several streams, superclusters, and halo objects, in which we identified 12 members of \G\ Enceladus. We present the first direct parallaxes of five objects in multiple systems within 10\,pc of the Sun.

The newly accessible mid-infrared (MIR) window offered by the James Webb Space Telescope (JWST) for exoplanet imaging is expected to provide valuable information to characterize their atmospheres. In particular, coronagraphs on board the JWST Mid-InfraRed instrument (MIRI) are capable of imaging the coldest directly imaged giant planets at the wavelengths where they emit most of their flux. The MIRI coronagraphs have been specially designed to detect the NH3 absorption around 10.5 microns, which has been predicted by atmospheric models. We aim to assess the presence of NH3 while refining the atmospheric parameters of one of the coldest companions detected by directly imaging GJ 504 b. Its mass is still a matter of debate and depending on the host star age estimate, the companion could either be placed in the brown dwarf regime or in the young Jovian planet regime. We present an analysis of MIRI coronagraphic observations of the GJ 504 system. We took advantage of previous observations of reference stars to build a library of images and to perform a more efficient subtraction of the stellar diffraction pattern. We detected the presence of NH3 at 12.5 sigma in the atmosphere, in line with atmospheric model expectations for a planetary-mass object and observed in brown dwarfs within a similar temperature range. The best-fit model with Exo-REM provides updated values of its atmospheric parameters, yielding a temperature of Teff = 512 K and radius of R = 1.08 RJup. These observations demonstrate the capability of MIRI coronagraphs to detect NH3 and to provide the first MIR observations of one of the coldest directly imaged companions. Overall, NH3 is a key molecule for characterizing the atmospheres of cold planets, offering valuable insights into their surface gravity. These observations provide valuable information for spectroscopic observations planned with JWST.

We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

Lensed quasars are key to many areas of study in astronomy, offering a unique probe into the intermediate and far universe. However, finding lensed quasars has proved difficult despite significant efforts from large collaborations. These challenges have limited catalogues of confirmed lensed quasars to the hundreds, despite theoretical predictions that they should be many times more numerous. We train machine learning classifiers to discover lensed quasar candidates. By using semi-supervised learning techniques we leverage the large number of potential candidates as unlabelled training data alongside the small number of known objects, greatly improving model performance. We present our two most successful models: (1) a variational autoencoder trained on millions of quasars to reduce the dimensionality of images for input to a gradient boosting classifier that can make accurate predictions and (2) a convolutional neural network trained on a mix of labelled and unlabelled data via virtual adversarial training. These models are both capable of producing high-quality candidates, as evidenced by our discovery of GRALJ140833.73+042229.98. The success of our classifier, which uses only images, is particularly exciting as it can be combined with existing classifiers, which use other data than images, to improve the classifications of both models and discover more lensed quasars.

The ARIEL (Atmospheric Remote-sensing Exoplanet Large-survey) mission concept is one of the three M4 mission candidates selected by the European Space Agency (ESA) for a Phase A study, competing for a launch in 2026. ARIEL has been designed to study the physical and chemical properties of a large and diverse sample of exoplanets, and through those understand how planets form and evolve in our galaxy. Here we describe the assumptions made to estimate an optimal sample of exoplanets - including both the already known exoplanets and the "expected" ones yet to be discovered - observable by ARIEL and define a realistic mission scenario. To achieve the mission objectives, the sample should include gaseous and rocky planets with a range of temperatures around stars of different spectral type and metallicity. The current ARIEL design enables the observation of ~1000 planets, covering a broad range of planetary and stellar parameters, during its four year mission lifetime. This nominal list of planets is expected to evolve over the years depending on the new exoplanet discoveries.

The long-term dynamical future of the Sun's planets has been simulated and statistically analyzed in great detail, but most prior work considers the solar system as completely isolated, neglecting the potential influence of field star passages. To understand the dynamical significance of field star encounters, we simulate several thousand realizations of the modern solar system in the presence of passing field stars for 5 Gyrs. We find that the impulse gradient of the strongest stellar encounter largely determines the net dynamical effect of field stars. Because the expected strength of such an encounter is uncertain by multiple orders of magnitude, the possible significance of field stars can be large. Our simulations indicate that isolated models of the solar system can underestimate the degree of our giant planets' future secular orbital changes by over an order of magnitude. In addition, our planets and Pluto are significantly less stable than previously thought. Field stars transform Pluto from a completely stable object over 5 Gyrs to one with a ~5% instability probability. Furthermore, field stars increase the odds of Mercury's instability by ~50-80%. We also find a ~0.3% chance that Mars will be lost through collision or ejection and a ~0.2% probability that Earth will be involved in a planetary collision or ejected. Compared to previously studied instabilities in isolated solar systems models, those induced by field stars are much more likely to involve the loss of multiple planets. In addition, they typically happen sooner in our solar system's future, making field star passages the most likely cause of instability for the next 4-4.5 Gyrs.

An instability among the giant planets' orbits can match many aspects of the Solar System's current orbital architecture. We explore the possibility that this dynamical instability was triggered by the close passage of a star or substellar object during the Sun's embedded cluster phase. We run N-body simulations starting with the giant planets in a resonant chain and an outer planetesimal disk, with a wide-enough planet-disk separation to preserve the planets' orbital stability for $>$100 Myr. We subject the system to a single flyby, testing a wide range in flyby mass, velocity and closest approach distance. We find a variety of outcomes, from flybys that over-excite the system (or strip the planets entirely) to flybys too weak to perturb the planets at all. An intermediate range of flybys triggers a dynamical instability that matches the present-day Solar System. Successful simulations -- that match the giant planets' orbits without over-exciting the cold classical Kuiper belt -- are characterized by the flyby of a substellar object ($3-30 M_{Jup}$) passing within 20 au of the Sun. We performed Monte Carlo simulations of the Sun's birth cluster phase, parameterized by the product of the stellar density $\eta$ and the cluster lifetime $T$. The balance between under- and over-excitation of the young Solar System is at $\eta T \approx 5 \times 10^4$~Myr pc$^{-3}$, in a range consistent with previous work. We find a probability of $\sim$1% that the Solar System's dynamical instability was triggered by a substellar flyby. The probability increases to $\sim$5% if the occurrence rate of free-floating planets and low-mass brown dwarfs is modestly higher than predicted by standard stellar initial mass functions.

Ground-based gamma-ray astronomy is now well established as a key observational approach to address critical topics at the frontiers of astroparticle physics and high-energy astrophysics. Whilst the field of TeV astronomy was once dominated by arrays of atmospheric Cherenkov Telescopes, ground-level particle detection has now been demonstrated to be an equally viable and strongly complementary approach. Ground-level particle detection provides continuous monitoring of the overhead sky, critical for the mapping of extended structures and capturing transient phenomena. As demonstrated by HAWC and LHAASO, the technique provides the best available sensitivity above a few tens of TeV, and for the first time access to the PeV energy range. Despite the success of this approach, there is so far no major ground-level particle-based observatory with access to the Southern sky. HESS, located in Namibia, is the only major gamma-ray instrument in the Southern Hemisphere, and has shown the extraordinary richness of the inner galaxy in the TeV band, but is limited in terms of field of view and energy reach. SWGO is an international effort to construct the first wide-field instrument in the south with deep sensitivity from 100s of GeV into the PeV domain. The project is now close to the end of its development phase and planning for construction of the array in Chile has begun. Here we describe the baseline design, expected sensitivity and resolution, and describe in detail the main scientific topics that will be addressed by this new facility and its initial phase SWGO-A. We show that SWGO will have a transformational impact on a wide range of topics from cosmic-ray acceleration and transport to the nature of dark matter. SWGO represents a key piece of infrastructure for multi-messenger astronomy in the next decade, with strong scientific synergies with the nearby CTA Observatory.

While small, Neptune-like planets are among the most abundant exoplanets, our understanding of their atmospheric structure and dynamics remains sparse. In particular, many unknowns remain on the way moist convection works in these atmospheres where condensable species are heavier than the non-condensable background gas. While it has been predicted that moist convection could shut-down above some threshold abundance of these condensable species, this prediction is based on simple linear analysis and relies on strong assumptions on the saturation of the atmosphere. To investigate this issue, we develop a 3D cloud resolving model for H2 atmospheres with large amounts of condensable species and apply this model to a prototypical temperate Neptune-like planet -- K2-18b. Our model confirms the shut-down of moist convection and the onset of a stably stratified layer in the atmosphere, leading to much hotter deep atmospheres and interiors. Our 3D simulations further provide quantitative estimates of the turbulent mixing in this stable layer, which is a key driver of the cycling of condensables in the atmosphere. This allows us to build a very simple, yet realistic 1D model that captures the most salient features of the structure of Neptune-like atmospheres. Our qualitative findings on the behavior of moist convection in hydrogen atmospheres go beyond temperate planets and should also apply to the regions where iron and silicates condense in the deep interior of H2-dominated planets. Finally, we use our model to investigate the likelihood of a liquid ocean beneath a H2 dominated atmosphere on K2-18b. We find that the planet would need to have a very high albedo (>0.5-0.6) to sustain a liquid ocean. However, due to the spectral type of the star, the amount of aerosol scattering that would be needed to provide such a high albedo is inconsistent with the latest observational data.

Exoplanet searches through space-based photometric time series have shown to be very efficient in recent years. However, follow-up efforts on the detected planet candidates have been demonstrated to be critical to uncover the true nature of the transiting objects. In this paper we show a detailed analysis of one of those false positives hidden as planetary signals. In this case, the candidate KOI-3886.01 showed clear evidence of a planetary nature from various techniques. Indeed, the properties of the fake planet set it among the most interesting and promising for the study of planetary evolution as the star leaves the main sequence. To unveil the true nature of this system, we present a complete set of observational techniques including high-spatial resolution imaging, high-precision photometric time series (showing eclipses, phase curve variations and asteroseismology signals), high-resolution spectroscopy and derived radial velocities, to unveil the true nature of this planet candidate. We find that KOI-3886.01 is an interesting false positive case: a hierarchical triple system composed by a $\sim$K2III giant star (KOI-3886A) accompanied by a close-in eclipsing binary formed by a subgiant $\sim$G4IV star (KOI-3886B) and a brown dwarf (KOI-3886C). In particular, KOI-3886C is one of the most irradiated brown dwarfs known to date, showing the largest radius in this substellar regime. It is also the first eclipsing brown dwarf known around an evolved star. In this paper we highlight the relevance of complete sets of follow-up observations to extrasolar planets detected by the transit technique using large-pixel photometers such as Kepler and TESS, and in the future, PLATO. In particular, multi-color high-spatial resolution imaging was the first hint toward ruling out the planet scenario in this system.

There has been an incredibly large investment in obtaining high-resolution stellar spectra for determining chemical abundances of stars. This information is crucial to answer fundamental questions in Astronomy by constraining the formation and evolution scenarios of the Milky Way as well as the stars and planets residing in it. We have just entered a new era, in which chemical abundances of FGK-type stars are being produced at industrial scales, where the observations, reduction, and analysis of the data are automatically performed by machines. Here we review the latest human efforts to assess the accuracy and precision of such industrial abundances by providing insights in the steps and uncertainties associated with the process of determining stellar abundances. To do so, we highlight key issues in the process of spectral analysis for abundance determination, with special effort in disentangling sources of uncertainties. We also provide a description of current and forthcoming spectroscopic surveys, focusing on their reported abundances and uncertainties. This allows us to identify which elements and spectral lines are best and why. Finally, we make a brief selection of main scientific questions the community is aiming to answer with abundances.

The improvements in the precision of the published data in \textit{Gaia} EDR3 with respect to \textit{Gaia} DR2, particularly for parallaxes and proper motions, offer the opportunity to increase the number of known open clusters in the Milky Way by detecting farther and fainter objects that have so far go unnoticed. Our aim is to keep completing the open cluster census in the Milky Way with the detection of new stellar groups in the Galactic disc. We use \textit{Gaia} EDR3 up to magnitude $G = 18$ mag, increasing in one unit the magnitude limit and therefore the search volume explored in our previous studies. We use the \texttt{OCfinder} method to search for new open clusters in \textit{Gaia} EDR3 using a Big Data environment. As a first step, \texttt{OCfinder} identifies stellar statistical overdensities in the five dimensional astrometric space (position, parallax and proper motions) using the \texttt{DBSCAN} clustering algorithm. Then, these overdensities are classified into random statistical overdensities or real physical open clusters using a deep artificial neural network trained on well-characterised $G$, $G_{\rm BP} - G_{\rm RP}$ colour-magnitude diagrams. We report the discovery of $664$ new open clusters within the Galactic disc, most of them located beyond $1$ kpc from the Sun. From the estimation of ages, distances and line-of-sight extinctions of these open clusters, we see that young clusters align following the Galactic spiral arms while older ones are dispersed in the Galactic disc. Furthermore, we find that most open clusters are located at low Galactic altitudes with the exception of a few groups older than $1$ Gyr. We show the success of the \texttt{OCfinder} method leading to the discovery of a total of $1\,310$ open clusters (joining the discoveries here with the previous ones based on \textit{Gaia} DR2)[abridged]

Context: Chamaeleon is the southernmost low-mass star-forming complex within 200 pc from the Sun. Its stellar population has been extensively studied in the past, but the current census of the stellar content is not complete yet and deserves further investigation. Aims: We take advantage of the second data release of the \textit{Gaia} space mission to expand the census of stars in Chamaeleon and to revisit the properties of the stellar populations associated to the Chamaeleon I (Cha I) and Chamaeleon II (Cha II) dark clouds. Methods: We perform a membership analysis of the sources in the \textit{Gaia} catalogue over a field of 100 deg$^{2}$ encompassing the Chamaeleon clouds, and use this new census of cluster members to investigate the 6D structure of the complex. Results: We identify 188 and 41 high-probability members of the stellar populations in Cha I and Cha II, respectively, including 19 and 7 new members. Our sample covers the magnitude range from $G=6$ to $G=20$ mag in Cha I, and from $G=12$ to $G=18$ mag in Cha II. We confirm that the northern and southern subgroups of Cha I are located at different distances ($191.4^{+0.8}_{-0.8}$ pc and $186.7^{+1.0}_{-1.0}$ pc), but they exhibit the same space motion within the reported uncertainties. Cha II is located at a distance of $197.5^{+1.0}_{-0.9}$ pc and exhibits a space motion that is consistent with Cha I within the admittedly large uncertainties on the spatial velocities of the stars that come from radial velocity data. The median age of the stars derived from the Hertzsprung-Russell diagram (HRD) and stellar models is about 1-2 Myr, suggesting that they are somewhat younger than previously thought. We do not detect significant age differences between the Chamaeleon subgroups, but we show that Cha II exhibits a higher fraction of disc-bearing stars compared to Cha I.

We aim to develop a new method to infer the sub-beam probability density function (PDF) of H2 column densities and the dense gas mass within molecular clouds using spatially unresolved observations of molecular emission lines in the 3 mm band. We model spatially unresolved line integrated intensity measurements as the average of an emission function weighted by the sub-beam column density PDF. The emission function, which expresses the line integrated intensity as a function of the gas column density, is an empirical fit to high resolution (< 0.05 pc) multi-line observations of the Orion B molecular cloud. The column density PDF is assumed to be parametric, composed of a lognormal distribution at moderate column densities and a power law distribution at higher column densities. To estimate the sub-beam column density PDF, the emission model is combined with a Bayesian inversion algorithm (the Beetroots code), which takes account of thermal noise and calibration errors. We validate our method by demonstrating that it recovers the true column density PDF of the Orion B cloud, reproducing the observed emission line integrated intensities. We apply the method to 12CO(J=1-0), 13CO(J=1-0), C18O(J=1-0), HCN(J=1-0), HCO+(J=1-0) and N2H+(J=1-0) observations of a 700 x 700 pc2 field of view (FoV) in the nearby galaxy M51. On average, the model reproduces the observed intensities within 30%. The column density PDFs obtained for the spiral arm region within our test FoV are dominated by a power-law tail at high column densities, with slopes that are consistent with gravitational collapse. Outside the spiral arm, the column density PDFs are predominantly lognormal, consistent with supersonic isothermal turbulence. We calculate the mass associated with the powerlaw tail of the column density PDFs and observe a strong, linear correlation between this mass and the 24$\mu$m surface brightness.

We introduce a novel technique to mitigate the adverse effects of atmospheric turbulence on astronomical imaging. Utilizing a video-to-image neural network trained on simulated data, our method processes a sliding sequence of short-exposure ($\sim$0.2s) stellar field images to reconstruct an image devoid of both turbulence and noise. We demonstrate the method with simulated and observed stellar fields, and show that the brief exposure sequence allows the network to accurately associate speckles to their originating stars and effectively disentangle light from adjacent sources across a range of seeing conditions, all while preserving flux to a lower signal-to-noise ratio than an average stack. This approach results in a marked improvement in angular resolution without compromising the astrometric stability of the final image.