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.     

On the relation between orbital structure and observed bar morphology

We investigate the morphological relation between the orbits of the central family of periodic orbits ( x ₁ family) and the bar itself using models of test particles moving in a barred potential. We show that different bar morphologies may have as a backbone the same set of x ₁ periodic orbits. We point out that by populating initially axisymmetric stellar discs exponentially with test particles in circular, or almost circular motion, we may end up with a response bar which reveals a shape different in crucial details from that of the individual stable x ₁ orbits. For example, a bar model in which the x ₁ orbits are pure ellipses may have a much more complicated response morphology. This depends on the particular invariant curves around x ₁, which are populated in each model.     

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.     

Shock Propagation Through Multiphase Media*

This paper presents two and three dimensional simulations of the interaction of shocks with media with large numbers of dense inclusions. An approximate model of the interaction of a starburst wind with the surrounding galactic ISM illustrates issues which must be addressed in global models of ISM dynamics. As a step towards developing the sub-grid model of multiphase turbulence, we define and study a form of ‘multiphase Riemann problem’. This allows us to develop macroscopic characteristics of the flows which may be compared to such subgrid models.