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.     

Cyclic variations in the solar lower atmosphere

The Ca K line has been measured regularly nearly every month since 1974 at Kitt Peak. It is well known that the K₁ component of the Ca K line is formed in the temperature minimum region (TMR) of the solar atmosphere. Our study of the data of CaII K profiles over two solar cycles indicates that both in full disc integrated spectra and in center disc spectra, the distance between the red K₁ and the blue K₁ of the profiles and its average intensity show periodic variations. But the variation for the full disc integrated spectra fluctuates in the same way as the sunspot number does, while that for the center disc spectra has a time delay with respect to sunspot number. Non-LTE computations yield a cyclic temperature variation of about 17 K of the TMR in the quiet-Sun atmosphere and a cyclic variation of about 15–20 km in the height position of the TMR.     

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.     

A 3NM-64_He added to LARC for Solar Extreme Event studies during solar cycle 24

The Antarctic Laboratory for Cosmic Rays (LARC, acronym for Laboratorio Antartico per i Raggi Cosmici or Laboratorio Antártico para Rayos Cósmicos) operates on King George Island (South Shetlands). Since January 1991 a standard 6NM-64 detector has been recording continuous cosmic ray measurements and several Ground-Level Enhancements have been registered. Here we describe the different phases performed in Italy for the realization of a 3NM-64_³He detector, which started its measurements during the Italian XXII Antarctic Summer Campaign. Data recorded during solar activity cycle 24 will furnish an useful research tool for the next Solar Extreme Events.