Showing 1–50 of 51 results for author: Gray, M D

Maser Flares Driven by Isothermal Shock Waves

Authors: M. D. Gray, S. Etoka, B. Pimpanuwat, A. M. S. Richards

Abstract: We use 3D computer modelling to investigate the timescales and radiative output from maser flares generated by the impact of shock-waves on astronomical unit-scale clouds in interstellar and star-forming regions, and in circumstellar regions in some circumstances. Physical conditions are derived from simple models of isothermal hydrodynamic (single-fluid) and C-type (ionic and neutral fluid) shock… ▽ More We use 3D computer modelling to investigate the timescales and radiative output from maser flares generated by the impact of shock-waves on astronomical unit-scale clouds in interstellar and star-forming regions, and in circumstellar regions in some circumstances. Physical conditions are derived from simple models of isothermal hydrodynamic (single-fluid) and C-type (ionic and neutral fluid) shock-waves, and based on the ortho-H$_2$O 22-GHz transition. Maser saturation is comprehensively included, and we find that the most saturated maser inversions are found predominantly in the shocked material. We study the effect on the intensity, flux density and duration of flares of the following parameters: the pre-shock level of saturation, the observer's viewpoint, and the shock speed. Our models are able to reproduce observed flare rise times of a few times 10 days, specific intensities of up to 10$^5$ times the saturation intensity and flux densities of order $100(R/d)^2$Jy from a source of radius $R$ astronomical units at a distance of $d$ kiloparsec. We found that flares from C-type shocks are approximately 5 times more likely to be seen by a randomly placed observer than flares from hydrodynamically shocked clouds of similar dimensions. We computed intrinsic beaming patterns of the maser emission, finding substantial extension of the pattern parallel to the shock front in the hydrodynamic models. Beaming solid angles for hydrodynamic models can be as small as $1.3\times 10^{-5}$sr, but are an order of magnitude larger for C-type models. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS - 27 pages, 21 figures, 2 Tables

Resolved ALMA observations of water in the inner astronomical units of the HL Tau disk

Authors: Stefano Facchini, Leonardo Testi, Elizabeth Humphreys, Mathieu Vander Donckt, Andrea Isella, Ramon Wrzosek, Alain Baudry, Malcom D. Gray, Anita M. S. Richards, Wouter Vlemmmings

Abstract: The water molecule is a key ingredient in the formation of planetary systems, with the water snowline being a favourable location for the growth of massive planetary cores. Here we present Atacama Large Millimeter/ submillimeter Array data of the ringed protoplanetary disk orbiting the young star HL Tauri that show centrally peaked, bright emission arising from three distinct transitions of the ma… ▽ More The water molecule is a key ingredient in the formation of planetary systems, with the water snowline being a favourable location for the growth of massive planetary cores. Here we present Atacama Large Millimeter/ submillimeter Array data of the ringed protoplanetary disk orbiting the young star HL Tauri that show centrally peaked, bright emission arising from three distinct transitions of the main water isotopologue. The spatially and spectrally resolved water content probes gas in a thermal range down to the water sublimation temperature. Our analysis implies a stringent lower limit of 3.7 Earth oceans of water vapour available within the inner 17 astronomical units of the system. We show that our observations are limited to probing the water content in the atmosphere of the disk, due to the high dust column density and absorption, and indicate that the main water isotopologue is the best tracer to spatially resolve water vapour in protoplanetary disks. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: Published in Nature Astronomy (2024), arXiv version presents text and figures before proofing

ATOMIUM: Probing the inner wind of evolved O-rich stars with new, highly excited H$_2$O and OH lines

Authors: A. Baudry, K. T. Wong, S. Etoka, A. M. S. Richards, H. S. P. Muller, F. Herpin, T. Danilovich, M. D. Gray, S. Wallstrom, D. Gobrecht, T. Khouri, L. Decin, C. A. Gottlieb, K. M. Menten, W. Homan, T. J. Millar, M. Montarges, B. Pimpanuwat, J. M. C. Plane, P. Kervella

Abstract: Water and the hydroxyl radical are major constituents of the envelope of O-rich late-type stars. Transitions involving energy levels that are highly excited have been observed in both H$_2$O and OH. These and more recently discovered transitions can now be observed at a high sensitivity and angular resolution with the ALMA Array. Spectra and maps of H$_2$O and OH observed with an angular resolutio… ▽ More Water and the hydroxyl radical are major constituents of the envelope of O-rich late-type stars. Transitions involving energy levels that are highly excited have been observed in both H$_2$O and OH. These and more recently discovered transitions can now be observed at a high sensitivity and angular resolution with the ALMA Array. Spectra and maps of H$_2$O and OH observed with an angular resolution of 20 to $\sim$200 mas were obtained at two epochs with the ALMA array. Observations with the Compact Array were also used to check for time variability of water transitions. Radiative transfer models of water were revisited to characterize masing conditions and up-to-date chemical models were used for comparison with our observations. Ten rotational transitions of H$_2$O with energies up to 9000 K were observed in various vibrational states. All but one are new detections in space, and from these we have derived accurate rest frequencies. Hyperfine split $Λ$-doubling transitions in v = 0, J = 27/2 and 29/2 levels of the $^2Π_{3/2}$ state and, $J = 33/2$ and 35/2 of the $^2Π_{1/2}$ state of OH with excitation energies up to 8900 K were also observed. Four of these transitions are new detections in space. Combining our measurements with earlier observations of OH, the v = 0 and v = 1 $Λ$-doubling frequencies have been improved. Our H$_2$O maps show compact emission and extensions up to twelve stellar radii or more. The 268.149 GHz emission line of water in the v$_2$ = 2 state is time variable, tends to be masing with dominant radiative pumping, and is widely excited. The widespread but weaker 262.898 GHz water line in v$_2$ = 1 also shows signs of maser emission. Emission and absorption of both H$_2$O and OH reveal an infall of matter and complex kinematics influenced by binarity. From our observed column densities, we derived OH/H$_2$O abundance ratios in a few stars. △ Less

Submitted 4 May, 2023; originally announced May 2023.

Comments: 41 pages with references and 25 figures in main text. 4 Tables in Appendix A. 10 figures in Appendix B. 2 figures in Appendix C. 14 figures in Appendix D. 5 figures in Appendix E. 2 figures in Appendix F

Journal ref: A&A 674, A125 (2023)

The VLT/SPHERE view of the ATOMIUM cool evolved star sample. I. Overview: Sample characterization through polarization analysis

Authors: M. Montargès, E. Cannon, A. de Koter, T. Khouri, E. Lagadec, P. Kervella, L. Decin, I. McDonald, W. Homan, L. B. F. M. Waters, R. Sahai, C. A. Gottlieb, J. Malfait, S. Maes, B. Pimpanuwat, M. Jeste, T. Danilovich, F. De Ceuster, M. Van de Sande, D. Gobrecht, S. H. J. Wallström, K. T. Wong, I. El Mellah, J. Bolte, F. Herpin , et al. (10 additional authors not shown)

Abstract: Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars). Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution dat… ▽ More Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars). Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution data. To help interpret the polarized signal, we produced synthetic maps of light scattering by dust, through 3D radiative transfer simulations with the RADMC3D code. Results. The degree of linear polarization (DoLP) observed by ZIMPOL spreads across several optical filters. We infer that it primarily probes dust located just outside of the point spread function, and in or near the plane of the sky, with a total optical depth close to unity in the line of sight, representing only a fraction of the total circumstellar dust. The maximum DoLP ranges from 0.03-0.38 depending on the source, fractions that can be reproduced by our 3D pilot models for grains composed of common dust species. The spatial structure of the DoLP shows a diverse set of shapes. Only for three sources do we note a correlation between the ALMA CO and SiO lines, which trace the gas density, and the DoLP, which traces the dust. Conclusion. The clumpiness of the DoLP and the lack of a consistent correlation between the gas and the dust location show that, in the inner circumstellar environment (CSE), dust formation occurs at very specific sites. This has potential consequences for the derived mass-loss rates and dust-to-gas ratio in the inner region of the CSE. Except for $π^1$~Gru and perhaps GY Aql, we do not detect interactions between the circumstellar wind and the hypothesized companions that shape the wind at larger scales. This suggests that the orbits of any other companions are tilted out of the plane of the sky. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 22 pages, 15 figures, 5 tables

Journal ref: A&A 671, A96 (2023)

Polarized Maser Emission with In-Source Faraday Rotation

Authors: Taylor L. Tobin, Malcolm D. Gray, Athol J. Kemball

Abstract: We discuss studies of polarization in astrophysical masers with particular emphasis on the case where the Zeeman splitting is small compared to the Doppler profile, resulting in a blend of the transitions between magnetic substates. A semi-classical theory of the molecular response is derived, and coupled to radiative transfer solutions for 1 and 2-beam linear masers, resulting in a set of non-lin… ▽ More We discuss studies of polarization in astrophysical masers with particular emphasis on the case where the Zeeman splitting is small compared to the Doppler profile, resulting in a blend of the transitions between magnetic substates. A semi-classical theory of the molecular response is derived, and coupled to radiative transfer solutions for 1 and 2-beam linear masers, resulting in a set of non-linear, algebraic equations for elements of the molecular density matrix. The new code, PRISM, implements numerical methods to compute these solutions. Using PRISM, we demonstrate a smooth transfer between this case and that of wider splitting. For a J=1-0 system, with parameters based on the $v=1, J=1-0$ transition of SiO, we investigate the behaviour of linear and circular polarization as a function of the angle between the propagation axis and the magnetic field, and with the optical depth, or saturation state, of the model. We demonstrate how solutions are modified by the presence of Faraday rotation, generated by various abundances of free electrons, and that strong Faraday rotation leads to additional angles where Stokes-Q changes sign. We compare our results to a number of previous models, from the analytical limits derived by Goldreich, Keeley and Kwan in 1973, through computational results by W. Watson and co-authors, to the recent work by Lankhaar and Vlemmings in 2019. We find that our results are generally consistent with those of other authors given the differences of approach and the approximations made. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: 36 pages, 15 figures. Accepted for publication in ApJ

Sciences with Thai National Radio Telescope

Authors: Phrudth Jaroenjittichai, Koichiro Sugiyama, Busaba H. Kramer, Boonrucksar Soonthornthum, Takuya Akahori, Kitiyanee Asanok, Willem Baan, Sherin Hassan Bran, Shari L. Breen, Se-Hyung Cho, Thanapol Chanapote, Richard Dodson, Simon P. Ellingsen, Sandra Etoka, Malcolm D. Gray, James A. Green, Kazuhiro Hada, Marcus Halson, Tomoya Hirota, Mareki Honma, Hiroshi Imai, Simon Johnston, Kee-Tae Kim, Michael Kramer, Di Li , et al. (22 additional authors not shown)

Abstract: This White Paper summarises potential key science topics to be achieved with Thai National Radio Telescope (TNRT). The commissioning phase has started in mid 2022. The key science topics consist of "Pulsars and Fast Radio Bursts (FRBs)", "Star Forming Regions (SFRs)", "Galaxy and Active Galactic Nuclei (AGNs)", "Evolved Stars", "Radio Emission of Chemically Peculiar (CP) Stars", and "Geodesy", cov… ▽ More This White Paper summarises potential key science topics to be achieved with Thai National Radio Telescope (TNRT). The commissioning phase has started in mid 2022. The key science topics consist of "Pulsars and Fast Radio Bursts (FRBs)", "Star Forming Regions (SFRs)", "Galaxy and Active Galactic Nuclei (AGNs)", "Evolved Stars", "Radio Emission of Chemically Peculiar (CP) Stars", and "Geodesy", covering a wide range of observing frequencies in L/C/X/Ku/K/Q/W-bands (1-115 GHz). As a single-dish instrument, TNRT is a perfect tool to explore time domain astronomy with its agile observing systems and flexible operation. Due to its ideal geographical location, TNRT will significantly enhance Very Long Baseline Interferometry (VLBI) arrays, such as East Asian VLBI Network (EAVN), Australia Long Baseline Array (LBA), European VLBI Network (EVN), in particular via providing a unique coverage of the sky resulting in a better complete "uv" coverage, improving synthesized-beam and imaging quality with reducing side-lobes. This document highlights key science topics achievable with TNRT in single-dish mode and in collaboration with VLBI arrays. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: 75 pages, 16 figures, 15 tables. White Paper for Potential Key Sciences to be Achieved with Thai National Radio Telescope (TNRT)

Stringent limits on $^{28}$SiO maser emission from the recurrent nova T Coronae Borealis

Authors: A. Evans, B. Pimpanuwat, A. M. S. Richards, D. P. K. Banerjee, U. Munari, M. D. Gray, B. Hutawarakorn Kramer, A. Kraus

Abstract: There are indications that the third known eruption of the recurrent nova T CrB is imminent, and multi-wavelength observations prior to the eruption are important to characterise the system before it erupts. T CrB is known to display the SiO fundamental vibrational feature at 8$\,μ$m. When the anticipated eruption occurs, it is possible that the shock produced when the ejected material runs into t… ▽ More There are indications that the third known eruption of the recurrent nova T CrB is imminent, and multi-wavelength observations prior to the eruption are important to characterise the system before it erupts. T CrB is known to display the SiO fundamental vibrational feature at 8$\,μ$m. When the anticipated eruption occurs, it is possible that the shock produced when the ejected material runs into the wind of the red giant in the system may be traced using SiO maser emission. We have used the 100m Effelsberg Radio Telescope to search for $^{28}$SiO emission in the $\upsilon=1$, $\upsilon=2$, $J=1\rightarrow0$ transitions, at 43.122 GHz and 42.820~GHz respectively, while the system is in quiescence. We find no evidence for such emission. △ Less

Submitted 29 July, 2022; originally announced July 2022.

Comments: Published in MNRAS, 514, 4894 - 4897 (2022)

Journal ref: MNRAS, 514, 4894 - 4897 (2022)

Atomium: The astounding complexity of the near circumstellar environment of the M-type AGB star R Hydrae. I. Morpho-kinematical interpretation of CO and SiO emission

Authors: Ward Homan, Bannawit Pimpanuwat, Fabrice Herpin, Taissa Danilovich, Iain McDonald, Sofia H. J. Wallström, Anita M. S. Richards, Alain Baudry, Raghvendra Sahai, Tom J. Millar, Alex de Koter, C. A. Gottlieb, Pierre Kervella, Miguel Montargès, Marie Van de Sande, Leen Decin, Albert Zijlstra, Sandra Etoka, Manali Jeste, Holger S. P. Müller, Silke Maes, Jolien Malfait, Karl Menten, John Plane, Kelvin Lee , et al. (14 additional authors not shown)

Abstract: Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we… ▽ More Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we observed the M-type asymptotic giant branch (AGB) star R Hydrae with ALMA. The morphology of the inner wind of R Hya, which has a known companion at ~3500 au, was determined from maps of CO and SiO obtained at high angular resolution. A map of the CO emission reveals a multi-layered structure consisting of a large elliptical feature at an angular scale of ~10'' that is oriented along the north-south axis. The wind morphology within the elliptical feature is dominated by two hollow bubbles. The bubbles are on opposite sides of the AGB star and lie along an axis with a position angle of ~115 deg. Both bubbles are offset from the central star, and their appearance in the SiO channel maps indicates that they might be shock waves travelling through the AGB wind. An estimate of the dynamical age of the bubbles yields an age of the order of 100 yr, which is in agreement with the previously proposed elapsed time since the star last underwent a thermal pulse. When the CO and SiO emission is examined on subarcsecond angular scales, there is evidence for an inclined, differentially rotating equatorial density enhancement, strongly suggesting the presence of a second nearby companion. The position angle of the major axis of this disc is ~70 deg in the plane of the sky. We tentatively estimate that a lower limit on the mass of the nearby companion is ~0.65 Msol on the basis of the highest measured speeds in the disc and the location of its inner rim at ~6 au from the AGB star. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: 21 pages, 23 figures

Journal ref: A&A 651, A82 (2021)

Nitric Oxide and other molecules: Molecular Modelling and Low FrequencyExploration using the Murchison Widefield Array

Authors: Chenoa D. Tremblay, Malcolm D. Gray, Natasha Hurley-Walker, James A. Green, Joanne R. Dawson R. Dawson, John M. Dickey, Paul A. Jones, Steven J. Tingay, O. Ivy Wong

Abstract: We present new molecular modelling for 14NO and 15NO and a deep, blind molecular line survey at low radio frequencies (99-129 MHz). This survey is the third in a series completed with the Murchison Widefield Array (MWA), but in comparison with the previous surveys, uses four times more data (17 hours vs. 4 hours) and is three times better in angular resolution (1' vs. 3'). The new molecular modell… ▽ More We present new molecular modelling for 14NO and 15NO and a deep, blind molecular line survey at low radio frequencies (99-129 MHz). This survey is the third in a series completed with the Murchison Widefield Array (MWA), but in comparison with the previous surveys, uses four times more data (17 hours vs. 4 hours) and is three times better in angular resolution (1' vs. 3'). The new molecular modelling for nitric oxide and its main isotopologue has seven transitions within the MWA frequency band (although we also present the higher frequency transitions). Although we did not detect any new molecular lines at a limit of 0.21 Jy beam^-1, this work is an important step in understanding the data processing challenges for the future Square Kilometre Array (SKA) and places solid limits on what is expected in the future of low-frequency surveys. The modelling can be utilised for future searches of nitric oxide. △ Less

Submitted 19 October, 2020; originally announced October 2020.

Comments: 13 pages, 5 figures, 1 table, accepted by Astrophysical Journals (ApJ)

Analysis of Methanol Maser Flares in G107.298+5.63 andS255-NIRS3

Authors: M. D. Gray, S. Etoka, B. Pimpanuwat

Abstract: A 3D maser model has been used to perform an inverse problem on the light curves from three high-amplitude maser flares, selected on the basis of contemporaneous infra-red observations. Plots derived from the model recover the size of the maser cloud, and two parameters linked to saturation, from three observational properties of the light curve. Recovered sizes are consistent with independent int… ▽ More A 3D maser model has been used to perform an inverse problem on the light curves from three high-amplitude maser flares, selected on the basis of contemporaneous infra-red observations. Plots derived from the model recover the size of the maser cloud, and two parameters linked to saturation, from three observational properties of the light curve. Recovered sizes are consistent with independent interferometric measurements. Maser objects transition between weak and moderate saturation during a flare. △ Less

Submitted 23 June, 2020; originally announced June 2020.

Comments: 5 pages, 5 figures, accepted for publication in MNRAS

Maser Flares Driven by Variations in Pumping and Background Radiation

Authors: M. D. Gray, S. Etoka, A. Travis, B. Pimpanuwat

Abstract: We simulate maser flares by varying either the pump rate or the background level of radiation in a 3D model of a maser cloud. We investigate the effect of different cloud shapes, saturation levels and viewpoints. Results are considered for clouds with both uniform and internally variable unsaturated inversion. Pumping and background variations are represented by several different driving functions… ▽ More We simulate maser flares by varying either the pump rate or the background level of radiation in a 3D model of a maser cloud. We investigate the effect of different cloud shapes, saturation levels and viewpoints. Results are considered for clouds with both uniform and internally variable unsaturated inversion. Pumping and background variations are represented by several different driving functions, some of which are light curves drawn from observations. We summarise the pumping variability results in terms of three observable parameters, the maximum flux density achieved, a variability index and duty cycle. We demonstrate typical ranges of the flux density that may result from viewing an aspherical object from random viewpoints. The best object for a flare is a prolate cloud, viewed close to its long axis and driven from unsaturated conditions to at least modest saturation. Results for variation of the background level are qualitatively different from the variable pumping results in that they tend to produce short intervals of low flux density under conditions of moderate saturation and sufficient variability to be consistent with strong flaring. Variable background models typically have a significantly higher duty cycle than those with variable pumping. △ Less

Submitted 11 February, 2020; originally announced February 2020.

Comments: 19 pages, 15 figures, accepted for publication in MNRAS

Maser Flare Simulations from Oblate and Prolate Clouds

Authors: M. D. Gray, J. Baggott, J. Westlake, S. Etoka

Abstract: We investigated, through numerical models, the flaring variability that may arise from the rotation of maser clouds of approximately spheroidal geometry, ranging from strongly oblate to strongly prolate examples. Inversion solutions were obtained for each of these examples over a range of saturation levels from unsaturated to highly saturated. Formal solutions were computed for rotating clouds wit… ▽ More We investigated, through numerical models, the flaring variability that may arise from the rotation of maser clouds of approximately spheroidal geometry, ranging from strongly oblate to strongly prolate examples. Inversion solutions were obtained for each of these examples over a range of saturation levels from unsaturated to highly saturated. Formal solutions were computed for rotating clouds with many randomly chosen rotation axes, and corresponding averaged maser light curves plotted with statistical information. The dependence of results on the level of saturation and on the degree of deformation from the spherical case were investigated in terms of a variability index and duty cycle. It may be possible to distinguish observationally between flares from oblate and prolate objects. Maser flares from rotation are limited to long timescales (at least a few years) and modest values of the variability index ($\lesssim 100$), and can be aperiodic or quasi-periodic. Rotation is therefore not a good model for H$_2$O variability on timescales of weeks to months, or of truly periodic flares. △ Less

Submitted 23 April, 2019; originally announced April 2019.

Comments: 11 pages, 12 figures, accepted for publication in MNRAS

Constraining Theories of Polarized SiO Maser Transport: Multi-Epoch Analysis of a $π/2$ Electric Vector Rotation Feature

Authors: Taylor L. Tobin, Athol J. Kemball, Malcolm D. Gray

Abstract: The detailed polarization mechanisms of SiO masers originating from the near circum\-stellar environment of Asymptotic Giant Branch stars are not yet definitively known. Prevailing theories are broadly classified as either Zeeman or non-Zeeman in origin, the latter including effects such as anisotropic pumping or anisotropic resonant scattering. The predicted behavior of the linear and circular po… ▽ More The detailed polarization mechanisms of SiO masers originating from the near circum\-stellar environment of Asymptotic Giant Branch stars are not yet definitively known. Prevailing theories are broadly classified as either Zeeman or non-Zeeman in origin, the latter including effects such as anisotropic pumping or anisotropic resonant scattering. The predicted behavior of the linear and circular polarization fractions and electric vector position angle vary by theory. In particular, individual maser features that exhibit a rotation in linear polarization of $\sim π/2$ as a function of frequency over their extent can be utilized as a test of several maser polarization transport theories. In this paper, we analyze one SiO ($ν=1$, $J=1-0$) maser feature toward the Mira variable, TX Cam that exhibits this internal polarization rotation and persists across five epochs (spanning $\sim3$ months). We compare our results to the predictions by several maser polarization theories and find that the linear polarization across the feature is consistent with a geometric effect for a saturated maser originating when the angle between the projected magnetic field and the line of sight ($θ$) crosses the Van Vleck angle $θ_F \sim 55^{\circ}$. However, the electric vector position angle (EVPA) exhibits a smooth rotation across the spatial extent of the feature rather than the expected abrupt $π/2$ flip. We discuss possible explanations for this discrepancy and alternative theories. Circular polarization across the feature is also analyzed and it is the most accurately described by Zeeman effects giving rise to a circular polarization fraction of the form $m_c$ is approximately proportional to $\cos θ$. △ Less

Submitted 22 December, 2018; originally announced December 2018.

Comments: Figure Sets included as Appendix B

A New 3D Maser Code Applied to Flaring Events

Authors: M. D. Gray, L. Mason, S. Etoka

Abstract: We set out the theory and discretization scheme for a new finite-element computer code, written specifically for the simulation of maser sources. The code was used to compute fractional inversions at each node of a 3-D domain for a range of optical thicknesses. Saturation behaviour of the nodes with regard to location and optical depth were broadly as expected. We have demonstrated via formal solu… ▽ More We set out the theory and discretization scheme for a new finite-element computer code, written specifically for the simulation of maser sources. The code was used to compute fractional inversions at each node of a 3-D domain for a range of optical thicknesses. Saturation behaviour of the nodes with regard to location and optical depth were broadly as expected. We have demonstrated via formal solutions of the radiative transfer equation that the apparent size of the model maser cloud decreases as expected with optical depth as viewed by a distant observer. Simulations of rotation of the cloud allowed the construction of light-curves for a number of observable quantities. Rotation of the model cloud may be a reasonable model for quasi-periodic variability, but cannot explain periodic flaring. △ Less

Submitted 5 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS, 2018 February 23, 12 pages, 6 figures

Sun-Sized Water Vapor Masers in Cepheus A

Authors: A. M. Sobolev, J. M. Moran, M. D. Gray, A. Alakoz, H. Imai, W. A. Baan, A. M. Tolmachev, V. A. Samodurov, D. A. Ladeyshchikov

Abstract: We present the first VLBI observations of a Galactic water maser (in Chepeus A) made with a very long baseline interferometric array involving the RadioAstron Earth-orbiting satellite station as one of its elements. We detected two distinct components at -16.9 and 0.6 km/s with a fringe spacing of 66 microarcseconds. In total power, the 0.6 km/s component appears to be a single Gaussian component… ▽ More We present the first VLBI observations of a Galactic water maser (in Chepeus A) made with a very long baseline interferometric array involving the RadioAstron Earth-orbiting satellite station as one of its elements. We detected two distinct components at -16.9 and 0.6 km/s with a fringe spacing of 66 microarcseconds. In total power, the 0.6 km/s component appears to be a single Gaussian component of strength 580 Jy and width of 0.7 km/s. Single-telescope monitoring showed that its lifetime was only 8~months. The absence of a Zeeman pattern implies the longitudinal magnetic field component is weaker than 120 mG. The space-Earth cross power spectrum shows two unresolved components smaller than 15 microarcseconds, corresponding to a linear scale of 1.6 x 10^11 cm, about the diameter of the Sun, for a distance of 700 pc, separated by 0.54 km/s in velocity and by 160 +/-35 microarcseconds in angle. This is the smallest angular structure ever observed in a Galactic maser. The brightness temperatures are greater than 2 x 10^14K, and the line widths are 0.5 km/s. Most of the flux (about 87%) is contained in a halo of angular size of 400 +/- 150 microarcseconds. This structure is associated with the compact HII region HW3diii. We have probably picked up the most prominent peaks in the angular size range of our interferometer. We discuss three dynamical models: (1) Keplerian motion around a central object, (2) two chance overlapping clouds, and (3) vortices caused by flow around an obstacle (i.e., von Karman vortex street) with Strouhal number of about~0.3. △ Less

Submitted 19 February, 2018; originally announced February 2018.

Comments: 15 pages, 9 figures. Accepted for publication in ApJ, February 16, 2018

Vibrationally excited water emission at 658 GHz from evolved stars

Authors: A. Baudry, E. M. L. Humphreys, F. Herpin, K. Torstensson, W. H. T. Vlemmings, A. M. S. Richards, M. D. Gray, C. De Breuck, M. Olberg

Abstract: Several rotational transitions of water have been identified toward evolved stars in the ground vibrational state as well as in the first excited state of the bending mode. In the latter vibrational state of water, the 658 GHz J = 1_1,0-1_0,1 rotational transition is often strong and seems to be widespread in late-type stars. Our main goals are to better characterize the nature of the 658 GHz emis… ▽ More Several rotational transitions of water have been identified toward evolved stars in the ground vibrational state as well as in the first excited state of the bending mode. In the latter vibrational state of water, the 658 GHz J = 1_1,0-1_0,1 rotational transition is often strong and seems to be widespread in late-type stars. Our main goals are to better characterize the nature of the 658 GHz emission, compare the velocity extent of the 658 GHz emission with SiO maser emission to help locate the water layers and, more generally, investigate the physical conditions prevailing in the excited water layers of evolved stars. Another goal is to identify new 658 GHz emission sources and contribute in showing that this emission is widespread in evolved stars. Eleven evolved stars were extracted from our mini-catalog of existing and potential 658 GHz sources for observations with the APEX telescope equipped with the SEPIA receiver. The 13CO J=6-5 line was placed in the same receiver sideband for simultaneous observation with the 658 GHz line of water. We have compared the 658 GHz line properties with our H2O radiative transfer models in stars and we have compared the velocity ranges of the 658 GHz and SiO J=2-1, v=1 maser lines. All stars show 658 GHz emission with a peak flux density in the range 50-70 Jy to 2000-3000 Jy. We have shown that the 658 GHz line is masing and we found that the 658 GHz velocity extent tends to be correlated with that of the SiO maser suggesting that both emission lines are excited in circumstellar layers close to the central star. Broad and stable line profiles are observed at 658 GHz. This could indicate maser saturation although we have tentatively provided first information on time variability at 658 GHz. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: 9 pages, 5 figures. Accepted in A&A

Journal ref: A&A 609, A25 (2018)

VLBA imaging of the 3mm SiO maser emission in the disk-wind from the massive protostellar system Orion Source I

Authors: S. Issaoun, C. Goddi, L. D. Matthews, L. J. Greenhill, M. D. Gray, E. M. L. Humphreys, C. J. Chandler, M. Krumholz, H. Falcke

Abstract: We present the first images of the 28SiO v=1, J=2-1 maser emission around the closest known massive young stellar object Orion Source I observed at 86 GHz (3mm) with the VLBA. These images have high spatial (~0.3 mas) and spectral (~0.054 km/s) resolutions. We find that the 3mm masers lie in an X-shaped locus consisting of four arms, with blue-shifted emission in the south and east arms and red-sh… ▽ More We present the first images of the 28SiO v=1, J=2-1 maser emission around the closest known massive young stellar object Orion Source I observed at 86 GHz (3mm) with the VLBA. These images have high spatial (~0.3 mas) and spectral (~0.054 km/s) resolutions. We find that the 3mm masers lie in an X-shaped locus consisting of four arms, with blue-shifted emission in the south and east arms and red-shifted emission in the north and west arms. Comparisons with previous images of the 28SiO v=1,2, J=1-0 transitions at 7mm (observed in 2001-2002) show that the bulk of the J=2-1 transition emission follows the streamlines of the J=1-0 emission and exhibits an overall velocity gradient consistent with the gradient at 7mm. While there is spatial overlap between the 3mm and 7mm transitions, the 3mm emission, on average, lies at larger projected distances from Source I (~44 AU compared with ~35 AU for 7mm). The spatial overlap between the v=1, J=1-0 and J=2-1 transitions is suggestive of a range of temperatures and densities where physical conditions are favorable for both transitions of a same vibrational state. However, the observed spatial offset between the bulk of emission at 3mm and 7mm possibly indicates different ranges of temperatures and densities for optimal excitation of the masers. We discuss different maser pumping models that may explain the observed offset. We interpret the 3mm and 7mm masers as being part of a single wide-angle outflow arising from the surface of an edge-on disk rotating about a northeast-southwest axis, with a continuous velocity gradient indicative of differential rotation consistent with a Keplerian profile in a high-mass proto-binary. △ Less

Submitted 24 July, 2017; originally announced July 2017.

Comments: 11 pages, 12 figures; accepted for publication in A&A

Journal ref: A&A 606, A126 (2017)

First detection of THz water maser in NGC7538-IRS1 with SOFIA and new 22 GHz e-MERLIN maps

Authors: F. Herpin, A. Baudry, A. M. S. Richards, M. D. Gray, N. Schneider, K. M. Menten, F. Wyrowski, S. Bontemps, R. Simon, H. Wiesemeyer

Abstract: The formation of massive stars is still not well understood. Accumulating a large amount of mass infalling within a single entity in spite of radiation pressure is possible if, among several other conditions, enough thermal energy is released. Despite numerous water line observations, with the Herschel Space Observatory, in most of the sources observations were not able to trace the emission from… ▽ More The formation of massive stars is still not well understood. Accumulating a large amount of mass infalling within a single entity in spite of radiation pressure is possible if, among several other conditions, enough thermal energy is released. Despite numerous water line observations, with the Herschel Space Observatory, in most of the sources observations were not able to trace the emission from the hot core around the newly forming protostellar object. We want to probe the physical conditions and water abundance in the inner layers of the host protostellar object NGC7538-IRS1 using a highly excited H2O line. Water maser models predict that several THz water masers should be detectable in these objects. We present SOFIA observations of the o-H2O 8(2,7)-7(3,4) line at 1296.41106 GHz and a 6(1,6)-5(2,3) 22 GHz e-MERLIN map of the region (first-ever 22 GHz images made after the e-MERLIN upgrade). In order to be able to constrain the nature of the emission - thermal or maser - we use near-simultaneous observations of the 22 GHz water maser performed with the Effelsberg radiotelescope and e-MERLIN. A thermal water model using the RATRAN radiative transfer code is presented based on HIFI pointed observations. Molecular water abundances are derived for the hot core. The H2O 8(2,7)- 7(3,4) line is detected toward NGC7538-IRS1 with one feature at the source velocity (-57.7 km/s) and another one at -48.4 km/s. We propose that the emission at the source velocity is consistent with thermal excitation and is excited in the innermost part of the IRS1a massive protostellar object's closest circumstellar environment. The other emission is very likely the first detection of a water THz maser line, pumped by shocks due to IRS1b outflow, in a star-forming region. Assuming thermal excitation of the THz line, the water abundance in NGC7538-IRS1's hot core is estimated to be 5.2x10^{-5} with respect to H2. △ Less

Submitted 20 June, 2017; originally announced June 2017.

Comments: accepted to A&A, 15 pages

Journal ref: A&A 606, A52 (2017)

Simultaneous 183 GHz H2O Maser and SiO Observations Towards Evolved Stars Using APEX SEPIA Band 5

Authors: E. M. L. Humphreys, K. Immer, M. D. Gray, E. De Beck, W. H. T. Vlemmings, A. Baudry, A. M. S. Richards, M. Wittkowski, K. Torstensson, C. De Breuck, P. Moller, S. Etoka, M. Olberg

Abstract: We investigate the use of 183 GHz H2O masers for characterization of the physical conditions and mass loss process in the circumstellar envelopes of evolved stars. We used APEX SEPIA Band 5 to observe the 183 GHz H2O line towards 2 Red Supergiant and 3 Asymptotic Giant Branch stars. Simultaneously, we observed lines in 28SiO v0, 1, 2 and 3, and for 29SiO v0 and 1. We detected the 183 GHz H2O line… ▽ More We investigate the use of 183 GHz H2O masers for characterization of the physical conditions and mass loss process in the circumstellar envelopes of evolved stars. We used APEX SEPIA Band 5 to observe the 183 GHz H2O line towards 2 Red Supergiant and 3 Asymptotic Giant Branch stars. Simultaneously, we observed lines in 28SiO v0, 1, 2 and 3, and for 29SiO v0 and 1. We detected the 183 GHz H2O line towards all the stars with peak flux densities greater than 100 Jy, including a new detection from VY CMa. Towards all 5 targets, the water line had indications of being due to maser emission and had higher peak flux densities than for the SiO lines. The SiO lines appear to originate from both thermal and maser processes. Comparison with simulations and models indicate that 183 GHz maser emission is likely to extend to greater radii in the circumstellar envelopes than SiO maser emission and to similar or greater radii than water masers at 22, 321 and 325 GHz. We speculate that a prominent blue-shifted feature in the W Hya 183 GHz spectrum is amplifying the stellar continuum, and is located at a similar distance from the star as mainline OH maser emission. From a comparison of the individual polarizations, we find that the SiO maser linear polarization fraction of several features exceeds the maximum fraction allowed under standard maser assumptions and requires strong anisotropic pumping of the maser transition and strongly saturated maser emission. The low polarization fraction of the H2O maser however, fits with the expectation for a non-saturated maser. 183 GHz H2O masers can provide strong probes of the mass loss process of evolved stars. Higher angular resolution observations of this line using ALMA Band 5 will enable detailed investigation of the emission location in circumstellar envelopes and can also provide information on magnetic field strength and structure. △ Less

Submitted 7 April, 2017; originally announced April 2017.

Comments: 11 pages, accepted by A&A. The abstract here has been shortened

Towards detecting methanol emission in low-mass protoplanetary discs with ALMA: The role of non-LTE excitation

Authors: S. Yu. Parfenov, D. A. Semenov, A. M. Sobolev, M. D. Gray

Abstract: The understanding of organic content of protoplanetary discs is one of the main goals of the planet formation studies. As an attempt to guide the observational searches for weak lines of complex species in discs, we modelled the (sub-)millimetre spectrum of gaseous methanol (CH$_3$OH), one of the simplest organic molecules, in the representative T Tauri system. We used 1+1D disc physical model cou… ▽ More The understanding of organic content of protoplanetary discs is one of the main goals of the planet formation studies. As an attempt to guide the observational searches for weak lines of complex species in discs, we modelled the (sub-)millimetre spectrum of gaseous methanol (CH$_3$OH), one of the simplest organic molecules, in the representative T Tauri system. We used 1+1D disc physical model coupled to the gas-grain ALCHEMIC chemical model with and without 2D-turbulent mixing. The computed CH$_3$OH abundances along with the CH$_3$OH scheme of energy levels of ground and excited torsional states were used to produce model spectra obtained with the non-local thermodynamic equilibrium (non-LTE) 3D line radiative transfer code LIME. We found that the modelled non-LTE intensities of the CH$_3$OH lines can be lower by factor of $>10$--$100$ than those calculated under assumption of LTE. Though population inversion occurs in the model calculations for many (sub-)millimetre transitions, it does not lead to the strong maser amplification and noticeably high line intensities. We identify the strongest CH$_3$OH (sub-)millimetre lines that could be searched for with the Atacama Large Millimeter Array (ALMA) in nearby discs. The two best candidates are the CH$_{3}$OH $5_0-4_0~A^+$ (241.791 GHz) and $5_{-1}-4_{-1}~E$ (241.767 GHz) lines, which could possibly be detected with the $\sim5σ$ signal-to-noise ratio after $\sim3$ hours of integration with the full ALMA array. △ Less

Submitted 12 May, 2016; originally announced May 2016.

Comments: Accepted for publication in MNRAS, 15 figures, 3 tables

Excited-state hydroxyl maser catalogue from the methanol multibeam survey -- I. Positions and Variability

Authors: A. Avison, L. J. Quinn, G. A. Fuller, J. L. Caswell, J. A. Green, S. L. Breen, S. P. Ellingsen, M. D. Gray, M. Pestalozzi, M. A. Thompson, M. A. Voronkov

Abstract: We present the results of the first complete unbaised survey of the Galactic Plane for 6035-MHz excited-state hydroxyl masers undertaken as part of the Methanol Multibeam Survey. These observations cover the Galactic longitude ranges $186^{\circ}< l < 60^{\circ}$ including the Galactic Centre. We report the detection of 127 excited-state hydroxyl masers within the survey region, 47 being new sourc… ▽ More We present the results of the first complete unbaised survey of the Galactic Plane for 6035-MHz excited-state hydroxyl masers undertaken as part of the Methanol Multibeam Survey. These observations cover the Galactic longitude ranges $186^{\circ}< l < 60^{\circ}$ including the Galactic Centre. We report the detection of 127 excited-state hydroxyl masers within the survey region, 47 being new sources. The positions of new detections were determined from interferometric observations with the Australia Telescope Compact Array. We discuss the association of 6035-MHz masers in our survey with the 6668-MHz masers from the MMB Survey, finding 37 likely methanol-excited-state hydroxyl masers maser pairs with physical separations of <=0.03pc and 55 pairings separated by <=0.1pc. Using these we calculate for the first time an excited-state hydroxyl maser life time of between 3.3x10^3 and 8.3x10^3 years. We also discuss the variability of the 6035-MHz masers and detection rates of counterpart 6030-MHz excited-state hydroxyl masers (28% of our sample having detection at both frequencies). △ Less

Submitted 29 April, 2020; v1 submitted 9 May, 2016; originally announced May 2016.

Comments: 32 pages, 11 figures, 4 tables, accepted for publication in MNRAS; Typo in figure 5 caption corrected in 2nd version; Transposition error of words 'right' and 'left' in Figure 11 caption corrected in 3rd version

First Detection of $^3$He$^+$ in the Planetary Nebula IC$\,$418

Authors: L. Guzman-Ramirez, J. R. Rizzo, A. A. Zijlstra, C. García-Miró, C. Morisset, M. D. Gray

Abstract: The $^3$He isotope is important to many fields of astrophysics, including stellar evolution, chemical evolution, and cosmology. The isotope is produced in low-mass stars which evolve through the planetary nebula (PN) phase. $^3$He abundances in PNe can help test models of the chemical evolution of the Galaxy. We present the detection of the $^3$He$^+$ emission line using the single dish Deep Space… ▽ More The $^3$He isotope is important to many fields of astrophysics, including stellar evolution, chemical evolution, and cosmology. The isotope is produced in low-mass stars which evolve through the planetary nebula (PN) phase. $^3$He abundances in PNe can help test models of the chemical evolution of the Galaxy. We present the detection of the $^3$He$^+$ emission line using the single dish Deep Space Station 63, towards the PN IC$\,$418. We derived a $^3$He/H abundance in the range 1.74$\pm$0.8$\times$10$^{-3}$ to 5.8$\pm$1.7$\times$10$^{-3}$, depending on whether part of the line arises in an outer ionized halo. The lower value for $^3$He/H ratio approaches values predicted by stellar models which include thermohaline mixing, but requires that large amounts of $^3$He are produced inside low-mass stars which enrich the interstellar medium (ISM). However, this over-predicts the $^3$He abundance in HII regions, the ISM, and proto-solar grains, which is known to be of the order of 10$^{-5}$. This discrepancy questions our understanding of the evolution of the $^3$He, from circumstellar environments to the ISM. △ Less

Submitted 10 April, 2016; originally announced April 2016.

Comments: 5 pages, 4 figures, accepted for publication in the Monthly Notices of the Royal Astronomical Society Letters

The Physics of Water Masers observable with ALMA and SOFIA: Model Predictions for Evolved Stars

Authors: M. D. Gray, A. Baudry, A. M. S. Richards, E. M. L. Humphreys, A. M. Sobolev, J. A. Yates

Abstract: We present the results of models that were designed to study all possible water maser transitions in the frequency range 0-1.91THz, with particular emphasis on maser transitions that may be generated in evolved-star envelopes and observed with the ALMA and SOFIA telescopes. We used tens of thousands of radiative transfer models of both spin species of H2O, spanning a considerable parameter space i… ▽ More We present the results of models that were designed to study all possible water maser transitions in the frequency range 0-1.91THz, with particular emphasis on maser transitions that may be generated in evolved-star envelopes and observed with the ALMA and SOFIA telescopes. We used tens of thousands of radiative transfer models of both spin species of H2O, spanning a considerable parameter space in number density, kinetic temperature and dust temperature. Results, in the form of maser optical depths, have been summarized in a master table, Table 6. Maser transitions identified in these models were grouped according to loci of inverted regions in the density/kinetic temperature plane, a property clearly related to the dominant mode of pumping. A more detailed study of the effect of dust temperature on maser optical depth enabled us to divide the maser transitions into three groups: those with both collisional and radiative pumping schemes (22,96,209,321,325,395,941 and 1486\,GHz), a much larger set that are predominantly radiatively pumped, and another large group with a predominantly collisional pump. The effect of accelerative and decelerative velocity shifts of up to 5km/s was found to be generally modest, with the primary effect of reducing computed maser optical depths. More subtle asymmetric effects, dependent on line overlap, include maximum gains offset from zero shift by >1km/s, but these effects were predominantly found under conditions of weak amplification. These models will allow astronomers to use multi-transition water maser observations to constrain physical conditions down to the size of individual masing clouds (size of a few astronomical units). △ Less

Submitted 21 October, 2015; originally announced October 2015.

SKA studies of nearby galaxies: star-formation, accretion processes and molecular gas across all environments

Authors: R. J. Beswick, E. Brinks, M. A. Perez-Torres, A. M. S. Richards, S. Aalto, A. Alberdi, M. K. Argo, I. van Bemmel, J. E. Conway, C. Dickinson, D. M. Fenech, M. D. Gray, H-R. Klockner, E. J. Murphy, T. W. B. Muxlow, M. Peel, A. P. Rushton, E. Schinnerer

Abstract: The SKA will be a transformational instrument in the study of our local Universe. In particular, by virtue of its high sensitivity (both to point sources and diffuse low surface brightness emission), angular resolution and the frequency ranges covered, the SKA will undertake a very wide range of astrophysical research in the field of nearby galaxies. By surveying vast numbers of nearby galaxies of… ▽ More The SKA will be a transformational instrument in the study of our local Universe. In particular, by virtue of its high sensitivity (both to point sources and diffuse low surface brightness emission), angular resolution and the frequency ranges covered, the SKA will undertake a very wide range of astrophysical research in the field of nearby galaxies. By surveying vast numbers of nearby galaxies of all types with $μ$Jy sensitivity and sub-arcsecond angular resolutions at radio wavelengths, the SKA will provide the cornerstone of our understanding of star-formation and accretion activity in the local Universe. In this chapter we outline the key continuum and molecular line science areas where the SKA, both during phase-1 and when it becomes the full SKA, will have a significant scientific impact. △ Less

Submitted 18 December, 2014; originally announced December 2014.

Comments: 23 pages, 4 figures, to appear as part of 'Continuum Science' in Proceedings 'Advancing Astrophysics with the Square Kilometre Array' PoS(AASKA14)070

ALMA sub-mm maser and dust distribution of VY Canis Majoris

Authors: A. M. S. Richards, C. M. V. Impellizzeri, E. M. Humphreys, C. Vlahakis, W. Vlemmings, A. Baudry, E. De Beck, L. Decin, S. Etoka, M. D. Gray, G. M. Harper, T. R. Hunter, P. Kervella, F. Kerschbaum, I. McDonald, G. Melnick, S. Muller, D. Neufeld, E. O'Gorman, S. Yu. Parfenov, A. B. Peck, H. Shinnaga, A. M. Sobolev, L. Testi, L. Uscanga , et al. (3 additional authors not shown)

Abstract: Cool, evolved stars have copious, enriched winds. The structure of these winds and the way they are accelerated is not well known. We need to improve our understanding by studying the dynamics from the pulsating stellar surface to about 10 stellar radii, where radiation pressure on dust is fully effective. Some red supergiants have highly asymmetric nebulae, implicating additional forces. We retri… ▽ More Cool, evolved stars have copious, enriched winds. The structure of these winds and the way they are accelerated is not well known. We need to improve our understanding by studying the dynamics from the pulsating stellar surface to about 10 stellar radii, where radiation pressure on dust is fully effective. Some red supergiants have highly asymmetric nebulae, implicating additional forces. We retrieved ALMA Science Verification data providing images of sub-mm line and continuum emission from VY CMa. This enables us to locate water masers with milli-arcsec precision and resolve the dusty continuum. The 658-, 321- and 325-GHz masers lie in irregular, thick shells at increasing distances from the centre of expansion. For the first time this is confirmed as the stellar position, coinciding with a compact peak offset to the NW of the brightest continuum emission. The maser shells (and dust formation zone) overlap but avoid each other on tens-au scales. Their distribution is broadly consistent with excitation models but the conditions and kinematics appear to be complicated by wind collisions, clumping and asymmetries. △ Less

Submitted 5 November, 2014; v1 submitted 18 September, 2014; originally announced September 2014.

Comments: Letter 4 pages, 5 figures plus appendix with 3 figures. Accepted by Astronomy and Astrophysics Letters

Journal ref: A&A 572, L9 (2014)

A quantum mechanical approach to establishing the magnetic field orientation from a maser Zeeman profile

Authors: J. A. Green, M. D. Gray, T. Robishaw, J. L. Caswell, N. M. McClure-Griffiths

Abstract: Recent comparisons of magnetic field directions derived from maser Zeeman splitting with those derived from continuum source rotation measures have prompted new analysis of the propagation of the Zeeman split components, and the inferred field orientation. In order to do this, we first review differing electric field polarization conventions used in past studies. With these clearly and consistentl… ▽ More Recent comparisons of magnetic field directions derived from maser Zeeman splitting with those derived from continuum source rotation measures have prompted new analysis of the propagation of the Zeeman split components, and the inferred field orientation. In order to do this, we first review differing electric field polarization conventions used in past studies. With these clearly and consistently defined, we then show that for a given Zeeman splitting spectrum, the magnetic field direction is fully determined and predictable on theoretical grounds: when a magnetic field is oriented away from the observer, the left-hand circular polarization is observed at higher frequency and the right-hand polarization at lower frequency. This is consistent with classical Lorentzian derivations. The consequent interpretation of recent measurements then raises the possibility of a reversal between the large-scale field (traced by rotation measures) and the small-scale field (traced by maser Zeeman splitting). △ Less

Submitted 23 April, 2014; originally announced April 2014.

Comments: 10 pages, 5 Figures, accepted for publication in MNRAS

Silicon isotopic abundance toward evolved stars and its application for presolar grains

Authors: T. -C. Peng, E. M. L. Humphreys, L. Testi, A. Baudry, M. Wittkowski, M. G. Rawlings, I. de Gregorio-Monsalvo, W. Vlemmings, L. -A. Nyman, M. D. Gray, C. de Breuck

Abstract: Galactic chemical evolution (GCE) is important for understanding the composition of the present-day interstellar medium (ISM) and of our solar system. In this paper, we aim to track the GCE by using the 29Si/30Si ratios in evolved stars and tentatively relate this to presolar grain composition. We used the APEX telescope to detect thermal SiO isotopologue emission toward four oxygen-rich M-type st… ▽ More Galactic chemical evolution (GCE) is important for understanding the composition of the present-day interstellar medium (ISM) and of our solar system. In this paper, we aim to track the GCE by using the 29Si/30Si ratios in evolved stars and tentatively relate this to presolar grain composition. We used the APEX telescope to detect thermal SiO isotopologue emission toward four oxygen-rich M-type stars. Together with the data retrieved from the Herschel science archive and from the literature, we were able to obtain the 29Si/30Si ratios for a total of 15 evolved stars inferred from their optically thin 29SiO and 30SiO emission. These stars cover a range of masses and ages, and because they do not significantly alter 29Si/30Si during their lifetimes, they provide excellent probes of the ISM metallicity (or 29Si/30Si ratio) as a function of time. The 29Si/30Si ratios inferred from the thermal SiO emission tend to be lower toward low-mass oxygen-rich stars (e.g., down to about unity for W Hya), and close to an interstellar or solar value of 1.5 for the higher-mass carbon star IRC+10216 and two red supergiants. There is a tentative correlation between the 29Si/30Si ratios and the mass-loss rates of evolved stars, where we take the mass-loss rate as a proxy for the initial stellar mass or current stellar age. This is consistent with the different abundance ratios found in presolar grains. We found that older objects (up to possibly 10 Gyr old) in our sample trace a previous, lower 29Si/30Si value of about 1. Material with this isotopic ratio is present in two subclasses of presolar grains, providing independent evidence of the lower ratio. Therefore, the 29Si/30Si ratio derived from the SiO emission of evolved stars is a useful diagnostic tool for the study of the GCE and presolar grains. △ Less

Submitted 14 November, 2013; originally announced November 2013.

Comments: 7 pages, 4 figures

Journal ref: A&A, 2013, 559, L8

e-MERLIN resolves Betelgeuse at wavelength 5 cm

Authors: A. M. S. Richards, R. J. Davis, L. Decin, S. Etoka, G. M. Harper, J. J. Lim, S. T. Garrington, M . D. Gray, I. McDonald, E. O'Gorman, M. Wittkowski

Abstract: Convection, pulsation and magnetic fields have all been suggested as mechanisms for the transport of mass and energy from the optical photosphere of red supergiants, out to the region where the stellar wind is launched. We imaged the red supergiant Betelgeuse at 0.06-0.18 arcsec resolution, using e-MERLIN at 5.5--6.0 GHz, with a sensitivity of ~0.01 mJy/beam. Most of the radio emission comes from… ▽ More Convection, pulsation and magnetic fields have all been suggested as mechanisms for the transport of mass and energy from the optical photosphere of red supergiants, out to the region where the stellar wind is launched. We imaged the red supergiant Betelgeuse at 0.06-0.18 arcsec resolution, using e-MERLIN at 5.5--6.0 GHz, with a sensitivity of ~0.01 mJy/beam. Most of the radio emission comes from within an ellipse (0.235x0.218) arcsec^2 (~5x the optical radius), with a flux density of 1.62 mJy, giving an average brightness temperature ~1250 K. This radio photosphere contains two hotspots of 0.53 and 0.79 mJy/beam, separated by 90 milli-arcsec, with brightness temperatures 5400+/-600 K and 3800+/-500 K. Similar hotspots, at more than double the distance from the photosphere of those seen in any other regime, were detected by the less-sensitive `old' MERLIN in 1992, 1995 and 1996 and many exceed the photospheric temperature of 3600 K. Such brightness temperatures are high enough to emanate from pockets of chromospheric plasma. Other possibilities include local shock heating, the convective dredge-up of hot material or exceptionally cool, low density regions, transparent down to the hottest layer at ~40 milliarcsec radius. We also detect an arc 0.2--0.3 arcsec to the SW, brightness temperature ~150 K, in a similar direction to extensions seen on both smaller and larger scales in the infra-red and in CO at mm wavelengths. These preliminary results will be followed by further e-MERLIN, VLA and ALMA observations to help resolve the problem of mass elevation from 1 to 10 R* in red supergiants. △ Less

Submitted 12 March, 2013; originally announced March 2013.

Comments: 6pp (incl. appendix), 3 figs Accepted by MNRAS letters 2013 Mar 12

Polarization morphology of SiO masers in the circumstellar envelope of the AGB star R Cassiopeiae

Authors: K. A. Assaf, P. J. Diamond, A. M. S. Richards, M. D. Gray

Abstract: Silicon monoxide maser emission has been detected in the circumstellar envelopes of many evolved stars in various vibrationally-excited rotational transitions. It is considered a good tracer of the wind dynamics close to the photosphere of the star. We have investigated the polarization morphology in the circumstellar envelope of an AGB star, R Cas. We mapped the linear and circular polarization o… ▽ More Silicon monoxide maser emission has been detected in the circumstellar envelopes of many evolved stars in various vibrationally-excited rotational transitions. It is considered a good tracer of the wind dynamics close to the photosphere of the star. We have investigated the polarization morphology in the circumstellar envelope of an AGB star, R Cas. We mapped the linear and circular polarization of SiO masers in the v=1, J=1-0 transition. The linear polarization is typically a few tens of percent while the circular polarization is a few percent. The fractional polarization tends to be higher for emission of lower total intensity. We found that, in some isolated features the fractional linear polarization appears to exceed 100%. We found the Faraday rotation is not negligible but is ~15 deg., which could produce small scale structure in polarized emission whilst total intensity is smoother and partly resolved out. The polarization angles vary considerably from feature to feature but there is a tendency to favour the directions parallel or perpendicular to the radial direction with respect to the star. In some features, the polarization angle abruptly flips 90 deg. We found that our data are in the regime where the model of Goldreich et al (1973) can be applied and the polarization angle flip is caused when the magnetic field is at close to 55 deg. to the line of sight. The polarization angle configuration is consistent with a radial magnetic field although other configurations are not excluded. △ Less

Submitted 11 February, 2013; originally announced February 2013.

Comments: 14 pages, 15 figures. Accepted for publication in MNRAS

Evolved star water maser cloud size determined by star size

Authors: A. M. S. Richards, S. Etoka, M. D. Gray, E. E. Lekht, J. E. Mendoza-Torres, K. Murakawa, G. Rudnitskij, J. A. Yates

Abstract: Cool, evolved stars undergo copious mass loss but the details of how the matter is returned to the ISM are still under debate. We investigated the structure and evolution of the wind at 5 to 50 stellar radii from Asymptotic Giant Branch and Red Supergiant stars. 22-GHz water masers around seven evolved stars were imaged using MERLIN, at sub-AU resolution. Each source was observed at between 2 and… ▽ More Cool, evolved stars undergo copious mass loss but the details of how the matter is returned to the ISM are still under debate. We investigated the structure and evolution of the wind at 5 to 50 stellar radii from Asymptotic Giant Branch and Red Supergiant stars. 22-GHz water masers around seven evolved stars were imaged using MERLIN, at sub-AU resolution. Each source was observed at between 2 and 7 epochs (several stellar periods). We compared our results with long-term Pushchino single dish monitoring. The 22-GHz emission is located in ~spherical, thick, unevenly filled shells. The outflow velocity doubles between the inner and outer shell limits. Water maser clumps could be matched at successive epochs separated by <2 years for AGB stars, or at least 5 years for RSG. This is much shorter than the decades taken for the wind to cross the maser shell, and comparison with spectral monitoring shows that some features fade and reappear. In 5 sources, most of the matched features brighten or dim in concert from one epoch to the next. One cloud in W Hya was caught in the act of passing in front of a background cloud leading to 50-fold, transient amplification. The masing clouds are 1-2 orders of magnitude denser than the wind average and contain a substantial fraction of the mass loss in this region, with a filling factor <1%. The RSG clouds are ~10x bigger than those round the AGB stars. Proper motions are dominated by expansion, with no systematic rotation. The maser clouds survive for decades (the shell crossing time) but the masers are not always beamed in our direction. Radiative effects cause changes in flux density throughout the maser shells on short timescales. Cloud size is proportional to parent star size; clouds have a similar radius to the star in the 22-GHz maser shell. Stellar properties such as convection cells must determine the clumping scale. △ Less

Submitted 11 July, 2012; originally announced July 2012.

Comments: Accepted by A&A 2012 July 10 Main text 29 pages, 62 figures Appendix 44 pages, 23 figures

Radio and IR interferometry of SiO maser stars

Authors: M. Wittkowski, D. A. Boboltz, M. D. Gray, E. M. L. Humphreys, I. Karovicova, M. Scholz

Abstract: Radio and infrared interferometry of SiO maser stars provide complementary information on the atmosphere and circumstellar environment at comparable spatial resolution. Here, we present the latest results on the atmospheric structure and the dust condensation region of AGB stars based on our recent infrared spectro-interferometric observations, which represent the environment of SiO masers. We dis… ▽ More Radio and infrared interferometry of SiO maser stars provide complementary information on the atmosphere and circumstellar environment at comparable spatial resolution. Here, we present the latest results on the atmospheric structure and the dust condensation region of AGB stars based on our recent infrared spectro-interferometric observations, which represent the environment of SiO masers. We discuss, as an example, new results from simultaneous VLTI and VLBA observations of the Mira variable AGB star R Cnc, including VLTI near- and mid-infrared interferometry, as well as VLBA observations of the SiO maser emission toward this source. We present preliminary results from a monitoring campaign of high-frequency SiO maser emission toward evolved stars obtained with the APEX telescope, which also serves as a precursor of ALMA images of the SiO emitting region. We speculate that large-scale long-period chaotic motion in the extended molecular atmosphere may be the physical reason for observed deviations from point symmetry of atmospheric molecular layers, and for the observed erratic variability of high-frequency SiO maser emission △ Less

Submitted 16 April, 2012; originally announced April 2012.

Comments: 8 pages, 4 figures, submitted to Proc. IAU Symp. 287 "Cosmic masers - from OH to H_0", R.S. Booth, E.M.L. Humphreys, W.H.T. Vlemmings (eds.), invited paper

Methanol and excited OH masers towards W51: I - Main and South

Authors: Sandra Etoka, Malcolm D. Gray, Gary A. Fuller

Abstract: MERLIN phase-referenced polarimetric observations towards the W51 complex were carried out in March 2006 in the Class II methanol maser transition at 6.668 GHz and three of the four excited OH maser hyperfine transitions at 6 GHz. Methanol maser emission is found towards both W51 Main and South. We did not detect any emission in the excited OH maser lines at 6.030 and 6.049 GHz down to a 3 sigma l… ▽ More MERLIN phase-referenced polarimetric observations towards the W51 complex were carried out in March 2006 in the Class II methanol maser transition at 6.668 GHz and three of the four excited OH maser hyperfine transitions at 6 GHz. Methanol maser emission is found towards both W51 Main and South. We did not detect any emission in the excited OH maser lines at 6.030 and 6.049 GHz down to a 3 sigma limit of ~20 mJy per beam. Excited OH maser emission at 6.035 GHz is only found towards W51 Main. This emission is highly circularly polarised (typically 45% and up to 87%). Seven Zeeman pairs were identified in this transition, one of which contains detectable linear polarisation. The magnetic field strength derived from these Zeeman pairs ranges from +1.6 to +6.8 mG, consistent with the previously published magnetic field strengths inferred from the OH ground-state lines. The bulk of the methanol maser emission is associated with W51 Main, sampling a total area of ~3"x2.2" (i.e., ~16200x11900 AU), while only two maser components, separated by ~2.5", are found in the W51 South region. The astrometric distributions of both 6.668-GHz methanol and 6.035-GHz excited-OH maser emission in the W51 Main/South region are presented here. The methanol masers in W51 Main show a clear coherent velocity and spatial structure with the bulk of the maser components distributed into 2 regions showing a similar conical opening angle with of a central velocity of ~+55.5 km/s and an expansion velocity of =<5 km/s. The mass contained in this structure is estimated to be at least 22 solar masses. The location of maser emission in the two afore-mentioned lines is compared with that of previously published OH ground-state emission. Association with the UCHII regions in the W51 Main/South complex and relationship of the masers to infall or outflow in the region are discussed. △ Less

Submitted 12 March, 2012; originally announced March 2012.

Comments: 19 pages, 16 figures and 4 tables, accepted for publication in MNRAS

Radiation Transfer in the Cavity and Shell of Planetary Nebulae

Authors: M. D. Gray, M. Matsuura, A. A. Zijlstra

Abstract: We develop an approximate analytical solution for the transfer of line-averaged radiation in the hydrogen recombination lines for the ionized cavity and molecular shell of a spherically symmetric planetary nebula. The scattering problem is treated as a perturbation, using a mean intensity derived from a scattering-free solution. The analytical function was fitted to Halpha and Hbeta data from the… ▽ More We develop an approximate analytical solution for the transfer of line-averaged radiation in the hydrogen recombination lines for the ionized cavity and molecular shell of a spherically symmetric planetary nebula. The scattering problem is treated as a perturbation, using a mean intensity derived from a scattering-free solution. The analytical function was fitted to Halpha and Hbeta data from the planetary nebula NGC6537. The position of the maximum in the intensity profile produced consistent values for the radius of the cavity as a fraction of the radius of the dusty nebula: 0.21 for Halpha and 0.20 for Hbeta. Recovered optical depths were broadly consistent with observed optical extinction in the nebula, but the range of fit parameters in this case is evidence for a clumpy distribution of dust. △ Less

Submitted 18 January, 2012; originally announced January 2012.

Comments: MNRAS accepted; 10 Figs

The 6-GHz methanol multibeam maser catalogue IV: Galactic longitudes 186 to 330 including the Orion-Monoceros region

Authors: J. A. Green, J. L. Caswell, G. A. Fuller, A. Avison, S. L. Breen, S. P. Ellingsen, M. D. Gray, M. Pestalozzi, L. Quinn, M. A. Thompson, M. A. Voronkov

Abstract: We present the fourth portion of a Galactic Plane survey of methanol masers at 6668 MHz, spanning the longitude range 186 degrees to 330 degrees. We report 207 maser detections, 89 new to the survey. This completes the southern sky part of the Methanol Multibeam survey and includes a large proportion of new sources, 43%. We also include results from blind observations of the Orion-Monoceros star f… ▽ More We present the fourth portion of a Galactic Plane survey of methanol masers at 6668 MHz, spanning the longitude range 186 degrees to 330 degrees. We report 207 maser detections, 89 new to the survey. This completes the southern sky part of the Methanol Multibeam survey and includes a large proportion of new sources, 43%. We also include results from blind observations of the Orion-Monoceros star forming region, formally outside the latitude range of the Methanol Multibeam survey; only the four previously known methanol emitting sites were detected, of which we present new positions and spectra for masers at Orion-A (south) and Orion-B, obtained with the MERLIN array. △ Less

Submitted 3 January, 2012; originally announced January 2012.

Comments: 39 Pages, 8 Figures, accepted for publication in MNRAS

The 43GHz SiO maser in the circumstellar envelope of the AGB star R Cassiopeiae

Authors: K. A. Assaf, P. J. Diamond, A. M. S. Richards, M. D. Gray

Abstract: We present multi-epoch, total intensity, high-resolution images of 43GHz, v=1, J=1-0 SiO maser emission toward the Mira variable R Cas. In total we have 23 epochs of data for R Cas at approximate monthly intervals over an optical pulsation phase range from 0.158 to 1.78. These maps show a ring-like distribution of the maser features in a shell, which is assumed to be centred on the star at a radiu… ▽ More We present multi-epoch, total intensity, high-resolution images of 43GHz, v=1, J=1-0 SiO maser emission toward the Mira variable R Cas. In total we have 23 epochs of data for R Cas at approximate monthly intervals over an optical pulsation phase range from 0.158 to 1.78. These maps show a ring-like distribution of the maser features in a shell, which is assumed to be centred on the star at a radius of 1.6 to 2.3 times the stellar radii. It is clear from these images that the maser emission is significantly extended around the star. At some epochs a faint outer arc can be seen at 2.2 stellar radii. The intensity of the emission waxes and wanes during the stellar phase. Some maser features are seen infalling as well as outflowing. We have made initial comparisons of our data with models by Gray et. al. (2009). △ Less

Submitted 5 May, 2011; originally announced May 2011.

Comments: 12 pages, 14 figures

Major structures of the inner Galaxy delineated by 6.7-GHz methanol masers

Authors: J. A. Green, J. L. Caswell, N. M. McClure-Griffiths, A. Avison, S. L. Breen, M. G. Burton, S. P. Ellingsen, G. A. Fuller, M. D. Gray, M. Pestalozzi, M. A. Thompson, M. A. Voronkov

Abstract: We explore the longitude-velocity distribution of 6.7-GHz methanol masers in the context of the inner structure of our Galaxy. We analyse the correlation in velocities within this distribution and identify density enhancements indicating large-scale regions of enhanced star formation. These are interpreted as the starting points of the spiral arms and the interaction of the Galactic bar with the 3… ▽ More We explore the longitude-velocity distribution of 6.7-GHz methanol masers in the context of the inner structure of our Galaxy. We analyse the correlation in velocities within this distribution and identify density enhancements indicating large-scale regions of enhanced star formation. These are interpreted as the starting points of the spiral arms and the interaction of the Galactic bar with the 3-kpc arms. The methanol masers support the presence of a long thin bar with a 45 degree orientation. Signatures of the full 3-kpc arm structure are seen, including a prominent tangent at approximately -22 degrees Galactic longitude. We compare this distribution with existing models of the gas dynamics of our Galaxy. The 3-kpc arm structure appears likely to correspond to the radius of corotation resonance of the bar, with the bar on its inner surface and the starting points of the spiral arms on its outer surface. △ Less

Submitted 20 March, 2011; originally announced March 2011.

Comments: 35 pages, 11 figures. Accepted for publication in ApJ

The 6-GHz methanol multibeam maser catalogue II: Galactic longitudes 6 to 20

Authors: J. A. Green, J. L. Caswell, G. A. Fuller, A. Avison, S. L. Breen, S. P. Ellingsen, M. D. Gray, M. R. Pestalozzi, L. Quinn, M. A. Thompson, M. A. Voronkov

Abstract: We present the second portion of an unbiased survey of the Galactic plane for 6668-MHz methanol masers. This section of the survey spans the longitude range 6 degrees to 20 degrees. We report the detection of 119 maser sources, of which 42 are new discoveries. The masers are tightly constrained to the Galactic plane, with only four outside a latitude range of +/- 1 degree. This longitude region in… ▽ More We present the second portion of an unbiased survey of the Galactic plane for 6668-MHz methanol masers. This section of the survey spans the longitude range 6 degrees to 20 degrees. We report the detection of 119 maser sources, of which 42 are new discoveries. The masers are tightly constrained to the Galactic plane, with only four outside a latitude range of +/- 1 degree. This longitude region includes the brightest known 6668-MHz methanol maser, 9.621+0.196, as well as the two brightest newly discovered sources in the southern survey as a whole. We list all the sources associated with the 3-kpc arms within +/- 15 degrees longitude and consider further candidates beyond 15 degrees longitude. We identify three new sources associated with the Galactic bar and comment on the density of masers in relation to the bar orientation. △ Less

Submitted 18 July, 2010; originally announced July 2010.

Comments: 25 pages, 4 figures, accepted for publication in MNRAS

The 6-GHz methanol multibeam maser catalogue I: Galactic Centre region, longitudes 345 to 6

Authors: J. L. Caswell, G. A. Fuller, J. A. Green, A. Avison, S. L. Breen, K. J. Brooks, M. G. Burton, A. Chrysostomou, J. Cox, P. J. Diamond, S. P. Ellingsen, M. D. Gray, M. G. Hoare, M. R. W. Masheder, N. M. McClure-Griffiths, M. R. Pestalozzi, C. J. Phillips, L. Quinn, M. A. Thompson, M. A. Voronkov, A. J. Walsh, D. Ward-Thompson, D. Wong-McSweeney, J. A. Yates, R. J. Cohen

Abstract: We have conducted a Galactic plane survey of methanol masers at 6668 MHz using a 7-beam receiver on the Parkes telescope. Here we present results from the first part, which provides sensitive unbiased coverage of a large region around the Galactic Centre. Details are given for 183 methanol maser sites in the longitude range 345$^{\circ}$ through the Galactic Centre to 6$^{\circ}$. Within 6… ▽ More We have conducted a Galactic plane survey of methanol masers at 6668 MHz using a 7-beam receiver on the Parkes telescope. Here we present results from the first part, which provides sensitive unbiased coverage of a large region around the Galactic Centre. Details are given for 183 methanol maser sites in the longitude range 345$^{\circ}$ through the Galactic Centre to 6$^{\circ}$. Within 6$^{\circ}$ of the Centre, we found 88 maser sites, of which more than half (48) are new discoveries. The masers are confined to a narrow Galactic latitude range, indicative of many sources at the Galactic Centre distance and beyond, and confined to a thin disk population; there is no high latitude population that might be ascribed to the Galactic Bulge. Within 2$^{\circ}$ of the Galactic Centre the maser velocities all lie between -60 and +77 \kms, a range much smaller than the 540 \kms range observed in CO. Elsewhere, the maser with highest positive velocity (+107 \kms) occurs, surprisingly, near longitude 355$^{\circ}$ and is probably attributable to the Galactic bar. The maser with the most negative velocity (-127 \kms) is near longitude 346$^{\circ}$, within the longitude-velocity locus of the near side of the `3-kpc arm'. It has the most extreme velocity of a clear population of masers associated with the near and far sides of the 3-kpc arm. Closer to the Galactic Centre the maser space density is generally low, except within 0.25 kpc of the Centre itself, the `Galactic Centre Zone', where it is 50 times higher, which is hinted at by the longitude distribution, and confirmed by the unusual velocities. △ Less

Submitted 12 February, 2010; originally announced February 2010.

Comments: 35 pages, 5 figures

Maser emission from SiO isotopologues traces the innermost 100 AU around Radio Source I in Orion BN/KL

Authors: C. Goddi, L. J. Greenhill, C. J. Chandler, E. M. L. Humphreys, L. D. Matthews, M. D. Gray

Abstract: We have used the Very Large Array (VLA) at 7 mm wavelength to image five rotational transitions (J=1-0) from three SiO isotopologues towards Orion BN/KL: 28SiO v=0,1,2; and 29SiO and 30SiO v=0. For the first time, we have mapped the 29SiO and 30SiO J=1-0 emission, established the maser nature of the emission, and confirmed association with the deeply embedded high-mass young stellar object commo… ▽ More We have used the Very Large Array (VLA) at 7 mm wavelength to image five rotational transitions (J=1-0) from three SiO isotopologues towards Orion BN/KL: 28SiO v=0,1,2; and 29SiO and 30SiO v=0. For the first time, we have mapped the 29SiO and 30SiO J=1-0 emission, established the maser nature of the emission, and confirmed association with the deeply embedded high-mass young stellar object commonly denoted radio Source I. The 28SiO v=0 maser emission shows a bipolar structure that extends over ~700 AU along a northeast-southwest axis, and we propose that it traces a bipolar outflow driven by Source I. The high-brightness isotopic SiO maser emission imaged with a ~0.2 arcsec resolution has a more compact distribution, generally similar to that of the 28SiO v=1,2 emission, and it probably traces bulk gas flows in a region of diameter <100 AU centered on Source I. On small scales of <10 AU, however, compact 29SiO/30SiO v=0 and 28SiO v=1,2 emission features may be offset from one another in position and line-of-sight velocity. From a radiative transfer analysis based on a large velocity gradient (LVG) pumping model, we derive similar temperatures and densities for the optimum excitation of both 29SiO/30SiO v=0 and 28SiO v=1,2 masers, significantly higher than required for 28SiO v=0 maser excitation. In order to account for the small-scale differences among the isotopologues (v=0) and the main species (v=1,2), follow-up radiative transfer modeling that incorporates non-local line overlap among transitions of all SiO isotopic species may be required. △ Less

Submitted 8 April, 2009; originally announced April 2009.

Comments: 10 pages, 3 figures, accepted for publication by ApJ

Journal ref: Astrophys.J.698:1165-1173,2009

An Improved Model of SiO Maser Emission in Miras

Authors: M. D. Gray, M. Wittkowski, M. Scholz, E. M. L. Humphreys, K. Ohnaka, D. Boboltz

Abstract: We describe a combined dynamic atmosphere and maser propagation model of SiO maser emission in Mira variables. This model rectifies many of the defects of an earlier model of this type, particularly in relation to the infra-red (IR) radiation field generated by dust and various wavelength-dependent, optically thick layers. Modelled masers form in rings with radii consistent with those found in V… ▽ More We describe a combined dynamic atmosphere and maser propagation model of SiO maser emission in Mira variables. This model rectifies many of the defects of an earlier model of this type, particularly in relation to the infra-red (IR) radiation field generated by dust and various wavelength-dependent, optically thick layers. Modelled masers form in rings with radii consistent with those found in VLBI observations and with earlier models. This agreement requires the adoption of a radio photosphere of radius approximately twice that of the stellar photosphere, in agreement with observations. A radio photosphere of this size renders invisible certain maser sites with high amplification at low radii, and conceals high-velocity shocks, which are absent in radio continuum observations. The SiO masers are brightest at an optical phase of 0.1 to 0.25, which is consistent with observed phase-lags. Dust can have both mild and profound effects on the maser emission. Maser rings, a shock and the optically thick layer in the SiO pumping band at 8.13\micron appear to be closely associated in three out of four phase samples. △ Less

Submitted 17 November, 2008; originally announced November 2008.

Comments: accepted by MNRAS, MN-08-1478-MJ.R1

Journal ref: 2009MNRAS.394...51G

The 6-GHz Multibeam Maser Survey I. Techniques

Authors: J. A. Green, J. L. Caswell, G. A. Fuller, A. Avison, S. L. Breen, K. Brooks, M. G. Burton, A. Chrysostomou, J. Cox, P. J. Diamond, S. P. Ellingsen, M. D. Gray, M. G. Hoare, M. R. W. Masheder, N. M. McClure-Griffiths, M. Pestalozzi, C. Phillips, L. Quinn, M. A. Thompson, M. Voronkov, A. Walsh, D. Ward-Thompson, D. Wong-McSweeney, J. A. Yates, R. J. Cohen

Abstract: A new 7-beam 6-7 GHz receiver has been built to survey the Galaxy and the Magellanic Clouds for newly forming high-mass stars that are pinpointed by strong methanol maser emission at 6668 MHz. The receiver was jointly constructed by Jodrell Bank Observatory (JBO) and the Australia Telescope National Facility (ATNF) and allows simultaneous coverage at 6668 and 6035 MHz. It was successfully commis… ▽ More A new 7-beam 6-7 GHz receiver has been built to survey the Galaxy and the Magellanic Clouds for newly forming high-mass stars that are pinpointed by strong methanol maser emission at 6668 MHz. The receiver was jointly constructed by Jodrell Bank Observatory (JBO) and the Australia Telescope National Facility (ATNF) and allows simultaneous coverage at 6668 and 6035 MHz. It was successfully commissioned at Parkes in January 2006 and is now being used to conduct the Parkes-Jodrell multibeam maser survey of the Milky Way. This will be the first systematic survey of the entire Galactic plane for masers of not only 6668-MHz methanol, but also 6035-MHz excited-state hydroxyl. The survey is two orders of magnitude faster than most previous systematic surveys and has an rms noise level of ~0.17 Jy.This paper describes the observational strategy, techniques and reduction procedures of the Galactic and Magellanic Cloud surveys, together with deeper, pointed, follow-up observations and complementary observations with other instruments. It also includes an estimate of the survey detection efficiency. The 111 days of observations with the Parkes telescope have so far yielded >800 methanol sources, of which ~350 are new discoveries. The whole project will provide the first comprehensive Galaxy-wide catalogue of 6668-MHz and 6035-MHz masers. △ Less

Submitted 29 October, 2008; originally announced October 2008.

Comments: 13 pages, 13 figures, accepted for publication by MNRAS

Multibeam Maser Survey of methanol and excited OH in the Magellanic Clouds: new detections and maser abundance estimates

Authors: J. A. Green, J. L. Caswell, G. A. Fuller, S. L. Breen, K. Brooks, M. G. Burton, A. Chrysostomou, J. Cox, P. J. Diamond, S. P. Ellingsen, M. D. Gray, M. G. Hoare, M. R. W. Masheder, N. McClure-Griffiths, M. Pestalozzi, C. Phillips, L. Quinn, M. A. Thompson, M. Voronkov, A. Walsh, D. Ward-Thompson, D. Wong-McSweeney, J. A. Yates, R. J. Cohen

Abstract: We present the results of the first complete survey of the Large and Small Magellanic Clouds for 6668-MHz methanol and 6035-MHz excited-state hydroxyl masers. In addition to the survey, higher-sensitivity targeted searches towards known star-formation regions were conducted. The observations yielded the discovery of a fourth 6668-MHz methanol maser in the Large Magellanic Cloud (LMC), found towa… ▽ More We present the results of the first complete survey of the Large and Small Magellanic Clouds for 6668-MHz methanol and 6035-MHz excited-state hydroxyl masers. In addition to the survey, higher-sensitivity targeted searches towards known star-formation regions were conducted. The observations yielded the discovery of a fourth 6668-MHz methanol maser in the Large Magellanic Cloud (LMC), found towards the star-forming region N160a, and a second 6035-MHz excited-state hydroxyl maser, found towards N157a. We have also re-observed the three previously known 6668-MHz methanol masers and the single 6035-MHz hydroxyl maser. We failed to detect emission from either transition in the Small Magellanic Cloud. All observations were initially made using the Methanol Multibeam (MMB) survey receiver on the 64-m Parkes telescope as part of the MMB project and accurate positions have been measured with the Australia Telescope Compact Array (ATCA). We compare the maser populations in the Magellanic Clouds with those of our Galaxy and discuss their implications for the relative rates of massive star-formation, heavy metal abundance, and the abundance of complex molecules. The LMC maser populations are demonstrated to be smaller than their Milky Way counterparts. Methanol masers are under-abundant by a factor of ~45, whilst hydroxyl and water masers are a factor of ~10 less abundant than our Galaxy. △ Less

Submitted 2 January, 2008; originally announced January 2008.

Comments: 10 pages, 5 figures, accepted for publication by MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.385:948-956,2008

Pumping of OH Main-Line Masers in Star-Forming Regions

Authors: M. D. Gray

Abstract: Pumping routes of masers can in principle be recovered from a small matrix of master equations, at an advanced stage of elimination, by tracing back the coefficients to a set of unmodified all-process rate coefficients, drawn from those which appeared in the original set of master equations, prior to any elimination operations. The traceback is achieved by logging the operations carried out on e… ▽ More Pumping routes of masers can in principle be recovered from a small matrix of master equations, at an advanced stage of elimination, by tracing back the coefficients to a set of unmodified all-process rate coefficients, drawn from those which appeared in the original set of master equations, prior to any elimination operations. The traceback is achieved by logging the operations carried out on each coefficient. There is no guarantee that a pumping scheme can be represented as a small set of important routes in this way. In the present work, the traceback method is applied to a model which is typical of a large volume of parameter space which produces very strong inversions in the main lines of the rotational ground state of OH, at 1665- and 1667-MHz. For both lines, the pumping scheme is largely restricted to the $^{2}Π_{3/2}$ stack of rotational levels, and it is possible to list a comparatively small set of routes (less than ten) which provide more than 80 per cent of the inversion. In both cases, the strongest, and simplest, route consists of a radiative upward stage, to the $^{2}Π_{3/2}, J=5/2$ rotational level, followed by a collisional de-excitation to the rotational ground state. △ Less

Submitted 22 November, 2006; originally announced November 2006.

Comments: 13 pages, 3 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.375:477-488,2007

OH megamasers, starburst and AGN activity in Markarian 231

Authors: A. M. S. Richards, J. H. Knapen, J. A. Yates, R. J. Cohen, J. L. Collett, M. M. Wright, M. D. Gray, D. Field

Abstract: We present MERLIN observations of OH maser and radio continuum emission from the Ultra Luminous IR Galaxy Markarian 231. The 1665- and 1667-MHz transitions have a combined velocity extent of 720 km/s and show a similar position-velocity structure including a gradient of 1.7 km/s/pc from NW to SE along the 420-pc major axis, steeper in the inner few tens of pc. The maser distribution is modelled… ▽ More We present MERLIN observations of OH maser and radio continuum emission from the Ultra Luminous IR Galaxy Markarian 231. The 1665- and 1667-MHz transitions have a combined velocity extent of 720 km/s and show a similar position-velocity structure including a gradient of 1.7 km/s/pc from NW to SE along the 420-pc major axis, steeper in the inner few tens of pc. The maser distribution is modelled as a torus rotating about an axis inclined at ~45deg. We estimate the enclosed mass density to be 320(90) Msun in a flattened distribution, including a central unresolved mass of </=8E+06 Msun. All the maser emission is projected against a region with a radio continuum brightness temperature >/=1E+05 K, giving a maser gain of </=2.2. The 1667:1665-MHz line ratio is close to the LTE ratio of 1.8 consistent with radiatively pumped, unsaturated masers. The size of individual masing regions is in the range 0.25-4 pc with a covering factor close to unity. There are no very bright compact masers, in contrast to galaxies such as the Seyfert 2 Markarian 273 where the masing torus is viewed nearer edge-on. The comparatively modest maser amplification seen from Markarian 231 is consistent with its classification as a Seyfert 1. Most of the radio continuum emission on 50-500 pc scales is probably of starburst origin but the compact peak is 0.4 per cent polarized by a magnetic field running north-south, similar to the jet direction on these scales. There is no close correlation between maser and continuum intensity. Comparisons with other data show that the jet changes direction close the nucleus and suggest that the sub-kpc disc hosting the masers and starburst activity is severely warped. △ Less

Submitted 13 September, 2005; originally announced September 2005.

Comments: 16 pages, 7 figures

Journal ref: Mon.Not.Roy.Astron.Soc.364:353-366,2005

The symmetric dust shell and the central star of the bipolar planetary nebula NGC 6537

Authors: M. Matsuura, A. A. Zijlstra, M. D. Gray, F. J. Molster, L. B. F. M. Waters

Abstract: We present high-resolution images of the strongly bipolar planetary nebula NGC 6537, obtained with Hubble Space Telescope and with the infrared adaptive optics system on the Very Large Telescope. The central star is detected for the first time. Using the multi-band photometry and constraints from the dynamical age of the nebula, we derive a temperature in the range 1.5-2.5 10^5 K, a luminosity~1… ▽ More We present high-resolution images of the strongly bipolar planetary nebula NGC 6537, obtained with Hubble Space Telescope and with the infrared adaptive optics system on the Very Large Telescope. The central star is detected for the first time. Using the multi-band photometry and constraints from the dynamical age of the nebula, we derive a temperature in the range 1.5-2.5 10^5 K, a luminosity~10^3 L_sun and a core mass M_c~0.7-0.9 M_sun. The progenitor mass is probably in the range M_i = 3-7 M_sun. The extinction map shows a largely symmetric, and compact dust structure, which is most likely a shell, located at the neck of the bipolar flow, only 4 arcsec from the star. The dust shell traces a short-lived phase of very high mass loss at the end of the AGB. The dynamical age of the shell and bipolar lobes are very similar but the morphologies are very different. The data suggests that the mass loss during the ejection of the compact shell was largely spherically symmetric, and the pronounced bipolarity formed afterwards. The dynamical ages of the bipolar lobes and dust shell are similar, which is consistent with suggestions that bipolar structures form in a run-away event at the very last stages of the AGB mass loss. The inner edge of the dust shell is ionized, and PAH emission is seen just outside the ionized gas. We associate the PAH emission with the photo-dissociation region of the molecular shell. △ Less

Submitted 27 July, 2005; originally announced July 2005.

Comments: 10 pages, accepted by MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.363:628-640,2005

Hidden bipolarity in red supergiant winds

Authors: A. M. S. Richards, M. R. W. Masheder, H. J. van Langevelde, R. J. Cohen, M. D. Gray, J. A. Yates, P. J. Diamond, W. H. T. Vlemmings, M. Szymczak, K. Murakawa

Abstract: Many observations of late-type M stars show large near-spherical circumstellar envelopes, yet planetary nebulae and supernova remnants are frequently axisymmetric. We present VLBI and MERLIN observations of masers around the red supergiant S Per which show varying degrees of axisymmetry and a dynamically significant magnetic field. There is no evidence for rotation here or in most similar object… ▽ More Many observations of late-type M stars show large near-spherical circumstellar envelopes, yet planetary nebulae and supernova remnants are frequently axisymmetric. We present VLBI and MERLIN observations of masers around the red supergiant S Per which show varying degrees of axisymmetry and a dynamically significant magnetic field. There is no evidence for rotation here or in most similar objects. We examine possible origins of the magnetic field. △ Less

Submitted 3 January, 2005; originally announced January 2005.

Comments: 4 pages. 4 figures. Proceedings of the 7th European VLBI Network Symposium held in Toledo, Spain on October 12-15, 2004. Editors: R. Bachiller, F. Colomer, J.-F. Desmurs, P. de Vicente (Observatorio Astronomico Nacional), p. 209-212. Needs evn2004.cls

MERLIN 6-GHz maser emission from W3(OH)

Authors: S. Etoka, R. J. Cohen, M. D. Gray

Abstract: We present the preliminary results of the first alignment to milliarsecond accuracy of the 6 GHz maser emission in OH and methanol maser lines toward W3(OH). The identifications of Zeeman pairs allowed us to infer the actual velocity of the OH material and the magnetic field strength at the location where 6.7 GHz methanol maser emission arises. We present the preliminary results of the first alignment to milliarsecond accuracy of the 6 GHz maser emission in OH and methanol maser lines toward W3(OH). The identifications of Zeeman pairs allowed us to infer the actual velocity of the OH material and the magnetic field strength at the location where 6.7 GHz methanol maser emission arises. △ Less

Submitted 30 December, 2004; originally announced December 2004.

Comments: 2 pages. 2 figures. Proceedings of the 7th European VLBI Network Symposium held in Toledo, Spain on October 12-15, 2004. Editors: R. Bachiller, F. Colomer, J.-F. Desmurs, P. de Vicente (Observatorio Astronomico Nacional), p. 199-200. Needs evn2004.cls

Angular diameters, fluxes and extinction of compact planetary nebulae: further evidence for steeper extinction towards the Bulge

Authors: P. M. E. Ruffle, A. A. Zijlstra, J. R. Walsh, M. D. Gray, K. Gesicki, D. Minniti, F. Comeron, .

Abstract: We present values for angular diameter, flux and extinction for 70 Galactic planetary nebulae observed using narrow band filters. Angular diameters are derived using constant emissivity shell and photoionization line emission models. The mean of the results from these two models are presented as our best estimate. Contour plots of 36 fully resolved objects are included and the low intensity cont… ▽ More We present values for angular diameter, flux and extinction for 70 Galactic planetary nebulae observed using narrow band filters. Angular diameters are derived using constant emissivity shell and photoionization line emission models. The mean of the results from these two models are presented as our best estimate. Contour plots of 36 fully resolved objects are included and the low intensity contours often reveal an elliptical structure that is not always apparent from FWHM measurements. Flux densities are determined, and for both H-alpha and O[III] there is little evidence of any systematic differences between observed and catalogued values. Observed H-alpha extinction values are determined using observed H-alpha and catalogued radio fluxes. H-alpha extinction values are also derived from catalogued H-alpha and H-beta flux values by means of an Rv dependent extinction law. Rv is then calculated in terms of observed extinction values and catalogued H-alpha and H-beta flux values. Comparing observed and catalogue extinction values for a subset of Bulge objects, observed values tend to be lower than catalogue values calculated with Rv = 3.1. For the same subset we calculate <Rv> = 2.0, confirming that toward the Bulge interstellar extinction is steeper than Rv = 3.1. For the inner Galaxy a relation with the higher supernova rate is suggested, and that the low-density warm ionized medium is the site of the anomalous extinction. Lowvalues of extinction are also derived using dust models with a turnover radius of 0.08 microns. △ Less

Submitted 8 July, 2004; originally announced July 2004.

Comments: Accepted by MNRAS. 17 pages, 9 figures (including 36 contour plots of PNe), 5 Tables (including 2 large tables of angular diameters, fluxes and extinction)

Modification of Dust Grain Structure by Sputtering

Authors: M. D. Gray, M. G. Edmunds

Abstract: We have applied the SRIM computer code to study the sputtering of some likely astrophysical grain materials, and we have shown that selective embedding of metallic projectiles offers a partial explanation of gas-phase depletions. We show that supernova shockwaves sweep a significantly larger mass of interstellar gas per unit time than the shockwaves generated by outflows in star-forming regions.… ▽ More We have applied the SRIM computer code to study the sputtering of some likely astrophysical grain materials, and we have shown that selective embedding of metallic projectiles offers a partial explanation of gas-phase depletions. We show that supernova shockwaves sweep a significantly larger mass of interstellar gas per unit time than the shockwaves generated by outflows in star-forming regions. We apply our sputtering model to the bombardment levels expected in a supernova shock, and show that net embedding may dominate over net sputtering, leading to grain growth under some circumstances, particularly when the bombarding gas is enriched with metals from the supernova progenitor star. A combination of short cooling times and net embedding mean that it is possible for a type II supernova to generate more dust that it destroys, and we conclude that, in general, the sputtering process often leads to a compositional change in the grain material rather than simply to grain erosion. △ Less

Submitted 12 December, 2003; originally announced December 2003.

Comments: tex file and nine postscript figure in tar.gz archive

Journal ref: Mon.Not.Roy.Astron.Soc. 349 (2004) 491

Numerical simulations of stellar SiO maser variability. Investigation of the effect of shocks

Authors: E. M. L. Humphreys, M. D. Gray, J. A. Yates, D. Field, G. H. Bowen, P. J. Diamond

Abstract: A stellar hydrodynamic pulsation model has been combined with a SiO maser model in an attempt to calculate the temporal variability of SiO maser emission in the circumstellar envelope (CE) of a model AGB star. This study investigates whether the variations in local physical conditions brought about by shocks are the predominant contributing factor to SiO maser variability because, in this work,… ▽ More A stellar hydrodynamic pulsation model has been combined with a SiO maser model in an attempt to calculate the temporal variability of SiO maser emission in the circumstellar envelope (CE) of a model AGB star. This study investigates whether the variations in local physical conditions brought about by shocks are the predominant contributing factor to SiO maser variability because, in this work, the radiative part of the pump is constant. We find that some aspects of the variability are not consistent with a pump provided by shock-enhanced collisions alone. In these simulations, gas parcels of relatively enhanced SiO abundance are distributed in a model CE by a Monte Carlo method, at a single epoch of the stellar cycle. From this epoch on, Lagrangian motions of individual parcels are calculated according to the velocity fields encountered in the model CE during the stellar pulsation cycle. The potentially masing gas parcels therefore experience different densities and temperatures, and have varying line-of-sight velocity gradients throughout the stellar cycle, which may or may not be suitable to produce maser emission. At each epoch (separated by 16.6 days), emission lines from the parcels are combined to produce synthetic spectra and VLBI-type images. We report here the results for v=1, J=1-0 (43-GHz) and J=2-1 (86-GHz) masers. △ Less

Submitted 22 February, 2002; originally announced February 2002.

Comments: 16 pages, 8 figures, accepted by A&A