Advances in numerical modeling of astrophysical and space plasmas

Plasma science is rich in distinguishable scales ranging from the atomic to the galactic to the meta-galactic, i.e., themesoscale. Thus plasma science has an important contribution to make in understanding the connection between microscopic and macroscopic phenomena. Plasma is a system composed of a large number of particles which interact primarily, but not exclusively, through the electromagnetic field. The problem of understanding the linkages and couplings in multi-scale processes is a frontier problem of modern science involving fields as diverse as plasma phenomena in the laboratory to galactic dynamics.Unlike the first three states of matter, plasma, often called the fourth state of matter, involves the mesoscale and its interdisciplinary founding have drawn upon various subfields of physics including engineering, astronomy, and chemistry. Basic plasma research is now posed to provide, with major developments in instrumentation and large-scale computational resources, fundamental insights into the properties of matter on scales ranging from the atomic to the galactic. In all cases, these are treated as mesoscale systems. Thus, basic plasma research, when applied to the study of astrophysical and space plasmas, recognizes that the behavior of the near-earth plasma environment may depend to some extent on the behavior of the stellar plasma, that may in turn be governed by galactic plasmas. However, unlike laboratory plasmas, astrophysical plasmas will forever be inaccessible to in situ observation. The inability to test concepts and theories of large-scale plasmas leaves only virtual testing as a means to understand the universe. Advances in in computer technology and the capability of performing physics first principles, fully three-dimensional, particle-in-cell simulations, are making virtual testing a viable alternative to verify our predictions about the far universe.The first part of this paper explores the dynamical and fluid properties of the plasma state, plasma kinetics, and the radiation emitted from plasmas. The second part of this paper outlines the formulation for the particle-in-cell simulation of astrophysical plasmas and advances in simulational techniques and algorithms, as-well-as the advances that may be expected as the computational resource grows to petaflop speed/memory capabilities.Dedicated to the memories of Hannes Alfvén and Oscar Buneman; Founders of the Subject.     

Spectral measurements of galactic and extragalactic sources in the hard X-ray region of 20-200 keV using balloon borne telescope

The spectral and temporal measurements in the hard X-ray region between 20-200 keV not only determines the extended behaviour of thermal X-ray spectrum below 10 keV but also provide a unique insight into the non-thermal processes in relativistic astrophysical plasma. From our present understanding of the X-ray sources, a significant fluxin the 20-200 keV band is expected from a variety of astrophysical phenomena, however, the available spectral data on the galactic and extragalactic X-ray source is very limited. This is mainly due to the fact that sensitivity of the detector systems used for earlier measurements was relatively poor. Since 1997, we have been carrying out a programme of hard X-ray observations galactic and extragalactic sources, in the 20-200 keV energy band using a highly sensitive balloon borne experiment. The X-ray telescope consists of three modules of large area scintillation counters specially configured in the back-to-back geometry and have a combined sensitivity of ∼ 10^-6 ph cm^-2 s^-1 keV^-1 for an on-source observations of 3 hrs. A total of 30 hours of ceiling data above an altitude of 3 mbar has been collected in 4 successful balloon flights from Hyderabad, India. Almost a dozen galactic and extragalactic X-ray sources were targeted and tracked during these observations. A positive detection was made in each case and in some cases the observed spectra extended right up to 150 keV. A brief account of the observed spectral and temporal features on some of the sources along with accurate measurement of diffuse background spectrum and a weak gamma ray burst will be presented in the paper. This revised version was published online in July 2006 with corrections to the Cover Date.     

The space environment of Mercury at the times of the second and third MESSENGER flybys

The second and third flybys of Mercury by the MESSENGER spacecraft occurred, respectively, on 6 October 2008 and on 29 September 2009. In order to provide contextual information about the solar wind properties and the interplanetary magnetic field (IMF) near the planet at those times, we have used an empirical modeling technique combined with a numerical physics-based solar wind model. The Wang–Sheeley–Arge (WSA) method uses solar photospheric magnetic field observations (from Earth-based instruments) in order to estimate the inner heliospheric radial flow speed and radial magnetic field out to 21.5 solar radii from the Sun. This information is then used as input to the global numerical magnetohydrodynamic model, ENLIL, which calculates solar wind velocity, density, temperature, and magnetic field strength and polarity throughout the inner heliosphere. WSA-ENLIL calculations are presented for the several-week period encompassing the second and third flybys. This information, in conjunction with available MESSENGER data, aid in understanding the Mercury flyby observations and provide a basis for global magnetospheric modeling. We find that during both flybys, the solar wind conditions were very quiescent and would have provided only modest dynamic driving forces for Mercury's magnetospheric system.     

Evolution of current loops in space

Ampere's law requires that every magnetic field have an associated current. The analysis of magnetic fields in this paper begins with that current in a circular loop and calculates the forces that make the loop evolve. A circular current generates a dipole field; and a second-order, ordinary differential equation represents the evolving magetic field. The theory describes cases where the conductor shrinks as the loop increases in size. The temperature of the conducting ions and electrons then decreases. The theory also describes cases where the conductor grows as the loop grows. Then the conducting particles heat up.Analysis shows that the magnetic clouds in the solar wind belong to the first type. In the provisional model adopted, the Klein-Burlaga clouds at one astronomical unit have a toroidal shape, centered on the sun, with a conductor radius of .125 au, and temperature (same for conducting electrons and protons) of 10⁵ K. After 26 days the toroid has a radius of 7.1 au, the conductor radius is .025 au, and the temperature is 2600 K.     

分形和混沌理论在太阳物理学中的应用

本文对非线性科学的两个重要分支－分形和混沌－在太阳物理学中的应用情况作了综述，主要内容包括：太阳活动混沌性的揭示；对太阳活动混沌性的可能解释－太阳非线性发电机理论；一些太阳现象的分形描述；耀斑的自组织临界行为研究。最后给出了作者对这一领域工作前景的展望。     

Ice lines, planetesimal composition and solid surface density in the solar nebula

To date, there is no core accretion simulation that can successfully account for the formation of Uranus or Neptune within the observed 2–3 Myr lifetimes of protoplanetary disks. Since solid accretion rate is directly proportional to the available planetesimal surface density, one way to speed up planet formation is to take a full accounting of all the planetesimal-forming solids present in the solar nebula. By combining a viscously evolving protostellar disk with a kinetic model of ice formation, which includes not just water but methane, ammonia, CO and 54 minor ices, we calculate the solid surface density of a possible giant planet-forming solar nebula as a function of heliocentric distance and time. Our results can be used to provide the starting planetesimal surface density and evolving solar nebula conditions for core accretion simulations, or to predict the composition of planetesimals as a function of radius. We find three effects that favor giant planet formation by the core accretion mechanism: (1) a decretion flow that brings mass from the inner solar nebula to the giant planet-forming region, (2) the fact that the ammonia and water ice lines should coincide, according to recent lab results from Collings et al. [Collings, M.P., Anderson, M.A., Chen, R., Dever, J.W., Viti, S., Williams, D.A., McCoustra, M.R.S., 2004. Mon. Not. R. Astron. Soc. 354, 1133–1140], and (3) the presence of a substantial amount of methane ice in the trans-saturnian region. Our results show higher solid surface densities than assumed in the core accretion models of Pollack et al. [Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y., 1996. Icarus 124, 62–85] by a factor of 3–4 throughout the trans-saturnian region. We also discuss the location of ice lines and their movement through the solar nebula, and provide new constraints on the possible initial disk configurations from gravitational stability arguments.     

Search for Phobos and Deimos gas/dust tori using in situ observations from Mars Global Surveyor MAG/ER

More than 490 elliptical aerobraking and science phasing orbits made by Mars Global Surveyor (MGS) in 1997 and 1998 provide unprecedented coverage of the solar wind in the vicinity of the orbits of the martian moons Phobos and Deimos. We have performed a comprehensive survey of magnetic field perturbations in the solar wind to search for possible signatures of solar wind interaction with dust or gas escaping from the moons. A total of 1246 solar wind disturbance events were identified and their distribution was examined relative to Phobos, the Phobos orbit, and the Deimos orbit. We find that the spatial distribution of solar wind perturbations does not increase near or downstream of Phobos, Phobos’ orbit, or Deimos’ orbit, which would have been expected if there is significant outgassing or dust escape from the martian moons. Of the 1246 magnetic field perturbation events found in the MGS data set, 11 events were found within 2000 km of the Phobos orbit, while three events were found within 2000 km of the Deimos orbit. These events were analyzed in detail and found to likely have other causes than outgassing/dust escape from the martian moons. Thus we conclude that the amount of gas/dust escaping the martian moons is not significant enough to induce detectable magnetic field perturbations in the solar wind. In essence we have not found any clear evidence in the MGS magnetic field data for outgassing or dust escape from the martian moons.     

Particle acceleration and nonthermal radiation in space plasmas

Acceleration processes for fast particles in astrophysical and space plasmas are reviewed with emphasis on stochastic acceleration by MHD turbulence and on acceleration by shock waves. Radiation processes in astrophysical and space plasmas are reviewed with emphasis on plasma emission from the solar corona and electron cyclotron maser emission from planets and stars.     

Making the Earth: Combining dynamics and chemistry in the Solar System

No terrestrial planet formation simulation completed to date has considered the detailed chemical composition of the planets produced. While many have considered possible water contents and late veneer compositions, none have examined the bulk elemental abundances of the planets produced as an important check of formation models. Here we report on the first study of this type. Bulk elemental abundances based on disk equilibrium studies have been determined for the simulated terrestrial planets of O’Brien et al. [O’Brien, D.P., Morbidelli, A., Levison, H.F., 2006. Icarus 184, 39-58]. These abundances are in excellent agreement with observed planetary values, indicating that the models of O’Brien et al. [O’Brien, D.P., Morbidelli, A., Levison, H.F., 2006. Icarus 184, 39-58] are successfully producing planets comparable to those of the Solar System in terms of both their dynamical and chemical properties. Significant amounts of water are accreted in the present simulations, implying that the terrestrial planets form “wet” and do not need significant water delivery from other sources. Under the assumption of equilibrium controlled chemistry, the biogenic species N and C still need to be delivered to the Earth as they are not accreted in significant proportions during the formation process. Negligible solar photospheric pollution is produced by the planetary formation process. Assuming similar levels of pollution in other planetary systems, this in turn implies that the high metallicity trend observed in extrasolar planetary systems is in fact primordial.     

Formation of terrestrial planets in a dissipating gas disk with Jupiter and Saturn

We have performed N-body simulations on final accretion stage of terrestrial planets, including the eccentricity and inclination damping effect due to tidal interaction with a gas disk. We investigated the dependence on a depletion time scale of the disk, and the effect of secular perturbations by Jupiter and Saturn. In the final stage, terrestrial planets are formed through coagulation of protoplanets of about the size of Mars. They would collide and grow in a decaying gas disk. Kominami and Ida [Icarus 157 (2002) 43-56] showed that it is plausible that Earth-sized, low-eccentricity planets are formed in a mostly depleted gas disk. In this paper, we investigate the formation of planets in a decaying gas disk with various depletion time scales, assuming disk surface density of gas component decays exponentially with time scale of τ_gas. Fifteen protoplanets with are initially distributed in the terrestrial planet regions. We found that Earth-sized planets with low eccentricities are formed, independent of initial gas surface density, when the condition (τ_cross+τ_growth)/2?τ_gas?τ_cross is satisfied, where τ_cross is the time scale for initial protoplanets to start orbit crossing in a gas-free case and τ_growth is the time scale for Earth-sized planets to accrete during the orbit crossing stage. In the cases satisfying the above condition, the final masses and eccentricities of the largest planets are consistent with those of Earth and Venus. However, four or five protoplanets with the initial mass remain. In the final stage of terrestrial planetary formation, it is likely that Jupiter and Saturn have already been formed. When Jupiter and Saturn are included, their secular perturbations pump up eccentricities of protoplanets and tend to reduce the number of final planets in the terrestrial planet regions. However, we found that the reduction is not significant. The perturbations also shorten τ_cross. If the eccentricities of Jupiter and Saturn are comparable to or larger than present values (∼0.05), τ_cross become too short to satisfy the above condition. As a result, eccentricities of the planets cannot be damped to the observed value of Earth and Venus. Hence, for the formation of terrestrial planets, it is preferable that the secular perturbations from Jupiter and Saturn do not have significant effect upon the evolution. Such situation may be reproduced by Jupiter and Saturn not being fully grown, or their eccentricities being smaller than the present values during the terrestrial planets' formation. However, in such cases, we need some other mechanism to eliminate the problem that numerous Mars-sized planets remain uncollided.     

Spectroscopic diagnostics in the VUV for solar and stellar plasmas

Summary The VUV emission spectra from the solar atmosphere and stellar atmospheres have been intensively studied during the past 25 years with several major space programs. In this review we discuss the spectroscopic diagnostic techniques used to study astrophysical plasmas, the atomic processes involved, the recent observations and the plans for future space missions.     

The thermal structure and the location of the snow line in the protosolar nebula: Axisymmetric models with full 3-D radiative transfer

The precise location of the water ice condensation front (‘snow line’) in the protosolar nebula has been a debate for a long time. Its importance stems from the expected substantial jump in the abundance of solids beyond the snow line, which is conducive to planet formation, and from the higher ‘stickiness’ in collisions of ice-coated dust grains, which may help the process of coagulation of dust and the formation of planetesimals. In an optically thin nebula, the location of the snow line is easily calculated to be around 3 AU, subject to brightness variations of the young Sun. However, in its first 5-10 myr, the solar nebula was optically thick, implying a smaller snowline radius due to shielding from direct sunlight, but also a larger radius because of viscous heating. Several models have attempted to treat these opposing effects. However, until recently treatments beyond an approximate 1 + 1D radiative transfer were unfeasible. We revisit the problem with a fully self-consistent 3D treatment in an axisymmetric disk model, including a density-dependent treatment of the dust and ice sublimation. We find that the location of the snow line is very sensitive to the opacities of the dust grains and the mass accretion rate of the disk. We show that previous approximate treatments are quite efficient at determining the location of the snow line if the energy budget is locally dominated by viscous accretion. Using this result we derive an analytic estimate of the location of the snow line that compares very well with results from this and previous studies. Using solar abundances of the elements we compute the abundance of dust and ice and find that the expected jump in solid surface density at the snow line is smaller than previously assumed. We further show that in the inner few AU the refractory species are also partly evaporated, leading to a significantly smaller solid state surface density in the regions where the rocky planets were formed.     

MHD compressive turbulence in the solar wind and the nearly incompressible approach

At MHD scales density fluctuation in the solar wind generally have a relative amplitude less than 0.1. The nearly incompressible MHD theory would seem then appropriate to describe a major part of the compressive turbulence at these scales. As a test of the theory, we focus on the scaling properties of density fluctuations with turbulent Mach numbers and on the level of correlation between density and temperature fluctuations. Our findings do not appear in favour of an extended applicability of the nearly incompressible theory to MHD compressive turbulence in the solar wind.     

Stochastic forces in space plasmas with ion-acoustic and lower-hybrid-drift turbulence

A Langevin equation for charged particles in a plasma with electrostatic turbulence is developed from first principles and in consistency with the kinetic theory in polarization approximation. For the case of ion-acoustic and electrostatic lower-hybrid-drift turbulence approximate expressions for the space-time spectral density of the wave energy are given and estimates of the intensities of the stochastic wave forces are made. The application is done for the plasmas of the earth's magnetosphere, the solar wind and solar flares. It seems, that ion-acoustic and electrostatic lower-hybrid-drift waves can contribute to electron chaotization in different regions of the space plasma.     

The dynamics of Jupiter and Saturn in the gaseous protoplanetary disk

We study the possibility that the mutual interactions between Jupiter and Saturn prevented Type II migration from driving these planets much closer to the Sun. Our work extends previous results by Masset and Snellgrove [Masset, F., Snellgrove, M., 2001. Mon. Not. R. Astron. Soc. 320, L55-L59], by exploring a wider set of initial conditions and disk parameters, and by using a new hydrodynamical code that properly describes for the global viscous evolution of the disk. Initially both planets migrate towards the Sun, and Saturn's migration tends to be faster. As a consequence, they eventually end up locked in a mean motion resonance. If this happens in the 2:3 resonance, the resonant motion is particularly stable, and the gaps opened by the planets in the disk may overlap. This causes a drastic change in the torque balance for the two planets, which substantially slows down the planets' inward migration. If the gap overlap is substantial, planet migration may even be stopped or reversed. As the widths of the gaps depend on disk viscosity and scale height, this mechanism is particularly efficient in low viscosity, cool disks. The initial locking of the planets in the 2:3 resonance is a likely outcome if Saturn formed at the edge of Jupiter's gap, but also if Saturn initially migrated rapidly from further away. We also explore the possibility of trapping in other resonances, and the subsequent evolutions. We discuss the compatibility of our results with the initial conditions adopted in Tsiganis et al. [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461] and Gomes et al. [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469] to explain the current orbital architecture of the giant planets and the origin of the Late Heavy Bombardment of the Moon.     

Ion flow and momentum transfer in the Venus plasma environment

An analysis of ion data from 390 Venus Express, VEX, orbits demonstrates that the flow of solar wind- and ionospheric ions near Venus is characterized by a marked asymmetry. The flow asymmetry of solar wind H⁺ and ionospheric O⁺ points steadily in the opposite direction to the planet’s orbital motion, and is most pronounced near the Pole and in the tail/nightside region. The flow asymmetry is consistent with aberration forcing, here defined as lateral forcing induced by the planet’s orbital motion. In addition to solar wind forcing by the radial solar wind expansion, Venus is also subject a lateral/aberration forcing induced by the planet’s orbital motion transverse to the solar wind flow.The ionospheric response to lateral solar wind forcing is analyzed from altitude profiles of the ion density, ion velocity and ion mass-flux. The close connection between decreasing solar wind H⁺ mass-flux and increasing ionospheric O⁺ mass-flux, is suggestive of a direct/local solar wind energy and momentum transfer to ionospheric plasma. The bulk O⁺ ion flow is accelerated to velocities less than 10 km/s inside the dayside/flank Ionopause, and up to 6000 km in the tail. Consequently, the bulk O⁺ outflow does not escape, but remains near Venus as a fast (km/s) O⁺ zonal wind in the Venus polar and nightside upper ionosphere. Furthermore, the total O⁺ mass-flux in the Venus induced magnetosphere, increases steadily downward to a maximum of 2 × 10⁻¹⁴ kg/(m² s) at ≈400 km altitude, suggesting a downward transport of energy and momentum. The O⁺, and total mass-flux, decay rapidly below 400 km. With no other plasma mass-flux as replacement, we argue that the reduction of ion mass-flux is caused by ion-neutral drag, a transfer of ion energy and momentum to neutrals, implying that the O⁺ plasma wind is converted to a neutral (thermosphere) wind at Venus. Incidentally, such a neutral wind would go in the same direction as the Venus atmosphere superrotation.     

The formation of Uranus and Neptune in solid-rich feeding zones: Connecting chemistry and dynamics

The core accretion theory of planet formation has at least two fundamental problems explaining the origins of Uranus and Neptune: (1) dynamical times in the trans-saturnian solar nebula are so long that core growth can take >15 Myr and (2) the onset of runaway gas accretion that begins when cores reach ∼10M_⊕ necessitates a sudden gas accretion cutoff just as Uranus and Neptune’s cores reach critical mass. Both problems may be resolved by allowing the ice giants to migrate outward after their formation in solid-rich feeding zones with planetesimal surface densities well above the minimum-mass solar nebula. We present new simulations of the formation of Uranus and Neptune in the solid-rich disk of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) using the initial semimajor axis distribution of the Nice model (Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A. [2005]. Nature 435, 466-469; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R. [2005]. Nature 435, 462-465; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F. [2005]. Nature 435, 459-461), with one ice giant forming at 12 AU and the other at 15 AU. The innermost ice giant reaches its present mass after 3.8-4.0 Myr and the outermost after 5.3-6 Myr, a considerable time decrease from previous one-dimensional simulations (e.g. Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y. [1996]. Icarus 124, 62-85). The core masses stay subcritical, eliminating the need for a sudden gas accretion cutoff.Our calculated carbon mass fractions of 22% are in excellent agreement with the ice giant interior models of Podolak et al. (Podolak, M., Weizman, A., Marley, M. [1995]. Planet. Space Sci. 43, 1517-1522) and Marley et al. (Marley, M.S., Gómez, P., Podolak, M. [1995]. J. Geophys. Res. 100, 23349-23354). Based on the requirement that the ice giant-forming planetesimals contain >10% mass fractions of methane ice, we can reject any Solar System formation model that initially places Uranus and Neptune inside of Saturn’s orbit. We also demonstrate that a large population of planetesimals must be present in both ice giant feeding zones throughout the lifetime of the gaseous nebula. This research marks a substantial step forward in connecting both the dynamical and chemical aspects of planet formation. Although we cannot say that the solid-rich solar nebula model of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) gives exactly the appropriate initial conditions for planet formation, rigorous chemical and dynamical tests have at least revealed it to be a viable model of the early Solar System.     

Use of multi-point analysis and modelling to address cross-scale coupling in space plasmas: Lessons from Cluster

The properties of plasmas (in space) are fundamentally governed by both ‘cross-scale’ coupling and comparative temporal behaviour operating over the micro-, meso-, and (MHD-) fluid regimes: for example, under conditions of turbulence, during magnetic reconnection and in shocks and other plasma boundaries. These themes map to a number of related and overlapping, phenomena, where known phenomena play different roles in each theme. Detailed understanding of fundamental plasma processes therefore requires analysis of both theoretical models (to distinguish the collisionless from the collisional regimes) and multi-scale measurements (suitable to address issues of stationarity). In particular, the investigation of phenomena requires analysis techniques which can distinguish and quantify temporal behaviour and the multi-scale spatial behaviour. The analysis of existing, multi-point data sets has led to a number of data co-ordination methods, such as the four spacecraft analysis tools developed for cluster, and we consider examples here. Advanced analysis concepts may be investigated with suitable considerations of measurement quality:adequate sampling of phenomena (for example, to extract the necessary information on the mechanisms operating) requires suitable spacecraft configurations and directly relates to the measurement quality achievable. A particular issue is how to resolve temporal behaviour across the spatial regimes, so that the data set is suitably coordinated. With the addition of theoretical modelling (in the context of particular phenomena) both the space and laboratory plasma regimes may be compared and we give an example of nonlinear wave coupling across spatial scales in this context.