This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

Ring-like distributions of the 6.7 GHz methanol maser spots at milliarcsecond scales represent a family of molecular structures of unknown origin associated with high-mass young stellar objects (HMYSOs). We aim to study G23.657-00.127, which has a nearly circular ring of the 6.7 GHz methanol masers, and is the most suitable target to test hypotheses on the origin of the maser rings. The European Very Long Baseline Interferometry Network (EVN) was used at three epochs spanning 10.3 yr to derive the spatio-kinematical structure of the 6.7 GHz methanol maser emission in the target. The maser cloudlets, lying in a nearly symmetric ring, expand mainly in the radial direction with a mean velocity of 3.2 km s$^{-1}$. There is an indication that the radial component of the velocity increases with cloudlet's distance from the ring centre. The tangential component does not show any clear evidence for rotation of the cloudlets or any relationship with distance from the ring centre. The blue-shifted masers may hint at an anticlockwise rotation of cloudlets in the southern part of the ring. The nearly circular structure of the ring clearly persisted for more than 10 yr. Interferometric data demonstrated that about one quarter of cloudlets show significant variability in their brightness, although the overall spectrum was non-variable in single-dish studies. Taking into account the three-dimensional motion of the maser cloudlets and their spatial distribution along a small ring, we speculate about two possible scenarios where the methanol masers trace either a spherical outflow arising from an (almost) edge-on disc, or a wide angle wind at the base of a protostellar jet. The latter is associated with near- and mid-infrared emission detected towards the ring. High angular resolution images of complementary (thermal) tracers are needed to interpret the environment of methanol masers.

GX 9+9 (4U 1728-16) is a low mass X-ray binary (LMXB) source harboring a neutron star. Although it belongs to the subclass of the bright Atoll sources together with GX 9+1, GX 3+1, and GX 13+1, its broadband spectrum is poorly studied and apparently does not show reflection features in the spectrum. To constrain the continuum well and verify whether a relativistic smeared reflection component is present, we analyze the broadband spectrum of GX 9+9 using {\it BeppoSAX} and \textit{XMM-Newton} spectra covering the 0.3-40 keV energy band. We fit the spectrum adopting a model composed of a disk-blackbody plus a Comptonized component whose seed photons have a blackbody spectrum (Eastern Model). A statistically equivalent model is composed of a Comptonized component whose seed photons have a disk-blackbody distribution plus a blackbody that mimics a saturated Comptonization likely associated with a boundary layer (Western model). Other trials did not return a good fit. The spectrum of GX 9+9 was observed in a soft state and its luminosity is $2.3 \times 10^{37}$ erg s$^{-1}$ assuming a distance to the source of 5 kpc. In the Eastern Model scenario, we find the seed-photon temperature and electron temperature of the Comptonized component to be $1.14^{+0.10}_{-0.07}$ keV and $2.80^{+0.09}_{-0.04}$ keV, respectively, while the optical depth of the Comptonizing corona is $8.9\pm0.4$. The color temperature of the inner accretion disk is $0.86^{+0.08}_{-0.02}$ keV and $0.82 \pm 0.02$ keV for the {\it BeppoSAX} and \textit{XMM-Newton} spectrum, respectively. In the Western Model scenario, instead, we find that the seed-photon temperature is $0.87 \pm 0.07$ keV and $1.01 \pm 0.08$ keV for the {\it BeppoSAX} and \textit{XMM-Newton} spectrum, respectively. (Abridged abstract)

We report on the phase-coherent timing analysis of the accreting millisecond X-ray pulsar IGR J17591-2342, using Neutron Star Interior Composition Explorer (NICER) data taken during the outburst of the source between 2018 August 15 and 2018 October 17. We obtain an updated orbital solution of the binary system. We investigate the evolution of the neutron star spin frequency during the outburst, reporting a refined estimate of the spin frequency and the first estimate of the spin frequency derivative ($\dot{\nu} \sim -7\times 10^{-14}$ Hz s$^{-1}$), confirmed independently from the modelling of the fundamental frequency and its first harmonic. We further investigate the evolution of the X-ray pulse phases adopting a physical model that accounts for the accretion material torque as well as the magnetic threading of the accretion disc in regions where the Keplerian velocity is slower than the magnetosphere velocity. From this analysis we estimate the neutron star magnetic field $B_{eq} = 2.8(3)\times10^{8}$ G. Finally, we investigate the pulse profile dependence on energy finding that the observed behaviour of the pulse fractional amplitude and lags as a function of energy are compatible with a thermal Comptonisation of the soft photons emitted from the neutron star caps.

We use anisotropic fluid cosmology to describe the present, dark energy-dominated, universe. Similarly to what has been proposed for galactic dynamics, the anisotropic fluid gives an effective description of baryonic matter, dark energy and their possible interaction, without assuming the presence of dark matter. The resulting anisotropic fluid spacetime naturally generates inhomogeneities at small scales, triggered by an anisotropic stress, and could therefore be responsible for structure formation at these scales. Solving the cosmological equations, we show that the dynamics of the scale factor $a$ is described by usual FLRW cosmology and decouples completely from that describing inhomogeneities. We assume that the cosmological anisotropic fluid inherits the equation of state from that used to describe galaxy rotation curves. We show that, in the large scale regime, the fluid can be described as a generalized Chaplygin gas and fits well the distance modulus experimental data of type Ia supernovae, thus correctly modelling the observed accelerated expansion of the universe. Conversely, in the small scale regime, we use cosmological perturbation theory to derive the power spectrum $P(k)$ for mass density distribution. At short wavelengths, we find a $1/k^4$ behaviour, in good accordance with the observed correlation function for matter distribution at small scales.

We wish to perform a statistical study of the location and motion of individual 22 GHz water maser cloudlets, characterized by sizes that are within a few au, with respect to the radio thermal emission from young stellar objects (YSO). For this purpose, we have been carrying out the Protostellar Outflows at the EarliesT Stages (POETS) survey of a sample (38) of high-mass YSOs. The water maser positions and three-dimensional (3D) velocities were determined through Very Long Baseline Array observations with accuracies of a few milliarcsec (mas) and a few km/s, respectively. The position of the ionized core of the protostellar wind, marking the YSO, was determined through sensitive continuum Jansky Very Large Array observations with a typical error of 20 mas. The statistic of the separation of the water masers from the radio continuum shows that 84% of the masers are found within 1000 au from the YSO and 45% of them are within 200 au. Therefore, we can conclude that the 22 GHz water masers are a reliable proxy for the YSO position. The distribution of maser luminosity is strongly peaked towards low values, indicating that about half of the maser population is still undetected with the current Very Long Baseline Interferometry detection thresholds of 50-100 mJy/beam. Next-generation, sensitive radio interferometers will exploit these weak masers for an improved sampling of the velocity and magnetic fields around the YSOs. The average direction of the water maser proper motions provides a statistically-significant estimate for the orientation of the jet emitted by the YSO: 55% of the maser proper motions are directed on the sky within an angle of 30 deg from the jet axis. Finally, we show that our measurements of 3D maser velocities statistically support models in which water maser emission arises from planar shocks with propagation direction close to the plane of the sky.

SAX J1748.9-2021 is a transient accretion powered millisecond X-ray pulsar located in the Globular cluster NGC 6440. We report on the spectral and timing analysis of SAX J1748.9-2021 performed on AstroSat data taken during its faint and short outburst of 2017. We derived the best-fitting orbital solution for the 2017 outburst and obtained an average local spin frequency of 442.361098(3) Hz. The pulse profile obtained from 3-7 keV and 7-20 keV energy bands suggest constant fractional amplitude ~0.5% for fundamental component, contrary to previously observed energy pulse profile dependence. Our AstroSat observations revealed the source to be in a hard spectral state. The 1-50 keV spectrum from SXT and LAXPC on-board AstroSat can be well described with a single temperature blackbody and thermal Comptonization. Moreover, we found that the combined spectra from XMM-Newton (EPIC-PN) and AstroSat (SXT+LAXPC) indicated the presence of reflection features in the form of iron (Fe K${\alpha}$) line that we modeled with the reflection model xillvercp. One of the two X-ray burst observed during the AstroSat/LAXPC observation showed hard X-ray emission (>30 keV) due to Compton up-scattering of thermal photons by the hot corona. Time resolved analysis performed on the bursts revealed complex evolution in emission radius of blackbody for second burst suggestive of mild photospheric radius expansion.

GrailQuest (Gamma-ray Astronomy International Laboratory for Quantum Exploration of Space-Time) is an ambitious astrophysical mission concept that uses a fleet of small satellites, whose scientific objectives are discussed below. Within Quantum Gravity theories, different models for space-time quantisation predict an energy dependent speed for photons. Although the predicted discrepancies are minuscule, Gamma-Ray Bursts, occurring at cosmological distances, could be used to detect this signature of space-time granularity with a new concept of modular observatory of huge overall collecting area consisting in a fleet of small satellites in low orbits, with sub-microsecond time resolution and wide energy band (keV-MeV). The enormous number of collected photons will allow to effectively search these energy dependent delays. Moreover, GrailQuest will allow to perform temporal triangulation of high signal-to-noise impulsive events with arc-second positional accuracies: an extraordinary sensitive X-ray/Gamma all-sky monitor crucial for hunting the elusive electromagnetic counterparts of Gravitational Waves. A pathfinder of GrailQuest is already under development through the HERMES (High Energy Rapid Modular Ensemble of Satellites) project: a fleet of six 3U cube-sats to be launched by the end of 2021.

We compile and analyze ~200 trigonometric parallaxes and proper motions of molecular masers associated with very young high-mass stars. These measurements strongly suggest that the Milky Way is a four-arm spiral. Fitting log-periodic spirals to the locations of the masers, allows us to significantly expand our view of the structure of the Milky Way. We present an updated model for its spiral structure and incorporate it into our previously published parallax-based distance-estimation program for sources associated with spiral arms. Modeling the three-dimensional space motions yields estimates of the distance to the Galactic center, Ro = 8.15 +/- 0.15 kpc, the circular rotation speed at the Sun's position, To = 236 +/- 7 km/s, and the nature of the rotation curve. Our data strongly constrain the full circular velocity of the Sun, To + Vsun = 247 +/- 4 km/s, and its angular velocity, (To + Vsun)/Ro = 30.32 +/- 0.27 km/s/kpc. Transforming the measured space motions to a Galactocentric frame which rotates with the Galaxy, we find non-circular velocity components typically about 10 km/s. However, near the Galactic bar and in a portion of the Perseus arm, we find significantly larger non-circular motions. Young high-mass stars within 7 kpc of the Galactic center have a scale height of only 19 pc and, thus, are well suited to define the Galactic plane. We find that the orientation of the plane is consistent with the IAU-defined plane to within +/-0.1 deg., and that the Sun is offset toward the north Galactic pole by Zsun = 5.5 +/- 5.8 pc. Accounting for this offset places the central supermassive black hole, Sgr A*, in the midplane of the Galaxy. Using our improved Galactic parameters, we predict the Hulse-Taylor binary pulsar to be at a distance of 6.54 +/- 0.24 kpc, assuming its orbital decay from gravitational radiation follows general relativity.

The goal of the Protostellar Outflows at the EarliesT Stages (POETS) survey is to image the disk-outflow interface on scales of 10-100 au in a statistically significant sample (36) of luminous young stellar objects (YSO), targeting both the molecular and ionized components of the outflows. The outflow kinematics is studied at milliarcsecond scales through VLBI observations of the 22 GHz water masers. We employed the JVLA at 6, 13, and 22 GHz in the A- and B-Array configurations to determine the spatial structure and the spectral index of the radio continuum emission. In about half of the targets, the water masers observed at separation <= 1000 au from the YSOs trace either or both of these kinematic structures: 1) a spatially elongated distribution oriented at close angle with the direction of collimation of the maser proper motions (PM), and 2) a linear LSR velocity (Vlsr) gradient across the YSO position. The kinematic structure (1) is readily interpreted in terms of a protostellar jet, as confirmed in some targets via the comparison with independent observations of the YSO jets, in thermal (continuum and line) emissions, reported in the literature. The kinematic structure (2) is interpreted in terms of a disk-wind (DW) seen almost edge-on on the basis of several pieces of evidence: first, it is invariably directed perpendicular to the YSO jet; second, it agrees in orientation and polarity with the Vlsr gradient in thermal emissions (when reported in the literature) identifying the YSO disk at scales of <= 1000~au; third, the PMs of the masers delineating the Vlsr gradients hint at flow motions at a speed of 10-20 km/s directed at large angles with the disk midplane. In the remaining targets, the maser PMs are not collimated but rather tend to align along two almost perpendicular directions, and could originate in DW-jet systems slightly inclined (<= 30 deg) with respect to edge-on.

The X-ray transient eclipsing source MXB 1659-298 went into outburst in 1999 and 2015. During these two outbursts the source was observed by XMM-Newton, nuSTAR, and Swift/XRT. Using these observations, we studied the broadband spectrum of the source to constrain the continuum components and to verify whether it had a reflection component, as is observed in other X-ray eclipsing transient sources. We combined the available spectra to study the soft and hard state of the source in the 0.45-55 keV energy range. We report a reflection component in the soft and hard state. The direct emission in the soft state can be modeled with a thermal component originating from the inner accretion disk plus a Comptonized component associated with an optically thick corona surrounding the neutron star. On the other hand, the direct emission in the hard state is described only by a Comptonized component with a temperature higher than 130 keV; this component is associated with an optically thin corona. We observed narrow absorption lines from highly ionized ions of oxygen, neon, and iron in the soft spectral state. We investigated where the narrow absorption lines form in the ionized absorber. The equivalent hydrogen column density associated with the absorber is close to $6 \times 10^{23}$ cm$^{-2}$ and $1.3 \times 10^{23}$ cm$^{-2}$ in the soft and hard state, respectively.

It is important to determine if massive stars form via disc accretion, like their low-mass counterparts. Theory and observation indicate that protostellar jets are a natural consequence of accretion discs and are likely to be crucial for removing angular momentum during the collapse. However, massive protostars are typically rarer, more distant and more dust enshrouded, making observational studies of their jets more challenging. A fundamental question is whether the degree of ionisation in jets is similar across the mass spectrum. Here we determine an ionisation fraction of $\sim5-12\%$ in the jet from the massive protostar G35.20-0.74N, based on spatially coincident infrared and radio emission. This is similar to the values found in jets from lower-mass young stars, implying a unified mechanism of shock ionisation applies in jets across most of the protostellar mass spectrum, up to at least $\sim10$ solar masses.

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

The X-ray transient eclipsing source MXB 1659-298 went in outburst in 1999 and 2015, respectively, during which it was observed by XMM-Newton, NuSTAR and Swift. Using these observations we studied the broadband spectrum of the source to constrain the continuum components and to verify the presence of a reflection component. We analysed the soft and hard state of the source, finding that the soft state can be modelled with a thermal component associated with the inner accretion disc plus a Comptonised component. A smeared reflection component and the presence of an ionised absorber are also requested in the best-fit model. On the other hand, the direct continuum emission in the hard state can be described by a Comptonised component with a temperature larger than 150 keV. Also in this case a reflection component and a ionised absorber are observed.

PSR J1023+0038 is the first millisecond pulsar discovered to pulsate in the visible band; such a detection took place when the pulsar was surrounded by an accretion disk and also showed X-ray pulsations. We report on the first high time resolution observational campaign of this transitional pulsar in the disk state, using simultaneous observations in the optical (TNG, NOT, TJO), X-ray (XMM-Newton, NuSTAR, NICER), infrared (GTC) and UV (Swift) bands. Optical and X-ray pulsations were detected simultaneously in the X-ray high intensity mode in which the source spends $\sim$ 70% of the time, and both disappeared in the low mode, indicating a common underlying physical mechanism. In addition, optical and X-ray pulses were emitted within a few km, had similar pulse shape and distribution of the pulsed flux density compatible with a power-law relation $F_{\nu} \propto \nu^{-0.7}$ connecting the optical and the 0.3-45 keV X-ray band. Optical pulses were detected also during flares with a pulsed flux reduced by one third with respect to the high mode; the lack of a simultaneous detection of X-ray pulses is compatible with the lower photon statistics. We show that magnetically channeled accretion of plasma onto the surface of the neutron star cannot account for the optical pulsed luminosity ($\sim 10^{31}$ erg/s). On the other hand, magnetospheric rotation-powered pulsar emission would require an extremely efficient conversion of spin-down power into pulsed optical and X-ray emission. We then propose that optical and X-ray pulses are instead produced by synchrotron emission from the intrabinary shock that forms where a striped pulsar wind meets the accretion disk, within a few light cylinder radii away, $\sim$ 100 km, from the pulsar.

MAXI J1820+070, a black hole candidate first detected in early March 2018, was observed by XMM-Newton during the outburst rise. In this letter we report on the spectral and timing analysis of the XMM-Newton X-ray and UV data, as well as contemporaneous X-ray data from the Swift satellite. The X-ray spectrum is well described by a hard thermal Comptonization continuum. The XMM-Newton X-ray light curve shows a pronounced dipping interval, and spectral analysis indicates that it is caused by a moderately ionized partial covering absorber. The XMM-Newton/OM U-filter data does not reveal any signs of the 17 hr orbital modulation that was seen later on during the outburst decay. The UV/X-ray cross correlation function shows a complex shape, with a peak at positive lags of about 4 seconds and a pre-cognition dip at negative lags, which is absent during the X-ray dipping episode. Such shape could arise if the UV emission comes partially from synchrotron self-Compton emission near the black hole, as well as from reprocessing of the X-rays in the colder accretion disc further out.

Context. Since the discovery of the first Accreting Millisecond X-ray Pulsar SAX J1808.4-3658 in 1998, the family of these sources kept growing on. Currently, it counts 22 members. All AMXPs are transients with usually very long quiescence periods, implying that mass accretion rate in these systems is quite low and not constant. Moreover, for at least three sources, a non-conservative evolution was also proposed. Aims. Our purpose is to study the long term averaged mass-accretion rates in all the Accreting Millisecond X-ray Pulsars discovered so far, to investigate a non-conservative mass-transfer scenario. Methods. We calculated the expected mass-transfer rate under the hypothesis of a conservative evolution based on their orbital periods and on the (minimum) mass of the secondary (as derived from the mass function), driven by gravitational radiation and/or magnetic braking. Using this theoretical mass-transfer, we determined the expected accretion luminosity of the systems. Thus, we achieved the lower limit to the distance of the sources by comparing the computed theoretical luminosity and the observed flux averaged over a time period of 20 years. Then, the lower limit to the distance of the sources has been compared to the value of the distance reported in literature to evaluate how reasonable is the hypothesis of a conservative mass-transfer. Results. Based on a sample of 18 sources, we found strong evidences of a non-conservative mass-transfer for five sources, for which the estimated distance lower limits are higher than their known distances. We also report hints for mass outflows in other six sources. The discrepancy can be fixed under the hypothesis of a non-conservative mass-transfer in which a fraction of the mass transferred onto the compact object is swept away from the system, likely due to the (rotating magnetic dipole) radiation pressure of the pulsar.

The source X 1822-371 is an eclipsing compact binary system with a period close to 5.57 hr and an orbital period derivative $\dot{P}_{\rm orb}$ of 1.51(7)$\times 10^{-10}$ s s$^{-1}$. The very large value of $\dot{P}_{\rm orb}$ is compatible with a super-Eddington mass transfer rate from the companion star, as suggested by X-ray and optical data. The XMM-Newton observation taken in 2017 allows us to update the orbital ephemeris and verify whether the orbital period derivative has been stable over the last 40 yr. We added to the X-ray eclipse arrival times from 1977 to 2008 two new values obtained from the RXTE and XMM-Newton observations performed in 2011 and 2017, respectively. We estimated the number of orbital cycles and the delays of our eclipse arrival times spanning 40 yr using as reference time the eclipse arrival time obtained from the Rossi-XTE observation taken in 1996. Fitting the delays with a quadratic model, we found an orbital period $P_{\rm orb}=5.57062957(20)$ hr and a $\dot{P}_{\rm orb}$ value of $1.475(54) \times 10^{-10}$ s s$^{-1}$. The addition of a cubic term to the model does not significantly improve the quality of the fit. We also determined a spin-period value of $P_{\rm spin}=0.5915669(4)$ s and its first derivative $\dot{P}_{\rm spin}= -2.595(11) \times 10^{-12}$ s s$^{-1}$. The obtained results confirm the scenario of a super-Eddington mass transfer rate; we also exclude a gravitational coupling between the orbit and the change in the oblateness of the companion star triggered by the nuclear luminosity of the companion star.

Context: XB 1916-053 is a low mass X-ray binary system (LMXB) hosting a neutron star (NS) and showing periodic dips. The spectrum of the persistent emission was modeled with a blackbody component having a temperature between 1.31 and 1.67 keV and with a Comptonization component with an electron temperature of 9.4 keV and a photon index $\Gamma$ between 2.5 and 2.9. The presence of absorption features associated with highly ionized elements suggested the presence of partially ionized plasma in the system. Aims: In this work we performed a study of the spectrum of XB 1916-053, which aims to shed light on the nature of the seed photons that contribute to the Comptonization component. Methods: We analyzed three Suzaku observations of XB 1916-053: the first was performed in November 2006 and the others were carried out in October 2014. We extracted the persistent spectra from each observation and combined the spectra of the most recent observations, obtaining a single spectrum with a higher statistic. We also extracted and combined the spectra of the dips observed during the same observations. Results: On the basis of the available data statistics, we infer that the scenario in which the corona Comptonizes photons emitted both by the innermost region of the accretion disk and the NS surface is not statistically relevant with respect to the case in which only photons emitted by the NS surface are Comptonized. We find that the source is in a soft spectral state in all the analyzed observations. We detect the K$\alpha$ absorption lines of \ion{Fe}{xxv} and \ion{Fe}{xxvi}, which have already been reported in literature, and for the first time the K$\beta$ absorption lines of the same ions. We also detect an edge at 0.876 keV, which is consistent with a \ion{O}{viii} K absorption edge. (Abridged)

The bright maser emission produced by several molecular species at centimeter to long millimeter wavelengths provides an essential tool for understanding the process of massive star formation. Unimpeded by the high dust optical depths that affect shorter wavelength observations, the high brightness temperature of these emission lines offers a way to resolve accretion and outflow motions down to scales below $\sim$1 au in deeply embedded Galactic star-forming regions at kiloparsec distances. The recent identification of extraordinary accretion outbursts in two high-mass protostars, both of which were heralded by maser flares, has rapidly impacted the traditional view of massive protostellar evolution, leading to new hydrodynamic simulations that can produce such episodic outbursts. In order to understand how these massive protostars evolve in response to such events, larger, more sensitive ground-based centimeter wavelength interferometers are needed that can simultaneously image multiple maser species in the molecular gas along with faint continuum from the central ionized gas. Fiducial observations of a large sample of massive protostars will be essential in order to pinpoint the progenitors of future accretion outbursts, and to quantify the outburst-induced changes in their protostellar photospheres and outflow and accretion structures. Knowledge gained from these studies will have broader impact on the general topic of accretion onto massive objects.

We present the results of a monitoring campaign using the KAT-7 and HartRAO 26m telescopes, of hydroxyl, methanol and water vapour masers associated with the high-mass star forming region G9.62+0.20E. Periodic flaring of the main line hydroxyl masers were found, similar to that seen in the 6.7 and 12.2 GHz methanol masers. The 1667 MHz flares are characterized by a rapid decrease in flux density which is coincident with the start of the 12.2 GHz methanol maser flare. The decrease in the OH maser flux density is followed by a slow increase till a maximum is reached after which the maser decays to its pre-flare level. A possible interpretation of the rapid decrease in the maser flux density is presented. Considering the projected separation between the periodic methanol and OH masers, we conclude that the periodic 12.2 methanol masing region is located about 1600 AU deeper into the molecular envelope compared to the location of the periodic OH masers. A single water maser flare was also detected which seems not to be associated with the same event that gives rise to the periodic methanol and OH maser flares.

IGR J17503-2636 is a hard X-ray transient discovered by INTEGRAL on 2018 August 11. This was the first ever reported X-ray emission from this source. Following the discovery, follow-up observations were carried out with Swift, Chandra, NICER, and NuSTAR. We report in this paper the analysis and results obtained from all these X-ray data. Based on the fast variability in the X-ray domain, the spectral energy distribution in the 0.5-80 keV energy range, and the reported association with a highly reddened OB supergiant at ~10 kpc, we conclude that IGR J17503-2636 is most likely a relatively faint new member of the supergiant fast X-ray transients. Spectral analysis of the NuSTAR data revealed a broad feature in addition to the typical power-law with exponential roll-over at high energy. This can be modeled either in emission or as a cyclotron scattering feature in absorption. If confirmed by future observations, this feature would indicate that IGR J17503-2636 hosts a strongly magnetized neutron star with B~2e12 G.

Centimeter continuum observations of protostellar jets have revealed the presence of knots of shocked gas where the flux density decreases with frequency. This spectrum is characteristic of nonthermal synchrotron radiation and implies the presence of both magnetic fields and relativistic electrons in protostellar jets. Here, we report on one of the few detections of nonthermal jet driven by a young massive star in the star-forming region G035.02$+$0.35. We made use of the NSF's Karl G. Jansky Very Large Array (VLA) to observe this region at C, Ku, and K bands with the A- and B-array configurations, and obtained sensitive radio continuum maps down to a rms of 10 $\mu$Jy beam$^{-1}$. These observations allow for a detailed spectral index analysis of the radio continuum emission in the region, which we interpret as a protostellar jet with a number of knots aligned with extended 4.5 $\mu$m emission. Two knots clearly emit nonthermal radiation and are found at similar distances, of approximately 10,000 au, each side of the central young star, from which they expand at velocities of hundreds km s$^{-1}$. We estimate both the mechanical force and the magnetic field associated with the radio jet, and infer a lower limit of $0.4\times10^{-4} $M$_{\odot}$ yr$^{-1}$ km s$^{-1}$ and values in the range $0.7-1.3 $mG, respectively.

We studied five XMM-Newton observations of the neutron-star binary 4U 1728$-$34 covering the hard, intermediate and soft spectral states. By jointly fitting the spectra with several reflection models, we obtained an inclination angle of 25$-$53$\deg$ and an iron abundance up to 10 times the solar. From the fits with reflection models, we found that the fluxes of the reflection and the Comptonised components vary inconsistently; since the latter is assumed to be the illuminating source, this result possibly indicates the contribution of the neutron star surface/boundary layer to the disc reflection. As the source evolved from the relatively soft to the intermediate state, the disc inner radius decreased, opposite to the prediction of the standard accretion disc model. We also explore the possible reasons why the supersolar iron abundance is required by the data and found that this high value is probably caused by the absence of the hard photons in the XMM-Newton data.

It is well established that Solar-mass stars gain mass via disk accretion, until the mass reservoir of the disk is exhausted and dispersed, or condenses into planetesimals. Accretion disks are intimately coupled with mass ejection via polar cavities, in the form of jets and less collimated winds, which allow mass accretion through the disk by removing a substantial fraction of its angular momentum. Whether disk accretion is the mechanism leading to the formation of stars with much higher masses is still unclear. Here, we are able to build a comprehensive picture for the formation of an O-type star, by directly imaging a molecular disk which rotates and undergoes infall around the central star, and drives a molecular jet which arises from the inner disk regions. The accretion disk is truncated between 2000-3000au, it has a mass of about a tenth of the central star mass, and is infalling towards the central star at a high rate (6x10^-4 Msun/yr), as to build up a very massive object. These findings, obtained with the Atacama Large Millimeter/submillimeter Array at 700au resolution, provide observational proof that young massive stars can form via disk accretion much like Solar-mass stars.

We previously observed the star-forming region G16.59$-$0.05 through interferometric observations of both thermal and maser lines, and identified a high-mass young stellar object (YSO) which is surrounded by an accretion disk and drives a nonthermal radio jet. We performed high-angular-resolution (beam FWHM ~0.15") 1.2-mm continuum and line observations towards G16.59$-$0.05 with the Atacama Large Millimeter Array (ALMA). The main dust clump, with size ~10$^4$ au, is resolved into four relatively compact (diameter ~2000 au) millimeter (mm) sources. The source harboring the high-mass YSO is the most prominent in molecular emission. By fitting the emission profiles of several unblended and optically thin transitions of CH$_3$OCH$_3$ and CH$_3$OH, we derived gas temperatures inside the mm-sources in the range 42--131 K, and calculated masses of 1--5 $M_{\odot}$. A well-defined Local Standard of Rest velocity (Vlsr) gradient is detected in most of the high-density molecular tracers at the position of the high-mass YSO, pinpointed by compact 22-GHz free-free emission. This gradient is oriented along a direction forming a large (~70 degree) angle with the radio jet, traced by elongated 13-GHz continuum emission. The butterfly-like shapes of the P-V plots and the linear pattern of the emission peaks of the molecular lines at high velocity confirm that this Vlsr gradient is due to rotation of the gas in the disk surrounding the high-mass YSO. The disk radius is ~500 au, and the Vlsr distribution along the major axis of the disk is well reproduced by a Keplerian profile around a central mass of 10$\pm$2 $M_{\odot}$. The position of the YSO is offset by >~ 0.1" from the axis of the radio jet and the dust emission peak. To explain this displacement we argue that the high-mass YSO could have moved from the center of the parental mm source owing to dynamical interaction with one or more companions.

The High Energy Modular Ensemble of Satellites (HERMES) project is aimed to realize a modular X/gamma-ray monitor for transient events, to be placed on-board of a CubeSat bus. This expandable platform will achieve a significant impact on Gamma Ray Burst (GRB) science and on the detection of Gravitational Wave (GW) electromagnetic counterparts: the recent LIGO/VIRGO discoveries demonstrated that the high-energy transient sky is still a field of extreme interest. The very complex temporal variability of GRBs (up to the millisecond scale) combined with the spatial and temporal coincidence between GWs and their electromagnetic counterparts suggest that upcoming instruments require sub-ms time resolution combined with a transient localization accuracy lower than a degree. The current phase of the ongoing HERMES project is focused on the realization of a technological pathfinder with a small network (3 units) of nano-satellites to be launched in mid 2020. We will show the potential and prospects for short and medium-term development of the project, demonstrating the disrupting possibilities for scientific investigations provided by the innovative concept of a new "modular astronomy" with nano-satellites (e.g. low developing costs, very short realization time). Finally, we will illustrate the characteristics of the HERMES Technological Pathfinder project, demonstrating how the scientific goals discussed are actually already reachable with the first nano-satellites of this constellation. The detector architecture will be described in detail, showing that the new generation of scintillators (e.g. GAGG:Ce) coupled with very performing Silicon Drift Detectors (SDD) and low noise Front-End-Electronics (FEE) are able to extend down to few keV the sensitivity band of the detector. The technical solutions for FEE, Back-End-Electronics (BEE) and Data Handling will be also described.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

Most of the X-ray binary systems containing neutron stars classified as Atoll sources show two different spectral states, called soft and hard. Moreover, a large number of these systems show a reflection component relativistically smeared in their spectra, which gives information on the innermost region of the system. Our aim is to investigate the poorly studied broadband spectrum of the low mass X-ray binary system 4U 1702-429, which was recently analysed combining XMM-Newton and INTEGRAL data. The peculiar value of the reflection fraction brought us to analyse further broadband spectra of 4U 1702-429. We re-analysed the spectrum of the XMM-Newton/INTEGRAL observation of 4U 1702-429 in the 0.3-60 keV energy range and we extracted three 0.1-100 keV spectra of the source analysing three observations collected with the BeppoSAX satellite. We find that the XMM-Newton/INTEGRAL spectrum is well fitted using a model composed of a disc blackbody plus a Comptonised component and a smeared reflection component. We used the same spectral model for the BeppoSAX spectra, finding out that the addition of a smeared reflection component is statistically significant. The best-fit values of the parameters are compatible to each other for the BeppoSAX spectra. We find that the reflection fraction is $0.05^{+0.03}_{-0.01}$ for the XMM-Newton/INTEGRAL spectrum and between 0.15 and 0.4 for BeppoSAX ones. The relative reflection fraction and the ionisation parameter are incompatible between the XMM-Newton/INTEGRAL and the BeppoSAX observations and the characteristics of the Comptonising corona suggest that the source was in a soft state in the former observation and in a hard state in the latter.

The first discovered accreting millisecond pulsar, SAX J1808.4-3658, went into X-ray outburst in April 2015. We triggered a 100 ks XMM-Newton ToO, taken at the peak of the outburst, and a 55 ks NuSTAR ToO, performed four days apart. We report here the results of a detailed spectral analysis of both the XMM-Newton and NuSTAR spectra. While the XMM-Newton spectrum appears much softer than in previous observations, the NuSTAR spectrum confirms the results obtained with XMM-Newton during the 2008 outburst. We find clear evidence of a broad iron line that we interpret as produced by reflection from the inner accretion disk. For the first time, we use a self-consistent reflection model to fit the reflection features in the NuSTAR spectrum; in this case we find a statistically significant improvement of the fit with respect to a simple Gaussian or diskline model to fit the iron line, implying that the reflection continuum is also significantly detected. Despite the differences evident between the XMM-Newton and NuSTAR spectra, the smearing best-fit parameters found for these spectra are consistent with each other and are compatible with previous results. In particular, we find an upper limit to the inner disk radius of $\sim 12~R_g$. In all the cases, a high inclination angle ($>50^\circ$) of the system is required. This inclination angle, combined with measurements of the radial velocity of the optical companion, results in a low value for the neutron star mass ($<0.8\,M_\odot$), a result that deserves further investigation.

Imaging the bright maser emission produced by several molecular species at centimeter wavelengths is an essential tool for understanding the process of massive star formation because it provides a way to probe the kinematics of dense molecular gas at high angular resolution. Unimpeded by the high dust optical depths that affect shorter wavelength observations, the high brightness temperature of these emission lines offers a way to resolve accretion and outflow motions down to scales as fine as $\sim$1-10 au in deeply embedded Galactic star-forming regions, and at sub-pc scales in nearby galaxies. The Next Generation Very Large Array will provide the capabilities needed to fully exploit these powerful tracers.

We report on the discovery by the Nuclear Spectroscopic Telescope Array (NuSTAR) and the Neutron Star Interior Composition Explorer (NICER) of the accreting millisecond X-ray pulsar IGR J17591-2342, detecting coherent X-ray pulsations around 527.4 Hz (1.9 ms) with a clear Doppler modulation. This implies an orbital period of ~8.8 hours and a projected semi-major axis of ~1.23 lt-s. From the binary mass function, we estimate a minimum companion mass of 0.42 solar masses, obtained assuming a neutron star mass of 1.4 solar masses and an inclination angle lower than 60 degrees, as suggested by the absence of eclipses or dips in the light-curve of the source. The broad-band energy spectrum is dominated by Comptonisation of soft thermal seed photons with a temperature of ~0.7 keV by electrons heated to 21 keV. We also detect black-body-like thermal direct emission compatible with an emission region of a few kilometers and temperature compatible with the seed source of Comptonisation. A weak Gaussian line centered on the iron K-alpha; complex can be interpreted as a signature of disc reflection. A similar spectrum characterises the NICER spectra, measured during the outburst fading.

We discuss the spectral and timing properties of the accreting millisecond X-ray pulsar SWIFT J1756.9-2508 observed by XMM-Newton, NICER and NuSTAR during the X-ray outburst occurred in April 2018. The spectral properties of the source are consistent with a hard state dominated at high energies by a non-thermal power-law component with a cut-off at ~70 keV. No evidence of iron emission lines or reflection humps has been found. From the coherent timing analysis of the pulse profiles, we derived an updated set of orbital ephemerides. Combining the parameters measured from the three outbursts shown by the source in the last ~11 years, we investigated the secular evolution of the spin frequency and the orbital period. We estimated a neutron magnetic field of 3.1E+8 G < B_pc< 4.5E+8 G and measured an orbital period derivative of -4.1E-12 s/s < P_dot_orb < 7.1E-12 s/s. We also studied the energy dependence of the pulse profile by characterising the behaviour of the pulse fractional amplitude in the energy range 0.3-80 keV. These results are compared with those obtained from the previous outbursts of SWIFT J1756.9-2508 and other previously known accreting millisecond X-ray pulsars.

GRO J1744-28, commonly known as the `Bursting Pulsar', is a low mass X-ray binary containing a neutron star and an evolved giant star. This system, together with the Rapid Burster (MXB 1730-33), are the only two systems that display the so-called Type II X-ray bursts. These type of bursts, which last for 10s of seconds, are thought to be caused by viscous instabilities in the disk; however the Type II bursts seen in GRO J1744-28 are qualitatively very different from those seen in the archetypal Type II bursting source the Rapid Burster. To understand these differences and to create a framework for future study, we perform a study of all X-ray observations of all 3 known outbursts of the Bursting Pulsar which contained Type II bursts, including a population study of all Type II X-ray bursts seen by RXTE. We find that the bursts from this source are best described in four distinct phenomena or `classes' and that the characteristics of the bursts evolve in a predictable way. We compare our results with what is known for the Rapid Burster and put out results in the context of models that try to explain this phenomena.

We report on the detection of X-ray pulsations at 2.1 ms from the known X-ray burster IGR J17379-3747 using XMM-Newton. The coherent signal shows a clear Doppler modulation from which we estimate an orbital period of ~1.9 hours and a projected semi-major axis of ~8 lt-ms. Taking into account the lack of eclipses (inclination angle of < 75 deg) and assuming a neutron star mass of 1.4 Msun, we estimated a minimum companion star of ~0.06 Msun. Considerations on the probability distribution of the binary inclination angle make less likely the hypothesis of a main-sequence companion star. On the other hand, the close correspondence with the orbital parameters of the accreting millisecond pulsar SAX J1808.4-3658 suggests the presence of a bloated brown dwarf. The energy spectrum of the source is well described by a soft disk black-body component (kT ~0.45 keV) plus a Comptonisation spectrum with photon index ~1.9. No sign of emission lines or reflection components is significantly detected. Finally, combining the source ephemerides estimated from the observed outbursts, we obtained a first constraint on the long-term orbital evolution of the order of dP_orb/dt = (-2.5 +/- 2.3)E-12 s/s.

Abridged. Methods. We made use of the Karl G. Jansky Very Large Array (VLA) in the B configuration at K band, and in the A configuration at both Ku and C bands, in order to image the radio continuum emission towards 25 H$_2$O maser sites with an angular resolution and thermal rms of the order of $0.1''$ and 10 $\mu$Jy beam$^{-1}$, respectively. These targets add to our pilot study of 11 maser sites presented in Moscadelli et al. (2016). The sample of H$_2$O maser sites was selected among those regions having an accurate distance measurement, obtained through maser trigonometric parallaxes, and H$_2$O maser luminosities in excess of 10$^{-6}$ L$_{\odot}$. Results. We present high-resolution radio continuum images of 33 sources belonging to 25 star-forming regions. In each region, we detect radio continuum emission within a few 1000 au of the H$_2$O masers' position; 50% of the radio continuum sources are associated with bolometric luminosities exceeding 5 $\times$ 10$^{3}$ L$_{\odot}$, including W33A and G240.32$+$0.07. We provide a detailed spectral index analysis for each radio continuum source, based on the integrated fluxes at each frequency, and produce spectral index maps with the multi-frequency-synthesis deconvolution algorithm of CASA. The radio continuum emission traces thermal bremsstrahlung in (proto)stellar winds and jets, with flux densities at 22 GHz below 3 mJy, and spectral index values between $-0.1$ and 1.3. We prove a strong correlation ($r>$0.8) between the radio continuum luminosity (L$_{\rm rad}$) and the H$_2$O maser luminosity (L$_{\rm H_2O}$) of ($\rm L_{8GHz}/mJy\times kpc^2)=10^{3.8}\times(L_{H_2O}/L_{\odot})^{0.74}$. Since H$_2$O masers are excited through shocks driven by (proto)stellar winds and jets, these results provide support to the idea that the radio continuum emission around young stars is dominated by shock-ionization.

SAX J1748.9-2021 is an accreting X-ray millisecond pulsar observed in outburst five times since its discovery in 1998. In early October 2017, the source started its sixth outburst, which lasted only ~13 days, significantly shorter than the typical 30 days duration of the previous outbursts. It reached a 0.3-70 keV unabsorbed peak luminosity of $\sim3\times10^{36}$ erg/s. This is the weakest outburst ever reported for this source to date. We analyzed almost simultaneous XMM-Newton, NuSTAR and INTEGRAL observations taken during the decaying phase of its 2017 outburst. We found that the spectral properties of SAX J1748.9-2021 are consistent with an absorbed Comptonization plus a blackbody component. The former, characterized by an electron temperature of ~20 keV, a photon index of ~1.6-1.7 keV and seed photon temperature of 0.44 keV, can be associated to a hot corona or the accretion column, while the latter is more likely originating from the neutron star surface (kT$_{bb}\sim0.6$ keV, R$_{bb}\sim2.5$ km). These findings suggest that SAX J1748.9-2021 was observed in a $hard$ spectral state, as it is typically the case for accreting millisecond pulsars in outburst.

The evolutionary status of the low mass X-ray binary SAX J1808.4-3658 is simulated by following the binary evolution of its possible progenitor system through mass transfer, starting at a period of $\sim$6.6 hr. The evolution includes angular momentum losses via magnetic braking and gravitational radiation. It also takes into account the effects of illumination of the donor by both the X-ray emission and the spin down luminosity of the pulsar. The system goes through stages of mass transfer and stages during which it is detached, where only the rotationally powered pulsar irradiates the donor. We show that the pulsar irradiation is a necessary ingredient to reach SAX J1808.4-3658 orbital period when the donor mass is reduced to 0.04 - 0.06 M$_\odot$. We also show that the models reproduce important properties of the system, including the orbital period derivative, which is shown to be directly linked to the evolution through mass transfer cycles. Moreover we find that the effects of the irradiation on the internal structure of the donor are non negligible, causing the companion star to be non completely convective at the values of mass observed for the system and significantly altering its long term evolution, as the magnetic braking remains active along the whole evolution.

We report on quasi-simultaneous observations from radio to X-ray frequencies of the neutron star X-ray binary Aql X-1 over accretion state transitions during its 2016 outburst. All the observations show radio to millimetre spectra consistent with emission from a jet, with a spectral break from optically thick to optically thin synchrotron emission that decreases from ~100 GHz to <5.5 GHz during the transition from a hard to a soft accretion state. The 5.5 GHz radio flux density as the source reaches the soft state, 0.82$\pm$0.03 mJy, is the highest recorded to date for this source. During the decay of the outburst, the jet spectral break is detected again at a frequency of ~30-100 GHz. The flux density is 0.75$\pm$0.03 mJy at 97.5 GHz at this stage. This is the first time that a change in the frequency of the jet break of a neutron star X-ray binary has been measured, indicating that the processes at play in black holes are also present in neutron stars, thus supporting the idea that the internal properties of the jet rely most critically on the conditions of the accretion disc and corona around the compact object, rather than the black hole mass or spin or the neutron star surface or magnetic field.

GRO J1744-28 (the Bursting Pulsar) is a neutron star LMXB which shows highly structured X-ray variability near the end of its X-ray outbursts. In this letter we show that this variability is analogous to that seen in Transitional Millisecond Pulsars such as PSR J1023+0038: 'missing link' systems consisting of a pulsar nearing the end of its recycling phase. As such, we show that the Bursting Pulsar may also be associated with this class of objects. We discuss the implications of this scenario; in particular, we discuss the fact that the Bursting Pulsar has a significantly higher spin period and magnetic field than any other known Transitional Pulsar. If the Bursting Pulsar is indeed transitional, then this source opens a new window of opportunity to test our understanding of these systems in an entirely unexplored physical regime.

We report on the results of the multiwavelength campaign carried out after the discovery of the INTEGRAL transient IGR J17329-2731. The optical data collected with the SOAR telescope allowed us to identify the donor star in this system as a late M giant at a distance of 2.7$^{+3.4}_{-1.2}$ kpc. The data collected quasi-simultaneously with XMM-Newton and NuSTAR showed the presence of a modulation with a period of 6680$\pm$3 s in the X-ray light curves of the source. This unveils that the compact object hosted in this system is a slowly rotating neutron star. The broadband X-ray spectrum showed the presence of a strong absorption ($\gg$10$^{23}$ cm$^{-2}$) and prominent emission lines at 6.4 keV, and 7.1 keV. These features are usually found in wind-fed systems, in which the emission lines result from the fluorescence of the X-rays from the accreting compact object on the surrounding stellar wind. The presence of a strong absorption line around $\sim$21 keV in the NuSTAR spectrum suggests a cyclotron origin, thus allowing us to estimate the neutron star magnetic field as $\sim$2.4$\times$10$^{12}$ G. All evidence thus suggests IGR J17329-2731 is a symbiotic X-ray binary. As no X-ray emission was ever observed from the location of IGR J17329-2731 by INTEGRAL (or other X-ray facilities) during the past 15 yr in orbit and considering that symbiotic X-ray binaries are known to be variable but persistent X-ray sources, we concluded that INTEGRAL caught the first detectable X-ray emission from IGR J17329-2731 when the source shined as a symbiotic X-ray binary. The Swift/XRT monitoring performed up to $\sim$3 months after the discovery of the source, showed that it maintained a relatively stable X-ray flux and spectral properties.

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.

IRAS 16293-2422 is a very well studied young stellar system seen in projection towards the L1689N cloud in the Ophiuchus complex. However, its distance is still uncertain with a range of values from 120 pc to 180 pc. Our goal is to measure the trigonometric parallax of this young star by means of H$_2$O maser emission. We use archival data from 15 epochs of VLBA observations of the 22.2 GHz water maser line. By modeling the displacement on the sky of the H$_2$O maser spots, we derived a trigonometric parallax of $7.1\pm1.3$ mas, corresponding to a distance of $141_{-21}^{+30}$ pc. This new distance is in good agreement with recent values obtained for other magnetically active young stars in the L1689 cloud. We relate the kinematics of these masers with the outflows and the recent ejections powered by source A in the system.

We report the discovery of X-ray pulsations at 105.2 Hz (9.5 ms) from the transient X-ray binary IGR J16597-3704 using NuSTAR and Swift. The source was discovered by INTEGRAL in the globular cluster NGC 6256 at a distance of 9.1 kpc. The X-ray pulsations show a clear Doppler modulation implying an orbital period of ~46 minutes and a projected semi-major axis of ~5 lt-ms, which makes IGR J16597-3704 an ultra-compact X-ray binary system. We estimated a minimum companion mass of 0.0065 solar masses, assuming a neutron star mass of 1.4 solar masses, and an inclination angle of <75 degrees (suggested by the absence of eclipses or dips in its light-curve). The broad-band energy spectrum of the source is well described by a disk blackbody component (kT ~1.4 keV) plus a comptonised power-law with photon index ~2.3 and an electron temperature of ~30 keV. Radio pulsations from the source were searched for with the Parkes observatory and not detected.

Methanol and water masers indicate young stellar objects. They often exhibit flares, and a fraction shows periodic activity. Several mechanisms might explain this behavior but the lack of concurrent infrared (IR) data complicates to identify the cause. Recently, 6.7 GHz methanol maser flares were observed, triggered by accretion bursts of high-mass YSOs which confirmed the IR-pumping of these masers. This suggests that regular IR changes might lead to maser periodicity. Hence, we scrutinized space-based IR imaging of YSOs associated with periodic methanol masers. We succeeded to extract the IR light curve from NEOWISE data for the intermediate mass YSO G107.298+5.639. Thus, for the first time a relationship between the maser and IR variability could be established. While the IR light curve shows the same period of ~34.6 days as the masers, its shape is distinct from that of the maser flares. Possible reasons for the IR periodicity are discussed.

MXB 1659-298 is a transient neutron star Low-Mass X-ray binary system that shows eclipses with a periodicity of 7.1 hr. The source went to outburst in August 2015 after 14 years of quiescence. We investigate the orbital properties of this source with a baseline of 40 years obtained combining the eight eclipse arrival times present in literature with 51 eclipse arrival times collected during the last two outbursts. A quadratic ephemeris does not fit the delays associated with the eclipse arrival times and the addition of a sinusoidal term with a period of $2.31 \pm 0.02$ yr is required. We infer a binary orbital period of $P=7.1161099(3)$ hr and an orbital period derivative of $\dot{P}=-8.5(1.2) \times 10^{-12}$ s s$^{-1}$. We show that the large orbital period derivative can be explained with a highly non conservative mass transfer scenario in which more than 98\% of the mass provided by the companion star leaves the binary system. We predict an orbital period derivative value of $\dot{P}=-6(3) \times 10^{-12}$ s s$^{-1}$ and constrain the companion star mass between $\sim$0.3 and $ 0.9 \pm 0.3$ M$_{\odot}$. Assuming that the companion star is in thermal equilibrium the periodic modulation can be due to either a gravitational quadrupole coupling due to variations of the oblateness of the companion star or with the presence of a third body of mass M$_3 >21 $ Jovian masses.

Little is known about the portion of the Milky Way lying beyond the Galactic center at distances of more than 9 kilo-parsec from the Sun. These regions are opaque at optical wavelengths due to absorption by interstellar dust, and distances are very large and hard to measure. We report a direct trigonometric parallax distance of 20.4_{-2.2}^{+2.8} kilo-parsec obtained with the Very Long Baseline Array to a water maser source in a region of active star formation. These measurements allow us to shed light on Galactic spiral structure by locating the Scutum-Centaurus spiral arm as it passes through the far side of the Milky Way, and to validate a kinematic method for determining distances in this region based on transverse motions.

We perform a reflection study on a new observation of the neutron star low-mass X-ray binary Aquila X-1 taken with NuSTAR during the August 2016 outburst and compare with the July 2014 outburst. The source was captured at $\sim32\%\ L_{\mathrm{Edd}}$, which is over four times more luminous than the previous observation during the 2014 outburst. Both observations exhibit a broadened Fe line profile. Through reflection modeling, we determine that the inner disk is truncated $R_{in,\ 2016}=11_{-1}^{+2}\ R_{g}$ (where $R_{g}=GM/c^{2}$) and $R_{in,\ 2014}=14\pm2\ R_{g}$ (errors quoted at the 90% confidence level). Fiducial neutron star parameters (M$_{NS}=1.4$ M$_{\odot}$, $R_{NS}=10$ km) give a stellar radius of $R_{NS}=4.85\ R_{g}$; our measurements rule out a disk extending to that radius at more than the $6\sigma$ level of confidence. We are able to place an upper limit on the magnetic field strength of $B\leq3.0-4.5\times10^{9}$ G at the magnetic poles, assuming that the disk is truncated at the magnetospheric radius in each case. This is consistent with previous estimates of the magnetic field strength for Aquila X-1. However, if the magnetosphere is not responsible for truncating the disk prior to the neutron star surface, we estimate a boundary layer with a maximum extent of $R_{BL,\ 2016}\sim10\ R_{g}$ and $R_{BL,\ 2014}\sim6\ R_{g}$. Additionally, we compare the magnetic field strength inferred from the Fe line profile of Aquila X-1 and other neutron star low-mass X-ray binaries to known accreting millisecond X-ray pulsars.

Context. In 1998 the first accreting millisecond pulsar, SAX J1808.4-3658, was discovered and to date 18 systems showing coherent, high frequency (> 100 Hz) pulsations in low mass X-ray binaries are known. Since their discovery, this class of sources has shown interesting and sometimes puzzling behaviours. In particular, apart from a few exceptions, they are all transient with very long X-ray quiescent periods implying a quite low averaged mass accretion rate onto the neutron star. Among these sources, XTE J0929-314 has been detected in outburst just once in about 15 years of continuous monitoring of the X-ray sky. Aims. We aim to demonstrate that a conservative mass transfer in this system will result in an X-ray luminosity that is higher than the observed, long-term averaged X-ray luminosity. Methods. Under the hypothesis of a conservative mass transfer driven by gravitational radiation, as expected for this system given the short orbital period of about 43.6 min and the low mass of the companion implied by the mass function derived from timing techniques, we calculate the expected mass transfer rate in this system and predict the long-term averaged X-ray luminosity. This is compared with the averaged, over 15 years, X-ray flux observed from the system, and a lower limit of the distance to the source is inferred. Results. This distance is shown to be > 7.4 kpc in the direction of the Galactic anticentre, implying a large height, > 1.8 kpc, of the source with respect to the Galactic plane, placing the source in an empty region of the Galaxy. We suggest that the inferred value of the distance is unlikely. (abridged)

XB 1254-690 is a dipping low mass X-ray binary system hosting a neutron star and showing type I X-ray bursts. We aim at obtaining more accurate orbital ephemeris and at constraining the orbital period derivative of the system for the first time. In addition, we want to better constrain the distance to the source in order to locate the system in a well defined evolutive scenario. We apply for the first time an orbital timing technique to XB 1254-690, using the arrival times of the dips present in the light curves that have been collected during 26 years of X-ray pointed observations performed from different space missions. We estimate the dip arrival times using a statistical method that weights the count-rate inside the dip with respect to the level of the persistent emission outside the dip. We fit the obtained delays as a function of the orbital cycles both with a linear and a quadratic function. We infer the orbital ephemeris of XB 1254-690 improving the accuracy of the orbital period with respect to previous estimates. We infer a mass of M$_{2}=0.42\pm 0.04$ M$_{\odot}$ for the donor star, in agreement with the estimations already present in literature, assuming that the star is in thermal equilibrium while it transfers part of its mass via the inner Lagrangian point, and assuming a neutron star mass of 1.4 M$_{\odot}$. Using these assumptions, we also constrain the distance to the source, finding a value of 7.6$\pm 0.8$ kpc. Finally, we discuss the evolution of the system suggesting that it is compatible with a conservative mass transfer driven by magnetic braking.

We investigate the relation between the parameters of the energy spectrum and the frequency and amplitude of the kilohertz quasi-periodic oscillations (kHz QPOs) in the low-mass X-ray binary 4U 1636$-$53. We fit the $3-180$-keV spectrum of this source with a model that includes a thermal Comptonisation component. We show that the frequencies of both kHz QPOs follow the same relation as a function of the parameters of this spectral component, except for a systematic frequency shift, whereas the rms fractional amplitude of each QPO follows a different relation with respect to those same parameters. This implies that, while the dynamical mechanism that sets the frequencies of the QPO can be the same for both kHz QPOs, the radiative mechanisms that set the amplitudes of the lower and the upper kHz QPO are likely different. We discuss the implications of these results to the modelling of the kHz QPOs and the possibility that the lower kHz QPO reflects a resonance between the Comptonising medium and the photons from the accretion disc and/or the neutron star surface.

The black-hole binary V404 Cyg entered the outburst phase in June 2015 after 26 years of X-ray quiescence, and with its behaviour broke the outburst evolution pattern typical of most black-hole binaries. We observed the entire outburst with the Swift satellite and performed time-resolved spectroscopy of its most active phase, obtaining over a thousand spectra with exposures from tens to hundreds of seconds. All the spectra can be fitted with an absorbed power law model, which most of the time required the presence of a partial covering. A blue-shifted iron-Kalpha line appears in 10% of the spectra together with the signature of high column densities, and about 20% of the spectra seem to show signatures of reflection. None of the spectra showed the unambiguous presence of soft disk-blackbody emission, while the observed bolometric flux exceeded the Eddington value in 3% of the spectra. Our results can be explained assuming that the inner part of the accretion flow is inflated into a slim disk that both hides the innermost (and brightest) regions of the flow, and produces a cold, clumpy, high-density outflow that introduces the high-absorption and fast spectral variability observed. We argue that the black hole in V404 Cyg might have been accreting erratically or even continuously at Eddington/Super-Eddington rates - thus sustaining a surrounding slim disk - while being partly or completely obscured by the inflated disk and its outflow. Hence, the largest flares produced by the source might not be accretion-driven events, but instead the effects of the unveiling of the extremely bright source hidden within the system.

We present a timing analysis of the 2015 outburst of the accreting millisecond X-ray pulsar SAX J1808.4-3658, using non-simultaneous XMM-Newton and NuStar observations. We estimate the pulsar spin frequency and update the system orbital solution. Combining the average spin frequency from the previous observed, we confirm the long-term spin down at an average rate $\dot{\nu}_{\text{SD}}=1.5(2)\times 10^{-15}$ Hz s$^{-1}$. We also discuss possible corrections to the spin down rate accounting for mass accretion onto the compact object when the system is X-ray active. Finally, combining the updated ephemerides with those of the previous outbursts, we find a long-term orbital evolution compatible with a binary expansion at a mean rate $\dot{P}_{orb}=3.6(4)\times 10^{-12}$ s s$^{-1}$, in agreement with previously reported values. This fast evolution is incompatible with an evolution driven by angular momentum losses caused by gravitational radiation under the hypothesis of conservative mass transfer. We discuss the observed orbital expansion in terms of non-conservative mass transfer and gravitational quadrupole coupling mechanism. We find that the latter can explain, under certain conditions, small fluctuations (of the order of few seconds) of the orbital period around a global parabolic trend. At the same time, a non-conservative mass transfer is required to explain the observed fast orbital evolution, which likely reflects ejection of a large fraction of mass from the inner Lagrangian point caused by the irradiation of the donor by the magneto-dipole rotator during quiescence (radio-ejection model). This strong outflow may power tidal dissipation in the companion star and be responsible of the gravitational quadrupole change oscillations.

The BeSSeL Survey is mapping the spiral structure of the Milky Way by measuring trigonometric parallaxes of hundreds of maser sources associated with high-mass star formation. While parallax techniques for water masers at high frequency (22 GHz) have been well documented, recent observations of methanol masers at lower frequency (6.7 GHz) have revealed astrometric issues associated with signal propagation through the ionosphere that could significantly limit parallax accuracy. These problems displayed as a "parallax gradient" on the sky when measured against different background quasars. We present an analysis method in which we generate position data relative to an "artificial quasar" at the target maser position at each epoch. Fitting parallax to these data can significantly mitigate the problems and improve parallax accuracy.

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

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.

We present 3-79 keV NuSTAR observations of the neutron star low-mass X-ray binary 4U 1636-53 in the soft, transitional and hard state. The spectra display a broad emission line at 5-10 keV. We applied several models to fit this line: A GAUSSIAN line, a relativistically broadened emission line model, KYRLINE, and two models including relativistically smeared and ionized reflection off the accretion disc with different coronal heights, RELXILL and RELXILLLP. All models fit the spectra well, however, the KYRLINE and RELXILL models yield an inclination of the accretion disc of $\sim88\degree$ with respect to the line of sight, which is at odds with the fact that this source shows no dips or eclipses. The RELXILLLP model, on the other hand, gives a reasonable inclination of $\sim56\degree$. We discuss our results for these models in this source and the possible primary source of the hard X-rays.

We present the spectral and timing analysis of the X-ray pulsar GRO J1744-28 during its 2014 outburst using data collected with the X-ray satellites Swift, INTEGRAL, Chandra, and XMM-Newton. We derived, by phase-connected timing analysis of the observed pulses, an updated set of the source ephemeris. We were also able to investigate the spin-up of the X-ray pulsar as a consequence of the accretion torque during the outburst. Relating the spin-up rate and the mass accretion rate as $\dot{\nu}\propto\dot{M}^{\beta}$, we fitted the pulse phase delays obtaining a value of $\beta=0.96(3)$. Combining the results from the source spin-up frequency derivative and the flux estimation, we constrained the source distance to be between 3.4-4.1 kpc, assuming a disc viscous parameter $\alpha$ to be in the range 0.1-1. Finally, we investigated the presence of a possible spin-down torque by adding a quadratic component to the pulse phase delay model. The marginal statistical improvement of the updated model does not allow us to firmly confirm the presence of this component.

We analysed the Rossi X-ray Timing Explorer (RXTE) data from a sample of bright accreting neutron star (NS) low-mass X-ray binaries (LMXBs). With the aim of studying the quasi-periodic variability as a function of the accretion regime, we carried out a systematic search of the quasi-periodic oscillations (QPOs) in the X-ray time series of these systems, using the integrated fractional variability as a tracker for the accretion states. We found that the three QPO types originally identified in the '80s for the brightest LMXBs, the so-called Z-sources, i.e., horizontal, normal and flaring branch oscillations (HBOs, NBOs and FBOs, respectively), are also identified in the slightly less bright NS LMXBs, the so-called Atoll sources, where we see QPOs with a behaviour consistent with the HBOs and FBOs. We compared the quasi-periodic variability properties of our NS sample with those of a sample of black hole (BH) LMXBs. We confirm the association between HBOs, NBOs and FBOs observed in Z-sources, with the type-C, type-B and type-A QPOs, respectively, observed in BH systems, and we extended the comparison to the HBO-like and FBO-like QPOs seen in Atoll sources. We conclude that the variability properties of BH and weakly-magnetized NS LMXBs show strong similarities, with QPOs only weakly sensitive to the nature of the central compact object in both classes of systems. We find that the historical association between kHz QPOs and high-frequency QPOs, seen around NSs and BHs, respectively, is not obvious when comparing similar accretion states in the two kinds of systems.

The progenitors of high-mass stars and clusters are still challenging to recognise. Only unbiased surveys, sensitive to compact regions of high dust column density, can unambiguously reveal such a small population of particularly massive and cold clumps. Here we study a flux limited sample of compact sources from the ATLASGAL survey to identify a sample of candidate progenitors of massive clusters in the inner Galaxy. Sensitive mid-infrared data at 21-24 $\mu$m from the WISE and MIPSGAL surveys were explored to search for embedded objects, and complementary spectroscopic data were used to investigate their stability and star formation activity. Based on such ancillary data we identify an unbiased sample of infrared-quiet massive clumps in the Galaxy that potentially represent the earliest stages of massive cluster formation. An important fraction of this sample consists of sources that have not been studied in detail before. Comparing their properties to clumps hosting more evolved embedded objects, we find that they exhibit similar physical properties in terms of mass and size, suggesting that infrared-quiet massive clumps are not only capable of forming high-mass stars, but likely also follow a single evolutionary track leading to the formation of massive clusters. The majority of the sources are not in virial-equilibrium, suggesting collapse on the clump scale. This is in line with the low number of infrared-quiet massive clumps and earlier findings that star formation, in particular for high-mass objects is a fast, dynamic process. We propose a scenario in which massive clumps start to fragment and collapse before their final mass is accumulated indicating that strong self-gravity and global collapse is needed to build up rich clusters and the most massive stars.

Context: Ser X-1 is a well studied LMXB which clearly shows a broad iron line. Recently, Miller et al. (2103) have presented broad-band, high quality NuSTAR data of SerX-1.Using relativistically smeared self-consistent reflection models, they find a value of R_in close to 1.0 R_ISCO (corresponding to 6 R_g), and a low inclination angle, less than 10 deg. Aims: The aim of this paper is to probe to what extent the choice of reflection and continuum models (and uncertainties therein) can affect the conclusions about the disk parameters inferred from the reflection component. To this aim we re-analyze all the available public NuSTAR and XMM-Newton. Ser X-1 is a well studied source, its spectrum has been observed by several instruments, and is therefore one of the best sources for this study. Methods: We use slightly different continuum and reflection models with respect to those adopted in literature for this source. In particular we fit the iron line and other reflection features with self-consistent reflection models as reflionx (with a power-law illuminating continuum modified with a high energy cutoff to mimic the shape of the incident Comptonization spectrum) and rfxconv. With these models we fit NuSTAR and XMM-Newton spectra yielding consistent spectral results. Results: Our results are in line with those already found by Miller et al. (2013) but less extreme. In particular, we find the inner disk radius at about 13 R_g and an inclination angle with respect to the line of sight of about 27 deg. We conclude that, while the choice of the reflection model has little impact on the disk parameters, as soon as a self-consistent model is used, the choice of the continuum model can be important in the precise determination of the disk parameters from the reflection component. Hence broad-band X-ray spectra are highly preferable to constrain the continuum and disk parameters.

We report on the spectral and timing properties of the accreting millisecond X-ray pulsar IGR J00291+5934 observed by XMM-Newton and NuSTAR during its 2015 outburst. The source is in a hard state dominated at high energies by a comptonization of soft photons ($\sim0.9$ keV) by an electron population with kT$_e\sim30$ keV, and at lower energies by a blackbody component with kT$\sim0.5$ keV. A moderately broad, neutral Fe emission line and four narrow absorption lines are also found. By investigating the pulse phase evolution, we derived the best-fitting orbital solution for the 2015 outburst. Comparing the updated ephemeris with those of the previous outbursts, we set a $3\sigma$ confidence level interval $-6.6\times 10^{-13}$ s/s $< \dot{P}_{orb} < 6.5 \times 10^{-13}$ s/s on the orbital period derivative. Moreover, we investigated the pulse profile dependence on energy finding a peculiar behaviour of the pulse fractional amplitude and lags as a function of energy. We performed a phase-resolved spectroscopy showing that the blackbody component tracks remarkably well the pulse-profile, indicating that this component resides at the neutron star surface (hot-spot).

We fitted the $3-180$-keV spectrum of all the observations of the neutron-star low-mass X-ray binary 4U 1636$-$53 taken with the {\it Rossi X-ray Timing Explorer} using a model that includes a thermal Comptonisation component. We found that in the low-hard state the power-law index of this component, $\Gamma$, gradually increases as the source moves in the colour-colour diagram. When the source undergoes a transition from the hard to the soft state $\Gamma$ drops abruptly; once the source is in the soft state $\Gamma$ increases again and then decreases gradually as the source spectrum softens further. The changes in $\Gamma$, together with changes of the electron temperature, reflect changes of the optical depth in the corona. The lower kilohertz quasi-periodic oscillation (kHz QPO) in this source appears only in observations during the transition from the hard to the soft state, when the optical depth of the corona is high and changes depends strongly upon the position of the source in the colour-colour diagram. Our results are consistent with a scenario in which the lower kHz QPO reflects a global mode in the system that results from the resonance between, the disc and/or the neutron-star surface, and the Comptonising corona.

In this paper, we report on the analysis of the peculiar X-ray variability displayed by the accreting millisecond X-ray pulsar IGR J00291+5934 in a 80 ks-long joint NuSTAR and XMM-Newton observation performed during the source outburst in 2015. The light curve of the source was characterized by a flaring-like behavior, with typical rise and decay time scales of ~120 s. The flares are accompanied by a remarkable spectral variability, with the X-ray emission being generally softer at the peak of the flares. A strong quasi periodic oscillation (QPO) is detected at ~8 mHz in the power spectrum of the source and clearly associated with the flaring-like behavior. This feature has the strongest power at soft X-rays (<3 keV). We carried out a dedicated hardness-ratio resolved spectral analysis and a QPO phase-resolved spectral analysis, together with an in-depth study of the source timing properties, to investigate the origin of this behavior. We suggest that the unusual variability of IGR J00291+5934 observed by XMM-Newton and NuSTAR could be produced by an heartbeat-like mechanism, similar to that operating in black-hole X-ray binaries. The possibility that this variability, and the associated QPO, are triggered by phases of quasi-stable nuclear burning, as suggested in the literature for a number of other neutron star binaries displaying a similar behavior, cannot be solidly tested in the case of IGR J00291+5934 due to the paucity of type-I X-ray bursts observed from this source.

We report on the discovery of coherent pulsations at a period of 2.9 ms from the X-ray transient MAXI J0911-655 in the globular cluster NGC 2808. We observed X-ray pulsations at a frequency of $\sim339.97$ Hz in three different observations of the source performed with XMM-Newton and NuSTAR during the source outburst. This newly discovered accreting millisecond pulsar is part of an ultra-compact binary system characterised by an orbital period of $44.3$ minutes and a projected semi-major axis of $\sim17.6$ lt-ms. Based on the mass function we estimate a minimum companion mass of 0.024 M$_{\odot}$, which assumes a neutron star mass of 1.4 M$_{\odot}$ and a maximum inclination angle of $75^{\circ}$ (derived from the lack of eclipses and dips in the light-curve of the source). We find that the companion star's Roche-Lobe could either be filled by a hot ($5\times 10^{6}$ K) pure helium white dwarf with a 0.028 M$_{\odot}$ mass (implying $i\simeq58^{\circ}$) or an old (>5 Gyr) brown dwarf with metallicity abundances between solar/sub-solar and mass ranging in the interval 0.065$-$0.085 M$_{\odot}$ (16 < $i$ < 21). During the outburst the broad-band energy spectra are well described by a superposition of a weak black-body component (kT$\sim$ 0.5 keV) and a hard cutoff power-law with photon index $\Gamma \sim$ 1.7 and cut-off at a temperature kT$_e\sim$ 130 keV. Up to the latest Swift-XRT observation performed on 2016 July 19 the source has been observed in outburst for almost 150 days, which makes MAXI J0911-655 the second accreting millisecond X-ray pulsar with outburst duration longer than 100 days.

Aims: We want to study the physical properties of the ionized jet emission in the vicinity of an O-type young stellar object (YSO), and estimate how efficient is the transfer of energy and momentum from small- to large-scale outflows. Methods: We conducted Karl G. Jansky Very Large Array (VLA) observations, at both 22 and 45 GHz, of the compact and faint radio continuum emission in the high-mass star-forming region G023.01-00.41, with an angular resolution between 0.3" and 0.1", and a thermal rms of the order of 10 uJy/beam. Results: We discovered a collimated thermal (bremsstrahlung) jet emission, with a radio luminosity (L_rad) of 24 mJy kpc^2 at 45 GHz, in the inner 1000 AU from an O-type YSO. The radio thermal jet has an opening angle of 44 degrees and brings a momentum rate of 8 10^-3 M_sun yr^-1 km/s. By combining the new data with previous observations of the molecular outflow and water maser shocks, we can trace the outflow emission from its driving source through the molecular clump, across more than two order of magnitude in length (500 AU-0.2 pc). We find that the momentum-transfer efficiency, between the inner jet emission and the extended outflow of entrained ambient gas, is near unity. This result suggests that the large-scale flow is swept-up by the mechanical force of the radio jet emission, which originates in the inner 1000 AU from the high-mass YSO.

The source 4U 1702-429 (Ara X-1) is a low-mass X-ray binary system hosting a neutron star. Albeit the source is quite bright ( $\sim10^{37}$ erg s$^{-1}$) its broadband spectrum has never been studied. Neither dips nor eclipses have been observed in the light curve suggesting that its inclination angle is smaller than 60$^{\circ}$.We analysed the broadband spectrum of 4U 1702-429 in the 0.3-60 keV energy range, using XMM-Newton and INTEGRAL data, to constrain its Compton reflection component if it is present. After excluding the three time intervals in which three type-I X-ray bursts occurred, we fitted the joint XMM-Newton and INTEGRAL spectra obtained from simultaneous observations. A broad emission line at 6.7 keV and two absorption edges at 0.87 and 8.82 keV were detected. We found that a self-consistent reflection model fits the 0.3-60 keV spectrum well. The broadband continuum is composed of an emission component originating from the inner region of the accretion disc, a Comptonised direct emission coming from a corona with an electron temperature of $2.63 \pm 0.06$ keV and an optical depth $\tau=13.6 \pm 0.2$, and, finally, a reflection component. The best-fit indicates that the broad emission line and the absorption edge at 8.82 keV, both associated with the presence of \ion{Fe}{xxv} ions, are produced by reflection in the region above the disc with a ionisation parameter of $Log(\xi) \simeq 2.7$. We have inferred that the inner radius, where the broad emission line originates, is $64^{+52}_{-15}$ km, and the inner radius of the accretion disc is $39^{+6}_{-8}$ km. (Abridged)

We report on the discovery and energy dependence of hard phase lags in the 2.14 Hz pulsed profiles of GRO J1744-28. We used data from XMM-Newton and NuSTAR. We were able to well constrain the lag spectrum with respect to the softest (0.3--2.3 keV) band: the delay shows increasing lag values reaching a maximum delay of $\sim$ 12 ms, between 6 and 6.4 keV. After this maximum, the value of the hard lag drops to 7 ms, followed by a recovery to a plateau at 9 ms for energies above 8 keV. NuSTAR data confirm this trend up to 30 keV, but the measurements are statistically poorer, and therefore, less constraining. The lag-energy pattern up to the discontinuity is well described by a logarithmic function. Assuming this is due to a Compton reverberation mechanism, we derive a size for the Compton cloud $R_{\rm{cc}}$ $\sim$ 120 $R_{\rm g}$, consistent with previous estimates on the magnetospheric radius. In this scenario, the sharp discontinuity at $\sim$ 6.5 keV appears difficult to interpret and suggests the possible influence of the reflected component in this energy range. We therefore propose the possible coexistence of both Compton and disk reverberation to explain the scale of the lags and its energy dependence.

We report full polarimetric VLBA observations of water masers towards the Turner-Welch Object in the W3(OH) high-mass star forming complex. This object drives a synchrotron jet, which is quite exceptional for a high-mass protostar, and is associated with a strongly polarized water maser source, W3(H$_2$O), making it an optimal target to investigate the role of magnetic fields on the innermost scales of protostellar disk-jet systems. The linearly polarized emission from water masers provides clues on the orientation of the local magnetic field, while the measurement of the Zeeman splitting from circular polarization provides its strength. The water masers trace a bipolar, biconical outflow at the center of the synchrotron jet. Although on scales of a few thousand AU the magnetic field inferred from the masers is on average orientated along the flow axis, on smaller scales (10s to 100s of AU), we have revealed a misalignment between the magnetic field and the velocity vectors, which arises from the compression of the field component along the shock front. Our measurements support a scenario where the magnetic field would evolve from having a dominant component parallel to the outflow velocity in the pre-shock gas, with field strengths of the order of a few tens of mG (at densities of $10^7$ cm$^{-3}$), to being mainly dominated by the perpendicular component of order of a few hundred of mG in the post-shock gas where the water masers are excited (at densities of $10^9$ cm$^{-3}$). The general implication is that in the undisturbed (i.e. not-shocked) circumstellar gas, the flow velocities would follow closely the magnetic field lines, while in the gas shocked by the prostostellar jet the magnetic field would be re-configured to be parallel to the shock front.

We present the results of VLA, ATCA, and Swift XRT observations of the 2015 outburst of the transient neutron star X-ray binary (NSXB), EXO 1745$-$248, located in the globular cluster Terzan 5. Combining (near-) simultaneous radio and X-ray measurements we measure a correlation between the radio and X-ray luminosities of $L_R\propto L_X^\beta$ with $\beta=1.68^{+0.10}_{-0.09}$, linking the accretion flow (probed by X-ray luminosity) and the compact jet (probed by radio luminosity). While such a relationship has been studied in multiple black hole X-ray binaries (BHXBs), this work marks only the third NSXB with such a measurement. Constraints on this relationship in NSXBs are strongly needed, as comparing this correlation between different classes of XB systems is key in understanding the properties that affect the jet production process in accreting objects. Our best fit disc-jet coupling index for EXO 1745$-$248 is consistent with the measured correlation in NSXB 4U 1728$-$34 ($\beta=1.5\pm 0.2$) but inconsistent with the correlation we fit using the most recent measurements from the literature of NSXB Aql X-1 ($\beta=0.76^{+0.14}_{-0.15}$). While a similar disc-jet coupling index appears to hold across multiple BHXBs in the hard accretion state, this does not appear to be the case with the three NSXBs measured so far. Additionally, the normalization of the EXO 1745$-$248 correlation is lower than the other two NSXBs, making it one of the most radio faint XBs ever detected in the hard state. We also report the detection of a type-I X-ray burst during this outburst, where the decay timescale is consistent with hydrogen burning.

While there are many dynamical mechanisms and models that try to explain the origin and phenomenology of the quasi-periodic oscillations (QPOs) seen in the X-ray light curves of low-mass X-ray binaries, few of them address how the radiative processes occurring in these extreme environments give rise to the rich set of variability features actually observed in these light curves. A step towards this end comes from the study of the energy and frequency dependence of the phase lags of these QPOs. Here we used a methodology that allowed us to study, for the first time, the dependence of the phase lags of all QPOs in the range of 1 Hz to 1300 Hz detected in the low-mass X-ray binary 4U 1636-53 upon energy and frequency as the source changes its states as it moves through the colour-colour diagram. Our results suggest that within the context of models of up-scattering Comptonization, the phase lags dependencies upon frequency and energy can be used to extract size scales and physical conditions of the medium that produces the lags.

We report on the timing analysis of the 2015 outburst of the intermittent accreting millisecond X-ray pulsar SAX J1748.9-2021 observed on March 4 by the X-ray satellite XMM-Newton. By phase-connecting the time of arrivals of the observed pulses, we derived the best-fit orbital solution for the 2015 outburst. We investigated the energy pulse profile dependence finding that the pulse fractional amplitude increases with energy while no significant time lags are detected. Moreover, we investigated the previous outbursts from this source, finding previously undetected pulsations in some intervals during the 2010 outburst of the source. Comparing the updated set of orbital parameters, in particular the value of the time of passage from the ascending node, with the orbital solutions reported from the previous outbursts, we estimated for the first time the orbital period derivative corresponding with $\dot{P}_{orb}=(1.1\pm0.3)\times 10^{-10}$ s/s. We note that this value is significant at 3.5 sigma confidence level, because of significant fluctuations with respect to the parabolic trend and more observations are needed in order to confirm the finding. Assuming the reliability of the result, we suggest that the large value of the orbital-period derivative can be explained as a result of an highly non-conservative mass transfer driven by emission of gravitational waves, which implies the ejection of matter from a region close to the inner Lagrangian point. We also discuss possible alternative explanations.

To probe the circumstellar environment of IRAS 13481-6124, a 20 M_sun high-mass young stellar object (HMYSO) with a parsec-scale jet and accretion disc, we investigate the origin of its Br\gamma-emission line through NIR interferometry. We present the first AMBER/VLTI observations of the Br\gamma-emitting region in an HMYSO at R~1500. Our AMBER/VLTI observations reveal a spatially and spectrally resolved Br\gamma-line in emission with a strong P Cygni profile, indicating outflowing matter with a terminal velocity of ~500 km/s. Visibilities, differential phases, and closure phases are detected in our observations within the spectral line and in the adjacent continuum. Both total visibilities (continuum plus line emitting region) and pure-line visibilities indicate that the Br\gamma-emitting region is more compact (2-4 mas in diameter or ~6-13 au at 3.2 kpc) than the continuum-emitting region (~5.4 mas or ~17 au). The absorption feature is also spatially resolved at the longest baselines (81 and 85 m) and has a visibility that is slightly smaller than the continuum-emitting region. The differential phases at the four longest baselines display an \u2018S\u2019-shaped structure across the line, peaking in the blue- and red-shifted high-velocity components. The calibrated photocentre shifts are aligned with the known jet axis, i.e they are probably tracing an ionised jet. The high-velocity components (v_r~100-500 km/s) are located far from the source, whereas the low-velocity components (0-100 km/s) are observed to be closer, indicating a strong acceleration of the gas flow in the inner 10 au. Finally, a non-zero closure phase along the continuum is detected. By comparing our observations with the synthetic images of the continuum around 2.16 um, we confirm that this feature originates from the asymmetric brightness distribution of the continuum owing to the inclination of the inner disc.

4U 1323-619 is a low mass X-ray binary system that shows type I X-ray bursts and dips. The most accurate estimation of the orbital period is 2.941923(36) hrs and a distance from the source that is lower than 11 kpc has been proposed. We aim to obtain the orbital ephemeris, the orbital period of the system, as well as its derivative to compare the observed luminosity with that predicted by the theory of secular evolution. We took the advantage of about 26 years of X-ray data and grouped the selected observations when close in time. We folded the light curves and used the timing technique, obtaining 12 dip arrival times. We fit the delays of the dip arrival times both with a linear and a quadratic function. We locate 4U 1323-619 within a circular area centred at RA (J2000)= 201.6543\degree and DEC (J2000)= -62.1358\degree with an associated error of 0.0002\degree, and confirm the detection of the IR counterpart already discussed in literature. We estimate an orbital period of P=2.9419156(6) hrs compatible with the estimations that are present in the literature, but with an accuracy ten times higher. We also obtain a constraint on the orbital period derivative for the first time, estimating $\dot{P}=(8\pm 13)\times 10^{-12}$ s/s. Assuming that the companion star is in thermal equilibrium in the lower main sequence, and is a neutron star of 1.4 M$_{\odot}$, we infer a mass of 0.28$\pm$0.03 M$_{\odot}$ for the companion star. Assuming a distance of 10 kpc, we obtained a luminosity of (4.3$\pm$0.5)$\times 10^{36}$ erg s$^{-1}$, which is not in agreement with what is predicted by the theory of secular evolution. Using a 3D extinction map of the K$_{s}$ radiation in our Galaxy, we obtain a distance of 4.2$^{+0.8}_{-0.7}$ kpc at 68\% confidence level. (Abridged)

We analyzed a 115 ks XMM-Newton observation and the stacking of 8 days of INTEGRAL observations, taken during the raise of the 2015 outburst of the accreting millisecond X-ray pulsar SAX J1748.9-2021. The source showed numerous type-I burst episodes during the XMM-Newton observation, and for this reason we studied separately the persistent and burst epochs. We described the persistent emission with a combination of two soft thermal components, a cold thermal Comptonization component (~2 keV) and an additional hard X-ray emission described by a power-law (photon index ~2.3). The continuum components can be associated with an accretion disc, the neutron star (NS) surface and a thermal Comptonization emission coming out of an optically thick plasma region, while the origin of the high energy tail is still under debate. In addition, a number of broad (~0.1-0.4 keV) emission features likely associated to reflection processes have been observed in the XMM-Newton data. The estimated 1.0-50 keV unabsorbed luminosity of the source is ~5x10^37 erg/s, about 25% of the Eddington limit assuming a 1.4 solar mass NS. We suggest that the spectral properties of SAX J1748.9-2021 are consistent with a soft state, differently from many other accreting X-ray millisecond pulsars which are usually found in the hard state. Moreover, none of the observed type-I burst reached the Eddington luminosity. Assuming that the burst ignition and emission are produced above the whole NS surface, we estimate a neutron star radius of ~7-8 km, consistent with previous results.

We report the detection of a possible gamma-ray counterpart of the accreting millisecond pulsar SAX J1808.4-3658. The analysis of ~6 years of data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope (Fermi-LAT) within a region of 15deg radius around the position of the pulsar reveals a point gamma-ray source detected at a significance of ~6 sigma (Test Statistic TS = 32), with position compatible with that of SAX J1808.4-3658 within 95% Confidence Level. The energy flux in the energy range between 0.6 GeV and 10 GeV amounts to (2.1 +- 0.5) x 10-12 erg cm-2 s-1 and the spectrum is well-represented by a power-law function with photon index 2.1 +- 0.1. We searched for significant variation of the flux at the spin frequency of the pulsar and for orbital modulation, taking into account the trials due to the uncertainties in the position, the orbital motion of the pulsar and the intrinsic evolution of the pulsar spin. No significant deviation from a constant flux at any time scale was found, preventing a firm identification via time variability. Nonetheless, the association of the LAT source as the gamma-ray counterpart of SAX J1808.4-3658 would match the emission expected from the millisecond pulsar, if it switches on as a rotation-powered source during X-ray quiescence.

We want to study the velocity and magnetic field morphology in the vicinity (<1000 AU) of a massive young stellar object (YSO), at very high spatial resolution (10-100 AU). We performed milli-arcsecond polarimetric observations of the strong CH3OH maser emission observed in the vicinity of an O-type YSO, in G023.01-00.41. We have combined this information with the velocity field of the CH3OH masing gas previously measured at the same angular resolution. We analyse the velocity and magnetic fields in the reference system defined by the direction of the molecular outflow and the equatorial plane of the hot molecular core at its base, as recently observed on sub-arcsecond scales. We provide a first detailed picture of the gas dynamics and magnetic field configuration within a radius of 2000 AU from a massive YSO. We have been able to reproduce the magnetic field lines for the outer regions (>600 AU) of the molecular envelope, where the magnetic field orientation shows a smooth change with the maser cloudlets position (0.2 degree/AU). Overall, the velocity field vectors well accommodate with the local, magnetic field direction, but still show an average misalignment of 30 degrees. We interpret this finding as the contribution of a turbulent velocity field of about 3.5 km/s, responsible for braking up the alignment between the velocity and magnetic field vectors. We do resolve different gas flows which develop both along the outflow axis and across the disk plane, with an average speed of 7 km/s. In the direction of the outflow axis, we establish a collimation of the gas flow, at a distance of about 1000 AU from the disk plane. In the disk region, gas appears to stream outward along the disk plane for radii greater than 500-600 AU, and inward for shorter radii.

The ultra-compact dipping source \object{XB 1916-053} has an orbital period of close to 50 min and a companion star with a very low mass (less than 0.1 M$_{\odot}$). The orbital period derivative of the source was estimated to be $1.5(3) \times 10^{-11}$ s/s through analysing the delays associated with the dip arrival times obtained from observations spanning 25 years, from 1978 to 2002. The known orbital period derivative is extremely large and can be explained by invoking an extreme, non-conservative mass transfer rate that is not easily justifiable. We extended the analysed data from 1978 to 2014, by spanning 37 years, to verify whether a larger sample of data can be fitted with a quadratic term or a different scenario has to be considered. We obtained 27 delays associated with the dip arrival times from data covering 37 years and used different models to fit the time delays with respect to a constant period model.We find that the quadratic form alone does not fit the data. The data are well fitted using a sinusoidal term plus a quadratic function or, alternatively, with a series of sinusoidal terms that can be associated with a modulation of the dip arrival times due to the presence of a third body that has an elliptical orbit. We infer that for a conservative mass transfer scenario the modulation of the delays can be explained by invoking the presence of a third body with mass between 0.10-0.14 M$_{\odot}$, orbital period around the X-ray binary system of close to 51 yr and an eccentricity of $0.28 \pm 0.15$. In a non-conservative mass transfer scenario we estimate that the fraction of matter yielded by the degenerate companion star and accreted onto the neutron star is $\beta = 0.08$, the neutron star mass is $\ge 2.2$ M$_{\odot}$, and the companion star mass is 0.028 M$_{\odot}$. (Abridged)

We present an analysis of the energy and frequency dependence of the Fourier time lags of the hectoHertz quasi-periodic oscillations (QPOs) and of the QPOs at the frequency at which the power density spectrum shows a break in the neutron-star low-mass X-ray binary 4U 1636-53, using a large data set obtained with the Rossi X-ray Timing Explorer. We found that: (i) For the break frequency QPO: for low frequencies, in general the time lag is positive, but it is decreasing with increasing frequency, reaching zero lag at 20 Hz. Between 20 and 35 Hz there is a small fluctuation around zero, from where the time lags become positive again and increase slightly above zero up to 65 Hz. (ii) For the hHz QPO: we see that when the frequency is 100 Hz the time lag is negative, but it increases to zero already at 110 Hz, being consistent with this value up to 130 Hz from where it increases to 0.5 msec at around 140 Hz. From 140 Hz the time lag decreases sharply, being strongly negative for hHz greater than 220 Hz. We compare these lags with our previous results for the lags of the kiloHertz QPOs in this same source and discuss possible scenarios for producing the lags in this system in the context of reflection off the accretion disc or (up-)down-scattering in a hot medium close to the neutron star.

The high-mass star-forming site G009.62+00.20E hosts the 6.7 GHz methanol maser source with the greatest flux density in the Galaxy which has been flaring periodically over the last ten years. We performed high-resolution astrometric measurements of the CH3OH, H2O, and OH maser emission and 7 mm continuum in the region. The radio continuum emission was resolved in two sources separated by 1300 AU. The CH3OH maser cloudlets are distributed along two north-south ridges of emission to the east and west of the strongest radio continuum component. This component likely pinpoints a massive young stellar object which heats up its dusty envelope, providing a constant IR pumping for the Class II CH3OH maser transitions. We suggest that the periodic maser activity may be accounted for by an independent, pulsating, IR radiation field provided by a bloated protostar in the vicinity of the brightest masers. We also report about the discovery of an elliptical distribution of CH3OH maser emission in the region of periodic variability.

Broad emission features of abundant chemical elements, such as Iron, are commonly seen in the X-ray spectra of accreting compact objects and their studies can provide useful information about the geometry of the accretion processes. In this work, we focus our attention on GX 3+1, a bright, persistent accreting low mass X-ray binary, classified as an atoll source. Its spectrum is well described by an accretion disc plus a stable comptonizing, optically thick corona which dominates the X-ray emission in the 0.3-20 keV energy band. In addition, four broad emission lines are found and we associate them with reflection of hard photons from the inner regions of the accretion disc where doppler and relativistic effects are important. We used self-consistent reflection models to fit the spectra of the 2010 XMM-Newton observation and the stacking of the whole datasets of 2010 INTEGRAL observations. We conclude that the spectra are consistent with reflection produced at ~10 gravitational radii by an accretion disc with an ionization parameter of xi~600 erg cm/s and viewed under an inclination angle of the system of ~35{\deg}. Furthermore, we detected for the first time for GX 3+1, the presence of a powerlaw component dominant at energies higher than 20 keV, possibly associated with an optically thin component of non-thermal electrons.

We have conducted the first parallax and proper motion measurements of 6.7 GHz methanol maser emission using the Australian Long Baseline Array (LBA). The parallax of G339.884$-$1.259 measured from five epochs of observations is 0.48$\pm $0.08 mas, corresponding to a distance of $2.1^{+0.4}_{-0.3}$ kpc, placing it in the Scutum spiral arm. This is consistent (within the combined uncertainty) with the kinematic distance estimate for this source at 2.5$\pm $0.5 kpc using the latest Solar and Galactic rotation parameters. We find from the Lyman continuum photon flux that the embedded core of the young star is of spectral type B1, demonstrating that luminous 6.7 GHz methanol masers can be associated with high-mass stars towards the lower end of the mass range.

We analyse all available X-ray observations of X1822-371 made with XMM-Newton, Chandra, Suzaku and INTEGRAL satellites. The observations were not simultaneous. The Suzaku and INTEGRAL broad band energy coverage allows us to constrain the spectral shape of the continuum emission well. We use the model already proposed for this source, consisting of a Comptonised component absorbed by interstellar matter and partially absorbed by local neutral matter, and we added a Gaussian feature in absorption at $\sim 0.7$ keV. This addition significantly improves the fit and flattens the residuals between 0.6 and 0.8 keV. We interpret the Gaussian feature in absorption as a cyclotron resonant scattering feature (CRSF) produced close to the neutron star surface and derive the magnetic field strength at the surface of the neutron star, $(8.8 \pm 0.3) \times 10^{10}$ G for a radius of 10 km. We derive the pulse period in the EPIC-pn data to be 0.5928850(6) s and estimate that the spin period derivative of X1822-371 is $(-2.55 \pm 0.03) \times 10^{-12}$ s/s using all available pulse period measurements. Assuming that the intrinsic luminosity of X1822-371is at the Eddington limit and using the values of spin period and spin period derivative of the source, we constrain the neutron star and companion star masses. We find the neutron star and the companion star masses to be $1.69 \pm 0.13$ M$_{\odot}$ and $0.46 \pm 0.02$ M$_{\odot}$, respectively, for a neutron star radius of 10 km.In a self-consistent scenario in which X1822-371 is spinning-up and accretes at the Eddington limit, we estimate that the magnetic field of the neutron star is $(8.8 \pm 0.3) \times 10^{10}$ G for a neutron star radius of 10 km. If our interpretation is correct, the Gaussian absorption feature near 0.7 keV is the very first detection of a CRSF below 1 keV in a LMXB. (abridged)

The bursting pulsar, GRO J1744-28, went again in outburst after $\sim$18 years of quiescence in mid-January 2014. We studied the broad-band, persistent, X-ray spectrum using X-ray data from a XMM-Newton observation, performed almost at the peak of the outburst, and from a close INTEGRAL observation, performed 3 days later, thus covering the 1.3-70.0 keV band. The spectrum shows a complex continuum shape that cannot be modelled with standard high-mass X-ray pulsar models, nor by two-components models. We observe broadband and peaked residuals from 4 to 15 keV, and we propose a self-consistent interpretation of these residuals, assuming they are produced by cyclotron absorption features and by a moderately smeared, highly ionized, reflection component. We identify the cyclotron fundamental at $\sim$ 4.7 keV, with hints for two possible harmonics at 10.4 keV and 15.8 keV. The position of the cyclotron fundamental allows an estimate for the pulsar magnetic field of (5.27 $\pm$ 0.06) $\times$ 10$^{11}$ G, if the feature is produced at its surface. From the dynamical and relativistic smearing of the disk reflected component, we obtain a lower limit estimate for the truncated accretion disk inner radius, ($\gtrsim$ 100 R$_g$), and for the inclination angle (18$^{\circ}$-48$^{\circ}$). We also detect the presence of a softer thermal component, that we associate with the emission from an accretion disk truncated at a distance from the pulsar of 50-115 R$_g$. From these estimates, we derive the magneto-spheric radius for disk accretion to be $\sim$ 0.2 times the classical Alfv\'en radius for radial accretion.

We report parallaxes and proper motions of three water maser sources in high-mass star-forming regions in the Outer Spiral Arm of the Milky Way. The observations were conducted with the Very Long Baseline Array as part of Bar and Spiral Structure Legacy Survey and double the number of such measurements in the literature. The Outer Arm has a pitch angle of 14.9 +/- 2.7 deg and a Galactocentric distance of 14.1 +/- 0.6 kpc toward the Galactic anticenter. The average motion of these sources toward the Galactic center is 10.7 +/- 2.1 km/s and we see no sign of a significant fall in the rotation curve out to 15 kpc from the Galactic center. The three-dimensional locations of these star-forming regions are consistent with a Galactic warp of several hundred parsecs from the plane.

We present a systematic study of the orbital inclination effects on black-hole transients fast time-variability properties. We have considered all the black-hole binaries that have been densely monitored by the Rossi XTE satellite. We find that the amplitude of low-frequency quasi periodic oscillations (QPOs) depends on the orbital inclination. Type-C QPOs are stronger for nearly edge-on systems (high inclination), while type-B QPOs are stronger when the accretion disk is closer to face-on (low inclination). Our results also suggest that the noise associated with type-C QPOs is consistent with being stronger for low-inclination sources, while the noise associated to type-B QPOs seems inclination independent. These results are consistent with a geometric origin of the type-C QPOs - for instance arising from relativistic precession of the inner flow within a truncated disk - while the noise would correspond to intrinsic brightness variability from mass accretion rate fluctuations in the accretion flow. The opposite behavior of type-B QPOs - stronger in low inclinations sources - supports the hypothesis that type-B QPOs are related to the jet, the power of which is the most obvious measurable parameter expected to be stronger in nearly face-on sources.

We aim at deriving the main physical properties of massive jets from near-IR observations, comparing them to those of a large sample of jets from low-mass YSOs, and relating them to the main features of their driving sources. We present a NIR imaging (H2 and Ks) and low-resolution spectroscopic (0.95-2.50 um) survey of 18 massive jets towards GLIMPSE extended green objects, driven by intermediate- and high-mass YSOs, which have Lbol between 4x10^2 and 10^5 Lsun. As in low-mass jets, H2 is the primary NIR coolant, detected in all the analysed flows, whereas the most important ionic tracer is [FeII], detected in half of the sampled jets. Our analysis indicates that the emission lines originate from shocks at high temperatures and densities. No fluorescent emission is detected along the flows, regardless of the source Lbol. On average, the physical parameters of these massive jets (i.e. Av, temperature, column density, mass, and luminosity) have higher values than those measured in their low-mass counterparts. The morphology of the H2 flows is varied, mostly depending on the complex, dynamic, and inhomogeneous environment in which these massive jets form and propagate. All flows and jets in our sample are collimated, showing large precession angles. Additionally, the presence of both knots and jets suggests that the ejection process is continuous with burst episodes, as in low-mass YSOs. We compare the flow H2 luminosity with the source Lbol confirming the tight correlation between these two quantities. Five sources, however, display a lower L(H2)/Lbol efficiency, which might be related to YSO evolution. Most important, the inferred L(H2) vs Lbol relationship agrees well with the correlation between the momentum flux of the CO outflows and the bolometric luminosities of high-mass YSOs indicating that outflows from high-mass YSOs are momentum driven, as are their low-mass counterparts.

When the EPIC-pn instrument on board XMM-Newton is operated in Timing mode, high count rates (>100 cts/s) of bright sources may affect the calibration of the energy scale, resulting in a modification of the real spectral shape. The corrections related to this effect are then strongly important in the study of the spectral properties. Tests of these calibrations are more suitable in sources which spectra are characterised by a large number of discrete features. Therefore, in this work, we carried out a spectral analysis of the accreting Neutron Star GX 13+1, which is a dipping source with several narrow absorption lines and a broad emission line in its spectrum. We tested two different correction approaches on an XMM-Newton EPIC-pn observation taken in Timing mode: the standard Rate Dependent CTI (RDCTI or epfast) and the new, Rate Dependent Pulse Height Amplitude (RDPHA) corrections. We found that, in general, the two corrections marginally affect the properties of the overall broadband continuum, while hints of differences in the broad emission line spectral shape are seen. On the other hand, they are dramatically important for the centroid energy of the absorption lines. In particular, the RDPHA corrections provide a better estimate of the spectral properties of these features than the RDCTI corrections. Indeed the discrete features observed in the data, applying the former method, are physically more consistent with those already found in other Chandra and XMM-Newton observations of GX 13+1.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

We report measurements of trigonometric parallaxes for six high-mass star-forming regions in the Scutum spiral arm of the Milky Way as part of the BeSSeL Survey. Combining our measurements with 10 previous measurements from the BeSSeL Survey yields a total sample of 16 sources in the Scutum arm with trigonometric parallaxes in the Galactic longitude range from 5 deg to 32 deg. Assuming a logarithmic spiral model, we estimate a pitch angle of 19.8 (+/-3.1) deg for the Scutum arm, which is larger than pitch angles reported for other spiral arms. The high pitch angle of the arm may be due to the arm's proximity to the Galactic bar. The Scutum arm sources show an average peculiar motion of 4 km/s slower than the Galactic rotation and 8 km/s toward the Galactic center. While the direction of this non-circular motion has the same sign as determined for sources in other spiral arms, the motion toward the Galactic center is greater for the Scutum arm sources.

Both the broad iron (Fe) line and the frequency of the kilohertz quasi-periodic oscillations (kHz QPOs) in neutron star low-mass X-ray binaries (LMXBs) can potentially provide independent measures of the inner radius of the accretion disc. We use XMM-Newton and simultaneous Rossi X-ray Timing Explorer observations of the LMXB 4U 1636-53 to test this hypothesis. We study the properties of the Fe-K emission line as a function of the spectral state of the source and the frequency of the kHz QPOs. We find that the inner radius of the accretion disc deduced from the frequency of the upper kHz QPO varies as a function of the position of the source in the colour-colour diagram, in accordance with previous work and with the standard scenario of accretion disc geometry. On the contrary, the inner disc radius deduced from the profile of the Fe line is not correlated with the spectral state of the source. The values of the inner radius inferred from kHz QPOs and Fe lines, in four observations, do not lead to a consistent value of the neutron star mass, regardless of the model used to fit the Fe line. Our results suggest that either the kHz QPO or the standard relativistic Fe line interpretation does not apply for this system. Furthermore, the simultaneous detection of kHz QPOs and broad Fe lines is difficult to reconcile with models in which the broadening of the Fe line is due to the reprocessing of photons in an outflowing wind.

(Abridged) METHODS: We performed SMA observations at 1.3 mm with both the most extended and compact array configurations, providing sub-arcsecond and high sensitivity maps of various molecular lines, including both hot-core and outflow tracers. We also reconstruct the spectral energy distribution of the region from millimeter to near infrared wavelengths, using the Herschel/Hi-GAL maps, as well as archival data. RESULTS: From the spectral energy distribution, we derive a bolometric luminosity of about 4x10^4 Lsun. Our interferometric observations reveal that the distribution of dense gas and dust in the HMC is significantly flattened and extends up to a radius of 8000 AU from the center of radio continuum and maser emission in the region. The equatorial plane of this HMC is strictly perpendicular to the elongation of the collimated bipolar outflow, as imaged on scales of about 0.1-0.5 pc in the main CO isotopomers as well as in the SiO(5-4) line. In the innermost HMC regions (ca. 1000 AU), the velocity field traced by the CH3CN(12_K-11_K) line emission shows that molecular gas is both expanding along the outflow direction following a Hubble-law, and rotating about the outflow axis, in agreement with the (3-D) velocity field traced by methanol masers. The velocity field associated with rotation indicates a dynamical mass of 19 Msun at the center of the core. The latter is likely to be concentrated in a single O9.5 ZAMS star, consistent with the estimated bolometric luminosity of G023.01-00.41. The physical properties of the CO(2-1) outflow emission, such as its momentum rate 6x10^-3 Msun km/s /yr and its outflow rate 2x10^-4 Msun/yr, support our estimates of the luminosity (and mass) of the embedded young stellar object.

We used six simultaneous XMM-Newton and Rossi X-ray Timing Explorer plus five Suzaku observations to study the continuum spectrum and the iron emission line in the neutron-star low-mass X-ray binary 4U 1636-53. We modelled the spectra with two thermal components (representing the accretion disc and boundary layer), a Comptonised component (representing a hot corona), and either a Gaussian or a relativistic line component to model an iron emission line at about 6.5 keV. For the relativistic line component we used either the diskline, laor or kyrline model, the latter for three different values of the spin parameter. The fitting results for the continuum are consistent with the standard truncated disc scenario. We also find that the flux and equivalent width of the iron line first increase and then decrease as the flux of the Comptonised component increases. This could be explained either by changes in the ionisation state of the accretion disc where the line is produced by reflection, or by light bending of the emission from the Comptonised component if the height at which this component is produced changes with mass accretion rate.

Over 100 trigonometric parallaxes and proper motions for masers associated with young, high-mass stars have been measured with the BeSSeL Survey, a VLBA key science project, the EVN, and the Japanese VERA project. These measurements provide strong evidence for the existence of spiral arms in the Milky Way, accurately locating many arm segments and yielding spiral pitch angles ranging from 7 to 20 degrees. The widths of spiral arms increase with distance from the Galactic center. Fitting axially symmetric models of the Milky Way with the 3-D position and velocity information and conservative priors for the solar and average source peculiar motions, we estimate the distance to the Galactic center, Ro, to be 8.34 +/- 0.16 kpc, a circular rotation speed at the Sun, To, to be 240 +/- 8 km/s, and a rotation curve that is nearly flat (a slope of -0.2 +/- 0.4 km/s/kpc) between Galactocentric radii of 5 and 16 kpc. Assuming a "universal" spiral galaxy form for the rotation curve, we estimate the thin disk scale length to be 2.44 +/- 0.16 kpc. The parameters Ro and To are not highly correlated and are relatively insensitive to different forms of the rotation curve. Adopting a theoretically motivated prior that high-mass star forming regions are in nearly circular Galactic orbits, we estimate a global solar motion component in the direction of Galactic rotation, Vsun = 14.6 +/- 5.0 km/s. While To and Vsun are significantly correlated, the sum of these parameters is well constrained, To + Vsun = 255.2 +/- 5.1 km/s, as is the angular speed of the Sun in its orbit about the Galactic center, (To + Vsun)/Ro = 30.57 +/- 0.43 km/s/kpc. These parameters improve the accuracy of estimates of the accelerations of the Sun and the Hulse-Taylor binary pulsar in their Galactic orbits, significantly reducing the uncertainty in tests of gravitational radiation predicted by general relativity.

We present a systematic analysis of the complete set of observations of the black hole (BH) binary XTE J1550-564 obtained by the Rossi X-ray Timing Explorer. We study the fast time variability properties of the source and determine the spin of the black hole through the relativistic precession model. Similarly to what is observed in the BH binary GRO J1655-40, the frequencies of the QPOs and broad band noise components match the general relativistic frequencies of particle motion close to the compact object predicted by the relativistic precession model. The combination of two simultaneously observed quasi-periodic oscillation (QPO) frequencies together with the dynamical BH mass from optical/infrared observations yields a spin equal to a = 0.34 +/- 0.01, consistent with previous determinations from X-ray spectroscopy. Based on the derived BH parameters, the low frequency QPO emission radii vary from about 30 gravitational radii to the innermost stable orbit for this spin (about 5 gravitational radii), where they sharply disappear as observed for the case of GRO J1655-40.

We report four trigonometric parallaxes for high-mass star-forming regions within 4 kpc of the Galactic center. These measurements were made with the VLBA as part of the BeSSeL Survey. By associating these sources kinematically with large-scale features in CO and HI longitude-velocity diagrams, we begin to outline some major features of the inner Milky Way: the Connecting arm, the near and far 3 kpc arms, and the Norma arm. The Connecting arm in the first Galactic quadrant lies closer to the Galactic center than the far 3 kpc arm and is offset by the long-bar's major axis near its leading edge, supporting the presence of an inner Lindblad resonance. Assuming the 3 kpc arms are a continuous physical structure, the relative Galactocentric distance of its near and far sides suggests highly elliptical streamlines of gas around the bar(s) and a bar corotation radius, r_CR > 3.6 kpc. At a Galactic longitude near 10{\deg} and a heliocentric distance of about 5 kpc, the near 3 kpc arm and the Norma arm intersect on a face-on view of our Galaxy, while passing at different Galactic latitudes. We provide an accurate distance measurement to the W31 star-forming complex of 4.95 (+0.5;-0.43) kpc from the Sun, which associates it with a bright CO feature belonging to the near 3 kpc arm.

We wish to investigate the structure, velocity field, and stellar content of the G16.59-0.05 high-mass star-forming region, where previous studies have established the presence of two almost perpendicular (NE-SW and SE-NW), massive outflows, and a rotating disk traced by methanol maser emission. We performed Very Large Array observations of the radio continuum and ammonia line emission, complemented by COMICS/Subaru and Hi-GAL/Herschel images in the mid- and far-infrared (IR). Our centimeter continuum maps reveal a collimated radio jet that is oriented E-W and centered on the methanol maser disk, placed at the SE border of a compact molecular core. The spectral index of the jet is negative, indicating non-thermal emission over most of the jet, except the peak close to the maser disk, where thermal free-free emission is observed. We find that the ammonia emission presents a bipolar structure consistent (on a smaller scale) in direction and velocity with that of the NE-SW bipolar outflow detected in previous CO observations. After analyzing our previous N2H+(1-0) observations again, we conclude that two scenarios are possible. In one case both the radio jet and the ammonia emission would trace the root of the large-scale CO bipolar outflow. The different orientation of the jet and the ammonia flow could be explained by precession and/or a non-isotropic density distribution around the star. In the other case, the N2H+(1-0) and ammonia bipolarity is interpreted as two overlapping clumps moving with different velocities along the line of sight. The ammonia gas also seems to undergo rotation consistent with the maser disk. Our IR images complemented by archival data allow us to derive a bolometric luminosity of about 10^4 L_sun and to conclude that most of the luminosity is due to the young stellar object associated with the maser disk.

We measured the phase-lag spectrum of the high-frequency quasi-periodic oscillations (QPO) in the black hole systems (at QPO frequencies) GRS 1915+105 (35 Hz and 67 Hz), GRO J1655-40 (300 Hz and 450 Hz), XTE J1550-564 (180 Hz and 280 Hz), and IGR J17091-3624 (67 Hz). The lag spectra of the 67-Hz QPO in, respectively, GRS 1915+105 and IGR J17091-3624, and the 450-Hz QPO in GRO J1655-40 are hard (hard photons lag the soft ones) and consistent with each other, with the hard lags increasing with energy. On the contrary, the lags of the 35-Hz QPO in GRS 1915+105 are soft, with the lags becoming softer as the energy increases; the lag spectrum of the 35-Hz QPO is inconsistent with that of the 67-Hz QPO. The lags of the 300-Hz QPO in GRO J1655-40, and the 180-Hz and the 280-Hz QPO in XTE J1550-564 are independent of energy, consistent with each other and with being zero or slightly positive (hard lags). For GRO J1655-40 the lag spectrum of the 300-Hz QPO differs significantly from that of the 450-Hz QPOs. The similarity of the lag spectra of the 180-Hz and 280-Hz QPO in XTE J1550-564 suggests that these two are the same QPO seen at a different frequency in different observations. The lag spectrum of the 67-Hz QPO in GRS 1915+105 is significantly different from that of the $2.7 \times 10^{-4}$ Hz QPO in the narrow-line Seyfert 1 galaxy RE J1034+396, which disproves the suggestion that the two QPOs are the same physical phenomenon with their frequencies scaled only by the black-hole mass. The lag spectrum of the QPO in RE J1034+396 is similar to that of the 35-Hz QPO in GRS 1915+105, although identifying these two QPOs as being the same physical feature remains problematic. We compare our results with those for the lags of the kilohertz QPOs in neutron-star systems and the broadband noise component in active galactic nuclei, and discuss possible scenarios for producing the lags in these systems.

We analysed the X-ray spectra of six observations, simultaneously taken with XMM-Newton and Rossi X-ray Timing Explorer (RXTE), of the neutron star low-mass X-ray binary 4U 1636-53. The observations cover several states of the source, and therefore a large range of inferred mass accretion rate. These six observations show a broad emission line in the spectrum at around 6.5 keV, likely due to iron. We fitted this line with a set of phenomenological models of a relativistically broadened line, plus a model that accounts for relativistically smeared and ionised reflection from the accretion disc. The latter model includes the incident emission from both the neutron-star surface or boundary layer and the corona that is responsible for the high-energy emission in these systems. From the fits with the reflection model we found that in four out of the six observations the main contribution to the reflected spectrum comes from the neutron-star surface or boundary layer, whereas in the other two observations the main contribution to the reflected spectrum comes from the corona. We found that the relative contribution of these two components is not correlated to the state of the source. From the phenomenological models we found that the iron line profile is better described by a symmetric, albeit broad, profile. The width of the line cannot be explained only by Compton broadening, and we therefore explored the case of relativistic broadening. We further found that the direct emission from the disc, boundary layer, and corona generally evolved in a manner consistent with the standard accretion disc model, with the disc and boundary layer becoming hotter and the disc moving inwards as the source changed from the hard in to the soft state. The iron line, however, did not appear to follow the same trend.

(abridged) We studied the energy and frequency dependence of the Fourier time lags and intrinsic coherence of the kHz QPOs in the NS LMXBs 4U 1608-52 and 4U 1636-53 using RXTE data. In both sources we confirmed energy-dependent soft lags of 10-100 \mu s for the lower kHz QPO. We also found that the time lags of the upper kHz QPO are independent of energy and inconsistent with the soft lags of the lower kHz QPO. The intrinsic coherence of the lower kHz QPO remains constant at 0.6 from 5 to 12 keV, and then drops to zero, while for the upper kHz QPO the intrinsic coherence is consistent with zero across the full energy range. The intrinsic coherence of the upper kHz QPO is consistent with zero over the full frequency range of the QPO, except in 4U 1636-53 at ~780 Hz where it increases to 0.13. In 4U 1636-53, for the lower kHz QPO the 4-12 keV photons lag the 12-20 keV ones by 25 \mu s in the QPO frequency range 500-850 Hz, with the lags decreasing to 15 \mu s at higher frequencies. In 4U 1608-52 the soft lags of the lower kHz QPO remain constant at 40 \mu s. In 4U 1636-53, for the upper kHz QPO the 12-20 keV photons lag the 4-12 keV ones by 11 +/- 3 \mu s, independent of QPO frequency; we found consistent results for the time lags of the upper kHz QPO in 4U 1608-52. The intrinsic coherence of the lower kHz QPO increases from ~0-0.4 at 600 Hz to 1 and 0.6 at 800 Hz in 4U 1636-53 and 4U 1608-52, respectively. In 4U 1636-53 it decreases to 0.5 at 920 Hz, while in 4U 1608-52 we do not have data above 800 Hz. We discuss our results in the context of scenarios in which the soft lags are either due to reflection off the accretion disc or up-/down-scattering in a hot medium close to the neutron star. We finally explore the connection between, on one hand the time lags and the intrinsic coherence of the kHz QPOs, and on the other the QPOs' amplitude and quality factor in these two sources.

This work presents a detailed study of the gas kinematics towards the "Hot Molecular Core" (HMC) G31.41+0.31 via multi-epoch VLBI observations of the H2O 22 GHz and CH3OH 6.7 GHz masers, and single-epoch VLBI of the OH 1.6 GHz masers. Water masers present a symmetric spatial distribution with respect to the HMC center, where two nearby (0.2" apart), compact, VLA sources (labeled "A" and "B") are previously detected. The spatial distribution of a first group of water masers, named "J1", is well fit with an elliptical profile, and the maser proper motions mainly diverge from the ellipse center, with average speed of 36 km s-1. These findings strongly suggest that the "J1" water maser group traces the heads of a young (dynamical time of 1.3 10^3 yr), powerful (momentum rate of ~0.2 M_sun yr-1 km s-1), collimated (semi-opening angle ~10 deg) jet emerging from a MYSO located close (within 0.15") to the VLA source "B". Most of the water features not belonging to "J1" present an elongated (about 2" in size), NE--SW oriented (PA = 70 deg), S-shape distribution, which we denote with the label "J2". The elongated distribution of the "J2" group and the direction of motion, approximately parallel to the direction of elongation, of most "J2" water masers suggests the presence of another collimated outflow, emitted from a MYSO near the VLA source "A". The orientation of the "J2" jet agrees well with that (PA = 68 deg) of the well-defined V_LSR gradient across the HMC revealed by previous interferometric, thermal line observations. Furthermore, the "J2" jet is powerful enough to sustain the large momentum rate, 0.3 M_sun yr-1 km s-1, estimated assuming that the V_LSR gradient represents a collimated outflow. These two facts lead us to favour the interpretation of the V_LSR gradient across the G31.41+0.31 HMC in terms of a compact and collimated outflow.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

We present the results of the analysis of a large database of X-ray observations of 22 galactic black-hole transients with the Rossi X-Ray timing explorer throughout its operative life for a total exposure time of ~12 Ms. We excluded persistent systems and the peculiar source GRS 1915+105, as well as the most recently discovered sources. The semi-automatic homogeneous analysis was aimed at the detection of high-frequency (100-1000 Hz) quasi-periodic oscillations (QPO), of which several cases were previously reported in the literature. After taking into account the number of independent trials, we obtained 11 detections from two sources only: XTE J1550-564 and GRO J1655-40. For the former, the detected frequencies are clustered around 180 Hz and 280 Hz, as previously found. For the latter, the previously-reported dichotomy 300-450 Hz is found to be less sharp. We discuss our results in comparison with kHz QPO in neutron-star X-ray binaries and the prospects for future timing X-ray missions.

We analysed all archival RXTE observations of the neutron-star low-mass X-ray binary 4U 1636-53 up to May 2010. In 528 out of 1280 observations we detected kilohertz quasi-periodic oscillations (kHz QPOs), with ~ 65% of these detections corresponding to the so-called lower kHz QPO. Using this QPO we measured, for the first time, the rate at which the QPO frequency changes as a function of QPO frequency. For this we used the spread of the QPO frequency over groups of 10 consecutive measurements, sampling timescales between 320 and 1600 s, and the time derivative of the QPO frequency over timescales of 32 to 160 s. We found that: (i) Both the QPO-frequency spread and the QPO time derivative decrease by a factor ~ 3 as the QPO frequency increases. (ii) The average value of the QPO time derivative decreases by a factor of ~ 2 as the timescale over which the derivative is measured increases from less than 64 s to 160 s. (iii) The relation between the absolute value of the QPO time derivative and the QPO frequency is consistent with being the same both for the positive and negative QPO-frequency derivative. We show that, if either the lower or the upper kHz QPO reflects the Keplerian frequency at the inner edge of the accretion disc, these results support a scenario in which the inner part of the accretion disc is truncated at a radius that is set by the combined effect of viscosity and radiation drag.

Whether the Cygnus X complex consists of one physically connected region of star formation or of multiple independent regions projected close together on the sky has been debated for decades. The main reason for this puzzling scenario is the lack of trustworthy distance measurements. We aim to understand the structure and dynamics of the star-forming regions toward Cygnus X by accurate distance and proper motion measurements. To measure trigonometric parallaxes, we observed 6.7 GHz methanol and 22 GHz water masers with the European VLBI Network and the Very Long Baseline Array. We measured the trigonometric parallaxes and proper motions of five massive star-forming regions toward the Cygnus X complex and report the following distances within a 10% accuracy: 1.30+-0.07 kpc for W 75N, 1.46^{+0.09}_{-0.08} kpc for DR 20, 1.50^{+0.08}_{-0.07} kpc for DR 21, 1.36^{+0.12}_{-0.11} kpc for IRAS20290+4052, and 3.33+-0.11kpc for AFGL 2591. While the distances of W 75N, DR 20, DR 21, and IRAS 20290+4052 are consistent with a single distance of 1.40+-0.08 kpc for the Cygnus X complex, AFGL 2591 is located at a much greater distance than previously assumed. The space velocities of the four star-forming regions in the Cygnus X complex do not suggest an expanding Stroemgren sphere.

We report on a detailed study of the water maser kinematics and radio continuum emission toward the most massive and young object in the star-forming region AFGL 2591. Our analysis shows at least two spatial scales of multiple star formation, one projected across 0.1 pc on the sky and another one at about 2000 AU from a ZAMS star of about 38 Solar masses. This young stellar object drives a powerful jet- and wind-driven outflow system with the water masers associated to the outflow walls, previously detected as a limb-brightened cavity in the NIR band. At about 1300 AU to the north of this object a younger protostar drives two bow shocks, outlined by arc-like water maser emission, at 200 AU either side of the source. We have traced the velocity profile of the gas that expands along these arc-like maser structures and compared it with the jet-driven outflow model. This analysis suggests that the ambient medium around the northern protostar is swept up by a jet-driven shock (>66 km/s) and perhaps a lower-velocity (~10 km/s) wind with an opening angle of about 20 degrees from the jet axis.

We report a trigonometric parallax measurement with the Very Long Baseline Array for the water maser in the distant high-mass star-forming region G75.30+1.32. This source has a heliocentric distance of 9.25+-0.45 kpc, which places it in the Outer arm in the first Galactic quadrant. It lies 200 pc above the Galactic plane and is associated with a substantial HI enhancement at the border of a large molecular cloud. At a Galactocentric radius of 10.7 kpc, G75.30+1.32 is in a region of the Galaxy where the disk is significantly warped toward the North Galactic Pole. While the star-forming region has an instantaneous Galactic orbit that is nearly circular, it displays a significant motion of 18 km/s toward the Galactic plane. The present results, when combined with two previous maser studies in the Outer arm, yield a pitch angle of about 12 degrees for a large section of the arm extending from the first quadrant to the third.

Observational signatures of infalling envelopes and outflowing material in early stages of protostellar evolution, and at small radii from the protostar, are essential to progress in the understanding of the mass-accretion process in star formation. In this letter, we report a detailed study of the accretion and outflow structure around a protostar in the well-known high-mass star-forming region AFGL 5142. We focus on the mm source MM-1, which exhibits hot-core chemistry, radio continuum emission, and strong water and methanol masers. Remarkably, our Very Long Baseline Interferometry (VLBI) observations of molecular masers over six years provided us with the 3-D velocity field of circumstellar molecular gas with a resolution of 0.001-0.005 arcseconds and at radii <0.23 arcseconds (or 400 AU) from the protostar. In particular, our measurements of methanol maser emission enabled, for the first time, a direct measurement of infall of a molecular envelope (radius of 300 AU and velocity of 5 km/s) onto an intermediate- to high-mass protostar. We estimate an infall rate of 0.0006 n_8 Msun/year, where n_8 is the ambient volume density in units of 10^8 cm-3 (required for maser excitation). In addition, our measurements of water maser (and radio continuum) emission identify a collimated bipolar molecular outflow (and ionized jet) from MM-1. The evidence of simultaneous accretion and outflow at small spatial scales, makes AFGL 5142 an extremely compelling target for high-angular resolution studies of high-mass star formation.

This work presents a study of the milliarcsecond structure of the 6.7 GHz methanol masers at high-velocity resolution (0.09 km s^(-1)) in four high-mass star-forming regions: G16.59-0.05, G23.01-0.41, IRAS20126+4104, and AFGL5142. We studied these sources by means of multi-epoch VLBI observations in the 22 GHz water and 6.7 GHz methanol masers, to determine the 3-D gas kinematics within a few thousand AU from the (proto)star. The present work exploits the 6.7 GHz maser data collected so far to investigate the milliarcsecond structure of this maser emission at high-velocity resolution. Most of the detected 6.7 GHz maser features present an ordered (linear, or arc-like) distribution of maser spots on the plane of the sky, together with a regular variation in the spot LSR velocity (V_LSR) with position. Typical values for the amplitude of the V_LSR gradients (defined in terms of the derivative of the spot V_LSR with position) are found to be 0.1-0.2 km s^(-1) mas^(-1). In each of the four target sources, the orientation and the amplitude of most of the feature V_LSR gradients remain remarkably stable in time, on timescales of (at least) several years. We also find that the data are consistent with having the V_LSR gradients and proper motion vectors in the same direction on the sky, considered the measurement uncertainties. The time persistency, the ordered angular and spatial distribution, and the orientation generally similar to the proper motions, altogether suggest a kinematical interpretation for the origin of the 6.7 GHz maser V_LSR gradients. This work shows that the organized motions (outflow, infall, and rotation) revealed by the (22 GHz water and 6.7 GHz methanol) masers on large scales (~100-1000 AU) also persist to very small (~10 AU) scales.

We present the results of a survey of radio continuum sources near the Galactic plane using the Very Long Baseline Array (VLBA). Our observations are designed to identify compact extragalactic sources of milliarcsecond size that can be used for parallax measurements in the Bar and Spiral Structure Legacy Survey. We selected point sources from the NVSS and CORNISH catalogs with flux densities above 30 mJy and within $1.5\degr$ of known maser targets. Of the 1529 sources observed, 199 were detected. For sources detected on 3 or more baselines, we determined accurate positions and evaluated their quality as potential calibrators. Most of the 1330 sources that were not detected with the VLBA are probably of extragalactic origin.

Astrometric Very Long Baseline Interferometry (VLBI) observations of maser sources in the Milky Way are used to map the spiral structure of our Galaxy and to determine fundamental parameters such as the rotation velocity ($\Theta_0$) and curve and the distance to the Galactic center (R$_0$). Here, we present an update on our first results, implementing a recent change in the knowledge about the Solar motion. It seems unavoidable that the IAU recommended values for R$_0$ and $\Theta_0$ need a substantial revision. In particular the combination of 8.5 kpc and 220 \kms\, can be ruled out with high confidence. Combining the maser data with the distance to the Galactic center from stellar orbits and the proper motion of Sgr\,A* gives best values of R$_0$ = 8.3 $\pm$ 0.23 kpc and $\Theta_0$ = 239 or 246 $\pm$ 7 \kms, for Solar motions of V$_ \odot$ = 12.23 and 5.25 \kms, respectively. Finally, we give an outlook to future observations in the Bar and Spiral Structure Legacy (BeSSeL) Survey.

Cygnus X is one of the closest giant molecular cloud complexes and therefore an extensively studied region of ongoing high mass star formation. However, the distance to this region has been a long-standing issue, since sources at galactic longitude of ~80 degrees could be in the Local Arm nearby (1-2 kpc), in the Perseus Arm at ~5 kpc, or even in the outer arm (~10 kpc). We use combined observations of the EVN plus two Japanese stations to measure very accurate parallaxes of methanol masers in five star-forming regions in Cygnus X to understand if they belong to one large star-forming complex or if they are separate entities located at different distances. Here we report our preliminary result for W75N based on six epochs of VLBI observations: we find that W75N is at a distance of 1.32^{+0.11}_{-0.09} kpc, which is significantly closer than the reported values in the literature (1.5-2 kpc).

We analysed 866 observations of the neutron-star low-mass X-ray binary XTE J1701-462 during its 2006-2007 outburst. XTE J1701-462 is the only example so far of a source that during an outburst showed, beyond any doubt, spectral and timing characteristics both of the Z and atoll type. We found that the lower kHz QPO in the atoll phase has a significantly higher coherence and fractional rms amplitude than any of the kHz QPOs seen during the Z phase, and that in the same frequency range, atoll lower kHz QPOs show coherence and fractional rms amplitude, respectively, 2 and 3 times larger than the Z kHz QPOs. Out of the 707 observations in the Z phase, there is no single observation in which the kHz QPOs have a coherence or rms amplitude similar to those seen when XTE J1701-462 was in the atoll phase, even though the total exposure time was about 5 times longer in the Z than in the atoll phase. Since it is observed in the same source, the difference in QPO coherence and rms amplitude between the Z and atoll phase cannot be due to neutron-star mass, magnetic field, spin, inclination of the accretion disk, etc. If the QPO frequency is a function of the radius in the accretion disk in which it is produced, our results suggest that in XTE J1701-462 the coherence and rms amplitude are not uniquely related to this radius. Here we argue that this difference is instead due to a change in the properties of the accretion flow around the neutron star. Regardless of the precise mechanism, our result shows that effects other than the geometry of space time around the neutron star have a strong influence on the coherence and rms amplitude of the kHz QPOs, and therefore the coherence and rms amplitude of the kHz QPOs cannot be simply used to deduce the existence of the innermost stable circular orbit around a neutron star.

The present paper focuses on the high-mass star-forming region G23.01-0.41. Methods: Using the VLBA and the EVN arrays, we conducted phase-referenced observations of the three most powerful maser species in G23.01-0.41: H2O at 22.2 GHz (4 epochs), CH3OH at 6.7 GHz (3 epochs), and OH at 1.665 GHz (1 epoch). In addition, we performed high-resolution (> 0".1), high-sensitivity (< 0.1 mJy) VLA observations of the radio continuum emission from the HMC at 1.3 and 3.6 cm. Results: We have detected H2O, CH3OH, and OH maser emission clustered within 2000 AU from the center of a flattened HMC, oriented SE-NW, from which emerges a massive 12CO outflow, elongated NE-SW, extended up to the pc-scale. Although the three maser species show a clearly different spatial and velocity distribution and sample distinct environments around the massive YSO, the spatial symmetry and velocity field of each maser specie can be explained in terms of expansion from a common center, which possibly denotes the position of the YSO driving the maser motion. Water masers trace both a fast shock (up to 50 km/s) closer to the YSO, powered by a wide-angle wind, and a slower (20 km/s) bipolar jet, at the base of the large-scale outflow. Since the compact free-free emission is found offset from the putative location of the YSO along a direction consistent with that of the maser jet axis, we interpret the radio continuum in terms of a thermal jet. The velocity field of methanol masers can be explained in terms of a composition of slow (4 km/s in amplitude) motions of radial expansion and rotation about an axis approximately parallel to the maser jet. Finally, the distribution of line of sight velocities of the hydroxyl masers suggests that they can trace gas less dense (n(H2) < 10^6 cm^-3) and more distant from the YSO than that traced by the water and methanol masers, which is expanding toward the observer. (Abridged)

The present paper focuses on the high-mass star-forming region G16.59-0.05. Methods: Using the VLBA and the EVN arrays, we conducted phase-referenced observations of the three most powerful maser species in G16.59-0.05: H2O at 22.2 GHz (4 epochs), CH3OH at 6.7 GHz (3 epochs), and OH at 1.665 GHz (1 epoch). In addition, we performed high-resolution (> 0".1), high-sensitivity (< 0.1 mJy) VLA observations of the radio continuum emission from the star-forming region at 1.3 and 3.6 cm. Results: This is the first work to report accurate measurements of the "relative" proper motions of the 6.7 GHz CH3OH masers. The different spatial and 3-D velocity distribution clearly indicate that the 22 GHz water and 6.7 GHz methanol masers are tracing different kinematic environments. The bipolar distribution of 6.7 GHz maser l.o.s. velocities and the regular pattern of observed proper motions suggest that these masers are tracing rotation around a central mass of about 35 solar masses. The flattened spatial distribution of the 6.7 GHz masers, oriented NW-SE, suggests that they can originate in a disk/toroid rotating around the massive YSO which drives the 12CO(2-1) outflow, oriented NE-SW, observed on arcsec scale. The extended, radio continuum source observed close to the 6.7 GHz masers could be excited by a wide-angle wind emitted from the YSO associated with the methanol masers, and such a wind is proven to be sufficiently energetic to drive the NE-SW 12CO(2-1) outflow. The H2O masers distribute across a region offset about 0".5 to the NW of the CH3OH masers, in the same area where emission of high-density molecular tracers, typical of HMCs, was detected. We postulate that a distinct YSO, possibly in an earlier evolutionary phase than that exciting the methanol masers, is responsible for the excitation of the water masers and the HMC molecular lines. (Abridged)

We report a trigonometric parallax of 12 GHz methanol masers associated with the massive star forming region G9.62+0.20, corresponding to a distance of 5.2+-0.6 kpc. With a local standard of rest velocity of about 2 km/s, the region's kinematic distances of 0.5 and 16 kpc differ greatly from the distance derived here. Our measurement of the peculiar motion of the star forming region shows a very large deviation from a circular Galactic orbit: 41 km/s radially outward from the Galactic center and 60 km/s counter to Galactic rotation. The combination of its radial velocity and distance places G9.62+0.20 in the inner region of the Galaxy close to the Expanding Near 3 kpc-Arm, where the bulge/bar potential has strong gravitational influence. We also map the distribution of 12 GHz methanol masers, locate them with respect to a hypercompact H II region, and compare our data with the periodic flare phenomenon reported previously for this source.

Most previous high-angular (<0.1 arcsec) resolution studies of molecular masers in high-mass star forming regions (SFRs) have concentrated mainly on either water or methanol masers. While high-angular resolution observations have clarified that water masers originate from shocks associated with protostellar jets, different environments have been proposed in several sources to explain the origin of methanol masers. Tha aim of the paper is to investigate the nature of the methanol maser birthplace in SFRs and the association between the water and methanol maser emission in the same young stellar object. We have conducted phase-reference Very Long Baseline Interferometry (VLBI) observations of water and methanol masers toward two high-mass SFRs, Sh 2-255 IR and AFGL 5142. In Sh 2-255 IR water masers are aligned along a direction close to the orientation of the molecular outflow observed on angular scales of 1-10 arcsec, tracing possibly the disk-wind emerging from the disk atmosphere. In AFGL 5142 water masers trace expansion at the base of a protostellar jet, whilst methanol masers are more probably tracing infalling than outflowing gas. The results for AFGL 5142 suggest that water and methanol masers trace different kinematic structures in the circumstellar gas.