Search

(Abridged) Massive stellar embryos are embedded in warm envelopes that provide mass reservoirs for the accretion process onto final stars. Feedback from star formation activities in return impacts the properties of the envelopes, offering us a unique opportunity to investigate star formation processes. We aim to characterise the warm envelopes of proto- or young stellar objects in different evolutionary stages based on the morphology and kinematics of the $^{13}$CO(6-5) emission and to examine their relations with star formation processes. Using the APEX telescope, we obtained maps of the mid-$J$ $^{13}$CO emission with an angular size of 80" x 80" towards 99 massive clumps from the ATLASGAL survey. Our maps are classified based on morphological complexities, and the radial structure of the emission is characterised for simple single-core sources. The velocity centroids of the emission are compared to small- and large-scale gas kinematics, aiming to shed light on the origin of envelope kinematics. $^{13}$CO(6-5) emission is detected towards sources in all stages of high-mass star formation, with a detection rate of 83% for the whole sample. The detection rate, line width, and line peak increase with evolution, and the line luminosity is strongly correlated with $L_{\mathrm{bol}}$ and $M_{\mathrm{clump}}$, indicating that the excitation of $^{13}$CO(6-5) emission is closely related to star formation processes. In addition, the radial distributions of the emission for single-core sources can be well fitted by power-law functions, suggesting a relatively simple envelope structure for many sources. As for the envelope kinematics, linear velocity gradients are common among the single-core sources. Our comparison of kinematics on different scales suggests that the origin of the linear velocity gradients in the warm envelopes is complex and unclear for many sources., Comment: 18 pages, 23 figures, accepted for publication in A&A

Classical Wolf-Rayet stars are descendants of massive OB-type stars that have lost their hydrogen-rich envelopes, and are in the final stages of stellar evolution, possibly exploding as type Ib/c supernovae. It is understood that the mechanisms driving this mass-loss are either strong stellar winds and or binary interactions, so intense studies of these binaries including their evolution can tell us about the importance of the two pathways in WR formation. WR 138 (HD 193077) has a period of just over 4 years and was previously reported to be resolved through interferometry. We report on new interferometric data combined with spectroscopic radial velocities in order to provide a three-dimensional orbit of the system. The precision on our parameters tend to be about an order of magnitude better than previous spectroscopic techniques. These measurements provide masses of the stars, namely $M_{\rm WR} = 13.93\pm1.49M_{\odot}$ and $M_{\rm O} = 26.28\pm1.71M_{\odot}$. The derived orbital parallax agrees with the parallax from \textit{Gaia}, namely with a distance of 2.13 kpc. We compare the system's orbit to models from BPASS, showing that the system likely may have been formed with little interaction but could have formed through some binary interactions either following or at the start of a red supergiant phase, but with the most likely scenario occurring as the red supergiant phase starts for a $\sim 40M_\odot$ star., Comment: accepted to ApJ

The Southern Photometric Local Universe Survey (S-PLUS) is a project to map $\sim9300$ sq deg of the sky using twelve bands (seven narrow and five broadbands). Observations are performed with the T80-South telescope, a robotic telescope located at the Cerro Tololo Observatory in Chile. The survey footprint consists of several large contiguous areas, including fields at high and low galactic latitudes, and towards the Magellanic Clouds. S-PLUS uses fixed exposure times to reach point source depths of about $21$ mag in the $griz$ and $20$ mag in the $u$ and the narrow filters. This paper describes the S-PLUS Data Release 4 (DR4), which includes calibrated images and derived catalogues for over 3000 sq deg, covering the aforementioned area. The catalogues provide multi-band photometry performed with the tools \texttt{DoPHOT} and \texttt{SExtractor} -- point spread function (\PSF) and aperture photometry, respectively. In addition to the characterization, we also present the scientific potential of the data. We use statistical tools to present and compare the photometry obtained through different methods. Overall we find good agreement between the different methods, with a slight systematic offset of 0.05\,mag between our \PSF and aperture photometry. We show that the astrometry accuracy is equivalent to that obtained in previous S-PLUS data releases, even in very crowded fields where photometric extraction is challenging. The depths of main survey (MS) photometry for a minimum signal-to-noise ratio $S/N = 3$ reach from $\sim19.5$ for the bluer bands to $\sim21.5$ mag on the red. The range of magnitudes over which accurate \PSF photometry is obtained is shallower, reaching $\sim19$ to $\sim20.5$ mag depending on the filter. Based on these photometric data, we provide star-galaxy-quasar classification and photometric redshift for millions of objects., Comment: 26 pages, 17 figures, 14 tables, accepted for A&A

We present extensive photometric and spectroscopic observations of the nearby Type Ia supernova (SN) 2023wrk at a distance of about 40 Mpc. The earliest detection of this SN can be traced back to a few hours after the explosion. Within the first few days the light curve shows a bump feature, while the B - V color is blue and remains nearly constant. The overall spectral evolution is similar to that of an SN 1991T/SN 1999aa-like SN Ia, while the C II $\lambda6580$ absorption line appears to be unusually strong in the first spectrum taken at $t \approx -$15.4 days after the maximum light. This carbon feature disappears quickly in subsequent evolution but it reappears at around the time of peak brightness. The complex evolution of the carbon line and the possible detection of Ni III absorption around 4700 {\AA} and 5300 {\AA} in the earliest spectra indicate macroscopic mixing of fuel and ash. The strong carbon lines is likely related to collision of SN ejecta with unbound carbon, consistent with the predictions of pulsational delayed-detonation or carbon-rich circumstellar-matter interaction models. Among those carbon-rich SNe Ia with strong C II $\lambda6580$ absorption at very early times, the line-strength ratio of C II to Si II and the B-V color evolution are found to exhibit large diversity, which may be attributed to different properties of unbound carbon and outward-mixing $^{56}$Ni., Comment: Accepted for publication in the Astrophysical Journal (27 pages, 14 figures, 7 tables)

H$\alpha$ emission is a clear indicator of circumstellar activity in Be stars, historically employed to assess the classical Be star (CBe) population in young open clusters (YOCs). The YOC NGC330 in the Small Magellanic Cloud exhibits a large known fraction of {CBe} stars and was selected for a pilot study to establish a comprehensive methodology for identifying H$\alpha$ emitters in the Magellanic Clouds, encompassing the entire B-type spectral range. Using the SOAR Adaptative Module Imager (SAMI), we investigated the stellar population of NGC330 using multi-band BVRI+H$\alpha$ imaging. We identified H$\alpha$ emitters within the entire V-band range covered by SAMI/SOAR observations ($V\lesssim22$), comprising the complete B-type stellar population and offering a unique opportunity to explore the Be phenomenon across all spectral sub-classes. The stellar radial distribution shows a clear bimodal pattern between the most massive (B5 or earlier) and the lower-mass main-sequence objects (later than B6) within the cluster. The former is concentrated towards the cluster center (showing a dispersion of $\sigma=4.26\pm0.20$ pc), whereas the latter extends across larger radii ($\sigma=10.83\pm0.65$ pc), indicating mass stratification within NGC330. The total fraction of emitters is $4.4\pm0.5\%$, notably smaller than previous estimates from flux- or seeing-limited observations. However, a higher fraction of H$\alpha$ emitters is observed among higher-mass stars ($32.8\pm3.4\%$) than within lower-mass ($4.4\pm0.9\%$). Consequently, the putative CBe population exhibits distinct dynamical characteristics compared to the bulk of the stellar population in NGC330. These findings highlight the significance of the current observations in providing a complete picture of the CBe population in NGC330., Comment: Accepted for publication in the Astrophysical Journal. 24 pages, 17 figures, 3 tables

Changes in the flux and spectrum of Eta Carinae since 1900 have been attributed to the evolution of the central binary by some. Others suggest evolution in the occulting ejecta. The brightness jump in the 1940s, which coincided with the appearance of narrow forbidden emission lines, may have been caused by the clearing and ionization of intervening circumstellar ejecta. The brightening changed at a slower pace up through forty years later. Here we continue earlier studies focused on the long-term showing that the forbidden line emission increased in the early 1990s with no noticeable increase in the brightness of the Homunculus. We interpret that the increase in narrow line emission is due to decreased extinction in the LOS from the central binary to the Weigelt clumps. In 2000, the central stellar core increased in brightness at a faster rate without associated changes in the Homunculus. By 2018, hundreds of narrow-line absorptions from singly-ionized metals in our LOS from Eta Carinae disappeared, thought to be caused by increased ionization of metals. These three events (1990, 2000, and 2018) are explained by the dissipation of circumstellar material within the Homunculus close to the binary. Combining these changes with the steadiness of the Homunculus and the primary winds over the past four decades indicates that circumstellar ejecta in our direction have been cleared., Comment: 26 pages, 8 figures, to appear in ApJ

We report the serendipitous discovery of VVV-WIT-12, an unusual variable source that seems to induce variability in its surrounding nebula. The source belongs to the rare objects that we call WITs (short for What Is This?) discovered within the VISTA Variables in the V\'ia L\'actea (VVV) survey. VVV-WIT-12 was discovered during a pilot search for light echoes from distant Supernovae (SNe) in the Milky Way using the near-IR images of the VVV survey. This source has an extremely red spectral energy distribution, consistent with a very reddened ($A_V \sim 100$ mag) long period variable star ($P\sim1525$ days). Furthermore, it is enshrouded in a nebula that changes brightness and color with time, apparently in synch with the central source variations. The near-IR light curve and complementary follow-up spectroscopy observations are consistent with a variable Young Stellar Object (YSO) illuminating its surrounding nebula. In this case the source periodic variation along the cycles produces an unprecedented light echo in the different regions of the nebula., Comment: Accepted for publication in The Astrophysical Journal Letters

Kilonovae are a rare class of astrophysical transients powered by the radioactive decay of nuclei heavier than iron, synthesized in the merger of two compact objects. Over the first few days, the kilonova evolution is dominated by a large number of radioactive isotopes contributing to the heating rate. On timescales of weeks to months, its behavior is predicted to differ depending on the ejecta composition and merger remnant. However, late-time observations of known kilonovae are either missing or limited. Here we report observations of a luminous red transient with a quasi-thermal spectrum, following an unusual gamma-ray burst of long duration. We classify this thermal emission as a kilonova and track its evolution up to two months after the burst. At these late times, the recession of the photospheric radius and the rapidly-decaying bolometric luminosity ($L_{\rm bol}\propto t^{-2.7\pm 0.4}$) support the recombination of lanthanide-rich ejecta as they cool., Comment: 47 pages, 14 figures, 9 tables; submitted; a minor typo fixed

The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe., Comment: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

W51A is the most active star formation region of the Giant \ion{H}{2} region W51. It harbors the two massive proto-clusters W51e and W51\,IRS2, which are very rare in the Galaxy. We aim to identify the new born massive stars and UCHII regions to derive its distance and age. We performed IFU observations with NIFS+ALTAIR of nine targets in the W51A sub-region. The distance modulus was obtained using the spectral classification in the $K$-band and a reddening law appropriate to the inner Galactic plane. We derived the distance and the spectral types for five of them, ranging from O8 to O9.5, similar to those derived from radio continuum data, except for two sources that we assigned somewhat a later spectral type. We included another seven objects with precise spectral classification from other works, which allowed us to better constrain the distance estimate. Our spectrophotometric distance d=4.80\,$\pm$\,1.27\,kpc is in good agreement with those derived from the Galactic rotation model and trigonometric parallaxes, placing the region near the tangent point of the Sagittarius arm. We conclude that the stars studied in this work have an age spread of 1.5-4 Myr, substantially older than thought to date., Comment: 18 pages, 12 figures and 3 tables. Accepted to ApJ for publication

Westerlund 1 (Wd 1) is one of the most relevant star clusters in the Milky Way to study massive star formation, although it is still poorly known. Here, we used photometric and spectroscopic data to model the eclipsing binary W36, showing that its spectral type is O6.5 III + O9.5 IV, hotter and more luminous than thought before. Its distance $d_{\rm W36}$ $=$ 4.03$\pm$0.25 kpc agrees, within the errors, with three recent Gaia-EDR3-based distances reported in Paper I, Beasor & Davies, and by Negueruela's group. However, they follow different approaches to fix the zero-points for red sources such as those in Wd 1 and to select the best approach, we used an accurate modelling of W36. The weighted mean distance of our parallax (Paper I) and binary distances results in $d_{\rm wd1}$ = 4.05 $\pm$0.20 kpc, with an unprecedented accuracy of 5%. We adopted isochrones based on the Geneva code with supersolar abundances to infer the age of W36B as 6.4 $\pm$ 0.7 Myr. This object seems to be part of the prolific star formation burst represented by OB giants and supergiants that occurred at 7.1 $\pm$ 0.5 Myr ago, which coincides with the recently published PMS isochrone with age 7.2 Myr. Other BA-type luminous evolved stars and Yellow Hypergiants spread in the age range of 8--11 Myr. The four Red Supergiants discussed in paper I represent the oldest population of the cluster with an age of 10.7 $\pm$ 1 Myr. The multiple episodes of star formation in Wd 1 are reminiscent of that reported for the R136/30 Dor LMC cluster., Comment: 15 pages, 6 figures, 2 appendix (with one table)

The Na D absorption doublet in the spectrum of $\eta$ Carinae is complex, with multiple absorption features associated with the Great Eruption (1840s), the Lesser Eruption (1890s), and interstellar clouds. The velocity profile is further complicated by the P Cygni profile originating in the system's stellar winds and blending with the He I $\lambda$5876 profile. The Na D profile contains a multitude of absorption components, including those at velocities of $-$145 km s$^{-1}$, $-$168 km s$^{-1}$, and $+$87 km s$^{-1}$ that we concentrate on in this analysis. Ground-based spectra recorded from 2008 to 2021 show significant variability of the $-$145 km s$^{-1}$ absorption throughout long-term observations. In the high ionization phases of $\eta$ Carinae prior to the 2020 periastron passage, this feature disappeared completely but briefly reappeared across the 2020 periastron, along with a second absorption at $-$168 km s$^{-1}$. Over the past few decades, $\eta$ Car has been gradually brightening demonstrated to be caused by a dissipating occulter. The decreasing absorption of the $-$145 km s$^{-1}$ component, coupled with similar trends seen in absorptions of ultraviolet resonant lines, indicate that this central occulter was possibly a large clump associated with the Little Homunculus or another clump between the Little Homunculus and the star. We also report on a foreground absorption component at $+$87 km s$^{-1}$. Comparison of Na D absorption in the spectra of nearby systems demonstrates that this red-shifted component likely originates in an extended foreground structure consistent with a previous ultraviolet spectral survey in the Carina Nebula., Comment: 8 figures, 13 pages plus online data, accepted to ApJ

We report high-precision X-ray monitoring observations in the 0.4-10 keV band of the luminous, long-period colliding-wind binary Eta Carinae up to and through its most recent X-ray minimum/periastron passage in February 2020. Eta Carinae reached its observed maximum X-ray flux on 7 January 2020, at a flux level of $3.30 \times 10^{-10}$ ergs s$^{-1}$ cm$^{-2}$, followed by a rapid plunge to its observed minimum flux, $0.03 \times 10^{-10}$ ergs s$^{-1}$ cm$^{-2}$ near 17 February 2020. The NICER observations show an X-ray recovery from minimum of only $\sim$16 days, the shortest X-ray minimum observed so far. We provide new constraints of the "deep" and "shallow" minimum intervals. Variations in the characteristic X-ray temperature of the hottest observed X-ray emission indicate that the apex of the wind-wind "bow shock" enters the companion's wind acceleration zone about 81 days before the start of the X-ray minimum. There is a step-like increase in column density just before the X-ray minimum, probably associated with the presence of dense clumps near the shock apex. During recovery and after, the column density shows a smooth decline, which agrees with previous $N_{H}$ measurements made by SWIFT at the same orbital phase, indicating that changes in mass-loss rate are only a few percent over the two cycles. Finally, we use the variations in the X-ray flux of the outer ejecta seen by NICER to derive a kinetic X-ray luminosity of the ejecta of $\sim 10^{41}$ ergs s$^{-1}$ near the time of the "Great Eruption'.

FUV spectra of Eta Car, recorded across two decades with HST/STIS, document multiple changes in resonant lines caused by dissipating extinction in our line of sight. The FUV flux has increased nearly ten-fold which has led to increased ionization of the multiple shells within the Homunculus and photo-destruction of molecular hydrogen. Comparison of observed resonant line profiles with CMFGEN model profiles allows separation of wind-wind collision and shell absorptions from the primary wind, P Cygni profiles.The dissipating occulter preferentially obscured the central binary and interacting winds relative to the very extended primary wind. We are now able to monitor changes in the colliding winds with orbital phase. High velocity transient absorptions occurred across the most recent periastron passage, indicating acceleration of the primary wind by the secondary wind which leads to a downstream, high velocity bowshock that is newly generated every orbital period. There is no evidence of changes in the properties of the binary winds., Comment: 36 pages, 22 figures, accepted Astrophysical Journal

Westerlund 1 (Wd 1) is one of the most massive young star clusters in the Milky Way. Although relevant for star formation and evolution, its fundamental parameters are not yet very well constrained. We aim to derive an accurate distance and provide constraints on the cluster age. We used the photometric and astrometric information available in the Gaia Early Data Release 3 (Gaia-EDR3) to infer its distance of 4.06$^{+0.36}_{-0.34}$ kpc. Modelling of the eclipsing binary system W36, reported in Paper II, led to the distance of 4.03$\pm$0.25 kpc, in agreement with the Gaia-EDR3 distance and, therefore, validating the parallax zero-point correction approach appropriate for red objects. The weighted average distance based on these two methods results in d$_{\rm wd1}$=4.05$\pm$0.20 kpc ($m-M$=13.04$^{+0.11}_{-0.12}$ mag), which has an unprecedented accuracy of 5%. Using the Binary Population and Spectral Synthesis (BPASS) models for the Red Supergiants with solar abundance, we derived an age of 10.7$\pm$1 Myr, in excellent agreement with recent work by Beasor & Davies (10.4$^{+1.3}_{-1.2}$ Myr) based on MIST evolutionary models. In Paper II, the age of W36B was reported to be 6.4$\pm$0.5 Myr, supporting recent claims of a temporal spread of several Myrs for the star-forming process within Wd 1 instead of a single monolithic starburst episode scenario., Comment: 16 pages, 9 figures, accepted for publication at MNRAS on August 18th, 2022

W33A is a well-known example of a high-mass young stellar object showing evidence of a circumstellar disc. We revisited the $K$-band NIFS/Gemini North observations of the W33A protostar using principal components analysis tomography and additional post-processing routines. Our results indicate the presence of a compact rotating disc based on the kinematics of the CO absorption features. The position-velocity diagram shows that the disc exhibits a rotation curve with velocities that rapidly decrease for radii larger than 0\farcs1 ($\sim$250 AU) from the central source, suggesting a structure about four times more compact than previously reported. We derived a dynamical mass of 10.0$^{+4.1}_{-2.2}$ M$_\odot$ for the "disc+protostar" system, about $\sim$33% smaller than previously reported, but still compatible with high-mass protostar status. A relatively compact H$_2$ wind was identified at the base of the large-scale outflow of W33A, with a mean visual extinction of $\sim$63 mag. By taking advantage of supplementary near-infrared maps, we identified at least two other point-like objects driving extended structures in the vicinity of W33A, suggesting that multiple active protostars are located within the cloud. The closest object (Source B) was also identified in the NIFS field of view as a faint point-like object at a projected distance of $\sim$7,000 AU from W33A, powering extended $K$-band continuum emission detected in the same field. Another source (Source C) is driving a bipolar H$_2$ jet aligned perpendicular to the rotation axis of W33A., Comment: to appear on MNRAS. 24 pages, 24 figures and supplementary data available at CDS

Since the start of ALMA observatory operation, new and important chemistry of infrared cold core was revealed. Molecular transitions at millimeter range are being used to identify and to characterize these sources. We have investigated the 231 GHz ALMA archive observations of the infrared dark cloud region C9, focusing on the brighter source that we called as IRDC-C9 Main. We report the existence of two sub-structures on the continuum map of this source: a compact bright spot with high chemistry diversity that we labelled as core, and a weaker and extended one, that we labelled as tail. In the core, we have identified lines of the molecules OCS(19-18), $^{13}$CS(5-4) and CH$_{3}$CH$_{2}$CN, several lines of CH$_{3}$CHO and the k-ladder emission of $^{13}$CH$_{3}$CN.We report two different temperature regions: while the rotation diagram of CH$_{3}$CHO indicates a temperature of 25 K, the rotation diagram of $^{13}$CH$_{3}$CN indicates a warmer phase at temperature of $\sim450$K. In the tail, only the OCS(19-18) and $^{13}$CS(5-4) lines were detected. We used the $Nautilus$ and the \textsc{Radex} codes to estimate the column densities and the abundances. The existence of hot gas in the core of IRDC-C9 Main suggests the presence of a protostar, which is not present in the tail., Comment: 13 pages, 11 figures. Accept in MNRAS

The Perseus Arm is the closest Galactic spiral arm from the Sun, offering an excellent opportunity to study in detail its stellar population. However, its distance has been controversial with discrepancies by a factor of two. Kinematic distances are in the range 3.9-4.2 kpc as compared to 1.9-2.3 kpc from spectrophotometric and trigonometric parallaxes, reinforcing previous claims that this arm exhibits peculiar velocities. We used the astrometric information of a sample of 31 OB stars from the star-forming W3 Complex to identify another 37 W3 members and to derive its distance from their Gaia-DR2 parallaxes with improved accuracy. The Gaia-DR2 distance to the W3 Complex,2.14$^{+0.08}_{-0.07}$ kpc, coincides with the previous stellar distances of $\sim$ 2 kpc. The Gaia-DR2 parallaxes tentatively show differential distances for different parts of the W3 Complex: W3 Main, located to the NE direction, is at 2.30$^{+0.19}_{-0.16}$ kpc, the W3 Cluster (IC 1795), in the central region of the complex, is at 2.17$^{+0.12}_{-0.11}$ kpc, and W3(OH) is at 2.00$^{+0.29}_{-0.23}$ kpc to the SW direction. The W3 Cluster is the oldest region, indicating that it triggered the formation of the other two star-forming regions located at the edges of an expanding shell around the cluster., Comment: 14 pages, 8 figures, accepted for publication in Monthly Notes of the Royal Astronomical Society in May 21th, 2019

The interacting binary Eta Carinae remains one of the most enigmatic massive stars in our Galaxy despite over four centuries of observations. In this work, its light curve from the ultraviolet to the near-infrared is analysed using spatially resolved HST observations and intense monitoring at the La Plata Observatory, combined with previously published photometry. We have developed a method to separate the central stellar object in the ground-based images using HST photometry and applying it to the more numerous ground-based data, which supports the hypothesis that the central source is brightening faster than the almost-constant Homunculus. After detrending from long-term brightening, the light curve shows periodic orbital modulation ($\Delta V$ $\sim$ 0.6 mag) attributed to the wind-wind collision cavity as it sweeps around the primary star and it shows variable projected area to our line-of-sight. Two quasi-periodic components with time scales of 2-3 and 8-10 yr and low amplitude, $\Delta V$ $<$ 0.2 mag, are superimposed on the brightening light-curve, being the only stellar component of variability found, which indicates minimal stellar instability. Moreover, the light curve analysis shows no evidence of `shell ejections' at periastron. We propose that the long-term brightening of the stellar core is due to the dissipation of a dusty clump in front of the central star, which works like a natural coronagraph. Thus, the central stars appear to be more stable than previously thought since the dominant variability originates from a changing circumstellar medium. We predict that the brightening phase, due mainly to dust dissipation, will be completed around 2032 $\pm$ 4 yr, when the star will be brighter than in the 1600's by up to $\Delta V$ $\sim$ 1 mag., Comment: 26 pages, 19 figures. Accepted to appear on MNRAS

High-mass stars are formed within massive molecular clumps, where a large number of stars form close together. The evolution of the clumps with different masses and luminosities is mainly regulated by its high-mass stellar content and the formation of such objects is still not well understood. In this work, we characterise the mid-J CO emission in a statistical sample of 99 clumps (Top100) selected from the ATLASGAL survey that are representative of the Galactic proto-cluster population. High-spatial resolution APEX-CHAMP+ maps of the CO(6-5) and CO(7-6) transitions were obtained and combined with additional single-pointing APEX-FLASH+ spectra of the CO(4-3) line. We study the correlations of the CO line luminosities and profiles for the three CO transitions with the clump properties and investigate if and how they change as a function of the evolution. All sources were detected above 3-$\sigma$ in all three CO transitions and most of the sources exhibit broad CO emission likely associated with molecular outflows. We found that the extension of the mid-J CO emission is correlated with the size of the dust emission traced by the Herschel-PACS 70 $\mu$m maps. The CO line luminosity is correlated with the luminosity and mass of the clumps. However, it does not correlate with the L/M ratio. The dependency of the CO luminosity with the properties of the clumps is steeper for higher-J transitions. Our data seem to exclude that this trend is biased by self-absorption features in the CO emission, but rather suggest that different J transitions arise from different regions of the inner envelope. Moreover, high-mass clumps show similar trends in CO luminosity as lower mass clumps, but are systematically offset towards larger values, suggesting that higher column density and/or temperature (of unresolved) CO emitters are found inside high-mass clumps., Comment: Accepted by A&A: 62 pages, 230 figures. ArXiV: 37 pages, 40 figures

The young stellar cluster Westerlund 1 (Wd 1: $l$=339.6$^\circ$, b=$-$0.4$^\circ$) is one of the most massive in the local Universe, but accurate parameters are pending on better determination of its extinction and distance. Based on our photometry and data collected from other sources, we have derived a reddening law for the cluster line-of-sight representative of the Galactic Plane (-5$^\circ<$b$<$+5$^\circ$) in the window 0.4-4.8 $\mu$m: The power law exponent $\alpha$=2.13$\pm$0.08 is much steeper than those published a decade ago (1.6-1.8) and our index $R_V$=2.50$\pm$0.04 also differs from them, but in very good agreement with recent works based on deep surveys in the inner Galaxy. As a consequence, the total extinction $A_{Ks}$=0.74$ \pm $0.08 ($A_V$=11.40$ \pm$ 2.40) is substantially smaller than previous results(0.91-1.13), part of which ($A_{Ks}$=0.63 or $A_V$=9.66) is from the ISM. The extinction in front of the cluster spans a range of $\Delta A_V\sim$8.7\,mag with a gradient increasing from SW to NE across the cluster face, following the same general trend of warm dust distribution. The map of the $J-Ks$ colour index also shows a trend of reddening in this direction. We measured the equivalent width of the diffuse interstellar band at 8620 $\AA$ (the "GAIA DIB") for Wd 1 cluster members and derived the relation $A_{Ks}$=0.612 $EW$ $-$ 0.191 $EW^2$. This extends the Munari et al. (2008) relation, valid for $E_{B-V}$ $<$ 1, to the non-linear regime ($A_V$ $>$ 4)., Comment: to appear in MNRAS, 16 pages, 22 figures

We present the results from a survey, designed to investigate the accretion process of massive young stellar objects (MYSOs) through near infrared narrow band imaging using the H$_2$ $\nu$=1-0 S(1) transition filter. A sample of 353 Massive Young Stellar Object (MYSO) candidates was selected from the Red MSX Source survey using photometric criteria at longer wavelengths (infrared and submillimeter) and chosen with positions throughout the Galactic Plane. Our survey was carried out at the SOAR Telescope in Chile and CFHT in Hawaii covering both hemispheres. The data reveal that extended H$_2$ emission is a good tracer of outflow activity, which is a signpost of accretion process on young massive stars. Almost half of the sample exhibit extended H$_2$ emission and 74 sources (21\%) have polar morphology, suggesting collimated outflows. The polar-like structures are more likely to appear on radio-quiet sources, indicating these structures occur during the pre-UCHII phase. We also found an important fraction of sources associated with fluorescent H$_2$ diffuse emission that could be due to a more evolved phase. The images also indicate only $\sim$23\% (80) of the sample is associated with extant (young) stellar clusters. These results support the scenario in which massive stars are formed by accretion disks, since the merging of low mass stars would not produce outflow structures., Comment: to appear on MNRAS. 18 pages, 86 figures and supplementary data available at CDS