Hydrodynamic instability of the solar nebula in the presence of a planetary core

When a planetary core composed of condensed matter is accumulated in the primitive solar nebula, the gas of the nebula becomes gravitationally concentrated as an envelope surrounding the planetary core. Models of such gaseous envelopes have been constructed subject to the assumption that the gas everywhere is on the same adiabat as that in the surrounding nebula. The gaseous envelope extends from the surface of the core to the distance at which the gravitational attraction of core plus envelope becomes equal to the gradient of the gravitational potential in the solar nebula; at this point the pressure and temperature of the gas in the envelope are required to attain the background values characteristics of the solar nebula. In general, as the mass of the condensed core increases, increasing amounts of gas became concentrated in the envelope, and these envelopes are stable against hydrodynamic instabilities. However, the core mass then goes through a maximum and starts to decrease. In most of the models tested, the envelopes were hydrodynamically unstable beyond the peak in the core mass. An unstable situation was always created if it was insisted that the core mass contain a larger amount of matter than given by these solutions. For an initial adiabat characterized by a temperature of 450°K and a pressure of 5 × 10^?6 atm, the maximum core mass at which instability occurs is approximately 115 earth masses; this value is rather insensitive to the position in the solar nebula or to the background pressure of the solar nebula. However, if the adiabat is lowered, then the core mass corresponding to instability is decreased. Since the core masses found by Podolak and Cameron for the giant planets are significantly less than the critical core mass corresponding to the initial solar nebula adiabat, we conclude that the giant planets obtained their large amounts of hydrogen and helium by a hydrodynamic collapse process in the solar nebula only after the nebula had been subjected to a considerable period of cooling.     

Propagation effects in the hydrogen-to-helium ratio in the solar cosmic rays

Taking as basis my dimensional solution of the propagation equation of solar cosmic rays in interplanetary space and using the equivalent diffusion coefficient found from observed data on solar protons, I discuss the effects of propagation on the hydrogen-to-hellum ratio Including its variation with the speed of the solar wind and with location in space. After eliminating the variations with solar distance and energy from the HEOS and PIONEER data collected by Perron et al., this ratio will be seen to increase with the magnetic longitude of the parent flare. The Initial value of the hydrogen-to-helium ratio at the Sun's surface, obtained after applying the corrections for propagation, is close to the value in the solar wind.     

Volatile element chemistry in the solar nebula: Na,K, F,Cl, Br,and P

The results of the most extensive set to date of thermodynamic calculations of the equilibrium chemistry of several hundred compounds of the elements Na, K, F, Cl, Br, and P in a solar composition system are reported. The calculations are carried out over a wide range of temperatures and pressures and along an adiabat in the primitive solar nebula. Two extreme models of accretion are investigated. In one extreme complete chemical equilibrium between condensates and gases is maintained because the time scale for accretion is long compared to the time scale for cooling or dissipation of the nebula. Condensates formed in this homogeneous accretion model include several phases such as whitlockite, alkali feldspars, and apatite minerals which are found in chondrites. In the other extreme complete isolation of newly formed condensates from prior condensates and gases occurs due to a time scale for accretion that is short relative to the time required for nebular cooling or dissipation. The condensates produced in this heterogeneous accretion model include alkali sulfides, ammonium halides, and ammonium phosphates. None of these phases are found in chondrites. Available observations of the Na, K, F, Cl, Br, and P elemental abundances in the terrestrial planets are found to be compatible with the predictions of the homogeneous accretion model.     

Helioseismic Constraints on Rotation and Magnetic Fields in the Solar Core

In recent years, more and more precise measurements have been made of solar oscillation frequencies and line widths. From space, the Solar and Heliospheric Observatory/Michelson Doppler Imager (MDI) data has led to much progress. From the ground, networks, like Global Oscillation Network Group (GONG), Taiwanese Oscillation Network (TON), and Birmingham Solar Oscillations Network (BiSON) have also led to much progress. The sharpened and enriched oscillation spectrum of data have been critically complemented by advances in the treatments of the opacities and the equation of state. All of this has led to a significantly more precise probing of the solar core. Here we discuss the progress made and suggest how the core may be better probed with seismic data on-hand. In particular, we review our knowledge of the rotation and structure of the core. We further argue that much may be learned about the core by exploiting the line width data from the aforementioned sources. Line-width data can be used to place sharper constraints on core properties, like the degree to which the Sun rotates on a single axis and the upper limit on magnetic fields that may be buried in the core.     

Constrained estimates of low-degree mode frequencies and the determination of the interior structure of the Sun

Low-degreep-modes penetrate to the solar centre and provide direct information about the core. However, the high observational accuracy that is required to resolve the details of the structure of the core is difficult to achieve because the oscillation power spectrum is significantly distorted by stochastic forcing of the oscillations, which appears as multiplicative noise. Here, an attempt is reported to reduce uncertainties of spectral parameter estimation by incorporating constraints imposed by smooth behaviour of some of the parameters (e.g., linewidths, background noise, rotational splitting) over a group of lines. Instead of estimating these parameters independently for each line, we determine them as smooth functions of frequency. It is expected that this procedure gives more accurate estimates of the average frequencies of any multiplet in the power spectrum, to which we have applied no constraints. We give some examples of the procedure for whole-disk measurements by the IPHIR space experiment. It is shown that the additional constraints do not result in significant changes in the frequency estimates, except for one mode whose peak in the power spectrum has the lowest signal-to-noise ratio. However, the uncertainty in the frequency of that mode does not influence substantially the results of the structure inversion in the core. Inversions of the IPHIR datasets are compared with corresponding inversions of data from the Birmingham Solar Oscillation Network (BISON). The IPHIR data indicate a sharp increase towards the centre of the deviation of the squared sound speed of the sun from that of a standard solar model, whereas the BISON data show a decrease. The difference between the IPHIR and BISON inversions is significant, preventing any definite conclusion about the deviation of the structure of the solar core from that of the model.     

A mixed solar core, solar neutrinos and helioseismology

We consider a wide class of solar models with mixed core. Most of these models can be excluded as the sound speed profile that they predict is in sharp disagreement with helioseismic constraints. All the remaining models predict ⁸B and/or ⁷Be neutrino fluxes of at least as large as those of SSMs. In conclusion, helioseismology shows that a mixed solar core cannot account for the neutrino deficit implied by the current solar neutrino experiments.     

Source regions and timescales for the delivery of water to the Earth

Abstract— In the primordial solar system, the most plausible sources of the water accreted by the Earth were in the outer asteroid belt, in the giant planet regions, and in the Kuiper Belt. We investigate the implications on the origin of Earth's water of dynamical models of primordial evolution of solar system bodies and check them with respect to chemical constraints. We find that it is plausible that the Earth accreted water all along its formation, from the early phases when the solar nebula was still present to the late stages of gas‐free sweepup of scattered planetesimals. Asteroids and the comets from the Jupiter‐Saturn region were the first water deliverers, when the Earth was less than half its present mass. The bulk of the water presently on Earth was carried by a few planetary embryos, originally formed in the outer asteroid belt and accreted by the Earth at the final stage of its formation. Finally, a late veneer, accounting for at most 10% of the present water mass, occurred due to comets from the Uranus‐Neptune region and from the Kuiper Belt. The net result of accretion from these several reservoirs is that the water on Earth had essentially the D/H ratio typical of the water condensed in the outer asteroid belt. This is in agreement with the observation that the D/H ratio in the oceans is very close to the mean value of the D/H ratio of the water inclusions in carbonaceous chondrites.     

Methods for computing giant planet formation and evolution

We present a numerical code for computing all stages of the formation and evolution of giant planets in the framework of the core instability mechanism. This code is a non-trivial adaption of the stellar binary evolution code and is based on a standard Henyey technique. To investigate the performance of this code we applied it to the computation of the formation and evolution of a Jupiter mass object from a half Earth core mass to ages in excess of the age of the Universe.
We also present a new smoothed linear interpolation algorithm devised especially for the purpose of circumventing some problems found when some physical data (e.g. opacities, equation of state, etc.) are introduced into an implicit algorithm like the one employed in this work.     

Models of the giant planets

Models of the giant planets were constructed based on the assumption that the hydrogen to helium ratio is solar in these planets. This assumption, together with arguments about the condensation sequence in the primitive solar nebula, yields models with a central core of rock and possibly ice surrounded by an envelope of hydrogen, helium, methane, ammonia, and water. These last three volatiles may be individually enhanced due to condensation at the period of core formation. Jupiter was found to have a core of about 40 earth masses and a water enhancement in the atmosphere of about 7.5 times the solar value. Saturn was found to have a core of 20 earth masses and a water enhancement in the atmosphere of about 25 times the solar value. Rock plus ice constitute 75–85% of the mass of Uranus and Neptune. Temperatures in the interiors of these planets are probably above the melting points, if there is an adiabatic relation throughout the interiors. Some aspects of the sensitivities of these results to uncertainties in rotational flattening are discussed.     

Helioseismology and the solar interior dynamics

The inversion of helioseismic modes leads to the sound velocity inside the Sun with a precision of about 0.1 per cent. Comparisoons of solar models with the “seismic sun” represent powerful tools to test the physics: depth of the convection zone, equation of state, opacities, element diffusion processes and mixing inside the radiative zone. We now have evidence that microscopic diffusion (element segregation) does occur below the convection zone, leading to a mild helium depletion in the solar outer layers. Meanwhile this process must be slowed down by some macroscopic effect, presumably rotation-induced mixing. The same mixing is also responsible for the observed lithium depletion. On the other hand, the observations of beryllium and helium 3 impose specific constraints on the depth of this mildly mixed zone. Helioseismology also gives information on the internal solar rotation: while differential rotation exists in the convection zone, solid rotation prevails in the radiative zone, and the transition layer (the so-called “tachocline”) is very small. These effects are discussed, together with the astrophysical constraints on the solar neutrino fluxes.     

Solar wind isotopic abundance ratios of ne,mg and si measured by SOHO/CELIAS/MTOF as diagnostic tool for the inner corona

Using the high-resolution mass spectrometer MTOF on board SOHO we have measured the solar wind isotopic abundance ratios of Ne, Mg, and Si in different solar wind regimes with bulk velocities ranging from 350 km/s to 650 km/s. Data indicate a systematic depletion of the heavier isotopes in the slow solar wind compared to their abundances in the fast solar wind from coronal holes. These variations in the solar wind isotopic composition represent a pure mass-dependent effect because the different isotopes of an element pass the inner corona with the same charge state distribution. The influence of particle mass on the acceleration of minor solar wind ions is discussed in the context of theoretical models and recent optical observations with other SOHO instruments.     

A comparison of the interiors of Jupiter and Saturn

Interior models of Jupiter and Saturn are calculated and compared in the framework of the three-layer assumption, which rely on the perception that both planets consist of three globally homogeneous regions: a dense core, a metallic hydrogen envelope, and a molecular hydrogen envelope. Within this framework, constraints on the core mass and abundance of heavy elements (i.e. elements other than hydrogen and helium) are given by accounting for uncertainties on the measured gravitational moments, surface temperature, surface helium abundance, and on the inferred protosolar helium abundance, equations of state, temperature profile and solid/differential interior rotation. Results obtained solely from static models matching the measured gravitational fields indicate that the mass of Jupiter’s dense core is less than 14 M_⊕ (Earth masses), but that models with no core are possible given the current uncertainties on the hydrogen–helium equation of state. Similarly, Saturn’s core mass is less than 22 M_⊕ but no lower limit can be inferred. The total mass of heavy elements (including that in the core) is constrained to lie between 11 and 42 M_⊕ in Jupiter, and between 19 and 31 M_⊕ in Saturn. The enrichment in heavy elements of their molecular envelopes is 1–6.5, and 0.5–12 times the solar value, respectively. Additional constraints from evolution models accounting for the progressive differentiation of helium (Hubbard WB, Guillot T, Marley MS, Burrows A, Lunine JI, Saumon D, 1999. Comparative evolution of Jupiter and Saturn. Planet. Space Sci. 47, 1175–1182) are used to obtain tighter, albeit less robust, constraints. The resulting core masses are then expected to be in the range 0–10 M_⊕, and 6–17 M_⊕ for Jupiter and Saturn, respectively. Furthermore, it is shown that Saturn’s atmospheric helium mass mixing ratio, as derived from Voyager, Y=0.06±0.05, is probably too low. Static and evolution models favor a value of Y=0.11−0.25. Using, Y=0.16±0.05, Saturn’s molecular region is found to be enriched in heavy elements by 3.5 to 10 times the solar value, in relatively good agreement with the measured methane abundance. Finally, in all cases, the gravitational moment J₆ of models matching all the constraints are found to lie between 0.35 and 0.38×10⁻⁴ for Jupiter, and between 0.90 and 0.98×10⁻⁴ for Saturn, assuming solid rotation. For comparison, the uncertainties on the measured J₆ are about 10 times larger. More accurate measurements of J₆ (as expected from the Cassini orbiter for Saturn) will therefore permit to test the validity of interior models calculations and the magnitude of differential rotation in the planetary interior.     

Are Non-Magnetic Mechanisms Such As Temporal Solar Diameter Variations Conceivable for an Irradiance Variability?

Irradiance variability has been monitored from space for more than two decades. Even though data come from different sources, it is well established that a temporal variability exists ≈0.1%, in phase with the solar cycle. Today, one of the best explanations for such an irradiance variability is provided by the evolution of the solar surface magnetic fields. But if some 90–95% can be reproduced, what would be the origin of the 10–5% left? Non-magnetic effects are conceivable. In this paper we will consider temporal variations of the diameter of the Sun as a possible contributor for the remaining part. Such an approach imposes strong constraints on the solar radius variability. We will show that over a solar cycle, variations of no more than 20 mas of amplitude can be considered. Such a variability – far from what is reported by observers conducting measurements by means of ground-based solar astrolabes – may explain a little part of the irradiance changes not explained by magnetic features. Further, requirements are needed that may help to reach a conclusion. Dedicated space missions are necessary (for example PICARD, GOLF-NG or SDO, scheduled for a launch around 2008); it is also proposed to reactivate SDS flights for such a purpose.     

Influence of Low-Degree High-Order p-Mode Splittings on the Solar Rotation Profile

The solar rotation profile is well constrained down to about 0.25R _⊙ thanks to the study of acoustic modes. Since the radius of the inner turning point of a resonant acoustic mode is inversely proportional to the ratio of its frequency to its degree, only the low-degree p modes reach the core. The higher the order of these modes, the deeper they penetrate into the Sun and thus they carry more diagnostic information on the inner regions. Unfortunately, the estimates of frequency splittings at high frequency from Sun-as-a-star measurements have higher observational errors because of mode blending, resulting in weaker constraints on the rotation profile in the inner core. Therefore inversions for the solar internal rotation use only modes below 2.4 mHz for ?≤3. In the work presented here, we used an 11.5-year-long time series to compute the rotational frequency splittings for modes ?≤3 using velocities measured with the GOLF instrument. We carried out a theoretical study of the influence of the low-degree modes in the region from 2 to 3.5 mHz on the inferred rotation profile as a function of their error bars.     

Effects of non-Newtonian gravity on the properties of strange stars

The properties of strange star matter are studied in the equivparticle model with inclusion of non-Newtonian gravity. It is found that the inclusion of non-Newtonian gravity makes the equation of state stiffer if Witten's conjecture is true. Correspondingly, the maximum mass of strange stars becomes as large as two times the solar mass, and the maximum radius also becomes bigger. The coupling to boson mass ratio has been constrained within the stability range of strange quark matter.     

On the collapse of neutron stars and stellar cores to pion-condensed stars

The transition from a neutron star to a pion-condensed star is investigated in Newtonian hydrodynamics. It is shown that in a certain range of ultradense equations of state, there occurs a mass ejection with energies comparable with usual supernova outputs. But the ejected mass is only in the order of 0.02M _⊙. Therefore, the observable consequences of this transition are not so dramatic as conjectured recently. In a realistic scenario including a stiff ultradense equation of state and a weak effect of pion condensation the mass ejection disappears. Additionally the collapse of a stellar core to a neutron star with pion-condensed core is considered. In comparison with a standard supernova scenario we find only a slightly reduced explosion energy. Further, the possible consequence of pion condensation during the secular evolution of the bounced core of a collapsing star to the cool final neutron star is discussed.     

Piecewise mass flows within a solar prominence observed by the New Vacuum Solar Telescope

The material of solar prominences is often observed in a state of flowing. These mass flows (MF) are important and useful for us to understand the internal structure and dynamics of prominences. In this paper, we present a high resolution H\(\alpha \) observation of MFs within a quiescent solar prominence. From the observation, we find that the plasma primarily has a circular motion and a downward motion separately in the middle section and legs of the prominence, which creates a piecewise mass flow along the observed prominence. Moreover, the observation also shows a clear displacement of MF’s velocity peaks in the middle section of the prominence. All of these provide us with a detailed record of MFs within a solar prominence and show a new approach to detecting the physical properties of prominence.     

Short‐lived radioactivity in the early solar system: The Super‐AGB star hypothesis

Abstract– The composition of the most primitive solar system condensates, such as calcium‐aluminum‐rich inclusions (CAIs) and micron‐sized corundum grains, show that short‐lived radionuclides (SLR), e.g., ²⁶Al, were present in the early solar system. Their abundances require a local or stellar origin, which, however, is far from being understood. We present for the first time the abundances of several SLR up to ⁶⁰Fe predicted from stars with initial mass in the range approximately 7–11 M_⊙. These stars evolve through core H, He, and C burning. After core C burning they go through a “Super”‐asymptotic giant branch (Super‐AGB) phase, with the H and He shells activated alternately, episodic thermal pulses in the He shell, a very hot temperature at the base of the convective envelope (approximately 10⁸ K), and strong stellar winds driving the H‐rich envelope into the surrounding interstellar medium. The final remnants of the evolution of Super‐AGB stars are mostly O–Ne white dwarfs. Our Super‐AGB models produce ²⁶Al/²⁷Al yield ratios approximately 0.02–0.26. These models can account for the canonical value of the ²⁶Al/²⁷Al ratio using dilutions with the solar nebula of the order of 1 part of Super‐AGB mass per several 10² to several 10³ of solar nebula mass, resulting in associated changes in the O‐isotope composition in the range Δ¹⁷O from 3 to 20‰. This is in agreement with observations of the O isotopic ratios in primitive solar system condensates, which do not carry the signature of a stellar polluter. The radionuclides ⁴¹Ca and ⁶⁰Fe are produced by neutron captures in Super‐AGB stars and their meteoritic abundances are also matched by some of our models, depending on the nuclear and stellar physics uncertainties as well as the meteoritic experimental data. We also expect and are currently investigating Super‐AGB production of SLR heavier than iron, such as ¹⁰⁷Pd.     

Geodesy constraints on the interior structure and composition of Mars

Knowledge of the interior structure of Mars is of fundamental importance to the understanding of its past and present state as well as its future evolution. The most prominent interior structure properties are the state of the core, solid or liquid, its radius, and its composition in terms of light elements, the thickness of the mantle, its composition, the presence of a lower mantle, and the density of the crust. In the absence of seismic sounding only geodesy data allow reliably constraining the deep interior of Mars. Those data are the mass, moment of inertia, and tides. They are related to Mars’ composition, to its internal mass distribution, and to its deformational response to principally the tidal forcing of the Sun. Here we use the most recent estimates of the moment of inertia and tidal Love number k₂ in order to infer knowledge about the interior structure of the Mars.We have built precise models of the interior structure of Mars that are parameterized by the crust density and thickness, the volume fractions of upper mantle mineral phases, the bulk mantle iron concentration, and the size and the sulfur concentration of the core. From the bulk mantle iron concentration and from the volume fractions of the upper mantle mineral phases, the depth dependent mineralogy is deduced by using experimentally determined phase diagrams. The thermoelastic properties at each depth inside the mantle are calculated by using equations of state. Since it is difficult to determine the temperature inside the mantle of Mars we here use two end-member temperature profiles that have been deduced from studies dedicated to the thermal evolution of Mars. We calculate the pressure and temperature dependent thermoelastic properties of the core constituents by using equations state and recent data about reference thermoelastic properties of liquid iron, liquid iron-sulfur, and solid iron. To determine the size of a possible inner core we use recent data on the melting temperature of iron-sulfur.Within our model assumptions the geodesy data imply that Mars has no solid inner core and that the liquid core contains a large fraction of sulfur. The absence of a solid inner is in agreement with the absence of a global magnetic field. We estimate the radius of the core to be 1794 ± 65 km and its core sulfur concentration to be 16 ± 2 wt%. We also show that it is possible for Mars to have a thin layer of perovskite at the bottom of the mantle if it has a hot mantle temperature. Moreover a chondritic Fe/Si ratio is shown to be consistent with the geodesy data, although significantly different value are also possible. Our results demonstrate that geodesy data alone, even if a mantle temperature is assumed, can almost not constrain the mineralogy of the mantle and the crust. In order to obtain stronger constraints on the mantle mineralogy bulk properties, like a fixed Fe/Si ratio, have to be assumed.     

Coupled orbit-attitude dynamics of high area-to-mass ratio (HAMR) objects: influence of solar radiation pressure, Earth’s shadow and the visibility in light curves

The orbital and attitude dynamics of uncontrolled Earth orbiting objects are perturbed by a variety of sources. In research, emphasis has been put on operational space vehicles. Operational satellites typically have a relatively compact shape, and hence, a low area-to-mass ratio (AMR), and are in most cases actively or passively attitude stabilized. This enables one to treat the orbit and attitude propagation as decoupled problems, and in many cases the attitude dynamics can be neglected completely. The situation is different for space debris objects, which are in an uncontrolled attitude state. Furthermore, the assumption that a steady-state attitude motion can be averaged over data reduction intervals may no longer be valid. Additionally, a subset of the debris objects have significantly high area-to-mass ratio (HAMR) values, resulting in highly perturbed orbits, e.g. by solar radiation pressure, even if a stable AMR value is assumed. Note, this assumption implies a steady-state attitude such that the average cross-sectional area exposed to the sun is close to constant. Time-varying solar radiation pressure accelerations due to attitude variations will result in un-modeled errors in the state propagation. This work investigates the evolution of the coupled attitude and orbit motion of HAMR objects. Standardized pieces of multilayer insulation (MLI) are simulated in a near geosynchronous orbits. It is assumed that the objects are rigid bodies and are in uncontrolled attitude states. The integrated effects of the Earth gravitational field and solar radiation pressure on the attitude motion are investigated. The light curves that represent the observed brightness variations over time in a specific viewing direction are extracted. A sensor model is utilized to generate light curves with visibility constraints and magnitude uncertainties as observed by a standard ground based telescope. The photometric models will be needed when combining photometric and astrometric observations for estimation of orbit and attitude dynamics of non-resolved space objects.