Characterizing the velocity field in hydrodynamical simulations of low-mass star formation using spectral line profiles

When low-mass stars form, the collapsing cloud of gas and dust goes through several stages which are usually characterized by the shape of their spectral energy distributions. Such classification is based on the cloud morphology only and does not address the dynamical state of the object. In this paper we investigate the initial cloud collapse and subsequent disk formation through the dynamical behavior as reflected in the sub-millimeter spectral emission line profiles. If a young stellar object is to be characterized by its dynamical structure it is important to know how accurately information about the velocity field can be extracted and which observables provide the best description of the kinematics. Of particular interest is the transition from infalling envelope to rotating disk, because this provides the initial conditions for the protoplanetary disk, such as mass and size. We use a hydrodynamical model, describing the collapse of a core and formation of a disk, to produce synthetic observables which we compare to calculated line profiles of a simple parameterized model. Because we know the velocity field from the hydrodynamical simulation we can determine in a quantitative way how well our best-fit parameterized velocity field reproduces the original. We use a molecular line excitation and radiation transfer code to produce spectra of both our hydro dynamical simulation as well as our parameterized model. We find that information about the velocity field can reasonably well be derived by fitting a simple model to either single-dish lines or interferometric data, but preferentially by using a combination of the two. Our result shows that it is possible to establish relative ages of a sample of young stellar objects using this method, independently of the details of the hydrodynamical model.Comment: 12 pages, 11 figures, accepted for publication in A&A on June 1

Chemistry of a newly detected circumbinary disk in Ophiuchus

(Abridged) Astronomers recently started discovering exoplanets around binary systems. Therefore, understanding the formation and evolution of circumbinary disks is crucial for a complete scenario of planet formation. The aim of this paper is to present the detection of a circumbinary disk around Oph-IRS67 and analyse its structure. We present high-angular-resolution (0.4", 60 AU) observations of C17O, H13CO+ , C34S, SO2, C2H and c-C3H2 molecular transitions with ALMA at 0.8 mm. The spectrally and spatially resolved maps reveal the kinematics of the circumbinary disk as well as its chemistry. Molecular abundances are estimated using RADEX. The continuum emission reveals the presence of a circumbinary disk around the two sources. This disk has a diameter of ~620 AU and is well traced by C17O and H13CO+ emission. C2H and c-C3H2 trace a higher-density region which is spatially offset from the sources (~430 AU). Finally, SO2 shows compact emission around one of the sources, Oph-IRS67 B. The molecular transitions which trace the circumbinary disk are consistent with a Keplerian profile on disk scales (< 200 AU) and an infalling profile for envelope scales (> 200 AU). The Keplerian fit leads to a mass of 2.2 Msun. Inferred CO abundances w.r.t. H2 are comparable to the canonical ISM value of 2.7e-4. This study proves the first detection of the circumbinary disk associated with Oph-IRS67. The disk is chemically differentiated from the nearby high-density region. The lack of methanol emission suggests the extended disk dominates the mass budget in the inner- most regions of the protostellar envelope, generating a flat density profile where less material is exposed to high temperatures. Thus, complex organic molecules would be associated with lower column densities. Finally, Oph-IRS67 is a promising candidate for the detection of both circumstellar disks with higher-angular-resolution observations.Comment: 19 pages, 14 figures, 6 table

Global Hydromagnetic Simulations of Protoplanetary Disks with Stellar Irradiation and Simplified Thermochemistry

Outflows driven by large-scale magnetic fields likely play an important role in the evolution and dispersal of protoplanetary disks, and in setting the conditions for planet formation. We extend our 2-D axisymmetric non-ideal MHD model of these outflows by incorporating radiative transfer and simplified thermochemistry, with the twin aims of exploring how heating influences wind launching, and illustrating how such models can be tested through observations of diagnostic spectral lines. Our model disks launch magnetocentrifugal outflows primarily through magnetic tension forces, so the mass-loss rate increases only moderately when thermochemical effects are switched on. For typical field strengths, thermochemical and irradiation heating are more important than magnetic dissipation. We furthermore find that the entrained vertical magnetic flux diffuses out of the disk on secular timescales as a result of non-ideal MHD. Through post-processing line radiative transfer, we demonstrate that spectral line intensities and moment-1 maps of atomic oxygen, the HCN molecule, and other species show potentially observable differences between a model with a magnetically driven outflow and one with a weaker, photoevaporative outflow. In particular, the line shapes and velocity asymmetries in the moment-1 maps could enable the identification of outflows emanating from the disk surface.Comment: 35 pages, 20 figures, accepted for publication in Ap

Profiling Cold New Early Dark Energy

Recent interest in New Early Dark Energy (NEDE), a cosmological model with a vacuum energy component decaying in a triggered phase transition around recombination, has been sparked by its impact on the Hubble tension. Previous constraints on the model parameters were derived in a Bayesian framework with Markov-chain Monte Carlo (MCMC) methods. In this work, we instead perform a frequentist analysis using the profile likelihood in order to assess the impact of prior volume effects on the constraints. We constrain the maximal fraction of NEDE $f_\mathrm{NEDE}$, finding $f_\mathrm{NEDE}=0.076^{+0.040}_{-0.035}$ at $68 \%$ CL with our baseline dataset and similar constraints using either data from SPT-3G, ACT or full-shape large-scale structure, showing a preference over $\Lambda$CDM even in the absence of a SH0ES prior on $H_0$. While this is stronger evidence for NEDE than obtained with the corresponding Bayesian analysis, our constraints broadly match those obtained by fixing the NEDE trigger mass. Including the SH0ES prior on $H_0$, we obtain $f_\mathrm{NEDE}=0.136^{+0.024}_{-0.026}$ at $68 \%$ CL. Furthermore, we compare NEDE with the Early Dark Energy (EDE) model, finding similar constraints on the maximal energy density fractions and $H_0$ in the two models. At $68 \%$ CL in the NEDE model, we find $H_0 = 69.56^{+1.16}_{-1.29} \text{ km s}^{-1}\text{ Mpc}^{-1}$ with our baseline and $H_0 = 71.62^{+0.78}_{-0.76} \text{ km s}^{-1}\text{ Mpc}^{-1}$ when including the SH0ES measurement of $H_0$, thus corroborating previous conclusions that the NEDE model provides a considerable alleviation of the $H_0$ tension.Comment: 14 pages, 7 figures, 1 tabl

Resolving the shocked gas in HH54 with Herschel: CO line mapping at high spatial and spectral resolution

The HH54 shock is a Herbig-Haro object, located in the nearby Chamaeleon II cloud. Observed CO line profiles are due to a complex distribution in density, temperature, velocity, and geometry. Resolving the HH54 shock wave in the far-infrared cooling lines of CO constrain the kinematics, morphology, and physical conditions of the shocked region. We used the PACS and SPIRE instruments on board the Herschel space observatory to map the full FIR spectrum in a region covering the HH54 shock wave. Complementary Herschel-HIFI, APEX, and Spitzer data are used in the analysis as well. The observed features in the line profiles are reproduced using a 3D radiative transfer model of a bow-shock, constructed with the Line Modeling Engine code (LIME). The FIR emission is confined to the HH54 region and a coherent displacement of the location of the emission maximum of CO with increasing J is observed. The peak positions of the high-J CO lines are shifted upstream from the lower J CO lines and coincide with the position of the spectral feature identified previously in CO(10-9) profiles with HIFI. This indicates a hotter molecular component in the upstream gas with distinct dynamics. The coherent displacement with increasing J for CO is consistent with a scenario where IRAS12500-7658 is the exciting source of the flow, and the 180 K bow-shock is accompanied by a hot (800 K) molecular component located upstream from the apex of the shock and blueshifted by -7 km s$^{-1}$. The spatial proximity of this knot to the peaks of the atomic fine-structure emission lines observed with Spitzer and PACS ([OI]63, 145 $\mu$m) suggests that it may be associated with the dissociative shock as the jet impacts slower moving gas in the HH54 bow-shock.Comment: 6 pages, 5 figure

Refueling Infrastructure for Alternative Fuel Vehicles: Lessons Learned for Hydrogen; Workshop Proceedings

DOE sponsored the Refueling Infrastructure for Alternative Fuel Vehicles: Lessons Learned for Hydrogen workshop to understand how lessons from past experiences can inform future efforts to commercialize hydrogen vehicles. This report contains the proceedings from the workshop

Dimethyl ether in its ground state, v=0, and lowest two torsionally excited states, v11=1 and v15=1, in the high-mass star-forming region G327.3-0.6

The goal of this paper is to determine the respective importance of solid state vs. gas phase reactions for the formation of dimethyl ether. This is done by a detailed analysis of the excitation properties of the ground state and the torsionally excited states, v11=1 and v15=1, toward the high-mass star-forming region G327.3-0.6. With the Atacama Pathfinder EXperiment 12 m submillimeter telescope, we performed a spectral line survey. The observed spectrum is modeled assuming local thermal equilibrium. CH3OCH3 has been detected in the ground state, and in the torsionally excited states v11=1 and v15=1, for which lines have been detected here for the first time. The emission is modeled with an isothermal source structure as well as with a non-uniform spherical structure. For non-uniform source models one abundance jump for dimethyl ether is sufficient to fit the emission, but two components are needed for the isothermal models. This suggests that dimethyl ether is present in an extended region of the envelope and traces a non-uniform density and temperature structure. Both types of models furthermore suggest that most dimethyl ether is present in gas that is warmer than 100 K, but a smaller fraction of 5%-28% is present at temperatures between 70 and 100 K. The dimethyl ether present in this cooler gas is likely formed in the solid state, while gas phase formation probably is dominant above 100 K. Finally, the v11=1 and v15=1 torsionally excited states are easily excited under the density and temperature conditions in G327.3-0.6 and will thus very likely be detectable in other hot cores as well.Comment: 12 pages (excluding appendices), 8 figures, A&A in pres

Determining the Parameters of Massive Protostellar Clouds via Radiative Transfer Modeling

A one-dimensional method for reconstructing the structure of prestellar and protostellar clouds is presented. The method is based on radiative transfer computations and a comparison of theoretical and observed intensity distributions at both millimeter and infrared wavelengths. The radiative transfer of dust emission is modeled for specified parameters of the density distribution, central star, and external background, and the theoretical distribution of the dust temperature inside the cloud is determined. The intensity distributions at millimeter and IR wavelengths are computed and quantitatively compared with observational data. The best-fit model parameters are determined using a genetic minimization algorithm, which makes it possible to reveal the ranges of parameter degeneracy as well. The method is illustrated by modeling the structure of the two infrared dark clouds IRDC-320.27+029 (P2) and IRDC-321.73+005 (P2). The derived density and temperature distributions can be used to model the chemical structure and spectral maps in molecular lines.Comment: Accepted for publication in Astronomy Report

Modelling Herschel observations of hot molecular gas emission from embedded low-mass protostars

Aims. Young stars interact vigorously with their surroundings, as evident from the highly rotationally excited CO (up to Eup=4000 K) and H2O emission (up to 600 K) detected by the Herschel Space Observatory in embedded low-mass protostars. Our aim is to construct a model that reproduces the observations quantitatively, to investigate the origin of the emission, and to use the lines as probes of the various heating mechanisms. Methods. The model consists of a spherical envelope with a bipolar outflow cavity. Three heating mechanisms are considered: passive heating by the protostellar luminosity, UV irradiation of the outflow cavity walls, and C-type shocks along the cavity walls. Line fluxes are calculated for CO and H2O and compared to Herschel data and complementary ground-based data for the protostars NGC1333 IRAS2A, HH 46 and DK Cha. The three sources are selected to span a range of evolutionary phases and physical characteristics. Results. The passively heated gas in the envelope accounts for 3-10% of the CO luminosity summed over all rotational lines up to J=40-39; it is best probed by low-J CO isotopologue lines such as C18O 2-1 and 3-2. The UV-heated gas and the C-type shocks, probed by 12CO 10-9 and higher-J lines, contribute 20-80% each. The model fits show a tentative evolutionary trend: the CO emission is dominated by shocks in the youngest source and by UV-heated gas in the oldest one. This trend is mainly driven by the lower envelope density in more evolved sources. The total H2O line luminosity in all cases is dominated by shocks (>99%). The exact percentages for both species are uncertain by at least a factor of 2 due to uncertainties in the gas temperature as function of the incident UV flux. However, on a qualitative level, both UV-heated gas and C-type shocks are needed to reproduce the emission in far-infrared rotational lines of CO and H2O.Comment: 15 pages (+4 pages appendix), 20 figures, accepted by A&