Simulation Techniques for Cosmological Simulations

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.     

Particle Acceleration Mechanisms

In this paper we review the possible mechanisms for production of non-thermal electrons which are responsible for the observed non-thermal radiation in clusters of galaxies. Our primary focus is on non-thermal Bremsstrahlung and inverse Compton scattering, that produce hard X-ray emission. We first give a brief review of acceleration mechanisms and point out that in most astrophysical situations, and in particular for the intracluster medium, shocks, turbulence and plasma waves play a crucial role. We also outline how the effects of the turbulence can be accounted for. Using a generic model for turbulence and acceleration, we then consider two scenarios for production of non-thermal radiation. The first is motivated by the possibility that hard X-ray emission is due to non-thermal Bremsstrahlung by nonrelativistic particles and attempts to produce non-thermal tails by accelerating the electrons from the background plasma with an initial Maxwellian distribution. For acceleration rates smaller than the Coulomb energy loss rate, the effect of energising the plasma is to primarily heat the plasma with little sign of a distinct non-thermal tail. Such tails are discernible only for acceleration rates comparable or larger than the Coulomb loss rate. However, these tails are accompanied by significant heating and they are present for a short time of <10⁶ years, which is also the time that the tail will be thermalised. A longer period of acceleration at such rates will result in a runaway situation with most particles being accelerated to very high energies. These more exact treatments confirm the difficulty with this model, first pointed out by Petrosian (Astrophys. J. 557:560, 2001). Such non-thermal tails, even if possible, can only explain the hard X-ray but not the radio emission which needs GeV or higher energy electrons. For these and for production of hard X-rays by the inverse Compton model, we need the second scenario where there is injection and subsequent acceleration of relativistic electrons. It is shown that a steady state situation, for example arising from secondary electrons produced from cosmic ray proton scattering by background protons, will most likely lead to flatter than required electron spectra or it requires a short escape time of the electrons from the cluster. An episodic injection of relativistic electrons, presumably from galaxies or AGN, and/or episodic generation of turbulence and shocks by mergers can result in an electron spectrum consistent with observations but for only a short period of less than one billion years.     

Soft X-Ray and Extreme Ultraviolet Excess Emission from Clusters of Galaxies

An excess over the extrapolation to the extreme ultraviolet and soft X-ray ranges of the thermal emission from the hot intracluster medium has been detected in a number of clusters of galaxies. We briefly present each of the satellites (EUVE, ROSAT PSPC and BeppoSAX, and presently XMM-Newton, Chandra and Suzaku) and their corresponding instrumental issues, which are responsible for the fact that this soft excess remains controversial in a number of cases. We then review the evidence for this soft X-ray excess and discuss the possible mechanisms (thermal and non-thermal) which could be responsible for this emission.     

Orbital phase spectroscopy of four high mass X-ray binary pulsars to study the stellar wind of the companion

Our work focuses on a comprehensive orbital phase-dependent spectroscopy of the four High Mass X-ray Binary Pulsars (HMXBPs) 4U 1538-52, GX 301-2, OAO 1657-415 and Vela X-1. We hereby report the measurements of the variation of the absorption column density and iron-line flux along with other spectral parameters over the binary orbit for the above-mentioned HMXBPs in elliptical orbits, as observed with the Rossi X-ray Timing Explorer (RXTE) and the BeppoSAX satellites. A spherically symmetric wind profile was used as a model to compare the observed column density variations. Out of the four pulsars, only in 4U 1538-52, we find the model having a reasonable corroboration with the observations, whereas in the remaining three the stellar wind seems to be clumpy and a smooth symmetric stellar wind model appears to be quite inadequate in explaining the data. Moreover, in GX 301-2, neither the presence of a disk nor a gas stream from the companion was validated. Furthermore, the spectral results obtained in the case of OAO 1657-415 and Vela X-1 were more or less similar to that of GX 301-2.     

Coronal fast wave trains of the decimetric type IV radio event observed during the decay phase of the June 6, 2000 flare

The 22 min long decimetric type IV radio event observed during the decay phase of the June 6, 2000 flare simultaneously by the Brazilian Solar Spectroscope (BSS) and the Ond?ejov radiospectrograph in frequency range 1200–4500 MHz has been analyzed. We have found that the characteristic periods of about 60 s belong to the long-period spectral component of the fast wave trains with a tadpole pattern in their wavelet power spectra. We have detected these trains in the whole frequency range 1200–4500 MHz. The behavior of individual wave trains at lower frequencies is different from that at higher frequencies. These individual wave trains have some common as well as different properties. In this paper, we focus on two examples of wave trains in a loop segment and the main statistical parameters in their wavelet power and global spectra are studied and discussed.     

星载单频GPS数据的电离层延迟改正方法分析

利用C/A码单点定位对LEO(Low Earth Orbit)卫星上的电离层延迟改正方法——"电离层比例因子法"进行了分析研究.计算的CHAMP卫星的轨道结果表明:采用电子密度峰值高度(hmF2,F2 region maximum electron density height)平均值和瞬时值计算的电离层比例因子α变化范围分别为0.3~0.4和0.2~0.65之间,两者最大差异可达0.3,相比较而言,hmF2瞬时值的结果更加合理,并且相应的大地高H方向的系统偏差要降低0.05~0.3m左右;与双频无电离层组合的普通单点定位结果相比表明该方法能较好地消除电离层一阶项所引入的H方向上的系统偏差;该方法适用的LEO卫星轨道高度范围大致在200~ 600km之间,当轨道高度超过700km时,该方法并不适用.     

Response of Hainan GPS ionospheric scintillations to the different strong magnetic storm conditions

Using the GPS ionospheric scintillation data at Hainan station (19.5°N, 109.1°E) in the eastern Asia equatorial regions and relevant ionospheric and geomagnetic data from July 2003 to June 2005, we investigate the response of L-band ionospheric scintillation activity over this region to different strong magnetic storm conditions (Dst < −100 nT) during the descending phase of the solar cycle. These strong storms and corresponding scintillations mainly took place in winter and summer seasons. When the main phase developed rapidly and reached the maximum near 20–21 LT (LT = UT + 8) after sunset, scintillations might occur in the following recovery phase. When the main phase maximum occurred shortly after midnight near 01–02 LT, following the strong scintillations in the pre-midnight main phase, scintillations might also occur in the post-midnight recovery phase. When the main phase maximum took place after 03 LT to the early morning hours no any scintillation could be observed in the latter of the night. Moreover, when the main phase maximum occurred during the daytime hours, scintillations could also hardly be observed in the following nighttime recovery phase, which might last until the end of recovery phase. Occasionally, scintillations also took place in the initial phase of the storm. During those scintillations associated with the nighttime magnetic storms, the height of F layer base (h’F) was evidently increased. However, the increase of F layer base height does not always cause the occurrence of scintillations, which indicates the complex interaction of various disturbance processes in ionosphere and thermosphere systems during the storms.     

Spectroscopic analysis of the solar flare event on 2002 August 3 with the use of RHESSI and RESIK data

We use simultaneous observations from RESIK and RHESSI instruments to compare plasma properties of a major solar flare in its rise and gradual phase. This event occurred on 2002 August 3 (peak time at 19:06 UT). The flare had a very good coverage with RESIK data and well-resolved soft and hard X-ray sources were seen in RHESSI images. Spectra of X-ray radiation from RHESSI images are studied and compared with RESIK measurements in different flare phases. Result shows large differences in flare morphology and spectra between flare rise and gradual phase.