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.     

Sharp boundaries of solar wind plasma structures and their relationship to solar wind turbulence

Sharp (<10 min) and large (>20%) solar wind ion flux changes are common phenomena in turbulent solar wind plasma. These changes are the boundaries of small- and middle-scale solar wind plasma structures which can have a significant influence on Earth’s magnetosphere. These solar wind ion flux changes are typically accompanied by only a small change in the bulk solar wind velocity, hence, the flux changes are driven mainly by plasma density variations. We show that these events occur more frequently in high-density solar wind. A characteristic of solar wind turbulence, intermittency, is determined for time periods with and without these flux changes. The probability distribution functions (PDF) of solar wind ion flux variations for different time scales are calculated for each of these periods and compared. For large time scales, the PDFs are Gaussian for both data sets. For small time scales, the PDFs from both data set are more flat than Gaussian, but the degree of flatness is much larger for the data near the sharp flux change boundaries.     

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.     

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.     

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.