While considerable advance has been made to account for statistical uncertainties in astronomical analyses, systematic instrumental uncertainties have been generally ignored. This can be crucial to a proper interpretation of analysis results because instrumental calibration uncertainty is a form of systematic uncertainty. Ignoring it can underestimate error bars and introduce bias into the fitted values of model parameters. Accounting for such uncertainties currently requires extensive case-specific simulations if using existing analysis packages. Here we present general statistical methods that incorporate calibration uncertainties into spectral analysis of high-energy data. We first present a method based on multiple imputation that can be applied with any fitting method, but is necessarily approximate. We then describe a more exact Bayesian approach that works in conjunction with a Markov chain Monte Carlo based fitting. We explore methods for improving computational efficiency, and in particular detail a method of summarizing calibration uncertainties with a principal component analysis of samples of plausible calibration files. This method is implemented using recently codified Chandra effective area uncertainties for low-resolution spectral analysis and is verified using both simulated and actual Chandra data. Our procedure for incorporating effective area uncertainty is easily generalized to other types of calibration uncertainties.

We summarize and critically evaluate the available data on nuclear fusion cross sections important to energy generation in the Sun and other hydrogen-burning stars and to solar neutrino production. Recommended values and uncertainties are provided for key cross sections, and a recommended spectrum is given for 8B solar neutrinos. We also discuss opportunities for further increasing the precision of key rates, including new facilities, new experimental techniques, and improvements in theory. This review, which summarizes the conclusions of a workshop held at the Institute for Nuclear Theory, Seattle, in January 2009, is intended as a 10-year update and supplement to Reviews of Modern Physics 70 (1998) 1265.

The detection and quantification of narrow emission lines in X-ray spectra is a challenging statistical task. The Poisson nature of the photon counts leads to local random fluctuations in the observed spectrum that often results in excess emission in a narrow band of energy resembling a weak narrow line. From a formal statistical perspective, this leads to a (sometimes highly) multimodal likelihood. Many standard statistical procedures are based on (asymptotic) Gaussian approximations to the likelihood and simply cannot be used in such settings. Bayesian methods offer a more direct paradigm for accounting for such complicated likelihood functions but even here multimodal likelihoods pose significant computational challenges. The new Markov chain Monte Carlo (MCMC) methods developed in 2008 by van Dyk and Park, however, are able to fully explore the complex posterior distribution of the location of a narrow line, and thus provide valid statistical inference. Even with these computational tools, standard statistical quantities such as means and standard deviations cannot adequately summarize inference and standard testing procedures cannot be used to test for emission lines. In this paper, we use new efficient MCMC algorithms to fit the location of narrow emission lines, we develop new statistical strategies for summarizing highly multimodal distributions and quantifying valid statistical inference, and we extend the method of posterior predictive p-values proposed by Protassov et al. (2002) to test for the presence of narrow emission lines in X-ray spectra. We illustrate and validate our methods using simulation studies and apply them to the Chandra observations of the high redshift quasar PG1634+706.

A commonly used measure to summarize the nature of a photon spectrum is the so-called Hardness Ratio, which compares the number of counts observed in different passbands. The hardness ratio is especially useful to distinguish between and categorize weak sources as a proxy for detailed spectral fitting. However, in this regime classical methods of error propagation fail, and the estimates of spectral hardness become unreliable. Here we develop a rigorous statistical treatment of hardness ratios that properly deals with detected photons as independent Poisson random variables and correctly deals with the non-Gaussian nature of the error propagation. The method is Bayesian in nature, and thus can be generalized to carry out a multitude of source-population--based analyses. We verify our method with simulation studies, and compare it with the classical method. We apply this method to real world examples, such as the identification of candidate quiescent Low-mass X-ray binaries in globular clusters, and tracking the time evolution of a flare on a low-mass star.

An effective field theory technique that combines the standard nuclear physics approach and chiral perturbation theory is applied to the $hen$ process, ${}^{3}{He}+n\to {}^{4}{He}+ \gamma$. For the initial and final nuclear states, high-precision wave functions are generated via the variational Monte Carlo method using the Argonne $v_{14}$ potential and Urbana VIII trinucleon interactions, while the relevant transition operators are calculated up to ${\cal O}(Q^4)$ in HB$\chi$PT. The imposition of the renormalization condition that the magnetic moments of ${}^{3}{He}$ and ${}^{3}{H}$ be reproduced allows us to carry out a parameter-free calculation of the $hen$ cross section. The result, $\sigma=(60\pm 3\pm 1) \mu b$, is in reasonable agreement with the experimental values, $(54\pm 6) \mu b$ and $(55\pm 3) \mu b$. This agreement demonstrates the validity of the calculational method previously used for estimating the reaction rate of the solar $hep$ process.

In interpreting the SNO experiments, accurate estimates of the \nu d reaction cross sections are of great importance. In our recent work, we have improved our previous calculation by updating some of its inputs and by incorporating the results of a recent effective-field-theoretical calculation. The new cross sections are slightly (\sim 1%) larger than the previously reported values. It is reasonable to assign 1% uncertainty to the \nu d cross sections reported here; this error estimate does not include radiative corrections.

Spurred by the recent complete determination of the weak currents in two-nucleon systems up to ${\cal O}(Q^3)$ in heavy-baryon chiral perturbation theory, we carry out a parameter-free calculation of the threshold $S$-factors for the solar $pp$ (proton-fusion) and $hep$ processes in an effective field theory that {\it combines} the merits of the standard nuclear physics method and systematic chiral expansion. The power of the EFT adopted here is that one can correlate in a unified formalism the weak-current matrix elements of two-, three- and four-nucleon systems. Using the tritium $\beta$-decay rate as an input to fix the only unknown parameter in the theory, we can evaluate the threshold $S$ factors with drastically improved precision; the results are $S_{pp}(0) = 3.94\times(1 \pm 0.004) \times 10^{-25} {MeV-b}$ and $S_{hep}(0) = (8.6\pm 1.3)\times 10^{-20} {keV-b}$. The dependence of the calculated $S$-factors on the momentum cutoff parameter $\Lambda$ has been examined for a physically reasonable range of $\Lambda$. This dependence is found to be extremely small for the $pp$ process, and to be within acceptable levels for the $hep$ process, substantiating the consistency of our calculational scheme.

In interpreting the SNO experiments, very accurate estimates of the \nu d reaction cross sections are of great importance. We improve the previous estimates of our group by updating some of its inputs and by taking into account the results of a recent effective-field-theoretical calculation. The new cross sections are slightly (\sim 1%) larger than the previously reported values. We present arguments that lead to the conclusion that it is reasonable to assign 1% uncertainty to the \nu d cross sections reported here.

The strategy of modern effective field theory is exploited to pin down accurately the flux $S$ factors for the $pp$ and $hep$ processes in the Sun. The technique used is to combine the high accuracy established in few-nucleon systems of the "standard nuclear physics approach" (SNPA) and the systematic power counting of chiral perturbation theory (ChPT) into a consistent effective field theory framework. Using highly accurate wave functions obtained in the SNPA and working to \nlo3 in the chiral counting for the current, we make totally parameter-free and error-controlled predictions for the $pp$ and $hep$ processes in the Sun.

Using effective field theory, we calculate the S-factor for the hep process in a totally parameter-free formulation. The transition operators are organized according to chiral counting, and their matrix elements are evaluated using the realistic nuclear wave functions obtained in the correlated-hyperspherical-harmonics method. Terms of up to next-to-next-to-next-to-leading order in heavy-baryon chiral perturbation theory are considered. Fixing the only parameter in the theory by fitting the tritium \beta-decay rate, we predict the hep S-factor with accuracy better than \sim 20 %.

Spurred by the recent complete determination of the weak currents in two-nucleon systems up to ${\cal O}(Q^3)$ in heavy-baryon chiral perturbation theory, we carry out a parameter-free calculation of the solar proton fusion rate in an effective field theory that combines the merits of the standard nuclear physics method and systematic chiral expansion. Using the tritium beta-decay rate as an input to fix the only unknown parameter in the effective Lagrangian, we can evaluate with drastically improved precision the ratio of the two-body contribution to the well established one-body contribution; the ratio is determined to be (0.86\pm 0.05) %. This result is essentially independent of the cutoff parameter for a wide range of its variation (500 MeV \le \Lambda \le 800 MeV), a feature that substantiates the consistency of the calculation.

The proton burning process p + p -> d + e(+) + \nu(e), important for the stellar evolution of main-sequence stars of mass equal to or less than that of the Sun, is computed in effective field theory using chiral perturbation expansion to the next-to-next-to leading chiral order. This represents a model-independent calculation consistent with low-energy effective theory of QCD comparable in accuracy to the radiative np capture at thermal energy previously calculated by first using very accurate two-nucleon wavefunctions backed up by an effective field theory technique with a finite cut-off. The result obtained thereby is found to support within theoretical uncertainties the previous calculation of the same process by Bahcall and his co-workers.