From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets

We present results from a suite of N-body simulations that follow the formation and accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5 and 4.0 AU and the collisional growth is followed until the completion of planetary accretion (>100 Myr). A total of 64 simulations were carried out to explore sensitivity to the key parameters and initial conditions. All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas (which dissipates uniformly in space and exponentially in time), and orbits of Jupiter and Saturn. We end up with 1-5 planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today’s eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts, as long as we include all the above effects of gas and assume initial disks similar to the minimum mass solar nebular.     

Alfven Wave Reflection model of field-aligned currents at Mercury

An Alfven Wave Reflection (AWR) model is proposed that provides closure for strong field-aligned currents (FACs) driven by the magnetopause reconnection in the magnetospheres of planets having no significant ionospheric and surface electrical conductance. The model is based on properties of the Alfven waves, generated at high altitudes and reflected from the low-conductivity surface of the planet. When magnetospheric convection is very slow, the incident and reflected Alfven waves propagate along approximately the same path. In this case, the net field-aligned currents will be small. However, as the convection speed increases, the reflected wave is displaced relatively to the incident wave so that the incident and reflected waves no longer compensate each other. In this case, the net field-aligned current may be large despite the lack of significant ionospheric and surface conductivity. Our estimate shows that for typical solar wind conditions at Mercury, the magnitude of Region 1-type FACs in Mercury’s magnetosphere may reach hundreds of kilo-Amperes. This AWR model of field-aligned currents may provide a solution to the long-standing problem of the closure of FACs in the Mercury’s magnetosphere.     

The invention that opened the solar system to exploration

The invention of gravity-propelled interplanetary space travel (also known as “gravity-assist trajectories”) in the early 1960s broke the high-energy barrier of classical space travel based on reaction propulsion, and made possible the exploration of the entire solar system with instrumented spacecraft. In this concept, a free-fall spacecraft is launched from a launch planet P₁ to a nearby planet P₂ such that its gravitational field (superimposed on the gravitational field of the Sun) catapults the vehicle to another planet P₃, which in turn is used to repeat the process. Thus, through a series of planetary encounters, a gravity-propelled trajectory P₁-P₂-P₃-P₄-…-P_N is generated. This paper describes how the invention was conceived and how the difficult mathematical problem of computing the trajectories was solved in order to numerically investigate and use the invention in actual missions. The crucial roles played by the UCLA Computing Facility and the Departments of Mathematics and Physics are also described.     

Scientific Objectives of China-Russia Joint Mars Exploration Program YH-1

Compared with other planets, Mars is a planet most similar with the earth and most possible to find the extraterrestrial life on it, and therefore especially concerned about by human beings. In recent years, some countries have launched Mars probes and announced their manned Mars exploration programs. China has become the fifth country in the world to launch independently artificial satellites, and the third country able to carry out an independent manned space program. However, China is just at the beginning of deep space explorations. In 2007, China and Russia signed an agreement on a joint Mars exploration program by sending a Chinese micro-satellite Yinghuo-1 (YH-1) to the Mars orbit. Once YH-1 enters its orbit, it will carry out its own exploration, as well as the joint exploration with the Russian Phobos-Grunt probe. This paper summarizes the scientific background and objectives of YH-1 and describes briefly its payloads for realizing these scientific objectives. In addition, the main exploration tasks of YH-1 and a preliminary prospect on its exploration results are also given.     

Resonance sticking in the scattered disk

We investigate the dynamical evolution of trans-neptunian objects (TNOs) in typical scattered disk orbits (scattered TNOs) by performing simulations using several thousand particles lying initially on Neptune-encountering orbits. We explore the role of resonance sticking in the scattered disk, a phenomenon characterized by multiple temporary resonance captures (‘resonances’ refers to external mean motion resonances with Neptune, which can be described in the form r:s, where the arguments r and s are integers). First, all scattered TNOs evolve through intermittent temporary resonance capture events and gravitational scattering by Neptune. Each scattered TNO experiences tens to hundreds of resonance captures over a period of 4 Gyr, which represents about 38% of the object's lifetime (mean value). Second, resonance sticking plays an important role at semimajor axes , where the great majority of such captures occurred. It is noteworthy that the stickiest (i.e., dominant) resonances in the scattered disk are located within this distance range and are those possessing the lowest argument s. This was evinced by r:1, r:2 and r:3 resonances, which played the greatest role during resonance sticking evolution, often leading to captures in several of their neighboring resonances. Finally, the timescales and likelihood of temporary resonance captures are roughly proportional to resonance strength. The dominance of low s resonances is also related to the latter. In sum, resonance sticking has an important impact on the evolution of scattered TNOs, contributing significantly to the longevity of these objects.     

Comparison of the Characteristics of Magnetic Clouds and Magnetic Cloud-Like Structures for the Events of 1995 – 2003

Using nine years of solar wind plasma and magnetic field data from the Wind mission, we investigated the characteristics of both magnetic clouds (MCs) and magnetic cloud-like structures (MCLs) during 1995 – 2003. A MCL structure is an event that is identified by an automatic scheme (Lepping, Wu, and Berdichevsky, Ann. Geophys. 23, 2687, 2005) with the same criteria as for a MC, but it is not usually identifiable as a flux rope by using the MC (Burlaga et al., J. Geophys. Res. 86, 6673, 1981) fitting model developed by Lepping, Jones, and Burlaga (Geophys. Res. Lett. 95(11), 957, 1990). The average occurrence rate is 9.5 for MCs and 13.6 for MCLs per year for the overall period of interest, and there were 82 MCs and 122 MCLs identified during this period. The characteristics of MCs and MCL structures are as follows: (1) The average duration, Δt, of MCs is 21.1 h, which is 40% longer than that for MCLs (Δt=15 h); (2) the average (minimum B _z found in MC/MCL measured in geocentric solar ecliptic coordinates) is −10.2 nT for MCs and −6 nT for MCLs; (3) the average Dst_min (minimum Dst caused by MCs/MCLs) is −82 nT for MCs and −37 nT for MCLs; (4) the average solar wind velocity is 453 km s⁻¹ for MCs and 413 km s⁻¹ for MCLs; (5) the average thermal speed is 24.6 km s⁻¹ for MCs and 27.7 km s⁻¹ for MCLs; (6) the average magnetic field intensity is 12.7 nT for MCs and 9.8 nT for MCLs; (7) the average solar wind density is 9.4 cm⁻³ for MCs and 6.3 cm⁻³ for MCLs; and (8) a MC is one of the most important interplanetary structures capable of causing severe geomagnetic storms. The longer duration, more intense magnetic field and higher solar wind speed of MCs, compared to those properties of the MCLs, are very likely the major reasons for MCs generally causing more severe geomagnetic storms than MCLs. But the fact that a MC is an important interplanetary structure with respect to geomagnetic storms is not new (e.g., Zhang and Burlaga, J. Geophys. Res. 93, 2511, 1988; Bothmer, ESA SP-535, 419, 2003).     

中国太阳能供热采暖技术的现状与发展

本文介绍了我国太阳能供热采暖技术的发展现状及工程实例,指出了该项技术与国际先进水平的差距及存在的主要问题,提出了太阳能供热采暖今后在我国的发展趋势、发展目标和重点任务。     

Measurement of the radial density gradient of cosmic ray in the heliosphere by the GRAPES-3 experiment

A radial anisotropy in the flux of cosmic rays in heliosphere was theoretically predicted by Parker and others within the framework of the diffusion–convection mechanism. The solar wind is responsible for sweeping out the galactic cosmic rays, creating a radial density gradient within the heliosphere. This gradient coupled with the interplanetary magnetic field induces a flow of charged particles perpendicular to the ecliptic plane which was measured and correctly explained by Swinson, and is hereafter referred as ‘Swinson flow’. The large area GRAPES-3 tracking muon telescope offers a powerful probe to measure the Swinson flow and the underlying radial density gradient of the galactic cosmic rays at a relatively high rigidity of ∼100 GV. The GRAPES-3 data collected over a period of six years (2000–2005) were analyzed and the amplitude of the Swinson flow was estimated to be (0.0644 ± 0.0008)% of cosmic ray flux which was an ∼80σ effect. The phase of the maximum flow was at a sidereal time of (17.70 ± 0.05) h which was 18 min earlier than the expected value of 18 h. This small 18 min phase difference had a significance of ∼6σ indicating the inherent precision of the GRAPES-3 measurement. The radial density gradient of the galactic cosmic rays at a median rigidity of 77 GV was found to be 0.65% AU⁻¹.     

无力磁场的稳定性

无力场被广泛用来模拟太阳活动区的强磁场,本文从Bernstein能量原理出发,导出了无力场能量原理的普遍形式,并给出了若干稳定性的充分条件,它们可方便而有效地对无力场进行稳定性判断。