用粒子模拟方法研究离子声波

<正>向传播朗缪尔波被离子声波散射是太阳射电Ⅲ型暴基波和谐波激发的重要过程.使用粒子模拟方法对电子束流激发朗缪尔波的过程进行了模拟,同时对产生的反向朗缪尔波、朗缪尔波2次谐波和朗缪尔波通过非线性过程产生的离子声波的性质进行了分析研究.为了更好地研究离子声波,模拟时单独计算了由离子扰动引起的电场.模拟计算了不同初始参数下产生的离子声波强度,发现离子的温度和质量对离子声波的产生有重要作用,验证了反向朗缪尔波与离子声波的相关性.同时在模拟中验证了朗缪尔波的衰变过程,确认了离子声波对反向朗缪尔波的放大作用.     

Magnetohydrodynamics in superconducting-superfluid neutron stars

Magnetohydrodynamic (MHD) equations are presented for the mixture of superfluid neutrons, superconducting protons and normal electrons believed to exist in the outer cores of neutron stars. The dissipative effects of electron viscosity and mutual friction resulting from electron-vortex scattering are also included. It is shown that Alfvén waves are replaced by cyclotron-vortex waves that have not been previously derived from MHD theory. The cyclotron-vortex waves are analogous to Alfvén waves with the tension arising from the magnetic energy density replaced by the vortex energy density. The equations are then put into a simplified form useful for studying the effect of the interior magnetic field on the dynamics. Of particular interest is the crust–core coupling time, which can be inferred from pulsar glitch observations. The hypothesis that cyclotron-vortex waves play a significant role in the core spin-up during a glitch is used to place limits on the interior magnetic field. The results are compared with those of other studies.     

Solar coronal heating by magnetosonic waves

Solar coronal heating by magnetohydrodynamic (MHD) waves is investigated. ultraviolet (UV) and X-ray emission lines of the corona show non-thermal broadenings. The wave rms velocities inferred from these observations are of the order of 25–60 km s⁻¹ . Assuming that these values are not negligible, we solved MHD equations in a quasi-linear approximation, by retaining the lowest order non-linear term in rms velocity. Plasma density distribution in the solar corona is assumed to be inhomogeneous. This plasma is also assumed to be permeated by dipole-like magnetic loops. Wave propagation is considered along the magnetic field lines. As dissipative processes, only the viscosity and parallel (to the local magnetic field lines) heat conduction are assumed to be important. Two wave modes emerged from the solution of the dispersion relation. The fast mode magneto-acoustic wave, if originated from the coronal base can propagate upwards into the corona and dissipate its mechanical energy as heat. The damping length-scale of the fast mode is of the order of 500 km. The wave energy flux associated with these waves turned out to be of the order of 2.5×10⁵ ergs cm⁻² s⁻¹ which is high enough to replace the energy lost by thermal conduction to the transition region and by optically thin coronal emission. The fast magneto-acoustic waves prove to be a likely candidate to heat the solar corona. The slow mode is absent, in other words cannot propagate in the solar corona.     

Energy, angular momentum and wave action associated with density waves in a rotating magnetized gas disc

Both fast and slow magnetohydrodynamic (MHD) density waves propagating in a thin rotating magnetized gas disc are investigated. In the tight-winding or WKBJ regime, the radial variation of MHD density-wave amplitude during wave propagation is governed by the conservation of wave action surface density which travels at a relevant radial group speed C _g. The wave energy surface density and the wave angular momentum surface density are related to by = and = m respectively, where is the angular frequency in an inertial frame of reference and the integer m , proportional to the azimuthal wavenumber, corresponds to the number of spiral arms. Consequently, both wave energy and angular momentum are conserved for spiral MHD density waves. For both fast and slow MHD density waves, net wave energy and angular momentum are carried outward or inward for trailing or leading spirals, respectively. The wave angular momentum flux contains separate contributions from gravity torque, advective transport and magnetic torque. While the gravity torque plays an important role, the latter two can be of comparable magnitudes to the former. Similar to the role of gravity torque, the part of MHD wave angular momentum flux by magnetic torque (in the case of either fast or slow MHD density waves) propagates outward or inward for trailing or leading spirals, respectively. From the perspective of global energetics in a magnetized gas sheet in rotation, trailing spiral structures of MHD density waves are preferred over leading ones. With proper qualifications, the generation and maintenance as well as transport properties of MHD density waves in magnetized spiral galaxies are discussed.     

The impact of correlated noise on SuperWASP detection rates for transiting extrasolar planets

The evolution of the Alfvén turbulence due to three-wave interactions is discussed using kinetic theory for a collisionless, thermal plasma. There are three low-frequency modes, analogous to the three modes of compressible magnetohydrodynamics (MHD). When only Alfvén waves are considered, the known anisotropy of turbulence in incompressible MHD theory is reproduced. Inclusion of a fast mode wave leads to the separation of turbulence into two regimes: small wave numbers where three-wave processes involving a fast mode are dominant, and large wave numbers where the three Alfvén wave process is dominant. Possible application of the anisotropic Alfvén turbulence to the interstellar medium and dissipation of magnetic energy in magnetars are discussed.     

Nature of microstructure in pulsar radio emission

We present a model for microstructure in pulsar radio emission. We propose that micropulses result from alteration of the radio wave generation region by nearly transverse drift waves propagating across the pulsar magnetic field and encircling the bundle of the open magnetic field lines. It is demonstrated that such waves can modify the curvature of the field lines significantly. This, in turn, affects strongly fulfilment of the resonance conditions necessary for the excitation of radio waves. The time-scale of micropulses is therefore determined by the wavelength of the drift waves. The main features of the microstructure are naturally explained within the framework of this model.     

The excitation, propagation and dissipation of waves in accretion discs: the non-linear axisymmetric case

We analyse the non-linear propagation and dissipation of axisymmetric waves in accretion discs using the ZEUS-2D hydrodynamics code. The waves are numerically resolved in the vertical and radial directions. Both vertically isothermal and thermally stratified accretion discs are considered. The waves are generated by means of resonant forcing, and several forms of forcing are considered. Compressional motions are taken to be locally adiabatic  ( γ =5/3)  . Prior to non-linear dissipation, the numerical results are in excellent agreement with the linear theory of wave channelling in predicting the types of modes that are excited, the energy flux by carried by each mode, and the vertical wave energy distribution as a function of radius. In all cases, waves are excited that propagate on both sides of the resonance (inwards and outwards). For vertically isothermal discs, non-linear dissipation occurs primarily through shocks that result from the classical steepening of acoustic waves. For discs that are substantially thermally stratified, wave channelling is the primary mechanism for shock generation. Wave channelling boosts the Mach number of the wave by vertically confining the wave to a small cool region at the base of the disc atmosphere. In general, outwardly propagating waves with Mach numbers near resonance  ℳ_r≳0.01  undergo shocks within a distance of order the resonance radius.     

Sound waves in the intracluster medium

We present an analysis of the behaviour of a perturbed radio cocoon. Comparisons with observations of sound waves detected in the Perseus and Virgo clusters suggest the separations of observed ripples correspond to the natural oscillation frequency of the cocoon. An energy injection rate consistent with active galactic nucleus power is required to offset the strong acoustic damping of cocoon oscillations, suggesting the sources are in equilibrium with the intracluster medium (ICM), and the oscillations are effectively undamped. Viscous dissipation of sound waves provides ICM heating that can quench cooling flows on time-scales greatly exceeding the oscillation time-scale. Thermal conductivity is likely to be heavily suppressed.     

AGN feedback through sound wave dissipation

We investigate the possibility of explaining the observed ripples in the X-ray gas in the Perseus and Virgo clusters through natural oscillations of a perturbed radio cocoon. Such a perturbation would result from an expanding overpressured cocoon of radio plasma overshooting its pressure equilibrium point with the cluster gas. The oscillations are heavily acoustically damped, and energy injection rates required to sustain them are consistent with observed AGN powers. Viscous dissipation of sound waves generated by these oscillations heats the cluster gas. By comparing our model with observations in Perseus and Virgo, we reproduce the observed ripple separations and amplitudes. Spitzer viscosity is largely sufficient in explaining the gas density profile, suggesting that thermal conductivity is likely to be heavily suppressed. In the central regions, viscous heating can suppress cooling flows on timescales exceeding the radio source lifetime.     

The stochastic background of gravitational waves from neutron star formation at cosmological distances

We investigate the stochastic gravitational wave background that results from neutron star birth throughout the Universe. The neutron star birth rate, as a function of redshift, is calculated using an observation-based model for the evolving star formation rate, together with an estimate of the rate of core-collapse supernovae in the nearby Universe and an estimate of the neutron star/black hole branching ratio. Using three sample waveforms, based on numerical models of stellar core collapse by Zwerger & Müller, the spectral flux density, spectral strain, spectral energy density and duty cycle of the background have been computed. Our results show, contrary to recent claims, that the spectrum of the stochastic background clearly reflects the different physics in the core-collapse models. For a star formation model that is corrected for dust extinction, the neutron star formation rate throughout the Universe is high enough to result in a nearly continuous background of gravitational waves, with spectral features that can be related to emission mechanisms.     

A new mechanism for pulsar gamma-ray emission

Gamma-ray emission in pulsar magnetospheres is attributed to synchrotron radiation, which tends to decrease the pitch angle of the particle, being balanced by plasma processes tending to increase the pitch angle. The plasma processes are non-resonant instabilities that drive non-resonant quasilinear diffusion (NQD), thereby pumping energy from waves and the parallel motion of the particle into the perpendicular motion of the particle. It is shown that NQD can maintain the pitch angles for particles near the light-cylinder such that they radiate synchrotron radiation at MeV energies. Compared to conventional emission mechanisms (such as polar cap or outer gap models), the resulting spectrum has a relatively low upper cut-off from about a few to 100 MeV. Possible observational consequences of this mechanism are discussed.     

Radio signature of cosmological structure formation shocks

In the course of the formation of cosmological structures, large shock waves are generated in the intracluster medium (ICM). In analogy to processes in supernova remnants, these shock waves may generate a significant population of relativistic electrons which, in turn, produce observable synchrotron emission. The extended radio relics found at the periphery of several clusters and possibly also a fraction of radio halo emission may have this origin. Here, we derive an analytic expression for (i) the total radio power in the downstream region of a cosmological shock wave, and (ii) the width of the radio-emitting region. These expressions predict a spectral slope close to −1 for strong shocks. Moderate shocks, such as those produced in mergers between clusters of galaxies, lead to a somewhat steeper spectrum. Moreover, we predict an upper limit for the radio power of cosmological shocks. Comparing our results to the radio relics in Abell 115, 2256 and 3667, we conclude that the magnetic field in these relics is typically at a level of 0.1 μG. Magnetic fields in the ICM are presumably generated by the shocks themselves; this allows us to calculate the radio emission as a function of the cluster temperature. The resulting emissions agree very well with the radio power–temperature relation found for cluster haloes. Finally, we show that cosmic accretion shocks generate less radio emission than merger shock waves. The latter may, however, be detected with upcoming radio telescopes.     

Possible Contribution of Mature γ-ray Pulsars to Cosmic-ray Positrons

We restudy the possible contribution of mature gamma-ray pulsars to cosmic ray positrons based on the new version of outer gap model. In this model, the inclination angle and average properties of the outer gap are taken into account, and more mature pulsars can have the outer gap and emit high energy photons. Half of the primary particles in the outer gaps will flow back toward the star surface and emit synchrotron photons, which can produce electron/positron pairs by the cascade of pair production. Some of these pairs will escape from the light cylinder and be accelerated to relativistic energies in the pulsar wind driven by low-frequency electromagnetic waves. Using a Monte Carlo method, we obtain a sample of mature gamma-ray pulsars and then calculate the production of the positrons from these pulsars. The observed excess of cosmic positrons can be well explained by this model.     

Investigating ultra-long gravitational waves with measurements of pulsar rotational parameters

A method is suggested with which to explore the gravitational wave background (GWB) in the frequency range 10⁻¹²–10⁻⁸ Hz. This method is based on the precise measurements of pulsar rotational parameters: the influence of gravitational waves (GWs) in this frequency range will affect these parameters and therefore some conclusions about the energy density of the GWB can be made using analysis of the derivatives of pulsar rotational frequency. The calculated values of the second derivative from a number of pulsars limit the density of the GWB, Ω_gw, as follows:  Ω_gw < 2 × 10⁻⁶  . Also, the time series of the frequency ν of different pulsars in a pulsar array can be cross-correlated pairwise in the same manner as in anomalous residuals analysis, thus providing the possibility of GWB detection in the ultra-low-frequency range.     

On viscosity, conduction and sound waves in the intracluster medium

Recent X-ray and optical observations of the Perseus cluster indicate that a combination of weak shocks at small radii  (≳20  kpc)  and viscous and conductive dissipation of sound waves at larger radii is responsible for heating the intracluster medium and can balance radiative cooling of cluster cores. We discuss this mechanism more generally and show how the specific heating and cooling rates vary with temperature and radius. It appears that this heating mechanism is most effective above  10⁷  K  , which allows for radiative cooling to proceed within normal galaxy formation but stifles the growth of very massive galaxies. The scaling of the wavelength of sound waves with cluster temperature and feedback in the system are investigated.     

The entropy and energy of intergalactic gas in galaxy clusters

Studies of the X-ray surface brightness profiles of clusters, coupled with theoretical considerations, suggest that the breaking of self-similarity in the hot gas results from an 'entropy floor', established by some heating process, which affects the structure of the intracluster gas strongly in lower-mass systems. By fitting analytical models for the radial variation in gas density and temperature to X-ray spectral images from the ROSAT PSPC and ASCA GIS, we have derived gas entropy profiles for 20 galaxy clusters and groups. We show that, when these profiles are scaled such that they should lie on top of one another in the case of self-similarity, the lowest-mass systems have higher-scaled entropy profiles than more massive systems. This appears to be due to a baseline entropy of depending on the extent to which shocks have been suppressed in low-mass systems. The extra entropy may be present in all systems, but is detectable only in poor clusters, where it is significant compared with the entropy generated by gravitational collapse. This excess entropy appears to be distributed uniformly with radius outside the central cooling regions.
We determine the energy associated with this entropy floor, by studying the net reduction in binding energy of the gas in low-mass systems, and find that it corresponds to a pre-heating temperature of 0.3 keV. Since the relationship between entropy and energy injection depends upon gas density, we are able to combine the excesses of 70140 keV cm² and 0.3 keV to derive the typical electron density of the gas into which the energy was injected. The resulting value of implies that the heating must have happened prior to cluster collapse but after a redshift z 710. The energy requirement is well matched to the energy from supernova explosions responsible for the metals which now pollute the intracluster gas.     

Sound waves excitation by jet-inflated bubbles in clusters of galaxies

We show that repeated sound waves in the intracluster medium (ICM) can be excited by a single inflation episode of an opposite bubble pair. To reproduce this behaviour in numerical simulations, the bubbles should be inflated by jets, rather than being injected artificially as already full-blown bubbles. The multiple sound waves are excited by the motion of the bubble–ICM boundary that is caused by vortices inside the inflated bubbles and the backflow ('cocoon') of the ICM around the bubble. These sound waves form a structure that can account for the ripples observed in the Perseus cooling flow cluster. We inflate the bubbles using slow massive jets either with a very wide opening angle or that are narrow and precessing. The wide jets (or collimated fast winds) are slow in the sense that they are highly subrelativistic,   v _j∼ 0.01 c – 0.1 c   , and they are massive in the sense that the pair of bubbles carries back to the ICM a large fraction of the cooling mass, i.e.  ∼1–50 M_⊙ yr⁻¹  . We use a two-dimensional axisymmetric (referred to as 2.5D) hydrodynamical numerical code ( vh-1 ).     

Using X‐ray observations to identify the particle acceleration mechanisms in VHE SNRs and “dark” VHE sources

Very high energy (VHE) γ‐ray observations have proven to be very successful in localizing Galactic acceleration sites of VHE particles. Observations of shell‐type supernova remnants have confirmed that particles are accelerated to VHE energies in supernova blast waves; the interpretation of the γ‐ray data in terms of hadronic or leptonic particle components in these objects relies nevertheless strongly on input from X‐ray observations. The largest identified Galactic VHE source class consists of pulsar wind nebulae, as detected in X‐rays. Many of the remaining VHE sources remain however unidentified until now. With X‐ray observations of these enigmatic “dark” objects one hopes to solve the following questions: What is the astrophysical nature of these sources? Are they predominantly electron or hadron accelerators? And what is their contribution to the overall cosmic ray energy budget? The paper aims to provide an overview over the identification status of the Galactic VHE source population. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

tempo2: a new pulsar timing package – III. Gravitational wave simulation

Analysis of pulsar timing data sets may provide the first direct detection of gravitational waves. This paper, the third in a series describing the mathematical framework implemented into the tempo2 pulsar timing package, reports on using tempo2 to simulate the timing residuals induced by gravitational waves. The tempo2 simulations can be used to provide upper bounds on the amplitude of an isotropic, stochastic, gravitational wave background in our Galaxy and to determine the sensitivity of a given pulsar timing experiment to individual, supermassive, binary black hole systems.     

On the formation of cluster radio relics

In several merging clusters of galaxies so-called cluster radio relics have been observed. These are extended radio sources which do not seem to be associated with any radio galaxies. Two competing physical mechanisms to accelerate the radio-emitting electrons have been proposed: (i) diffusive shock acceleration and (ii) adiabatic compression of fossil radio plasma by merger shock waves. Here the second scenario is investigated. We present detailed three-dimensional magneto-hydrodynamical simulations of the passage of a radio plasma cocoon filled with turbulent magnetic fields through a shock wave. Taking into account synchrotron, inverse Compton and adiabatic energy losses and gains, we evolved the relativistic electron population to produce synthetic polarization radio maps. On contact with the shock wave the radio cocoons are first compressed and finally torn into filamentary structures, as is observed in several cluster radio relics. In the synthetic radio maps the electric polarization vectors are mostly perpendicular to the filamentary radio structures. If the magnetic field inside the cocoon is not too strong, the initially spherical radio cocoon is transformed into a torus after the passage of the shock wave. Very recent, high-resolution radio maps of cluster radio relics seem to exhibit such toroidal geometries in some cases. This supports the hypothesis that cluster radio relics are fossil radio cocoons that have been revived by a shock wave. For a late-stage relic the ratio of its global diameter to the filament diameter should correlate with the shock strength. Finally, we argue that the total radio polarization of a radio relic should be well correlated with the three-dimensional orientation of the shock wave that produced the relic.