Dense cores in the L1630 molecular cloud: discovering new protostars with SCUBA

Maps of the 450- and 850-μm dust continuum emission from three star-forming condensations within the Lynds 1630 molecular cloud, made with the SCUBA bolometer array, reveal the presence of four new submillimetre sources, each of a few solar masses (two of which are probably class I and two of which are class 0), as well as several sources the existence of which was previously known. The sources are located in filaments and appear elongated when observed at 450 μm. They probably have dust temperatures in the range 10 to 20 K, in good agreement with previous ammonia temperature estimates. Attempts to fit their structures with power-law and Gaussian density distributions suggest that the central distribution is flatter than expected for a simple singular isothermal sphere.
Although the statistics are poor, our results suggest that the ratio of 'protostellar core' mass to total virial mass may be similar for both large and small condensations.     

Study of star formation in RCW 106 using far-infrared observations

High-resolution far-infrared observations of a large area of the star-forming complex RCW 106 obtained using the TIFR 1-m balloon-borne telescope are presented. Intensity maps have been obtained simultaneously in two bands centred around 150 and 210 μm. Intensity maps have also been obtained in the four IRAS bands using HIRES-processed IRAS data. From the 150- and 210-μm maps, reliable maps of dust temperature and optical depth have been generated. The star formation in this complex has occurred in five linear sub-clumps. Using the map at 210 μm, which has a spatial resolution superior to that of IRAS at 100 μm, 23 sources have been identified. The spectral energy distribution (SED) and luminosity of these sources have been determined using the associations with the IRAS maps. The luminosity distribution of these sources has been obtained. Assuming these embedded sources to be zero-age main-sequence stars and using the mass–luminosity relation for these, the power-law slope of the initial mass function is found to be −1.73±0.5 . This index for this very young complex is about the same as that for more evolved complexes and clusters. Radiation transfer calculations in spherically symmetric geometry have been undertaken to fit the SEDs of 13 sources with fluxes in both the TIFR and the IRAS bands. From this, the r ⁻² density distribution in the envelopes is ruled out. Finally, a correlation is seen between the luminosity of embedded sources and the computed dust masses of the envelopes.     

The mid-infrared colours of Galactic bulge, disc and Magellanic planetary nebulae

We present mid-infrared (MIR) photometry for 367 Galactic disc, bulge and Large Magellanic Cloud (LMC) planetary nebulae (PNe), determined using data acquired with the Spitzer Space Telescope , and through the Legacy Programs GLIMPSE II (Galactic Legacy Infrared Mid-plane Survey Extraordinaire II) and SAGE (Surveying the Agents of the Galaxy's Evolution). This has permitted us to make a comparison between the luminosity functions of bulge and LMC PNe, and between the MIR colours of all three categories of source. It is determined that whilst the  3.6 μm  luminosity functions of the LMC and bulge sources are likely to be closely similar, the [3.6]–[5.8] and [5.8]–[8-0] indices of LMC nebulae are different from those of their disc and bulge counterparts. This may arise because of enhanced  6.2 μm  polycyclic aromatic hydrocarbon emission within the LMC sources, and/or as a result of further, and more radical differences between the spectra of LMC and Galactic PNe. We also determine that the more evolved disc sources listed in the Macquarie/AAO/Strasbourg (MASH) catalogues of Parker et al. and Miszalski et al. have similar colours to those of the less evolved (and higher surface brightness) sources in the catalogue of Acker et al., a result which appears at variance with previous studies of these sources.     

On the origin of double main-sequence turn-offs in star clusters of the Magellanic Clouds

Recent observational studies of intermediate-age star clusters (SCs) in the Large Magellanic Cloud (LMC) have reported that a significant number of these objects show double main-sequence turn-offs (DMSTOs) in their colour-magnitude diagrams (CMDs). One plausible explanation for the origin of these DMSTOs is that the SCs are composed of two different stellar populations with age differences of ∼300 Myr. Based on analytical methods and numerical simulations, we explore a new scenario in which SCs interact and merge with star-forming giant molecular clouds (GMCs) to form new composite SCs with two distinct component populations. In this new scenario, the possible age differences between the two different stellar populations responsible for the DMSTOs are due largely to secondary star formation within GMCs interacting and merging with already-existing SCs in the LMC disc. The total gas masses being converted into new stars (i.e. the second generation of stars) during GMC-SC interaction and merging can be comparable to or larger than the masses of the original SCs (i.e. the first generation of stars) in this scenario. Our simulations show that the spatial distributions of new stars in composite SCs formed from GMC-SC merging are more compact than those of stars initially in the SCs. We discuss both advantages and disadvantages of the new scenario in explaining fundamental properties of SCs with DMSTOs in the LMC and in the Small Magellanic Cloud (SMC). We also discuss the merits of various alternative scenarios for the origin of the DMSTOs.     

On the prospect of inferring the halo structure and the masses of dark objects through parallax microlensing

We study the proposed use of parallax microlensing in the direction of the Large Magellanic Cloud (LMC) to separate the effects of the mass function of dark massive halo objects (MHOs or 'machos') on the one hand, and their spatial distribution and kinematics on the other. This disentanglement is supposed to allow a much better determination of the two than could be achieved entirely on the basis of the durations of events. We restrict our treatment to the same class of power-law spherical models for the halo of MHOs studied in a previous paper by Marković 38 Sommer-Larsen, and assume that one can eliminate microlensing events caused by massive objects outside the halo (e.g., the LMC halo). Whereas the duration-based error in the average MHO mass, μ¯ ≡  M ¯/M, exceeds (at N  = 100 events) μ¯ by a factor of 2 or more, parallax microlensing remarkably brings it down to 15–20 per cent of μ¯, regardless of the shape of the mass function. In addition, the slope α of the mass function, d n /dμ ∝ μ^α, can be inferred relatively accurately (σ_α < 0.4) for a broader range, −3 < α < 0. The improvement in the inference of the halo structure is also significant: the index γ of the density profile ( ρ ∼  R ^−γ) can be obtained with the error σ_γ < 0.4. While in a typical situation the errors for the parameters specifying the velocity dispersion profile are of about the same magnitude as the parameters themselves, virtually all the uncertainty is 'concentrated' in linear combinations of the parameters that may have little influence on the profile, thus allowing its reasonably accurate inference.     

An Australia Telescope Compact Array 20-cm radio continuum study of the Large Magellanic Cloud

We present a mosaic image of the 1.4-GHz radio continuum emission from the Large Magellanic Cloud (LMC) observed with the Australia Telescope Compact Array (ATCA) and the Parkes Telescope. The mosaic covers     with an angular resolution of 40 arcsec, corresponding to a spatial scale of ∼10 pc in the LMC. The final image is suitable for studying emission on all scales between 40 arcsec and the surveyed area. In this paper, we discuss (i) the characteristics of the LMC's diffuse and compact radio continuum emission, (ii) the fraction of the emission produced by thermal processes and the implied star formation rate in the LMC and (iii) variations in the radio spectral index across the LMC. Two non-standard reduction techniques that we used to process the ATCA visibility data may be of interest for future wide-field radio continuum surveys. The data are open to the astronomical community and should be a rich resource for studies of individual objects such as supernova remnants, H  ii regions and planetary nebulae as well as extended features such as the diffuse emission from synchrotron radiation.     

Measuring angular diameters of extended sources

When measuring diameters of partially resolved sources like planetary nebulae, H  ii regions or galaxies, often a technique called Gaussian deconvolution is used. This technique yields a Gaussian diameter, which subsequently has to be multiplied by a conversion factor to obtain the true angular diameter of the source. This conversion factor is a function of the FWHM of the beam or point spread function, and also depends on the intrinsic surface brightness distribution of the source.
In this paper, conversion factors are presented for a number of simple geometries: a circular constant surface brightness disc and a spherical constant emissivity shell, using a range of values for the inner radius. Also, more realistic geometries are studied, based on a spherically symmetric photoionization model of a planetary nebula. This enables a study of optical depth effects, a comparison between images in various emission lines, and the use of power-law density distributions. It is found that the conversion factor depends quite critically on the intrinsic surface brightness distribution, which is usually unknown. The uncertainty is particularly large if extended regions of low surface brightness are present in the nebula. In such cases the use of Gaussian or second-moment deconvolution is not recommended.
As an alternative, a new algorithm is presented which allows the determination of the intrinsic FWHM of the source using only the observed surface brightness distribution and the FWHM of the beam. Hence no assumptions concerning the intrinsic surface brightness distribution are needed. Tests show that this implicit deconvolution method works well in realistic conditions, even when the signal-to-noise ratio is low, provided that the beamsize is less than roughly 2/3 of the observed FWHM and the beam profile can be approximated by a Gaussian. A code implementing this algorithm is available.     

The far-infrared–submillimetre spectral energy distribution of high-redshift quasars

We combine photometric observations of high-redshift     quasars, obtained at submillimetre to millimetre wavelengths, to obtain a mean far-infrared (rest-frame) spectral energy distribution (SED) of the thermal emission from dust, parametrized by a single temperature ( T ) and power-law emissivity index ( β ). The best-fitting values are     and     . Our method exploits the redshift spread of this set of quasars, which allows us to sample the SED at a larger number of rest wavelengths than is possible for a single object: the wavelength range extends down to ∼60 μm, and therefore samples the turnover in the greybody curve for these temperatures. This parametrization is of use to any studies that extrapolate from a flux at a single wavelength, for example to infer dust masses and far-infrared luminosities.
We interpret the cool, submillimetre component as arising from dust heated by star formation in the host galaxy of the quasar, although we do not exclude the presence of dust heated directly by the active galactic nucleus (AGN). Applying the mean SED to the data, we derive consistent star formation rates ∼1000 M_⊙ yr⁻¹ and dust masses ∼10⁹ M_⊙, and investigate a simple scheme of AGN and host galaxy co-evolution to account for these quantities. The time-scale for formation of the host galaxy is     , and the luminous quasar phase occurs towards the end of this period, just before the reservoir of cold gas is depleted. Given the youth of the Universe at     (1.6 Gyr), the coexistence of a massive black hole and a luminous starburst at high redshifts is a powerful constraint on models of quasar host galaxy formation.     

New constraints on turbulent transport in accretion discs

Many objects studied in astronomy follow a power-law distribution function (DF), for example the masses of stars or star clusters. A still used method by which such data is analysed is to generate a histogram and fit a straight line to it. The parameters obtained in this way can be severely biased, and the properties of the underlying DF, such as its shape or a possible upper limit, are difficult to extract. In this work, we review techniques available in the literature and present newly developed (effectively) bias-free estimators for the exponent and the upper limit. Furthermore, we discuss various graphical representations of the data and powerful goodness-of-fit tests to assess the validity of a power law for describing the distribution of data. As an example, we apply the presented methods to the data set of massive stars in R136 and the young star clusters in the Large Magellanic Cloud. For R136 we confirm the result of Koen of a truncated power law with a bias-free estimate for the exponent of  2.20 ± 0.78/2.87 ± 0.98  (where the Salpeter–Massey value is 2.35) and for the upper limit of  143 ± 9/163 ± 9 M_⊙  , depending on the stellar models used. The star clusters in the Large Magellanic Cloud (with ages up to  10^7.5 yr  ) follow a truncated power-law distribution with exponent  1.62 ± 0.06  and upper limit  68 ± 12 × 10³ M_⊙  . Using the graphical data representation, a significant change in the form of the mass function below  10^2.5 M_⊙  can be detected, which is likely caused by incompleteness in the data.     

The initial mass function of the rich young cluster NGC 1818 in the Large Magellanic Cloud

We use deep Hubble Space Telescope photometry of the rich, young (∼20- to 45-Myr old) star cluster NGC 1818 in the Large Magellanic Cloud to derive its stellar mass function (MF) down to  ∼0.15 M_⊙  . This represents the deepest robust MF thus far obtained for a stellar system in an extragalactic, low-metallicity  ([Fe/H]≃−0.4 dex)  environment. Combining our results with the published MF for masses above  1.0 M_⊙  , we obtain a complete present-day MF. This is a good representation of the cluster's initial MF (IMF), particularly at low masses, because our observations are centred on the cluster's uncrowded half-mass radius. Therefore, stellar and dynamical evolution of the cluster will not have affected the low-mass stars significantly. The NGC 1818 IMF is well described by both a lognormal and a broken power-law distribution with slopes of  Γ= 0.46 ± 0.10  and  Γ≃−1.35  (Salpeter-like) for masses in the range from 0.15 to  0.8 M_⊙  and greater than  0.8 M_⊙  , respectively. Within the uncertainties, the NGC 1818 IMF is fully consistent with both the Kroupa solar neighbourhood and the Chabrier lognormal mass distributions.     

Higher order moments of the density field in a parametrized sequence of non-Gaussian theories

We calculate the higher order moments in a sequence of models where the initial density fluctuations are drawn from a     distribution with a power-law power spectrum. For large values of , the distribution is approximately Gaussian, and we reproduce the values known from perturbation theory. As is lowered the distribution becomes progressively more non-Gaussian, approximating models with rare, high-amplitude peaks. The limit   =1  is a realization of recently proposed isocurvature models for producing early galaxy formation, where the density perturbations are quadratic in a Gaussian field.     

Hyperfine structure and interstellar curves of growth

The effects of the hyperfine structure (HFS) that is present in some interstellar absorption lines are investigated in the case of a single absorbing cloud. If the respective total equivalent widths of two or more unresolved HFS multiplets measured in relatively low-resolution spectra are analysed specifically by means of a curve of growth, the column density N (X) and the linewidth parameter b (X) inferred for absorbing species X in the cloud will generally be in error if the HFS is ignored. The fundamental physical effect is the reduced line saturation that arises because the total column density is divided among the HFS levels of the ground atomic level, each of which generally gives rise to an HFS line at a different wavelength. For nuclear spins   I = 3/2  and   I = 5/2  , theoretical curves of growth are calculated for some of the resonance lines of some alkali atoms, for each of four illustrative choices of the parameter  α=Δ/ b   , the ratio of the HFS splitting in the ground atomic level to the linewidth. Applications of the results to interstellar absorption by Na  i , K  i and Al  iii are emphasized. HFS is, fortunately, unimportant for most interstellar lines, however. Among the 35 elements that have been detected in diffuse clouds via interstellar absorption in the ultraviolet/optical spectral region, the most abundant isotope of each of 25 of these shows no HFS, because   I = 0  or, in the relevant ground atomic level,   J = 0  .     

Long-term scintillation observations of five pulsars at 1540 MHz

From 2001 January to 2002 June, we monitored PSRs B0329+54, B0823+26, B1929+10, B2020+28 and B2021+51 using the Nanshan 25-m radio telescope of the Urumqi Observatory to study their diffractive interstellar scintillation (DISS). The average interval between observations was about 9 d and the observation duration ranged between 2 and 6 h depending on the pulsar. Wide variations in the DISS parameters were observed over the 18-month data span. Despite this, the average scintillation velocities are in excellent agreement with the proper motion velocities. The average two-dimensional autocorrelation function for PSR B0329+54 is well described by a thin-screen Kolmogorov model, at least along the time and frequency axes. Observed modulation indices for the DISS time-scale and bandwidth and the pulsar flux density are greater than values predicted for a Kolmogorov spectrum of electron density fluctuations. Correlated variations over times that are long compared to the nominal refractive scintillation time are observed, suggesting that larger scale density fluctuations are important. For these pulsars, the scintillation bandwidth as a function of frequency has a power-law index  (∼3.6)  much less than that expected for Kolmogorov turbulence (∼4.4). Sloping fringes are commonly observed in the dynamic spectra, especially for PSR B0329+54. The detected range of fringe slopes are limited by our observing resolution. Our observations are sensitive to larger-scale fringes and hence smaller refractive angles, corresponding to the central part of the scattering disc.     

Period–luminosity relations for type II Cepheids and their application

JHK _s magnitudes corrected to mean intensity are estimated for Large Magellanic Cloud (LMC) type II Cepheids in the OGLE-III survey the third phase of the Optical Gravitational Lensing Experiment (OGLE). Period–luminosity (PL) relations are derived in JHK _s as well as in a reddening-free VI parameter. Within the uncertainties, the BL Her stars  ( P < 4 d)  and the W Vir stars (   P = 4  to 20 d) are colinear in these PL relations. The slopes of the infrared relations agree with those found previously for type II Cepheids in globular clusters within the uncertainties. Using the pulsation parallaxes of V553 Cen and SW Tau, the data lead to an LMC modulus uncorrected for any metallicity effects of  18.46 ± 0.10  mag. The type II Cepheids in the second-parameter globular cluster, NGC 6441, show a PL( VI ) relation of the same slope as that in the LMC, and this leads to a cluster distance modulus of  15.46 ± 0.11  mag, confirming the hypothesis that the RR Lyrae variables in this cluster are overluminous for their metallicity. It is suggested that the Galactic variable κ Pavonis is a member of the peculiar W Vir class found by the OGLE-III group in the LMC. Low-resolution spectra of OGLE-III type II Cepheids with   P > 20  d (RV Tau stars) show that a high proportion have TiO bands; only one has been found showing C₂. The LMC RV Tau stars, as a group, are not colinear with the shorter period type II Cepheids in the infrared PL relations in marked contrast to such stars in globular clusters. Other differences between LMC, globular cluster and Galactic field type II Cepheids are noted in period distribution and infrared colours.     

Properties of warm absorbers in active galaxies: a systematic stability curve analysis

Signatures of warm absorbers are seen in soft X-ray spectra of about half of all type 1 Seyfert galaxies observed and in some quasars and blazars. We use the thermal equilibrium curve to study the influence of the shape of the ionizing continuum, density and the chemical composition of the absorbing gas on the existence and nature of the warm absorbers. We describe circumstances in which a stable warm absorber can exist as a multiphase medium or one with continuous variation in pressure. In particular, we find the following results: (i) the warm absorber exists only if the spectral index of the X-ray power-law ionizing continuum  α > 0.2  and has a multiphase nature if  α∼ 0.8  , which interestingly is the spectral index for most of the observed type 1 Seyfert galaxies; (ii) thermal and ionization states of highly dense warm absorbers are sensitive to their density if the ionizing continuum is sufficiently soft, i.e. dominated by the ultraviolet; (iii) absorbing gas with super-solar metallicity is more likely to have a multiphase nature and (iv) the nature of the warm absorber is significantly influenced by the absence of iron and associated elements which are produced in the later stages of star formation history in Type Ia supernovae.     

Automated analysis of eclipsing binary light curves – II. Statistical analysis of OGLE LMC eclipsing binaries

In the first paper of this series, we presented EBAS – Eclipsing Binary Automated Solver, a new fully automated algorithm to analyse the light curves of eclipsing binaries, based on the ebop code. Here, we apply the new algorithm to the whole sample of 2580 binaries found in the Optical Gravitational Lensing Experiment (OGLE) Large Magellanic Cloud (LMC) photometric survey and derive the orbital elements for 1931 systems. To obtain the statistical properties of the short-period binaries of the LMC, we construct a well-defined subsample of 938 eclipsing binaries with main-sequence B-type primaries. Correcting for observational selection effects, we derive the distributions of the fractional radii of the two components and their sum, the brightness ratios and the periods of the short-period binaries. Somewhat surprisingly, the results are consistent with a flat distribution in log P between 2 and 10 d. We also estimate the total number of binaries in the LMC with the same characteristics, and not only the eclipsing binaries, to be about 5000. This figure leads us to suggest that  (0.7 ± 0.4)  per cent of the main-sequence B-type stars in the LMC are found in binaries with periods shorter than 10 d. This frequency is substantially smaller than the fraction of binaries found by small Galactic radial-velocity surveys of B stars. On the other hand, the binary frequency found by Hubble Space Telescope ( HST ) photometric searches within the late main-sequence stars of 47 Tuc is only slightly higher and still consistent with the frequency we deduced for the B stars in the LMC.     

A study of interstellar medium of dwarf galaxies using H i power spectrum analysis

We estimate the power spectrum of H  i intensity fluctuations for a sample of eight galaxies (seven dwarf and one spiral). The power spectrum can be fitted to a power-law     for six of these galaxies, indicating turbulence is operational. The estimated best-fitting value for the slope ranges from  ∼−1.5  (AND IV, NGC 628, UGC 4459 and GR 8) to  ∼−2.6  (DDO 210 and NGC 3741). We interpret this bi-modality as being due to having effectively 2D turbulence on length-scales much larger than the scale-height of the galaxy disc and 3D otherwise. This allows us to use the estimated slope to set bounds on the scale-heights of the face-on galaxies in our sample. We also find that the power-law slope remains constant as we increase the channel thickness for all these galaxies, suggesting that the fluctuations in H  i intensity are due to density fluctuations and not velocity fluctuations, or that the slope of the velocity structure function is ∼0. Finally, for the four galaxies with '2D turbulence' we find that the slope α correlates with the star formation rate (SFR) per unit area, with larger SFRs leading to steeper power laws. Given our small sample size, this result needs to be confirmed with a larger sample.     

Molecular line intensities as measures of cloud masses – I. Sensitivity of CO emissions to physical parameter variations

A reliable estimate of the molecular gas content in galaxies plays a crucial role in determining their dynamical and star-forming properties. However, H₂, the dominant molecular species, is difficult to observe directly, particularly in the regions where most molecular gas is thought to reside. Its mass is therefore commonly inferred by assuming a direct proportionality with the integrated intensity of the  ¹²CO( J = 1 → 0)  emission line, using a CO-to-H₂ conversion factor, X . Although a canonical value for X is used extensively in such estimates, there is increasing evidence, both theoretical and observational, that the conversion factor may vary by over an order of magnitude under conditions different from those of the local neighbourhood. In an effort to understand the influence of changing environmental conditions on the conversion factor, we derive theoretical estimates of X for a wide range of physical parameters using a photon-dominated region (PDR) time-dependent chemical model, benchmarking key results against those of an independent PDR code to ensure reliability. Based on these results, the sensitivity of the X factor to change in each physical parameter is interpreted in terms of the chemistry and physical processes within the cloud. In addition to confirming previous observationally derived trends, we find that the time-dependence of the chemistry, often neglected in such models, has a considerable influence on the value of the conversion factor.     

The Galactic disc distribution of planetary nebulae with warm dust emission features – I

We investigate the Galactic disc distribution of a sample of planetary nebulae characterized in terms of their mid-infrared spectral features. The total number of Galactic disc PNe with 8–13 μm spectra is brought up to 74 with the inclusion of 24 new objects, the spectra of which we present for the first time. 54 PNe have clearly identified warm dust emission features, and form a sample that we use to construct the distribution of the C/O chemical balance in Galactic disc PNe. The dust emission features complement the information on the progenitor masses brought by the gas-phase N/O ratios: PNe with unidentified infrared emission bands have the highest N/O ratios, while PNe with the silicate signature have either very high N enrichment or close to none. We find a trend for a decreasing proportion of O-rich PNe towards the third and fourth Galactic quadrants. Two independent distance scales confirm that the proportion of O-rich PNe decreases from     per cent inside the solar circle to     per cent outside. PNe with warm dust are also the youngest. PNe with no warm dust are uniformly distributed in C/O and N/O ratios, and do not appear to be confined to     They also have higher 6-cm fluxes, as expected from more evolved PNe. We show that the IRAS fluxes are a good representation of the bolometric flux for compact and IR-bright PNe, which are probably optically thick. Selection of objects with     should probe a good portion of the Galactic disc for these young, dense and compact nebulae, and the dominant selection effects are rooted in the PN catalogues.