Search

An established method measuring the hydrogen ionisation fraction in shock excited gas is the BE99 method, which utilises six bright forbidden emission lines of [SII]6716, 6731, [NII]6548, 6583, and [OI]6300, 6363. We aim to extent the BE99 method by including more emission lines in the blue and near-infrared part of the spectrum ($\lambda$ = 3500-11000A), and considering higher hydrogen ionisation fractions ($x_e > 0.3$). In addition, we investigate how a non-equilibrium state of the gas and the presence of extinction influence the BE99 technique. We find that plenty additional emission line ratios can in principle be exploited as extended curves (or stripes) in the ($x_e, T_e$)-diagram. If the BE99 equilibrium is reached and extinction is corrected for, all stripes overlap in one location in the ($x_e, T_e$)-diagram indicating the existing gas parameters. The application to the Par Lup 3-4 outflow shows that the classical BE99 lines together with the [NI]5198+5200 lines do not meet in one locationin the ($x_e, T_e$)-diagram. This indicates that the gas parameters derived from the classical BE99 method are not fully consistent with other observed line ratios. A multi-line approach is necessary to determine the gas parameters. From our analysis we derive $n_e \sim$ 45 000 cm^-3 - 53000 cm^-3 , $T_e$ = 7600K - 8000K, and $x_e \sim$ 0.027 - 0.036 for the Par Lup 3-4 outflow. For the 244-440 Proplyd we were able to use the line ratios of [SII]6716+6731, [OI]6300+6363, and [OII]7320, 7330 in the BE99 diagram to estimate the ionisation fraction at knot E3 ($x_e = 0.58 \pm 0.05$). In conclusion, exploiting new line ratios reveals more insights on the state of the gas. Our analysis indicates, however, that a multi-line approach is more robust in deriving gas parameters, especially for high density gas., Comment: 12 pages, 10 figures, accepted to A&A

Understanding the interstellar and potentially circumstellar extinction in the sight lines of classical T Tauri stars is an important ingredient for constructing reliable spectral energy distributions, which catalyze protoplanetary disk chemistry, for example. Therefore, some attempts of measuring $A_{V}$ toward individual stars have been made using partly different wavelength regimes and different underlying assumptions. We used strong lines of Ly{\alpha} fluorescent H2 and derived the extinction based on the assumption of optically thin transitions. We investigated a sample of 72 classical T Tauri stars observed with the Hubble Space Telescope in the framework of the ULLYSES program. We computed $A_{V}$ and $R_{V}$ values for the 34 objects with sufficient data quality and an additionally $A_{V}$ value for the canonical $R_{V}$ = 3.1 value. Our results agree largely with values obtained from optical data. Moreover, we confirm the degeneracy between $A_{V}$ and $R_{V}$ and present possibilities to break this. Finally, we discuss whether the assumption of optical thin lines is valid., Comment: 12 pages, 6 figures, accepted to A&A

Characterising stellar jet asymmetries is key to setting robust constraints on jet launching models and improving our understanding of the underlying mechanisms behind jet launching. We aim to characterise the asymmetric properties of the bipolar jet coming from the Classical T Tauri Star Th 28. We combined data from integral field spectroscopy with VLT/MUSE and high-resolution spectra from VLT/X-shooter to map the optical emission line ratios in both jet lobes. We carried out a diagnostic analysis of these ratios to compare the density, electron temperature, and ionisation fraction within both lobes. The mass accretion rate was derived from the emission lines at the source and compared with the mass outflow rate derived for both lobes, using the estimated densities and measured [O I]6300 and [S II ]6731 luminosities. The blue-shifted jet exhibits a significantly higher electron temperature and moderately higher ionisation fraction than the red-shifted jet. In contrast to previous studies, we also estimated higher densities in the blue-shifted jet by a factor of ~2. These asymmetries are traced to within 160 au of the source in the line ratio maps. We estimate the mass accretion rate onto the central star and compare this with estimates of the mass outflow rates through each side of the jet. The emission line maps and diagnostic results suggest that the jet asymmetries originate close to the source and are likely to be intrinsic to the jet. Furthermore, the combined dataset offers access to a broad array of accretion tracers. In turn, this enables a more accurate estimation of the mass accretion rate, revealing a value of Macc that is higher by a factor >350 than would otherwise be determined. Supplementary figures and tables are available via a public Zenodo repository (doi:10.5281/zenodo.13373809)., Comment: 19 pages, 16 figures, 4 tables

The main goal of this study is to screen the PENELLOPE/UVES targets for outflow activity and find microjets via spectro-astrometry in, e.g., the [OI]$\lambda$6300 line. In total, 34 T\,Tauri stars of the PENELLOPE survey have been observed with the high resolution slit spectrograph UVES in three different slit positions rotated by $120^\text{o}$. Our spectro-astrometric analysis in the [OI]$\lambda$6300 wind line reveals two newly discovered microjets associated with Sz\,103 and XX\,Cha. Both microjets have an extent of about $0.04$ arcseconds, that is, $<10\,\text{au}$, and we confined their orientation by the three slit observations. Furthermore, we confirm the binary nature of VW\,Cha and CVSO\,109. We present (further) evidence that DK\,Tau\,B and CVSO\,104\,A are spectroscopic binaries. Sz\,115 is tentatively a spectroscopic binary. The origin of the LVC, that is, MHD winds versus photoevaporative winds, of the Sz\,103 and XX\,Cha microjets remains unclear., Comment: 12 pages, 8 figures, accepted for publication in A&A

Accretion bursts from low-mass young stellar objects (YSOs) are known for many decades. In recent years, the first accretion bursts of massive YSOs (MYSOs) have been observed. These phases of intense protostellar growth are of particular importance for studying massive star formation. Bursts of MYSOs are accompanied by flares of Class II methanol masers (hereafter masers), caused by an increase in exciting mid-infrared (MIR) emission. The G323.46$-$0.08 (hereafter G323) event extends the small sample of known MYSO bursts. Maser observations of the MYSO G323 show evidence of a flare, which was presumed to be caused by an accretion burst. This should be verified with IR data. We used time-dependent radiative transfer (TDRT) to characterize the heating and cooling timescales for eruptive MYSOs and to infer the main burst parameters. The G323 accretion burst is confirmed. It reached its peak in late 2013/early 2014 with a Ks-band increase of 2.5mag. TDRT indicates that the duration of the thermal afterglow in the far-infrared (FIR) can exceed the burst duration by years. The latter was proved by SOFIA observations, which indicate a flux increase of $(14.2\pm4.6)$% at $70\, \rm \mu m$ and $(8.5\pm6.1)$% at $160\, \mu$m in 2022 (2 years after the burst end). A one-sided light echo emerged that was propagating into the interstellar medium. The G323 burst is probably the most energetic MYSO burst observed so far. Within $8.4 \rm \, yrs$, an energy of $(0.9\pm_{0.8}^{2.5}) \times 10^{47}\,\rm erg$ was released. The short timescale points to the accretion of a compact body, while the burst energy corresponds to an accumulated mass of at least $(7\pm_{6}^{20})\,M_{Jup}$ and possibly even more if the protostar is bloated. In this case, the accretion event might have triggered protostellar pulsations, which give rise to the observed maser periodicity.

We present the first results of the JWST program PROJECT-J (PROtostellar JEts Cradle Tested with JWST ), designed to study the Class I source HH46 IRS and its outflow through NIRSpec and MIRI spectroscopy (1.66 to 28 micron). The data provide line-images (~ 6.6" in length with NIRSpec, and up to 20" with MIRI) revealing unprecedented details within the jet, the molecular outflow and the cavity. We detect, for the first time, the red-shifted jet within ~ 90 au from the source. Dozens of shock-excited forbidden lines are observed, including highly ionized species such as [Ne III] 15.5 micron, suggesting that the gas is excited by high velocity (> 80 km/s) shocks in a relatively high density medium. Images of H2 lines at different excitations outline a complex molecular flow, where a bright cavity, molecular shells, and a jet-driven bow-shock interact with and are shaped by the ambient conditions. Additional NIRCam 2 micron images resolve the HH46 IRS ~ 110 au binary system and suggest that the large asymmetries observed between the jet and the H2 wide angle emission could be due to two separate outflows being driven by the two sources. The spectra of the unresolved binary show deep ice bands and plenty of gaseous lines in absorption, likely originating in a cold envelope or disk. In conclusion, JWST has unraveled for the first time the origin of the HH46 IRS complex outflow demonstrating its capability to investigate embedded regions around young stars, which remain elusive even at near-IR wavelengths., Comment: 28 pages, 15 figures, Accepted for publication on The Astrophysical Journal (9 April 2024)

We report new Very Large Array high-resolution observations of the radio jet from the outbursting high-mass star S255IR~NIRS3. The images at 6, 10, and 22.2 GHz confirm the existence of a new lobe emerging to the SW and expanding at a mean speed of ~285 km/s, about half as fast as the NE lobe. The new data allow us to reproduce both the morphology and the continuum spectrum of the two lobes with the model already adopted in our previous studies. We conclude that in all likelihood both lobes are powered by the same accretion outburst. We also find that the jet is currently fading down, recollimating, and recombining.

Observations indicate that the accretion process in star formation may occur through accretion outbursts. This phenomenon has also now been detected in a few young massive (proto)stars (>8 Msun). The recent outburst at radio wavelengths of the massive (proto)star S255 NIRS3 has been interpreted by us as expansion of a thermal jet, fed by the infalling material. To follow up on our previous study and confirm our interpretation, we monitored the source for more than 1 yr in six bands from 1.5 GHz to 45.5 GHz and, after ~1.5 yr, with the Atacama Large Millimeter/submillimeter Array at two epochs, which made it possible to detect the proper motions of the jet lobes. The prediction of our previous study is confirmed by the new results. The radio jet is found to expand, while the flux, after an initial exponential increase, appears to stabilise and eventually decline. The radio flux measured during our monitoring is attributed to a single NE lobe, However, from 2019 a second lobe has been emerging to the SW, probably powered by the same accretion outburst, although with a delay of at least a couple of years. Flux densities at >6 GHz were satisfactorily fitted with a jet model, whereas those below 6 GHz are clearly underestimated by the model. This indicates that non-thermal emission becomes dominant at long wavelengths. Our results suggest that thermal jets can be a direct consequence of accretion events, when yearly flux variations are detected. The end of the accretion outburst is mirrored in the radio jet, as ~1 yr after the onset of the radio outburst, the inner radius of the jet began to increase while the jet mass stopped growing, as expected if the powering mechanism of the jet is quenched. Our findings support a tight connection between accretion and ejection in massive stars, consistent with a formation process involving a disk-jet system similar to that of low-mass stars.

Using observations with e-MERLIN and the VLA, together with archival data from ALMA, we obtain high-resolution radio images of two binary YSOs: L1551 IRS 5 and L1551 NE, covering a wide range of frequencies from 5 - 336 GHz, and resolving emission from the radio jet on scales of only ~15 au. By comparing these observations to those from a previous epoch, it is shown that there is a high degree of variability in the free-free emission from the jets of these sources. In particular, the northern component of L1551 IRS 5 shows a remarkable decline in flux density of a factor of ~5, suggesting that the free-free emission of this source has almost disappeared. By fitting the spectra of the sources, the ionised mass-loss rates of the jets are derived and it is shown that there is significant variability of up to a factor of ~6 on timescales of ~20 years. Using radiative transfer modelling, we also obtained a model image for the jet of the southern component of L1551 IRS 5 to help study the inner region of the ionised high-density jet. The findings favour the X-wind model launched from a very small innermost region., Comment: 13 pages, 7 figures, accepted for publication in A&A

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

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

$Context.$ The YSO Th 28 possesses a highly collimated jet, which clearly exhibits an asymmetric brightness of its jet lobes at optical and NIR wavelengths. There may be asymmetry in the jet plasma parameters in opposite jet lobes (e.g. electron density, temperature, and outflow velocity). $Aims.$ We examined the Th 28 jet in a 3"x3" where the jet material is collimated and accelerated. Our goal is to map the morphology and determine its physical parameters to determine the physical origin of such asymmetries. $Methods.$ We present $JHK$-spectra of Th 28 obtained with the SINFONI on the (VLT, ESO) in June-July 2015. $Results.$ The [Fe II] emission originates in collimated jet lobes. Two new axial knots are detected at 1" in the blue lobe and 1".2 in the red lobe. The H$_2$ radiation is emitted from an extended region with a radius of $\gtrsim270$ au, which is perpendicular to the jet. The PV diagrams of the bright H$_2$ lines reveal faint H$_2$ emission along both jet lobes as well. The compact and faint H I emission (Pa$\beta$ and Br$\gamma$) comes from two regions, namely from a spherical region around the star and from the jet lobes. The size of the jet launching region is derived as 0".015 ($\sim$3 au at 185 pc), and the initial opening angle of the Th 28 jet is $\sim28^0$, which makes this jet substantially less collimated than most jets from other CTTs. $Conclusions.$ The emission in [Fe II], H$_2$, and H I lines suggests a morphology in which the ionised gas in the disc appears to be disrupted by the jet. The resolved disc-like H$_2$ emission most likely arises in the disc atmosphere from shocks caused by a radial uncollimated wind. The asymmetry of the [Fe II] photocentre shifts with respect to the jet source arises in the immediate vicinity of the driving source of Th28 and suggests that the observed brightness asymmetry is intrinsic as well., Comment: 18 pages, 14 figures

Recent studies have shown that 6.7 GHz methanol maser flares can be a powerful tool for verifying the mechanisms of maser production and even the specific signatures of accretion rate changes in the early stages of high-mass star formation. We characterize the spatial structure and evolution of methanol and water masers during a flare of methanol maser emission at 6.7 GHz in the HMYSO G24.33$+$0.14. VLBA was used to image the 6.7 and 12.2 GHz methanol and 22.2 GHz water vapor masers at three epochs guided by monitoring the methanol line with the Torun 32m telescope. The 6.7 GHz maser maps were also obtained with the EVN and LBA during the flare. WISE data were used to find correlations between the 6.7 GHz maser and IR fluxes. The 6.7 GHz methanol maser cloudlets are distributed over $\sim$3500 au, and the morphology of most of them is stable although their brightness varies following the course of the total flux density on a timescale of two months. The 12.2 GHz methanol maser cloudlets cover an area an order of magnitude smaller than that of 6.7 GHz emission, and both transitions emerge from the same masing gas. The 22.2 GHz maser cloudlets lie in the central region and show a systematic increase in brightness and moderate changes in size and orientation, together with the velocity drift of the strongest cloudlet during two months of the VLBI observing period. Time lag estimates imply the propagation of changes in the physical conditions of the masing region with a subluminal speed ($\sim$0.3c). A tight correlation of IR (4.6$\mu$m) and 6.7 GHz flux densities is found, supporting the radiative pumping model. Comparison with the 230 GHz ALMA data indicates that the methanol masers are distributed in the inner part of the rotating disk, whereas the 22.2 GHz emission traces the compact inner component of the bipolar outflow or a jet structure., Comment: Accepted by A\&A

HD163296 is a Herbig Ae star which drives a bipolar knotty jet with a total length of ~6000au. Strong evidence exists that the disk of HD163296 harbors planets. Studies have shown that the presence of companions around jet-driving stars could affect the morphology of the jets. This includes a `wiggling' of the jet axis and a periodicity in the positions of the jet knots. In this study we investigate the morphology (including the jet width and axis position) and proper motions of the HD163296 jets, and use our results to better understand the whole system.This study is based on optical integral-field spectroscopy observations obtained with VLT/MUSE in 2017. Using spectro-images and position velocity diagrams extracted from the MUSE data cube, we investigated the number and positions of the jet knots. A comparison was made to X-Shooter data collected in 2012 and the knot proper motions were estimated. The jet width and jet axis position with distance from the star were studied from the extracted spectro-images. We observe the merging of knots and identify two previously undetected knots. Measurements of the jet axis position reveal a similar pattern of deviation in all forbidden emission lines along the first 20 arc seconds of the jets. This result is interpreted as being due to asymmetric shocks and not due to a wiggling of the jet axis. The number of new knots detected and their positions challenge the 16-year knot ejection periodicity proposed in prior studies, arguing for a more complicated jet system than was previously assumed. We use the non-detection of a jet axis wiggling to rule out companions with a mass $>$0.1~\Msun\ and orbits between 1~au and 35~au. Any object inferred at these distances using other methods must be a brown dwarf or planet, otherwise it would have impacted the jet axis position. Both the precession and orbital motion scenarios are considered.

The Orion molecular cloud complex harbours the nearest GMCs and site of high-mass star formation. Its YSO populations are thoroughly characterized. The region is therefore a prime target for the study of star formation. Here, we verify the performance of the SuperCAM 64 pixel heterodyne array on APEX. We give a descriptive overview of a set of wide-field CO(3-2) spectral cubes obtained towards the Orion GMC complex, aimed at characterizing the dynamics and structure of the extended molecular gas in diverse regions of the clouds, ranging from very active sites of clustered star formation in Orion B to comparatively quiet regions in southern Orion A. We present a 2.7 square degree (130pc$^2$) mapping survey in the CO(3-2) transition, obtained using SuperCAM on APEX at an angular resolution of 19'' (7600AU or 0.037pc at a distance of 400pc), covering L1622, NGC2071, NGC2068, OriB9, NGC2024, and NGC2023 in Orion B, and the southern part of the L1641 cloud in Orion A. We describe CO integrated emission and line moment maps and position-velocity diagrams and discuss a few sub-regions in some detail. Evidence for expanding bubbles is seen with lines splitting into double components, most prominently in NGC2024, where we argue that the bulk of the molecular gas is in the foreground of the HII region. High CO(3-2)/CO(1-0) line ratios reveal warm CO along the western edge of Orion B in the NGC2023/NGC2024 region facing the IC434 HII region. Multiple, well separated radial velocity components seen in L1641-S suggest that it consists of a sequence of clouds at increasingly larger distances. We find a small, spherical cloud - the 'Cow Nebula' globule - north of NGC2071. We trace high velocity line wings for the NGC2071-IR outflow and the NGC2024 CO jet. The protostellar dust core FIR4 (rather than FIR5) is the true driving source of the NGC2024 monopolar outflow., Comment: Accepted for publication in Astronomy and Astrophysics

[abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $\mu$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $\mu$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR. 24 pages, 16 figures

The Low-Frequency Array (LOFAR) Long-Baseline Calibrator Survey (LBCS) was conducted between 2014 and 2019 in order to obtain a set of suitable calibrators for the LOFAR array. In this paper we present the complete survey, building on the preliminary analysis published in 2016 which covered approximately half the survey area. The final catalogue consists of 30006 observations of 24713 sources in the northern sky, selected for a combination of high low-frequency radio flux density and flat spectral index using existing surveys (WENSS, NVSS, VLSS, and MSSS). Approximately one calibrator per square degree, suitable for calibration of $\geq$ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 degrees and away from the Galactic plane, with a considerably lower density south of this point due to relative difficulty in selecting flat-spectrum candidate sources in this area of the sky. Use of the VLBA calibrator list, together with statistical arguments by comparison with flux densities from lower-resolution catalogues, allow us to establish a rough flux density scale for the LBCS observations, so that LBCS statistics can be used to estimate compact flux densities on scales between 300 mas and 2 arcsec, for sources observed in the survey. The LBCS can be used to assess the structures of point sources in lower-resolution surveys, with significant reductions in the degree of coherence in these sources on scales between 2 arcsec and 300 mas. The LBCS survey sources show a greater incidence of compact flux density in quasars than in radio galaxies, consistent with unified schemes of radio sources. Comparison with samples of sources from interplanetary scintillation (IPS) studies with the Murchison Widefield Array (MWA) shows consistent patterns of detection of compact structure in sources observed both interferometrically with LOFAR and using IPS., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR

The INT Galactic Plane Survey (IGAPS) is the merger of the optical photometric surveys, IPHAS and UVEX, based on data from the Isaac Newton Telescope (INT) obtained between 2003 and 2018. These capture the entire northern Galactic plane within the Galactic coordinate range, -5

Context: Th 28 is a Classical T Tauri star in the Lupus 3 cloud which drives an extended bipolar jet. Previous studies of the inner jet identified signatures of rotation around the outflow axis, a key result for theories of jet launching. Thus this is an important source in which to investigate the poorly understood jet launching mechanism. We investigate the morphology and kinematics of the Th 28 micro-jets with the aim of characterizing their structure and outflow activity, using optical integral-field spectroscopy observations obtained with VLT/MUSE. We use spectro-imaging and position-velocity maps to investigate the kinematic and morphological features of the jet, and obtain a catalogue of emission lines in which the jet is visible. A Lucy-Richardson deconvolution procedure is used to differentiate the structure of the inner micro-jet region. Spatial profiles extracted perpendicular to the jet axis are fitted to investigate the jet width, opening angle and the evolution of the jet axis. We confirm the previously identified knot HHW$_{2}$ within the red-shifted jet and identify three additional knots in each lobe for the first time. We also find [O III]$\lambda$5007 emission from the blue-shifted micro-jet including the knot closest to the star. Proper motions for the innermost knots on each side are estimated and we show that new knots are ejected on an approximate timescale of 10-15 years. The jet axis centroids show a point-symmetric wiggle within the inner portion of both micro-jets indicating precession. We use the jet shape to measure a precession period of 8 years, with a half-opening angle < 0.6$^{\circ}$. This may provide an alternative explanation for the rotation signatures previously reported. We find the jet shape to be compatible with precession due to a brown dwarf companion orbiting at a separation $\leq$ 0.3 au., Comment: Accepted by Astronomy & Astrophysics; 21 pages, 23 figures

In recent years, thanks to new facilities such as LOFAR capable of sensitive observations, much work has been done on the detection of stellar radio emission at low frequencies. Such emission has commonly been shown to be coherent emission, generally attributed to electron-cyclotron maser emission, and has usually been detected from main-sequence M dwarfs. Here we report the first detection of coherent emission at low frequencies from T Tauri stars, which are known to be associated with high levels of stellar activity. Using LOFAR, we have detected several bright radio bursts at 150 MHz from two weak-line T Tauri stars: KPNO-Tau 14 and LkCa 4. All of the bursts have high brightness temperatures ($10^{13} - 10^{14}\ \mathrm{K}$) and high circular polarization fractions (60 - 90 \%), indicating that they must be due to a coherent emission mechanism. This could be either plasma emission or electron-cyclotron maser (ECM) emission. Due to the exceptionally high brightness temperatures seen in at least one of the bursts ($\geq 10^{14}\ \mathrm{K}$), and the high circular polarization levels, it seems unlikely that plasma emission could be the source and so ECM is favoured as the most likely emission mechanism. Assuming this is the case, the required magnetic field in the emission regions would be 40 - 70 G. We determine that the most likely method of generating ECM emission is plasma co-rotation breakdown in the stellar magnetosphere. There remains the possibility, however, it could be due to an interaction with an orbiting exoplanet., Comment: 16 pages, 9 figures, accepted for publication in A&A

We present SOFIA/FIFI-LS observations of three Class 0 and one Class I outflows (Cep E, HH 1, HH 212, and L1551 IRS5) in the far-infrared [O I]63mum and [O I]145mum transitions. Spectroscopic [O I]63mum maps enabled us to infer the spatial extent of warm, low-excitation atomic gas within these protostellar outflows. If proper shock conditions prevail, the instantaneous mass-ejection rate is directly connected to the [O I]63mum luminosity. In order to unravel evolutionary trends, we analysed a set of 14 Class 0/I outflow sources that were spatially resolved in the [O I]63mum emission. We compared these data with a sample of 72 Class 0/I/II outflow sources that have been observed with Herschel (WISH, DIGIT, WILL, GASPS surveys) without spatially resolving the [O I]63mum line., Comment: 14 pages, 7 figures

The evolution of young stars and disks is driven by the interplay of several processes, notably accretion and ejection of material. Critical to correctly describe the conditions of planet formation, these processes are best probed spectroscopically. About five-hundred orbits of the Hubble Space Telescope (HST) are being devoted in 2020-2022 to the ULLYSES public survey of about 70 low-mass (M<2Msun) young (age<10 Myr) stars at UV wavelengths. Here we present the PENELLOPE Large Program that is being carried out at the ESO Very Large Telescope (VLT) to acquire, contemporaneous to HST, optical ESPRESSO/UVES high-resolution spectra to investigate the kinematics of the emitting gas, and UV-to-NIR X-Shooter medium-resolution flux-calibrated spectra to provide the fundamental parameters that HST data alone cannot provide, such as extinction and stellar properties. The data obtained by PENELLOPE have no proprietary time, and the fully reduced spectra are made available to the whole community. Here, we describe the data and the first scientific analysis of the accretion properties for the sample of thirteen targets located in the Orion OB1 association and in the sigma-Orionis cluster, observed in Nov-Dec 2020. We find that the accretion rates are in line with those observed previously in similarly young star-forming regions, with a variability on a timescale of days of <3. The comparison of the fits to the continuum excess emission obtained with a slab model on the X-Shooter spectra and the HST/STIS spectra shows a shortcoming in the X-Shooter estimates of <10%, well within the assumed uncertainty. Its origin can be either a wrong UV extinction curve or due to the simplicity of this modelling, and will be investigated in the course of the PENELLOPE program. The combined ULLYSES and PENELLOPE data will be key for a better understanding of the accretion/ejection mechanisms in young stars., Comment: Accepted for publication on Astronomy & Astrophysics. 15 pages + appendix, language edited version

Class II methanol masers are signs of massive young stellar objects (MYSOs). Recent findings show that MYSO accretion bursts cause flares of these masers. Thus, maser monitoring can be used to identify such bursts. Burst-induced SED changes provide valuable information on a very intense phase of high-mass star formation. In mid-January 2019, a maser flare of the MYSO G358.93-0.03 was reported. ALMA and SMA imaging resolved the core of the star forming region and proved the association of the masers with the brightest continuum source MM1. However, no significant flux rise of the (sub)mm dust continuum was found. Thus, we performed NIR imaging with GROND and IFU spectroscopy with FIFI-LS aboard SOFIA to detect possible counterparts to the (sub)mm sources, and compare their photometry to archival measurements. The comparison of pre-burst and burst SEDs is of crucial importance to judge whether a luminosity increase due to the burst is present and if it triggered the maser flare. The FIR fluxes of MM1 measured with FIFI-LS exceed those from Herschel significantly, which clearly confirms the presence of an accretion burst. The second epoch data, taken about 16 months later, still show increased fluxes. Our RT modeling yielded major burst parameters and suggests that the MYSO features a circumstellar disk which might be transient. From the multi-epoch SEDs, conclusions on heating and cooling time-scales could be drawn. Circumstances of the burst-induced maser relocation have been explored. The verification of the accretion burst from G358 is another confirmation that Class II methanol maser flares represent an alert for such events. The few events known to date already indicate that there is a broad range in burst strength and duration as well as environmental characteristics. The G358 event is the shortest and least luminous MYSO accretion burst so far., Comment: 22 pages, 11 figures, accepted by A&A. Abstract abridged for arxiv submission

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

We present SOFIA/FIFI-LS observations of five prototypical, low-mass Class I outflows (HH111, SVS13, HH26, HH34, HH30) in the far-infrared [OI]63mum and [OI]145mum transitions. The obtained spectroscopic [OI]63mum and [OI]145mum maps enable us to study the spatial extent of warm, low-excitation atomic gas within outflows driven by Class I protostars. These [OI] maps may potentially allow us to measure the mass-loss rates ($\dot{M}_\text{jet}$) of this warm component of the atomic jet., Comment: 24 pages, 41 figures, accepted for publication in A&A

Young accreting stars drive outflows that collimate into jets, which can be seen hundreds of au from their driving sources. Accretion and outflow activity cease with system age, and it is believed that magneto-centrifugally launched disk winds are critical agents in regulating accretion through the protoplanetary disk. Protostellar jets are well studied in classical T~Tauri stars ($M_\star\lesssim2\,M_\odot$), while few nearby ($d\lesssim150\,$pc) intermediate-mass stars ($M_\star=2-10\,M_\odot$), known as Herbig Ae/Be stars, have detected jets. We report VLT/MUSE observations of the Herbig~Ae/Be star HD~100546 and the discovery of a protostellar jet. The jet is similar in appearance to jets driven by low-mass stars and compares well with the jet of HD~163296, the only other known optical jet from a nearby Herbig~Ae/Be star. We derive a (one-sided) mass-loss rate in the jet of $\log \dot{M}_{jet} \sim -9.5$ (in $M_\odot$\,yr$^{-1}$) and a ratio of outflow to accretion of roughly $3\times10^{-3}$, which is lower than that of CTTS jets. The discovery of the HD~100546 jet is particularly interesting because the protoplanetary disk around HD~100546 shows a large radial gap, spiral structure, and might host a protoplanetary system. A bar-like structure previously seen in H$\alpha$ with VLT/SPHERE shares the jet position angle, likely represents the base of the jet, and suggests a jet-launching region within about 2\,au. We conclude that the evolution of the disk at radii beyond a few au does not affect the ability of the system to launch jets., Comment: 6 pages, 6 figures, accepted for publication in A&A letters

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$\sigma$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $\alpha^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations., Comment: 9 pages, 2 figures, accepted for publication in MNRAS

The INT Galactic Plane Survey (IGAPS) is the merger of the optical photometric surveys, IPHAS and UVEX, based on data from the Isaac Newton Telescope (INT) obtained between 2003 and 2018. Here, we present the IGAPS point source catalogue. It contains 295.4 million rows providing photometry in the filters, i, r, narrow-band Halpha, g and U_RGO. The IGAPS footprint fills the Galactic coordinate range, |b| < 5deg and 30deg < l < 215deg. A uniform calibration, referred to the Pan-STARRS system, is applied to g, r and i, while the Halpha calibration is linked to r and then is reconciled via field overlaps. The astrometry in all 5 bands has been recalculated on the Gaia DR2 frame. Down to i ~ 20 mag (Vega system), most stars are also detected in g, r and Halpha. As exposures in the r band were obtained within the IPHAS and UVEX surveys a few years apart, typically, the catalogue includes two distinct r measures, r_I and r_U. The r 10sigma limiting magnitude is ~21, with median seeing 1.1 arcsec. Between ~13th and ~19th magnitudes in all bands, the photometry is internally reproducible to within 0.02 magnitudes. Stars brighter than r=19.5 have been tested for narrow-band Halpha excess signalling line emission, and for variation exceeding |r_I-r_U| = 0.2 mag. We find and flag 8292 candidate emission line stars and over 53000 variables (both at >5sigma confidence). The 174-column catalogue will be available via CDS Strasbourg., Comment: 28 pages, 22 figures

Massive young stellar object (MYSOs) form during the collapse of high-mass pre-stellar cores, where infalling molecular material is accreted through a centrifugally-balanced accretion disc that is subject to efficient gravitational instabilities. In the resulting fragmented accretion disc of the MYSO, gaseous clumps and low-mass stellar companions can form, which will influence the future evolution of massive protostars in the Hertzsprung-Russell diagram. We perform dust continuum radiative transfer calculations and compute synthetic images of disc structures modelled by the gravito-radiation-hydrodynamics simulation of a forming MYSO, in order to investigate the Atacama Large Millimeter/submillimeter Array (ALMA) observability of circumstellar gaseous clumps and forming multiple systems. Both spiral arms and gaseous clumps located at ~a few 100 au from the protostar can be resolved by interferometric ALMA Cycle 7 C43-8 and C43-10 observations at band 6 (1.2 mm), using a maximal 0.015" beam angular resolution and at least 10-30 min exposure time for sources at distances of 1-2 kpc. Our study shows that substructures are observable regardless of their viewing geometry or can be inferred in the case of an edge-viewed disc. The observation probability of the clumps increases with the gradually increasing efficiency of gravitational instability at work as the disc evolves. As a consequence, large discs around MYSOs close to the zero-age-main-sequence line exhibit more substructures than at the end of the gravitational collapse. Our results motivate further observational campaigns devoted to the close surroundings of the massive protostars S255IR-NIRS3 and NGC 6334I-MM1, whose recent outbursts are a probable signature of disc fragmentation and accretion variability., Comment: Published at MNRAS, link: http://adsabs.harvard.edu/abs/2019MNRAS.487.4473M Previous version (2019arXiv190602015M) was identical but had an incomplete title

Variability in Young Stellar Objects (YSOs) is one of their primary characteristics. Long-term, multi-filter, high-cadence monitoring of large YSO samples is the key to understand the partly unusual light-curves that many of these objects show. Here we introduce and present the first results of the HOYS-CAPS citizen science project which aims to perform such monitoring for nearby (d<1kpc) and young (age<10Myr) clusters and star forming regions, visible from the northern hemisphere, with small telescopes. We have identified and characterised 466 variable (413 confirmed young) stars in 8 young, nearby clusters. All sources vary by at least 0.2mag in V, have been observed at least 15 times in V, R and I in the same night over a period of about 2yrs and have a Stetson index of larger than 1. This is one of the largest samples of variable YSOs observed over such a time-span and cadence in multiple filters. About two thirds of our sample are classical T-Tauri stars, while the rest are objects with depleted or transition disks. Objects characterised as bursters show by far the highest variability. Dippers and objects whose variability is dominated by occultations from normal interstellar dust or dust with larger grains (or opaque material) have smaller amplitudes. We have established a hierarchical clustering algorithm based on the light-curve properties which allows the identification of the YSOs with the most unusual behaviour, and to group sources with similar properties. We discuss in detail the light-curves of the unusual objects V2492Cyg, V350Cep and 2MASSJ21383981+5708470., Comment: 13 pages, 9 figures, 3 tables, accepted for publication in MNRAS, full version with full appendix available at http://astro.kent.ac.uk/~df/papers.html

With their wide fields of view and often relatively long coverage of any position in the sky in imaging survey mode, modern radio telescopes provide a data stream that is naturally suited to searching for rare transients. However, Radio Frequency Interference (RFI) can show up in the data stream in similar ways to such transients, and thus the normal pre-treatment of filtering RFI (flagging) may also remove astrophysical transients from the data stream before imaging. In this paper we investigate how standard flagging affects the detectability of such transients by examining the case of transient detection in an observing mode used for Low Frequency Array (LOFAR; \citep{LOFAR}) surveys. We quantify the fluence range of transients that would be detected, and the reduction of their SNR due to partial flagging. We find that transients with a duration close to the integration sampling time, as well as bright transients with durations on the order of tens of seconds, are completely flagged. For longer transients on the order of several tens of seconds to minutes, the flagging effects are not as severe, although part of the signal is lost. For these transients, we present a modified flagging strategy which mitigates the effect of flagging on transient signals. We also present a script which uses the differences between the two strategies, and known differences between transient RFI and astrophysical transients, to notify the observer when a potential transient is in the data stream., Comment: 31 pages, 6 figures, accepted for publication in Astronomy and Computing

Context. Recent observations of the massive young stellar object S255 NIRS 3 have revealed a large increase in both methanol maser flux density and IR emission, which have been interpreted as the result of an accretion outburst, possibly due to instabilities in a circumstellar disk. This indicates that this type of accretion event could be common in young/forming early-type stars and in their lower mass siblings, and supports the idea that accretion onto the star may occur in a non-continuous way. Aims. As accretion and ejection are believed to be tightly associated phenomena, we wanted to confirm the accretion interpretation of the outburst in S255 NIRS 3 by detecting the corresponding burst of the associated thermal jet. Methods. We monitored the radio continuum emission from S255 NIRS 3 at four bands using the Karl G. Jansky Very Large Array. The millimetre continuum emission was also observed with both the Northern Extended Millimeter Array of IRAM and the Atacama Large Millimeter/submillimeter Array. Results. We have detected an exponential increase in the radio flux density from 6 to 45 GHz starting right after July 10, 2016, namely about 13 months after the estimated onset of the IR outburst. This is the first ever detection of a radio burst associated with an IR accretion outburst from a young stellar object. The flux density at all observed centimetre bands can be reproduced with a simple expanding jet model. At millimetre wavelengths we infer a marginal flux increase with respect to the literature values and we show this is due to free-free emission from the radio jet. Abridged., Comment: 14 pages, 9 figures, 3 tables, accepted by Astronomy & Astrophysics

We present LOFAR high-band data over the frequency range 115-189 MHz for the X-ray binary SS 433, obtained in an observing campaign from 2013 February - 2014 May. Our results include a deep, wide-field map, allowing a detailed view of the surrounding supernova remnant W 50 at low radio frequencies, as well as a light curve for SS 433 determined from shorter monitoring runs. The complex morphology of W 50 is in excellent agreement with previously published higher-frequency maps; we find additional evidence for a spectral turnover in the eastern wing, potentially due to foreground free-free absorption. Furthermore, SS 433 is tentatively variable at 150 MHz, with both a debiased modulation index of 11 per cent and a $\chi^2$ probability of a flat light curve of $8.2 \times 10^{-3}$. By comparing the LOFAR flux densities with contemporaneous observations carried out at 4800 MHz with the RATAN-600 telescope, we suggest that an observed $\sim$0.5-1 Jy rise in the 150-MHz flux density may correspond to sustained flaring activity over a period of approximately six months at 4800 MHz. However, the increase is too large to be explained with a standard synchrotron bubble model. We also detect a wealth of structure along the nearby Galactic plane, including the most complete detection to date of the radio shell of the candidate supernova remnant G 38.7-1.4. This further demonstrates the potential of supernova remnant studies with the current generation of low-frequency radio telescopes., Comment: 19 pages, 8 figures, 4 tables, accepted for publication in MNRAS

Solar-mass stars form via circumstellar disk accretion (disk-mediated accretion). Recent findings indicate that this process is likely episodic in the form of accretion bursts, possibly caused by disk fragmentation. Although it cannot be ruled out that high-mass young stellar objects (HMYSOs; $M>$8 M$_\odot$, $L_{bol}>$5$\times$10$^3$ L$_\odot$) arise from the coalescence of their low-mass brethren, latest results suggest that they more likely form via disks. Accordingly, disk-mediated accretion bursts should occur. Here we report on the discovery of the first disk-mediated accretion burst from a $\sim$20 M$_\odot$ HMYSO. Our near-infrared images show the brightening of the central source and its outflow cavities. Near-infrared spectroscopy reveals emission lines typical of accretion bursts in low-mass protostars, but orders of magnitude more luminous. Moreover, the energy released and the inferred mass-accretion rate are also orders of magnitude larger. Our results identify disk accretion as the common mechanism of star formation across the entire stellar mass spectrum., Comment: 27 pages, 3 figures. Published in Nature Physics, 2016, Vol. 13, pag. 276

Massive OB stars are critical to the ecology of galaxies, and yet our knowledge of OB stars in the Milky Way, fainter than $V \sim 12$, remains patchy. Data from the VST Photometric H$\alpha$ Survey (VPHAS+) permit the construction of the first deep catalogues of blue excess-selected OB stars, without neglecting the stellar field. A total of 14900 candidates with 2MASS cross-matches are blue-selected from a 42 square-degree region in the Galactic Plane, capturing the Carina Arm over the Galactic longitude range $282^{\circ} \lesssim \ell \lesssim 293^{\circ}$. Spectral energy distribution fitting is performed on these candidates' combined VPHAS+ $u,g,r,i$ and 2MASS $J,H,K$ magnitudes. This delivers: effective temperature constraints, statistically separating O from early-B stars; high-quality extinction parameters, $A_0$ and $R_V$ (random errors typically $< 0.1$). The high-confidence O-B2 candidates number 5915 and a further 5170 fit to later B spectral type. Spectroscopy of 276 of the former confirms 97% of them. The fraction of emission line stars among all candidate B stars is 7--8% . Greyer ($R_V > 3.5$) extinction laws are ubiquitous in the region, over the distance range 2.5--3 kpc to $\sim$10~kpc. Near prominent massive clusters, $R_V$ tends to rise, with particularly large and chaotic excursions to $R_V \sim 5$ seen in the Carina Nebula. The data reveal a hitherto unnoticed association of 108 O-B2 stars around the O5If$+$ star LSS 2063 ($\ell = 289.77^{\circ}$, $b = -1.22^{\circ}$). Treating the OB star scale-height as a constant within the thin disk, we find an orderly mean relation between extinction ($A_0$) and distance in the Galactic longitude range, $287.6^{\circ} < \ell < 293.5^{\circ}$, and infer the subtle onset of thin-disk warping. A halo around NGC 3603, roughly a degree in diameter, of $\sim$500 O-B2 stars with $4 < A_0 (\rm{mag}) < 7$ is noted., Comment: Accepted for publication in MNRAS. 26 pages, 20 figures, 4 tables, 1 appendix

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree., Comment: Accepted by Astronomy & Astrophysics. Error in figure 6 corrected

The eclipses of certain types of binary millisecond pulsars (i.e. `black widows' and `redbacks') are often studied using high-time-resolution, `beamformed' radio observations. However, they may also be detected in images generated from interferometric data. As part of a larger imaging project to characterize the variable and transient sky at radio frequencies <200 MHz, we have blindly detected the redback system PSR J2215+5135 as a variable source of interest with the Low-Frequency Array (LOFAR). Using observations with cadences of 2 weeks - 6 months, we find preliminary evidence that the eclipse duration is frequency dependent ($\propto \nu^{-0.4}$), such that the pulsar is eclipsed for longer at lower frequencies, in broad agreement with beamformed studies of other similar sources. Furthermore, the detection of the eclipses in imaging data suggests an eclipsing medium that absorbs the pulsed emission, rather than scattering it. Our study is also a demonstration of the prospects of finding pulsars in wide-field imaging surveys with the current generation of low-frequency radio telescopes., Comment: 9 pages, 4 figures, 4 tables, accepted for publication in MNRAS

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range., Comment: 35 pages, 11 figures, updated version: Pierre Auger Observatory data ICRC 2015 added to Fig 2

We present simultaneous photometric and spectroscopic observations of seven young and highly variable M dwarfs in star forming regions in Orion, conducted in 4 observing nights with FORS2 at ESO/VLT. All seven targets show significant photometric variability in the I-band, with amplitudes between 0.1-0.8 mag, The spectra, however, remain remarkably constant, with spectral type changes less than 0.5 subtypes. Thus, the brightness changes are not caused by veiling that 'fills in' absorption features. Three objects in the $\sigma$ Ori cluster (age $\sim$3 Myr) exhibit strong H$\alpha$ emission and H$\alpha$ variability, in addition to the continuum variations. Their behaviour is mostly consistent with the presence of spots with temperature of $\sim300$ K above the photosphere and filling factors between 0.2-0.4, in contrast to typical hot spots observed in more massive stars. The remaining targets near $\epsilon$ Ori, likely to be older, show eclipse-like lightcurves, no significant H$\alpha$ activity and are better represented by variable extinction due to circumstellar material. Interestingly, two of them show no evidence of infrared excess emission. Our study shows that high-amplitude variability in young very low mass stars can be caused by different phenomena than in more massive T Tauri stars and can persist when the disk has disappeared and accretion has ceased., Comment: 19 pages, 8 figures, accepted

Only a few well characterized very low-mass M dwarfs are known today. Our understanding of M dwarfs is vital as these are the most common stars in our solar neighborhood. We aim to characterize the properties of a rare F+dM stellar system for a better understanding of the low-mass end of the Hertzsprung-Russel diagram. We used photometric light curves and radial velocity follow-up measurements to study the binary. Spectro- scopic analysis was used in combination with isochrone fitting to characterize the primary star. The primary star is an early F-type main-sequence star with a mass of (1.493 +- 0.073) Msun and a radius of (1.474 +- 0.040) Rsun. The companion is an M dwarf with a mass of (0.188 +- 0.014) Msun and a radius of (0.234 +- 0.009) Rsun. The orbital period is (1.35121 +- 0:00001)d. The secondary star is among the lowest-mass M dwarfs known to date. The binary has not reached a 1:1 spin-orbit synchronization. This indicates a young main-sequence binary with an age below ~250 Myrs. The mass-radius relation of both components are in agreement with this finding., Comment: 7 pages, 5 figures, accepted by AJ

LS I +61 303 is a gamma-ray binary that exhibits an outburst at GHz frequencies each orbital cycle of $\approx$ 26.5 d and a superorbital modulation with a period of $\approx$ 4.6 yr. We have performed a detailed study of the low-frequency radio emission of LS I +61 303 by analysing all the archival GMRT data at 150, 235 and 610 MHz, and conducting regular LOFAR observations within the Radio Sky Monitor (RSM) at 150 MHz. We have detected the source for the first time at 150 MHz, which is also the first detection of a gamma-ray binary at such a low frequency. We have obtained the light-curves of the source at 150, 235 and 610 MHz, all of them showing orbital modulation. The light-curves at 235 and 610 MHz also show the existence of superorbital variability. A comparison with contemporaneous 15-GHz data shows remarkable differences with these light-curves. At 15 GHz we see clear outbursts, whereas at low frequencies we see variability with wide maxima. The light-curve at 235 MHz seems to be anticorrelated with the one at 610 MHz, implying a shift of $\sim$ 0.5 orbital phases in the maxima. We model the shifts between the maxima at different frequencies as due to changes in the physical parameters of the emitting region assuming either free-free absorption or synchrotron self-absorption, obtaining expansion velocities for this region close to the stellar wind velocity with both mechanisms., Comment: 12 pages, 10 figures, accepted for publication in MNRAS

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration., Comment: 23 pages, 16 figures, 5 tables, accepted for publication in MNRAS

Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 {\deg} $\pm$ 25 {\deg}. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles., Comment: 10 pages, 12 figures, accepted for publication in A&A (27/11/2015) - abstract edited because of limited characters

We present first results from a LOFAR census of non-recycled pulsars. The census includes almost all such pulsars known (194 sources) at declinations ${\rm Dec}> 8^\circ$ and Galactic latitudes $|{\rm Gb}| > 3^\circ$, regardless of their expected flux densities and scattering times. Each pulsar was observed for $\geq 20$ minutes in the contiguous frequency range of 110--188 MHz. Full-Stokes data were recorded. We present the dispersion measures, flux densities, and calibrated total intensity profiles for the 158 pulsars detected in the sample. The median uncertainty in census dispersion measures ($1.5 \times 10^{-3}$ pc cm$^{-3}$) is ten times smaller, on average, than in the ATNF pulsar catalogue. We combined census flux densities with those in the literature and fitted the resulting broadband spectra with single or broken power-law functions. For 48 census pulsars such fits are being published for the first time. Typically, the choice between single and broken power-laws, as well as the location of the spectral break, were highly influenced by the spectral coverage of the available flux density measurements. In particular, the inclusion of measurements below 100 MHz appears essential for investigating the low-frequency turnover in the spectra for most of the census pulsars. For several pulsars, we compared the spectral indices from different works and found the typical spread of values to be within 0.5--1.5, suggesting a prevailing underestimation of spectral index errors in the literature. The census observations yielded some unexpected individual source results, as we describe in the paper. Lastly, we will provide this unique sample of wide-band, low-frequency pulse profiles via the European Pulsar Network Database., Comment: 39 pages, 26 figures, accepted by A&A. Fixed a number of typos, rewritten DM variations section (conclusions unchanged). Arxiv pdf contains tables that will be published as electronic only by A&A

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/)., Comment: 38 pages, 11 figures, 5 tables, A&A in press, updated with editorial corrections

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs., Comment: 23 pages, 19 figures. Accepted for publication in Astronomy & Astrophysics. MSSS Verification Field images and catalog data may be downloaded from http://vo.astron.nl