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.     

Great solar bursts of October 19, 22 and 23, 1989

Itapetinga measurements at 48 GHz with the multibeam technique are used to determine the relative position of solar burst centroid of emission with high spatial accuracy and time resolution. For the Great Bursts of October 19,22, 1989, with a large production of relativistic particles, and October 23, it is suggested that, at 48 GHz, the bursts might have originated in more then one source in space and time. Additionally the October 19 and 22 Ground Level Events exhibited very unusual intensity-time profiles including double component structures for the onset phase. The Bern observatory spectral radio emission data show a strong spectral flattening typical for large source inhomogeneties. The interpretation for this is that large solar flares are a superposition of a few strong bursts (separated both in space and time) in the same flaring region.     

Models of P/Borrelly: Activity and dust mantle formation

Comet 19P/Borrelly was observed by Deep Space One spacecraft on September 22, 2001 (Soderblom et al., 2002).The DS1 images show a very dark and elongate nucleus with a complex topography; the IR spectra show a strong red-ward slope consistent with a very hot and dry surface (345K to 300K). During DS1 encounter the comet coma was dominated by a prominent jet but most of the comet was inactive, confirming the Earth-based observations that <10% of the surface is actively sublimating. We have developed a thermal evolution model of comet PBorrelly, using a numerical code that is able to solve the heat conduction and gas diffusion equations at the same time across an idealized spherical nucleus ( De Sanctis et al., 1999, 2000; Capria et al., 2000; Coradini et al., 1997a,b). The comet nucleus is composed by water, volatiles ices and dust in different proportions. The refractory component is made by grains that are embedded in the icy matrix. The code is able to account for the dust release, contributing to the dust flux, and the formation of dust mantles on the comet surface. The model was applied to a cometary nucleus with the estimated physical and dynamical characteristics of P/Borrelly in order to infer the status and activity level of a body on such an orbit during the DS1 observation. The comet gas flux, differentiation and thermal behavior were simulated and reproduced. The model results are in good agreement with the DS1 flyby results and the ground based observations, in terms of activity, dust coverage and temperatures of the surface.     

Measurements of bidirectional reflectance of ordinary chondrite for muses-C in-situ detection

MUSES-C, a Japanese sample return mission, is targeting a small near Earth asteroid, 1998SF36, which is considered an S-type asteroid and is similar in spectroscopy to LL class ordinary chondrite meteorite ([Binzel et al., 2001]). Although this mission will bring us detailed photometric data, that is, disk-resolved bidirectional reflectance data of the asteroid, there were few bidirectional reflectance data of ordinary chondrite meteorites. For the purpose of comparison with the data obtained by the in-situ observation, we have performed measurements of bidirectional reflectance of ordinary chondrite samples.

Here we summarize the result of our laboratory measurements of the bidirectional reflectance and compare them with the scattering property of 1998SF36. Although the geometric albedo of 1998SF36 is higher than the typical value of S-type asteroids, however, the laboratory data of ordinary chondrite are similar to or brighter than the model disk-resolved reflectance of 1998SF36 derived from disk-integrated ground-based data. We found in our laboratory data that there is a positive correlation between the surface roughness and the strength of the opposition effect. A stronger opposition effect in the reflectance of 1998SF36 than the laboratory data may therefore indicate that the surface of the asteroid has rougher structure than our laboratory samples.     

The acceleration of Anomalous Cosmic Rays by stochastic acceleration in the heliosheath

Stochastic acceleration in the heliosheath appears to be a likely mechanism by which Anomalous Cosmic Rays (ACRs) are accelerated. However, most stochastic acceleration mechanisms are not appropriate. The energy density in the ACRs and in the interstellar pickup ions out of which the ACRs are accelerated greatly exceeds the energy density in the turbulence in the heliosheath. Thus, a traditional stochastic acceleration mechanism in which particles are accelerated by damping the turbulence will not work. A stochastic acceleration mechanism has been developed in which the total energy of the pickup ions and the ACRs is conserved. Energy is redistributed from the core pickup ions into a suprathermal tail to create the ACRs. A model for the acceleration of the ACRs in the heliosheath, based on this stochastic acceleration mechanism, is presented. The model provides reasonable fits to the spectra of suprathermal particles and ACRs observed by Voyager.     

磁暴对海南地区电离层扩展F的影响

利用2002年2月至2007年12月海南地区DPS-4测高仪观测数据, 用统计分析方法研究了磁暴对电离层扩展F的影响. 结果认为磁暴从整体上抑制了扩展F现象的发生. 但若把扩展F分为不同类型, 则结果却有所不同. 对于频率型(FSF), 在2002年和2003年磁暴对其有促发作用, 在2004---2007年有抑制作用; 对于区域型(RSF), 在2002---2005年磁暴对其有抑制作用, 在2006年和2007年, 对其有弱促发作用; 对于混合型(MSF), 在2002年磁暴对其有抑制作用, 在2003年和2004年有促发作用, 从2005年开始, 磁暴对其有抑制作用; 对于强区域型 (SSF), 在2002---2004年磁暴对其有抑制作用, 在2005年和2006年有促发作用, 2007年有弱抑制作用.      

Understanding Substorms in the Magnetotail: Early Development and Recent Progress

This is a concise review of physics of the substorm in the magnetotail. It consists of two parts. The first part summarizes historical developments in the early days of the space age (1960--1975) when the basic concepts such as magnetotail and reconnection were established and the leading model of the substorm was introduced. The second part is an overview of the research conducted in recent years (1995--2010) when very significant advances have been achieved in understanding the substorm physics by virtue of several major satellites missions that addressed the magnetotail physics intensively.      

Ionospheric response of equatorial and low latitude F-region during the intense geomagnetic storm on 24–25 August 2005

In this investigation, we present and discuss the response of the ionospheric F-region in the South American and East Asian sectors during an intense geomagnetic storm in August 2005. The geomagnetic storm studied reached a minimum Dst of −216 nT at 12:00 UT on 24 August. In this work ionospheric sounding data obtained of 24, 25, and 26 August 2005 at Palmas (PAL; 10.2° S, 48.2° W; dip latitude 6.6° S), São José dos Campos (SJC, 23.2° S, 45.9° W; dip latitude 17.6° S), Brazil, Ho Chi Minh City, (HCM; 10.5° N, 106.3° E; dip latitude 2.9° N), Vietnam, Okinawa (OKI; 26.3° N, 127.8° E; dip latitude 21.2° N), Japan, are presented. Also, the GPS observations obtained at different stations in the equatorial and low-latitude regions in the Brazilian sector are presented. On the night of 24–25 August 2005, the h′F variations show traveling ionospheric disturbances associated with Joule heating in the auroral zone from SJC to PAL. The foF2 variations show a positive storm phase on the night of 24–25 August at PAL and SJC during the recovery phase. Also, the GPS-VTEC observations at several stations in the Brazilian sector show a fairly similar positive storm phase on 24 August. During the fast decrease of Dst (between 10:00 and 11:00 UT) on 24 August, there is a prompt penetration of electric field of magnetospheric origin that result in abrupt increase (∼12:00 UT) in foF2 at PAL, SJC (Brazil) and OKI (Japan) and in VTEC at IMPZ, BOMJ, PARA and SMAR (Brazil). OKI showed strong oscillations of the F-region on the night 24 August resulted to the propagation of traveling atmospheric disturbances (TADs) by Joule heating in the auroral region. These effects result a strong positive observed at OKI station. During the daytime on 25 August, in the recovery phase, the foF2 observations showed positive ionospheric storm at HCM station. Some differences in the latitudinal response of the F-region is also observed in the South American and East Asian sectors.