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.     

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.     

WIND/3DP and Nancay radioheliographobservations of solar energetic electron events

Initial results of a combined study of electron events using the 3DP experiment on the WIND spacecraftand the Nançay Radioheliograph (NRH) are presented. A total of 57 electron events whose solar release time could be inferred from WIND/3DP observations occurred during NRH observing times. In 40 of them a distinct signature was detected in maps at decimetric and metric wavelengths (dm-m-λ) taken by the NRH. These events are equally distributed among two categories: (1) Electron release together with dm-m-λ bursts of a few minutes duration: these events are also accompanied by decametric-hectometric type III bursts seen by WAVES/WIND. They correspond to the well-known impulsive electron events. (2) Electron release during long duration (several tens of minutes) dm-m-λ emission: the electrons are most often released more than ten minutes after the start of the radio event. In the majority of cases the dm-m-λ radio source changes position, size, and/or intensity near the time of electron release.     

An empirical model of the plasma environment around Mercury

Since the flyby observations by Mariner 10 in 1974 and 1975, Mercury has been one of the most interesting objects for space physics and planetary exploration. The MESSENGER and BepiColombo missions now plan to revisit this planet. In order to design plasma instruments for the BepiColombo mission, we have estimated electron and ion fluxes around Mercury with an empirical model, which has been developed for the Earth’s magnetotail. The solar wind data needed as input parameters are derived from Helios observations. The result shows that our predicted electron fluxes at aphelion agree well with the Mariner-10 data. It is also noted that ion instruments must cover a very wide dynamic range of proton fluxes. However, the applicability of the Earth’s magnetospheric model to Mercury is, in itself, an important issue for comparative magnetospheric studies.     

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.     

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.     

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.     

Future of Space Astronomy: A global Road Map for the next decades

The use of space techniques continues to play a key role in the advance of astrophysics by providing access to the entire electromagnetic spectrum from radio to high energy γ rays. The increasing size, complexity and cost of large space observatories places a growing emphasis on international collaboration. Furthermore, combining existing and future datasets from space and “ground based” observatories is an emerging mode of powerful and relatively inexpensive research to address problems that can only be tackled by the application of large multi-wavelength observations. While the present set of astronomical facilities is impressive and covers the entire electromagnetic spectrum, with complementary space and “ground based” telescopes, the situation in the next 10–20 years is of critical concern. The James Webb Space Telescope (JWST), to be launched not earlier than 2018, is the only approved future major space astronomy mission. Other major highly recommended space astronomy missions, such as the Wide-field Infrared Survey Telescope (WFIRST), the International X-ray Observatory (IXO), Large Interferometer Space Antenna (LISA) and the Space Infrared Telescope for Cosmology and Astrophysics (SPICA), have yet to be approved for development.     

Digital signal processing and numerical analysis for radar in geophysical applications

Numerical solutions for signal processing are described in this work as a contribution to study of echo detection methods for ionospheric sounder design. The ionospheric sounder is a high frequency radar for geophysical applications. The main detection approach has been done by implementing the spread-spectrum techniques using coding methods to improve the radar’s range resolution by transmitting low power. Digital signal processing has been performed and the numerical methods were checked. An algorithm was proposed and its computational complexity was calculated.