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.     

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.     

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.     

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.     

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.     

Modulation of galactic cosmic ray anisotropy in heliomagnetosphere: Influence of cosmic ray scattering on sidereal daily variation

In previous papers, the present authors have shown that the galactic anisotropy is modulated due to cosmic ray orbital deflection in the heliomagnetosphere, and that the sidereal time daily variations of galactic origin can be expressed using the basic vectors, which have been obtained by calculating trajectories of cosmic rays in a model magnetosphere having Parker's Archimedian spiral structure with a flat or a wavy neutral sheet. In the present paper, the magnetic irregularities superposed on the Parker's spiral field have been taken into account, which cause the scattering of cosmic rays and disturb their orbits. We examined the fluctuations of asymptotic directions calculating their orbits by the Monte-Carlo simulation, based on the theory of the multiple scattering process. It is shown that the dispersion of the projected deviation angle is determined mainly by the scattering mean free path and by the structure of the order magnetic field, e.g. the polarity state of the heliomagnetosphere and the extent of the neutral sheet. We investigated also the influence of the fluctuations of asymptotic directions on the sidereal daily variation. It is found that, under some conditions, the scattering causes only the attenuation of the amplitude of the basic vector, and does not change its phase. The attenuation is negligibly small at high rigidities larger than ~ 1000 GV, but becomes more serious with decreasing rigidity. The rigidity dependence curve of the attenuation rate was calculated for various cases. A simple and approximate method is also presented for the derivation of those curves for any value of the magnitude of the mean free path and for various model magnetospheres. It is noted, however, that the lower limiting rigidity below which the present method is not applicable is relatively high in the Positive polarity state.     

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.     

The Citation Impact of Digital Preprint Archives for Solar Physics Papers

Papers that are posted to a digital preprint archive are typically cited twice as often as papers that are not posted. This has been demonstrated for papers published in a wide variety of journals, and in many different subfields of astronomy. Most astronomers now use the arXiv.org server (astro-ph) to distribute preprints, but the solar physics community has an independent archive hosted at Montana State University. For several samples of solar physics papers published in 2003, I quantify the boost in citation rates for preprints posted to each of these servers. I show that papers on the MSU archive typically have citation rates 1.7 times higher than the average of similar papers that are not posted as preprints, while those posted to astro-ph get 2.6 times the average. A comparable boost is found for papers published in conference proceedings, suggesting that the higher citation rates are not the result of self-selection of above-average papers. Editors’ Note: This paper lies outside the normal purview of Solar Physics papers, however the editors feel that the content is of sufficient importance for all Solar Physics authors and readers to merit its publication.