Source author record

W. H. T. Vlemmings

W. H. T. Vlemmings appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

60works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

60 published item(s)

preprint2022arXiv

First detection of AlF line emission towards M-type AGB stars

The nucleosynthesis production of fluorine (F) is still a matter of debate. Asymptotic giant branch (AGB) stars are one of the main candidates for F production. However, their contribution to the total F budget is not fully known due to the lack of observations. In this paper, we report the detection of AlF line emission, one of the two main carriers of F in the gas-phase in the outflow of evolved stars, towards five nearby oxygen-rich AGB stars, $o$ Ceti, R Leo, IK Tau, R Dor, and W Hya. From spatially resolved observations, we estimated the AlF emitting region with a radius $\sim11R_{\star}$ for $o$ Ceti and $\sim9R_{\star}$ for R Leo. From population diagram analysis, we report the AlF column densities of $\sim 5.8\times10^{15}$ cm$^{-2}$ and $\sim 3\times10^{15}$ cm$^{-2}$ for $o$ Ceti and R Leo, respectively, within these regions. For $o$ Ceti, we used the C$^{18}$O ($v=0$, $J=3-2$) observations to estimate the H$_2$ column density of the emitting region. We found a fractional abundance of $f_{\rm AlF/H_2}\sim(2.5\pm1.7)\times10^{-8}$. This gives a lower limit on the F budget in $o$ Ceti and is compatible with the solar F budget $f_{\rm F/H_2}=(5\pm2)\times10^{-8}$. For R Leo, a fractional abundance $f_{\rm AlF/H_2}=(1.2\pm0.5)\times10^{-8}$ is estimated. For other sources, we cannot precisely determine the emitting region based on the available data. Assuming an emitting region with a radius of $\sim 11R_{\star}$ and the rotational temperatures derived for $o$ Ceti and R Leo, we crudely approximated the AlF column density to be $\sim(1.2-1.5)\times10^{15}$ cm$^{-2}$ in W Hya, $\sim(2.5-3.0)\times10^{14}$ cm$^{-2}$ in R Dor, and $\sim(0.6-1.0)\times10^{16}$ cm$^{-2}$ in IK Tau. These result in fractional abundances within a range of $f_{\rm AlF/H_2}\sim(0.1-4)\times10^{-8}$ in W Hya, R Dor, and IK Tau.

preprint2021arXiv

EVN observations of 6.7 GHz methanol maser polarization in massive star-forming regions. V. Completion of the flux-limited sample

Although the role of magnetic fields in launching molecular outflows in massive YSOs has been convincingly demonstrated by theoretical arguments, observationally, the alignment of the magnetic field lines with the molecular outflows is still under debate. We aim to complete the measurements of the direction of the magnetic fields at mas resolution around a sample of massive star-forming regions (MSFRs) to determine whether the magnetic field and outflows are aligned. In 2012, we started a large VLBI campaign with the EVN to measure the magnetic field orientation and strength toward a sample of 31 MSFRs (the flux-limited sample) by analyzing the polarized emission of 6.7GHz CH3OH masers. In the previous papers of the series, we have presented 80% of the sample. Here, we report the linearly and circularly polarized emission of 6.7GHz CH3OH masers toward the last five MSFRs of the flux-limited sample. The sources are G30.70-0.07, G30.76-0.05, G31.28+0.06, G32.03+0.06, and G69.52-0.97. We detected a total of 209 masers, 15% of which show linearly polarized emission (0.07%-16.7%), and 2% of which show circularly polarized emission (0.2%-4.2%). Zeeman splitting was measured toward G30.70-0.07, G32.03+0.06, and G69.52-0.97. The statistical analysis of the entire flux-limited sample shows that the observations are consistent with a bimodal distribution in the difference between the 3D magnetic field direction and the outflow axis, with half the magnetic field directions being perpendicular and the other half being parallel to the outflow. In addition, we determined that typical values of the linear and circular polarization fractions for 6.7 GHz CH3OH masers are Pl=1.0%-2.5% and Pv=0.5%-0.75%, respectively. We found that a typical Zeeman splitting is in the range between 0.5 m/s and 2.0 m/s. This would correspond to 9 mG<$|B_{||}|$<40 mG if F=3->4 is the most favored hyperfine transition.

preprint2020arXiv

ALMA and VLA reveal the lukewarm chromospheres of the nearby red supergiants Antares and Betelgeuse

We first present spatially resolved ALMA and VLA continuum observations of the early-M red supergiant Antares to search for the presence of a chromosphere at radio wavelengths. We resolve the free-free emission of the Antares atmosphere at 11 unique wavelengths between 0.7 mm (ALMA band 8) and 10 cm (VLA S band). The projected angular diameter is found to continually increase with increasing wavelength, from a low of 50.7 mas at 0.7 mm up to a diameter of 431 mas at 10 cm, which corresponds to 1.35 and 11.6 times the photospheric angular diameter, respectively. All four ALMA measurements show that the shape of the atmosphere is elongated, with a flattening of 15% at a similar position angle. The disk-averaged gas temperature of the atmosphere initially rises from a value of 2700 K at 1.35 $R_{\star}$ (i.e., 0.35 $R_{\star}$ above the photosphere) to a peak value of 3800 K at $\sim$2.5 $R_{\star}$, after which it then more gradually decreases to 1650 K at 11.6 $R_{\star}$. The rise in gas temperature between 1.35 $R_{\star}$ and $\sim$2.5 $R_{\star}$ is evidence for a chromospheric temperature rise above the photosphere of a red supergiant. We detect a clear change in the spectral index across the sampled wavelength range, with the flux density ${S_ν} \propto ν^{1.42}$ between 0.7 mm and 1.4 cm, which we associate with chromosphere-dominated emission, while the flux density ${S_ν} \propto ν^{0.8}$ between 4.3 cm and 10 cm, which we associate with wind-dominated emission. We then perform nonlocal thermal equilibrium modeling of the far-ultraviolet radiation field of another early-M red supergiant, Betelgeuse, and find that an additional hot (i.e., $>7000$ K) chromospheric photoionization component with a much smaller filling factor must also exist throughout the chromospheres of these stars.

preprint2020arXiv

ALMA reveals the coherence of the magnetic field geometry in OH 231.8+4.2

In a continuing effort to investigate the role of magnetic fields in evolved low and intermediate mass stars (principally regarding the shaping of their envelopes), we present new ALMA high resolution polarization data obtained for the nebula OH 231.8+4.2. We found that the polarized emission likely arises from aligned grains in the presence of magnetic fields rather than radiative alignment and self scattering. The ALMA data show well organized electric field orientations in most of the nebula and the inferred magnetic field vectors (rotated by 90 degrees) trace an hourglass morphology centred on the central system of the nebula. One region in the southern part of OH 231.8+4.2 shows a less organized distribution probably due to the shocked environment. These findings, in conjunction with earlier investigations (maser studies and dust emission analysis at other scales and wavelengths) suggest an overall magnetic hourglass located inside a toroidal field. We propose the idea that the magnetic field structure is closely related to the architecture of a magnetic tower and that the outflows were therefore magnetically launched. While the current dynamical effect of the fields might be weak in the equatorial plane principally due to the evolution of the envelope, it would still be affecting the outflows. In that regard, the measurement of the magnetic field at the stellar surface, which is still missing, combined with a full MHD treatment are required to better understand and constrain the events occurring in OH 231.8+4.2.

preprint2020arXiv

Characterizing the radio continuum nature of sources in the massive star-forming region W75N (B)

The massive star-forming region W75N~(B) is thought to host a cluster of massive protostars (VLA~1, VLA~2, and VLA~3) undergoing different evolutionary stages. In this work, we present radio continuum data with the highest sensitivity and angular resolution obtained to date in this region, using the VLA-A and covering a wide range of frequencies (4-48~GHz), which allowed us to study the morphology and the nature of the emission of the different radio continuum sources. We also performed complementary studies with multi-epoch VLA data and ALMA archive data at 1.3 mm wavelength. We find that VLA~1 is driving a thermal radio jet at scales of $\approx$0.1 arcsec ($\approx$130 au), but also shows signs of an incipient hyper-compact HII region at scales of $\lesssim$ 1 arcsec ($\lesssim$ 1300~au). VLA~3 is also driving a thermal radio jet at scales of a few tenths of arcsec (few hundred of au). We conclude that this jet is shock-exciting the radio continuum sources Bc and VLA~4 (obscured HH objects), which show proper motions moving outward from VLA~3 at velocities of $\approx$112--118~km/s. We have also detected three new weak radio continuum sources, two of them associated with millimeter continuum cores observed with ALMA, suggesting that these two sources are also embedded YSOs in this massive star-forming region.

preprint2020arXiv

CO and HCN isotopologue ratios in the outflows of AGB stars

Isotopologue line intensity ratios of circumstellar molecules have been widely used to trace the photospheric elemental isotopic ratios of evolved stars. However, depending on the molecular species and the physical conditions of the environment, the circumstellar isotopologue ratio may deviate considerably from the stellar atmospheric value. In this paper, we aim to examine how the CO and HCN abundance ratios vary radially due to chemical reactions in the outflows of AGB stars and the effect of excitation and optical depth on the resulting line intensity ratios. We find that the circumstellar 12CO/13CO can deviate from its atmospheric value by up to 25-94% and 6-60% for C- and O-type CSEs, respectively. We show that variations of the intensity of the ISRF and the gas kinetic temperature can significantly influence the CO isotopologue ratio in the outer CSEs. On the contrary, the H12CN/H13CN ratio is stable for all tested mass-loss rates. The RT modeling shows that the integrated line intensity ratio of CO of different rotational transitions varies significantly for stars with intermediate mass-loss rates due to combined chemical and excitation effects. In contrast, the excitation conditions for the both HCN isotopologues are the same. We demonstrate the importance of using the isotopologue abundance profiles from chemical models as inputs to RT models in the interpretation of isotopologue observations. Previous studies of CO isotopologue ratios are based on multi-transition data for individual sources and it is difficult to estimate the errors in the reported values due to assumptions that are not entirely correct according to this study. If anything, previous studies may have overestimated the circumstellar 12CO/13CO abundance ratio. The use of the HCN as a tracer of C isotope ratios is affected by fewer complicating problems, provided one accounts corrections for high optical depths.

preprint2020arXiv

DEATHSTAR: Nearby AGB stars with the Atacama Compact Array I. CO envelope sizes and asymmetries: A new hope for accurate mass-loss-rate estimates

This is the first publication of the DEATHSTAR project. The goal of the project is to reduce the uncertainties of observational estimates of mass-loss rates from Asymptotic Giant Branch (AGB) stars. Line emission from 12CO J=2-1 and 3-2 were mapped using the ACA. In this initial analysis, the emission distribution was fit to a Gaussian distribution in the uv-plane. Detailed radiative transfer analysis will be presented in the future. The axes of the best-fit Gaussian at the line center of the 12CO J=2-1 emission gives a first indication of the size of the emitting region. Furthermore, the fitting results, such as the major and minor axis, center position, and the goodness of fit across both lines, constrain the symmetry of the emission distribution. We find that the CO envelope sizes are, in general, larger for C-type than for M-type AGB stars, which is expected if the CO/H2 ratio is larger in C-type stars. Furthermore, a relation between the 12CO J=2-1 size and circumstellar density is shown that, while in broad agreement with photodissociation calculations, reveals large scatter and systematic differences between the stellar types. The majority of the sources have CO envelopes that are consistent with a spherically symmetric, smooth outflow. For about a third of the sources, indications of strong asymmetries are found. This is consistent with previous interferometric investigations of northern sources. Smaller scale asymmetries are found in a larger fraction of sources. These results for CO envelope radii and shapes can be used to constrain detailed radiative transfer modeling of the same stars so as to determine mass-loss rates that are independent of photodissociation models. For a large fraction of the sources, observations at higher spatial resolution will be necessary to further investigate the complex circumstellar dynamics revealed by our ACA observations.

preprint2020arXiv

Inner dusty envelope of the AGB stars W~Hydrae, SW\,Virginis, and R~Crateris using SPHERE/ZIMPOL

Stars with initial masses between $\sim0.8$ and 8~$M_\odot$ present copious mass loss during the asymptotic giant branch (AGB) at the end of their lives. The accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the extended AGB stellar atmospheres. The details of this process are not yet well understood, however. Using the extreme-adaptive-optics imager and polarimeter SPHERE/ZIMPOL, we observed light polarised by grains around W\,Hya, SW\,Vir, and R\,Crt, which have mass-loss rates between 10$^{-7}$ and 10$^{-6}~M_\odot~{\rm yr^{-1}}$. We find the distribution of dust to be asymmetric around the three targets. A biconical morphology is seen for R Crt, with a position angle that is very similar to those inferred from interferometric observations of maser emission and of mid-infrared continuum emission. The cause of the biconical outflow cannot be directly inferred from the ZIMPOL data. The dust grains polarise light more efficiently at 0.65~$μ$m for R\,Crt and SW\,Vir and at 0.82~$μ$m for W\,Hya. This indicates that at the time of the observations, the grains around SW\,Vir and R\,Crt had sizes $< 0.1~μ$m, while those around W\,Hya were larger, with sizes $\gtrsim 0.1~μ$m. The asymmetric distribution of dust around R\,Crt makes the interpretation more uncertain for this star, however. We find that polarised light is produced already from within the visible photosphere of W~Hya, which we reproduce using models with an inner dust shell that is optically thick to scattering. The radial profile of the polarised light observed around W\,Hya reveal a steep dust density profile. We find the wind-acceleration region of W\,Hya to extend to at least $\sim 7~R_\star$, in agreement with theoretical predictions of acceleration up to $\sim 12~R_\star$.

preprint2020arXiv

Mapping circumstellar magnetic fields of late-type evolved stars with the Goldreich-Kylafis effect: CARMA observations at $λ1.3$ mm of R Crt and R Leo

Mapping magnetic fields is the key to resolving what remains an unclear physical picture of circumstellar magnetic fields in late-type evolved stars. Observations of linearly polarized emission from thermal molecular line transitions due to the Goldreich-Kylafis (G-K) effect provides valuable insight into the magnetic field geometry in these sources that is complementary to other key studies. In this paper, we present the detection of spectral-line polarization from both the thermal $J=2-1$ CO line and the $v=1, J=5-4$ SiO maser line toward two thermal-pulsating (TP-) AGB stars, R Crt and R Leo. The observed fractional linear polarization in the CO emission is measured as $m_l\sim 3.1\%$ and $m_l\sim9.7\%$ for R Crt and R Leo respectively. A circumstellar envelope (CSE) model profile and the associated parameters are estimated and used as input to a more detailed modeling of the predicted linear polarization expected from the G-K effect. The observed thermal line polarization level is consistent with the predicted results from the G-K model for R Crt; additional effects need to be considered for R Leo.

preprint2019arXiv

The extended molecular envelope of the asymptotic giant branch star $π^{1}$ Gruis as seen by ALMA II. The spiral-outflow observed at high-angular resolution

The AGB star $π^{1}$ Gruis has a known companion (at a separation of ~400 AU) which cannot explain the strong deviations from the spherical symmetry of the CSE. Recently, hydrodynamic simulations of mass transfer in closer binary systems have successfully reproduced the spiral-shaped CSEs found around a handful of sources. There is growing evidence for an even closer, undetected companion complicating the case of $π^{1}$ Gruis further. The improved spatial resolution allows for the investigation of the complex circumstellar morphology and the search for imprints on the CSE of the third component. We have observed the 12CO J=3-2 line emission from $π^{1}$ Gruis using both the compact and extended array of Atacama Large Millimeter/submillimeter Array (ALMA). The interferometric data has furthermore been combined with data from the ALMA total power (TP) array. The imaged brightness distribution has been used to constrain a non-local, non-LTE 3D radiative transfer model of the CSE. The high-angular resolution ALMA data have revealed the first example of a source on the AGB where both a faster bipolar outflow and a spiral pattern along the orbital plane can be seen in the gas envelope. The spiral can be traced in the low- to intermediate velocity, from 13 to 25 km s$^{-1}$, equatorial torus. The largest spiral-arm separation is $\approx$5".5 and consistent with a companion with an orbital period of $\approx$330 yrs and a separation of less than 70 AU. The kinematics of the bipolar outflow is consistent with it being created during a mass-loss eruption where the mass-loss rate from the system increased by at least a factor of 5 during 10-15 yrs. The spiral pattern is the result of an undetected companion. The bipolar outflow is the result of a rather recent mass-loss eruption event.

preprint2016arXiv

ALMA observations of the vibrationally-excited rotational CO transition $v=1, J=3-2$ towards five AGB stars

We report the serendipitous detection with ALMA of the vibrationally-excited pure-rotational CO transition $v=1, J=3-2$ towards five asymptotic giant branch (AGB) stars, $o$ Cet, R Aqr, R Scl, W Aql, and $π^1$ Gru. The observed lines are formed in the poorly-understood region located between the stellar surface and the region where the wind starts, the so-called warm molecular layer. We successfully reproduce the observed lines profiles using a simple model. We constrain the extents, densities, and kinematics of the region where the lines are produced. R Aqr and R Scl show inverse P-Cygni line profiles which indicate infall of material onto the stars. The line profiles of $o$ Cet and R Scl show variability. The serendipitous detection towards these five sources shows that vibrationally-excited rotational lines can be observed towards a large number of nearby AGB stars using ALMA. This opens a new possibility for the study of the innermost regions of AGB circumstellar envelopes.

preprint2016arXiv

Dual differential polarimetry. A technique to recover polarimetric information from dual polarization observations

Current mm/submm interferometers, like the Atacama Large mm/submm Array (ALMA), use receivers that register the sky signal in a linear polarization basis. In the case of observations performed in full-polarization mode (where the cross-correlations are computed among all the polarization channels) it is possible to reconstruct the full-polarization brightness distribution of the observed sources, as long as a proper calibration of delay offsets and leakage among polarization channels can be performed. Observations of calibrators, preferably with some linear polarization, with a good parallactic angle coverage are usually needed for such a calibration. In principle, dual-polarization observations only allow us to recover the Stokes $I$ intensity distribution of the sources, regardless of the parallactic angle coverage of the observations. In this paper, we present a novel technique of dual differential polarimetry that makes it possible to obtain information related to the full-polarization brightness distribution of the observed sources from dual-polarization observations. This technique is inspired in the Earth-rotation polarization synthesis and can be applied even to sources with complex structures.

preprint2016arXiv

Measuring Magnetic Fields from Water Masers in the Synchrotron Protostellar Jet in W3(H$_2$O)

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.

preprint2016arXiv

Study of the inner dust envelope and stellar photosphere of the AGB star R Doradus using SPHERE/ZIMPOL

We use high-angular-resolution images obtained with SPHERE/ZIMPOL to study the photosphere, the warm molecular layer, and the inner wind of the close-by oxygen-rich AGB star R Doradus. We present observations in filters V, cntH$α$, and cnt820 and investigate the surface brightness distribution of the star and of the polarised light produced in the inner envelope. Thanks to second-epoch observations in cntH$α$, we are able to see variability on the stellar photosphere. We find that in the first epoch the surface brightness of R Dor is asymmetric in V and cntH$α$, the filters where molecular opacity is stronger, while in cnt820 the surface brightness is closer to being axisymmetric. The second-epoch observations in cntH$α$ show that the morphology of R Dor changes completely in a timespan of 48 days to a more axisymmetric and compact configuration. The polarised intensity is asymmetric in all epochs and varies by between a factor of 2.3 and 3.7 with azimuth for the different images. We fit the radial profile of the polarised intensity using a spherically symmetric model and a parametric description of the dust density profile, $ρ(r)=ρ_\circ r^{-n}$. On average, we find exponents of $- 4.5 \pm 0.5$ that correspond to a much steeper density profile than that of a wind expanding at constant velocity. The dust densities we derive imply an upper limit for the dust-to-gas ratio of $\sim 2\times10^{-4}$ at 5.0 $R_\star$. Given the uncertainties in observations and models, this value is consistent with the minimum values required by wind-driving models for the onset of a wind, of $\sim 3.3\times10^{-4}$. However, if the steep density profile we find extends to larger distances from the star, the dust-to-gas ratio will quickly become too small for the wind of R Dor to be driven by the grains that produce the scattered light.

preprint2016arXiv

The ALMA detection of CO rotational line emission in AGB stars in the Large Magellanic Cloud

Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB stars in the Large Magellanic Cloud. Methods: ALMA was used to observe 2 OH/IR stars and 4 carbon stars in the LMC in the CO J= 2-1 line. Results: We present the first measurement of expansion velocities in extragalactic carbon stars. All four C-stars are detected and wind expansion velocities and stellar velocities are directly measured. Mass-loss rates are derived from modelling the spectral energy distribution and Spitzer/IRS spectrum with the DUSTY code. Gas-to-dust ratios are derived that make the predicted velocities agree with the observed ones. The expansion velocities and MLRs are compared to a Galactic sample of well-studied relatively low MLRs stars supplemented with "extreme" C-stars that have properties more similar to the LMC targets. Gas MLRs derived from a simple formula are significantly smaller than derived from the dust modelling, indicating an order of magnitude underestimate of the estimated CO abundance, time-variable mass loss, or that the CO intensities in LMC stars are lower than predicted by the formula derived for Galactic objects. This could be related to a stronger interstellar radiation field in the LMC. Conclusions: Although the LMC sample is small and the comparison to Galactic stars is non-trivial because of uncertainties in their distances it appears that for C stars the wind expansion velocities in the LMC are lower than in the solar neighbourhood, while the MLRs appear similar. This is in agreement with dynamical dust-driven wind models.

preprint2015arXiv

A detailed view of the gas shell around R Sculptoris with ALMA

Thermal pulses are fundamental to the chemical evolution of AGB stars and their circumstellar envelopes. A further consequence of thermal pulses is the formation of detached shells of gas and dust around the star. We aim to determine the physical properties of the detached gas shell around R Sculptoris, in particular the shell mass and temperature, and to constrain the evolution of the mass-loss rate during and after a thermal pulse. We analyse CO(1-0), CO(2-1), and CO(3-2) emission, observed by. The spatial resolution of the ALMA data allows us to separate the detached shell emission from the extended emission inside the shell. We perform radiative transfer modelling of both components to determine the shell properties and the post-pulse mass-loss properties. The ALMA data show a gas shell with a radius of 19.5" expanding at 14.3km/s. The different scales probed by the ALMA Cycle 0 array show that the shell must be entirely filled with gas, contrary to the idea of a detached shell. The comparison to single-dish spectra and radiative transfer modelling confirms this. We derive a shell mass of 4.5e-3 Msun with a temperature of 50K. Typical timescales for thermal pulses imply a pulse mass-loss rate of 2.3e-5 Msun/yr. For the post-pulse mass-loss rate, we find evidence for a gradual decline of the mass-loss rate, with an average value of 1.6e-5 Msun/yr. The total amount of mass lost since the last thermal pulse is 0.03 Msun, a factor four higher compared to classical models, with a sharp decline in mass-loss rate immediately after the pulse. We find that the mass-loss rate after a thermal pulse has to decline more slowly than generally expected from models of thermal pulses. This may cause the star to lose significantly more mass during a thermal pulse cycle, which affects the chemical evolution of the star and the interstellar medium.

preprint2015arXiv

ALMA view of the circumstellar environment of the post-common-envelope-evolution binary system HD101584

We study the circumstellar evolution of the binary HD101584, consisting of a post-AGB star and a low-mass companion, which is most likely a post-common-envelope-evolution system. We used ALMA observations of the 12CO, 13CO, and C18O J=2-1 lines and the 1.3mm continuum to determine the morphology, kinematics, masses, and energetics of the circumstellar environment. The circumstellar medium has a bipolar hour-glass structure, seen almost pole-on, formed by an energetic jet, about 150 km/s. We conjecture that the circumstellar morphology is related to an event that took place about 500 year ago, possibly a capture event where the companion spiraled in towards the AGB star. However, the kinetic energy of the accelerated gas exceeds the released orbital energy, and, taking into account the expected energy transfer efficiency of the process, the observed phenomenon does not match current common-envelope scenarios. This suggests that another process must augment, or even dominate, the ejection process. A significant amount of material resides in an unresolved region, presumably in the equatorial plane of the binary system.

preprint2015arXiv

EVN observations of 6.7 GHz methanol maser polarization in massive star-forming regions III. The flux-limited sample

Theoretical simulations and observations at different angular resolutions have shown that magnetic fields have a central role in massive star formation. Like in low-mass star formation, the magnetic field in massive young stellar objects can either be oriented along the outflow axis or randomly. Measuring the magnetic field at milliarcsecond resolution (10-100 au) around a substantial number of massive young stellar objects permits determining with a high statistical significance whether the direction of the magnetic field is correlated with the orientation of the outflow axis or not. In late 2012, we started a large VLBI campaign with the European VLBI Network to measure the linearly and circularly polarized emission of 6.7 GHz methanol masers around a sample of massive star-forming regions. This paper focuses on the first seven observed sources, G24.78+0.08, G25.65+1.05, G29.86-0.04, G35.03+0.35, G37.43+1.51, G174.20-0.08, and G213.70-12.6. For all these sources, molecular outflows have been detected in the past. We detected a total of 176 methanol masing cloudlets toward the seven massive star-forming regions, 19% of which show linearly polarized emission. The methanol masers around the massive young stellar object MM1 in G174.20-0.08 show neither linearly nor circularly polarized emission. The linear polarization vectors are well ordered in all the other massive young stellar objects. We measured significant Zeeman splitting toward both A1 and A2 in G24.78+0.08, and toward G29.86-0.04 and G213.70-12.6. By considering all the 19 massive young stellar objects reported in the literature for which both the orientation of the magnetic field at milliarcsecond resolution and the orientation of outflow axes are known, we find evidence that the magnetic field (on scales 10-100 au) is preferentially oriented along the outflow axes.

preprint2015arXiv

Magnetic field measurements at milliarcsecond resolution around massive young stellar objects

Magnetic fields have only recently been included in theoretical simulations of high-mass star formation. The simulations show that magnetic fields can play a crucial role not only in the formation and dynamics of molecular outflows, but also in the evolution of circumstellar disks. Therefore, new measurements of magnetic fields at milliarcsecond resolution close to massive young stellar objects (YSOs) are fundamental for providing new input for numerical simulations and for understanding the formation process of massive stars. The polarized emission of 6.7 GHz CH3OH masers allows us to investigate the magnetic field close to the massive YSO where the outflows and disks are formed. Recently, we have detected with the EVN CH3OH maser polarized emission towards 10 massive YSOs. From a first statistical analysis we have found evidence that magnetic fields are primarily oriented along the molecular outflows. To improve our statistics we are carrying on a large observational EVN campaign for a total of 19 sources, the preliminary results of the first seven sources are presented in this contribution. Furthermore, we also describe our efforts to estimate the Lande' g-factors of the CH3OH maser transition to determine the magnetic field strength from our Zeeman-splitting measurements.

preprint2015arXiv

Observing the onset of outflow collimation in a massive protostar

The current paradigm of star formation through accretion disks, and magnetohydrodynamically driven gas ejections, predicts the development of collimated outflows, rather than expansion without any preferential direction. We present radio continuum observations of the massive protostar W75N(B)-VLA 2, showing that it is a thermal, collimated ionized wind and that it has evolved in 18 years from a compact source into an elongated one. This is consistent with the evolution of the associated expanding water-vapor maser shell, which changed from a nearly circular morphology, tracing an almost isotropic outflow, to an elliptical one outlining collimated motions. We model this behavior in terms of an episodic, short-lived, originally isotropic, ionized wind whose morphology evolves as it moves within a toroidal density stratification.

preprint2015arXiv

Resolving the stellar activity of the Mira AB binary with ALMA

We present the size, shape and flux densities at millimeter continuum wavelengths, based on ALMA science verification observations in Band 3 (~94.6 GHz) and Band 6 (~228.7 GHz), from the binary Mira A (o Ceti) and Mira B. The Mira AB system has been observed with ALMA at a spatial resolution of down to ~25 mas. The extended atmosphere of Mira A and the wind around Mira B sources are resolved and we derive the size of Mira A and of the ionized region around Mira B. The spectral indices within Band 3 (between 89-100 GHz) and between Band 3 and Band 6 are also derived. The spectral index of Mira A is found to change from 1.71+-0.05 within Band 3 to 1.54+-0.04 between Band 3 and 6. The spectral index of Mira B is 1.3+-0.2 in Band 3, in good agreement with measurements at longer wavelengths. However it rises to 1.72+-0.11 between the bands. For the first time the extended atmosphere of a star is resolved at these frequencies and for Mira A the diameter is ~3.8x3.2 AU in Band 3 (with brightness temperature Tb~5300 K) and ~4.0x3.6 AU in Band 6 (Tb~2500 K). Additionally, a bright hotspot of ~0.4 AU and with Tb~10000 K is found on the stellar disc of Mira A. The size of the ionized region around the accretion disk of Mira B is found to be ~2.4 AU. The emission around Mira B is consistent with that from a partially ionized wind of gravitationally bound material from Mira A close to the accretion disk of Mira B. The Mira A atmosphere does not fully match predictions, with brightness temperatures in Band 3 significantly higher than expected, potentially due to shock heating. The hotspot is likely due to magnetic activity and could be related to the previously observed X-ray flare of Mira A.

preprint2015arXiv

Velocity and magnetic fields within 1000 AU from a massive YSO

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.

preprint2014arXiv

From the ashes: JVLA observations of water fountain nebula candidates show the rebirth of IRAS 18455+0448

[abridged] The class of water fountain nebulae is thought to represent the stage of the earliest onset of collimated bipolar outflows during the post-Asymptotic Giant Branch phase. They thus play a crucial role in the study of the formation of bipolar Planetary Nebulae (PNe). To date, 14 water fountain nebulae have been identified. The identification of more sources in this unique stage of stellar evolution will enable us to study the origin of bipolar PNe morphologies in more detail. We present the results of seven sources observed with the JVLA that were identified as water fountain candidates in an Effelsberg 100m telescope survey of 74 AGB and early post-AGB stars. We find that our sample of water fountain candidates displays strong variability in their 22 GHz H2O maser spectra. The JVLA observations show an extended bipolar H2O maser outflow for one source, the OH/IR star IRAS 18455+0448. This source was previously classified as a 'dying' OH/IR star based on the exponential decrease of its 1612 MHz OH maser and the lack of H2O masers. We therefore also re-observed the 1612, 1665, and 1667 MHz OH masers. We confirm that the 1612 MHz masers have not reappeared and find that the 1665/1667 MHz masers have decreased in strength by several orders of magnitude during the last decade. The OH/IR star IRAS 18455+0448 is confirmed to be a new addition to the class of water fountain nebulae. Its kinematic age is approximately 70 yr, but could be lower, depending on the distance and inclination. Previous observations indicate, with significant uncertainty, that IRAS 18455+0448 has a surprisingly low mass compared to available estimates for other water fountain nebulae. The available historical OH maser observations make IRAS 18455+0448 unique for the study of water fountain nebulae and the launch of post-AGB bipolar outflows...

preprint2014arXiv

Rapidly increasing collimation and magnetic field changes of a H2O maser outflow

W75N(B) is a massive star-forming region that contains three radio continuum sources (VLA 1, VLA 2, and VLA 3), which are thought to be three massive young stellar objects at three different evolutionary stages. VLA 1 is the most evolved and VLA 2 the least evolved source. The 22 GHz H2O masers associated with VLA 1 and VLA 2 have been mapped at several epochs over eight years. While the H2O masers in VLA 1 show a persistent linear distribution along a radio jet, those in VLA 2 are distributed around an expanding shell. Furthermore, H2O maser polarimetric measurements revealed magnetic fields aligned with the two structures. Using new polarimetric observations of H2O masers, we aim to confirm the elliptical expansion of the shell-like structure around VLA 2 and, at the same time, to determine if the magnetic fields around the two sources have changed. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22 GHz H2O masers towards the massive star-forming region W75N(B). The H2O maser distribution around VLA 1 is unchanged from that previously observed. We made an elliptical fit of the H2O masers around VLA 2. We find that the shell-like structure is still expanding along the direction parallel to the thermal radio jet of VLA 1. While the magnetic field around VLA 1 has not changed in the past 7 years, the magnetic field around VLA 2 has changed its orientation according to the new direction of the major-axis of the shell-like structure and it is now aligned with the magnetic field in VLA 1.

preprint2014arXiv

Submillimeter H2O masers in water-fountain nebulae

We report the first detection of submillimeter water maser emission toward water-fountain nebulae, which are post-AGB stars that exhibit high-velocity water masers. Using APEX we found emission in the ortho-H2O (10_29-9_36) transition at 321.226 GHz toward three sources: IRAS 15445-5449, IRAS 18043-2116 and IRAS 18286-0959. Similarly to the 22 GHz masers, the submillimeter water masers are expanding with a velocity larger than that of the OH masers, suggesting that these masers also originate in fast bipolar outflows. In IRAS 18043-2116 and IRAS 18286-0959, which figure among the sources with the fastest water masers, the velocity range of the 321 GHz masers coincides with that of the 22 GHz masers, indicating that they likely coexist. Towards IRAS 15445-5449 the submillimeter masers appear in a different velocity range, indicating that they are tracing different regions. The intensity of the submillimeter masers is comparable to that of the 22 GHz masers, implying that the kinetic temperature of the region where the masers originate should be Tk > 1000 K. We propose that the passage of two shocks through the same gas can create the conditions necessary to explain the presence of strong high-velocity 321 GHz masers coexisting with the 22 GHz masers in the same region.

preprint2014arXiv

The magnetic field at milliarcsecond resolution around IRAS20126+4104

IRAS20126+4104 is a well studied B0.5 protostar that is surrounded by a ~1000 au Keplerian disk and is where a large-scale outflow originates. Both 6.7-GHz CH3OH masers and 22-GHz H2O masers have been detected toward this young stellar object. The CH3OH masers trace the Keplerian disk, while the H2O masers are associated with the surface of the conical jet. Recently, observations of dust polarized emission (350 um) at an angular resolution of 9 arcseconds (~15000 au) have revealed an S-shaped morphology of the magnetic field around IRAS20126+4104. The observations of polarized maser emissions at milliarcsecond resolution (~20 au) can make a crucial contribution to understanding the orientation of the magnetic field close to IRAS20126+4104. This will allow us to determine whether the magnetic field morphology changes from arcsecond resolution to milliarcsecond resolution. The European VLBI Network was used to measure the linear polarization and the Zeeman splitting of the 6.7-GHz CH3OH masers toward IRAS20126+4104. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman splitting of the 22-GHz H2O masers toward the same region. We detected 26 CH3OH masers and 5 H2O masers at high angular resolution. Linear polarization emission was observed toward three CH3OH masers and toward one H2O maser. Significant Zeeman splitting was measured in one CH3OH maser (ΔV_{Z}=-9.2 +/- 1.4 m/s). No significant (5 sigma) magnetic field strength was measured using the H2O masers. We found that in IRAS20126+4104 the rotational energy is less than the magnetic energy.

preprint2014arXiv

The wonderful complexity of the Mira AB system

We have mapped the CO(3-2) line emission around the Mira AB system at 0.5 resolution using the Atacama Large Millimeter/submillimeter Array (ALMA). The CO map shows amazing complexity. The circumstellar gas has been shaped by different dynamical actors during the evolution of the system and several morphological components can be identified. The companion is marginally resolved in continuum emission and is currently at 0.487$\pm$0.006 separation. In the main line component, centered on the stellar velocity, spiral arcs around Mira A are found. The spiral appears to be relatively flat and oriented in the orbital plane. An accretion wake behind the companion is clearly visible and the projected arc separation is of order 5''. In the blue wing of the line emission, offset from the main line, several large ($\sim$5-10''), opposing arcs are found. We tentatively suggest that this structure is created by the wind of Mira B blowing a bubble in the expanding envelope of Mira A.

preprint2014arXiv

UVMULTIFIT: A versatile tool for fitting astronomical radio interferometric data

The analysis of astronomical interferometric data is often performed on the images obtained after deconvolution of the interferometer's point spread function (PSF). This strategy can be understood (especially for cases of sparse arrays) as fitting models to models, since the deconvolved images are already non-unique model representations of the actual data (i.e., the visibilities). Indeed, the interferometric images may be affected by visibility gridding, weighting schemes (e.g., natural vs. uniform), and the particulars of the (non-linear) deconvolution algorithms. Fitting models to the direct interferometric observables (i.e., the visibilities) is preferable in the cases of simple (analytical) sky intensity distributions. In this paper, we present UVMULTIFIT, a versatile library for fitting visibility data, implemented in a Python-based framework. Our software is currently based on the CASA package, but can be easily adapted to other analysis packages, provided they have a Python API. We have tested the software with synthetic data, as well as with real observations. In some cases (e.g., sources with sizes smaller than the diffraction limit of the interferometer), the results from the fit to the visibilities (e.g., spectra of close by sources) are far superior to the output obtained from the mere analysis of the deconvolved images. UVMULTIFIT is a powerful improvement of existing tasks to extract the maximum amount of information from visibility data, especially in cases close to the sensitivity/resolution limits of interferometric observations.

preprint2013arXiv

A synchrotron jet from a post-AGB star

(Abridged) The evolution of low- and intermediate-initial-mass stars beyond the asymptotic giant branch (AGB) remains poorly understood. High-velocity outflows launched shortly after the AGB phase are thought to be the primary shaping mechanism of bipolar and multipolar planetary nebulae (PNe). However, little is known about the launching and driving mechanism for these jets, whose momentum and energy often far exceed the energy that can be provided by radiation pressure alone. Here, we report direct evidence of a magnetically collimated jet shaping the bipolar morphology of the circumstellar envelope of a post-AGB star. We present radio continuum observations of the post-AGB star IRAS 15445-5449 which has water masers tracing a fast bipolar outflow. Our observations confirm the earlier observed steep negative spectral index of the spectral energy distribution (SED) above 3 GHz, and resolve, for the first time, the emission to originate from a radio jet, proving the existence of such jets around a post-AGB star. The SED is consistent with a synchrotron jet embedded in a sheath of thermal electrons. We find a close correspondence between the extent and direction of the synchrotron jet and the bipolar shape of the object observed at other wavelenghts, suggesting that the jet is responsible for the source morphology. The jet is collimated by a magnetic field of the order of mG at almost 7000 AU from the central star. We recover observations from the ATCA archive that indicate that the emission measure of the thermal component has increased by a factor of three between 1998 and 2005 after which it has remained constant. The short timescale evolution of the radio emission suggests a short lifetime for the jet.

preprint2013arXiv

ALMA Observations of the Coldest Place in the Universe: The Boomerang Nebula

The Boomerang Nebula is the coldest known object in the Universe, and an extreme member of the class of Pre-Planetary Nebulae, objects which represent a short-lived transitional phase between the AGB and Planetary Nebula evolutionary stages. Previous single-dish CO (J=1-0) observations (with a 45 arcsec beam) showed that the high-speed outflow in this object has cooled to a temperature significantly below the temperature of the cosmic background radiation. Here we report the first observations of the Boomerang with ALMA in the CO J=2-1 and J=1-0 lines to resolve the structure of this ultra-cold nebula. We find a central hourglass-shaped nebula surrounded by a patchy, but roughly round, cold high-velocity outflow. We compare the ALMA data with visible-light images obtained with HST and confirm that the limb-brightened bipolar lobes seen in these data represent hollow cavities with dense walls of molecular gas and dust producing both the molecular-emission-line and scattered-light structures seen at millimeter and visible wavelengths. The large diffuse biconical shape of the nebula seen in the visible is likely due to preferential illumination of the cold high-velocity outflow. We find a compact source of millimeter-wave continuum in the nebular waist -- these data, together with sensitive upper limits on the radio continuum using observations with ATCA, indicate the presence of a substantial mass of very large (mm-sized) grains in the waist of the nebula. Another unanticipated result is the detection of CO emission regions beyond the ultracold region which indicate the re-warming of the cold gas, most likely due to photoelectric grain heating.

preprint2013arXiv

ALMA observations of the variable 12CO/13CO ratio around the asymptotic giant branch star R Sculptoris

[abridged] The 12CO/13CO ratio is often used as a measure of the 12C/13C ratio in the circumstellar environment, carrying important information about the stellar nucleosynthesis. External processes can change the 12CO and 13CO abundances, and spatially resolved studies of the 12CO/13CO ratio are needed to quantify the effect of these processes on the globally determined values. Additionally, such studies provide important information on the conditions in the circumstellar environment. The detached-shell source R Scl, displaying CO emission from recent mass loss, in a binary-induced spiral structure as well as in a clumpy shell produced during a thermal pulse, provides a unique laboratory for studying the differences in CO isotope abundances throughout its recent evolution. We observed both the 12CO(J=3-2) and the 13CO(J=3-2) line using ALMA. We find significant variations in the 12CO/13CO intensity ratios and consequently in the abundance ratios. The average CO isotope abundance ratio is at least a factor three lower in the shell (~19) than that in the present-day (<300 years) mass loss (>60). Additionally, variations in the ratio of more than an order of magnitude are found in the shell itself. We attribute these variations to the competition between selective dissociation and isotope fractionation in the shell, of which large parts cannot be warmer than ~35 K. However, we also find that the 12CO/13CO ratio in the present-day mass loss is significantly higher than the 12C/13C ratio determined in the stellar photosphere from molecular tracers (~19). The origin of this discrepancy is still unclear, but we speculate that it is due to an embedded source of UV-radiation that is primarily photo-dissociating 13CO. This radiation source could be the hitherto hidden companion. Alternatively, the UV-radiation could originate from an active chromosphere of R Scl itself....

preprint2013arXiv

EVN observations of 6.7-GHz methanol maser polarization in massive star-forming regions II. First statistical results

Magnetic fields have only recently been included in theoretical simulations of high-mass star formation. The simulations show that magnetic fields play an important role in the formation and dynamics of molecular outflows. Masers, in particular 6.7-GHz CH3OH masers, are the best probes of the magnetic field morphologies around massive young stellar objects on the smallest scales of 10-100 AU. This paper focuses on 4 massive young stellar objects, IRAS06058+2138-NIRS1, IRAS22272+6358A, S255-IR, and S231, which complement our previous 2012 sample (the first EVN group). From all these sources, molecular outflows have been detected in the past. Seven of the European VLBI Network antennas were used to measure the linear polarization and Zeeman-splitting of the 6.7-GHz CH3OH masers in the star-forming regions in this second EVN group. We detected a total of 128 CH3OH masing cloudlets. Fractional linear polarization (0.8%-11.3%) was detected towards 18% of the CH3OH masers in our sample. The linear polarization vectors are well ordered in all the massive young stellar objects. We measured significant Zeeman-splitting in IRAS06058+2138-NIRS1 (DVz=3.8+/-0.6 m/s) and S255-IR (DVz=3.2+/-0.7 m/s). By considering the 20 massive young stellar objects towards which the morphology of magnetic fields was determined by observing 6.7-GHz CH3OH masers in both hemispheres, we find no evident correlation between the linear distributions of CH3OH masers and the outflows or the linear polarization vectors. On the other hand, we present first statistical evidence that the magnetic field (on scales 10-100 AU) is primarily oriented along the large-scale outflow direction. Moreover, we empirically find that the linear polarization fraction of unsaturated CH3OH masers is P_l<4.5%.

preprint2013arXiv

Linear polarization of submillimetre masers. Tracing magnetic fields with ALMA

(Abridged) The goal of this work is to investigate if SiO, H2O and HCN maser emission within the ALMA frequency range can be detected with observable levels of fractional linear polarization in the regime where the Zeeman frequency (gOmega) is greater than the stimulated emission rate (R) and the decay rate of the molecular state (Gamma). We have used a radiative transfer code to calculate the fractional linear polarization as a function of the emerging brightness temperature for a number of rotational transition of SiO, H2O and HCN which have been observed to display maser emission at submillimetre wavelengths. We assume typical magnetic field strengths measured towards galactic star forming regions and circumstellar envelopes of late-type stars from previous VLBI observations. Since the Landé g-factors have not been reported for the different rotational transitions we have modeled, we performed our calculations assuming conservative values of the Zeeman frequency for the different molecular species. Setting a lower limit for the Zeeman frequency which still satisfies the criteria gOmega>R and gOmega>Gamma, we find fractional polarization levels of up to 13%, 14% and 19% for the higher J transitions analysed for SiO, H2O and HCN, respectively, without considering anisotropic pumping or any other non-Zeeman effect. Such upper limits were calculated assuming a magnetic field oriented perpendicular to the direction of propagation of the maser radiation. According to our results SiO, H2O and HCN maser emission within the ALMA frequency range can be detected with suitable linear polarization to trace the magnetic field structure towards star forming regions and late-type stars, even if the detected polarization has been enhanced by non-Zeeman effects.

preprint2013arXiv

Magnetic fields around evolved stars: further observations of H$_2$O maser polarization

We aim to detect the magnetic field and infer its properties around four AGB stars using H$_2$O maser observations. The sample we observed consists of the following sources: the semi-regular variable RT Vir and the Mira variables AP Lyn, IK Tau, and IRC+60370. We observed the 6$_{1,6}-5_{2,3}$ H$_2$O maser rotational transition, in full-polarization mode, to determine its linear and circular polarization. Based on the Zeeman effect, one can infer the properties of the magnetic field from the maser polarization analysis. We detected a total of 238 maser features, in three of the four observed sources. No masers were found toward AP Lyn. The observed masers are all located between 2.4 and 53.0 AU from the stars. Linear and circular polarization was found in 18 and 11 maser features, respectively. We more than doubled the number of AGB stars in which magnetic field has been detected from H$_2$O maser polarization, as our results confirm the presence of fields around IK Tau, RT Vir and IRC+60370. The strength of the field along the line of sight is found to be between 47 and 331 mG in the H$_2$O maser region. Extrapolating this result to the surface of the stars, assuming a toroidal field ($\propto$ r$^{-1}$), we find magnetic fields of 0.3-6.9 G on the stellar surfaces. If, instead of a toroidal field, we assume a poloidal field ($\propto$ r$^{-2}$), then the extrapolated magnetic field strength on the stellar surfaces are in the range between 2.2 and $\sim$115 G. Finally, if a dipole field ($\propto$ r$^{-3}$) is assumed, the field strength on the surface of the star is found to be between 15.8 and $\sim$1945 G. The magnetic energy of our sources is higher than the thermal and kinetic energy in the H$_2$O maser region of this class of objects. This leads us to conclude that, indeed, magnetic fields probably play an important role in shaping the outflows of evolved stars. (abridged)

preprint2013arXiv

Rapid Angular Expansion of the Ionized Core of CRL 618

CRL 618 is an object that exhibits characteristics of both AGB and post-AGB star. It also displays a spectacular array of bipolar lobes with a dense equatorial region, which makes it an excellent object to study the development of asymmetries in evolved stars. In the recent decades, an elliptical compact HII region located in the center of the nebula has been seen to be increasing in size and flux. This seems to be due to the ionization of the circumstellar envelope by the central star, and it would be indicating the beginning of the planetary nebula phase for CRL 618. We analyzed interferometric radio continuum data at ~5 and 22 GHz from observations carried out at seven epochs with the VLA. We traced the increase of the flux of the ionized region over a period of ~26 years. We measured the dimensions of the HII region directly from the brightness distribution images to determine the increase of its size over time. For one of the epochs we analyzed observations at six frequencies from which we estimated the electron density distribution. We carried out model calculations of the spectral energy distribution at two different epochs to corroborate our observational results. We found that the radio continuum flux and the size of the ionized region have been increasing monotonically in the last three decades. The size of the major axis of the HII region shows a dependance with frequency, which has been interpreted as a result of a gradient of the electron density in this direction. The growth of the HII region is due to the expansion of an ionized wind whose mass-loss rate increased continuously for a period of ~100 years until a few decades ago, when the mass-loss rate experienced a sudden decline. Our results indicate that the beginning of the ionization of the circumstellar envelope began around 1971, which marks the start of the planetary nebula phase of CRL 618.

preprint2012arXiv

Adaptable Radiative Transfer Innovations for Submillimetre Telescopes (ARTIST) - Dust polarisation module (DustPol)

We present a new publicly available tool (DustPol) aimed to model the polarised thermal dust emission. The module DustPol, which is publicly available, is part of the ARTIST (Adaptable Radiative Transfer Innovations for Submillimetre Telescopes) package, which also offers tools for modelling the polarisation of line emission together with a model library and a Python-based user interface. DustPol can easily manage analytical as well as pre-gridded models to generate synthetic maps of the Stokes I, Q, and U parameters. These maps are stored in FITS format which is straightforwardly read by the data reduction software used, e.g., by the Atacama Large Millimeter Array (ALMA). This turns DustPol into a powerful engine for the prediction of the expected polarisation features of a source observed with ALMA or the Planck satellite as well as for the interpretation of existing submillimetre observations obtained with other telescopes. DustPol allows the parameterisation of the maximum degree of polarisation and we find that, in a prestellar core, if there is depolarisation, this effect should happen at densities of 10^6 cm-3 or larger. We compare a model generated by DustPol with the observational polarisation data of the low-mass Class 0 object NGC 1333 IRAS 4A, finding that the total and the polarised emission are consistent.

preprint2012arXiv

EVN observations of 6.7 GHz methanol maser polarization in massive star-forming regions

The role of magnetic fields in the formation of high-mass stars is still under debate, and recent measurements of their orientation and strength by using polarized maser emissions are contributing new insights. Masers polarization, in particular of the 6.7-GHz methanol masers, are one of the best probes of the magnetic field morphologies around massive protostars. Determining the magnetic field morphology around an increasing number of massive protostars at milliarcsecond resolution by observing 6.7-GHz methanol masers is crucial to better understand the role of magnetic fields in massive star formation.The First EVN Group consists of 4 massive star-forming complexes: W51, W48, IRAS18556+0138, and W3(OH). These contain well-studied \hii ~regions from some of which molecular bipolar outflows were also detected (W51-e2, G35.20-0.74N). Nine of the European VLBI Network antennas were used to measure the linear polarization and Zeeman-splitting of the 6.7-GHz methanol masers in the star-forming regions of the First EVN Group. We detected a total of 154 CH3OH masers, one third of these towards W3(OH). Fractional linear polarization (1.2-11.5%) was detected towards 55 masers. The linear polarization vectors are well-ordered in all the massive star-forming regions. We measured significant Zeeman-splitting in 3 massive star-forming regions (W51, W48, and W3(OH)) revealing a range of separations -3.5 m/s<ΔV_{z}<3.8 m/s with the smallest |ΔV_{z}|=0.4m/s. We were also able to compare our magnetic field results with those obtained from submillimeter wavelength dust observation in W51 and show that the magnetic field at low and high resolutions are in perfect agreement.

preprint2012arXiv

No evidence for a central IMBH in M15

Intermediate mass black holes (IMBHs) with expected masses M_BH ~ 10^4 M_sun are thought to bridge the gap between stellar mass black holes (M_BH ~ 3 - 100 M_sun) and supermassive black holes found at the centre of galaxies (M_BH > 10^6 M_sun). Until today, no IMBH has been confirmed observationally. The most promising objects to host an IMBH as their central mass are globular clusters. Here, we present high sensitivity multi-epoch 1.6 GHz very long baseline interferometry observations of the globular cluster M15 that has been suggested to host an IMBH. Assuming the IMBH to be accreting matter from its surrounding we expect to detect it as a point source moving with the global motion of the cluster. However, we do not detect any such object within a radius of 6000 AU of the cluster centre in any of the five observations spread over more than one year. This rules out any variability of the putative IMBH on the time scale of one to two months. To get the most stringent upper limit for the flux density of the putative IMBH we concatenate the data of all five epochs. In this data we measure a 3σ upper flux limit of 10 μJy for a central source. We employ the fundamental plane of black hole activity to estimate the mass of the central IMBH candidate. Based on previous X-ray observations of M15 our measurements indicate a 3σ upper mass limit of ~500 M_sun.

preprint2012arXiv

Polarization of thermal molecular lines in the envelope of IK Tau

Molecular line polarization is a unique source of information about the magnetic fields and anisotropies in the circumstellar envelopes of evolved stars. Here we present the first detection of thermal CO(J=2-1) and SiO(J=5-4, v=0) polarization, in the envelope of the asymptotic giant branch star IK Tau. The observed polarization direction does not match predictions for circumstellar envelope polarization induced only by an anisotropic radiation field. Assuming that the polarization is purely due to the Goldreich-Kylafis effect, the linear polarization direction is defined by the magnetic field as even the small Zeeman splitting of CO and SiO dominates the molecular collisional and spontaneous emission rates. The polarization was mapped using the Submillimeter Array (SMA) and is predominantly north-south. There is close agreement between the CO and SiO observations, even though the CO polarization arises in the circumstellar envelope at ~800 AU and the SiO polarization at <250 AU. If the polarization indeed traces the magnetic field, we can thus conclude that it maintains a large-scale structure throughout the circumstellar envelope. We propose that the magnetic field, oriented either east-west or north-south is responsible for the east-west elongation of the CO distribution and asymmetries in the dust envelope. In the future, the Atacama Large Millimeter/submillimeter Array will be able to map the magnetic field using CO polarization for a large number of evolved stars.

preprint2012arXiv

Rotten Egg Nebula: The magnetic field of a binary evolved star

Most of PNe are not spherical. The loss of spherical symmetry occurs somewhere between the AGB and PN phase. The cause of this change of morphology is not yet well known, but magnetic fields are one of the possible agents. Its origin remains to be determined, and potentially requires the presence of a massive companion to the AGB star. Therefore, further detections of the magnetic field around evolved stars (in particular those thought to be part of a binary system) are crucial to improve our understanding of the origin and role of magnetism on evolved stars. One such binaries is the pre-PN OH231.8, around which a magnetic field was detected in the OH maser region of the outer circumstellar envelope. We aim to detect and infer the properties of the magnetic field of this source in the water maser region. We observed the 6_{1,6}-5_{2,3} water maser rotational transition to determine its linear and circular polarization. These emissions are located within the inner regions of OH231.8 (at few tens of AU). We detected 30 water maser features, which occur in two distinct regions that are moving apart with a velocity on the sky of 2.3 mas/year. Taking into account the inclination angle of the source, this corresponds to an average separation velocity of 21 km/s. Based on the velocity gradient of the maser emission, the masers appear to be dragged along the direction of the nebula jet. Linear polarization is present in 3 of the features, and circular polarization was detected in the 2 brightest ones. We found that the strength of the magnetic field is |B_{||}|~45 mG which, when assuming a toroidal magnetic field, implies B~2.5 G on the stellar surface. The morphology of the field is not yet determined, but the high scatter found on the directions of the linear polarization vectors could indicate that the masers occur near the tangent points of a toroidal field.

preprint2012arXiv

Searching for X-ray emission from AGB stars

Magnetic fields have been measured around Asymptotic Giant Branch (AGB) stars of all chemical types using maser polarization observations. If present, a large-scale magnetic field would lead to X-ray emission, which should be observable using current X-ray observatories. The aim of this work is to search the archival data for AGB stars that are intrinsic X-ray emitters. We have searched the ROSAT, CXO, and XMM archives for serendipitous X-ray observations of a sample of ~500 AGB stars. We specifically searched for the AGB stars detected with GALEX. The data is calibrated, analyzed and the X-ray luminosities and temperatures are estimated as functions of the circumstellar absorption. We identify 13 AGB stars as having either serendipitous or targeted observations in the X-ray data archives, however for a majority of the sources the detailed analysis show that the detections are questionable. Two new sources are detected by ROSAT: T Dra and R UMa. The spectral analysis suggests that the emission associated with these sources could be due to coronal activity or interaction across a binary system. Further observations of the detected sources are necessary to clearly determine the origin of the X-ray emission. Moreover, additional objects should be subject to targeted X-ray observations in order to achieve better constraints for the magnetic fields around AGB stars.

preprint2012arXiv

Unexpectedly large mass loss during the thermal pulse cycle of the red giant R Sculptoris!

The asymptotic giant branch star R Sculptoris is surrounded by a detached shell of dust and gas. The shell originates from a thermal pulse during which the star undergoes a brief period of increased mass loss. It has hitherto been impossible to constrain observationally the timescales and mass-loss properties during and after a thermal pulse - parameters that determine the lifetime on the asymptotic giant branch and the amount of elements returned by the star. Here we report observations of CO emission from the circumstellar envelope and shell around R Sculptoris with an angular resolution of 1.3 arcsec. What was hitherto thought to be only a thin, spherical shell with a clumpy structure, is revealed to contain a spiral structure. Spiral structures associated with circumstellar envelopes have been seen previously, from which it was concluded that the systems must be binaries. Using the data, combined with hydrodynamic simulations, we conclude that R Sculptoris is a binary system that underwent a thermal pulse approximately 1800 years ago, lasting approximately 200 years. About 0.003 Msun of mass was ejected at a velocity of 14.3 km s-1 and at a rate approximately 30 times higher than the prepulse mass-loss rate. This shows that approximately 3 times more mass is returned to the interstellar medium during and immediately after a pulse than previously thought.

preprint2011arXiv

H2O Maser Polarization of the Water Fountains IRAS 15445-5449 and IRAS 18043-2116

We present the morphology and linear polarization of the 22-GHz H2O masers in the high-velocity outflow of two post-AGB sources, d46 (IRAS 15445-5449) and b292 (IRAS 18043-2116). The observations were performed using The Australia Telescope Compact Array. Different levels of saturated maser emission have been detected for both sources. We also present the mid-infrared image of d46 overlaid with the distribution of the maser features that we have observed in the red-shifted lobe of the bipolar structure. The relative position of the observed masers and a previous radio continuum observation suggests that the continnum is produced along the blue-shifted lobe of the jet. It is likely due to synchrontron radiation, implying the presence of a strong magnetic field in the jet. The fractional polarization levels measured for the maser features of d46 indicate that the polarization vectors are tracing the poloidal component of the magnetic field in the emitting region. For the H2O masers of b292 we have measured low levels of fractional linear polarization. The linear polarization in the H2O maser region of this source likely indicates a dominant toroidal or poloidal magnetic field component. Since circular polarization was not detected it is not possible to determine the magnetic field strength. However, we present a 3-sigma evaluation of the upper limit intensity of the magnetic field in the maser emitting regions of both observed sources.

preprint2011arXiv

Magnetic Fields in Evolved Stars: Imaging the Polarized Emission of High-Frequency SiO Masers

We present Submillimeter Array observations of high frequency SiO masers around the supergiant VX Sgr and the semi-regular variable star W Hya. The J=5-4, v=1 28SiO and v=0 29SiO masers of VX Sgr are shown to be highly linearly polarized with a polarization from ~5-60%. Assuming the continuum emission peaks at the stellar position, the masers are found within ~60 mas of the star, corresponding to ~100 AU at a distance of 1.57 kpc. The linear polarization vectors are consistent with a large scale magnetic field, with position and inclination angles similar to that of the dipole magnetic field inferred in the H2O and OH maser regions at much larger distances from the star. We thus show for the first time that the magnetic field structure in a circumstellar envelope can remain stable from a few stellar radii out to ~1400 AU. This provides further evidence supporting the existence of large scale and dynamically important magnetic fields around evolved stars. Due to a lack of parallactic angle coverage, the linear polarization of masers around W Hya could not be determined. For both stars we observed the 28SiO and 29SiO isotopologues and find that they have a markedly different distribution and that they appear to avoid each other. Additionally, emission from the SO 5_5-4_4 line was imaged for both sources. Around W Hya we find a clear offset between the red- and blue-shifted SO emission. This indicates that W Hya is likely host to a slow bipolar outflow or a rotating disk-like structure.

preprint2011arXiv

The properties and polarization of the H2O and CH3OH maser environment of NGC7538-IRS1

NGC7538 is a complex massive star-forming region. The region is composed of several radio continuum sources, one of which is IRS1, a high-mass protostar, from which a 0.3 pc molecular bipolar outflow was detected. Several maser species have been detected around IRS1. The CH3OH masers have been suggested to trace a Keplerian-disk, while the H2O masers are almost aligned to the outflow. More recent results suggested that the region hosts a torus and potentially a disk, but with a different inclination than the Keplerian-disk that is supposed to be traced by the CH3OH masers. Tracing the magnetic field close to protostars is fundamental for determining the orientation of the disk/torus. Recent studies showed that during the protostellar phase of high-mass star formation the magnetic field is oriented along the outflows and around or on the surfaces of the disk/torus. The observations of polarized maser emissions at milliarcsecond resolution can make a crucial contribution to understanding the orientation of the magnetic field and, consequently, the orientation of the disk/torus in NGC7538-IRS1. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22GHz H2O masers toward NGC7538-IRS1. The European VLBI Network and the MERLIN telescopes were used to measure the linear polarization and the Zeeman-splitting of the 6.7GHz CH3OH masers toward the same region. We detected 17 H2O masers and 49 CH3OH masers at high angular resolution. We detected linear polarization emission toward two H2O masers and toward twenty CH3OH masers. The CH3OH masers, most of which only show a core structure, seem to trace rotating and potentially infalling gas in the inner part of a torus. Significant Zeeman-splitting was measured in three CH3OH masers. [...] We also propose a new description of the structure of the NGC7538-IRS1 maser region.

preprint2011arXiv

The structure of the magnetic field in the massive star-forming region W75N

A debated topic in star formation theory is the role of magnetic fields during the protostellar phase of high-mass stars. It is still unclear how magnetic fields influence the formation and dynamics of massive disks and outflows. Most current information on magnetic fields close to high-mass protostars comes from polarized maser emissions, which allows us to investigate the magnetic field on small scales by using very long-baseline interferometry. The massive star-forming region W75N contains three radio continuum sources (VLA1, VLA2, and VLA3), at three different evolutionary stages, and associated masers, while a large-scale molecular bipolar outflow is also present. Very recently, polarization observations of the 6.7 GHz methanol masers at milliarsecond resolution have been able to probe the strength and structure of the magnetic field over more than 2000 AU around VLA1. The magnetic field is parallel to the outflow, suggesting that VLA1 is its powering source. The observations of water masers at 22 GHz can give more information about the gas dynamics and the magnetic fields around VLA1 and VLA2. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22 GHz water masers in the star-forming region W75N. We detected 124 water masers, 36 around VLA1 and 88 around VLA2 of W75N, which indicate two different physical environments around the two sources, where VLA1 is in a more evolved state. The linear polarization of the masers confirms the tightly ordered magnetic field around VLA1, which is aligned with the large-scale molecular outflow, and also reveals an ordered magnetic field around VLA2, which is not parallel to the outflow. [abridged]

preprint2011arXiv

VLBA determination of the distance to nearby star-forming regions V. Dynamical mass, distance and radio structure of V773 Tau A

(ABRIDGED) We present multi-epoch Very Long Baseline Array (VLBA) observations of V773 Tau A, the 51-day binary subsystem in the multiple young stellar system V773 Tau. Combined with previous interferometric and radial velocity measurements, these new data enable us to improve the characterization of the physical orbit of the A subsystem. In particular, we infer updated dynamical masses for the primary and the secondary components of 1.55 pm 0.11 Msun, and 1.293 pm 0.068 Msun, respectively, and an updated orbital parallax distance to the system of 135.7 pm 3.2 pc, all consistent with previous estimates. Using the improved orbit, we can calculate the absolute coordinates of the barycenter of V773 Tau A at each epoch of our VLBA observations, and fit for its trigonometric parallax and proper motion. This provides a direct measurement of the distance to the system almost entirely independent of the orbit modeling. The best fit yields a distance of 129.9 pm 3.2 pc, in good agreement (i.e. within 1 sigma) with the distance estimate based on the orbital fit. Taking the mean value of the orbital and trigonometric parallaxes, we conclude that V773 Tau is located at d =132.8 pm 2.3 pc. The accuracy of this determination is nearly one order of magnitude better than that of previous estimates. In projection, V773 Tau and two other young stars (Hubble 4 and HDE 283572) recently observed with the VLBA are located toward the dark cloud Lynds 1495, in the central region of Taurus. These three stars appear to have similar trigonometric parallaxes, radial velocities, and proper motions, and we argue that the weighted mean and dispersion of their distances (d = 131.4 pc and sigma_d = 2.4 pc) provide a good estimate of the distance to and depth of Lynds 1495 and its associated stellar population.

preprint2011arXiv

VLBA SiO maser observations of the OH/IR star OH 44.8-2.3: magnetic field and morphology

SiO maser emission occurs in the extended atmosphere of evolved stars and can be studied at high angular resolution. Very long baseline interferometry (VLBI) observations of Mira variables indicate that SiO masers are significantly linearly polarized with linear polarization fraction up to 100%. However, no information is available at high angular resolution for SiO masers in higher mass loss OH/IR stars. We extend the VLBI SiO maser studies to OH/IR stars. The observations enable us to understand the SiO pumping mechanisms in higher mass loss evolved objects and compare those with Mira variables. Additionally, polarimetric observations of SiO masers help us to understand the magnetic field strength and morphology and to distinguish between conflicting polarization theories. The 43 GHz SiO maser observations of the OH/IR star OH 44.8-2.3 were performed with the VLBA in full polarization spectral line mode. Auxiliary EVLA observations were performed to allow for the absolute calibration of the polarization angle. The Zeeman splitting was measured by cross correlating the right and left circular polarization spectra as well as the S-curve fitting. Additionally, we analyzed the 1612 MHz OH maser observations of OH 44.8-2.3 from the VLA archive. The SiO masers of OH 44.8-2.2 form a ring located at ~5.4 AU around the star. The masers appear to be highly linearly polarized with fractional linear polarization up to 100%. The linear polarization vectors are consistent with a dipole field morphology in this star. We report a tentative detection of circular polarization of ~0.7% for the brightest maser feature. The magnetic field measured for this feature corresponds to 1.5\pm0.3 G. Additionally, the distribution of the 1612 MHz OH maser emission could indicate an elongated morphology.

preprint2010arXiv

Dynamics of the 6.7 and 12.2 GHz methanol masers around Cepheus A HW2

The 6.7 GHz methanol maser is exclusively associated with high-mass star formation. However, it remains unclear what structures harbour the methanol masers. Cepheus A is one of the closest regions of massive star formation, making it an excellent candidate for detailed studies. We determine the dynamics of maser spots in the high-mass star-forming region Cepheus A in order to infer where and when the maser emission occurs. Very long baseline interferometry (VLBI) observations of the 6.7 and 12.2 GHz methanol masers allows for mapping their spatial and velocity distribution. Phase-referencing is used to determine the astrometric positions of the maser emission, and multi-epoch observations can reveal 3D motions. The 6.7 GHz methanol masers are found in a filamentary structure over ~1350 AU, straddling the waist of the radio jet HW2. The positions agree well with previous observations of both the 6.7 and 12.2 GHz methanol masers. The velocity field of the maser spots does not show any sign of rotation, but is instead consistent with an infall signature. The 12.2 GHz methanol masers are closely associated with the 6.7 GHz methanol masers, and the parallax that we derive confirms previous measurements. We show that the methanol maser emission very likely arises in a shock interface in the equatorial region of Cepheus A HW2 and presents a model in which the maser emission occurs between the infalling gas and the accretion disk/process.

preprint2010arXiv

Magnetic field regulated infall on the disc around the massive protostar Cepheus A HW2

We present polarization observations of the 6.7-GHz methanol masers around the massive protostar Cepheus A HW2 and its associated disc. The data were taken with the Multi-Element Radio Linked Interferometer Network. The maser polarization is used to determine the full three-dimensional magnetic field structure around Cepheus A HW2. The observations suggest that the masers probe the large scale magnetic field and not isolated pockets of a compressed field. We find that the magnetic field is predominantly aligned along the protostellar outflow and perpendicular to the molecular and dust disc. From the three-dimensional magnetic field orientation and measurements of the magnetic field strength along the line of sight, we are able to determine that the high density material, in which the masers occurs, is threaded by a large scale magnetic field of ~23 mG. This indicates that the protostellar environment at ~1000 AU from Cepheus A HW2 is slightly supercritical (lambda approximately 1.7) and the relation between density and magnetic field is consistent with collapse along the magnetic field lines. Thus, the observations indicate that the magnetic field likely regulates accretion onto the disc. The magnetic field dominates the turbulent energies by approximately a factor of three and is sufficiently strong to be the crucial component stabilizing the massive accretion disc and sustaining the high accretion rates needed during massive star-formation.

preprint2010arXiv

Magnetic field structure in a high-mass outflow/disk system

To characterize the magnetic field structure of the outflow and core region within a prototypical high-mass star-forming region, we analyzed polarized CO(3-2) - for the first time observed with the Submillimeter Array -- as well as 880mum submm continuum emission from the high-mass outflow/disk system IRAS18089-1732. Both emission features with polarization degrees at a few percent level indicate that the magnetic field structure is largely aligned with the outflow/jet orientation from the small core scales to the larger outflow scales. Although quantitative estimates are crude, the analysis indicates that turbulent energy dominates over magnetic energy. The data also suggest a magnetic field strength increase from the lower-density envelope to the higher-density core.

preprint2009arXiv

Methanol masers probing the ordered magnetic field of W75N

The role of magnetic fields during the protostellar phase of high-mass star-formation is a debated topic. In particular, it is still unclear how magnetic fields influence the formation and dynamic of disks and outflows. Most current information on magnetic fields close to high-mass protostars comes from H2O and OH maser observations. Recently, the first 6.7 GHz methanol maser polarization observations were made, and they reveal strong and ordered magnetic fields. The morphology of the magnetic field during high-mass star-formation needs to be investigated on small scales, which can only be done using very long baseline interferometry observations. The massive star-forming regionW75N contains three radio sources and associated masers, while a large-scale molecular bipolar outflow is also present. Polarization observations of the 6.7 GHz methanol masers at high angular resolution probe the strength and structure of the magnetic field and determine its relation to the outflow. Eight of the European VLBI network antennas were used to measure the linear polarization and Zeeman-splitting of the 6.7 GHz methanol masers in the star-forming region W75N. We detected 10 methanol maser features, 4 of which were undetected in previous work. All arise near the source VLA1 of W75N. The linear polarization of the masers reveals a tightly ordered magnetic field over more than 2000 AU around VLA1 that is exactly aligned with the large-scale molecular outflow. This is consistent with the twisted magnetic field model proposed for explaining dust polarization observations. The Zeeman-splitting measured on 3 of the maser features indicates a dynamically important magnetic field in the maser region of the order of 50mG. We suggest VLA1 is the powering sources of the bipolar outflow.

preprint2009arXiv

Possible magnetic field variability during the 6.7 GHz methanol maser flares of G09.62+0.20

(Abridged) Recently, the magnetic field induced Zeeman splitting was measured for the strongest known 6.7 GHz methanol maser, which arises in the massive star forming region G09.62+0.20. This maser is one of a handful of periodically flaring methanol masers. The 100-m Effelsberg telescope was used to monitor the 6.7 GHz methanol masers of G09.62+0.20. With the exception of a two week period during the peak of the maser flare, we measure a constant magnetic field of B_||~11+-2 mG in the two strongest maser components of G09.62+0.20 that are separated by over 200 AU. In the two week period that coincides exactly with the peak of the maser flare of the strongest maser feature, we measure a sharp decrease and possible reversal of the Zeeman splitting. The exact cause of both maser and polarization variability is still unclear, but it could be related to either background amplification of polarized emission or the presence of a massive protostar with a close-by companion. Alternatively, the polarization variability could be caused by non-Zeeman effects related to the radiative transfer of polarized maser emission.

preprint2009arXiv

Precision Astrometry with the Very Long Baseline Array: Parallaxes and Proper Motions for 14 Pulsars

Astrometry can bring powerful constraints to bear on a variety of scientific questions about neutron stars, including their origins, astrophysics, evolution, and environments. Using phase-referenced observations at the VLBA, in conjunction with pulsar gating and in-beam calibration, we have measured the parallaxes and proper motions for 14 pulsars. The smallest measured parallax in our sample is 0.13+-0.02 mas for PSR B1541+09, which has a most probable distance of 7.2+1.3-1.1 kpc. We detail our methods, including initial VLA surveys to select candidates and find in-beam calibrators, VLBA phase-referencing, pulsar gating, calibration, and data reduction. The use of the bootstrap method to estimate astrometric uncertainties in the presence of unmodeled systematic errors is also described. Based on our new model-independent estimates for distance and transverse velocity, we investigate the kinematics and birth sites of the pulsars and revisit models of the Galactic electron density distribution. We find that young pulsars are moving away from the Galactic plane, as expected, and that age estimates from kinematics and pulsar spindown are generally in agreement, with certain notable exceptions. Given its present trajectory, the pulsar B2045-16 was plausibly born in the open cluster NGC 6604. For several high-latitude pulsars, the NE2001 electron density model underestimates the parallax distances by a factor of two, while in others the estimates agree with or are larger than the parallax distances, suggesting that the interstellar medium is irregular on relevant length scales. The VLBA astrometric results for the recycled pulsar J1713+0747 are consistent with two independent estimates from pulse timing, enabling a consistency check between the different reference frames.

preprint2009arXiv

The magnetic field of the evolved star W43A

The majority of the observed planetary nebulae exhibit elliptical or bipolar structures. Theoretical modeling has indicated that magnetically collimated jets may be responsible for the formation of the non-spherical planetary nebulae. The aim of this project is to measure the Zeeman splitting caused by the magnetic field in the OH and H2O maser regions occurring in the circumstellar envelope and bipolar outflow of the evolved star W43A. We report a measured magnetic field of approximately 100 micro-gauss in the OH maser region of the circumstellar envelope around W43A. The GBT observations reveal a magnetic field strength B|| of ~30 mG changing sign across the H2O masers at the tip of the red-shifted lobe of the bipolar outflow. We also find that the OH maser shell shows no sign of non-spherical expansion and that it probably has an expansion velocity that is typical for the shells of regular OH/IR stars. The GBT observations confirm that the magnetic field collimates the H2O maser jet, while the OH maser observations show that a strong large scale magnetic field is present in the envelope surrounding the W43A central star. The magnetic field in the OH maser envelope is consistent with the one extrapolated from the H2O measurements, confirming that magnetic fields play an important role in the entire circumstellar environment of W43A.

preprint2006arXiv

Intrinsic properties of the magnetically collimated water maser jet of W43A

Water maser polarization observations in the precessing jet of W43A have revealed that it is magnetically collimated. Here we present a detailed description of the physical properties of the water maser environment in the jet. We discuss the maser saturation level and beaming angle as well as the intrinsic temperatures and densities. Additionally, we show that the polarization angle of the strongest red-shifted maser feature undergoes a fast rotation of 90 degrees across the maser. Along with the variation of linear polarization fraction, this strongly supports the current theoretical description of maser linear polarization.

preprint2003arXiv

VLBI astrometry of circumstellar OH masers; proper motions and parallaxes of four AGB stars

The main-line OH masers around 4 AGB stars have been observed with the NRAO Very Long Baseline Array (VLBA) at 8 epochs over a period of 2.5 years. Using a phase referencing technique, the position of the most compact maser spot of each star was monitored with respect to two extragalactic reference sources. For U Her and W Hya, we observe the most blue-shifted maser spot, while for R Cas and S CrB we only detect a compact red-shifted maser spot. We managed to determine an accurate proper motion and parallax for U Her, R Cas and S CrB, while additional motion of the compact blue-shifted maser of W Hya is shown to possibly be related to the stellar pulsation. The motion and radio position are compared with the stellar trajectory and absolute optical position determined by the Hipparcos satellite. For U Her and W Hya, the most blue-shifted maser is consistent with the amplified stellar image. The new distances are compared with several published P--L relations, and in this respect the VLBI distances seem an improvement upon the Hipparcos distances.

preprint2002arXiv

Astrometry of circumstellar masers

The circumstellar masers around evolved stars offer an interesting possibility to measure stellar parameters through VLBI astrometry. In this paper the application of this technique is discussed, including the accuracy and the uncertainties of the method. The different maser species (OH, H_2O, SiO) have slightly different characteristics and applications. This paper does not concern astrometry of maser spots to study the kinematics of the envelope, but concentrates on attempting to measure the motion of the underlying star.

preprint2002arXiv

Magnetic Fields in the envelopes of late-type stars: Circular Polarization of H2O Masers

We present the first circular polarization measurements of circumstellar water masers around a sample of late-type stars. These observations are used to obtain the magnetic field strength in the water maser region with both an LTE and non-LTE analysis. We find fields from a few hundred milliGauss up to a few Gauss, indicating a solar-type r^(-2) dependence of the magnetic field on the distance. No linear polarization is detected to less than 1%.