Testing the circumstellar disc hypothesis: a search for H2 outflow signatures from massive young stellar objects with linearly distributed methanol masers

The results of a survey searching for outflows using near-infrared imaging are presented. Targets were chosen from a compiled list of massive young stellar objects associated with methanol masers in linear distributions. Presently, it is a widely held belief that these methanol masers are found in (and delineate) circumstellar accretion discs around massive stars. If this scenario is correct, one way to test the disc hypothesis is to search for outflows perpendicular to the methanol maser distributions. The main objective of the survey was to obtain wide-field near-infrared images of the sites of linearly distributed methanol masers using a narrow-band 2.12-μm filter. This filter is centred on the  H₂ v = 1–0 S(1)  line; a shock diagnostic that has been shown to successfully trace CO outflows from young stellar objects. 28 sources in total were imaged of which 18 sources display H₂ emission. Of these, only two sources showed emission found to be dominantly perpendicular to the methanol maser distribution. Surprisingly, the H₂ emission in these fields is not distributed randomly, but instead the majority of sources are found to have H₂ emission dominantly parallel to their distribution of methanol masers. These results seriously question the hypothesis that methanol masers exist in circumstellar discs. The possibility that linearly distributed methanol masers are instead directly associated with outflows is discussed.     

Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

A Quantitative Method to Optimise Magnetic Field Line Fitting of Observed Coronal Loops

Many authors use magnetic-field models to extrapolate the field in the solar corona from magnetic data in the photosphere. The accuracy of such extrapolations is usually judged qualitatively by eye, where a less judgemental quantitative approach would be more desirable. In this paper, a robust method for obtaining the best fit between a theoretical magnetic field and intensity observations of coronal loops on the solar disk will be presented. The method will be applied to Yohkoh data using a linear force-free field as an illustration. Any other theoretical model for the magnetic field can be used, provided there is enough freedom in the model to optimize the fit.     

Correcting systematic effects in a large set of photometric light curves

We suggest a new algorithm to remove systematic effects in a large set of light curves obtained by a photometric survey. The algorithm can remove systematic effects, such as those associated with atmospheric extinction, detector efficiency, or point spread function changes over the detector. The algorithm works without any prior knowledge of the effects, as long as they linearly appear in many stars of the sample. The approach, which was originally developed to remove atmospheric extinction effects, is based on a lower rank approximation of matrices, an approach which has already been suggested and used in chemometrics, for example. The proposed algorithm is especially useful in cases where the uncertainties of the measurements are unequal. For equal uncertainties, the algorithm reduces to the Principal Component Analysis (PCA) algorithm. We present a simulation to demonstrate the effectiveness of the proposed algorithm and we point out its potential, in the search for transit candidates in particular.     

The evolution of binary fractions in globular clusters

We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code ( startrack ) and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing three- and four-body interactions ( fewbody ). We find that the combination of stellar evolution and dynamical interactions (binary–single and binary–binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster such as 47 Tuc, assuming an initial binary fraction of 100 per cent, is only ∼ 5–10 per cent. We show that this is in good agreement with recent Hubble Space Telescope observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modelling of globular clusters, suggesting that 'realistic models' should incorporate much larger initial binary fractions than has usually been the case in the past.     

Richtmyer-Meshkov Experiments on the Omega Laser

Observations of the interstellar medium reveal a dynamic realm permeated by shocks. These shocks are generated on a large range of scales by galactic rotation, supernovae, stellar winds, and other processes. Whenever a shock encounters a density interface, Richtmyer-Meshkov instabilities may develop. Perturbations along the interface grow, leading to structure formation and material mixing. An understanding of the evolution of Richtmyer-Meshkov instabilities is essential for understanding galactic structure, molecular cloud morphology, and the early stages of star formation. An ongoing experimental campaign studies Richtmyer-Meshkov mixing in a convergent, compressible, miscible plasma at the Omega laser facility. Cylindrical targets, consisting of a low density foam core and an aluminum shell covered by an epoxy ablator, are directly driven by fifty laser beams. The aluminum shell is machined to produce different perturbation spectra. Surface types include unperturbed (smooth), single-mode sinusoids, multi-mode (rough), and multi-mode with particular modes accentuated (specified-rough). Experimental results are compared to theory and numerical simulations.