Search

We study the structure, interstellar absorption, color-magnitude diagrams, kinematics, and dynamical state of embedded star clusters in the star-forming region associated with the giant molecular cloud G174+2.5. Our investigation is based on photometric data from the UKIDSS Galactic Plane Survey catalog and astrometric data from the Gaia DR3 catalogs. First, we recover all the known embedded clusters and candidate clusters in the region using surface density maps. Then, for the detected clusters, we determine their general parameters: the center positions, radii, number of stars, and reddening. To evaluate the reddening, we use both the NICEST algorithm and the Q-method. Both methods produce consistent extinction maps in the regions of the four studied clusters. However, the Q-method yields a much smaller color scatter in the CMD. For four clusters in particular (S235~North-West, S235~A-B-C, S235~Central, and S235~East1+East2), we were able to compute individual membership probabilities, the cluster distances, the cluster masses, and their average proper motions. By building on these results, we have studied the clusters' kinematics and dynamics. Moreover, we estimate the mass of the gas component and the star formation efficiency (SFE) in the regions of these four clusters. Finally, we provide an estimate of the total energy of the stellar and gas components in the area of these four clusters to determine whether the clusters are bound (here we consider a `cluster' as the system `stars + gas'). The gravitational bound strongly depends on the region for which we estimate the gas mass. If we consider the mass of the entire cloud, all these four clusters turn out to be bound., Comment: 30 pages, 16 figures, accepted for publication in the Astronomical Journal

HMYSOs gain most of their mass in short bursts of accretion. Maser emission is an invaluable tool in discovering and probing accretion bursts. We observed the 22 GHz water maser response induced by the accretion burst in NGC6334I-MM1B and identified the underlying maser variability mechanisms. We report seven epochs of VLBI observations of 22 GHz water masers in NGC6334I with the VERA array, from 2014 to 2016, spanning the onset of the accretion burst in 2015.1. We also report 2019 ALMA observations of 321 GHz water masers and 22 GHz maser monitoring by HartRAO. We analyze variability patterns and use proper motions with the 22 GHz to 321 GHz line ratio to distinguish between masers in C-shocks and J-shocks. We also calculated the burst-to-quiescent variance ratio of the single-dish time series. The constant mean proper motion before and after the burst indicates that maser variability is due to excitation effects from variable radiation rather than jet ejecta. We find that the flux density variance ratio in the single-dish time series can identify maser efficiency variations in 22 GHz masers. The northern region, CM2-W2, is excited in C-shocks and showed long-term flaring with velocity-dependent excitation of new maser features. We propose that radiative heating of H2 due to high-energy radiation from the accretion burst be the mechanism for the flaring in CM2-W2. The southern regions are excited by J-shocks and have short-term flaring and dampening of water masers. We attributed the diverse variability patterns in the southern regions to the radiative transfer of the burst energy in the source. Our results indicate that the effects of source geometry, shock type, and incident radiation spectrum are fundamental factors affecting 22 GHz maser variability. Investigating water masers in irradiated shocks will improve their use as a diagnostic in time-variable radiation environments., Comment: 21 pages, 12 figures, Table 2 in electronic publication. Contact: jakobus.vorster@helsinki.fi

High-mass protostars (M$_{\star} >$ 8 M$_{\odot}$) are thought to gain the majority of their mass via short, intense bursts of growth. This episodic accretion is thought to be facilitated by gravitationally unstable and subsequently inhomogeneous accretion disks. Limitations of observational capabilities, paired with a lack of observed accretion burst events has withheld affirmative confirmation of the association between disk accretion, instability and the accretion burst phenomenon in high-mass protostars. Following its 2019 accretion burst, a heat-wave driven by a burst of radiation propagated outward from the high-mass protostar G358.93-0.03-MM1. Six VLBI (very long baseline interferometry) observations of the raditively pumped 6.7 GHz methanol maser were conducted during this period, tracing ever increasing disk radii as the heat-wave propagated outward. Concatenating the VLBI maps provided a sparsely sampled, milliarcsecond view of the spatio-kinematics of the accretion disk covering a physical range of $\sim$ 50 - 900 AU. We term this observational approach `heat-wave mapping'. We report the discovery of a Keplerian accretion disk with a spatially resolved four-arm spiral pattern around G358.93-0.03-MM1. This result positively implicates disk accretion and spiral arm instabilities into the episodic accretion high-mass star formation paradigm., Comment: Published in Nature Astronomy in 2023

High-mass stars are thought to accumulate much of their mass via short, infrequent bursts of disk-aided accretion. Such accretion events are rare and difficult to observe directly but are known to drive enhanced maser emission. In this Letter we report high-resolution, multi-epoch methanol maser observations toward G358.93-0.03 which reveal an interesting phenomenon; the sub-luminal propagation of a thermal radiation "heat-wave" emanating from an accreting high-mass proto-star. The extreme transformation of the maser emission implies a sudden intensification of thermal infrared radiation from within the inner (40 mas, 270 au) region. Subsequently, methanol masers trace the radial passage of thermal radiation through the environment at $\geq$ 4-8\% the speed of light. Such a high translocation rate contrasts with the $\leq$ 10 km s$^{-1}$ physical gas motions of methanol masers typically observed using very long baseline interferometry (VLBI). The observed scenario can readily be attributed to an accretion event in the high-mass proto-star G358.93-0.03-MM1. While being the third case in its class, G358.93-0.03-MM1 exhibits unique attributes hinting at a possible `zoo' of accretion burst types. These results promote the advantages of maser observations in understanding high-mass star formation, both through single-dish maser monitoring campaigns and via their international cooperation as VLBI arrays., Comment: Published in Nature Astronomy in 2020

We report the detection of ammonia masers in the non-metastable (6, 3), (7, 5) and (6, 5) transitions, the latter is the first unambiguous maser detection of that transition ever made. Our observations include the first VLBI detection of ammonia maser emission, which allowed effective constrain of the (6, 5) maser brightness temperature. The masers were detected towards G358.931-0.030, a site of 6.7-GHz class~II methanol maser emission that was recently reported to be undergoing a period of flaring activity. These ammonia masers appear to be flaring contemporaneously with the class~II methanol masers during the accretion burst event of G358.931-0.030. This newly detected site of ammonia maser emission is only the twelfth such site discovered in the Milky Way. We also report the results of an investigation into the maser pumping conditions, for all three detected masing transitions, through radiative transfer calculations constrained by our observational data. These calculations support the hypothesis that the ammonia (6, 5) maser transition is excited through high colour temperature infrared emission, with the (6, 5) and (7, 5) transition line-ratio implying dust temperatures >400K. Additionally, we detect significant linearly polarised emission from the ammonia (6, 3) maser line. Alongside our observational and radiative transfer calculation results, we also report newly derived rest frequencies for the ammonia (6, 3) and (6, 5) transitions., Comment: Accepted into MNRAS 2023 April 24. 13 pages, 5 figures, 3 tables

The burst mode of accretion in massive star formation is a scenario linking the initial gravitational collapse of parent pre-stellar cores to the properties of their gravitationally unstable discs and of their accretion-driven bursts. In this study, we present a series of high-resolution 3D radiation-hydrodynamics numerical simulations for young massive stars formed out of collapsing 100 Mo molecular cores spinning with several values of the ratio of rotational-to-gravitational energies beta=5%-9%. The models include the indirect gravitational potential caused by disc asymmetries. We find that this modifies the barycenter of the disc, causing significant excursions of the central star position, which we term stellar wobbling. The stellar wobbling slows down and protracts the development of gravitational instability in the disc, reducing the number and magnitude of the accretion-driven bursts undergone by the young massive stars, whose properties are in good agreement with that of the burst monitored from the massive protostar M17 MIR. Including stellar wobbling is therefore important for accurate modeling disc structures. Synthetic ALMA interferometric images in the millimeter waveband show that the outcomes of efficient gravitational instability such as spiral arms and gaseous clumps can be detected for as long as the disc is old enough and has already entered the burst mode of accretion., Comment: Accepted to MNRAS, see https://academic.oup.com/mnras/advance-article-abstract/doi/10.1093/mnras/stac2956/6761717?utm_source=advanceaccess&utm_campaign=mnras&utm_medium=email

We report VLBI monitoring observations of the 22 GHz water (H$_{2}$O) masers around the Mira variable BX Cam, which were carried out as a part of the EAVN Synthesis of Stellar Maser Animations (ESTEMA) project. Data of 37 epochs in total were obtained from 2018 May to 2021 June with a time interval of 3-4 weeks, spanning approximately three stellar pulsation periods ($P= \sim$440 d). In particular, the dual-beam system equipped on the VERA stations was used to measure the kinematics and parallaxes of the H$_{2}$O maser features. The measured parallax, $\pi=1.79\pm 0.08$ mas, is consistent with $Gaia$ EDR3 and previously measured VLBI parallaxes within a 1-$\sigma$ error level. The position of the central star was estimated, based on both the $Gaia$ EDR3 data and the center position of the ring-like 43 GHz silicon-monoxide (SiO) maser distribution imaged with the KVN. The three-dimensional H$_{2}$O maser kinematics indicates that the circumstellar envelope is expanding at a velocity of $13\pm4$ km s$^{-1}$, while there are asymmetries in both the spatial and velocity distributions of the maser features. Furthermore, the H$_{2}$O maser animation achieved by our dense monitoring program manifests the propagation of shock waves in the circumstellar envelope of BX Cam., Comment: 31 pages, 16 figures, accepted by ApJ

The massive young stellar object (MYSO) G358.93-0.03-MM1 showed an extraordinary near-infrared- to (sub-)millimetre-dark and far-infrared-loud accretion burst, which is closely associated with flares of several class II methanol maser transitions, and, later, a 22 GHz water maser flare. Water maser flares provide an invaluable insight into ejection events associated with accretion bursts. Although the short timescale of the 22 GHz water maser flare made it impossible to carry out a very long baseline interferometry observation, we could track it with the Karl G. Jansky Very Large Array (VLA). The evolution of the spatial structure of the 22 GHz water masers and their association with the continuum sources in the region is studied with the VLA during two epochs, pre- and post-H2O maser flare. A drastic change in the distribution of the water masers is revealed: in contrast to the four maser groups detected during epoch I, only two newly formed clusters are detected during epoch II. The 22 GHz water masers associated with the bursting source MM1 changed in morphology and emission velocity extent. Clear evidence of the influence of the accretion burst on the ejection from G358.93-0.03-MM1 is presented. The accretion event has also potentially affected a region with a radius of ~2'' (~13 500 AU at 6.75 kpc), suppressing water masers associated with other point sources in this region.

We present results of the investigations of the properties of the methanol $J_1 - J_0$ A$^{-+}$ line series motivated by the recent serendipitous detection of the maser emission in the $14_1 - 14_0$ A$^{-+}$ line at 349~GHz in S255IR-SMA1 soon after the accretion burst. The study includes further observations of several lines of this series in S255IR with the SMA, a mini-survey of methanol lines in the 0.8 mm range toward a sample of bright 6.7~GHz methanol maser sources with the IRAM-30m telescope, and theoretical modeling. We found that the maser component of the $14_1 - 14_0$ A$^{-+}$ line in S255IR decayed by more than order of magnitude in comparison with that in 2016. No clear sign of maser emission is observed in other lines of this series in the SMA observations except the $7_1 - 7_0$ A$^{-+}$ line where an additional bright component is detected at the velocity of the maser emission observed earlier in the $14_1 - 14_0$ A$^{-+}$ line. Our LVG model constrains the ranges of the physical parameters that matches the observed emission intensities. No obvious maser emission in the $J_1 - J_0$ A$^{-+}$ lines was detected in the mini-survey of the 6.7~GHz methanol maser sources, though one component in NGC7538 may represent a weak maser. In general, the maser effect in the $J_1 - J_0$ A$^{-+}$ lines may serve as a tracer of rather hot environments and in particular luminosity flaring events during high mass star formation., Comment: 16 pages, 12 figures

We present the most complete to date interferometric study of the centimeter wavelength methanol masers detected in G358.93-0.03 at the burst and post-burst epochs. A unique, NIR/(sub)mm-dark and FIR-loud MYSO accretion burst was recently discovered in G358.93-0.03. The event was accompanied by flares of an unprecedented number of rare methanol maser transitions. The first images of three of the newly-discovered methanol masers at 6.18, 12.23, and 20.97 GHz are presented in this work. The spatial structure evolution of the methanol masers at 6.67, 12.18, and 23.12 GHz is studied at two epochs. The maser emission in all detected transitions resides in a region of $\sim$0.2$^{\prime\prime}$ around the bursting source and shows a clear velocity gradient in the north-south direction, with red-shifted features to the north and blue-shifted features to the south. A drastic change in the spatial morphology of the masing region is found: a dense and compact "spiral" cluster detected at epoch I evolved into a disperse, "round" structure at epoch II. During the transition from the first epoch to the second, the region traced by masers expanded. The comparison of our results with the complementary VLA, VLBI, SMA, and ALMA maser data is conducted. The obtained methanol maser data support the hypothesis of the presence of spiral-arm structures within the accretion disk, which was suggested in previous studies of the source.

The paper aims to study relation between the distributions of the young stellar objects (YSOs) of different ages and the gas-dust constituents of the S254-S258 star-formation complex. This is necessary to study the time evolution of the YSO distribution with respect to the gas and dust compounds which are responsible for the birth of the young stars. For this purpose we use correlation analysis between different gas, dust and YSOs tracers. We compared the large-scale CO, HCO$^+$, near-IR extinction, and far-IR {\it Herschel} maps with the density of YSOs of the different evolutionary Classes. The direct correlation analysis between these maps was used together with the wavelet-based spatial correlation analysis. This analysis reveals a much tighter correlation of the gas-dust tracers with the distribution of Class I YSOs than with that of Class II YSOs. We argue that Class I YSOs which were initially born in the central bright cluster S255-IR (both N and S parts) during their evolution to Class II stage ($\sim$2 Myr) had enough time to travel through the whole S254-S258 star-formation region. Given that the region contains several isolated YSO clusters, the evolutionary link between these clusters and the bright central S255-IR (N and S) cluster can be considered. Despite the complexity of the YSO cluster formation in the non-uniform medium, the clusters of Class II YSOs in the S254-258 star-formation region can contain objects born in the different locations of the complex., Comment: Accepted for publication in MNRAS

We study molecular emission in a massive condensation at the border of the HII region RCW 120, paying particular attention to the Core 1 and Core 2 objects, the most massive fragments of the condensation found previously by ALMA. The latter fragment was previously suggested to host a high-mass analogue of Class 0 young stellar object. We present spectra of molecular emission in the 1 mm range made with the APEX telescope. We detect CH$_3$OH and C$^{34}$S lines in Core 1 and Core 2. The CH$_3$CN series and the SO$_2$ lines are only found in Core 2. We estimate gas physical parameters using methanol lines and obtain gas temperature less than 100 K in both regions. Molecular hydrogen number density in Core 2 is in the range of $10^5-10^7$ cm$^{-3}$ and is more uncertain in Core 1. However, the detection of the CH$_3$CN lines corresponding to highly excited transitions ($E_{\rm u}> 400$~K) in Core~2 indicates that the region contains hot gas, while the abundances of CH$_3$OH, CS, SO$_2$ and CH$_3$CN are quite low for a hot core stage. We propose that Core 2 is in the warm-up phase prior to the establishing of the hot gas chemistry. We suggest that Core 2 is in the beginning of the hot core stage. There are no detected CH$_3$CN lines in Core 1, therefore, it might be on an even less evolved evolutionary stage., Comment: 11 pages, 6 figures, accepted by MNRAS

It is now a widely held view that, in their formation and early evolution, stars build up mass in bursts. The burst mode of star formation scenario proposes that the stars grow in mass via episodic accretion of fragments migrating from their gravitationally-unstable circumstellar discs and it naturally explains the existence of observed pre-main-sequence bursts from high mass protostars. We present a parameter study of hydrodynamical models of massive young stellar objects (MYSOs) that explores the initial masses of the collapsing clouds (Mc = 60-200Mo) and ratio of rotational-to-gravitational energies (beta = 0:005-0:33). An increase in Mc and/or beta produces protostellar accretion discs that are more prone to develop gravitational instability and to experience bursts. We find that all MYSOs have bursts even if their pre-stellar core is such that beta <= 0.01. Within our assumptions, the lack of stable discs is therefore a major difference between low- and high-mass star formation mechanisms. All our disc masses and disk-to-star mass ratios Md=M* > 1 scale as a power-law with the stellar mass. Our results confirm that massive protostars accrete about 40-60% of their mass in the burst mode. The distribution of time periods between two consecutive bursts is bimodal: there is a short duration (~ 1-10 yr) peak corresponding to the short, faintest bursts and a long duration peak (at ~ 10^3-10^4 yr) corresponding to the long, FU-Orionis-type bursts appearing in later disc evolution, i.e., around 30 kyr after disc formation. We discuss this bimodality in the context of the structure of massive protostellar jets as potential signatures of accretion burst history., Comment: 22 pages, 11 figures. Accepted at MNRAS

We present a study of the association between class I methanol masers and cold dust clumps from the ATLASGAL survey. It was found that almost 100% of class I methanol masers are associated with objects listed in the ATLASGAL compact source catalog. We find a statistically significant difference in the flux density, luminosity, number and column density and temperature distributions of ATLASGAL sources associated with 95/44 GHz methanol masers compared with those ATLASGAL sources devoid of 95 GHz methanol masers. The masers tend to arise in clumps with higher densities, luminosities and temperatures compared with both the full sample of the ATLASGAL clumps, as well as the sample of ATLASGAL sources that were cross-matched with positions previously searched for methanol masers but with no detections. Comparison between the peak position of ATLASGAL clumps and the interferometric positions of the associated class I and II methanol masers reveals that class I masers are generally located at larger physical distances from the peak submillimetre emission than class II masers. We conclude that the tight association between ATLASGAL sources and class I methanol masers may be used as a link toward understanding the conditions of the pumping of these masers and evolutionary stages at which they appear., Comment: Accepted to the Astronomical Journal, 14 pages, 7 figures

The aim of the present work is to study structure and gas kinematics in the photodissociation regions (PDRs) around the compact HII regions S235A and S235C. We observe the [CII], [13CII] and [OI] line emission, using SOFIA/upGREAT and complement them by data of HCO+ and CO. We use the [13CII] line to measure the optical depth of the [CII] emission, and find that the [CII] line profiles are influenced by self-absorption, while the [13CII] line remains unaffected by these effects. Hence, for dense PDRs, [13CII] emission is a better tracer of gas kinematics. The optical depth of the [CII] line is up to 10 in S235A. We find an expanding motion of the [CII]-emitting layer of the PDRs into the front molecular layer in both regions. Comparison of the gas and dust columns shows that gas components visible neither in the [CII] nor in low-J CO lines may contribute to the total column across S235A. We test whether the observed properties of the PDRs match the predictions of spherical models of expanding HII region + PDR + molecular cloud. Integrated intensities of the [13CII], [CII] and [OI] lines are well-represented by the model, but the models do not reproduce the double-peaked [CII] line profiles due to an insufficient column density of C+. The model predicts that the [OI] line could be a more reliable tracer of gas kinematics, but the foreground self-absorbing material does not allow using it in the considered regions., Comment: accepted to MNRAS

We report a simultaneous 44 and 95 GHz class I methanol maser survey toward 144 sources from the 95 GHz class I methanol maser catalog. The observations were made with the three telescopes of the Korean very long baseline interferometry network operating in single-dish mode. The detection rates are 89% at 44 GHz and 77% at 95 GHz. There are 106 new discoveries at 44 GHz. Comparing the previous 95 GHz detections with new observations of the same transitions made using the Purple Mountain Observatory 13.7 m radio telescope shows no clear evidence of variability on a timescale of six years. Emission from the 44 and 95 GHz transitions shows strong correlations in peak velocity, peak flux density, and integrated flux density, indicating that they are likely cospatial. We found that the peak flux density ratio Spk,95/Spk,44 decreases as the 44 GHz peak flux density increases. We found that some class I methanol masers in our sample could be associated with infrared dark clouds, while others are associated with H II regions, indicating that some sources occur at an early stage of high-mass star formation, while others are located toward more evolved sources., Comment: 26 pages, including 10 figures and 8 tables. Accepted for publication in ApJS

This paper reports observations of a 22 GHz water maser `superburst' in the G25.65+1.05 massive star forming region, conducted in response to an alert from the Maser Monitoring Organisation (M2O). Very long baseline interferometry (VLBI) observations using the European VLBI Network (EVN) recorded a maser flux density of $1.2 \times 10^{4}$ Jy. The superburst was investigated in the spectral, structural and temporal domains and its cause was determined to be an increase in maser path length generated by the superposition of multiple maser emitting regions aligning in the line of sight to the observer. This conclusion was based on the location of the bursting maser in the context of the star forming region, its complex structure, and its rapid onset and decay.

We report the detection of new 12.178, 12.229, 20.347, and 23.121 GHz methanol masers in the massive star-forming region G358.93-0.03, which are flaring on similarly short timescales (days) as the 6.668 GHz methanol masers also associated with this source. The brightest 12.178 GHz channel increased by a factor of over 700 in just 50 d. The masers found in the 12.229 and 20.347 GHz methanol transitions are the first ever reported and this is only the fourth object to exhibit associated 23.121 GHz methanol masers. The 12.178 GHz methanol maser emission appears to have a higher flux density than that of the 6.668 GHz emission, which is unusual. No associated near-infrared flare counterpart was found, suggesting that the energy source of the flare is deeply embedded., Comment: 9 pages, 3 figures, 2 tables, Published in MNRAS

The H II region RCW120 is a well-known object, which is often considered as a target to verify theoretical models of gas and dust dynamics in the interstellar medium. However, the exact geometry of RCW120 is still a matter of debate. In this work, we analyse observational data on molecular emission in RCW120 and show that 13CO(2-1) and C18O(2-1) lines are fitted by a 2D model representing a ring-like face-on structure. The changing of the C18O(3-2) line profile from double-peaked to single-peaked from the dense molecular Condensation 1 might be a signature of stalled expansion in this direction. In order to explain a self-absorption dip of the 13CO(2-1) and 13CO(3-2) lines, we suggest that RCW120 is surrounded by a diffuse molecular cloud, and find confirmation of this cloud on a map of interstellar extinction. Optically thick 13CO(2-1) emission and the infrared 8 um PAH band form a neutral envelope of the H II region resembling a ring, while the envelope breaks into separate clumps on images made with optically thin C18O(2-1) line and far-infrared dust emission., Comment: accepted to MNRAS

Star formation is thought to be driven by two groups of mechanisms; spontaneous collapse and triggered collapse. Triggered star formation mechanisms further diverge into cloud-cloud collision (CCC), "collect and collapse" (C\&C) and shock induced collapse of pre-existing, gravitationally stable cores, or 'radiation driven implosion' (RDI). To evaluate the contributions of these mechanisms and establish whether these processes can occur together within the same star forming region we performed mapping observations of radio frequency ammonia, and water maser emission lines in the S235 massive star forming region. Via spectral analyses of main, hyperfine and multi-transitional ammonia lines we explored the distribution of temperature and column density in the dense gas in the S235 and S235AB star forming region. The most remarkable result of the mapping observations is the discovery of high density gas in inter-core bridges which physically link dense molecular cores that house young proto-stellar clusters. The presence of dense gas implies the potential for future star formation within the system of cores and gas bridges. Cluster formation implies collapse and the continuous physical links, also seen in re-imaged archival CS and $^{13}$CO maps, suggests a common origin to the molecular cores housing these clusters, i.e the structure condensed from a single, larger parent cloud, brought about by the influence of a local expanding H${\rm II}$ region. An ammonia absorption feature co-locating with the center of the extended H${\rm II}$ region may be attributed to an older gas component left over from the period prior to formation of the H${\rm II}$ region. Our observations also detail known and new sites of water maser emission, highlighting regions of active ongoing star formation., Comment: 13 pages, 9 figures, accepted for publication in PASJ

We present (sub)millimeter imaging at 0.5'' resolution of the massive star-forming region G358.93-0.03 acquired in multiple epochs at 2 and 3 months following the recent flaring of its 6.7 GHz methanol maser emission. Using SMA and ALMA, we have discovered 14 new Class II methanol maser lines ranging in frequency from 199 GHz to 361 GHz, which originate mostly from vt=1 torsionally-excited transitions and include one vt=2 transition. The latter detection provides the first observational evidence that Class II maser pumping involves levels in the vt=2 state. The masers are associated with the brightest continuum source (MM1), which hosts a line-rich hot core. The masers present a consistent curvilinear spatial velocity pattern that wraps around MM1, suggestive of a coherent physical structure 1200 au in extent. In contrast, the thermal lines exhibit a linear pattern that crosses MM1 but at progressive position angles that appear to be a function of either increasing temperature or decreasing optical depth. The maser spectral profiles evolved significantly over one month, and the intensities dropped by factors of 3.0 to 7.2, with the vt=2 line showing the largest change. A small area of maser emission from only the highest excitation lines closest to MM1 has disappeared. There are seven additional dust continuum sources in the protocluster, including another hot core (MM3). We do not find evidence for a significant change in (sub)millimeter continuum emission from any of the sources during the one month interval, and the total protocluster emission remains comparable to prior single dish measurements., Comment: 10 pages, 2 figures, 3 tables, accepted for publication in The Astrophysical Journal Letters. All maser transitions detected in this object can be displayed via an option on the advanced interface of http://splatalogue.net

Dynamics of slender magnetic flux tubes (MFT) in the accretion discs of T Tauri stars is investigated. We perform simulations taking into account buoyant, aerodynamic and turbulent drag forces, radiative heat exchange between MFT and ambient gas, magnetic field of the disc. The equations of MFT dynamics are solved using Runge-Kutta method of the fourth order. The simulations show that there are two regimes of MFT motion in absence of external magnetic field. In the region $r<0.2$ au, the MFT of radii $0.05 \leq a_0 \leq 0.16\,H$ ($H$ is the scale height of the disc) with initial plasma beta of 1 experience thermal oscillations above the disc. The oscillations decay over some time, and MFT continue upward motion afterwards. Thinner or thicker MFT do not oscillate. MFT velocity increases with initial radius and magnetic field strength. MFT rise periodically with velocities up to 5-15 km s$^{-1}$ and periods of $0.5-10$ yr determined by the toroidal magnetic field generation time. Approximately 20 % of disc mass and magnetic flux can escape to disc atmosphere via the magnetic buoyancy over characteristic time of disc evolution. MFT dispersal forms expanding magnetized corona of the disc. External magnetic field causes MFT oscillations near the disc surface. These magnetic oscillations have periods from several days to 1-3 months at $r < 0.6$ au. The magnetic oscillations decay over few periods. We simulate MFT dynamics in accretion discs in the Chameleon I cluster. The simulations demonstrate that MFT oscillations can produce observed IR-variability of T Tauri stars., Comment: 18 pages, 12 figures, 3 tables, accepted for publication in MNRAS

H$_2$O maser emission associated with the massive star formation region W49N were observed with the Space-VLBI mission RadioAstron. The procedure for processing of the maser spectral line data obtained in the RadioAstron observations is described. Ultra-fine spatial structures in the maser emission were detected on space-ground baselines of up to 9.6 Earth diameters. The correlated flux densities of these features range from 0.1% to 0.6% of the total flux density. These low values of correlated flux density are probably due to turbulence either in the maser itself or in the interstellar medium., Comment: Accepted for publication in Advances in Space Research

Although the 2P3/2-2P1/2 transition of [CII] at 158um is known to be an excellent tracer of active star formation, we still do not have a complete understanding of where within star formation regions the emission originates. Here, we use SOFIA upGREAT observations of [CII] emission toward the HII region complex Sh2-235 (S235) to better understand in detail the origin of [CII] emission. We complement these data with a fully-sampled Green Bank Telescope radio recombination line map tracing the ionized hydrogen gas. About half of the total [CII] emission associated with S235 is spatially coincident with ionized hydrogen gas, although spectroscopic analysis shows little evidence that this emission is coming from the ionized hydrogen volume. Velocity-integrated [CII] intensity is strongly correlated with WISE 12um intensity across the entire complex, indicating that both trace ultra-violet radiation fields. The 22um and radio continuum intensities are only correlated with [CII] intensity in the ionized hydrogen portion of the S235 region and the correlations between the [CII] and molecular gas tracers are poor across the region. We find similar results for emission averaged over a sample of external galaxies, although the strength of the correlations is weaker. Therefore, although many tracers are correlated with the strength of [CII] emission, only WISE 12um emission is correlated on small-scales of the individual HII region S235 and also has a decent correlation at the scale of entire galaxies. Future studies of a larger sample of Galactic HII regions would help to determine whether these results are truly representative., Comment: Accepted to ApJ

We present the unambiguous discovery of six new class II methanol maser transitions, three of which are torsionally excited (vt=1). The newly discovered 6.18-GHz 17_-2 -> 18_-3 E (vt=1), 7.68-GHz 12_4 -> 13_3 A- (vt=0), 7.83-GHz 12_4 -> 13_3 A+ (vt = 0), 20.9-GHz 10_1 -> 11_2 A+ (vt=1), 44.9-GHz 2_0 -> 3_1 E (vt=1) and 45.8-GHz 9_3 -> 10_2 E (vt=0) methanol masers were detected towards G358.931-0.030, where the known 6.68-GHz maser has recently been reported to be undergoing a period flaring. The detection of the vt=1 torsionally excited lines corroborates one of the missing puzzle pieces in class II maser pumping, but the intensity of the detected emission provides an additional challenge, especially in the case of the very highly excited 6.18-GHz line. Together with the newly detected vt=0 lines, these observations provide significant new information which can be utilised to improve class II methanol maser modelling. We additionally present detections of 6.68-, 19.9-, 23.1- and 37.7-GHz class II masers, as well as 36.2- and 44.1-GHz class I methanol masers, and provide upper limits for the 38.3- and 38.5-GHz class II lines. Near simultaneous Australia Telescope Compact Array (ATCA) observations confirm that all 10 of the class II methanol maser detections are co-spatial to ~0.2 arcsec, which is within the uncertainty of the observations. We find significant levels of linearly polarised emission in the 6.18-, 6.67-, 7.68-, 7.83-, 20.9-, 37.7-, 44.9- and 45.8-GHz transitions, and low levels of circular polarisation in the 6.68-, 37.7- and 45.8-GHz transitions., Comment: ApJL accepted 15 April 2019

The 20 cm radio continuum fluxes of 91 HII regions in a previously compiled catalog have been determined. The spectral types of the ionizing stars in 42 regions with known distances are estimated. These spectral types range from B0.5 to O7, corresponding to effective temperatures of 29 000-37 000 K. The dependences of the infrared (IR) fluxes at 8, 24, and 160 $\mu$m on the 20 cm flux are considered. The IR fluxes are used as a diagnostic of heating of the matter, and the radio fluxes as measurements of the number of ionizing photons. It is established that the IR fluxes grow approximately linearly with the radio flux. This growth of the IR fluxes probably indicates a growth of the mass of heated material in the envelope surrounding the HII region with increasing effective temperature of the star., Comment: 16, pages, 10 figures, published in Astronomy Reports

We report variation of K-band infrared (IR) emission in the vicinity of the G025.65+1.05 water and methanol maser source. New observational data were obtained with 2.5m telescope of the Caucasian Mountain Observatory (CMO) of Moscow State University on 2017-09-21 during the strong water maser flare. We found that the IR source situated close to the maser position had decreased brightness in comparison to archive data. This source is associated with a massive young stellar object (MYSO) corresponding to the compact infrared source IRAS 18316-0602 (RAFGL 7009S). Similar decrease in K-brightness of the IR source close to the maser position was observed in March~2011 when the water maser activity was increased. The dips in MYSO brightness can be related to the maser flare phases. Maser flares that are concurrent with dips of the IR emission can be explained if the lower IR radiation field enables more efficient sink of the pumping cycle by allowing IR photons to escape the maser region., Comment: Accepted for publication in RAA

We compare multi-epoch sub-arcsecond VLA imaging of the 22 GHz water masers toward the massive protocluster NGC6334I observed before and after the recent outburst of MM1B in (sub)millimeter continuum. Since the outburst, the water maser emission toward MM1 has substantially weakened. Simultaneously, the strong water masers associated with the synchrotron continuum point source CM2 have flared by a mean factor of 6.5 (to 4.2 kJy) with highly-blueshifted features (up to 70 km/s from LSR) becoming more prominent. The strongest flaring water masers reside 3000 au north of MM1B and form a remarkable bow shock pattern whose vertex coincides with CM2 and tail points back to MM1B. Excited OH masers trace a secondary bow shock located ~120 au downstream. ALMA images of CS (6-5) reveal a highly-collimated north-south structure encompassing the flaring masers to the north and the non-flaring masers to the south seen in projection toward the MM3-UCHII region. Proper motions of the southern water masers over 5.3 years indicate a bulk projected motion of 117 km/s southward from MM1B with a dynamical time of 170 yr. We conclude that CM2, the water masers, and many of the excited OH masers trace the interaction of the high velocity bipolar outflow from MM1B with ambient molecular gas. The previously-excavated outflow cavity has apparently allowed the radiative energy of the current outburst to propagate freely until terminating at the northern bow shock where it strengthened the masers. Additionally, water masers have been detected toward MM7 for the first time, and a highly-collimated CS (6-5) outflow has been detected toward MM4., Comment: 19 pages, 9 figures, accepted for publication by The Astrophysical Journal

We present the database of maser sources in H2O, OH and SiO lines that can be used to identify and study variable stars at evolved stages. Detecting the maser emission in H2O, OH and SiO molecules toward infrared-excess objects is one of the methods of identification long-period variables (LPVs, including Miras and Semi-Regular), because these stars exhibit maser activity in their circumstellar shells. Our sample contains 1803 known LPV objects. 46% of these stars (832 objects) manifest maser emission in the line of at least one molecule: H2O, OH or SiO. We use the database of circumstellar masers in order to search for long-periodic variables which are not included in the General Catalogue of Variable Stars (GCVS). Our database contains 4806 objects (3866 objects without associations in GCVS catalog) with maser detection in at least one molecule. Therefore it is possible to use the database in order to locate and study the large sample of long-period variable stars. Entry to the database at http://maserdb.net, Comment: Accepted for publication in RAA

We present results of a high resolution study of the filamentary infrared dark cloud G192.76+00.10 in the S254-S258 OB complex in several molecular species tracing different physical conditions. These include three isotopologues of carbon monoxide (CO), ammonia (NH$_3$), carbon monosulfide (CS). The aim of this work is to study the general structure and kinematics of the filamentary cloud, its fragmentation and physical parameters. The gas temperature is derived from the NH$_3 $ $(J,K) = (1,1), (2,2)$ and $^{12}$CO(2--1) lines and the $^{13}$CO(1--0), $^{13}$CO(2--1) emission is used to investigate the overall gas distribution and kinematics. Several dense clumps are identified from the CS(2--1) data. Values of the gas temperature lie in the ranges $10-35$ K, column density $N(\mathrm{H}_2)$ reaches the value 5.1 10$^{22}$ cm$^{-2}$. The width of the filament is of order 1 pc. The masses of the dense clumps range from $ \sim 30 $ M$_\odot$ to $ \sim 160 $ M$_\odot$. They appear to be gravitationally unstable. The molecular emission shows a gas dynamical coherence along the filament. The velocity pattern may indicate longitudinal collapse., Comment: 10 pages, 9 figures, accepted for publication in Research in Astronomy and Astrophysics

We present the first VLBI observations of a Galactic water maser (in Chepeus A) made with a very long baseline interferometric array involving the RadioAstron Earth-orbiting satellite station as one of its elements. We detected two distinct components at -16.9 and 0.6 km/s with a fringe spacing of 66 microarcseconds. In total power, the 0.6 km/s component appears to be a single Gaussian component of strength 580 Jy and width of 0.7 km/s. Single-telescope monitoring showed that its lifetime was only 8~months. The absence of a Zeeman pattern implies the longitudinal magnetic field component is weaker than 120 mG. The space-Earth cross power spectrum shows two unresolved components smaller than 15 microarcseconds, corresponding to a linear scale of 1.6 x 10^11 cm, about the diameter of the Sun, for a distance of 700 pc, separated by 0.54 km/s in velocity and by 160 +/-35 microarcseconds in angle. This is the smallest angular structure ever observed in a Galactic maser. The brightness temperatures are greater than 2 x 10^14K, and the line widths are 0.5 km/s. Most of the flux (about 87%) is contained in a halo of angular size of 400 +/- 150 microarcseconds. This structure is associated with the compact HII region HW3diii. We have probably picked up the most prominent peaks in the angular size range of our interferometer. We discuss three dynamical models: (1) Keplerian motion around a central object, (2) two chance overlapping clouds, and (3) vortices caused by flow around an obstacle (i.e., von Karman vortex street) with Strouhal number of about~0.3., Comment: 15 pages, 9 figures. Accepted for publication in ApJ, February 16, 2018

Observations of the masers in the course of RadioAstron mission yielded detections of fringes for a number of sources in both water and hydroxyl maser transitions. Several sources display numerous ultra-compact details. This proves that implementation of the space VLBI technique for maser studies is possible technically and is not always prevented by the interstellar scattering, maser beaming and other effects related to formation, transfer, and detection of the cosmic maser emission. For the first time, cosmic water maser emission was detected with projected baselines exceeding Earth Diameter. It was detected in a number of star-forming regions in the Galaxy and megamaser galaxies NGC 4258 and NGC 3079. RadioAstron observations provided the absolute record of the angular resolution in astronomy. Fringes from the NGC 4258 megamaser were detected on baseline exceeding 25 Earth Diameters. This means that the angular resolution sufficient to measure the parallax of the water maser source in the nearby galaxy LMC was directly achieved in the cosmic maser observations. Very compact features with angular sizes about 20 microarcsec have been detected in star-forming regions of our Galaxy. Corresponding linear sizes are about 5-10 million kilometers. So, the major step from milli- to micro-arcsecond resolution in maser studies is done in the RadioAstron mission. The existence of the features with extremely small angular sizes is established. Further implementations of the space-VLBI maser instrument for studies of the nature of cosmic objects, studies of the interaction of extremely high radiation field with molecular material and studies of the matter on the line of sight are planned., Comment: To be published in Astrophysical Masers: Unlocking the Mysteries of the Universe, IAU Symposium 336, 2017

We present estimates of brightness temperature for 5 galactic masers in star-forming regions detected at space baselines. Very compact features with angular sizes of about 23-60 micro arcsec were detected in these regions with corresponding linear sizes of about 4-10 million km. Brightness temperatures range from 1e+14 up to 1e+16 K., Comment: To be published in Astrophysical Masers: Unlocking the Mysteries of the Universe, IAU Symposium 336, 2017

We simulate the dynamics of slender magnetic flux tubes (MFTs) in the accretion disks of T Tauri stars. The dynamical equations of our model take into account the aerodynamic and turbulent drag forces, and the radiative heat exchange between the MFT and ambient gas. The structure of the disk is calculated with the help of our MHD model of the accretion disks. We consider the MFTs formed at the distances $0.027-0.8$ au from the star with various initial radii and plasma betas $\beta_0$. The simulations show that the MFT with weak magnetic field ($\beta_0=10$) rise slowly with speeds less than the sound speed. The MFTs with $\beta_0=1$ form outflowing magnetized corona of the disk. Strongly magnetized MFTs ($\beta_0=0.1$) can cause the outflows with velocities $20-50$ km s$^{-1}$. The tubes rise periodically over times from several days to several months according to our simulations. We propose that periodically rising MFTs can absorb stellar radiation and contribute to the IR-variability of Young Stellar Objects., Comment: accepted for publication in RAA

The relationships between atomic and molecular hydrogen and dust of various sizes in extragalactic star-forming regions are considered, based on observational data from the Spitzer and Herschel infrared space telescopes, the Very Large Array (atomic hydrogen emission) and IRAM (CO emission). The source sample consists of approximately 300 star-forming regions in 11 nearby galaxies. Aperture photometry has been applied to measure the fluxes in eight infrared bands (3.6, 4.5, 5.8, 8, 24, 70, 100, and 160$\mu$m), the atomic hydrogen (21cm) line and CO (2--1) lines. The parameters of the dust in the starforming regions were determined via synthetic-spectra fitting, such as the total dust mass, the fraction of polycyclic aromatic hydrocarbons (PAHs), etc. Comparison of the observed fluxes with the measured parameters shows that the relationships between atomic hydrogen, molecular hydrogen, and dust are different in low- and high-metallicity regions. Low-metallicity regions contain more atomic gas, but less molecular gas and dust, including PAHs. The mass of dust constitutes about $1\%$ of the mass of molecular gas in all regions considered. Fluxes produced by atomic and molecular gas do not correlate with the parameters of the stellar radiation, whereas the dust fluxes grow with increasing mean intensity of stellar radiation and the fraction of enhanced stellar radiation. The ratio of the fluxes at 8 and 24$\mu$m, which characterizes the PAH content, decreases with increasing intensity of the stellar radiation, possibly indicating evolutionary variations of the PAH content. The results confirm that the contribution of the 24$\mu$m emission to the total IR luminosity of extragalactic star-forming regions does not depend on the metallicity., Comment: Published in Astronomy Reports, 2017, vol. 61, issue 8

We studied the distribution of dense gas in a filamentary molecular cloud containing several dense clumps. The center of the filament is given by the dense clump WB673. The clumps are high-mass and intermediate-mass star-forming regions. We observed CS(2-1), 13CO(1-0), C18O(1-0) and methanol lines at 96GHz toward WB673 with the Onsala Space Observatory 20-m telescope. We found CS(2-1) emission in the inter-clump medium so the clumps are physically connected and the whole cloud is indeed a filament. Its total mass is $10^4$ M$_{\odot}$ and mass-to-length ratio is 360 M$_{\odot}$pc$^{-1}$ from 13CO(1-0) data. Mass-to-length ratio for the dense gas is $3.4-34$ M$_{\odot}$pc$^{-1}$ from CS(2-1) data. The PV-diagram of the filament is V-shaped. We estimated physical conditions in the molecular gas using methanol lines. Location of the filament on the sky between extended shells suggests that it could be a good example to test theoretical models of formation of the filaments via multiple compression of interstellar gas by supersonic waves.

Aims. We aimed at investigating the structure and kinematics of the gaseous disk and outflows around the massive YSO S255 NIRS3 in the S255IR-SMA1 dense clump. Methods. Observations of the S255IR region were carried out with ALMA at two epochs in the compact and extended configurations. Results. We serendipitously detected a new, never predicted, bright maser line at about 349.1 GHz, which most probably represents the CH$_3$OH $14_{1} - 14_{0}$ A$^{- +}$ transition. The emission covers most of the 6.7 GHz methanol maser emission area of almost 1$^{\prime\prime}$ in size and shows a velocity gradient in the same sense as the disk rotation. No variability was found on the time interval of several months. It is classified as Class II maser and probably originates in a ring at a distance of several hundreds AU from the central star., Comment: 4 pages, 4 figures, accepted by Astronomy and Astrophysics

We present results of a survey of 14 star-forming regions from the Perseus spiral arm in CS(2-1) and 13CO(1-0) lines with the Onsala Space Observatory 20 m telescope. Maps of 10 sources in both lines were obtained. For the remaining sources a map in just one line or a single-point spectrum were obtained. On the basis of newly obtained and published observational data we consider the relation between velocities of the "quasi-thermal" CS(2-1) line and 6.7 GHz methanol maser line in 24 high-mass star-forming regions in the Perseus arm. We show that, surprisingly, velocity ranges of 6.7 GHz methanol maser emission are predominantly red-shifted with respect to corresponding CS(2-1) line velocity ranges in the Perseus arm. We suggest that the predominance of the "red-shifted masers" in the Perseus arm could be related to the alignment of gas flows caused by the large-scale motions in the Galaxy. Large-scale galactic shock related to the spiral structure is supposed to affect the local kinematics of the star-forming regions. Part of the Perseus arm, between galactic longitudes from 85deg to 124deg, does not contain blue-shifted masers at all. Radial velocities of the sources are the greatest in this particular part of the arm, so the velocity difference is clearly pronounced. 13CO(1-0) and CS(2-1) velocity maps of G183.35-0.58 show gas velocity difference between the center and the periphery of the molecular clump up to 1.2 km/s. Similar situation is likely to occur in G85.40-0.00. This can correspond to the case when the large-scale shock wave entrains the outer parts of a molecular clump in motion while the dense central clump is less affected by the shock., Comment: accepted by Astronomy Reports

Aims: In order to test the nature of an (accretion) disk in the vicinity of Cepheus A HW2, we measured the three-dimensional velocity field of the CH3OH maser spots, which are projected within 1000au of the HW2 object, with an accuracy of the order of 0.1km/s. Methods: We made use of the European VLBI Network (EVN) to image the 6.7GHz CH3OH maser emission towards Cepheus A HW2 with 4.5 milli-arcsecond resolution (3au). We observed at three epochs spaced by one year between 2013 and 2015. During the last epoch, on mid-march 2015, we benefited from the new deployed Sardinia Radio Telescope. Results: We show that the CH3OH velocity vectors lie on a preferential plane for the gas motion with only small deviations of 12+/-9 degrees away from the plane. This plane is oriented at a position angle of 134 degrees east of north, and inclined by 26 degrees with the line-of-sight, closely matching the orientation of the disk-like structure previously reported by Patel et al.(2005). Knowing the orientation of the equatorial plane, we can reconstruct a face-on view of the CH3OH gas kinematics onto the plane. CH3OH maser emission is detected within a radius of 900au from HW2, and down to a radius of about 300au, the latter coincident with the extent of the dust emission at 0.9mm. The velocity field is dominated by an infall component of about 2km/s down to a radius of 300au, where a rotational component of 4km/s becomes dominant. We discuss the nature of this velocity field and the implications for the enclosed mass. Conclusions: These findings bring direct support to the interpretation that the high-density gas and dust emission, surrounding Cepheus A HW2, trace an accretion disk., Comment: 9 pages, 4 figures, 2 tables, accepted by Astronomy & Astrophysics

The recent detection of gas-phase methanol (CH$_3$OH) lines in the disc of TW Hya by Walsh et al. provided the first observational constraints on the complex O-bearing organic content in protoplanetary discs. The emission has a ring-like morphology, with a peak at $\sim 30-50$ au and an inferred column density of $\sim 3-6\times10^{12}$ cm$^{-2}$. A low CH$_3$OH fractional abundance of $\sim 0.3-4\times 10^{-11}$ (with respect to H$_2$) is derived, depending on the assumed vertical location of the CH$_3$OH molecular layer. In this study, we use a thermo-chemical model of the TW Hya disc, coupled with the ALCHEMIC gas-grain chemical model, assuming laboratory-motivated, fast diffusivities of the surface molecules to interpret the CH$_3$OH detection. Based on this disc model, we performed radiative transfer calculations with the LIME code and simulations of the observations with the CASA simulator. We found that our model allows to reproduce the observations well. The CH$_3$OH emission in our model appears as a ring with radius of $\sim60$ au. Synthetic and observed line flux densities are equal within the rms noise level of observations. The synthetic CH$_3$OH spectra calculated assuming local thermodynamic equilibrium (LTE) can differ by up to a factor of 3.5 from the non-LTE spectra. For the strongest lines, the differences between LTE and non-LTE flux densities are very small and practically negligible. Variations in the diffusivity of the surface molecules can lead to variations of the CH$_3$OH abundance and, therefore, line flux densities by an order of magnitude., Comment: Accepted for publication in MNRAS, 8 pages, 8 figures