Search

New observations are presented of millimetre line emissions of the circumstellar envelope (CSE) of AGB star $\chi$ Cygni using the recently upgraded NOEMA array. $\chi$ Cygni is an S-type Mira variable, at the border between oxygen-rich and carbon-rich, that has been observed for over 40 years to display features giving evidence for the strong role played by pulsation-associated shock waves in the generation of its wind. The new observations give evidence for a bright H$^{12}$CN(3-2) line emission confined to the very close neighbourhood of the star, however significantly more extended in 2024 than in 2023. Interpretation of such variability in terms of maser emission is considered but found to raise significant unanswered questions. Moreover, other unexpected features are observed in the very close neighbourhood of the star, including low Si$^{16}$O(6-5)/Si$^{17}$O(6-5), $^{28}$SiO(5-4)/$^{29}$SiO(5-4) and $^{12}$CO(2-1)/$^{13}$CO(2-1) line emission ratios. Possibly confirming the important role played by shocks, we discuss a measurement of the SiO(5-4)/SiO(6-5) emission ratio, evidence for a recent mass ejection particularly enhanced in the north-western red-shifted octant of the circumstellar envelope, leaving a depression of emission in its wake, and patterns of enhanced CO(2-1) line-emission suggesting an interpretation in terms of episodic outflows enhanced over solid angles associated with the surface of convective cells and on a time scale of a few decades. Unravelling the mechanisms underlying such newly observed features is very challenging and confirmation of the reported observations with improved sensitivity and angular resolution would be highly welcome. Observation of SiO maser emission in the ($\nu$=1,$J$=6-5) transition is reported for the first time., Comment: accepted for publication in A&A

Deuterated molecules and their molecular D/H-ratios (RD(D)) are important diagnostic tools to study the physical conditions of star-forming regions. The degree of deuteration, RD(D), can be significantly enhanced over the elemental D/H-ratio depending on physical parameters. Within the Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE), we aim to explore the large-scale distribution of deuterated molecules in the nearby Cygnus-X region. We focus on the analysis of large-scale structures of deuterated molecules in the filamentary region hosting the prominent Hii region DR21 and DR21(OH). Here we discuss the HCO+, HNC and HCN molecules and their deuterated isotopologues DCO+, DNC and DCN. The spatial distributions of integrated line emissions from DCO+, DNC, and DCN reveal morphological differences. DCO+ displays the most extended emission, characterized by several prominent peaks. Likewise, DNC exhibits multiple peaks, although its emission appears less extended compared to DCO+. In contrast to the extended emission of DCO+ and DNC, DCN appears the least extended, with distinct peaks. Focusing only on the regions where all three molecules are observed, the mean deuteration ratios for each species are 0.01 for both DNC and DCN, and = 0.005 for DCO+. Anti-correlations are found with deuterated molecules and dust temperature or N(H2). The strongest anti-correlation is found with RD(DCO+) and N(H2). The anti-correlation of RD(DCO+) and N(H2) is suggested to be a result of a combination of an increased photodissociation degree and shocks. A strong positive correlation between the ratio of integrated intensities of DCN and DNC with their 13C-isotopologues, are found in high column density regions. The positive relationship between the ratios implies that the D-isotopologue of the isomers could potentially serve as a tracer for the kinetic gas temperature., Comment: 24 pages, 21 figures, accepted to A&A

GRB~230812B is a bright and relatively nearby ($z =0.36$) long gamma-ray burst (GRB) that has generated significant interest in the community and has thus been observed over the entire electromagnetic spectrum. We report over 80 observations in X-ray, ultraviolet, optical, infrared, and sub-millimeter bands from the GRANDMA (Global Rapid Advanced Network for Multi-messenger Addicts) network of observatories and from observational partners. Adding complementary data from the literature, we then derive essential physical parameters associated with the ejecta and external properties (i.e. the geometry and environment) of the GRB and compare with other analyses of this event. We spectroscopically confirm the presence of an associated supernova, SN2023pel, and we derive a photospheric expansion velocity of v $\sim$ 17$\times10^3$ km s$^{-1}$. We analyze the photometric data first using empirical fits of the flux and then with full Bayesian Inference. We again strongly establish the presence of a supernova in the data, with a maximum (pseudo-)bolometric luminosity of $5.75 \times 10^{42}$ erg/s, at $15.76^{+0.81}_{-1.21}$ days (in the observer frame) after the trigger, with a half-max time width of 22.0 days. We compare these values with those of SN1998bw, SN2006aj, and SN2013dx. Our best-fit model favours a very low density environment ($\log_{10}({n_{\rm ISM}/{\rm cm}^{-3}}) = -2.38^{+1.45}_{-1.60}$) and small values for the jet's core angle $\theta_{\rm core} = 1.54^{+1.02}_{-0.81} \ \rm{deg}$ and viewing angle $\theta_{\rm obs} = 0.76^{+1.29}_{-0.76} \ \rm{deg}$. GRB 230812B is thus one of the best observed afterglows with a distinctive supernova bump.

Molecular outflows are believed to be a key ingredient in the process of star formation. The molecular outflow associated with DR21 Main in Cygnus-X is one of the most extreme, in mass and size, molecular outflows in the Milky Way. The outflow is suggested to belong to a rare class of explosive outflows which are formed by the disintegration of protostellar systems.We aim to explore the morphology, kinematics,and energetics of the DR21 Main outflow, and compare those properties to confirmed explosive outflows to unravel the underlying driving mechanism behind DR21. Line and continuum emission are studied at a wavelength of 3.6\,mm with IRAM 30 m and NOEMA telescopes as part of the Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE) program. The spectra include ($J= 1-0$) transitions of HCO$^+$, HCN, HNC, N$_2$H$^+$, H$_2$CO, CCH tracing different temperature and density regimes of the outflowing gas at high-velocity resolution ($\sim$ 0.8 km s$^{-1}$). The map encompasses the entire DR21 Main outflow and covers all spatial scales down to a resolution of ~3" ($\sim$ 0.02 pc). Integrated intensity maps of the HCO$^+$ emission reveal a strongly collimated bipolar outflow with significant overlap of the blue- and red-shifted emission. The opening angles of both outflow lobes decrease with velocity, from $\sim80$ to 20$^{\circ}$ for the velocity range from 5 to 45 km s$^{-1}$ relative to the source velocity. No evidence is found for the presence of elongated, "filament-like" structures expected in explosive outflows. N$_2$H$^+$ emission near the western outflow lobe reveals the presence of a dense molecular structure which appears to be interacting with the DR21 Main outflow. The overall morphology as well as the detailed kinematics of the DR21 Main outflow is more consistent with that of a typical bipolar outflow instead of an explosive counterpart., Comment: 26 pages, 37 figures. Accepted for publication in A&A

The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars., Comment: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

Context: While star formation on large molecular cloud scales and on small core and disk scales has been investigated intensely over the past decades, the connection of the large-scale interstellar material with the densest small-scale cores has been a largely neglected field. Methods: Using NOEMA and the IRAM 30\,m telescope, we mapped large areas (640\,arcmin$^2$) of the archetypical star formation complex Cygnus X at 3.6\,mm wavelengths in line and continuum emission. Results: The scope and outline of The Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE) is presented. We then focus on the first observed subregion in Cygnus X, namely the DR20 star formation site, which comprises sources in a range of evolutionary stages from cold pristine gas clumps to more evolved ultracompact H{\sc ii} regions. The data covering cloud to cores scales at a linear spatial resolution of $<5000$\,au reveal several kinematic cloud components that are likely part of several large-scale flows onto the central cores. The temperature structure of the region is investigated by means of the HCN/HNC intensity ratio and compared to dust-derived temperatures. We find that the deuterated DCO$^+$ emission is almost exclusively located toward regions at low temperatures below 20\,K. Investigating the slopes of spatial power spectra of dense gas tracer intensity distributions (HCO$^+$, H$^{13}$CO$^+$, and N$_2$H$^+$), we find comparatively flat slopes between $-2.9$ and $-2.6$, consistent with high Mach numbers and/or active star formation in DR20., Comment: 18 pages, 16 figures, accepted to A&A, also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

The latest evolutionary phases of low- and intermediate mass stars are characterized by complex physical processes like turbulence, convection, stellar pulsations, magnetic fields, condensation of solid particles, and the formation of massive outflows that inject freshly produced heavy elements and dust particles into the interstellar medium. We use the Northern Extended Millimeter Array (NOEMA) to obtain spatially and spectrally resolved observations of the semi-regular Asymptotic Giant Branch star RS Cancri to shed light on the morpho-kinematic structure of its inner, wind forming environment by applying detailed 3-D reconstruction modeling and LTE radiative transfer calculations. We detect 32 lines of 13 molecules and isotopologs (CO, SiO, SO, SO$_2$, H$_2$O, HCN, PN), including several transitions from vibrationally excited states. HCN, H$^{13}$CN, millimeter vibrationally excited H$_2$O, SO, $^{34}$SO, SO$_2$, and PN are detected for the first time in RS Cnc. Evidence for rotation is seen in HCN, SO, SO$_2$, and SiO(v=1). From CO and SiO channel maps, we find an inner, equatorial density enhancement, and a bipolar outflow structure with a mass loss rate of $1 \times 10^{-7}M_\odot{\rm yr}^{-1}$ for the equatorial region and of $2 \times 10^{-7}M_\odot{\rm yr}^{-1}$ for the polar outflows. The $^{12}$CO/$^{13}$CO ratio is measured to be $\sim20$ on average, $24\pm2$ in the polar outflows and $19\pm3$ in the equatorial region. We do not find direct evidence of a companion that might explain this kind of kinematic structure, and explore the possibility that a magnetic field might be the cause of it. The innermost molecular gas is influenced by stellar pulsation and possibly by convective cells that leave their imprint on broad wings of certain molecular lines, such as SiO and SO., Comment: 28 pages, 40 figures, accepted for publication in Astronomy & Astrophysics

We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from the 1.3 mm continuum visibility profiles and the radial temperature profiles are estimated from H2CO rotation temperature maps. Molecular column densities are derived with the line fitting tool XCLASS. The physical and chemical properties are combined by applying the physical-chemical model MUSCLE in order to constrain the chemical timescales of a few line-rich cores. The morphology and spatial correlations of the molecular emission are analyzed using the HOG method. The mid-infrared data show that both regions contain a cluster of young stellar objects. Bipolar molecular outflows are observed in the CO 2-1 transition toward the strong mm cores indicating protostellar activity. We find strong molecular emission of SO, SiO, H2CO, and CH3OH in locations which are not associated with the mm cores. These shocked knots can be either associated with the bipolar outflows or, in the case of ISOSS J23053+5953, with a colliding flow that creates a large shocked region between the mm cores. The mean chemical timescale of the cores is lower (20 000 yr) compared to that of the sources of the more evolved CORE sample (60 000 yr). With the HOG method, we find that the spatial emission of species tracing the extended emission and of shock-tracing molecules are well correlated within transitions of these groups., Comment: 34 pages, 17 figures, accepted for publication in A&A

Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources., Comment: 16 pages, 11 figures, Accepted for publication in A&A

Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we observed two young high-mass star-forming regions, ISOSS22478 and ISOSS23053, in the 1.3mm continuum and spectral line emission at a high angular resolution (~0.8''). Results: We resolved 29 cores that are mostly located along filament-like structures. Depending on the temperature assumption, these cores follow a mass-size relation of approximately M~r^2.0, corresponding to constant mean column densities. However, with different temperature assumptions, a steeper mass-size relation up to M~r^3.0, which would be more likely to correspond to constant mean volume densities, cannot be ruled out. The correlation of the core masses with their nearest neighbor separations is consistent with thermal Jeans fragmentation. We found hardly any core separations at the spatial resolution limit, indicating that the data resolve the large-scale fragmentation well. Although the kinematics of the two regions appear very different at first sight - multiple velocity components along filaments in ISOSS22478 versus a steep velocity gradient of more than 50km/s/pc in ISOSS23053 - the findings can be explained within the framework of a dynamical cloud collapse scenario. Conclusions: While our data are consistent with a dynamical cloud collapse scenario and subsequent thermal Jeans fragmentation, the importance of additional environmental properties, such as the magnetization of the gas or external shocks triggering converging gas flows, is nonetheless not as well constrained and would require future investigation., Comment: 16 pages, 16 figures, accepted for Astronomy and Astrophysics, a higher-resolution version can also be found at https://www.mpia.de/homes/beuther/papers.html

We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profiles are inferred from the 1.37 mm continuum visibilities. The column densities of 11 different species are determined by fitting the emission lines with XCLASS. Within the 18 observed regions, we identify 22 individual cores with associated 1.37 mm continuum emission and with a radially decreasing temperature profile. We find an average temperature power-law index of q = 0.4$\pm$0.1 and an average density power-law index of p = 2.0$\pm$0.2 on scales on the order of several 1 000 au. Comparing these results with values of p derived in the literature suggest that the density profiles remain unchanged from clump to core scales. The column densities relative to N(C18O) between pairs of dense gas tracers show tight correlations. We apply the physical-chemical model MUSCLE to the derived column densities of each core and find a mean chemical age of $\sim$60 000 yrs and an age spread of 20 000-100 000 yrs. With this paper we release all data products of the CORE project available at https://www.mpia.de/core. The CORE sample reveals well constrained density and temperature power-law distributions. Furthermore, we characterize a large variety in molecular richness that can be explained by an age spread confirmed by our physical-chemical modeling. The hot molecular cores show the most emission lines, but we also find evolved cores at an evolutionary stage, in which most molecules are destroyed and thus the spectra appear line-poor again., Comment: 29 pages, 14 figures, accepted for publication in A&A

In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively., Comment: 25 pages, 16 figures

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth., Comment: 22 pages, 13 figures, accepted for publication in A&A

The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions., Comment: 23 pages, 16 figures, 5 tables, accepted for publication in A&A

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

Using ALMA observations of $^{12}$CO(2-1), $^{28}$SiO(5-4) and $^{32}$SO$_2$(16$_{6,10}$-17$_{5,13}$) emissions of the circumstellar envelope of AGB star EP Aqr, we describe the morpho-kinematics governing the nascent wind. Main results are: 1) Two narrow polar structures, referred to as jets, launched from less than 25 au away from the star, build up between $\sim$ 20 au and $\sim$ 100 au to a velocity of $\sim$ 20 \kms. They fade away at larger distances and are barely visible in CO data. 2) SO$_2$, SiO and CO emissions explore radial ranges reaching respectively $\sim$30 au, 250 au and 1000 au from the star, preventing the jets to be detected in SO$_2$ data. 3) Close to the star photosphere, rotation (undetected in SiO and CO data) and isotropic radial expansion combine with probable turbulence to produce a broad SO$_2$ line profile ($\sim$ 7.5 \kms\ FWHM). 4) A same axis serves as axis of rotation close to the star, as jet axis and as axi-symmetry axis at large distances. 5) A radial wind builds up at distances up to $\sim$ 300 au from the star, with larger velocity near polar than equatorial latitudes. 6) A sharp depletion of SiO and CO emissions, starting near the star, rapidly broadens to cover the whole blue-western quadrant, introducing important asymmetry in the CO and particularly SiO observations. 7) The $^{12}$C/$^{13}$C abundance ratio is measured as 9$\pm$2. 8) Plausible interpretations are discussed, in particular assuming the presence of a companion., Comment: 18 pages, 16 figures, MNRAS accepted

Red supergiant stars are surrounded by a gaseous and dusty circumstellar environment created by their mass loss which spreads heavy elements into the interstellar medium. The structure and the dynamics of this envelope are crucial to understand the processes driving the red supergiant mass loss and the shaping of the pre-supernova ejecta. We have observed the emission from the CO $J = 2-1$ line from the red supergiant star $\mu$~Cep with the NOEMA interferometer. In the line the synthesized beam was $0.92 \times 0.72$~arcsec ($590 \times 462$~au at 641~pc). The continuum map shows only the unresolved contribution of the free-free emission of the star chromosphere. The continuum-subtracted channel maps reveal a very inhomogeneous and clumpy circumstellar environment. In particular, we detected a bright CO clump, as bright as the central source in the line, at 1.80~arcsec south-west from the star, in the blue channel maps. After a deprojection of the radial velocity assuming two different constant wind velocities, the observations were modelled using the 3D radiative transfer code \textsc{lime} to derive the characteristics of the different structures. We determine that the gaseous clumps observed around $\mu$~Cep are responsible for a mass loss rate of $(4.9 \pm 1.0) \times 10^{-7}~{\rm M}_\odot\,{\rm yr}^{-1}$, in addition to a spatially unresolved wind component with an estimated mass-loss rate of $2.0 \times 10^{-6}~{\rm M}_\odot\,{\rm yr}^{-1}$. Therefore, the clumps have a significant role in $\mu$~Cep's mass loss ($\ge 25 \%$). We cannot exclude that the unresolved central outflow may be made of smaller unresolved clumps., Comment: 15 pages, 4 tables, 9 figures. 2nd version : one co-author removed and acknowledgement updated (consistent with erratum https://doi.org/10.1093/mnras/stz1006)

We investigate the prompt emission and the afterglow properties of short duration gamma-ray burst (sGRB) 130603B and another eight sGRB events during 2012-2015, observed by several multi-wavelength facilities including the GTC 10.4m telescope. Prompt emission high energy data of the events were obtained by INTEGRAL/SPI/ACS, Swift/BAT and Fermi/GBM satellites. The prompt emission data by INTEGRAL in the energy range of 0.1-10 MeV for sGRB 130603B, sGRB 140606A, sGRB 140930B, sGRB 141212A and sGRB 151228A do not show any signature of the extended emission or precursor activity and their spectral and temporal properties are similar to those seen in case of other short bursts. For sGRB130603B, our new afterglow photometric data constraints the pre jet-break temporal decay due to denser temporal coverage. For sGRB 130603B, the afterglow light curve, containing both our new as well as previously published photometric data is broadly consistent with the ISM afterglow model. Modeling of the host galaxies of sGRB 130603B and sGRB 141212A using the LePHARE software supports a scenario in which the environment of the burst is undergoing moderate star formation activity. From the inclusion of our late-time data for 8 other sGRBs we are able to; place tight constraints on the non-detection of the afterglow, host galaxy or any underlying kilonova emission. Our late-time afterglow observations of the sGRB 170817A/GW170817 are also discussed and compared with the sub-set of sGRBs., Comment: 24 pages, 22 figures, Accepted to MNRAS, 2019 February 19. Received 2019 February 19; in original form 2018 August 30

ALMA observations of CO(1-0) and CO(2-1) emissions of the circumstellar envelope of EP Aqr, an oxygen-rich AGB star, are reported. A thorough analysis of their properties is presented using an original method based on the separation of the data-cube into a low velocity component associated with an equatorial outflow and a faster component associated with a bipolar outflow. A number of important and new results are obtained concerning the distribution in space of the effective emissivity, the temperature, the density and the flux of matter. A mass loss rate of (1.6$\pm$0.4)10$^{-7}$ solar masses per year is measured. The main parameters defining the morphology and kinematics of the envelope are evaluated and uncertainties inherent to de-projection are critically discussed. Detailed properties of the equatorial region of the envelope are presented including a measurement of the line width and a precise description of the observed inhomogeneity of both morphology and kinematics. In particular, in addition to the presence of a previously observed spiral enhancement of the morphology at very small Doppler velocities, a similarly significant but uncorrelated circular enhancement of the expansion velocity is revealed, both close to the limit of sensitivity. The results of the analysis place significant constraints on the parameters of models proposing descriptions of the mass loss mechanism, but cannot choose among them with confidence., Comment: 26 pages, 31 figures, accepted for publication in MNRAS

Atacama Large Millimetre/sub-millimetre Array (ALMA) observations of the CO(1-0) and CO(2-1) emissions from the circumstellar envelope of the Asymptotic Giant Branch (AGB) star EP Aqr have been made with four times better spatial resolution than previously available. They are analysed with emphasis on the de-projection in space of the effective emissivity and flux of matter using as input a prescribed configuration of the velocity field, assumed to be radial. The data are found to display an intrinsic axi-symmetry with respect to an axis making a small angle with the line of sight. A broad range of wind configurations, from prolate (bipolar) to oblate (equatorial) has been studied and found to be accompanied by significant equatorial emission. Qualitatively, the effective emissivity is enhanced near the equator to produce the central narrow component observed in the Doppler velocity spectra and its dependence on star latitude generally follows that of the wind velocity with the exception of an omni-present depression near the poles. In particular, large equatorial expansion velocities produce a flared disc or a ring of effective emissivity and mass loss. The effect on the determination of the orientation of the star axis of radial velocity gradients and possibly competing rotation and expansion in the equatorial disc is discussed. In general, the flux of matter is found to reach a broad maximum at distances of the order of 500 au from the star. Arguments are given that may be used to prefer one wind velocity distribution to another. As a result of the improved quality of the data, a deeper understanding of the constraints imposed on morphology and kinematics has been obtained., Comment: 17 pages, 11 figures, accepted for publication in RAA

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales., Comment: 23 pages, 26 figures, accepted for publication in Astronomy and Astrophysics

The problem of de-projection of radio line observations of axi-symmetric expanding circumstellar envelopes is studied with the aim of easing their analysis in terms of physics models. The arguments developed rest on the remark that, in principle, when the wind velocity distribution is known, the effective emissivity can be calculated at any point in space. The paper provides a detailed study of how much this is true in practice. The wind velocity distribution assumed to be axi-symmetric and in expansion, is described by four parameters: the angles defining the orientation of the symmetry axis, an overall velocity scale and a parameter measuring the elongation (prolateness) of the distribution. Tools are developed that allow for measuring, or at least constraining, each of the four parameters. The use of effective emissivity as relevant quantity, rather than temperature and density being considered separately, implies important assumptions and simplifications meaning that the approach being considered here is only a preliminary to, and by no means a replacement for, a physics analysis accounting for radiative transfer and hydrodynamics arguments. While most considerations are developed using simulated observations as examples, two case studies (EP Aqr, observed with ALMA, and RS Cnc, with NOEMA) are presented that illustrate their usefulness in practical cases., Comment: 16 pages, 16 figures, accepted for publication in MNRAS

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation., Comment: 29 pages, 17 figures, Astronomy & Astrophysics in press, for a higher-resolution version of the paper see http://www.mpia.de/homes/beuther/papers.html

Quasars are the most luminous non-transient objects known and as a result they enable studies of the Universe at the earliest cosmic epochs. Despite extensive efforts, however, the quasar ULAS J1120+0641 at z=7.09 has remained the only one known at z>7 for more than half a decade. Here we report observations of the quasar ULAS J134208.10+092838.61 (hereafter J1342+0928) at redshift z=7.54. This quasar has a bolometric luminosity of 4e13 times the luminosity of the Sun and a black hole mass of 8e8 solar masses. The existence of this supermassive black hole when the Universe was only 690 million years old---just five percent of its current age---reinforces models of early black-hole growth that allow black holes with initial masses of more than about 1e4 solar masses or episodic hyper-Eddington accretion. We see strong evidence of absorption of the spectrum of the quasar redwards of the Lyman alpha emission line (the Gunn-Peterson damping wing), as would be expected if a significant amount (more than 10 per cent) of the hydrogen in the intergalactic medium surrounding J1342+0928 is neutral. We derive a significant fraction of neutral hydrogen, although the exact fraction depends on the modelling. However, even in our most conservative analysis we find a fraction of more than 0.33 (0.11) at 68 per cent (95 per cent) probability, indicating that we are probing well within the reionization epoch of the Universe., Comment: Updated to match the final journal version

The formation of silicate dust in oxygen-rich envelopes of evolved stars is thought to be initiated by formation of seed particles that can withstand the high temperatures close to the stellar photosphere and act as condensation cores farther away from the star. Candidate species considered as first condensates are TiO and TiO$_2$. We aim to identify and characterize the circumstellar gas-phase chemistry of titanium that leads to the formation of solid titanium compounds in the envelope of $o$ Cet, the prototypical Mira, and seek an observational verification of whether titanium oxides play a major role in the onset of dust formation in M-type AGB stars. We present high angular-resolution ALMA observations at submillimeter (submm) wavelengths supplemented by APEX and Herschel spectra of the rotational features of TiO and TiO$_2$. In addition, circumstellar features of TiO and TiI are identified in optical spectra which cover multiple pulsation cycles of $o$ Cet. The submm ALMA data reveal TiO and TiO$_2$ bearing gas within the extended atmosphere of Mira. While TiO is traceable up to a radius (FWHM/2) of 4.0 R$_{\star}$, TiO$_2$ extends as far as 5.5 R$_{\star}$ and unlike TiO appears to be anisotropically distributed. Optical spectra display variable emission of TiI and TiO from inner parts of the extended atmosphere (<3 R$_{\star}$). Chemical models which include shocks are in general agreement with the observations of gas-phase titanium-bearing molecules. It is unlikely that substantial amounts of titanium is locked up in solids because the abundance of the gaseous titanium species is very high. In particular, formation of hot titanium-rich condensates is very improbable because we find no traces of their hot precursor species in the gas phase. It therefore appears unlikely that the formation of dust in Mira, and possibly other M-type AGB stars, is initiated by titanium oxides., Comment: 19 pages plus spectral atlas of the TiO gamma' band in 1998-2014 submitted to A&A

We report on observations obtained at IRAM on two semi-regular variable Asymptotic Giant Branch (AGB) stars, RS Cnc and EP Aqr, undergoing mass loss at an intermediate rate of ~ 10^-7 solar mass per year. Interferometric data obtained with the Plateau-de-Bure interferometer (NOEMA) have been combined with On-The-Fly maps obtained with the 30-m telescope in the CO(1-0) and (2-1) rotational lines. The spectral maps of spatially resolved sources reveal an axisymmetric morphology in which matter is flowing out at a low velocity (~ 2 km/s) in the equatorial planes, and at a larger velocity (~ 8 km/s) along the polar axes. There are indications that this kind of morpho-kinematics is relatively frequent among stars at the beginning of their evolution on the Thermally-Pulsing AGB, in particular among those that show composite CO line profiles, and that it might be caused by the presence of a companion. We discuss the progress that could be expected for our understanding of the mass loss mechanisms in this kind of sources by increasing the spatial resolution of the observations with ALMA or NOEMA., Comment: Proceedings 2016 SF2A days

Context: Dust is efficiently produced by cool giant stars, but the condensation of inorganic dust is poorly understood. Aims: Identify and characterize aluminum bearing species in the circumstellar gas of Mira ($o$ Ceti) in order to elucidate their role in the production of Al$_2$O$_3$ dust. Methods: Multiepoch spectral line observations at (sub-)millimeter, far-infrared, and optical wavelengths including: maps with ALMA which probe the gas distribution in the immediate vicinity of the star at ~30 mas; observations with ALMA, APEX, and Herschel in 2013-2015 for studying cycle and inter-cycle variability of the rotational lines of Al bearing molecules; optical records as far back as 1965 to examine variations in electronic transitions over time spans of days to decades; and velocity measurements and excitation analysis of the spectral features which constrain the physical parameters of the gas. Results: Three diatomic molecules AlO, AlOH, and AlH, and atomic Al I are the main observable aluminum species in Mira, although a significant fraction of aluminum might reside in other species that have not yet been identified. Strong irregular variability in the (sub-)millimeter and optical features of AlO (possibly the direct precursor of Al$_2$O$_3$) indicates substantial changes in the excitation conditions, or varying abundance that is likely related to shocks in the star. The inhomogeneous distribution of AlO might influence the spatial and temporal characteristics of dust production. Conclusions: We are unable to quantitatively trace aluminum depletion from the gas, but the rich observational material constrains time dependent chemical networks. Future improvements should include spectroscopic characterization of higher aluminum oxides, coordinated observations of dust and gas species at different variability phases, and tools to derive abundances in shock excited gas., Comment: submitted to A&A, 23 pages plus an atlas of spectral variability of Mira in Al-bearing species since 1965 (90 pages total)

Studies of the CO and HI radio emission of some evolved stars are presented using data collected by the IRAM Plateau de Bure interferometer and Pico Veleta telescope, the Nan\c{c}ay Radio Telescope and the JVLA and ALMA arrays. Approximate axial symmetry of the physical and kinematic properties of the circumstellar envelope (CSE) are observed in CO emission, in particular, from RS Cnc, EP Aqr and the Red Rectangle. A common feature is the presence of a bipolar outflow causing an enhanced wind velocity in the polar directions. HI emission extends to larger radial distances than probed by CO emission and displays features related to the interaction between the stellar outflow and interstellar matter. With its unprecedented sensitivity, FAST will open a new window on such studies. Its potential in this domain is briefly illustrated., Comment: 6 pages, 4 figures, conference (Frontiers in Radio Astronomy and FAST Early Sciences Symposium 2015)

The CO(1-0) and (2-1) emission of the circumstellar envelope of the AGB star EP Aqr has been observed using the IRAM PdBI and the IRAM 30-m telescope. The line profiles reveal the presence of two distinct components centered on the star velocity, a broad component extending up to ~10 km/s and a narrow component indicating an expansion velocity of ~2 km/s. An early analysis of these data was performed under the assumption of isotropic winds. The present study revisits this interpretation by assuming instead a bipolar outflow nearly aligned with the line of sight. A satisfactory description of the observed flux densities is obtained with a radial expansion velocity increasing from ~2 km/s at the equator to ~10 km/s near the poles. The angular aperture of the bipolar outflow is ~45 deg with respect to the star axis, which makes an angle of ~13 deg with the line of sight. A detailed study of the CO(1-0) to CO(2-1) flux ratio reveals a significant dependence of the temperature on the star latitude, smaller and steeper at the poles than at the equator at large distances from the star. Under the hypothesis of radial expansion and of rotation invariance about the star axis, the effective density has been evaluated in space as a function of star coordinates. Evidence is found for an enhancement of the effective density in the northern hemisphere of the star at angular distances in excess of ~3" and covering the whole longitudinal range. The peak velocity of the narrow component is observed to vary slightly with position on the sky, a variation consistent with the model and understood as the effect of the inclination of the star axis with respect to the line of sight. While the phenomenological model presented here reproduces well the general features of the observations, significant differences are also revealed, which would require a better spatial resolution to be properly described., Comment: accepted for publication in Astronomy & Astrophysics

We obtained high resolution (0."25 to 0."90) observations of HCN and HCO+ J=3-2 of the ultraluminous QSO galaxy Mrk231 with the IRAM Plateau de Bure Interferometer. We find luminous HCN and HCO+ 3-2 emission in the main disk and we detect compact (r<90 pc) vibrationally excited HCN 3-2, v2=1f emission centred on the nucleus. The velocity field of the vibrationally excited HCN is strongly inclined (PA=155 deg.) compared to the east-west rotation of the main disk. The nuclear molecular mass is estimated to 8e8 Msun with an average N(H2)of 1.2e24 cm-2. Prominent, spatially extended (>350 pc) line wings are found for HCN 3-2 with velocities +-750 km/s. Line ratios indicate that the emission is emerging in dense gas n=1e4 - 5e5 cm-3 of elevated HCN abundance X(HCN)=1e-8 to 1e-6. High X(HCN) also allows for the emission to originate in gas of more moderate density. We tentatively detect nuclear emission from the reactive ion HOC+ with HCO+/HOC+=10-20. The HCN v2=1f line emission is consistent with the notion of a hot, dusty, warped inner disk of Mrk231 where the v2=1f line is excited by bright mid-IR 14 micron continuum. We estimate the vibrational temperature T_vib to 200-400 K. We propose that 50% of the main HCN emission may have its excitation affected by the radiation field through IR pumping of the vibrational ground state. The HCN emission in the line wings, however, is more extended and thus likely not strongly affected by IR pumping. Our results reveal that dense clouds survive (and/or are formed) in the AGN outflow on scales of at least several hundred pc before evaporating or collapsing. The elevated HCN abundance in the outflow is consistent with warm chemistry possibly related to shocks and/or X-ray irradiated gas. An upper limit to the mass and momentum flux is 4e8 Msun and 12L_AGN/c, respectively, and we discuss possible driving mechanisms for the dense outflow., Comment: Accepted for publication in Astronomy and Astrophysics, 12 pages, 3 tables 8 figures

We study the dark nature of GRB 130528A through multi-wavelength observations and conclude that the main reason for the optical darkness is local extinction inside of the host galaxy. Automatic observations were performed at BOOTES-4/MET robotic telescope. We also triggered target of opportunity (ToO) observation at the OSN, IRAM PdBI and the GTC+OSIRIS. The host galaxy photometric observations in optical to near-infrared (nIR) wavelengths were achieved through large ground-based aperture telescopes, such as the 10.4m GTC, the 4.2m WHT, 6m BTA, and the 2m LT. Based on these observations, spectral energy distributions (SED) for the host galaxy and afterglow were constructed. Thanks to mm observations at PdBI, we confirm the presence of a mm source within the XRT error circle that faded over the course of our observations and identify the host galaxy. However, we do not find any credible optical source within early observations with BOOTES-4/MET and 1.5m OSN telescopes. Spectroscopic observation of this galaxy by GTC showed a single faint emission line that likely corresponds to [OII] 3727\{AA} at a redshift of 1.250+/-0.001 implying a SFR(M_sun/yr) > 6.18 M_sun/yr without correcting for dust extinction. The probable extinction was revealed through analysis of the afterglow SED, resulting in a value of AV >= ~ 0.9 at the rest frame, this is comparable to extinction levels found among other dark GRBs. The SED of the host galaxy is explained well (chi2/d.o.f.=0.564) by a luminous (MB=-21.16), low-extinction (AV =0, rest frame), and aged (2.6 Gyr) stellar population. We can explain this apparent contradiction in global and line-of-sight extinction if the GRB birth place happened to lie in a local dense environment. In light of having relatively small specific SFR (SSFR) ~ 5.3 M_sun/yr (L/L_star)-1, this also could explain the age of the old stellar population of host galaxy., Comment: 6 pages, 4 figures

We present a detailed study of the circumstellar gas distribution and kinematics of the semi-regular variable star RS Cnc on spatial scales ranging from ~1" (~150 AU) to ~6' (~0.25 pc). New modeling of CO1-0 and CO2-1 imaging observations leads to a revised characterization of RS Cnc's previously identified axisymmetric molecular outflow. Rather than a simple disk-outflow picture, we find that a gradient in velocity as a function of latitude is needed to fit the spatially resolved spectra, and in our preferred model, the density and the velocity vary smoothly from the equatorial plane to the polar axis. In terms of density, the source appears quasi-spherical, whereas in terms of velocity the source is axi-symmetric with a low expansion velocity in the equatorial plane and faster outflows in the polar directions. The flux of matter is also larger in the polar directions than in the equatorial plane. An implication of our model is that the stellar wind is still accelerated at radii larger than a few hundred AU, well beyond the radius where the terminal velocity is thought to be reached in an asymptotic giant branch star. The HI data show the previously detected head-tail morphology, but also supply additional detail about the atomic gas distribution and kinematics. We confirm that the `head' seen in HI is elongated in a direction consistent with the polar axis of the molecular outflow, suggesting that we are tracing an extension of the molecular outflow well beyond the molecular dissociation radius (up to ~0.05 pc). The 6'-long HI `tail' is oriented at a PA of 305{\deg}, consistent with the space motion of the star. We measure a total mass of atomic hydrogen ~0.0055 solar mass and estimate a lower limit to the timescale for the formation of the tail to be ~6.4x10^4 years. (abridged), Comment: accepted for publication in Astronomy & Astrophysics

At low redshift, a handful of gamma-ray bursts (GRBs) have been discovered with peak luminosities ($L_{\rm iso} < 10^{48.5}~\rm{erg\,s}^{-1}$) substantially lower than the average of the more distant ones ($L_{\rm iso} > 10^{49.5}~\rm{erg\,s}^{-1}$). The properties of several low-luminosity (low-$L$) GRBs indicate that they can be due to shock break-out, as opposed to the emission from ultrarelativistic jets. Owing to this, it is highly debated how both populations are connected, and whether there is a continuum between them. The burst at redshift $z=0.283$ from 2012 April 22 is one of the very few examples of intermediate-$L$ GRBs with a $\gamma$-ray luminosity of $L\sim10^{48.9}~\rm{erg\,s}^{-1}$ that have been detected up to now. Together with the robust detection of its accompanying supernova SN 2012bz, it has the potential to answer important questions on the origin of low- and high-$L$ GRBs and the GRB-SN connection. We carried out a spectroscopy campaign using medium- and low-resolution spectrographs at 6--10-m class telescopes, covering the time span of 37.3 days, and a multi-wavelength imaging campaign from radio to X-ray energies over a duration of $\sim270$ days. Furthermore, we used a tuneable filter centred at H$\alpha$ to map star formation in the host galaxy and the surrounding galaxies. We used these data to extract and model the properties of different radiation components and incorporate spectral-energy-distribution fitting techniques to extract the properties of the host galaxy. Modelling the light curve and spectral energy distribution from the radio to the X-rays revealed the blast-wave to expand with an initial Lorentz factor of $\Gamma_0\sim60$, low for a high-$L$ GRB, and that the afterglow had an exceptional low peak luminosity-density of $\lesssim2\times10^{30}~\rm{erg\,s}^{-1}\,\rm{Hz}^{-1}$ in the sub-mm. [Abridged], Comment: 30 pages, 17 figures, 9 tables; abstract is abridged; images are shown at reduced resolution; comments are welcome

Events such as GRB130606A at z=5.91, offer an exciting new window into pre-galactic metal enrichment in these very high redshift host galaxies. We study the environment and host galaxy of GRB 130606A, a high-z event, in the context of a high redshift population of GRBs. We have obtained multiwavelength observations from radio to gamma-ray, concentrating particularly on the X-ray evolution as well as the optical photometric and spectroscopic data analysis. With an initial Lorentz bulk factor in the range \Gamma_0 ~ 65-220, the X-ray afterglow evolution can be explained by a time-dependent photoionization of the local circumburst medium, within a compact and dense environment. The host galaxy is a sub-DLA (log N (HI) = 19.85+/-0.15), with a metallicity content in the range from ~1/7 to ~1/60 of solar. Highly ionized species (N V and Si IV) are also detected. This is the second highest redshift burst with a measured GRB-DLA metallicity and only the third GRB absorber with sub-DLA HI column density. GRB ' lighthouses' therefore offer enormous potential as backlighting sources to probe the ionization and metal enrichment state of the IGM at very high redshifts for the chemical signature of the first generation of stars., Comment: 10 pages, 9 figures, submitted to A&A. Typos corrected

We present multiwavelength observations of the afterglow of GRB 130427A, the brightest (in total fluence) gamma-ray burst of the past 29 years. Optical spectroscopy from Gemini-North reveals the redshift of the GRB to be z=0.340, indicating that its unprecedented brightness is primarily the result of its relatively close proximity to Earth; the intrinsic luminosities of both the GRB and its afterglow are not extreme in comparison to other bright GRBs. We present a large suite of multiwavelength observations spanning from 300 s to 130 d after the burst and demonstrate that the afterglow shows relatively simple, smooth evolution at all frequencies with no significant late-time flaring or rebrightening activity. The entire dataset from 1 GHz to 10 GeV can be modeled as synchrotron emission from a combination of reverse and forward shocks in good agreement with the standard afterglow model, providing strong support to the applicability of the underlying theory and clarifying the nature of the GeV emission observed to last for minutes to hours following other very bright GRBs. A tenuous, wind-stratified circumburst density profile is required by the observations, suggesting a massive-star progenitor with a low mass-loss rate, perhaps due to low metallicity. GRBs similar in nature to GRB 130427A, inhabiting low-density media and exhibiting strong reverse shocks, are probably not uncommon but may have been difficult to recognize in the past due to their relatively faint late-time radio emission; more such events should be found in abundance by the new generation of sensitive radio and millimeter instruments., Comment: Accepted to ApJ. Light curves and SEDs are available at http://www.astro.caltech.edu/~dperley/grb/130427a/data/

The formation of inorganic dust in circumstellar environments of evolved stars is poorly understood. Spectra of molecules thought to be most important for the nucleation, i.e. AlO, TiO, and TiO2, have been recently detected in the red supergiant VY CMa. These molecules are effectively formed in VY CMa and the observations suggest that non-equilibrium chemistry must be involved in their formation and nucleation into dust. In addition to exploring the recent observations of VY CMa, we briefly discuss the possibility of detecting these molecules in the dust-poor circumstellar environment of Betelgeuse., Comment: contribution to Betelgeuse Workshop 2012: "The physics of Red Supergiants: recent advances and open questions", 26-29 Nov 2012 Paris (France)

We report the first detection of pure rotational transitions of TiO and TiO_2 at (sub-)millimeter wavelengths towards the red supergiant VY CMa. A rotational temperature, T_rot, of about 250 K was derived for TiO_2. Although T_rot was not well constrained for TiO, it is likely somewhat higher than that of TiO_2. The detection of the Ti oxides confirms that they are formed in the circumstellar envelopes of cool oxygen-rich stars and may be the "seeds" of inorganic-dust formation, but alternative explanations for our observation of TiO and TiO_2 in the cooler regions of the envelope cannot be ruled out at this time. The observations suggest that a significant fraction of the oxides is not converted to dust, but instead remains in the gas phase throughout the outflow., Comment: to appear in Astronomy and Astrophysics

We report a search for linear polarization in the active galactic nucleus (AGN) 3C 84 (NGC 1275) at observed frequencies of 239 GHz and 348 GHz, corresponding to rest-frame frequencies of 243 GHz and 354 GHz. We collected polarization data with the IRAM Plateau de Bure Interferometer via Earth rotation polarimetry. We do not detect linear polarization. Our analysis finds 3-sigma upper limits on the degree of polarization of 0.5% and 1.9% at 239 GHz and 348 GHz, respectively. We regard the influence of Faraday conversion as marginal, leading to expected circular polarizations <0.3%. Assuming depolarization by a local Faraday screen, we constrain the rotation measure, as well as the fluctuations therein, to be 10^6 rad/m^2. From this we estimate line-of-sight magnetic field strengths of >100 microG. Given the physical dimensions of 3C 84 and its observed structure, the Faraday screen appears to show prominent small-scale structure, with \DeltaRM > 10^6 rad/m^2 on projected spatial scales <1 pc., Comment: 7 pages, 4 figures. Accepted by MNRAS

We present an analysis of the linear polarization of six active galactic nuclei - 0415+379 (3C~111), 0507+179, 0528+134 (OG+134), 0954+658, 1418+546 (OQ+530), and 1637+574 (OS+562). Our targets were monitored from 2007 to 2011 in the observatory-frame frequency range 80-253 GHz, corresponding to a rest-frame frequency range 88-705 GHz. We find average degrees of polarization m_L ~ 2-7%; this indicates that the polarization signals are effectively averaged out by the emitter geometries. We see indication for fairly strong shocks and/or complex, variable emission region geometries in our sources, with compression factors <0.9 and/or changes in viewing angles by >10 deg. An analysis of correlations between source fluxes and polarization parameter points out special cases: the presence of (at least) two distinct emission regions with different levels of polarization (for 0415+379) as well as emission from a single, predominant component (for 0507+179 and 1418+546). Regarding the evolution of flux and polarization, we find good agreement between observations and the signal predicted by "oblique shock in jet" scenarios in one source (1418+546). We attempt to derive rotation measures for all sources, leading to actual measurements for two AGN and upper limits for three sources. We derive values of RM = -39,000 +/- 1,000 (stat) +/- 13,000 (sys) rad/m^2 and RM = 420,000 +/- 10,000 (stat) +/- 110,000 (sys) rad/m^2 for 1418+546 and 1637+574, respectively; these are the highest values reported to date for AGN. These values indicate magnetic field strengths of the order ~0.0001 G. For 0415+379, 0507+179, and 0954+658 we derive upper limits |RM| < 17,000 rad/m^2. From the relation |RM| ~ nu^a we find a = 1.9 +/- 0.3 for 1418+546, in good agreement with a = 2 as expected for a spherical or conical outflow., Comment: 23 pages, 8 figures, 4 tables. Accepted by Astronomy and Astrophysics. Minor language editing, one missing reference (Macquart et al. 2006) added

We detect luminous emission from HCN, HCO+ and HNC 1--0 in the QSO ULIRG Mrk~231 with the IRAM Plateau de Bure Interferometer at 1."55 by 1."28 resolution. All three lines show broad line wings - which are particularly prominent for HCN. Velocities are found to be similar (750 km/s) to those found for CO 1-0. This is the first time bright HCN, HCO+ and HNC emission has been detected in a large-scale galactic outflow. We find that both the blue- and red-shifted line wings are spatially extended by at least 0."75 (700 pc) in a north-south direction. The line wings are brighter (relative to the line center intensity) in HCN than in CO 1-0 and line ratios suggest that the molecular outflow consists of dense (n>10E4 cmE-3) and clumpy gas with a high HCN abundance X(HCN)>10E-8. These properties are consistent with the molecular gas being compressed and fragmented by shocks in the outflow. Alternatively, HCN is instead pumped by mid-IR continuum, but we propose that this effect is not strong for the spatially extended outflowing gas. In addition, we find that the rotation of the main disk, in east-west direction, is also evident in the HCN, HCO+ and HNC line emission. An unexpectedly bright HC3N 10-9 line is detected inside the central 400 pc of Mrk231. This HC3N emission may emerge from a shielded, dust-enshrouded region within the inner 40-50 pc where the gas is heated to high temperatures (200 - 300 K) by the AGN., Comment: 8 pages, 7 figures, accepted for publication in Astronomy and Astrophysics

GRBs generate an afterglow emission that can be detected from radio to X-rays during days, or even weeks after the initial explosion. The peak of this emission crosses the mm/submm range during the first hours to days, making their study in this range crucial for constraining the models. Observations have been limited until now due to the low sensitivity of the observatories in this range. We present observations of 10 GRB afterglows obtained from APEX and SMA, as well as the first detection of a GRB with ALMA, and put them into context with all the observations that have been published until now in the spectral range that will be covered by ALMA. The catalogue of mm/submm observations collected here is the largest to date and is composed of 102 GRBs, of which 88 had afterglow observations, whereas the rest are host galaxy searches. With our programmes, we contributed with data of 11 GRBs and the discovery of 2 submm counterparts. In total, the full sample, including data from the literature, has 22 afterglow detections with redshift ranging from 0.168 to 8.2. GRBs have been detected in mm/submm wavelengths with peak luminosities spanning 2.5 orders of magnitude, the most luminous reaching 10^33erg s^-1 Hz^-1. We observe a correlation between the X-ray brightness at 0.5 days and the mm/submm peak brightness. Finally we give a rough estimate of the distribution of peak flux densities of GRB afterglows, based on the current mm/submm sample. Observations in the mm/submm bands have been shown to be crucial for our understanding of the physics of GRBs, but have until now been limited by the sensitivity of the observatories. With the start of the operations at ALMA, the sensitivity will be increased by more than an order of magnitude. Our estimates predict that, once completed, ALMA will detect up to 98% of the afterglows if observed during the passage of the peak synchrotron emission., Comment: 23 pages, 14 figures, 5 tables (one big one!), Accepted for publication in A&A. Includes the first observation of a GRB afterglow with ALMA

We analyze the long-term evolution of the fluxes of six active galactic nuclei (AGN) - 0923+392, 3C 111, 3C 273, 3C 345, 3C 454.3, and 3C 84 - in the frequency range 80-267 GHz using archival calibration data of the IRAM Plateau de Bure Interferometer. Our dataset spans a long timeline of ~14 years with 974-3027 flux measurements per source. We find strong (factors ~2-8) flux variability on timescales of years for all sources. The flux density distributions of five out of six sources show clear signatures of bi- or even multimodality. Our sources show mostly steep (alpha~0.5-1), variable spectral indices that indicate outflow dominated emission; the variability is most probably due to optical depth variations. The power spectra globally correspond to red-noise spectra with five sources being located between the cases of white and flicker noise and one source (3C 111) being closer to the case of random walk noise. For three sources the low-frequency ends of their power spectra appear to be upscaled in spectral power by factors ~2-3 with respect to the overall powerlaws. In two sources, 3C 454.3 and 3C 84, the 1.3-mm emission preceeds the 3-mm emission by ~55 and ~300 days, respectively, probably due to (combinations of) optical depth and emission region geometry effects. We conclude that the source emission cannot be described by uniform stochastic emission processes; instead, a distinction of "quiescent" and (maybe multiple) "flare" states of the source emission appears to be necessary., Comment: 20 pages, 11 figures, 2 tables. Accepted by Astronomy + Astrophysics. Error in affiliation of F. Fontani corrected

The detection of narrow SiO thermal emission toward young outflows has been proposed to be a signature of the magnetic precursor of C-shocks. Recent modeling of the SiO emission across C-shocks predicts variations in the SiO line intensity and line shape at the precursor and intermediate-velocity regimes in only few years. We present high-angular resolution (3.8"x3.3") images of the thermal SiO J=2-1 emission toward the L1448-mm outflow in two epochs (November 2004-February 2005, March-April 2009). Several SiO condensations have appeared at intermediate velocities (20-50 km/s) toward the red-shifted lobe of the outflow since 2005. Toward one of the condensations (clump D), systematic differences of the dirty beams between 2005 and 2009 could be responsible for the SiO variability. At higher velocities (50-80 km/s), SiO could also have experienced changes in its intensity. We propose that the SiO variability toward L1448-mm is due to a real SiO enhancement by young C-shocks at the internal working surface between the jet and the ambient gas. For the precursor regime (5.2-9.2 km/s), several narrow and faint SiO components are detected. Narrow SiO tends to be compact, transient and shows elongated (bow-shock) morphologies perpendicular to the jet. We speculate that these features are associated with the precursor of C-shocks appearing at the interface of the new SiO components seen at intermediate velocities., Comment: 10 pages, 6 figures, 2 tables, accepted for publication in Astrophysical Journal

Variable X-ray and gamma-ray emission is characteristic of the most extreme physical processes in the Universe, and studying the sources of these energetic photons has been a major driver in astronomy for the past 50 years. Here we present multiwavelength observations of a unique gamma-ray selected transient, discovered by Swift, which was accompanied by bright emission across the electromagnetic spectrum, and whose properties are unlike any previously observed source. We pinpoint the event to the center of a small, star-forming galaxy at redshift z=0.3534. Its high-energy emission has lasted much longer than any gamma-ray burst, while its peak luminosity was about 100 times higher than the brightest active galactic nuclei. The association of the outburst with the center of its host galaxy suggests that this phenomenon has its origin in a new, rare mechanism associated with a massive black hole in the nucleus of a galaxy.