An estimation of the oblateness of Sun from the motion of Icarus

From astrometric observations of minor planet (1566) Icarus from 1949 to 1987 were made solutions for improved orbital elements of Icarus and the quadrupole moment of the Sun. The formal result was J₂ = -0.6±5.8 &d 10^–6. From this we can conclude that J ₂ is very probably less than 2 · 10^–-5.     

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     

Solar quadrupole moment and purely relativistic gravitation contributions to Mercury's perihelion advance

The perihelion advance of the orbit of Mercury has long been one of the observational cornerstones for testing General Relativity (G.R.).The main goal of this paper is to discuss how, presently, observational and theoretical constraints may challenge Einstein's theory of gravitation characterized by β=γ=1. To achieve this purpose, we will first recall the experimental constraints upon the Eddington-Robertson parameters γ,β and the observational bounds for the perihelion advance of Mercury, Δω_obs. A second point will address the values given, up to now, to the solar quadrupole moment by several authors. Then, we will briefly comment why we use a recent theoretical determination of the solar quadrupole moment, J ₂=(2.0 ± 0.4) 10^-7, which takes into account both surfacic and internal differential rotation, in order to compute the solar contribution to Mercury's perihelion advance. Further on, combining bounds on γ and J ₂ contributions, and taking into account the observational data range for Δω_obs,we will be able to give a range of values for β. Alternatively, taking into account the observed value of Δω_obs, one can deduce a dynamical estimation of J ₂ in the setting of G.R. This point is important as it provides a solar model independent estimation that can be confronted with other determinations of J ₂ based upon solar theory and solar observations (oscillation data, oblateness...). Finally, a glimpse at future satellite experiments will help us to understand how stronger constraints upon the parameter space (γω J ₂) as well as a separation of the two contributions (from the quadrupole moment, J ₂, or purely relativistic, 2α²+2αγ–β) might be expected in the future. This revised version was published online in July 2006 with corrections to the Cover Date.     

Analytical Approach to the Motion of a Lunar Artificial Satellite

The canonical equations of motion of an artificial lunar satellite are formulated including the effects of the asphericity of the Moon comprising the harmonics J ₂, J ₂₂, J ₃, J ₃₁, J ₄ andJ ₅, the oblateness of the Earth up to the second zonal harmonic, as well as the disturbing function due to the attractions of the Earth and of the Sun (terms are retained up to order 10^-6 for the higher orbits and 10^-8 for the lower orbits). This revised version was published online in July 2006 with corrections to the Cover Date.     

Impact of the Quadrupole Moment of the Sun on the Dynamics of the Earth-Moon System

Range of values of the Sun's mass quadrupole moment of coefficient J₂ arising both from experimental and theoretical determinations enlarge across literature on two orders of magnitude, from around 10^-7 until to 10^-5. The accurate knowledge of the Moon's physical librations, for which the Lunar Laser Ranging data reach an outstanding precision level, prove to be appropriate to reduce the interval of J₂ values by giving an upper bound of J₂. A solar quadrupole moment as high as 1.1 10^-5 given either from the upper bounds of the error bars of the observations, or from the Roche's theory, is not compatible with the knowledge of the lunar librations accurately modeled and observed with the LLR experiment. The suitable values of J₂ have to be smaller than 3.0 10^-6. As a consequence, this upper bound of 3.0 10^-6 is accepted to study the impact of the Sun's quadrupole moment of mass on the dynamics of the Earth-Moon system. Such as effect (with J₂ = 5.5±1.3 × 10^-6) has been already tested in 1983 by Campbell & Moffat using analytical approximate equations, and thus for the orbits of Mercury, Venus, the Earth and Icarus. The approximate equations are no longer sufficient compared with present observational data and exact equations are required. As if to compute the effect on the lunar librations, we have used our BJV relativistic model of solar system integration including the spin-orbit coupled motion of the Moon. The model is solved by numerical integration. The BJV model stems from general relativity by using the DSX formalism for purposes of celestial mechanics when it is about to deal with a system of n extended, weakly self-gravitating, rotating and deformable bodies in mutual interactions. The resulting effects on the orbital elements of the Earth have been computed and plotted over 160 and 1600 years. The impact of the quadrupole moment of the Sun on the Earth's orbital motion is mainly characterized by variations of , , and . As a consequence, the Sun's quadrupole moment of mass could play a sensible role over long time periods of integration of solar system models. This revised version was published online in July 2006 with corrections to the Cover Date.     

The abundance of 3He in the solar wind - A constraint for models of solar evolution

³He is an intermediate product in the proton-proton chain, and standard models of the Sun predict a large bulge of enhanced ³He abundance near M _r /M ₀ = 0.6 in the contemporary Sun. The relatively low abundance of ³He at the solar surface, which is derived from solar wind observations, poses severe constraints to non-standard solar models.Direct measurements of the ³He abundance in the solar atmosphere are extremely difficult, whereas indirect measurements, e.g., in the solar wind, have been performed with considerable precision. The interpretation of solar wind observations with respect to solar surface abundances has been greatly improved in recent years. Abundance measurements have been performed under a large variety of solar wind conditions and refined models have been developed for the transport processes in the chromosphere and the transition region and for the processes occurring in the solar corona. From these measurements we estimate the present isotopic number ratio ³He/⁴He to be (4.1 ± 1.0) × 10^–4 at the solar surface, corresponding to the weight abundance X ₃ = (9.0 ± 2.4) × 10^–5. The zero-age Main-Sequence abundance of ³He (after burning of D) might have been slightly lower (by about 10 to 20%) than the present-day value.Non-standard solar models involving mild turbulent diffusion (Lebreton and Maeder, 1987) could account for a slow secular increase of the ³He/⁴He ratio in the solar atmosphere. On the other hand it is difficult to reconcile models with severe mass loss as proposed by Guzik, Willson, and Brunish (1987) with this constraint. The slowing down of the solar rotation during the early Main-Sequence evolution was accompanied by stronger differential rotation probably implying a more effective mixing of the inner parts. Again, the surface abundance of ³He imposes severe limits on the evolution of the distribution of momentum within the early Sun.     

Electron Temperatures and Flow Speeds of the Low Solar Corona: MACS Results from the Total Solar Eclipse of 29 March 2006 in Libya

An experiment was conducted in conjunction with the total solar eclipse on 29 March 2006 in Libya to measure both the electron temperature and its flow speed simultaneously at multiple locations in the low solar corona by measuring the visible K-coronal spectrum. Coronal model spectra incorporating the effects of electron temperature and its flow speed were matched with the measured K-coronal spectra to interpret the observations. Results show electron temperatures of (1.10±0.05) MK, (0.70±0.08) MK, and (0.98±0.12) MK, at 1.1 R _⊙ from Sun center in the solar north, east and west, respectively, and (0.93±0.12) MK, at 1.2 R _⊙ from Sun center in the solar west. The corresponding outflow speeds obtained from the spectral fit are (103±92) km s⁻¹, (0+10) km s⁻¹, (0+10) km s⁻¹, and (0+10) km s⁻¹. Since the observations were taken only at 1.1 R _⊙ and 1.2 R _⊙ from Sun center, these speeds, consistent with zero outflow, are in agreement with expectations and provide additional confirmation that the spectral fitting method is working. The electron temperature at 1.1 R _⊙ from Sun center is larger at the north (polar region) than the east and west (equatorial region).     

The use of solar eclipse timings to compare the reference systems of Newcomb's Tables of the Sun and of the Improved Lunar Ephemeris

Photoelectric, photographic, and visual observations of second and third contacts of the total solar eclipses of 12 November, 1966 and 7 March, 1970 are analyzed. The observations indicate a difference in the reference system of Newcomb's Tables of the Sun as compared to the reference system of the Improved Lunar Ephemeris (J = 2) of 0^s.911 ± 0^s.034 m.e. The correction is to be applied in the sense of reducing the value of T as tabulated in the American Ephemeris when interpolating the position of the Moon.     

A new outlook on the `Differential Theory' of the solar Quadrupole Moment and Oblateness

The modeling of the quadrupole moment J ₂ and of the oblateness , two key solar parameters, derives from the development in successive spherical harmonics of the gravitational potential. These harmonics are representative of the shape of the Sun, by studying the local distortion of the internal layers, under their distribution of mass and velocity. The first aim of this paper is to study, over the radius r and the colatitude , the structure of the internal layers of the Sun through a geometrical approach, considering J ₂ and under a differential form. The second aim is to determine their theoretical values, after integration over r and , taking the best available models of density and rotation into account constrained by helioseismic data. The novelty of our approach lies in the analysis of the profiles of the two above-mentioned solar parameters, under differential form, from the core to the surface. This analysis allows us to comply with the physical processes located in the transition regions, namely the tachocline and maybe a new subsurface layer which could be called the leptocline. The profiles of tildeJ ₂ show two sharp decreases, which are directly connected to the shear layers located at 0.7 R and beneath the surface. The profiles of tilde show five changes of curvature, which seem to be connected to solar processes, such as the matter circulation flows, seismic events or the storage of the magnetic field, phenomena taking place in the transition regions. These sets of profiles allow us to propose a configuration scenario composed of a double layer. Moreover, as a result of the theoretical determination of tildeJ ₂ and tilde, the values at the surface of the quadrupole moment and of the oblateness can be deduced, which are 1.60×10^–7 and 8.77×10^–6, respectively. As a result of an analysis of available data, we may admit J ₂=(2.0±0.4)×10^–7. The theoretical computations of J ₂ and at the surface will be confronted in the near future with the values measured in space by means of the Picard microsatellite. This mission should permit one to measure at the same time both the total solar irradiance and the latitudinal diameters in any position angle (after removing the passing spots or faculae at the limb).     

New values of gravitational moments J 2 and J 4 deduced from helioseismology

By applying the theory of slowly rotating stars to the Sun, the solar quadrupole and octopole moments J ₂ and J ₄ were computed using a solar model obtained from CESAM stellar evolution code (Morel, 1997) combined with a recent model of solar differential rotation deduced from helioseismology (Corbard et al., 2002). This model takes into account a near-surface radial gradient of rotation which was inferred and quantified from MDI f-mode observations by Corbard and Thompson (2002). The effect of this observational near-surface gradient on the theoretical values of the surface parameters J ₂, J ₄ is investigated. The results show that the octopole moment J ₄ is much more sensitive than the quadrupole moment J ₂ to the subsurface radial gradient of rotation.     

A New Pumping Mechanism for A Series of Class II Methanol Masers in the J 0 - J -1 E Transitions

A new pumping mechanism – methanol masers without population inversion is presented in this paper. It can be used to explain the formation of a series of J ₀ → J _-1 E methanol masers, while the 2₁ → 3₀ A ⁺ methanol masers are regarded as a driving coherent micrwave field. In the new mechanism, the intensities of J ₀ - J _-1 E methanol masers are increased with the decreasing transition frequencies (or with rotational number J, approximately). These results agree with Slysh et al. (1995) and Slysh et al. (1999) J ≤ 5 observations for G3345.01+1.79 and W48, in which both J ₀ → J _-1 E and 2₁ → 3₀ A ⁺ methanol masers are detected coincidentally. Other astronomical conditions, such as magnetic field, 2₁ → 3₀ A ⁺ coherent radition, incoherent pumping rate by thermal radition and so on are also discussed. The new mechanism can operate as a complement to other ordinary maser pumping mechanisms for some class II methonal maser sources. This revised version was published online in July 2006 with corrections to the Cover Date.     

160-min pulsation of the Sun: New observational results

The 1974–1988 Crimean measurements of the solar line-of-sight velocity continue to show the presence of a statistically significant periodicity P ₁ = 160.009 (±) min with an average harmonic amplitude of about 21 cm s^–1. The period is supposed to be that of the global pulsation of the Sun but with a little-known physical mechanism of excitation.The new observations give some evidence for the existence of a second periodicity, P ₁ = 160.014 (±) min. It is hypothesized that the appearance of P ₁ might be a sidelobe mode (of the P ₀-oscillation) induced by rapid rotation of the central solar core.It is also noted that the spacing, in frequency, between P ₀ and P ₁, corresponds to a beat period of 10 ± 3 yr, which happens to be in good agreement with the average duration of the 11 yr cycle of the magnetic activity of the Sun. Accordingly, we suppose that the phase shift of the P ₀-mode between the 1974–1982 and 1986–1988 time intervals reflects a remarkable change of the general magnetic field of the Sun in the course of the 22 yr solar cycle.     

Distribution of Matter in the Solar System

A new formula for the distribution of matter in the solar system is derived by assuming that the planets were formed from trapped particles of a cosmic dust disk attached to the Sun. Contrary to Boltzmann's distribution which predicts thermal collapse of this cloud on the Sun, it is found that if the primeval particles move on circular orbits according to Kepler's law, then their velocities obey a 2-D global Maxwellian and their distribution in space is given by p ₀ (r)=(α r ²)\exp(-α r) (Km^-1); α = 888.73 × 10⁶ Km. The form ofp ₀ (r) agrees with the observed mass distribution of the planets and explains their present large angular momentum. PACS numbers: 96.35.Cp, 96.35.Fs This revised version was published online in July 2006 with corrections to the Cover Date.     

The OSO-1 solar neutron experiment

BF₃ counters on OSO-1 were used to look for solar neutrons by trying to observe a diurnal variation in count rate. No effect was observed and an upper limit was placed on the solar neutron flux at the earth of J _n < 2 × 10^–3 neut/cm²/sec of 10 keV < E _n < 10 meV for the period March to May 1962. No proton-producing flares occurred during this time, so the most obvious source of solar neutrons could not be studied. The emulsion experiment of Apparao et al., flown during this period, which seemed to indicate solar neutrons has probably been misinterpreted.     

Magnetic reconnection as the cause of a photospheric canceling feature and mass flows in a filament

Magnetic reconnection in the temperature minimum region of the solar photosphere can account for the canceling magnetic features on the Sun. Litvinenko (1999a) showed that a reconnection model explains the quiet-Sun features with the magnetic flux cancelation rate of order 10¹⁷ Mx hr⁻¹. In this paper the model is applied to cancelation in solar active regions, which is characterized by a much larger rate of cancelation ∖ ge10¹⁹ Mx hr⁻¹. In particular, the evolution of a photospheric canceling feature observed in an active region on July 2, 1994 is studied. The theoretical predictions are demonstrated to be in reasonable agreement with the measured speed of approaching magnetic fragments, the magnetic field in the fragments, and the flux cancelation rate, deduced from the combined Big Bear Hα time-lapse images and videomagnetograms calibrated against the daily NSO/Kitt Peak magnetogram. Of particular interest is the prediction that photospheric reconnection should lead to a significant upward mass flux and the formation of a solar filament. Hα observations indeed showed a filament that had one of its ends spatially superposed with the canceling feature. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005284116353     

New concept for testing General Relativity: the Laser Astrometric Test Of Relativity (LATOR) mission

This paper discusses a Fundamental physics experiment that will test relativistic gravity at the accuracy better than the effects of the second order in the gravitational field strength, ∝ G². The Laser Astrometric Test Of Relativity (LATOR) mission uses laser interferometry between two micro‐spacecraft whose lines of sight pass close by the Sun to accurately measure deflection of light in the solar gravity. The key element of the experimental design is a redundant geometry optical truss provided by a long‐baseline (100 m) multi‐channel stellar optical interferometer placed on the International Space Station (ISS). The spatial interferometer is used for measuring the angles between the two spacecraft and for orbit determination purposes. In Euclidean geometry, determination of a triangle's three sides determines any angle therein; with gravity changing the optical lengths of sides passing close by the Sun and deflecting the light, the Euclidean relationships are overthrown. The geometric redundancy enables LATOR to measure the departure from Euclidean geometry caused by the solar gravity field to a very high accuracy. LATOR will not only improve the value of the parameterized post‐Newtonian (PPN) γ to unprecedented levels of accuracy of 1 part in 10⁸, it will also reach ability to measure effects of the next post‐Newtonian order (c⁻⁴) of light deflection resulting from gravity's intrinsic non‐linearity. The solar quadrupole moment parameter, J₂, will be measured with high precision, as well as a variety of other relativistic effects including Lense‐Thirring precession. LATOR will lead to very robust advances in the tests of Fundamental physics: this mission could discover a violation or extension of general relativity, or reveal the presence of an additional long range interaction in the physical law. There are no analogs to the LATOR experiment; it is unique and is a natural culmination of solar system gravity experiments. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Electron densities derived from line intensity ratios: Beryllium isoelectronic sequence

A direct method for determining electron densities from emission line intensities of ions in the beryllium isoelectronic sequence is described and then applied to the analysis of extreme ultraviolet Ciii and Ov spectra from both quiet and active areas in the solar transition region. The results are consistent with a value of N _e T _e = 6 × 10¹⁴ cm^-3K for the quiet Sun at temperatures of 5 × 10⁴ to 3 × 10⁵K. Electron densities are approximately five times greater in active regions than in the quiet Sun.     

SUMER - Solar Ultraviolet Measurements of Emitted Radiation

The instrument SUMER - Solar Ultraviolet Measurements of Emitted Radiation is designed to investigate structures and associated dynamical processes occurring in the solar atmosphere, from the chromosphere through the transition region to the inner corona, over a temperature range from 10⁴ to 2 × 10⁶ K and above. These observations will permit detailed spectroscopic diagnostics of plasma densities and temperatures in many solar features, and will support penetrating studies of underlying physical processes, including plasma flows, turbulence and wave motions, diffusion transport processes, events associated with solar magnetic activity, atmospheric heating, and solar wind acceleration in the inner corona. Specifically, SUMER will measure profiles and intensities of EUV lines; determine Doppler shifts and line broadenings with high accuracy; provide stigmatic images of the Sun in the EUV with high spatial, spectral, and temporal resolution; and obtain monochromatic maps of the full Sun and the inner corona or selected areas thereof. SUMER will be flown on the Solar and Heliospheric Observatory (SOHO), scheduled for launch in November, 1995. This paper has been written to familiarize solar physicists with SUMER and to demonstrate some command procedures for achieving certain scientific observations.     

The Influence of the Electron κ-Distribution in the Solar Corona on the Fe VIII – Fe XV Line Intensities

Changes in the excitation equilibrium of Feviii – Fe xv in the solar corona due to the electron non-thermal κ-distribution are studied. The shape of the distribution affects the electron excitation rate and thus the relative intensities of the spectral lines. Since the shape of the electron distribution function influences also the ionization equilibrium of Fe, both effects change the final intensities of the lines. Possibilities for diagnostics of the shape of the electron distribution in the solar corona are discussed. Synthetic spectrum of Fe for T = 1.58× 10⁶ K and n_e = 10⁸ cm⁻³ is shown together with the synthetic spectra computed with DEM for the quiet Sun.     

The effect of dense cores on the structure and evolution of Jupiter and Saturn

We have calculated evolutionary and static models of Jupiter and Saturn with homogeneous solar composition mantles and dense cores of material consisting of solar abundances of SiO₂, MgO, Fe, and Ni. Evolutionary sequences for Jupiter were calculated with cores of mass 2, 4, 6, and 8% of the Jovian mass. Evolutionary sequences for Saturn were calculated with cores of mass 16, 18, 20, and 22% of total mass. Two envelope mixtures, representative of the solar abundances were used: X (mass fraction of hydrogen) = 0.74, Y (mass fraction of helium) = 0.24 and X = 0.77 and Y = 0.21. For Jupiter, the observations of the temperature at 1 bar pressure (T_1bar), radius and internal luminosity were best fit by evolutionary models with a core mass of ～6.5% and chemical composition of X = 0.77, Y = 0.21. The calculated cooling time for Jupiter is approximately 4.9 × 10⁹ years, which is consistent, within our error bars, with the known age of the solar system. For Saturn, the observations of the radius, internal luminosity and T_1BAR can be best fit by evolutionary models with a core mass of ～21% and chemical composition of X = 0.77, Y = 0.21. The cooling time calculated for Saturn is approximately 2.6 × 10⁹ years, almost a factor 2 less than the present age of the solar system. Static models of Jupiter and Saturn were calculated for the above chemical compositions in order to investigate the sensitivity of the calculated gravitational moments, J₂ and J₄, to the mass of the dense core, T_1BAR and hydrogen/helium ratio. We find for Jupiter that a model having a core mass of approximately 7% gives values of J₂, J₄, and T_1BAR that are within observational limits, for the mixture X = 0.77, Y = 0.21. The static Jupiter models are completely consistent with the evolutionary results. For Saturn, the quantities J₂, J₄, and J₆ determined from the static models with the most probable T_1BAR of 140°K, using modeling procedures which result in consistent models for Jupiter, are considerably below the observed values.