A note on the algorithm of symplectic integrators

For a Hamiltonian that can be separated into N+1(N\geq 2) integrable parts, four algorithms can be built for a symplectic integrator. This research compares these algorithms for the first and second order integrators. We found that they have similar local truncation errors represented by error Hamiltonian but rather different numerical stability. When the computation of the main part of the Hamiltonian, H ₀, is not expensive, we recommend to use S ^* type algorithm, which cuts the calculation of the H ₀ system into several small time steps as Malhotra(1991) did. As to the order of the N+1 parts in one step calculation, we found that from the large to small would get a slower error accumulation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Spine–sheath polarization structures in four active galactic nuclei jets

We present the results of multifrequency (15 + 8 + 5 GHz) polarization Very Long Baseline Array (VLBA) observations of the three BL Lacertae objects 0745+241, 1418+546 and 1652+398 together with 5-GHz VLBI Space Observatory Programme (VSOP) observations of 1418+546 and 1.6- and 5-GHz VSOP observations of the blazar 1055+018. The jets of all these sources have polarization structure transverse to the jet axis, with the polarization E vectors aligned with the jet along the jet spine and 'sheaths' of orthogonal E vectors at one or both edges of the jet. The presence of polarization aligned with the jet near the 'spine' may indicate that the jets are associated with helical B fields that propogate outward with the jet flow; the presence of orthogonal polarization near the edges of the jet may likewise be a consequence of a helical jet B field, or may be owing to an interaction with the ambient medium on parsec scales. We have tentatively detected interknot polarization in 1055+018 with E aligned with the local jet direction, consistent with the possibility that the jet of this source is associated with a helical B field.     

Magnetic susceptibility variations of the Edicaran cap carbonates in the Yangtze platform and their implications for paleoclimate

1 INTRODUCTION One of the most intriguing phenomena in the late Neoproterozoic (~750 to 543 Ma) is the globa occurrence of thin carbonates that directly overlie glacial deposits in almost every continent (Kennedy 1996; Hoffman et al., 1998; Hoffman and Schrag 2002; Brasier and Shields, 2000; James et al., 2001 Jiang et al., 2003; Nogueira et al., 2003). These “cap carbonates”, commonly several to tens of meters thick, have attracted enormous interests because o their unusually negati…     

Faint galaxies: merging model of E/SO galaxies

We study a model of mergers affecting only the progenitors of present E/SO. We adopt the standard scenarios of star formation as prescribed by Guiderdoni & Rocca-Volmerange. The merging process is parametrized under the assumptions of(1) self-similarity of the Schechter MF and(2) mass conservation. Nine models are discussed. The predictions are compared with counts ofB _J ,U ⁺,F ⁺,N ⁺ bands. E/SO mergers account for the excess of the faintest blue galaxies without causing excess in redder bands. However, as we no longer have enough mergers at brighter magnitudes, a plain E/SO merging model fits less tightly for the redshift and the colour distributions. Detection effect, a steeper slope of LF may be ways to improve. Our models predict acceptable merger frequencies atz = 0.5 although some models predict more interacting galaxies than observation atz = 0.     

A kinetic model of gas flow in a porous cometary mantle

A numerical study of gas flow through a porous cometary mantle is presented. A kinetic model based on the well-known Test Particle Monte Carlo Method for the solution of rarefied gas dynamics problems is proposed. The physical model consists of two spatial plane regions: the condensed ice phase and a porous dust mantle. The structure of the porous dust layer is described as a bundle of cylindrical inclined channels not crossing each other. A vertical temperature gradient may exist across the dust mantle. The aim is to investigate how the characteristics of molecular flow depend on the capillary length, inclination angle, and temperature gradient. Examples illustrating a significant deviation of some results from equilibrium values are shown. In particular, the gas velocity distribution at both ends of the pore is strongly non-Maxwellian if there is an important temperature contrast across the pore. The emergent gas flow rate is found to vary with the pore length/radius ratio in excellent agreement with Clausing's empirical formula. The degree of collimation of the flow is quantitatively studied as a function of the length/radius ratio, and consequences for the jet force of outgassing through a dust mantle or, indeed, a rough surface are estimated.     

Computer simulation of the particle acceleration in pulsar magnetospheres

In the spherically-symmetric case, a computer simulation of the electron acceleration inside the outflow channel of the pulsar magnetosphere is produced. The stationary motion of electrons is shown to be unstable in the case of > _c, where is a parameter describing inhomogeneity of the background charge, and _c is its critical value. The arising non-stationary motion of electrons leads to a formation of electron bunches, which move chaotically. The mean electron energy appears to be much greater at the non-stationary motion, than at the stationary one. The time-averaged parameters of the non-stationary electron flow and their dependence upon have been investigated. Distributions of the mean values of parameters (charge density, electron velocity, electric field energy density, pressure, and internal energy of the gas composed of the electron bunches) over the magnetosphere altitude have been investigated. The mean spectra of the charge density have been obtained. The results of numerical investigation of the spherically-symmetric model are used for estimation of the electron energy and of the electron flux in the case of the more realistic model. The radioemission loss is estimated, and is shown to be large enough for explaining the radiopulsar phenomenon as a thermal radioemission of the pulsar magnetosphere. In particular, such common properties of the pulsar radioemission as the high bright temperature, the sharp radioemission directivity, and the characteristic turn-over of the radioemission spectrum at the frequency of the order 10⁸ Hz are found a natural explanation in frames of this model.     

Neutrino mass,luminosity variation,and spectrum of SN1987A

In this paper, the neutrino mass has been determined from SN1987a observation in a manner that the simultaneity of neutrino emission is not regarded as the starting point, but is itself defined through the analysis by Monte-Carlo simulation. The result is that the neutrino mass lies in 3–4 eV, possiblym _v 3.6 eV. Neutrino luminosity variation and neutrino spectrum are also obtained. Comparison with theories gives further support to the mass determination, and also predicts the mass of progenitor star of SN1987a to be in the range of 12–25M .The project supported by the National Natural Science Foundation of China.     

A model of the galactic centre with magnetic monopoles

A model of the galactic centre with magnetic monopole has been presented here. The positron can be produced continuously through magnetic monopoles to induce baryon decay (Rubakov catalytic reaction) and a lot of energy can be released as well. The calculation results show that even if the galactic center contains only a few magnetic monopoles (=N _M/N_B10^–24), this massive object can not collapse into a black hole. This model can explain the observed intensities of the annihilation line and higher energy photons ofE >511 keV from the galactic centre.Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984.     

The production of light elements in unconventional cosmologies with VMOs

The abundances of light elements (deuterium, helium and lithium) are calculated in unconventional cosmologies in which these elements are produced by interaction of fast -particles ejected by massive pregalactic stars with a purely hydrogenous cosmological substrate. The number density and energetics of such stars needed for suitable production of light elements are estimated. It is shown that the models always lead to overproduction of lithium. The intensity of ultraviolet, X-ray, and -ray emission caused by explosions of pregalactic massive stars is calculated. The distortions of cosmic background radiation by Compton losses of electrons heated by explosive stars are also calculated.