No Arabic abstract
We present a model aimed to reproduce the observed spectral energy distribution (SED) as well as the ammonia line emission of the G31.41+0.31 hot core. The core is modeled as an infalling envelope onto a massive star that is undergoing an intense accretion phase. We assume an envelope with a density and velocity structure resulting from the dynamical collapse of a singular logatropic sphere. The stellar and envelope physical properties are determined by fitting the SED. From these physical conditions, the ammonia line emission is calculated and compared with subarcsecond resolution VLA data of the (4,4) transition. The only free parameter in this line fitting is the ammonia abundance. The observed properties of the NH3(4,4) lines and their spatial distribution can be well reproduced provided it is taken into account the steep increase of the abundance in the hotter (> 100 K), inner regions of the core produced by the sublimation of icy mantles where ammonia molecules are trapped. The model predictions for the (2,2), (4,4), and (5,5) transitions are also in reasonably agreement with the single-dish spectra available in the literature. The best fit is obtained for a model with a star of 25 Msun, a mass accretion rate of 0.003 Msun/yr, and a total luminosity of 200,000 Lsun. The gas-phase ammonia abundance ranges from 2 times 10^{-8} in the outer region to 3 times 10^{-6} in the inner region. To our knowledge, this is the first time that the dust and molecular line data of a hot molecular core, including subarcsecond resolution data that spatially resolve the structure of the core, have been simultaneously explained by a physically self-consistent model. This modeling shows that massive protostars are able to excite high excitation ammonia transitions up to the outer edge (30,000 AU) of the large scale envelope.
As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22 (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.
G31.41+0.31 is a well known chemically rich hot molecular core (HMC). Using Band 3 observations of Atacama Large Millimeter Array (ALMA), we have analyzed the chemical and physical properties of the source. We have identified methyl isocyanate (CH3NCO), a precursor of prebiotic molecules, towards the source. In addition to this, we have reported complex organic molecules (COMs) like methanol (CH3OH), methanethiol (CH3SH), and methyl formate (CH3OCHO). Additionally, we have used transitions from molecules like HCN, HCO+, SiO to trace the presence of infall and outflow signatures around the star-forming region. For the COMs, we have estimated the column densities and kinetic temperatures, assuming molecular excitation under local thermodynamic equilibrium (LTE) conditions. From the estimated kinetic temperatures of certain COMs, we found that multiple temperature components may be present in the HMC environment. Comparing the obtained molecular column densities between the existing observational results toward other HMCs, it seems that the COMs are favourably produced in the hot-core environment ($sim 100$ K or higher). Though the spectral emissions towards G31.41+0.31 are not fully resolved, we find that CH$_3$NCO and other COMs are possibly formed on the grain/ice phase and populate the gas environment similar to other hot cores like Sgr B2, Orion KL, and G10.47+0.03, etc.
An inverse P-Cygni profile of H13CO+ (1-0) in G31.41+0.31 was recently observed, which indicates the presence of an infalling gas envelope. Also, an outflow tracer, SiO, was observed. Here, exclusive radiative transfer modelings have been implemented to generate synthetic spectra of some key species (H13 CO+, HCN, SiO, NH3, CH3 CN, CH3OH, CH3SH, and CH3NCO) and extract the physical features to infer the excitation conditions of the surroundings where they observed. The gas envelope is assumed to be accreting in a spherically symmetric system towards the central hot core region. Our principal intention was to reproduce the observed line profiles toward G31.41+0.31 and extract various physical parameters. The LTE calculation with CASSIS and non-LTE analysis with the RATRAN radiative transfer codes are considered for the modeling purpose. The best-fitted line parameters are derived, which represents the prevailing physical condition of the gas envelope. Our results suggest that an infalling gas could explain the observed line profiles of all the species mentioned above except SiO. An additional outflow component is required to confer the SiO line profile. Additionally, an astrochemical model is implemented to explain the observed abundancests various species in this source.
We study the origin of large abundances of complex organic molecules in the Galactic center (GC). We carried out a systematic study of the complex organic molecules CH3OH, C2H5OH, (CH3)2O, HCOOCH3, HCOOH, CH3COOH, H2CO, and CS toward 40 GC molecular clouds. Using the LTE approximation, we derived the physical properties of GC molecular clouds and the abundances of the complex molecules.The CH3OH abundance between clouds varies by nearly two orders of magnitude from 2.4x10^{-8} to 1.1x10^{-6}. The abundance of the other complex organic molecules relative to that of CH3OH is basically independent of the CH3OH abundance, with variations of only a factor 4-8. The abundances of complex organic molecules in the GC are compared with those measured in hot cores and hot corinos, in which these complex molecules are also abundant. We find that both the abundance and the abundance ratios of the complex molecules relative to CH3OH in hot cores are similar to those found in the GC clouds. However, hot corinos show different abundance ratios than observed in hot cores and in GC clouds. The rather constant abundance of all the complex molecules relative to CH3OH suggests that all complex molecules are ejected from grain mantles by shocks. Frequent (similar 10^{5}years) shocks with velocities >6km/s are required to explain the high abundances in gas phase of complex organic molecules in the GC molecular clouds. The rather uniform abundance ratios in the GC clouds and in Galactic hot cores indicate a similar average composition of grain mantles in both kinds of regions. The Sickle and the Thermal Radio Arches, affected by UV radiation, show different relative abundances in the complex organic molecules due to the differentially photodissociation of these molecules.
Context. ALMA observations at 1.4 mm and 0.2 (750au) angular resolution of the Main core in the high-mass star forming region G31.41+0.31 have revealed a puzzling scenario: on the one hand, the continuum emission looks very homogeneous and the core appears to undergo solid-body rotation, suggesting a monolithic core stabilized by the magnetic field; on the other hand, rotation and infall speed up toward the core center, where two massive embedded free-free continuum sources have been detected, pointing to an unstable core having undergone fragmentation. Aims. To establish whether the Main core is indeed monolithic or its homogeneous appearance is due to a combination of large dust opacity and low angular resolution, we carried out millimeter observations at higher angular resolution and different wavelengths. Methods. We carried out ALMA observations at 1.4 mm and 3.5 mm that achieved angular resolutions of 0.1(375 au) and 0.075 (280 au), respectively. VLA observations at 7 mm and 1.3 cm at even higher angular resolution, 0.05 (190 au) and 0.07 (260 au), respectively, were also carried out to better study the nature of the free-free continuum sources detected in the core. Results. The millimeter continuum emission of the Main core has been clearly resolved into at least four sources, A, B, C, and D, within 1, indicating that the core is not monolithic. The deconvolved radii of the dust emission of the sources, estimated at 3.5 mm, are 400-500au, their masses range from 15 to 26 Msun, and their number densities are several 1E9 cm-3. Sources A and B, located closer to the center of the core and separated by 750 au, are clearly associated with two free-free continuum sources, likely thermal radio jets, and are the brightest in the core. The spectral energy distribution of these two sources and their masses and sizes are similar and suggest a common origin.