Search

The young stellar object [BHB2007]-1 has been extensively studied in the past at radio, millimeter, and infrared wavelengths. It shows a gap in the disk and previous observations claimed the possible emission from a forming sub-stellar object, in correspondence to the disk gap. Here, we analyze a set of 8 Karl Jansky Very Large Array (VLA) observations at 15 GHz and spread over a month. We infer a slowly variable emission from the star, with a $\sim 15 \text{-} 20\%$ circular polarization detected in two of the eight observations. The latter can be related to the magnetic fields in the system, while the unpolarized and moderately varying component can be indicative of free-free emission associated with jet induced shocks or interaction of the stellar wind with dense surrounding material. We discard any relevant short flaring activities when sampling the radio light curves down to 10 seconds and find no clear evidence of emission from the sub-stellar object inferred from past observations, although deeper observations could shed further light on this., Comment: 9 pages, 6 figures, 2 tables, Accepted for publication in A&A

This paper presents the development of a deep-atmosphere, nonhydrostatic dynamical core (DyCore) targeted towards ground-thermosphere atmospheric prediction. This DyCore is based on a novel formulation of the specific internal energy equation (SIEE), which, unlike standard potential temperature formulations, is valid for variable composition atmospheres. Two versions of a SIEE are derived from basic principles. The first version, which uses a product-rule (PR) continuity equation, contains an additional compressible term and does not conserve mass. The second version, which does not use the product-rule (No-PR) in the continuity equation, contains two compressible terms and conserves mass to machine precision regardless of time truncation error. The pressure gradient and gravitational forces in the momentum balance equation are formulated in a manner appropriate for HA applications and the spectral element method (SEM) is used with Implicit-Explicit (IMEX) and Horizontally Explicit Vertically Implicit (HEVI) time-integration. These new equation sets were implemented in two U. S. Navy atmospheric models: the Nonhydrostatic Unified Model of the Atmosphere (NUMA) and the Navy Environmental Prediction sysTem Using a Nonhydrostatic Core (NEPTUNE). Numerical results using a nonhydrostatic and hydrostatic baroclinic instability, a balanced zonal flow, and HA mountain wave experiments are shown. These results are compared to existing deep-atmosphere dynamical cores, indicating that the proposed discretized IEE equation sets are viable next-generation ground-to-thermosphere DyCores.

Measurements have historically presented a problem for the consistent description of quantum theories, be it in non-relativistic quantum mechanics or in quantum field theory. Drawing on a recent surge of interest in the description of measurements in Algebraic Quantum Field theory, it was decided that this dissertation would be focused on trying to close the gap between the description of measurements proposed by K. Hepp in the 70's, considering decoherence of states in quasilocal algebras and the new framework of generally covariant measurement schemes proposed recently by C. Fewster and R. Verch. Another recent result that we shall also consider is the Frauchinger-Renner Gedankenexperiment, that has taken inspiration on Hepp's article about decoherence based measurements to arrive at a no-go result about the consistency of quantum descriptions of systems containing rational agents, we shall seek to provide a closure for the interpretation of this result. In doing so we naturally arrive at the study of the contextual properties of measurement setups., Comment: Dissertation submitted in partial fulfillment of the requirements for the degree of Master of Science in Physics at the university of S\~ao Paulo (Brazil). xvi+225 pages, 16 figures, 7 tables

Magnetic fields may play a crucial role in setting the initial conditions of massive star and star cluster formation. To investigate this, we report SOFIA-HAWC+ $214\:\mu$m observations of polarized thermal dust emission and high-resolution GBT-Argus C$^{18}$O(1-0) observations toward the massive Infrared Dark Cloud (IRDC) G28.37+0.07. Considering the local dispersion of $B$-field orientations, we produce a map of $B$-field strength of the IRDC, which exhibits values between $\sim0.03 - 1\:$mG based on a refined Davis-Chandrasekhar-Fermi (r-DCF) method proposed by Skalidis \& Tassis. Comparing to a map of inferred density, the IRDC exhibits a $B-n$ relation with a power law index of $0.51\pm0.02$, which is consistent with a scenario of magnetically-regulated anisotropic collapse. Consideration of the mass-to-flux ratio map indicates that magnetic fields are dynamically important in most regions of the IRDC. A virial analysis of a sample of massive, dense cores in the IRDC, including evaluation of magnetic and kinetic internal and surface terms, indicates consistency with virial equilibrium, sub-Alfv\'enic conditions and a dominant role for $B-$fields in regulating collapse. A clear alignment of magnetic field morphology with direction of steepest column density gradient is also detected. However, there is no preferred orientation of protostellar outflow directions with the $B-$field. Overall, these results indicate that magnetic fields play a crucial role in regulating massive star and star cluster formation and so need to be accounted for in theoretical models of these processes., Comment: Submitted to ApJ, comments welcome

We conduct a thorough study of different forms of horizontally explicit and vertically implicit (HEVI) time-integration strategies for the compressible Euler equations on spherical domains typical of nonhydrostatic global atmospheric applications. We compare the computational time and complexity of two nonlinear variants (NHEVI-GMRES and NHEVI-LU) and a linear variant (LHEVI). We report on the performance of these three variants for a number of additive Runge-Kutta Methods ranging in order of accuracy from second through fifth, and confirm the expected order of accuracy of the HEVI methods for each time-integrator. To gauge the maximum usable time-step of each HEVI method, we run simulations of a nonhydrostatic baroclinic instability for 100 days and then use this time-step to compare the time-to-solution of each method. The results show that NHEVI-LU is 2x faster than NHEVI-GMRES, and LHEVI is 5x faster than NHEVI-LU, for the idealized cases tested. The baroclinic instability and inertia-gravity wave simulations indicate that the optimal choice of time-integrator is LHEVI with either second or third order schemes, as both schemes yield similar time to solution and relative L2 error at their maximum usable time-steps. In the future, we will report on whether these results hold for more complex problems using, e.g., real atmospheric data and/or a higher model top typical of space weather applications., Comment: 36 pages, 10 figures, 2 tables

Context. Deuterium in H-bearing species is enhanced during the early stages of star formation, however, only a small number of high spatial resolution deuteration studies exist towards protostellar objects, leaving the small-scale structures unrevealed and understudied. Aims. We aim to constrain the deuterium fractionation ratios in a Class 0/I protostellar object in formaldehyde (H2CO), which has abundant deuterated isotopologues in this environment. Methods. We observed the Class 0/I protobinary system [BHB2007] 11, whose emission components are embedded in circumstellar disks that have radii of 2-3 au, using ALMA within the context of the Large Program FAUST. The system is surrounded by a complex filamentary structure connecting to the larger circumbinary disk. In this work we present the first study of formaldehyde D-fractionation towards this source with detections of H2CO 3(0,3)-2(0,2), combined with HDCO 4(2,2)-3(2,1), HDCO 4(1,4)-3(1,3) and D2CO 4(0,4)-3(0,3). These observations enable multiple velocity components associated with the methanol hotspots also uncovered by FAUST data, as well as the external envelope, to be resolved. In addition, based on the kinematics seen in the observations of the H2CO emission, we propose the presence of a second large scale outflow. Results. HDCO and D2CO are only found in the central regions of the core while H2CO is found more ubiquitously. From radiative transfer modelling, the column densities ranges found for H2CO, HDCO and D2CO are (3-8)x10$^{14}$ cm$^{-2}$, (0.8-2.9)x10$^{13}$ cm$^{-2}$ and (2.6-4.3)x10$^{12}$ cm$^{-2}$, respectively, yielding an average D/H ratio of 0.01-0.04. Following the results of kinematic modelling, the second large scale feature is inconsistent with a streamer-like nature and we thus tentatively conclude that the feature is an asymmetric molecular outflow launched by a wide-angle disk wind., Comment: 17 pages, 15 figures

This paper presents the results of full polarization observations of the massive star-forming region W75N, conducted with 3 arcsec spatial resolutions at 345 GHz using the Submillimeter Array (SMA). The magnetic field structures in the dense cores of the region are derived using the linearly polarized continuum emission. The overall magnetic field strength and orientation are found to agree with those from the previous observations. The plane-of-sky (POS) component of the magnetic field in the region was calculated to be 0.8 \pm 0.1 mG using the angular dispersion function (ADF) method. Further analyses involving the polarization-intensity gradient-local gravity method and H13CO+ (4-3) line data indicated that the cloud is undergoing global gravitational collapse and the magnetic field is shaped by gravity and outflows in the dense core regions., Comment: 16 pages, 11 figures, accepted for publication in ApJ

We have observed the low-mass protostellar source, IRAS 15398$-$3359, at a resolution of 0.$''$2-0.$''$3, as part of the Atacama Large Millimeter/Submillimeter Array Large Program FAUST, to examine the presence of a hot corino in the vicinity of the protostar. We detect nine CH$_3$OH lines including the high excitation lines with upper state energies up to 500 K. The CH$_3$OH rotational temperature and the column density are derived to be 119$^{+20}_{-26}$ K and 3.2$^{+2.5}_{-1.0}\times$10$^{18}$ cm$^{-2}$, respectively. The beam filling factor is derived to be 0.018$^{+0.005}_{-0.003}$, indicating that the emitting region of CH$_3$OH is much smaller than the synthesized beam size and is not resolved. The emitting region of three high excitation lines, 18$_{3,15}-18_{2,16}$, A ($E_u=$447 K), 19$_{3,16}-19_{2,17}$, A ($E_u=$491 K), and 20$_{3,17}-20_{2,18}$, A ($E_u=$537 K), is located within the 50 au area around the protostar, and seems to have a slight extension toward the northwest. Toward the continuum peak, we also detect one emission line from CH$_2$DOH and two features of multiple CH$_3$OCHO lines. These results, in combination with previous reports, indicate that IRAS 15398$-$3359 is a source with hybrid properties showing both hot corino chemistry rich in complex organic molecules on small scales $\sim$10 au) and warm carbon-chain chemistry (WCCC) rich in carbon-chain species on large scales ($\sim$100-1000 au). A possible implication of the small emitting region is further discussed in relation to the origin of the hot corino activity., Comment: 20pages, 3figures

Observations of protostellar envelopes are essential to understand better the process of gravitational collapse toward star and planet formation. From a theoretical perspective, magnetic fields are considered an important factor during the early stages of star formation, especially during the main accretion phase. We aim to study the relation between kinematics and magnetic fields at a very early stage of the star formation process by using data from the Atacama Pathfinder EXperiment (APEX) single dish antenna with the angular resolution of 28". We observed the two molecular lines C18O(2-1) and DCO+(3-2), toward the Class 0 young stellar object IRAS15398-3359. We implement a multi-component Gaussian fitting on the molecular data to study the kinematics. Also, we use previous polarization observations on this source to predict the influence of the magnetic field on the core. The velocity gradient along the central object can be explained as an ongoing outflow motion. We report flowing of material from the filament toward the central object, and of the merging of two velocity components in the C18O (2-1) emission around the protostar position, probably due to the merging of filamentary clouds. Our analysis shows that the large-scale magnetic field line observed previously is preferentially aligned to the rotation axis of the core.

This work presents IMplicit-EXplicit (IMEX) formulations for discontinuous Galerkin (DG) discretizations of the compressible Euler equations governing non-hydrostatic atmospheric flows. In particular, we show two different IMEX formulations that not only treat the stiffness due to the governing dynamics but also the domain discretization. We present these formulations for two different equation sets typically employed in atmospheric modeling. For both equation sets, efficient Schur complements are derived and the challenges and remedies for deriving them are discussed. The performance of these IMEX formulations of different orders are investigated on both 2D (box) and 3D (sphere) test problems and shown to achieve their theoretical rates of convergence and their efficiency with respect to both mesoscale and global applications are presented.

We present our current understanding of the formation and early evolution of protostars, protoplanetary disks, and the driving of outflows as dictated by the interplay of magnetic fields and partially ionized gas in molecular cloud cores. In recent years, the field has witnessed enormous development through sub-millimeter observations which in turn have constrained models of protostar formation. As a result of these observations % that the observations provided, the state-of-the-art theoretical understanding of the formation and evolution of young stellar objects is described. In particular, we emphasize the importance of the coupling, decoupling, and re-coupling between weakly ionized gas and the magnetic field on appropriate scales. This highlights the complex and intimate relationship between gravitational collapse and magnetic fields in young protostars., Comment: 32 pages, 15 figures, accepted as a chapter of Protostars and Planets VII

The diffusion motion of spin-bearing molecules is considerably affected in confined environments. In Nuclear Magnetic Resonance (NMR) experiments, under the presence of magnetic field gradients, this movement can be encoded in spin phase and the NMR signal attenuation due to diffusion can be evaluated. This paper considered this effect in both normal and anomalous diffusion by means of the Langevin stochastic model within the Gaussian phase approximation (GPA) for a constant gradient temporal profile. The phenomenological approach illustrates the emergence of fractional order exponential decay in the NMR signal and retrieves the classical result from Hanh echo., Comment: 8 pages, 2 figures

Numerical weather prediction is pushing the envelope of grid resolution at local and global scales alike. Aiming to model topography with higher precision, a handful of articles introduced unstructured vertical grids and tested them for dry atmospheres. The next step towards effective high-resolution unstructured grids for atmospheric modeling requires that also microphysics is independent of any vertical columns, in contrast to what is ubiquitous across operational and research models. In this paper, we present a non-column based continuous and discontinuous spectral element implementation of Kessler's microphysics with warm rain as a first step towards fully unstructured atmospheric models. We test the proposed algorithm against standard three-dimensional benchmarks for precipitating clouds and show that the results are comparable with those presented in the literature across all of the tested effective resolutions. While presented for both continuous and discontinuous spectral elements in this paper, the method that we propose can very easily be adapted to any numerical method utilized in other research and legacy codes.

This study aims to identify the properties, advantages, and drawbacks of some common (and some less common) sub-grid scale (SGS) models for large eddy simulation of low Mach compressible flows using high order spectral elements. The models investigated are the classical constant coefficient Smagorinsky-Lilly, the model by Vreman and two variants of a dynamic SGS (DSGS) model designed to stabilize finite and spectral elements for transport dominated problems. In particular, we compare one variant of DSGS that is based on a time-dependent residual version (R-DSGS) in contrast to a time-independent residual based scheme (T-DSGS). The SGS models are compared against the reference model by Smagorinsky and Lilly for their ability to: (i) stabilize the numerical solution, (ii) minimize undershoots and overshoots, (iii) capture/preserve discontinuities, and (iv) transfer energy across different length scales. These abilities are investigated on problems for: (1) passively advected tracers, (2) coupled, nonlinear system of equations exhibiting discontinuities, (3) gravity-driven flows in a stratified atmosphere, and (4) homogenous, isotropic turbulence. All models were able to preserve sharp discontinuities. Vreman and the R-DSGS models also reduce the undershoots and overshoots in the solution of linear and non-linear advection with sharp gradients. Our analysis shows that the R-DSGS and T-DSGS models are more robust than Vreman and Smagorinsky-Lilly for numerical stabilization of high-order spectral methods. The Smagorinsky and Vreman models are better able to resolve the finer flow structures in shear flows, while the nodal R-DSGS model shows better energy conservation. Overall, the nodal implementation of R-DSGS (in contrast to its element-based counterpart) is shown to outperform the other SGS models in most metrics listed above, and on par with respect to the remaining ones.

Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($<$10 K) allow the formation of thick icy dust mantles, which will be partially preserved in the future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particular CO, are heavily depleted in pre-stellar cores due to the efficient molecular freeze-out onto the surface of cold dust grains. However, N-bearing species such as NH$_3$ and, in particular, its deuterated isotopologues, appear to maintain high abundances where CO molecules are mainly in solid phase. Thanks to ALMA, we present here the first clear observational evidence of NH$_2$D freeze-out toward the L1544 pre-stellar core, suggestive of the presence of a"complete-depletion zone" within a $\simeq$1800 au radius, in agreement with astrochemical pre-stellar core model predictions. Our state-of-the-art chemical model coupled with a non-LTE radiative transfer code demonstrates that NH$_2$D becomes mainly incorporated in icy mantles in the central 2000 au and starts freezing-out already at $\simeq$7000 au. Radiative transfer effects within the pre-stellar core cause the NH$_2$D(1$_{11}$-1$_{01}$) emission to appear centrally concentrated, with a flattened distribution within the central $\simeq$3000 au, unlike the 1.3 mm dust continuum emission which shows a clear peak within the central $\simeq$1800 au. This prevented NH$_2$D freeze-out to be detected in previous observations, where the central 1000 au cannot be spatially resolved., Comment: accepted for publication in The Astrophysical Journal

We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed., Comment: 27 pages, 21 figures, 2 Tables, Accepted for publication in ApJ

Isotopic fractionation is an important tool to investigate the chemical history of our Solar System (SS). In particular, the isotopic fraction of nitrogen (14N/15N) is lower in comets and other pristine SS bodies with respect to the value measured for the protosolar nebula, suggesting a local chemical enrichment of 15N during the SS formation. Therefore, interferometric studies of nitrogen fractionation in SS precursors are imperative for us to obtain clues about our astrochemical origins. In this work, we investigated the variation of the 14N/15N ratio in one of the closest analogues of the environment in which the SS was born: the protocluster OMC-2 FIR4. We present the first comparison at high angular resolution between HCN and N2H+ using interferometric data. We analysed observations of the HCN isotopologues H13CN and HC15N in the OMC-2 FIR4 protocluster, specifically the transitions H13CN (1-0) and HC15N (1-0), from NOEMA within the context of the IRAM Seeds Of Life In Space Large Program. We combined our results with analysis of archival data obtained with ALMA of N2H+ and its 15N isotopologues. Our results show a small regional variation in the ratio for HCN from ~250 to 500. The ratios in the central regions of FIR4, where the candidate protostars are located, are largely consistent (~300). They also show little variation from the part of the protocluster known to harbour a high cosmic-ray ionisation rate, to the portion with lower rate. We also found a small variation in the ratio of N2H+ across different regions from ~200 to ~400. These results suggest that local changes in the physical parameters occurring on small linear scales probed by our observations do not seem to affect the 14N/15N ratio in either HCN or N2H+ and hence that this is independent of the molecule used. Moreover, the high level of irradiation due to cosmic rays does not affect the N-fractionation either.

Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of this conclusion, we use ammonia line data from the Green Bank Ammonia Survey and Planck-calibrated dust emission data from Herschel to estimate the masses and kinetic and gravitational energies for dense clumps in the Gould Belt clouds. We show that several types of systematic error can enhance the appearance of low kinetic-to-gravitational energy ratios: insufficient removal of foreground and background material; ignoring the kinetic energy associated with velocity differences across a resolved cloud; and over-correcting for stratification when evaluating the gravitational energy. Using an analysis designed to avoid these errors, we find that the most massive Gould Belt clumps harbor virial motions, rather than sub-virial ones. As a byproduct, we present a catalog of masses, energies, and virial energy ratios for 85 Gould Belt clumps., Comment: Submitted to ApJ

The presence of planets or sub-stellar objects still embedded in their native protoplanetary disks is indirectly suggested by disk sub-structures like gaps, cavities, and spirals. However, these companions are rarely detected. We present VLT/NACO high-contrast images in $J$, $H$, $K_S$, and $L^{\prime}$ band of the young star [BHB2007]-1 probing the inclined disk in scattered light and revealing the probable presence of a companion. The point source is detected in the $L^{\prime}$ band in spatial correspondence with complementary VLA observations. This object is constrained to have a mass in the range of 37-47 M$_{Jup}$ and is located at 50 au from the central star, inside the 70 au-large disk cavity recently imaged by ALMA, that is absent from our NACO data (down to 20 au). This mass range is compatible with the upper end derived from the size of the ALMA cavity. The NIR disk brightness is highly asymmetric around the minor axis, with the southern side 5.5 times brighter than the northern side. The constant amount of asymmetry across all wavelengths suggests that it is due to a shadow cast by a misaligned inner disk. The massive companion that we detect could, in principle, explain the possible disk misalignment, as well as the different cavity sizes inferred by the NACO and ALMA observations. The confirmation and characterization of the companion is entrusted to future observations., Comment: 9 pages, 5 figures, 3 tables. Accepted for publication on ApJ

We have observed the very low-mass Class 0 protostar IRAS 15398-3359 at scales ranging from 50 au to 1800 au, as part of the ALMA Large Program FAUST. We uncover a linear feature, visible in H2CO, SO, and C18O line emission, which extends from the source along a direction almost perpendicular to the known active outflow. Molecular line emission from H2CO, SO, SiO, and CH3OH further reveals an arc-like structure connected to the outer end of the linear feature and separated from the protostar, IRAS 15398-3359, by 1200 au. The arc-like structure is blue-shifted with respect to the systemic velocity. A velocity gradient of 1.2 km/s over 1200 au along the linear feature seen in the H2CO emission connects the protostar and the arc-like structure kinematically. SO, SiO, and CH3OH are known to trace shocks, and we interpret the arc-like structure as a relic shock region produced by an outflow previously launched by IRAS 15398-3359. The velocity gradient along the linear structure can be explained as relic outflow motion. The origins of the newly observed arc-like structure and extended linear feature are discussed in relation to turbulent motions within the protostellar core and episodic accretion events during the earliest stage of protostellar evolution.

While it is widely accepted that planets are formed in protoplanetary disks, there is still much debate on when this process happens. In a few cases protoplanets have been directly imaged, but for the vast majority of systems, disk gaps and cavities -- seen especially in dust continuum observations -- have been the strongest evidence of recent or on-going planet formation. We present ALMA observations of a nearly edge-on ($i = 75^{\circ}$) disk containing a giant gap seen in dust but not in $^{12}$CO gas. Inside the gap, the molecular gas has a warm (100 K) component coinciding in position with a tentative free-free emission excess observed with the VLA. Using 1D hydrodynamic models, we find the structure of the gap is consistent with being carved by a planet with 4-70 $M_{\rm Jup}$. The coincidence of free-free emission inside the planet-carved gap points to the planet being very young and/or still accreting. In addition, the $^{12}$CO observations reveal low-velocity large scale filaments aligned with the disk major axis and velocity coherent with the disk gas that we interpret as ongoing gas infall from the local ISM. This system appears to be an interesting case where both a star (from the environment and the disk) and a planet (from the disk) are growing in tandem., Comment: 7 pages, 4 figures, Accepted for publication in ApJ Letters

Observations indicate that molecular clouds are strongly magnetized, and that magnetic fields influence the formation of stars. A key observation supporting the conclusion that molecular clouds are significantly magnetized is that the orientation of their internal structure is closely related to that of the magnetic field. At low column densities the structure aligns parallel with the field, whereas at higher column densities, the gas structure is typically oriented perpendicular to magnetic fields, with a transition at visual extinctions $A_V\gtrsim{}3~\rm{}mag$. Here we use far-infrared polarimetric observations from the HAWC+ polarimeter on SOFIA to report the discovery of a further transition in relative orientation, i.e., a return to parallel alignment at $A_V\gtrsim{}21~\rm{}mag$ in parts of the Serpens South cloud. This transition appears to be caused by gas flow and indicates that magnetic supercriticality sets in near $A_V\gtrsim{}21~\rm{}mag$, allowing gravitational collapse and star cluster formation to occur even in the presence of relatively strong magnetic fields., Comment: 16 pages, 4 figures, Published in Nature Astronomy (August 2020). This is the authors' version before final edits. Link to the NA publication: https://www.nature.com/articles/s41550-020-1172-6

Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity towards the coherent cores in L1688 in order to study their kinematical properties in unprecedented detail. Methods : We used $\rm NH_3$ (1,1) and (2,2) data from the first data release (DR1) in the Green Bank Ammonia Survey (GAS). We first smoothed the data to a larger beam of 1' to obtain substantially more extended maps of velocity dispersion and kinetic temperature, compared to the DR1 maps. We then identified the coherent cores in the cloud and analysed the averaged line profiles towards the cores. Results : For the first time, we detected a faint (mean $\rm NH_3$(1,1) peak brightness $<$0.25 K in $T_{MB}$), supersonic component towards all the coherent cores in L1688. We fitted two components, one broad and one narrow, and derived the kinetic temperature and velocity dispersion of each component. The broad components towards all cores have supersonic linewidths ($\mathcal{M}_S \ge 1$). This component biases the estimate of the narrow dense core component's velocity dispersion by $\approx$28% and the kinetic temperature by $\approx$10%, on average, as compared to the results from single-component fits. Conclusions : Neglecting this ubiquitous presence of a broad component towards all coherent cores causes the typical single-component fit to overestimate the temperature and velocity dispersion. This affects the derived detailed physical structure and stability of the cores estimated from $\rm NH_3$ observations., Comment: Accepted for publication in Astronomy & Astrophysics on 06/07/2020. 15 pages, 16 figures, 1 table. Language edits from previous version

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., Comment: 18 pages, 11 figures, accepted for publication in MNRAS

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., Comment: Accepted for publication in A&A on 12/10/19

The majority of stars are part of gravitationally bound stellar systems, such as binaries. Observations of protobinary systems constrain the conditions that lead to stellar multiplicity and subsequent orbital evolution. We report high-angular resolution observations of the circumbinary disk around [BHB2007] 11, a young binary protostar system. The two protostars are embedded in circumstellar disks that have radii of 2 to 3 astronomical units and probably contain a few Jupiter masses. These systems are surrounded by a complex structure of filaments connecting to the larger circumbinary disk. We also observe accretion and radio jets associated with the protobinary system. The accretion is preferentially onto the lower-mass protostar, consistent with theoretical predictions., Comment: 28 pages, 6 figures, accepted for publication in Science

Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometries, independent from current methods of measuring these parameters. To accurately model these features and disentangle the various polarization mechanisms, multiwavelength observations at very high resolution and sensitivity are required. With significant upgrades to current interferometric facilities, we can understand how grains evolve in disks during the planet formation process., Comment: Science white paper submitted to the Astro2020 Decadal Survey

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., Comment: 35 pages, 8 tables, and 14 figures consisting of 16 .pdf files. Accepted for publication in the Astrophysical Journal

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., Comment: Accepted by ApJ in April, 2019

We use polarization observations of a circumbinary disk to investigate how the polarization properties change at distinct frequency bands. Our goal is to discern the main mechanism responsible for the polarization through comparison between our observations and model predictions. We used ALMA to perform full polarization observations at 97.5 GHz, 233 GHz and 343.5 GHz. The target is the Class I object BHB07-11, which is the youngest object in the Barnard 59 protocluster. Complementary VLA observations at 34.5 GHz revealed a binary system within the disk. We detect an extended and structured polarization pattern remarkably consistent among all three bands. The distribution of polarized intensity resembles a horseshoe shape with polarization angles following this morphology. From the spectral index between bands 3 and 7, we derive a dust opacity index $\beta \sim 1$ consistent with maximum grain sizes larger than expected to produce self-scattering polarization in each band. The polarization morphology do not match predictions from self-scattering. On the other hand, marginal correspondence is seen between our maps and predictions from radiation field assuming the brightest binary component as main radiation source. Molecular line data from BHB07-11 indicates disk rotation. We produced synthetic polarization maps from a rotating magnetized disk model assuming combined poloidal and toroidal magnetic field components. The magnetic field vectors (i. e., the polarization vectors rotated by 90$\degr$) are better represented by a model with poloidal magnetic field strength about 3 times the toroidal one. The similarity of our polarization patterns among the three bands provides a strong evidence against self-scattering and radiation fields. On the other hand, our data are reasonably well reproduced by a model of disk with toroidal magnetic field components slightly smaller than poloidal ones., Comment: 8 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Towards the pre-stellar core L1544, the methanol (CH$_3$OH) emission forms an asymmetric ring around the core centre, where CH$_3$OH is mostly in solid form, with a clear peak 4000~au to the north-east of the dust continuum peak. As part of the NOEMA Large Project SOLIS (Seeds of Life in Space), the CH$_3$OH peak has been spatially resolved to study its kinematics and physical structure and to investigate the cause behind the local enhancement. We find that methanol emission is distributed in a ridge parallel to the main axis of the dense core. The centroid velocity increases by about 0.2~km~s$^{-1}$ and the velocity dispersion increases from subsonic to transonic towards the central zone of the core, where the velocity field also shows complex structure. This could be indication of gentle accretion of material onto the core or interaction of two filaments, producing a slow shock. We measure the rotational temperature and show that methanol is in local thermodynamic equilibrium (LTE) only close to the dust peak, where it is significantly depleted. The CH$_3$OH column density, $N_{tot}({\rm CH_3OH})$, profile has been derived with non-LTE radiative transfer modelling and compared with chemical models of a static core. The measured $N_{tot}({\rm CH_3OH})$ profile is consistent with model predictions, but the total column densities are one order of magnitude lower than those predicted by models, suggesting that the efficiency of reactive desorption or atomic hydrogen tunnelling adopted in the model may be overestimated; or that an evolutionary model is needed to better reproduce methanol abundance., Comment: 19 pages, 18 figures. Accepted by ApJ