Hot stars in old stellar populations: a continuing need for intermediate ages

We present numerical investigations into the formation of massive stars from turbulent cores of density structure  ρ∝ r ^−1.5  . The results of five hydrodynamical simulations are described, following the collapse of the core, fragmentation and the formation of small clusters of protostars. We generate two different initial turbulent velocity fields corresponding to power-law spectra   P ∝ k ⁻⁴  and   P ∝ k ^−3.5  , and we apply two different initial core radii. Calculations are included for both completely isothermal collapse, and a non-isothermal equation of state above a critical density  (10⁻¹⁴ g cm⁻³)  . Our calculations reveal the preference of fragmentation over monolithic star formation in turbulent cores. Fragmentation was prevalent in all the isothermal cases. Although disc fragmentation was largely suppressed in the non-isothermal runs due to the small dynamic range between the initial density and the critical density, our results show that some fragmentation still persisted. This is inconsistent with previous suggestions that turbulent cores result in the formation of a single massive star. We conclude that turbulence cannot be measured as an isotropic pressure term.     

Effects of clumping on temperature – I. Externally heated clouds

We present a study of radiative transfer in dusty, clumpy star-forming regions. A series of self-consistent, 3D, continuum radiative transfer models are constructed for a grid of models parametrized by central luminosity, filling factor, clump radius and face-averaged optical depth. The temperature distribution within the clouds is studied as a function of this parametrization. Among our results, we find that: (i) the effective optical depth in clumpy regions is less than in equivalent homogeneous regions of the same average optical depth, leading to a deeper penetration of heating radiation in clumpy clouds, and temperatures higher by over 60 per cent; (ii) penetration of radiation is driven by the fraction of open sky (FOS) – which is a measure of the fraction of solid angle along which no clumps exist; (iii) FOS increases as clump radius increases and as filling factor decreases; (iv) for values of   FOS >0.6–0.8  the sky is sufficiently 'open' that the temperature distribution is relatively insensitive to FOS; (v) the physical process by which radiation penetrates is preferentially through streaming of radiation between clumps as opposed to diffusion through clumps; (vi) filling factor always dominates the determination of the temperature distribution for large optical depths, and for small clump radii at smaller optical depths; (vii) at lower face-averaged optical depths, the temperature distribution is most sensitive to filling factors of 1–10 per cent, in accordance with many observations; (viii) direct shadowing by clumps can be important for distances approximately one clump radius behind a clump.     

Spatially intermittent fields in photospheric magnetoconvection

Motivated by recent high-resolution observations of the solar surface, we investigate the problem of non-linear magnetoconvection in a three-dimensional compressible layer. We present results from a set of numerical simulations which model the situation in which there is a weak imposed magnetic field. This weak-field regime is characterized by vigorous granular convection and spatially intermittent magnetic field structures. When the imposed field is very weak, magnetic flux tends to accumulate at the edges of the convective cells, where it forms compact, almost 'point-like' structures which are reminiscent of those observed in the quiet Sun. If the imposed field is slightly stronger, there is a tendency for magnetic flux to become concentrated into 'ribbon-like' structures which are comparable to those observed in solar plages. The dependence of these simulations upon the strength of the imposed magnetic field is analysed in detail, and the concept of the fractal dimension is used to make a further, more quantitative comparison between these simulations and photospheric observations.     

Deep-ocean vertical noise directionality

The structure of beam noise measured at the output of a vertical array in a range dependent ocean basin was investigated using the modified wide-angle parabolic equation (PE). Noise sources were distributed throughout the basin, and the field due to each noise source at an array located in the midbasin was calculated. The response of the array to the superposition of the noise sources was found by beamforming. An efficient and direct approach that superimposes the noise sources on the PE field as the field is marched toward the array was developed. Downslope calculations of the midbasin vertical directionality were made between 50 and 400 Hz with this technique. Use of a geoacoustic model shows that the bottom behaves as a low-pass filter     

Inverted echo sounder measurement of dynamic height through an ENSOcycle in the central equatorial Pacific

A four-year record from an inverted echo sounder deployed near Palmyra Island at 6°N in the central Pacific Ocean is compared with a simultaneous record of subsurface pressure from this island lagoon. A factor m, converting round-trip acoustic travel time to surface dynamic height relative to a deep pressure level, was estimated from the ratio of the spectra of the two records in the energetic synoptic oscillation band. Year-to-year variation in m was not statistically significant. For the overall record, m was found to be -70±8 dynamic m/s, where the error bounds represent a 90% confidence interval. This is consistent with first-baroclinic-mode excitation     

Simulations of a giant hybrid air shower detector

We have used Monte Carlo simulations to investigate the capabilities of a giant air shower observatory designed to detect showers initiated by cosmic rays with energies exceeding 10¹⁹ eV. The observatory is to consist of an array of detectors that will characterise the air shower at ground level, and optical detectors to measure the fluorescence light emitted by the shower in the atmosphere. Using these detectors together in a ‘hybrid’ configuration, we find that precise geometrical reconstruction of the shower axis is possible, leading to excellent resolution in energy. The technique is also shown to provide very good reconstruction below 10¹⁹ eV, at energies where the ground array is not fully efficient.     

Derivation of frequency-dependentX- andY-functions for a non-coherent scattering problem

The equation of transfer for the case of non-coherent scattering (Hummer, 1968; Ivanov, 1973; McCormick and Siewert, 1970) has been considered. The correspondingX- andY-functions have been derived by a combination of eigenfunction method developed by Case, and from the principle of invariance as developed by Chandrasekhar (1960).