An iterative scheme for use in MT 2D forward modelling — to speed up production of the Jacobian and to improve accuracy of the inverse matrix

An important stage in two-dimensional magnetotelluric modelling is the calculation of the Earth's response functions for an assumed conductivity model and the calculation of the associated Jacobian relating those response functions to the model parameters. The efficiency of the calculation of the Jacobian will affect the efficiency of the inversion modelling. Rodi (1976) produced all the Jacobian elements by inverting a single matrix and using an approximate first-order algorithm. Since only one inverse matrix required calculation the procedure speeded up the inversion. An iterative scheme to improve the approximation to the Jacobian information is presented in this paper. While this scheme takes a little longer than Rodi's algorithm, it enables a more accurate determination of the Jacobian information. It is found that the Jacobian elements can be produced in 10% of the time required to calculate an inverse matrix or to calculate a 2D starting model. A modification of the algorithm can further be used to improve the accuracy of the original inverse matrix calculated in a 2D finite difference program and hence the solution this program produces. The convergence of the iteration scheme is found to be related both to the originally calculated inverse matrix and to the change in the newly formed matrix arising from perturbation of the model parameter. A ridge regression inverse algorithm is used in conjunction with the iterative scheme for forward modelling described in this paper to produce a 2D conductivity section from field data.     

Photometric Inference of the Corotation Radii of Ten Barred Galaxies

We have used photometric images of ten barred galaxies in the B and I bands to infer the geometrical and dynamical parameters of their bars: Their length, (R _bar), their strengths, (S _b), and their corotation radii, (R _CR).These parameters have been obtained studying azimuthal profiles from the B and I band images of the galaxies. We find that R _CR is in all cases slightly larger than R _bar, with mean values close to 1.2 R _bar. We have also found a dependence of the ratio R _CR/R _bar on S _bar. This revised version was published online in July 2006 with corrections to the Cover Date.     

Rosetta mission mutual impedance probe modelling: the short and long Debye length plasma cases

The Rosetta spacecraft (S/C), which is planned to meet comet 46P/Wirtanen in 2011, will carry a set of five wave and plasma instruments (i.e. the Rosetta Plasma Consortium). This is to measure the cometary plasma properties from the minimum value of activity of the comet to its maximum value at perihelion. The mutual impedance probe, MIP, is one of those (Trotignon et al., 1999) five. This revised version was published online in July 2006 with corrections to the Cover Date.     

A Negative Branch of Polarization for Comets and Atmosphereless Celestial Bodies and the Light Scattering by Aggregate Particles

Optical observations of comets and atmosphereless celestial bodies show that a change of sign of the linear polarization of scattered light from negative to positive at phase angles less than 20° is typical of the cometary coma, as well as of the regolith of Mercury, the Moon, planetary satellites, and asteroids. To explain a negative branch of polarization, this research suggests a unified approach to the treatment of cometary-dust particles and regolith grains as aggregate forms. A composite structure of aggregate particles resulting in the interaction of composing structural elements (monomers) in the light-scattering process is responsible for the negative polarization at small phase angles, if the monomer sizes are comparable to the wavelength. The characteristics of single scattering of light calculated for aggregates of this kind turned out to be close to the properties observed for cometary dust. Unlike the cometary coma, the regolith is an optically semi-infinite medium, where the interaction between particles is significant. To find the reflectance characteristics of regolith, the radiative-transfer equation should be solved for a regolith layer. In this case, the interaction between scatterers can be modeled to a certain extent by representing the regolith grains as aggregate structures consisting of several or many elements. Although real regolith grains are much larger than the particles considered here, laboratory measurements have shown that it is precisely the surface irregularities comparable to the wavelength that cause a negative branch of polarization. The main observed features of the phase and spectral dependence of the linear polarization of light scattered from comets and atmosphereless celestial bodies, which are due to the difference of the elementary scatterers in composition, size, and structure, can be successfully explained using the aggregate model of particles.     

Regolith Surface Reflectance: A New Attempt to Model

The reflectance coefficient of the regolith layer of celestial bodies has been studied in relation to the physical properties of regolith particles (size, refractive index, and packing density) on the basis of an accurate numerical radiative-transfer algorithm for a semi-infinite flat layer. Using the geometric-optics approximation, we have found that a shape mixture of randomly oriented spheroids can successfully model the single-scattering phase function of independent soil grains. In order to take into account the effect of packing density in a regolith layer, the concept of the so-called static structure factor was used. The main effect of increasing packing density is to suppress the forward-scattering peak of the phase function and to increase the albedo of the reflecting surface. We also investigated the influence of fine dust on the reflected light. An addition of small particles not only increases the surface albedo, but also changes the brightness profile and enhances the backscattering. Although the problem of unique solution, which is inherent in the retrieval of the properties of a medium from the measurements of the intensity of light scattered by this media, cannot be removed in the proposed model, the procedure used here, in contrast to widely used approximations, allows us to fit observational data with a set of real characteristics of the regolith. Semiempirical approaches are able to fit the measurements well with a small number of free parameters, but they do not explicitly contain crucial physical characteristics of the regolith such as grain sizes or the refractive index. We compared the numerical solution of the radiative-transfer equation with the Hapke approximation, which is most often used by investigators. The errors introduced by the Hapke model are small only for near-isotropic scattering by isolated particles. However, independent regolith grains are known to scatter light mainly in the forward direction.     

Multiple stellar populations in the Sextans dwarf spheroidal galaxy?

We present wide-field     multiband ( BVI ) CCD photometry (down to     of the very low surface brightness dwarf spheroidal galaxy Sextans. In the derived colour–magnitude diagrams we find evidence suggesting the presence of multiple stellar populations in this dwarf spheroidal. In particular, we discover (i) a blue horizontal branch tail that appears to lie on a brighter sequence with respect to the prominent red horizontal branch and the RR Lyrae stars, very similar to what was found by Majewski et al. for the Sculptor dwarf spheroidal, (ii) hints of a bimodal distribution in colour of the red giant branch stars, (iii) a double red giant branch bump. All of these features suggest that (at least) two components are present in the old stellar population of this galaxy: the main one with     and a minor component around     . The similarity to the Sculptor case may indicate that multiple star formation episodes are also common in the most nearby dwarf spheroidals that ceased their star formation activity at very early epochs.     

BRIGHTNESS VARIATIONS OF Mg Ib BRIGHT FEATURES

We study the temporal intensity variations of Mgib bright features and investigate the corresponding Hα velocity pattern. The network bright features are well visible in the continuum, in images averaged over the duration of the observations (130 min). We detected `flashing' bright features, which appear and disappear within two to five minutes, while the rest of the bright features undergo small variations of either their shape and/or their intensity. A power spectrum analysis reveals a 10-min oscillation for approximately half of the stable bright features. The 5-min oscillations are detected mainly at the network boundaries, where stable bright features are located, while 3-min oscillations coincide with few bright features, but are also quite intense inside the network cells. The majority of bright features are associated with Hα downflows. The downflow is very intense at the location of `flashing' bright features.     

Observations of GRS 1915+105 from the USA Experiment on ARGOS

We present first results from a monitoring campaign of GRS 1915+105 undertaken with the USA X-ray timing experiment on the ARGOS satellite. A variety of behaviour has been observed, ranging from low, steady X-ray emission to rapid quasi-periodic flaring on timescales of approximately 10–120 seconds.     

A unified theory of the formation of galactic bars

We give arguments for a basically unified formation mechanism of slow (Lynden-Bell) and fast (common) galactic bars. This mechanism is based on an instability that is akin to the well-known instability of radial orbits and is produced by the mutual attraction and alignment of precessing stellar orbits (so far, only the formation of slow bars has been explained in this way). We present a general theory of the low-frequency modes in a disk that consists of orbits precessing at different angular velocities. The problem of determining these modes is reduced to integral equations of moderately complex structure. The characteristic pattern angular velocities Ω_p of the low-frequency modes are of the order of the mean orbital precession angular velocity \(\bar \Omega _{pr}\). Bar modes are also among the low-frequency modes; while \(\Omega _p \approx \bar \Omega _{pr}\) for slow bars, Ω_p for fast bars can appreciably exceed even the maximum orbital precession angular velocity in the disk Ω _pr ^max (however, it remains of the order of these precession angular velocities). The possibility of such an excess of Ω_p over Ω _pr ^max is associated with the effect of “repelling” orbits. The latter tend to move in a direction opposite to the direction in which they are pushed. We analyze the pattern of orbital precession in potentials typical of galactic disks. We note that the maximum radius of an “attracting” circular orbit r_c can serve as a reasonable estimate of the bar length l_b. Such an estimate is in good agreement with the available results of N-body simulations.     

Chandra observations of neutron stars: an overview

We present an overview of Chandra X-ray Observatory observations of neutron stars. The outstanding spatial and spectral resolution of this great observatory have allowed for observations of unprecedented clarity and accuracy. Many of these observations have provided new insights into neutron star physics. We present an admittedly biased and overly brief review of these observations, highlighting some new discoveries made possible by the Observatory’s unique capabilities. This includes our analysis of recent multiwavelength observations of the putative pulsar and its pulsar-wind nebula in the IC443 SNR.