Our Sun. IV. The Standard Model and Helioseismology: Consequences of Uncertainties in Input Physics and in Observed Solar Parameters

Helioseismology provides a powerful tool to explore the deep interior of the Sun: for example, the adiabatic sound speed can be inferred with an accuracy of a few parts in 10,000. This has become a serious challenge to theoretical models of the Sun. Therefore, we have undertaken a self-consistent, systematic study of sources of uncertainties in the standard solar model, which must be understood before the helioseismic observations can be used as constraints on theory. We find that the largest uncertainty in the sound speed in the solar interior, namely, 3 parts in 1000, arises from uncertainties in the observed photospheric abundances of the elements; uncertainties of 1 part in 1000 arise from (1) the 4% uncertainty in the OPAL opacities, (2) the 5% uncertainty in the basic pp nuclear reaction rate, (3) the 15% uncertainty in the diffusion constants for the gravitational settling of helium, and (4) the 50% uncertainties in diffusion constants for the heavier elements. (Other investigators have shown that similar uncertainties arise from uncertainties in the interior equation of state and in rotation-induced turbulent mixing.) The predicted pre-main-sequence solar lithium depletion is a factor of order 20 (an order of magnitude larger than that predicted by earlier models that neglected gravitational settling and used older opacities), and is uncertain by a factor of 2. The predicted neutrino capture rate is uncertain by 30% for the Cl-37 experiment and by 3% for the Ga-71 experiments (not including uncertainties in the capture cross sections), while the B-8 neutrino flux is uncertain by 30%.Comment: LaTeX, 38 pages (including 8 figures); ApJ, in press. Added figures/color figurea available at http://www.cita.utoronto.ca/~boothroy/sun4.htm

Dissolution of coccolithophorid calcite by microzooplankton and copepod grazing

International audienceIndependent of the ongoing acidification of surface seawater, the majority of the calcium carbonate produced in the pelagial is dissolved by natural processes above the lysocline. We investigate to what extent grazing and passage of coccolithophorids through the guts of copepods and the food vacuoles of microzooplankton contribute to calcite dissolution. In laboratory experiments where the coccolithophorid Emiliania huxleyi was fed to the rotifer Brachionus plicatilis, the heterotrophic flagellate Oxyrrhis marina and the copepod Acartia tonsa, calcite dissolution rates of 45?55%, 37?53% and 5?22% of ingested calcite were found. We ascribe higher loss rates in microzooplankton food vacuoles as compared to copepod guts to the strongly acidic digestion and the individual packaging of algal cells. In further experiments, specific rates of calcification and calcite dissolution were also measured in natural populations during the PeECE III mesocosm study under differing ambient pCO2 concentrations. Microzooplankton grazing accounted for between 27 and 70% of the dynamic calcite stock being lost per day, with no measurable effect of CO2 treatment. These measured calcite dissolution rates indicate that dissolution of calcite in the guts of microzooplankton and copepods can account for the calcite losses calculated for the global ocean using budget and model estimates

Solar Cycle Related Changes at the Base of the Convection Zone

The frequencies of solar oscillations are known to change with solar activity. We use Principal Component Analysis to examine these changes with high precision. In addition to the well-documented changes in solar normal mode oscillations with activity as a function of frequency, which originate in the surface layers of the Sun, we find a small but statistically significant change in frequencies with an origin at and below the base of the convection zone. We find that at r=(0.712^{+0.0097}_{-0.0029})R_sun, the change in sound speed is \delta c^2 / c^2 = (7.23 +/- 2.08) x 10^{-5} between high and low activity. This change is very tightly correlated with solar activity. In addition, we use the splitting coefficients to examine the latitudinal structure of these changes. We find changes in sound speed correlated with surface activity for r >~ 0.9R_sun.Comment: 29 pages, 11 figures, accepted for publication in Ap

Non-LTE Model Atmospheres for Late-Type Stars II. Restricted NLTE Calculations for a Solar-Like Atmosphere

We test our knowledge of the atomic opacity in the solar UV spectrum. Using the atomic data compiled in Paper I from modern, publicly available, databases, we perform calculations that are confronted with space-based observations of the Sun. At wavelengths longer than about 260 nm, LTE modeling can reproduce quite closely the observed fluxes; uncertainties in the atomic line data account fully for the differences between calculated and observed fluxes. At shorter wavelengths, departures from LTE appear to be important, as our LTE and restricted NLTE calculations differ. Analysis of visible-near infrared Na I and O I lines, two species that produce a negligible absorption in the UV, shows that observed departures from LTE for theses species can be reproduced very accurately with restricted (fixed atmospheric structure) NLTE calculations.Comment: 13 pages, 11 figures, to appear in Ap

X-ray and UV spectroscopy of Galactic diffuse hot gas along the LMC X--3 sight line

We present Suzaku spectra of X-ray emission in the fields just off the LMC X-3 sight line. OVII, OVIII, and NeIX emission lines are clearly detected, suggesting the presence of an optically thin thermal plasma with an average temperature of 2.4E6. This temperature is significantly higher than that inferred from existing X-ray absorption line data obtained with Chandra grating observations of LMC X-3, strongly suggesting that the gas is not isothermal. We then jointly analyze these data to characterize the spatial and temperature distributions of the gas. Assuming a vertical exponential Galactic disk model, we estimate the gas temperature and density at the Galactic plane and their scale heights as 3.6(2.9, 4.7)E6 K and 1.4(0.3, 3.4)E-3 cm^{-3} and 1.4(0.2, 5.2) kpc and $2.8(1.0, 6.4)$ kpc, respectively. This characterization can account for all the \ovi line absorption, as observed in a FUSE spectrum of LMC X-3, but only predicts less than one tenth of the OVI line emission intensity typically detected at high Galactic latitudes. The bulk of the OVI emission most likely arises at interfaces between cool and hot gases.Comment: 10 pages, 7 figures, 3 tables, accepted for publication in ApJ, 200

Can three-flavor oscillations solve the solar neutrino problem?

The most promising solution to the solar neutrino problem are neutrino oscillations, which usually are analyzed within the reduced 2-flavor scheme, because the solutions found therein reasonably well reproduce the recent data of Super-Kamiokande about the recoil-electron energy spectrum, zenith-angle and seasonal variations, and the event rate data of all the neutrino detectors. In this work, however, a survey of the complete parameter space of 3-flavor oscillations is performed. Basically eight new additional solutions could be identified, where the best one with \Delta m(12)^2=2.7x10^(-10) eV^2, \Delta m(13)^2=1.0x10^(-5) eV^2, \Theta(12)=23, and \Theta(13)=1.3 (denoted SVO) is slightly more probable than any 2-flavor solution. However, including the results of the atmospheric neutrino problem excludes all 3-flavour solutions apart from the SLMA-solution (\Delta m(12)^2=7.9x10^(-6) eV^2, \Delta m(13)^2=2.5x10^(-4) eV^2, \Theta(12)=1.4, and \Theta(13)=20). Besides, the ability of SNO and Borexino to discriminate the various 2- and 3-flavor solutions is investigated. Only with very good statistics in these experiments the correct solution to the solar neutrino problem can be identified unambiguously.Comment: 22 pages, 19 figures, REVTeX, submitted to Phys.Rev.D, article with better resolved figures available under http://www.mpa-garching.mpg.de/~schlattl/public.htm

Solar Models: current epoch and time dependences, neutrinos, and helioseismological properties

We calculate accurate solar models and report the detailed time dependences of important solar quantities. We use helioseismology to constrain the luminosity evolution of the sun and report the discovery of semi-convection in evolved solar models that include diffusion. In addition, we compare the computed sound speeds with the results of p-mode observations by BiSON, GOLF, GONG, LOWL, and MDI instruments. We contrast the neutrino predictions from a set of eight standard-like solar models and four deviant (or deficient) solar models with the results of solar neutrino experiments. For solar neutrino and for helioseismological applications, we present present-epoch numerical tabulations of characteristics of the standard solar model as a function of solar radius, including the principal physical and composition variables, sound speeds, neutrino fluxes, and functions needed for calculating solar neutrino oscillations.Comment: Accepted ApJ. Have used refined satellite value for solar luminosity. Changes slightly best neutrino fluxes. Include new references, number density of scatterers of sterile neutrinos, some additional helioseismological predictions. 70 pages, 16 figures, additional material at http://www.sns.ias.edu/~jn

Equation of state of fully ionized electron-ion plasmas

Thermodynamic quantities of Coulomb plasmas consisting of point-like ions immersed in a compressible, polarizable electron background are calculated for ion charges Z=1 to 26 and for a wide domain of plasma parameters ranging from the Debye-Hueckel limit to the crystallization point and from the region of nondegenerate to fully degenerate nonrelativistic or relativistic electrons. The calculations are based on the linear-response theory for the electron-ion interaction, including the local-field corrections in the electronic dielectric function. The thermodynamic quantities are calculated in the framework of the N-body hypernetted-chain equations and fitted by analytic expressions. We present also accurate analytic approximations for the free energy of the ideal electron gas at arbitrary degeneracy and relativity and for the excess free energy of the one-component plasma of ions (OCP) derived from Monte Carlo simulations. The extension to multi-ionic mixtures is discussed within the framework of the linear mixing rule. These formulae provide a completely analytic, accurate description of the thermodynamic quantities of fully ionized electron-ion Coulomb plasmas, a useful tool for various applications from liquid state theory to dense stellar matter.Comment: 13 pages, 2 tables, 7 figures, REVTeX using epsf.sty. To be published in Phys. Rev. E, vol. 58 (1998

A Comparison of Solar Cycle Variations in the Equatorial Rotation Rates of the Sun's Subsurface, Surface, Corona, and Sunspot Groups

Using the Solar Optical Observing Network (SOON) sunspot-group data for the period 1985-2010, the variations in the annual mean equatorial-rotation rates of the sunspot groups are determined and compared with the known variations in the solar equatorial-rotation rates determined from the following data: i) the plasma rotation rates at 0.94Rsun, 0.95Rsun,...,1.0Rsun measured by Global Oscillation Network Group (GONG) during the period 1995-2010, ii) the data on the soft X-ray corona determined from Yohkoh/SXT full disk images for the years 1992-2001, iii) the data on small bright coronal structures (SBCS) which were traced in Solar and Heliospheric Observatory (SOHO)/EIT images during the period 1998-2006, and iv) the Mount Wilson Doppler-velocity measurements during the period 1986-2007. A large portion (up to approximate 30 deg latitude) of the mean differential-rotation profile of the sunspot groups lies between those of the internal differential-rotation rates at 0.94Rsun and 0.98Rsun.The variation in the yearly mean equatorial-rotation rate of the sunspot groups seems to be lagging that of the equatorial-rotation rate determined from the GONG measurements by one to two years.The amplitude of the latter is very small.The solar-cycle variation in the equatorial-rotation rate of the solar corona closely matches that determined from the sunspot-group data.The variation in the equatorial-rotation rate determined from the Mount Wilson Doppler-velocity data closely resembles the corresponding variation in the equatorial-rotation rate determined from the sunspot-group data that included the values of the abnormal angular motions (> 3 deg per day) of the sunspot groups. Implications of these results are pointed out.Comment: 22 pages, 10 figures, accepted by Solar Physic

Forbush decreases and turbulence levels at CME fronts

We seek to estimate the average level of MHD turbulence near coronal mass ejection (CME) fronts as they propagate from the Sun to the Earth. We examine the cosmic ray data from the GRAPES-3 tracking muon telescope at Ooty, together with the data from other sources for three well observed Forbush decrease events. Each of these events are associated with frontside halo Coronal Mass Ejections (CMEs) and near-Earth magnetic clouds. In each case, we estimate the magnitude of the Forbush decrease using a simple model for the diffusion of high energy protons through the largely closed field lines enclosing the CME as it expands and propagates from the Sun to the Earth. We use estimates of the cross-field diffusion coefficient $D_{\perp}$ derived from published results of extensive Monte Carlo simulations of cosmic rays propagating through turbulent magnetic fields. Our method helps constrain the ratio of energy density in the turbulent magnetic fields to that in the mean magnetic fields near the CME fronts. This ratio is found to be $\sim$ 2% for the 11 April 2001 Forbush decrease event, $\sim$ 6% for the 20 November 2003 Forbush decrease event and $\sim$ 249% for the much more energetic event of 29 October 2003.Comment: Accepted for publication in Astronomy and Astrophysics. (Abstract abridged) Typos correcte