An HLLC Riemann solver for relativistic flows – II. Magnetohydrodynamics

An approximate Riemann solver for the equations of relativistic magnetohydrodynamics (RMHD) is derived. The Harten–Lax–van Leer contact wave (HLLC) solver, originally developed by Toro, Spruce and Spears, generalizes the algorithm described in a previous paper to the case where magnetic fields are present. The solution to the Riemann problem is approximated by two constant states bounded by two fast shocks and separated by a tangential wave. The scheme is Jacobian-free, in the sense that it avoids the expensive characteristic decomposition of the RMHD equations and it improves over the HLL scheme by restoring the missing contact wave.
Multidimensional integration proceeds via the single step, corner transport upwind (CTU) method of Colella, combined with the constrained transport (CT) algorithm to preserve divergence-free magnetic fields. The resulting numerical scheme is simple to implement, efficient and suitable for a general equation of state. The robustness of the new algorithm is validated against one- and two-dimensional numerical test problems.     

Radiative Transfer in Inhomogeneous Atmospheres. III

The approach proposed in the previous parts of this series of papers is used to solve the radiative transfer problem in scattering and absorbing multicomponent atmospheres. Linear recurrence relations are obtained for both the reflectance and transmittance of these kinds of atmospheres, as well as for the emerging intensities when the atmosphere contains energy sources. Spectral line formation in a one-dimensional inhomogeneous atmosphere is examined as an illustration of the possibility of generalizing our approach to the matrix case. It is shown that, in this case as well, the question reduces to solving an initial value problem for linear differential equations. Some numerical calculations are presented.     

Optical monitoring of the quasar 4C 38.41

The results of a photometric monitoring of the quasar 4C 38.41, performed at the optical R and B bands in 2002 February–March, are presented. With a 60/90 cm Schmidt telescope at the Xinglong station of the National Astronomical Observatories of China, we observed the source exhibiting amplitude variations of up to 0.78 mag in both bands during the whole campaign. Intraday and even intranight variations are detected as well. A typical variability time-scale of about 36 d is derived from our 2-month observations at the optical bands, which is identical to that found at a radio wavelength of 92 cm, suggesting a common origin for the variations in 4C 38.41 from optical to radio bands.     

High-resolution calculations of merging neutron stars – II. Neutrino emission

The remnant resulting from the merger of two neutron stars produces neutrinos in copious amounts. In this paper we present the neutrino emission results obtained via Newtonian, high-resolution simulations of the coalescence event. These simulations use three-dimensional smoothed particle hydrodynamics together with a nuclear, temperature-dependent equation of state and a multiflavour neutrino leakage scheme. We present the details of our scheme, discuss the neutrino emission results from a neutron star coalescence and compare them with the core-collapse supernova case where neutrino emission has been studied for several decades. The average neutrino energies are similar to those in the supernova case, but contrary to the latter, the luminosities are dominated by electron-type antineutrinos that are produced in the hot, neutron-rich, thick disc of the merger remnant. The cooler parts of this disc contain substantial fractions of heavy nuclei, which, however, do not influence the overall neutrino emission results significantly. Our total neutrino luminosities from the merger event are considerably lower than those found in previous investigations. This imposes constraints on the ability of neutron star mergers to produce a gamma-ray burst via neutrino annihilation. The neutrinos are emitted preferentially along the initial binary rotation axis, an event seen 'pole-on' would appear much brighter in neutrinos than a similar event seen 'edge-on'.     

Bar dissolution in non-spherical halos

Dynamical evolution of N-body bars embedded in spherical and prolate dark matter halos is investigated. In particular, the configuration such that galactic disks are placed in the plane perpendicular to the equatorial plane of the prolate halos is considered. Such a configuration is frequently found in cosmological simulations. N-body disks embedded in a fixed external halo potential were simulated, so that the barred structure was formed via dynamical instability in initially cool disks. In the subsequent evolution, bars in prolate halos dissolved gradually with time, while the bar pattern in spherical halos remained almost unchanged until the end of simulations. The e-folding time of bars suggest that they could be destroyed in a time smaller than a Hubble time. This revised version was published online in August 2006 with corrections to the Cover Date.     

The efficiency of the spiral-in of a black hole to the Galactic Centre

We study the efficiency at which a black hole or dense star cluster spirals in to the Galactic Centre. This process takes place on a dynamical friction time-scale, which depends on the value of the Coulomb logarithm (ln Λ). We determine the accurate value of this parameter using the direct N -body method, a tree algorithm and a particle-mesh technique with up to two million plus one particles. The three different techniques are in excellent agreement. Our measurement for the Coulomb logarithm appears to be independent of the number of particles. We conclude that  ln Λ= 6.6 ± 0.6  for a massive point particle in the inner few parsec of the Galactic bulge. For an extended object, such as a dense star cluster, ln Λ is smaller, with a value of the logarithm argument Λ inversely proportional to the object size.     

物探方法在江西省茅排金矿的应用效果

茅排金矿床属中－低温硫化物金矿床，矿体的分布严格受断裂构造控制．通过应用不同的物探方法来寻找控矿断裂，从而达到间接地寻找矿体分布的目的，取得了较好的效果。