Trapped antiprotons produced by cosmic rays in the Earth's magnetosphere.

The existence of significant fluxes of antiparticles in the Earth magnetosphere has been predicted on theoretical considerations in this article. These antiparticles (positrons or antiprotons) at several hundred kilometers of altitudes, we believe are not of direct extraterrestrial origin, but are the natural products of nuclear reactions of the high energy primary cosmic rays (CR) and trapped protons (TP) confined in the terrestrial radiation belt, with the constituents of terrestrial atmosphere. Extraterrestrial positrons and antiprotons born in nuclear reactions of the same CR particles passing through only 5-7 g/cm2 of interstellar matter, exhibit lower fluxes compared to the antiprotons born at hundreds of g/cm2 in the atmosphere, which when confined in the magnetic field of the Earth (in any other planet), get accumulated. We present the results of the computations of the antiproton fluxes at 10 MeV to several GeV energies due to CR particle interactions with the matter in the interstellar space, and also with the residual atmosphere at altitudes of approximately 1000 km over the Earth's surface. The estimates show that the magnetospheric antiproton fluxes are greater by two orders of magnitude compared to the extraterrestrial fluxes measured at energies <1-2 GeV.     

Estimation of long-term cosmic ray intensity variation in near future and prediction of their contribution in expected global climate change

On the basis of results obtained in our paper [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], we determine: the dimension of the Heliosphere (modulation region), radial diffusion coefficient and other parameters of convection–diffusion; drift mechanisms of long-term variations of cosmic ray (CR) dependence on particle energy; level of solar activity (SA); and generally, the solar magnetic field. We obtain this important information on the basis of CR and SA data in the past, taking into account the theory of convection–diffusion and global drift modulation of galactic CR in the Heliosphere. By using these results and other regularly published predictions of expected SA variation in the near future, as well as predictions of the next SA cycle, we may make predictions of long-term cosmic ray intensity variation expected in the near future (up to 10–12 years). In [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], properties of connections between long-term variation in CR intensity and some part of a global climate change were estimated, controlled by solar activity through CR. We show that in this way we may make predictions of some part of a global climate change expected in the near future (up to 10–12 years and maybe more, depending upon the period during which definite predictions of SA can be made), controlled by solar activity through CR. In this case, estimations of expected long-term changes in the planetary distribution of cutoff rigidities, which also influence CR intensity, as well as CR-influenced effects on global climate variation, become important.     

Comparisons of several transport models in their predictions in typical space radiation environments

We have used several transport codes to calculate dose and dose equivalent values as well as the particle spectra behind a slab or inside a spherical shell shielding in typical space radiation environments. Two deterministic codes, HZETRN and UPROP, and two Monte Carlo codes, FLUKA and Geant4, are included. A soft solar particle event, a hard solar particle event, and a solar minimum galactic cosmic rays environment are considered; and the shielding material is either aluminum or polyethylene. We find that the dose values and particle spectra from HZETRN are in general rather consistent with Geant4 except for neutrons. The dose equivalent values from HZETRN and Geant4 are not far from each other, but the HZETRN values behind shielding are often lower than the Geant4 values. Results from FLUKA and Geant4 are mostly consistent for considered cases. However, results from the legacy code UPROP are often quite different from the other transport codes, partly due to its non-consideration of neutrons. Comparisons for the spherical shell geometry exhibit the same qualitative features as for the slab geometry. In addition, results from both deterministic and Monte Carlo transport codes show that the dose equivalent inside the spherical shell decreases from the center to the inner surface and this decrease is large for solar particle events; consistent with an earlier study based on deterministic radiation transport results. This study demonstrates both the consistency and inconsistency among these transport models in their typical space radiation predictions; further studies will be required to pinpoint the exact physics modules in these models that cause the differences and thus may be improved.     

Turbo码最优周期交织器

从码字重量分布的角度分析了影响Turbo码性能的原因,介绍了Turbo码最优周期交织器的概念,然后基于这一概念,为Turbo码提供了一个新的交织器设计方法, 并对采用这种交织器的Turbo码的性能进行仿真.与几种常用的交织器相比,这样的交织器使Turbo码的性能有了很大改善,这种改善在高信噪比的环境中更为明显.最优周期交织器的应用对于实际的Turbo码交织器设计有重要的意义,对于提高移动通信系统的性能有重要价值.     

Results of measurements of the response of a delta-doped CCD to neutral and charged particle beams

We report the measurements of the response of a delta-doped Charge Coupled Device (CCD) in imaging mode to beams of charged and neutral particles. That is, the detector imaged the incident beam over its 1024 × 1024 pixels, integrating the number of particles counted in each pixel during the exposure period. In order to count individual particles the exposure time would have had to be reduced considerably compared to the typical ?5 s used in these studies. Our CCD thus operated in a different manner than do conventional particle detectors such as the CEM and MCP that normally are used in a particle counting mode. The measurements were carried out over an energy range from 0.8 to 30 keV. The species investigated include H, H⁺, He⁺, N⁺, N₂⁺, and Ar⁺. The energy and ion mass covered wider ranges than previous measurements for the CCD. The results of these measurements show, as in the case of the previous measurement, for a given ion the CCD response increases with energy and for a given particle energy the response decreases with increasing mass of the particle. These results are in agreement with predictions of the theory of the range of ions in solids. The results also show the possibility for the application of the delta doped CCD as a detector for low energy particle measurements for space plasma physics applications.     

Ground-based measurements of galactic cosmic ray fragmentation in shielding.

The mean free path for nuclear interactions of galactic cosmic-rays is comparable to shielding and tissue thicknesses present in human interplanetary exploration, resulting in a significant fraction of nuclear reaction products at depth. In order to characterize the radiation field, the energy spectrum, the angular distribution, and the multiplicity of each type of secondary particles must also be known as a function of depth. Reactions can take place anywhere in a thick absorber; therefore, it is necessary to know these quantities as a function of particle energy for all particles produced. HZE transport methods are used to predict the radiation field; they are dependent on models of the interaction of man-made systems with the space environment to an even greater extent than methods used for other types of radiation. Hence, there is a major need to validate these transport codes by comparison with experimental data. The most cost-effective method of validation is a comparison with ground-based experimental measurements. A research program to provide such validation measurements using neon, iron and other accelerated heavy ion beams will be discussed and illustrated using results from ongoing experiments and their comparison with current transport codes. The extent to which physical measurements yield radiobiological predictions will be discussed.     

Gamma-ray astronomy and cosmic ray origin problem

We analyse the results of the nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) in order to describe their relevant properties: the remnant dynamics and the characteristics of nonthermal emission produced by CRs. It is shown that the theory fits the existing data in a satisfactory way and that the magnetic field in SNRs is significantly amplified due to efficient acceleration of the nuclear CR component. From the fact that magnetic field amplification occurs in all the young SNRs for which relevant data exist, and given the strong theoretical connection between magnetic field amplification and efficiently accelerating the nuclear CRs, we tentatively conclude that the Galactic SNRs are the source population of the Galactic CRs. Due to high interior magnetic fields in young SNRs the π⁰-decay γ-rays generated by the nuclear CR component as a rule dominate over γ-rays generated by the electron CR component, and the calculated γ-ray flux fits existing data.     

Propagation of energetic particles in the solar wind

The characteristics of solar energetic particle events as observed in interplanetary space depend on many physical processes acting at the source and during the transport from the source to the observer. These processes, such as acceleration at the Sun and the propagation near the Sun and in interplanetary space depend, in general, on both the particle velocity and rigidity. Thus, the evaluation of both the nuclear charge and/or atomic mass and the ionic charge of heavy ions turns out to be essential for the interpretation of the physical parameters observed, such as the energy spectra and the compositional variations during individual solar energetic particle events. In this paper recent results on the direct determination of the charge states of He, C, O, and Fe will be summarized. Using these results the compositional variations during individual solar particle events will be discussed. It will be shown that ratio changes by a factor of ～ 10 during the onset phase of solar particle events, as frequently observed, could be explained not only by rigidity dependent interplanetary propagation, but also by rigidity dependent diffusive propagation in the corona. However, there is now increasing experimental evidence that also other processes, such as compositional variations at the source and discontinuities of the interplanetary magnetic field, separating two different particle populations, may be important. Thus the picture emerges that these variations do not have a unique explanation but rather that each event has to be investigated individually.     

飞行时间成像系统的数字仿真技术

飞行时间传感器的仿真主要用于评估基于光相位信息的深度数据的处理算法，现有的方法包括基于对场景三维模型生成的特定噪声模式的仿真模型和基于光信号传播的物理过程仿真模型两大类，为飞行时间成像系统建立了特定的仿真平台.但现有的工作并未考虑环境光干扰对测量数据间的非线性耦合关系，而这种关系是附着/着陆敏感器的重要设计要素之一.本文针对小行星附着任务的需求，研究并建立了包含环境光矢量的飞行时间成像仿真环境，验证了包括环境光模型的飞行时间成像仿真方法，对于飞行时间成像技术在光照条件下的成像性能的评估具有重要意义和应用前景.     

Evolution of Dst index,cosmic rays and global temperature during solar cycles 20–23

We have studied conditions in interplanetary space, which can have an influence on galactic cosmic ray (CR) and climate change. In this connection the solar wind and interplanetary magnetic field parameters and cosmic ray variations have been compared with geomagnetic activity represented by the equatorial Dst index from the beginning 1965 to the end of 2012. Dst index is commonly used as the solar wind–magnetosphere–ionosphere interaction characteristic. The important drivers in interplanetary medium which have effect on cosmic rays as CMEs (coronal mass ejections) and CIRs (corotating interaction regions) undergo very strong changes during their propagation to the Earth. Because of this CMEs, coronal holes and the solar spot numbers (SSN) do not adequately reflect peculiarities concerned with the solar wind arrival to 1 AU. Therefore, the geomagnetic indices have some inestimable advantage as continuous series other the irregular solar wind measurements. We have compared the yearly average variations of Dst index and the solar wind parameters with cosmic ray data from Moscow, Climax, and Haleakala neutron monitors during the solar cycles 20–23. The descending phases of these solar cycles (CSs) had the long-lasting solar wind high speed streams occurred frequently and were the primary contributors to the recurrent Dst variations. They also had effects on cosmic rays variations. We show that long-term Dst variations in these solar cycles were correlated with the cosmic ray count rate and can be used for study of CR variations. Global temperature variations in connection with evolution of Dst index and CR variations is discussed.     

A comparison of total reaction cross section models used in particle and heavy ion transport codes

Understanding the interactions and propagations of high energy protons and heavy ions are essential when trying to estimate the biological effects of Galactic Cosmic Rays (GCR) and Solar Particle Events (SPE) on personnel in space. To be able to calculate the shielding properties of different materials and radiation risks, particle and heavy ion transport codes are needed. In all particle and heavy ion transport codes, the probability function that a projectile particle will collide within a certain distance x in the matter depends on the total reaction cross sections, and the calculated partial fragmentation cross sections scale with the total reaction cross sections. It is therefore crucial that accurate total reaction cross section models are used in the transport calculations. In this paper, different models for calculating nucleon–nucleus and nucleus–nucleus total reaction cross sections are compared with each other and with measurements. The uncertainties in the calculations with the different models are discussed, as well as their overall performances with respect to the available experimental data. Finally, a new compilation of experimental data is briefly presented.     

Propagation of secondary antiprotons and cosmic rays in the Galaxy

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratios (e.g., B/C) produce too few antiprotons. In the present paper, we discuss one possibility to overcome these difficulties. Using the measured antiproton flux and B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei as measured in the heliosphere may contain a fresh local “unprocessed” component at low energies perhaps associated with the Local Bubble, thus decreasing the measured secondary to primary nuclei ratio. The independent evidence for SN activity in the solar vicinity in the last few Myr supports this idea. The model reproduces antiprotons, B/C ratio, and elemental abundances up to Ni (Z ⩽ 28). Calculated isotopic distributions of Be and B are in perfect agreement with CR data. The abundances of three “radioactive clock” isotopes in CR, ¹⁰Be, ²⁶Al, ³⁶Cl, are all consistent and indicate a halo size z_h ∼ 4 kpc based on the most accurate data taken by the ACE spacecraft.     

The photometric study of light scattering from the surface of alumina powder and interpretations by Hapke formula

A laboratory experiment helps to understand the light scattering property of regolith like samples with known compositions and other physical parameters. The laboratory data so obtained can be compared with the existing in situ data on celestial objects like asteroids. Further, it may be analyzed with the help of various theoretical models to understand the light scattering processes from regolith more clearly. In this work we have performed laboratory based photometry of the light scattered from the surfaces of powdered alumina (Al₂O₃) at various tilt angles of the sample and at large phase angles, with the particles having diameter 0.3 μm. The wavelength of observation was 632.8 nm. These data have been fitted by a surface scattering model originally suggested by Hapke. Instead of using empirical Henyey–Greenstein phase function to fix the values of albedo and phase function to be used within Hapke formula, we have used Mie theory for the same. This approach helped us to determine the single particle properties such as particle diameter and complex refractive index from surface scattering phase curve alone. Mie theory depends only on the size parameter X(=2π(radius/wavelength)) and complex refractive index (n, k) of the material. Since the absorption coefficient (k) for alumina is known to be very low but not exactly zero, the best fit to the experimental data was obtained by least square technique with k as a free parameter, as the other parameters are known. Finally, we compare our results with other published results and discuss the scope of application of the method we adopted.     

Three-dimensional global simulation of multiple ICMEs’ interaction and propagation from the Sun to the heliosphere following the 25–28 October 2003 solar events

This study performs simulations of interplanetary coronal mass ejection (ICME) propagation in a realistic three-dimensional (3D) solar wind structure from the Sun to the Earth by using the newly developed hybrid code, HAFv.2+3DMHD. This model combines two simulation codes, Hakamada–Akasofu–Fry code version 2 (HAFv.2) and a fully 3D, time-dependent MHD simulation code. The solar wind structure is simulated out to 0.08 AU (18 Rs) from source surface maps using the HAFv.2 code. The outputs at 0.08 AU are then used to provide inputs for the lower boundary, at that location, of the 3D MHD code to calculate solar wind and its evolution to 1 AU and beyond. A dynamic disturbance, mimicking a particular flare’s energy output, is delivered to this non-uniform structure to model the evolution and interplanetary propagation of ICMEs (including their shocks). We then show the interaction between two ICMEs and the dynamic process during the overtaking of one shock by the other. The results show that both CMEs and heliosphere current sheet/plasma sheet were deformed by interacting with each other.     

Depth dependency of neutron density produced by cosmic rays in the lunar subsurface

Depth dependency of neutrons produced by cosmic rays (CRs) in the lunar subsurface was estimated using the three-dimensional Monte Carlo particle and heavy ion transport simulation code, PHITS, incorporating the latest high energy nuclear data, JENDL/HE-2007. The PHITS simulations of equilibrium neutron density profiles in the lunar subsurface were compared with the measurement by Apollo 17 Lunar Neutron Probe Experiment (LNPE). Our calculations reproduced the LNPE data except for the 350–400 mg/cm² region under the improved condition using the CR spectra model based on the latest observations, well-tested nuclear interaction models with systematic cross section data, and JENDL/HE-2007.     

The beginning of cosmic ray astrophysics

A short history of the beginning of cosmic ray (CR) astrophysics is considered: from the hypothesis on CR origin as a result of Supernova explosions in the Metagalaxy, to a model of solar origin of CR, galactic origin based on the stochastic mechanism of charged particle acceleration in interstellar space, to extragalactic and hierarchical models of CR origin, as well as galactic CR origin taking into account radio-astronomical data. We consider also the first balloon results on the chemical contents of primary CR (especially of the contents Li, Be, B), important for any model of CR origin. Investigations of the injection problem, CR drift and diffusion acceleration by shock waves, and CR generation in Supernova remnants were also important steps in the beginning of CR astrophysics.     

Klobuchar电离层模型误差分析及预测

电离层延迟误差是全球导航卫星系统(global navigation satellite system,GNSS)中的重要误差源之一.目前在电离层延迟改正模型中,应用最广泛的是Klobuchar参数模型,但是该模型的改正率仅能达到60％左右,无法满足日益增长的精度需求.将国际GNSS监测评估系统(internation...     

Modeling the charging of geosynchronous and interplanetary spacecraft using Nascap-2k

Nascap-2k is the updated version of the NASCAP/GEO spacecraft charging analysis code. In addition to packaging the physical content of NASCAP/GEO in a modern way, Nascap-2k incorporates other plasma analysis codes (in particular, the DynaPAC code) in order to extend its applicability to a wide variety of plasma environments. Nascap-2k also includes an interactive Object Toolkit for defining spacecraft surface models for analysis. In this paper we focus on the tenuous plasma charging capabilities of the code, with application to DSCS-III (geosynchronous environment), STEREO (solar wind environment) and MESSENGER (solar wind environment near 0.4 AU).     

The FLUKA code for space applications: recent developments.

The FLUKA Monte Carlo transport code is widely used for fundamental research, radioprotection and dosimetry, hybrid nuclear energy system and cosmic ray calculations. The validity of its physical models has been benchmarked against a variety of experimental data over a wide range of energies, ranging from accelerator data to cosmic ray showers in the earth atmosphere. The code is presently undergoing several developments in order to better fit the needs of space applications. The generation of particle spectra according to up-to-date cosmic ray data as well as the effect of the solar and geomagnetic modulation have been implemented and already successfully applied to a variety of problems. The implementation of suitable models for heavy ion nuclear interactions has reached an operational stage. At medium/high energy FLUKA is using the DPMJET model. The major task of incorporating heavy ion interactions from a few GeV/n down to the threshold for inelastic collisions is also progressing and promising results have been obtained using a modified version of the RQMD-2.4 code. This interim solution is now fully operational, while waiting for the development of new models based on the FLUKA hadron-nucleus interaction code, a newly developed QMD code, and the implementation of the Boltzmann master equation theory for low energy ion interactions.     

基于UPF的修正罗德里格参数卫星姿态估计

粒子滤波(PF)在卫星姿态估计中解决了非高斯分布的情况,在粒子滤波的基础上,利用修正罗德里格参数(MRPs)表示姿态,用动力学方程进行角速率的传播,提出将采样粒子滤波(UPF)算法应用到基于MRPs的卫星姿态估计中,该方法通过UKF得到粒子滤波的重要性采样函数,更充分地利用了最新的量测信息。仿真结果表明,这种基于UPF的姿态估计方法能够获得比UKF和标准PF更高的滤波精度,并且具有更好的滤波稳定性。