A Search for SiO Masers in Globular Clusters Using the VLA at 43 GHz
M. Colleen Gino
Abstract
Thirty globular clusters were observed with the Very Large Array in an effort to detect SiO maser emission from possible long period variable stars in the clusters. Data analysis was distributed among 10 individuals of the Swinburne University of Technologys Introductory Radio Astronomy course HET608, and this paper reports on the results from the measurements of M13, M12 and M72. No SiO maser emission was detected from these three objects.
Introduction
This paper presents a subset of a larger study of 30 globular clusters that were observed as part of a search for SiO maser emission from evolved stars in globular clusters. The importance of this search is that the detection of masers in globular clusters would be a key contribution in constraining the various models of galactic structure and formation presently in use.
Over the past 30 years there have been a number of unsuccessful searches for masers in globular clusters (Knapp and Kerr 1973, Cohen and Malkan 1979, Dickey and Malkan 1980), and at least one partially successful search (Frail and Beasley 1994). However, a search of the literature yields no record of an attempt to detect SiO masers in globular clusters.
SiO maser emission has long been observed in evolved stars, in particular Mira variables and long period variables (Buhl et al. 1974, Kaifu et al. 1975). Such stars are known to occur in globular clusters, although their frequency appears to be low with clusters containing such stars having on average only 2-4 members per cluster. This initial search for a detection of SiO masers was performed using the Very Large Array (VLA), which consists of 27 antennas of 25 meters in diameter arranged in a Y-shaped baseline. In the A-configuration the VLA has a synthesized aperture of 27 km.
This paper describes the results of the search for SiO maser emission from 3 globular clusters, M13, M12 and M72 using the VLA.
Background
A globular cluster is a densely packed association of stars characteristically arranged in a compact sphere, as can be seen in Figure 1A. These star clusters can be found in the galactic halo, the nearly spherical aggregation of dust, gas and stars surrounding our galaxy. Approximately 150 globular clusters are currently known to be associated with

Figure 1. Image of M13 taken using an unfiltered commercial CCD camera at Mt. Wilson Observatory. Panel A shows a 9 by 6 view of the object and panel B shows the central region with the 1 VLA pointings as an overlay in magenta. For reference the VLA pointings are identified in this paper by raster order beginning in the upper left-hand corner.
the Milky Way Galaxy, and they are believed to have formed around the same time as the Galaxy itself. In addition, a number of globular clusters in the galactic halo may have been removed from other nearby galaxies and captured by the Milky Ways gravitational influence. One important feature of globular clusters is that many of the globular clusters surrounding our galaxy are known to contain Mira-type and long period variables, which, at least in some cases, produce SiO masers (i.e. WX SER).
As low to medium mass (1-7 x solar mass) stars near the end of their lives, they experience a brief period of tremendous mass loss. Stellar masers are found in the circumstellar shells produced by this loss of mass from the dying star, and occur when the turbulent upper photosphere is exposed to the intense radiation from below. SiO masers exist quite close in to the photosphere of these evolved stars, which contain an abundance of heavier elements (see http://spacesun.rice.edu/~parviz/maserspots.htm).
Since globular clusters are located in the halo of the galaxy and orbit around the center of mass of the galactic bulge, the measurement of variations in their orbits is a sensitive indicator of the total enclosed mass and the local mass distribution, which may include dark matter. The relatively close proximity of the globular clusters to our solar system allows the position of the cluster to be measured using the parallax technique, provided that the position of the globular cluster can be measured to sufficient accuracy. In the simplest case, comparing the position of the globular cluster over time with that of a distant unmoving reference point such as a quasar or other bright distant astronomical object allows variations in the clusters orbit to be detected.
While the VLA is an excellent instrument to determine the presence of SiO masers (at a rest frequency of 43 GHz) follow-up observations with the Very Long Baseline Array, a system of 10 identical 25-meter antennas spread across the United States from Hawaii to St. Croix and controlled from a central headquarters in New Mexico, allows the position of a globular cluster to be determined with a very high accuracy. Such high angular resolution measurements, typically having 0.1 milliarcsecond beamwidth (Ulvestad and Goss 1999), provides fine enough resolution to determine the angular transverse motions and, for some of the nearest clusters, obtain accurate positions using parallax.
Observations
The observations were conducted over a five-day period beginning October 12, 2000 with the NRAO Very Large Array located outside of Socorro, New Mexico. Observations took place during the move from the D-configuration to the A-configuration, the period of time when the size of the array was being changed from the smallest baseline to the largest baseline of ~27 km. The measurements for this study were insensitive to this change, making it an ideal use of the non-standard and changing array configuration. Although a total of 23 of the 27 VLA antennas are equipped with 7mm (43 GHz) receivers, the number of antennas used varied from 19 to 21 over the course of the observations due to antenna maintenance and complex array reconfiguration issues.
The observations of the globular clusters were made at 43 GHz in spectral line mode 2AC (a two-IF mode consisting of an RCP-LCP pair) using 64 channels across a 6.25 MHz bandwidth. With the primary beam size of the VLA antennas (~1) being larger than the half-mass radius of nearly the entire sample, adequate coverage was accomplished with a single pointing toward the core of each cluster. A 3 by 3 grid was then produced by phase rotating the data to search directions adjacent to the phase center, enabling the detection of maser emission around the edges of the cluster (Beasley, 2000 private communication). This resulted in a grid of 1 arcminute diameter sampled areas spaced at 30 arcsecond offsets.
In the entire project, thirty clusters were ultimately observed, with a series of 10-second integrations equaling a total of 30 minutes integration time for each cluster. A total of 20.5 hours of observations were carried out over four separate observing sessions. The data from the first 2.5-hour observing session was rendered unusable by a high percentage of moisture in the atmosphere. All but one of the clusters (Pyxis) observed during the first session were re-observed in later sessions. A second cluster (NGC6402) was not observed due to scheduling problems. Also included was an observation of the known SiO maser WX SER in order to confirm the telescopes ability to detect maser emission.
Three criteria were met in selecting 27 of the 32 clusters chosen for this survey and are itemized below:
Five additional clusters that were chosen failed to meet one or more of the above criteria, but were included in the search since they may be captured globular clusters and as such are scientifically interesting objects. The complete list of the 32 globular clusters originally intended to be observed is listed in the Appendix.
Results and Discussion
The observational results for the three globular clusters M13, M12 and M72 are presented in this paper along with the results for the known SiO maser WX SER. As this paper presents a null result for the detection of an SiO maser, the analysis for M13 will be presented in detail with the results for M12 and M72 shown only as a summary in Table 1. Table 1 lists the name, raster position, mean, standard deviation, minimum detection limit, maximum detection value, possible detection and detection confirmation status for all three globular clusters and WX SER for easy comparison.

Figure 2. Amplitude vs Velocity profiles for M13 identified by raster order, and the known SiO maser WX SER as profile 0. The profiles include a unit offset in the plot for clarity, although all are on a comparable scale. The profiles have been corrected for the M13s barycentric velocity of 246.600 and WX SERs LSR velocity of +7.000 (Beasley 2000, private communication.)
The raw data for the globular clusters consist of 64 spectral channels measured with respect to the Earths velocity. The velocities were transformed using the barycentric velocities appropriate for the astronomical object. Position information within a cluster was obtained by phase shifting the data resulting in 9 measurements for each object located on a 30 arcsecond spaced grid. The position of the measurements with respect to the center of the cluster is shown in Figure 1B and the 9 positions are identified using raster ordering.
Table 1. Summary of results for the three globular clusters and the maser WX SER.
Object |
Raster |
Mean |
Standard |
Minimum |
Maximum |
Possible |
Confirmed |
Position |
Deviation |
Detection |
Value |
Detection |
|||
Limit |
|||||||
| M13 (NGC 6205) | 1 |
0.087 |
0.051 |
0.255 |
0.198 |
- |
|
2 |
0.089 |
0.045 |
0.226 |
0.222 |
- |
||
3 |
0.085 |
0.047 |
0.236 |
0.213 |
- |
||
4 |
0.089 |
0.058 |
0.288 |
0.285 |
- |
||
5 |
0.076 |
0.051 |
0.256 |
0.237 |
- |
||
6 |
0.093 |
0.054 |
0.270 |
0.295 |
Y |
N |
|
7 |
0.081 |
0.038 |
0.192 |
0.187 |
- |
||
8 |
0.080 |
0.042 |
0.210 |
0.220 |
Y |
N |
|
9 |
0.098 |
0.054 |
0.271 |
0.294 |
Y |
N |
|
| M12 (NGC 6218) | 1 |
0.058 |
0.032 |
0.158 |
0.117 |
- |
|
2 |
0.064 |
0.038 |
0.189 |
0.163 |
- |
||
3 |
0.067 |
0.038 |
0.192 |
0.144 |
- |
||
4 |
0.058 |
0.035 |
0.173 |
0.131 |
- |
||
5 |
0.063 |
0.035 |
0.177 |
0.168 |
- |
||
6 |
0.065 |
0.035 |
0.174 |
0.153 |
- |
||
7 |
0.067 |
0.036 |
0.178 |
0.180 |
Y |
N |
|
8 |
0.064 |
0.039 |
0.196 |
0.212 |
Y |
N |
|
9 |
0.060 |
0.036 |
0.179 |
0.173 |
- |
||
| M72 (NGC 6981) | 1 |
0.048 |
0.028 |
0.141 |
0.170 |
Y |
N |
2 |
0.049 |
0.025 |
0.125 |
0.118 |
- |
||
3 |
0.046 |
0.029 |
0.143 |
0.111 |
- |
||
4 |
0.053 |
0.027 |
0.133 |
0.125 |
- |
||
5 |
0.050 |
0.028 |
0.142 |
0.145 |
Y |
N |
|
6 |
0.043 |
0.024 |
0.120 |
0.098 |
- |
||
7 |
0.043 |
0.021 |
0.107 |
0.096 |
- |
||
8 |
0.046 |
0.027 |
0.134 |
0.111 |
- |
||
9 |
0.045 |
0.026 |
0.130 |
0.114 |
- |
||
| WXSER | 1 |
0.112 |
0.075 |
0.376 |
1.646 |
Y |
Y |
2 |
0.122 |
0.141 |
0.706 |
2.700 |
Y |
Y |
|
3 |
0.131 |
0.098 |
0.490 |
2.222 |
Y |
Y |
|
4 |
0.129 |
0.096 |
0.478 |
2.022 |
Y |
Y |
|
5 |
0.182 |
0.151 |
0.754 |
6.132 |
Y |
Y |
|
6 |
0.124 |
0.076 |
0.382 |
0.847 |
Y |
Y |
|
7 |
0.112 |
0.075 |
0.376 |
1.646 |
Y |
Y |
|
8 |
0.126 |
0.089 |
0.444 |
1.904 |
Y |
Y |
|
9 |
0.117 |
0.077 |
0.387 |
0.709 |
Y |
Y |
The data analysis consisted of taking the measured amplitude vs velocity profiles for each pointing position and determining the mean and standard deviation of the mean outside of the window where an SiO maser should be seen. A five-sigma detection limit was determined for each position, and those profiles where fluctuations exceeded the five-sigma limit were flagged as possible detections for further analysis. The phase vs velocity profiles for those possible detection events were examined to see if the phase indicated detection of a radio source (see Figure 3). Only those events where the five-sigma detection threshold was exceeded and the correct phase relationship was seen were considered as confirmed detections. This evaluation criteria showed that no masers were detected in the globular clusters, but clearly identified SiO maser detection in the WX SER source.
The amplitude information for each measurement on the grid for M13 is shown in Figure 2 overlaid with the WX SER SiO maser detection. It is clear from Figure 2 that the measured amplitude of the globular cluster is significantly below that of the WX SER maser. In addition the phase information shows no detection of a maser in the data as illustrated in Figure 3. This can be translated into a lower limit of approximately 300 mJy for the detection of a maser in M13.

Figure 3. Plot of the Phase vs Velocity profile of M13 (30,0) (blue squares) and WX SER (0,0) (red diamonds). In the case of the WX SER detection, which occurs around the zero velocity, the phase stabilizes at or near zero, while in the case of the null result for M13, the phase does not.
The data from M12 and M72 were reduced in a similar fashion using 42.200 and
345.100 respectively as the barycentric corrections for these objects. Again, no masers were detected in the profiles for these globular clusters and the reduced data is presented in Table 1 for comparison with the M13 and WX SER results.
The results of the study show that in an examination of the three globular clusters M13, M12 and M72, we were unable to detect SiO maser emission using the VLA down to a lower limit of approximately 0.3, 0.2 and 0.15 Jy for M13, M12 and M72 respectively.
Acknowledgements
I would like to thank Dr. A.J. Beasley, Dr. M.J. Claussen and Dr. J. Wrobel of the National Radio Astronomy Observatory, and Dr. S.W. Teare of New Mexico Tech for the many useful discussions. In addition, I would like to express my appreciation to the National Radio Astronomy Observatory for the observing time, as well as for the opportunity to act as the Array Operator of the Very Large Array for several hours of the observations. I would also like to thank G. Clark and B. Duff for the use of the 14" TIE telescope at Mt. Wilson Observatory. Finally, I would like to thank Dr. Tony Beasley for giving me and nine other students of the Swinburne Astronomy Online Astronomy Masters Program the opportunity to participate in this research program.
References
Buhl, D., Snyder, L.E., Lovas, F.J., Johnson, D.R. 1974, ApJ 192, 97.
Cohen, N.L. & Malkan, M.A. 1979, ApJ 84, 74.
Dickey, J.M. & Malkan, M.A. 1980, AJ 85, 145.
Frail, D.A. & Beasley, A.J. 1994, A&A 290, 796.
Kaifu, N., Buhl, D., Snyder, L.E. 1975, ApJ 195, 359.
Knapp, G.R. & Kerr, F.J. 1973, AJ 78, 458.
Ulvestad, J. & Goss, M. 1999, S&T V.98 No. 6, 36.
Appendix
ID Alt ID RA(2000) Dec(2000) l b D(kpc) Fe/H V(km/s) notes
======================================================================================================
*Pyxis 09 07 57.8 -37 13 17 261.32 7.00 39.4 -1.2 36.0
NGC 4147 12 10 06.2 +18 32 31 252.85 77.19 19.3 -1.83 183.2 opt ZY
NGC 5272 M 3 13 42 11.2 +28 22 32 42.21 78.71 10.4 -1.57 -148.5 opt +LPVs ZY
Pal 5 15 16 05.3 -00 06 41 0.85 45.86 23.2 -1.43 -55.0 opt ZY
NGC 5904 M 5 15 18 33.8 +02 04 58 3.86 46.80 7.5 -1.29 52.6 opt +LPVs
NGC 6171 M 107 16 32 31.9 -13 03 13 3.37 23.01 6.4 -1.04 -33.6 opt +LPV
NGC 6205 M 13 16 41 41.5 +36 27 37 59.01 40.91 7.7 -1.54 -246.6 opt +LPVs (many)
*NGC 6229 16 46 58.9 +47 31 40 73.64 40.31 30.7 -1.43 -154.2 +LPVs?
NGC 6218 M 12 16 47 14.5 -01 56 52 15.72 26.31 4.9 -1.48 -42.2 opt
NGC 6254 M 10 16 57 08.9 -04 05 58 15.14 23.08 4.4 -1.52 75.8 opt?
NGC 6333 M 9 17 19 11.8 -18 30 59 5.54 10.70 8.2 -1.72 229.1
NGC 6342 17 21 10.2 -19 35 14 4.90 9.73 8.6 -0.65 116.2
NGC 6356 17 23 35.0 -17 48 47 6.72 10.22 15.2 -0.50 27.0 +Miras
NGC 6366 17 27 44.3 -05 04 36 18.41 16.04 3.6 -0.82 -122.3
NGC 6402 M 14 17 37 36.1 -03 14 45 21.32 14.81 8.9 -1.39 -66.1 +LPV
NGC 6517 18 01 50.6 -08 57 32 19.23 6.76 10.8 -1.37 -39.6
NGC 6535 18 03 50.7 -00 17 49 27.18 10.44 6.7 -1.80 -215.1
NGC 6539 18 04 49.8 -07 35 09 20.80 6.78 8.4 -0.66 -45.6
IC 1276 Pal 7 18 10 44.2 -07 12 27 21.83 5.67 5.4 -0.73 155.7
Pal 8 18 41 29.9 -19 49 33 14.10 -6.80 12.9 -0.48 -43.0
NGC 6712 18 53 04.3 -08 42 22 25.35 -4.32 6.9 -1.01 -107.5 opt +LPVs
NGC 6715 M 54 18 55 03.3 -30 28 42 5.61 -14.09 27.2 -1.59 141.9 +LPVs
NGC 6749 Be42 19 05 15.3 +01 54 03 36.20 -2.20 7.9 -1.60 -61.7
NGC 6760 19 11 12.1 +01 01 50 36.11 -3.92 7.4 -0.52 -27.5
Pal 10 19 18 02.1 +18 34 18 52.44 2.72 5.9 -0.10 -31.7 +LPVs
Arp 2 19 28 44.1 -30 21 14 8.55 -20.78 28.6 -1.76 115.0 ZY
Pal 11 19 45 14.4 -08 00 26 31.81 -15.58 12.9 -0.39 -68.0
NGC 6838 M 71 19 53 46.1 +18 46 42 56.74 -4.56 3.9 -0.73 -22.9 +LPVs
NGC 6934 20 34 11.6 +07 24 15 52.10 -18.89 17.4 -1.54 -411.4 ZY
NGC 6981 M 72 20 53 27.9 -12 32 13 35.16 -32.68 17.0 -1.40 -345.1 ZY
*NGC 7006 21 01 29.5 +16 11 15 63.77 -19.41 41.5 -1.63 -384.1 +LPV
NGC 7089 M 2 21 33 29.3 -00 49 23 53.38 -35.78 11.5 -1.62 -5.3 opt
======================================================================================================
Notes:
opt: optical proper motion (opt? means not a good one).
LPVs: contains identified long period variables in Christine
Clements' table of variables (Clements 2000).
ZY: Zinn "Young" type, or second parameter globular cluster.
*: possible captured LMC cluster.
| Return to the Astrophys-Assist Educational Resources Page | Copyright ©2001 C. Gino |