  We examine contributions of second order physical processes to results of
stellar evolution calculations amenable to direct observational testing. In the
first paper in the series (Young et al. 2001) we established baseline results
using only physics which are common to modern stellar evolution codes. In the
current paper we establish how much of the discrepancy between observations and
baseline models is due to particular elements of new physics. We then consider
the impact of the observational uncertainties on the maximum predictive
accuracy achievable by a stellar evolution code. The sun is an optimal case
because of the precise and abundant observations and the relative simplicity of
the underlying stellar physics. The Standard Model is capable of matching the
structure of the sun as determined by helioseismology and gross surface
observables to better than a percent. Given an initial mass and surface
composition within the observational errors, and no additional constraints for
which the models can be optimized, it is not possible to predict the sun's
current state to better than ~7%. Convectively induced mixing in radiative
regions, seen in multidimensional hydrodynamic simulations, dramatically
improves the predictions for radii, luminosity, and apsidal motions of
eclipsing binaries while simultaneously maintaining consistency with observed
light element depletion and turnoff ages in young clusters (Young et al. 2003).
Systematic errors in core size for models of massive binaries disappear with
more complete mixing physics, and acceptable fits are achieved for all of the
binaries without calibration of free parameters. The lack of accurate abundance
determinations for binaries is now the main obstacle to improving stellar
models using this type of test.
