The bright star $\pi$ Men was chosen as the first target for a radial velocity follow-up to test the performance of ESPRESSO, the new high-resolution spectrograph at the ESO's Very-Large Telescope (VLT). The star hosts a multi-planet system (a transiting 4 M$_\oplus$ planet at $\sim$0.07 au, and a sub-stellar companion on a $\sim$2100-day eccentric orbit) which is particularly appealing for a precise multi-technique characterization. With the new ESPRESSO observations, that cover a time span of 200 days, we aim to improve the precision and accuracy of the planet parameters and search for additional low-mass companions. We also take advantage of new photometric transits of $\pi$ Men c observed by TESS over a time span that overlaps with that of the ESPRESSO follow-up campaign. We analyse the enlarged spectroscopic and photometric datasets and compare the results to those in the literature. We further characterize the system by means of absolute astrometry with Hipparcos and Gaia. We used the spectra of ESPRESSO for an independent determination of the stellar fundamental parameters. We present a precise characterization of the planetary system around $\pi$ Men. The ESPRESSO radial velocities alone (with typical uncertainty of 10 cm/s) allow for a precise retrieval of the Doppler signal induced by $\pi$ Men c. The residuals show an RMS of 1.2 m/s, and we can exclude companions with a minimum mass less than $\sim$2 M$_\oplus$ within the orbit of $\pi$ Men c). We improve the ephemeris of $\pi$ Men c using 18 additional TESS transits, and in combination with the astrometric measurements, we determine the inclination of the orbital plane of $\pi$ Men b with high precision ($i_{b}=45.8^{+1.4}_{-1.1}$ deg). This leads to the precise measurement of its absolute mass $m_{b}=14.1^{+0.5}_{-0.4}$ M$_{Jup}$, and shows that the planetary orbital planes are highly misaligned.

We characterized the transiting planetary system orbiting the G2V star K2-38 using the new-generation echelle spectrograph ESPRESSO. We carried out a photometric analysis of the available K2 photometric light curve of this star to measure the radius of its two known planets. Using 43 ESPRESSO high-precision radial velocity measurements taken over the course of 8 months along with the 14 previously published HIRES RV measurements, we modeled the orbits of the two planets through a MCMC analysis, significantly improving their mass measurements. Using ESPRESSO spectra, we derived the stellar parameters, $T_{\rm eff}$=5731$\pm$66, $\log g$=4.38$\pm$0.11~dex, and $[Fe/H]$=0.26$\pm$0.05~dex, and thus the mass and radius of K2-38, $M_{\star}$=1.03 $^{+0.04}_{-0.02}$~M$_{\oplus}$ and $R_{\star}$=1.06 $^{+0.09}_{-0.06}$~R$_{\oplus}$. We determine new values for the planetary properties of both planets. We characterize K2-38b as a super-Earth with $R_{\rm P}$=1.54$\pm$0.14~R$_{\rm \oplus}$ and $M_{\rm p}$=7.3$^{+1.1}_{-1.0}$~M$_{\oplus}$, and K2-38c as a sub-Neptune with $R_{\rm P}$=2.29$\pm$0.26~R$_{\rm \oplus}$ and $M_{\rm p}$=8.3$^{+1.3}_{-1.3}$~M$_{\oplus}$. We derived a mean density of $\rho_{\rm p}$=11.0$^{+4.1}_{-2.8}$~g cm$^{-3}$ for K2-38b and $\rho_{\rm p}$=3.8$^{+1.8}_{-1.1}$~g~cm$^{-3}$ for K2-38c, confirming K2-38b as one of the densest planets known to date. The best description for the composition of K2-38b comes from an iron-rich Mercury-like model, while K2-38c is better described by a rocky model with a H2 envelope. The maximum collision stripping boundary shows how giant impacts could be the cause for the high density of K2-38b. The irradiation received by each planet places them on opposite sides of the radius valley. We find evidence of a long-period signal in the radial velocity time-series whose origin could be linked to a 0.25-3~M$_{\rm J}$ planet or stellar activity.

We aim to confirm the presence of Proxima b using independent measurements obtained with the new ESPRESSO spectrograph, and refine the planetary parameters taking advantage of its improved precision. We analysed 63 spectroscopic ESPRESSO observations of Proxima taken during 2019. We obtained radial velocity measurements with a typical radial velocity photon noise of 26 cm/s. We ran a joint MCMC analysis on the time series of the radial velocity and full-width half maximum of the cross-correlation function to model the planetary and stellar signals present in the data, applying Gaussian process regression to deal with stellar activity. We confirm the presence of Proxima b independently in the ESPRESSO data. The ESPRESSO data on its own shows Proxima b at a period of 11.218 $\pm$ 0.029 days, with a minimum mass of 1.29 $\pm$ 0.13 Me. In the combined dataset we measure a period of 11.18427 $\pm$ 0.00070 days with a minimum mass of 1.173 $\pm$ 0.086 Me. We find no evidence of stellar activity as a potential cause for the 11.2 days signal. We find some evidence for the presence of a second short-period signal, at 5.15 days with a semi-amplitude of merely 40 cm/s. If caused by a planetary companion, it would correspond to a minimum mass of 0.29 $\pm$ 0.08 Me. We find that the FWHM of the CCF can be used as a proxy for the brightness changes and that its gradient with time can be used to successfully detrend the radial velocity data from part of the influence of stellar activity. The activity-induced radial velocity signal in the ESPRESSO data shows a trend in amplitude towards redder wavelengths. Velocities measured using the red end of the spectrograph are less affected by activity, suggesting that the stellar activity is spot-dominated. The data collected excludes the presence of extra companions with masses above 0.6 Me at periods shorter than 50 days.

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

Ultra-hot giant exoplanets receive thousands of times Earth's insolation. Their high-temperature atmospheres (>2,000 K) are ideal laboratories for studying extreme planetary climates and chemistry. Daysides are predicted to be cloud-free, dominated by atomic species and substantially hotter than nightsides. Atoms are expected to recombine into molecules over the nightside, resulting in different day-night chemistry. While metallic elements and a large temperature contrast have been observed, no chemical gradient has been measured across the surface of such an exoplanet. Different atmospheric chemistry between the day-to-night ("evening") and night-to-day ("morning") terminators could, however, be revealed as an asymmetric absorption signature during transit. Here, we report the detection of an asymmetric atmospheric signature in the ultra-hot exoplanet WASP-76b. We spectrally and temporally resolve this signature thanks to the combination of high-dispersion spectroscopy with a large photon-collecting area. The absorption signal, attributed to neutral iron, is blueshifted by -11+/-0.7 km s-1 on the trailing limb, which can be explained by a combination of planetary rotation and wind blowing from the hot dayside. In contrast, no signal arises from the nightside close to the morning terminator, showing that atomic iron is not absorbing starlight there. Iron must thus condense during its journey across the nightside.

Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficient, robust, and automatic method to fit ammonia lines with multiple velocity components, generalizable to other molecular species. We also present CRISPy, a Python package to identify filament spines in 3D images (e.g., position-position-velocity cubes), along with a complementary technique to sort fitted velocity components into velocity-coherent filaments. In NGC 1333, we find a wealth of velocity gradient structures on a beam-resolved scale of ~0.05 pc. Interestingly, these local velocity gradients are not randomly oriented with respect to filament spines and their perpendicular, i.e., radial, component decreases in magnitude towards the spine for many filaments. Together with remarkably constant velocity gradients on larger scales along many filaments, these results suggest a scenario in which gas falling onto filaments is progressively damped and redirected to flow along these filaments.

Phosphorus-bearing species are an essential key to form life on Earth, however they have barely been detected in the interstellar medium. Since only PN and PO have been identified so far towards star-forming regions, the chemical formation pathways of P-bearing molecules are not easy to constrain and are thus highly debatable. An important factor still missing in the chemical models is the initial elemental abundance of phosphorus, i.e. the depletion level of P. In order to overcome this problem, we study P-bearing species in diffuse/translucent clouds. In these objects phosphorus is mainly in the gas phase and therefore the elemental initial abundance needed in our chemical simulations corresponds to the cosmic one and is thus well constrained. An advanced chemical model with an updated P-chemistry network has been used. Single-pointing observations were performed with the IRAM 30m telescope towards the line of sight to the blazar B0355+508 aiming for the (2-1) transitions of PN, PO, HCP and CP. This line of sight incorporates five diffuse/translucent clouds. The (2-1) transitions of the PN, PO, HCP and CP were not detected. We report detections of the (1-0) lines of $\mathrm{^{13}CO}$, HNC and CN along with a first detection of $\mathrm{C^{34}S}$. We have reproduced the observations of HNC, CN, CS and CO in every cloud by applying typical conditions for diffuse/translucent clouds. According to our best-fit model, the most abundant P-bearing species are HCP and CP ($\sim10^{-10}$), followed by PN, PO and $\mathrm{PH_3}$ ($\sim10^{-11}$). We show that the production of P-bearing species is favoured towards translucent rather than diffuse clouds, where the environment provides a stronger shielding from the interstellar radiation. Based on our improved model, the (1-0) transitions of HCP, CP, PN and PO are expected to be detectable with estimated intensities up to $\sim200$ mK.

Context. Magnetic fields can affect significantly the star formation process. The theory of the magnetically-driven collapse in a uniform field predicts that initially the contraction happens along the field lines. When the gravitational pull grows strong enough, the magnetic field lines pinch inwards, giving rise to a characteristic hourglass shape. Aims. We investigate the magnetic field structure of a young Class 0 object, IRAS 15398-3359, embedded in the Lupus I cloud. Previous observations at large scales suggest that this source evolved in an highly magnetised environment. This object thus appears an ideal candidate to study the magnetically driven core collapse in the low-mass regime. Methods. We have performed polarisation observations of IRAS 15398-3359 at 214$\mu$m using the SOFIA/HAWC+ instrument, thus tracing the linearly polarised thermal emission of cold dust. Results. Our data unveil a significant bend of the magnetic field lines due to the gravitational pull. The magnetic field appears ordered and aligned with the large-scale B-field of the cloud and with the outflow direction. We estimate a magnetic field strength of $B= 78 \mu$G, expected to be accurate within a factor of two. The measured mass-to-flux parameter is $\lambda= 0.95$, indicating that the core is in a transcritical regime.

L1521E is a dense starless core in Taurus that was found to have relatively low molecular depletion by earlier studies, thus suggesting a recent formation. We aim to characterize the chemical structure of L1521E and compare it to the more evolved L1544 pre-stellar core. We have obtained $\sim$2.5$\times$2.5 arcminute maps toward L1521E using the IRAM-30m telescope in transitions of various species. We derived abundances for the species and compared them to those obtained toward L1544. We estimated CO depletion factors. Similarly to L1544, $c$-C$_3$H$_2$ and CH$_3$OH peak at different positions. Most species peak toward the $c$-C$_3$H$_2$ peak. The CO depletion factor derived toward the $Herschel$ dust peak is 4.3$\pm$1.6, which is about a factor of three lower than that toward L1544. The abundances of sulfur-bearing molecules are higher toward L1521E than toward L1544 by factors of $\sim$2-20. The abundance of methanol is similar toward the two cores. The higher abundances of sulfur-bearing species toward L1521E than toward L1544 suggest that significant sulfur depletion takes place during the dynamical evolution of dense cores, from the starless to pre-stellar stage. The CO depletion factor measured toward L1521E suggests that CO is more depleted than previously found. Similar CH$_3$OH abundances between L1521E and L1544 hint that methanol is forming at specific physical conditions in Taurus, characterized by densities of a few $\times$10$^4$ cm$^{-3}$ and $N$(H$_2$)$\gtrsim$10$^{22}$ cm$^{-2}$, when CO starts to catastrophically freeze-out, while water can still be significantly photodissociated, so that the surfaces of dust grains become rich in solid CO and CH$_3$OH, as already found toward L1544. Methanol can thus provide selective crucial information about the transition region between dense cores and the surrounding parent cloud.

Context. The deuterium fraction in low-mass prestellar cores is a good diagnostic indicator of the initial phases of star formations, and it is also a fundamental quantity to infer the ionisation degree in these objects. Aims. With the analysis of multiple transitions of $\rm N_2H^+$, $\rm N_2D^+$, $\rm HC^{18}O^+$ and $\rm DCO^+$ we are able to determine the molecular column density maps and the deuterium fraction in $\rm N_2H^+$ and $\rm HCO^+$ toward the prototypical prestellar core L1544. This is the preliminary step to derive the ionisation degree in the source. Methods. We use a non-local thermodynamic equilibrium (non-LTE) radiative transfer code, combined with the molecular abundances derived from a chemical model, to infer the excitation conditions of all the observed transitions, which allows us to derive reliable maps of each molecule's column density. The ratio between the column density of a deuterated species and its non-deuterated counterpart gives the searched deuteration level. Results. The non-LTE analysis confirms that, for the analysed molecules, higher-J transitions are characterised by excitation temperatures $\approx 1-2\,$K lower than the lower-J ones. The chemical model that provides the best fit to the observational data predict the depletion of $\rm N_2H^+$ and to a lesser extent of $\rm N_2D^+$ in the innermost region. The peak values for the deuterium fraction that we find are $\mathrm{D/H_{N_2H^+}} = 0.26^{+0.15}_{-0.14}$ and $\mathrm{D/H_{HCO^+}} = 0.035^{+0.015}_{-0.012}$, in good agreement with previous estimates in the source.

We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process.

We carried out a parameter-space exploration of the ammonia abundance in the pre-stellar core L1544, where it has been observed to increase toward the center of the core with no signs of freeze-out onto grain surfaces. We considered static and dynamical physical models coupled with elaborate chemical and radiative transfer calculations, and explored the effects of varying model parameters on the (ortho+para) ammonia abundance profile. None of our models are able to reproduce the inward-increasing tendency in the observed profile; ammonia depletion always occurs in the center of the core. In particular, our study shows that including the chemical desorption process, where exothermic association reactions on the grain surface can result in the immediate desorption of the product molecule, leads to ammonia abundances that are over an order of magnitude above the observed level in the innermost 15000 au of the core - at least when one employs a constant efficiency for the chemical desorption process irrespective of the ice composition. Our results seemingly constrain the chemical desorption efficiency of ammonia on water ice to below 1%. It is increasingly evident that time-dependent effects must be considered so that the results of chemical models can be reconciled with observations.

Adaptive optics (AO) instruments for the future extremely large telescopes (ELTs) are characterized by advanced optical systems with diffraction-limited optical quality. Low geometric distortion is also crucial for high accuracy astrometric applications. Optical alignment of such systems is a crucial step of the instrument integration. Due to relative inaccessibility of these giant instruments, automatic alignment methods are also favored to improve the instrument availability after major events, such as extraordinary maintenance. The proposed alignment concept for these systems is described: the notable example which is analyzed here is the case of the multiconjugate AO relay for the future ELT. The results of ray-tracing simulations carried out to validate the method are discussed in detail, covering the error sources, which could degrade the alignment performance.

We perform a virial analysis of starless dense cores in three nearby star-forming regions : L1688 in Ophiuchus, NGC 1333 in Perseus, and B18 in Taurus. Our analysis takes advantage of comprehensive kinematic information for the dense gas in all of these regions made publicly available through the Green Bank Ammonia Survey Data Release 1, which used to estimate internal support against collapse. We combine this information with ancillary data used to estimate other important properties of the cores, including continuum data from the James Clerk Maxwell Telescope Gould Belt Survey for core identification, core masses, and core sizes. Additionally, we used \textit{Planck} and \textit{Herschel}-based column density maps for external cloud weight pressure, and Five College Radio Astronomy Observatory $^{13}$CO observations for external turbulent pressure. Our self-consistent analysis suggests that many dense cores in all three star-forming regions are not bound by gravity alone, but rather require additional pressure confinement to remain bound. Unlike a recent, similar study in Orion~A, we find that turbulent pressure represents a significant portion of the external pressure budget. Our broad conclusion emphasizing the importance of pressure confinement in dense core evolution, however, agrees with earlier work.

Context: The $^{15}$N fractionation has been observed to show large variations among astrophysical sources, depending both on the type of target and on the molecular tracer used. These variations cannot be reproduced by the current chemical models. Aims: Until now, the $^{14}$N/$^{15}$N ratio in N$_2$H$^+$ has been accurately measured in only one prestellar source, L1544, where strong levels of fractionation, with depletion in $^{15}$N, are found ($^{14}$N/$^{15}$N$\: \approx 1000$). In this paper we extend the sample to three more bona fide prestellar cores, in order to understand if the antifractionation in N$_2$H$^+$ is a common feature of this kind of sources. Methods: We observed N$_2$H$^+$,N$^{15}$NH$^+$ and $^{15}$NNH$^+$ in L183, L429 and L694-2 with the IRAM 30m telescope. We modeled the emission with a non-local radiative transfer code in order to obtain accurate estimates of the molecular column densities, including the one for the optically thick N$_2$H$^+$. We used the most recent collisional rate coefficients available, and with these we also re-analysed the L1544 spectra previously published. Results: The obtained isotopic ratios are in the range $630-770$ and significantly differ with the value, predicted by the most recent chemical models, of $\approx 440$, close to the protosolar value. Our prestellar core sample shows high level of depletion of $^{15}$N in diazenylium, as previously found in L1544. A revision of the N chemical networks is needed in order to explain these results.

Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency.

We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores.

We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined.

We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, ${\sigma}_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of ${\sigma}_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available.