Recently developed analytic approximation for the equation of state of fully ionized nonideal electron-ion plasma mixtures [Potekhin et al., Phys. Rev. E, 79, 016411 (2009); arXiv:0812.4344], which covers the transition between the weak and strong Coulomb coupling regimes and reproduces numerical results obtained in the hypernetted chain (HNC) approximation, is modified in order to fit the small deviations from the linear mixing in the strong coupling regime, revealed by recent Monte Carlo simulations. In addition, a mixing rule is proposed for the regime of weak coupling, which generalizes post-Debye density corrections to the case of mixtures and numerically agrees with the HNC approximation in that regime.

We compare Path Integral Monte Carlo calculations by Militzer and Pollock (Phys. Rev. B 71, 134303, 2005) of Coulomb tunneling in nuclear reactions in dense matter to semiclassical calculations assuming WKB Coulomb barrier penetration through the radial mean-field potential. We find a very good agreement of two approaches at temperatures higher than ~1/5 of the ion plasma temperature. We obtain a simple parameterization of the mean field potential and of the respective reaction rates. We analyze Gamow-peak energies of reacting ions in various reaction regimes and discuss theoretical uncertainties of nuclear reaction rates taking carbon burning in dense stellar matter as an example.