Research on a complex CME event including Hα, LASCO, radio and MDI observations

We present our research on a fast and decelerating partial halo coronal mass ejection (CME) event detected in multi-wavelengths in the chromosphere and the corona on 14 October, 1999. The event involved a whole complex active area which spanned more than 40° of heliolongitude. It included a strong solar flare (XI/1N) and a complex eruptive filament within an active region of the entire complex. Especially, several radio sources were detected in the decimetric range prior to the CME by the Nançay Radioheliograph (NRH). A linear force-free field extrapolation of the Michelson Doppler Imager (MDI) magnetogram was performed to calculate the magnetic topology of the complex prior to the triggering of the event. The presence of a coronal null point combined with the occurrence of two distant and nearly simultaneous radio sources put strong arguments in favor of the generalized breakout model for the triggering of the eruption. The analysis of the subsequent development of the event suggests that large interconnecting loops were ejected together with the CME.     

Non-Thermal Processes in Cosmological Simulations

Non-thermal components are key ingredients for understanding clusters of galaxies. In the hierarchical model of structure formation, shocks and large-scale turbulence are unavoidable in the cluster formation processes. Understanding the amplification and evolution of the magnetic field in galaxy clusters is necessary for modelling both the heat transport and the dissipative processes in the hot intra-cluster plasma. The acceleration, transport and interactions of non-thermal energetic particles are essential for modelling the observed emissions. Therefore, the inclusion of the non-thermal components will be mandatory for simulating accurately the global dynamical processes in clusters. In this review, we summarise the results obtained with the simulations of the formation of galaxy clusters which address the issues of shocks, magnetic field, cosmic ray particles and turbulence.     

2000年和2001年52°N地区OI5577气辉强度与原子氧浓度峰值的时间变化特征

分析了2000年和2001年期间52°N地区OI5577气辉强度的夜间变化特征和季节变化特征．利用由 OI5577气辉强度反演原子氧浓度峰值的方法反演出原子氧浓度的峰值,分析了峰值的夜间变化特征和季节变化特征．结果表明, OI5577气辉强度的夜间变化特征随季节变化, 2000年春季的夜间强度最大值出现在0000LT 之后,夏季和秋季的出现在0000LT之前,冬季的出现在0000LT,2001年春季和秋季的夜间强度最大值出现在0000LT之前,夏季和冬季的出现在0000LT;OI5577气辉强度在2000年2月份,8月份和10月份出现最大值,在2001年9月份有最大值．就主要特征而言,反演出的原子氧浓度峰值的夜间变化特征和季节变化特征分别与OI5577气辉强度的一致．     

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.     

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.