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.     

Intrinsic time scale for reconnection on the dayside magnetopause

Recently much attention has been focused on the transient behavior of the magnetopause in response to pressure pulses and southward fluctuations of the interplanetary magnetic field. We examine the motion of the magnetopause behind the foreshock and conclude that this motion is affected by foreshock pressure variations but not by fluctuations in the direction of the magnetic field. Neither magnetopause erosion nor flux transfer event occurrence is controlled by the foreshock. On the contrary, flux transfer events occur at times of steady IMF and thier quasi-periodic behavior is controlled by the magnetopause or the magnetosphere and is not driven by the external boundary conditions. Since flux transfer events are clearly due to reconnection, this observation implies that the IMF must be southward some time perhaps as long as 7 minutes before flux transfer begins.     

Neutral Atmospheres

This paper summarizes the understanding of aeronomy of neutral atmospheres in the solar system, discussing most planets as well as Saturn’s moon Titan and comets. The thermal structure and energy balance is compared, highlighting the principal reasons for discrepancies amongst the atmospheres, a combination of atmospheric composition, heliocentric distance and other external energy sources not common to all. The composition of atmospheres is discussed in terms of vertical structure, chemistry and evolution. The final section compares dynamics in the upper atmospheres of most planets and highlights the importance of vertical dynamical coupling as well as magnetospheric forcing in auroral regions, where present. It is shown that a first order understanding of neutral atmospheres has emerged over the past decades, thanks to the combined effects of spacecraft and Earth-based observations as well as advances in theoretical modeling capabilities. Key gaps in our understanding are highlighted which ultimately call for a more comprehensive programme of observation and laboratory measurements.     

Development of thermal sensors and drilling systems for lunar and planetary regoliths

We consider some novel concepts for thermal properties experiments aboard lunar landers or rovers, that may lead to an improved understanding of both the structure of the lunar near surface layers and the lunar thermal history. The new instruments could be developed using the experience and heritage from recently developed systems, like the Rosetta Lander thermal conductivity experiment MUPUS and existing designs used for terrestrial measurements of thermal conductivity. We describe shortly the working principle of such sensors and the main challenges faced when using them in the airless regolith layers of the Moon or other airless bodies. In addition new concepts to create appropriate drill holes for thermal and other measurements in the lunar regolith are discussed.     

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.     

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.     

A steady-state model for the D- to F-region of the polar cap

For obvious reasons the ionosphere of the polar cap, surrounded by the auroral zone, is only poorly investigated. Even ionosonde data are very scant from geomagnetic latitudes beyond 70°. Since 1997 the European incoherent scatter radar facility EISCAT has an additional installation on Svalbard and has been providing electron density data nearly continuously ever since. These measurements which mainly cover the E- and F-regions are supplemented by rocket data from Heiss Island at a comparable magnetic latitude; these data are more sporadic, but cover lower altitudes and densities. A provisional, steady-state, neural network-based model is presented which uses the data of both sites.