A high-resolution radio study of neutral gas in the starburst galaxy NGC 520

We present subarcsec angular resolution observations of the neutral gas in the nearby starburst galaxy NGC 520. The central kpc region of NGC 520 contains an area of significantly enhanced star formation. The radio continuum structure of this region resolves into ∼10 continuum components. By comparing the flux densities of the brightest of these components at 1.4 GHz with published 15-GHz data we infer that these components detected at 1.4 and 1.6 GHz are related to the starburst and are most likely to be collections of several supernova remnants within the beam. None of these components is consistent with emission from an active galactic nuclei. Both neutral hydrogen (H  i ) and hydroxyl (OH) absorption lines are observed against the continuum emission, along with a weak OH maser feature probably related to the star formation activity in this galaxy. Strong H  i absorption  ( N _H∼ 10²² atoms cm⁻²)  traces a velocity gradient of 0.5 km s⁻¹ pc⁻¹ across the central kpc of NGC 520. The H  i absorption velocity structure is consistent with the velocity gradients observed in both the OH absorption and in CO emission observations. The neutral gas velocity structure observed within the central kpc of NGC 520 is attributed to a kpc-scale ring or disc. It is also noted that the velocity gradients observed for these neutral gas components appear to differ with the velocity gradients observed from optical ionized emission lines. This apparent disagreement is discussed and attributed to the extinction of the optical emission from the actual centre of this source hence implying that optical ionized emission lines are only detected from regions with significantly different radii to those sampled by the observations presented here.     

Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

The depolarization properties of powerful extragalactic radio sources as a function of cosmic epoch

We use the observed polarization properties of a sample of 26 powerful radio galaxies and radio-loud quasars to constrain the conditions in the Faraday screens local to the sources. We adopt the cosmological redshift, low-frequency radio luminosity and physical size of the large-scale radio structures as our 'fundamental' parameters. We find no correlation of the radio spectral index with any of the fundamental parameters. The observed rotation measure is also independent of these parameters, suggesting that most of the Faraday rotation occurs in the Galactic foreground. The difference between the rotation measures of the two lobes of an individual source, as well as the dispersion of the rotation measure, shows significant correlations with the source redshift, but not with the radio luminosity or source size. This is evidence that the small-scale structure observed in the rotation measure is caused by a Faraday screen local to the sources. The observed asymmetries between the lobes of our sources show no significant trends with each other or other source properties. Finally, we show that the commonly used model for the depolarization of synchrotron radio emission by foreground Faraday screens is inconsistent with our observations. We apply alternative models to our data and show that they require a strong increase of the dispersion of the rotation measure inside the Faraday screens with cosmological redshift. Correcting our observations with these models for redshift effects, we find a strong correlation of the depolarization measure with redshift and a significantly weaker correlation with radio luminosity. We do not find any (anti-)correlation of depolarization measure with source size. All our results are consistent with a decrease in the order of the magnetic field structure of the Faraday screen local to the sources for increasing cosmological redshift.     

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.     

Rayleigh-Taylor instability of two viscous superposed rotating and conducting fluids

Classical R-M and synthetic W-D analysis of V758 Centauri are presented. Two solutions (semi-detached and contact) were found from differential corrections approach. The semi-detached model is physically acceptable since the system is thermally decoupled. The solution for this case and the photometric data are consistent with a B9 primary and A9 secondary components having parameters close to Main-Sequence values. It is suggested that V758 Centauri is a B-type W UMa system at the brokencontact phase predicted by the Thermal Relaxation Oscillations theory.