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.     

Day-to-day variability of equatorial anomaly in GPS-TEC during low solar activity period

The ionospheric total electron content (TEC) in the northern hemispheric equatorial ionization anomaly (EIA) crest region is investigated by using dual-frequency signals of the Global Positioning System (GPS) acquired from Rajkot (Geog. Lat. 22.29°N, Geog. Long. 70.74°E; Geom. Lat. 14.21°N, Geom. Long. 144.90°E), India. The day-to-day variability of EIA characteristics is examined during low solar activity period (F10.7∼83 sfu). It is found that the daily maximum TEC at EIA crest exhibits a day-to-day and strong semi-annual variability. The seasonal anomaly and equinoctial asymmetry in TEC at EIA is found non-existent and weaker, respectively. We found a moderate and positive correlation of daily magnitude of crest, Ic with daily F10.7 and EUV fluxes with a correlation coefficient of 0.43 and 0.33, respectively indicating an existence of a short-term relation between TEC at EIA and the solar radiation even during low solar activity period. The correlation of daily Ic with Dst index is also moderate (r = −0.35), whereas no correlation is found with the daily Kp index (r = 0.14) respectively. We found that the magnitude of EIA crest is moderately correlated with solar flux in all seasons except winter where it is weakly related (0.27). The magnitude of EIA crest is also found highly related with EEJ strength in spring (r = 0.69) and summer (r = 0.65) than autumn (0.5) and winter (r = 0.47), though EEJ is stronger in autumn than spring.     

GNSS single frequency ionospheric range delay corrections: NeQuick data ingestion technique

A study of the performance of the NeQuick model and the Klobuchar model for GNSS single frequency range delay correction on a global scale was done using data for moderate solar activity. In this study NeQuick was used in the way intended for Galileo. This study is to assess the performance of the two models at each ionospheric geographic region during moderate solar activity as previously published studies were concentrated only on high solar activity. The results obtained showed that NeQuick outperformed Klobuchar for the whole year at the three geographical regions of the ionosphere. In terms of monthly root mean square of mismodeling, NeQuick outperformed Klobuchar by 15 TECU or more at low-latitudes, 5 TEC or more at mid-latitudes, and 1 TECU or more at high-latitudes.     

Relative amplitude of the total electron content variations depending on geomagnetic activity

We developed a method of estimation of a relative amplitude dI/I of the total electron content (TEC) variations in the ionosphere as deduced from the data of the global GPS receivers network. To obtain statistically significant results we picked out three latitudinal belts provided in the Internet by the maximum number of GPS sites. They are high-latitudinal belt (50–80°N, 200–300°E; 59 sites), mid latitude belt (20–50°N, 200–300°E; 817 sites), and equatorial belt (±20°N, 0–360°E; 76 sites). The results of the analysis of the diurnal and latitudinal dependencies of dI/I and dI/I distribution probability for 52 days with different levels of geomagnetic activity are presented. It was found that on average the relative amplitude of the TEC variations varies within the range 0–10% proportionally to the value of the K_p geomagnetic index. In quiet conditions the relative amplitude dI/I of the TEC variations at night significantly exceeds the daytime relative amplitude. At high levels of magnetic field disturbances, the geomagnetic control of the amplitude of TEC variations at high and middle latitudes is much more significant than the regular diurnal variations. At the equatorial belt, on average, the amplitude of TEC variations in quiet and disturbed periods almost does not differ. The obtained results may be useful for development of the theory of ionospheric irregularities.     

Report on the Kiso cometary dust trail survey

Cometary dust trails were first observed by IRAS; they are widely known to be the origins of meteoric showers. A new window has been opened for the study of dust trails, using ground-based observations. We succeeded in obtaining direct images of the 22P/Kopff dust trail with the Kiso 1.05-m Schmidt telescope. Following this initial success, we have continued to perform a dust trail survey at Kiso. As a result of this survey, we have detected dust trails along the orbit of six periodic comets, between February 2002 and March 2004. The optical depth of these dust trails are 10⁻⁹ to 10⁻⁸, which is consistent with IRAS measurements. In this paper, we describe the observations and data reduction procedures, and report the brief result obtained between February 2002 and March 2004.     

一种可抑制三次谐波的五相混合式步进电机空间矢量脉宽调制优化控制策略

五相混合式步进电机因其定位精度高、转矩脉动小以及转矩密度高等优点而应用愈加广泛。为进一步减小其运行噪声、提高系统整体效率,本文提出一种可抑制三次谐波的五相混合式步进电机空间矢量脉宽调制(Space vector pulse width modulation,SVPWM)优化控制策略。首先搭建了考虑互感条件下的五相混合式步进电机解耦控制数学模型。在此基础上,采用相邻两大矢量和两中矢量合成混合电压矢量的SVPWM控制策略,并给出了相应的五相混合式步进电机双闭环控制系统。本文所提的混合SVPWM算法通过利用大、中矢量在基波、三次谐波双坐标系对应关系以达到抑制三次谐波的效果,能够有效降低相电流纹波和谐波含量。最后,研制了一套基于FPGA的实验平台,对所提策略与传统滞环算法、大矢量SVPWM算法进行对比研究,重点分析了3种算法对相电流纹波和谐波含量的影响,从而验证了所提算法的正确性和优越性。