平行激波加速离子的数值研究

伴随耀斑和日冕物质抛射共生的日冕和行星际快激波作为一种粒子加速机制一直是理论研究关注的热点课题.在准平行激波传播条件下,首先建立数值求解一维输运方程的方法,然后探讨加速离子分布与激波和背景等离子参数之间的关系.取扩散系数分别为常数和能量的函数、有限自由逃逸边界的计算结果表明:(1)随着加速时间的增大,高能粒子近似呈双幂律分布,低能端(3～10 MeV)谱指数逐渐从10.2减小到2.4,能谱逐渐变硬,粒子被激波加速后能量逐渐增大;(2)随着激波压缩比从2增大到4,相同时间同一能量范围的粒子能谱谱指数逐渐从3.2减小到2.2,能谱逐渐变硬,表明激波强度的增大使得加速效率增大;(3)上下游逃逸边界由5减小到2后,粒子能谱的谱指数由2.4增大到3.3,粒子的加速效率减小;(4)当粒子注入能量增大时,粒子能谱的谱指数由2.4减小到0.9,加速效率增大;(5)当扩散系数与能量成正比时,粒子能谱指数由2.2增大到4.3,能谱变软.     

脉冲型耀斑中阿尔芬波湍流加速3He和重离子机制探讨

基于在^3He丰富事件中，高能^3He和重离子具有相似的幂律谱分布这一观测结果，通过数值求解Fokker-Planck方程，探讨经阿尔芬波湍动速后的离子分布随时间的演化特征。计算结果表明：加速源区的等离子体密度和阿尔芬波湍动能量密度对粒子能谱分布起主要作用，如果取加速源区等离子体密度n=(0.1-1)10^10cm^-3、磁场强度B=50-100Gs、湍动能量密度为0.4-2ergs cm^-3，则在1秒左右的时间内，湍动阿尔芬波能够将^3He和重离子加速到10MeV/nucleon量级，能谱指数为2.0-3.5。理论计算与观测结果基本一致。     

太阳风、高能暴粒子、太阳高能粒子和太阳宇宙线的谱指数研究

本文根据卫星提供的1963—1978年太阳风实验资料,将太阳风中的质子流作为极低能宇宙线,则能得到0.3—4kev的质子积分通量—动能曲线,使低能宇宙线的能谱向前推进了约三个数量级。所得的极低能宇宙线能谱亦呈幂律谱,即:J(>E)=A_sE~(-γ),具有双幂指数,约在1kev处发生转折,与低能太阳宇宙线能谱非常类似。 最近,卫星ISEE—3观测到46次与行星际激波相联系的高能暴粒子(ESP)事例,在能域35—53kev的各次质子峰值强度恰好绘于联结两能谱的虚线之中。这样,从太阳风、ESP、太阳高能粒子(SEP)到太阳低能宇宙线的能谱都被连接了起来,对于它们的起源,也能获得合理地很好地解释。     

射电喷流中湍动阿耳芬波对相对论电子的加速

本文通过数值求解带电粒子与Ａｌｆｖｅｎ波湍动相互作用的动力学方程，得到了相对论电子在射电喷流中被加速随时间演化的解．高能电子可以加速到Ｌｏｒｅｎｔｚ因子γｙ～１０６，且形成稳态的幂律谱，尽管其谱指数Ｓ≈ｌ比观测值小，但粒子加速时间约为１０１２－１０１４秒，小于射电斑的寿命１０７年．粒子能谱指数几乎与Ａｌｆｖｅｎ波谱指数和能量损失函数无关．能量损失对加速上限有较大影响．     

宇宙线的超新星遗迹起源

宇宙线的起源是高能天体物理的核心问题之一.一直以来,超新星爆发被认为是能谱膝区以下宇宙线的主要来源.多波段观测表明,超新星遗迹有能力加速带电粒子至亚PeV (10~(15)eV)能量.扩散激波加速被认为是最有效的天体高能粒子加速机制之一,而超新星遗迹的大尺度激波正好为这一机制提供平台.近年来,一系列较高精度的地面和空间实验极大地推动了对宇宙线以及超新星遗迹的研究.新的观测事实挑战着传统的扩散激波加速模型以及其在银河系宇宙线超新星遗迹起源学说上的应用,深化了人们对宇宙高能现象的认识.结合超新星遗迹辐射能谱的时间演化特性,构建的时间依赖的超新星遗迹粒子加速模型,不仅能够解释200 GV附近宇宙线的能谱反常,还自然地形成能谱膝区,甚至可以将超新星遗迹粒子加速对宇宙线能谱的贡献延伸至踝区.该模型预期超新星遗迹中粒子的输运行为表现为湍流扩散,这需要未来的观测以及与粒子输运相关的等离子体数值模拟工作来进一步验证.     

质子耀斑爆发相的粒子加速

本文从粒子分布函数所满足的带有规则电场的准线性方程出发,得到了包含有规则电场与湍动起伏场相互耦合作用在内的等效动量空间扩散系数,提出了太阳质子耀斑中性片中的规则电场与湍动起伏场的联合加速机制。根据太阳质子耀斑的物理条件,计算表明荷电粒子在中性片中可以被有效地加速,能量可以达到～20MeV,甚至～1GeV。本文证明了离子声湍动起伏场与规则电场的联合加速机制有效地使质子和其它重离子注入到Langmuir湍动加速区中去;并且表明,在Langmuir湍动起伏场与规则电场联合加速的情形下,可以得到与观测事实符合得较好的高能质子的谱以及高能电子的幂律分布的谱。     

回旋同步辐射自吸收的特性

详细研究了高能电子产生的回旋同步辐射自吸收的特性,并用磁偶极子场的 射电微波源模型计算了它的光学厚度.发现: (1)自吸率Kv随谐波数s的增加迅速地 下降,以致只有低次谐波(s<5)上的自吸收才会对微波爆发谱产生实质性的影响; (2) 自吸率κv随暗波数s下降的陡度还随高能电子的能谱指数δ的增加和其低能截止能量 E0的下降而迅速增加; (3)κv值还随传播角的增大而增加,并在70°至75°范围内达到 极值; (4)假设高能电子数密度为103cm-3,则在2≤s≤5范围内的自吸收的光学厚 度τvself在101-10-2之间变化,这些值约比回旋共振吸收的光学厚度τvgyro小3至4 个数量级.τvself在均匀磁场情况下可能被高估.只有当被加速的高能电子的数密度大于 104cm-3时,自吸收的光学厚度才开始能与回旋共振吸收的光学厚度相比较.     

相对论电子的物理效应及其在高能天体中的应用

本文对相对论电子的辐射性质、能谱演化和加速机制等进行了简要的介绍。同时，对相对论电子在高能天体中的辐射作用和特性进行了简要的综述。本文给出了相对论电子在Blazar天体的射电辐射机制、光变机制、BL Lac天体的辐射机制以及γ暴的辐射机制等方面的应用研究成果。1、提出了相对论电子的光学薄同步辐射模型：解释Blazar天体的射电平谱：Blazar中心体的剧烈活动，使射电辐射区处于等离子体湍动状态，其中的相对论电子在湍动等离子体波的二交费米加速、激波加速、辐射损失、粒子逃逸和辐射区的绝热减速等物理过程作用下，形成较平的能谱，产生射电平谱。2、提出了新的Blazar天体光变模型：当Blazar天体爆发时，中心天体产生大量的相对论电子，注入喷流中；相对论电子产生同步辐射，并不断损失能量和逃逸辐射区，使它们的能谱快速变化，引起辐射发生快速光变。3、给出了BL Lac天体的等离子体反应堆模型：大量相对论电子从中心天体注入周围的等离子体反应堆中，通过同步辐射快速损失能量，同时这些电子同步吸收反应堆中不透明的光子，产生一个稳定、各向同性的幂律分布，其谱指数为γ＝3；然后，这些相对论电子通过等离子体反应堆的爆发或其表面扩散过程逃逸出来，辐射低频的同步辐射。模型解释了BL Lac天体的高频辐射表现出快速的谱变化性质，即流量减小时谱变陡。4、提出了相对论电子的内激波加速模型，解释γ暴的尖峰光变特性：在γ暴产生的相对论运动的壳层中，有内激波产生；激波在壳层中传播，耗散壳层的运动能，使其中的部分电子加速成为相对论电子。然后，这些电子通过同步辐射产生观测到的γ辐射。模型认为，γ暴中的每个尖峰辐射是一对内激波加速相对论电子的辐射过程，复杂的γ暴光变曲线是多对内激波辐射过程的叠加。     

Mrk421中TeV和X射线爆发的相关性

用同步自康普顿机制（ＳＳＣ）来解释１９９４年５月在ＢＬＬａｃ天体Ｍｒｋ４２１中观测到的ＴｅＶＸ射线爆发．认为此ＴｅｖＸ射线爆发是由一个高能成分产生的,它具有平坦的电子能量分布,其低能段上能谱指数αｅ≈０２５．这个爆发成分独立于另一个低能成分（或宁静成分）．此宁静成分具有较陡的电子能量分布,它产生宁静态的辐射能谱．大多数射电爆发表明电子能量分布有坦谱,能谱指数αｅ～０２－０４．与这个事实一起,上述结果意味着活动星系核中粒子加速的基本机制产生平坦的电子能量分布．     

磁重联电流片的参数变化对电子加速的影响

通过采用试验粒子的方法,研究了在有引导磁场Bz存在的磁重联电流片中,电子被super-Dreicer电场Ez加速后的运动特征.首先,考虑了引导磁场恒定且与电场有不同方向时对粒子加速的影响.在这种情况下,Bz方向的改变直接改变了电子的运动轨迹,使其沿着不同的路径离开电流片.在Bz和Ez同向时,高能电子的pitch-angle接近于180°.然而,当2者反向时,高能电子的pitch-angle接近0°.引导磁场的取向只是使电场有选择地对不同区域的电子进行加速,不会最终影响电子的能量分布,最终得到的能谱是普遍的幂率谱E-γ.在典型的日冕条件下, γ大约等于2.9.进一步的研究表明γ的大小依赖于引导磁场及磁重联电场的强弱,以及电流片的尺度.随后,也研究了包含多个X-点和O-点电流片中被加速粒子的运动特征.结果表明X-点和O-点的存在使得粒子被束缚在加速区并获得最大的加速,而且最终的能谱具有多幂率谱的特征.     

On the origin of the solar system and the exceptional position of the Sun in the Galaxy

The solar system's position in the Galaxy is an exclusive one, since the Sun is close to the corotation circle, which is the place where the angular velocity of the galactic differential rotation is equal to that of density waves displaying as spiral arms. Each galaxy contains only one corotation circle; therefore, it is an exceptional place. In the Galaxy, the deviation of the Sun from the corotation is very small — it is equal to ΔR/R _⊙≈0.03, where ΔR=R _c ?R _⊙,R _c is the corotation distance from the galactic center andR _⊙ is the Sun's distance from the galactic center. The special conditions of the Sun's position in the Galaxy explain the origin of the fundamental cosmogony timescalesT ₁≈4.6×10⁹ yr,T ₂?10⁸ yr,T ₃?10⁶ yr detected by the radioactive decay of various nuclides. The timescaleT ₁ (the solar system's ‘lifetime’) is the protosolar cloud lifetime in a space between the galactic spiral arms. The timescaleT ₂ is the presolar cloud lifetime in a spiral arm.T ₃ is a timescale of hydrodynamical processes of a cloud-wave interaction. The possibility of the natural explanation of the cosmogony timescales by the unified process (on condition that the Sun is near the state of corotation) can become an argument in favour of the fact that the nearness to the corotation is necessary for the formation of systems similar to the Solar system. If the special position of the Sun is not incidental, then the corotation circles of our Galaxy, as well as those of other galaxies, are just regions where situations similar to ours are likely to be found.     

Perturbations by the oblateness of the Earth and by the planets in the motion of the Moon

Perturbations in the motion of the Moon are computed for the effect by the oblateness of the Earth and for the indirect effect of planets. Based on Delaunay's analytical solution of the main problem, the computations are performed by a method of Fourier series operation. The effect of the oblateness of the Earth is obtained to the second order, partly adopting an analytical evaluation. Both in longitude and latitude are found a few terms whose coefficient differs from the current lunar ephemeris based on Brown's theory by about 0.01. While, concerning the indirect effect of planets, several periodic terms in the current ephemeris seem to have errors reaching 0.05.As for the secular variations of and due to the figure of the Earth and the indirect effect of planets, the newly-computed values agree within 1/cy with Brown's results reduced to the same values of the parameters. Further, the accelerations in the mean longitude, and caused by the secular changes in the eccentricity of the Earth's orbite and in the obliquity of the ecliptic are obtained. The comparison with Brown shows an agreement within 0.3/cy² for the former cause and 0.02/cy² for the latter. An error is found in the argument of the principal term for the perturbations due to the ecliptic motion in the current ephemeris.Proceedings of the Conference on Analytical Methods and Ephemerides: Theory and Observations of the Moon and Planets. Facultés universitaires Notre Dame de la Paix, Namur, Belgium, 28–31 July, 1980.     

On the model of the accumulation of the moon compatible with the data on the composition and the age of lunar rocks

It is suggested that the overall early melting of the lunar surface is not necessary for the explanation of facts and that the structure of highlands is more complicated than a solidified anorthositic ‘plot’. The early heating of the interior of the Moon up to 1000K is really needed for the subsequent thermal history with the maximum melting 3.5 × 10⁹ yr ago, to give the observed ages for mare basalts. This may be considered as an indication that the Moon during the accumulation retained a portion of its gravitational energy converted into heat, which may occur only at rapid processes. A rapid (t < 10³ yr) accretion of the Moon from the circumterrestrial swarm of small particles would give necessary temperature, but it is not compatible with the characteristic time 10⁸ yr of the replenishment of this swarm which is the same as the time-scale of the accumulation of the Earth. It is shown that there were conditions in the circumterrestial swarm for the formation at a first stage of a few large protomoons. Their number and position is evaluated from the simple formal laws of the growth of satellites in the vicinity of a planet. Such ‘systems’ of protomoons are compared with the observed multiple systems, and the conclusion is reached that there could have been not more than 2–3 large protomoons with the Earth. The tidal evolution of protomoon orbits was short not only for the present value of the tidal phase-lag but also for a considerably smaller value. The coalescence of protomoons into a single Moon had to occur before the formation of the observed relief on the Moon. If we accept the age 3.9 × 10⁹ yr for the excavation of the Imbrium basin and ascribe the latter to the impact of an Earth satellite, this collision had to be roughly at 30R, whereR is the radius of the Earth, because the Moon at that time had to be somewhere at this distance. Therefore, the protomoons had to be orbiting inside 20–25R, and their coalescence had to occur more than 4.0x10⁹ yr ago. The energy release at coalescence is equivalent to several hundred degrees and even 1000 K. The process is very rapid (of the order of one hour). Therefore, the model is valid for the initial conditions of the Moon.     

On the theory of the oblateness of the Sun

In this paper we first emphasize why it is important to know the successive zonal harmonics of the Sun's figure with high accuracy: mainly fundamental astrometry, helioseismology, planetary motions and relativistic effects. Then we briefly comment why the Sun appears oblate, going back to primitive definitions in order to underline some discrepancies in theories and to emphasize again the relevant hypotheses. We propose a new theoretical approach entirely based on an expansion in terms of Legendre's functions, including the differential rotation of the Sun at the surface. This permits linking the two first spherical harmonic coefficients (J ₂ and J ₄) with the geometric parameters that can be measured on the Sun (equatorial and polar radii). We emphasize the difficulties in inferring gravitational oblateness from visual measurements of the geometric oblateness, and more generally a dynamical flattening. Results are given for different observed rotational laws. It is shown that the surface oblateness is surely upper bounded by 11 milliarcsecond. As a consequence of the observed surface and sub-surface differential rotation laws, we deduce a measure of the two first gravitational harmonics, the quadrupole and the octopole moment of the Sun: J ₂=−(6.13±2.52)×10⁻⁷ if all observed data are taken into account, and respectively, J ₂=−(6.84±3.75)×10⁻⁷ if only sunspot data are considered, and J ₂=−(3.49±1.86)×10⁻⁷ in the case of helioseismic data alone. The value deduced from all available data for the octopole is: J ₄=(2.8±2.1)×10⁻¹². These values are compared to some others found in the literature. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005238718479     

Accounting for the viscosity of the fluid core in the problem of the physical libration of the moon

A two-component theoretical model of the physical libration of the Moon in longitude is constructed with account taken of the viscosity of the core. In the new version, a hydrodynamic problem of motion of a fluid filling a solid rotating shell is solved. It is found that surfaces of equal angular velocity are spherical, and a velocity field of the fluid core of the Moon is described by elementary functions. A distribution of the internal pressure in the core is found. An angular momentum exchange between the fluid core and solid mantle is described by a third-order differential equation with a right-hand side. The roots of a characteristic equation are studied and the stability of rotation is proved. A libration angle as a function of time is found using the derived solution of the differential equation. Limiting cases of infinitely large and infinitely small viscosity are considered and an effect of lag of a libration phase from a phase of action of an external moment of forces is ascertained. This makes it possible to estimate the viscosity and sizes of the lunar fluid core from data of observations.     

On some problems concerning the calculation of the form of the magnetopause and the electric field on the magnetopause in the framework of the continuum medium model

In order to understand the reason of the existence of the electric field in the magnetosphere, and for the theoretical evaluation of its value, it is necessary to find the solution of the problem of determination of the magnetosphere boundary form in the frameworks of the continuum medium model which takes into account part of the magnetospheric plasma movement in supporting the magnetospheric boundary equilibrium. A number of problems for finding the distribution of the pressure, the density, the magnetic field and the electric field on the particular tangential discontinuity is considered in the case when the form of discontinuity is set (the direct problem) and a number of problems for finding the form of the discontinuity and the distribution of the above-mentioned physical quantities on the discontinuity is considered when the law of the change of the external pressure along the boundary is set (for example, with the help of the approximate Newton equation). The problem which is considered here, which deals with the calculation of the boundary form and with the calculation of the distribution of the corresponding physical quantities on the discontinuity of the 1st kind for the compressible fluid with the magnetic field with field lines which are perpendicular to the plane of the flow in question, concerns the last sort of problems. The comparison of the results of the calculation with the data in the equatorial cross-section of the magnetosphere demonstrates that the calculated form of the boundary, the value of the velocity of the return flow and the value of the electric field on the magnetopause, agree satisfactorily with the observational data.