State-of-the-art stellar structure and evolution codes fail to adequately describe turbulent convection. For stars with convective envelopes, such as red giants, this leads to an incomplete depiction of the surface layers. As a result, the predicted stellar oscillation frequencies are haunted by systematic errors, the so-called surface effect. Different empirically and theoretically motivated correction relations have been proposed to deal with this issue. In this paper, we compare the performance of these surface correction relations for red giant branch stars. For this purpose, we apply the different surface correction relations in asteroseismic analyses of eclipsing binaries and open clusters. In accordance with previous studies of main-sequence stars, we find that the use of different surface correction relations biases the derived global stellar properties, including stellar age, mass, and distance estimates. We furthermore demonstrate that the different relations lead to the same systematic errors for two different open clusters. Our results overall discourage from the use of surface correction relations that rely on reference stars to calibrate free parameters. Due to the demonstrated systematic biasing of the results, the use of appropriate surface correction relations is imperative to any asteroseismic analysis of red giants. Accurate mass, age, and distance estimates for red giants are fundamental when addressing questions that deal with the chemo-dynamical evolution of the Milky Way galaxy. In this way, our results also have implications for fields, such as galactic archaeology, that draw on findings from stellar physics.

Models of stellar structure and evolution are an indispensable tool in astrophysics, yet they are known to incorrectly reproduce the outer convective layers of stars. In the first paper of this series, we presented a novel procedure to include the mean structure of 3D hydrodynamical simulations on-the-fly in stellar models, and found it to significantly improve the outer stratification and oscillation frequencies of a standard solar model. In the present work, we extend the analysis of the method; specifically how the transition point between envelope and interior affects the models. We confirm the versatility of our method by successfully repeating the entire procedure for a different grid of 3D hydro-simulations. Furthermore, the applicability of the procedure was investigated across the HR diagram and an accuracy comparable to the solar case was found. Moreover, we explored the implications on stellar evolution and find that the red-giant branch is shifted about 40 K to higher effective temperatures. Finally, we present for the first time an asteroseismic analysis based on stellar models fully utilising the stratification of 3D simulations on-the-fly. These new models significantly reduce the asteroseismic surface term for the two selected stars in the Kepler field. We extend the analysis to red giants and characterise the shape of the surface effect in this regime. Lastly, we stress that the interpolation required by our method would benefit from new 3D simulations, resulting in a finer sampling of the grid.

Stellar models utilising one-dimensional (1D), heuristic theories of convection fail to adequately describe the energy transport in superadiabatic layers. The improper modelling leads to well-known discrepancies between observed and predicted oscillation frequencies for stars with convective envelopes. Recently, three-dimensional (3D) hydrodynamic simulations of stellar envelopes have been shown to facilitate a realistic depiction of superadiabatic convection in 1D stellar models. The resulting structural changes of the boundary layers have been demonstrated to impact not only the predicted oscillation spectra but evolution tracks as well. In this paper, we quantify the consequences that the change in boundary conditions has for stellar parameter estimates of main-sequence stars. For this purpose, we investigate two benchmark stars, Alpha Centauri A and B, using Bayesian inference. We show that the improved treatment of turbulent convection makes the obtained 1D stellar structures nearly insensitive to the mixing length parameter. By using 3D simulations in 1D stellar models, we hence overcome the degeneracy between the mixing length parameter and other stellar parameters. By lifting this degeneracy, the inclusion of 3D simulations has the potential to yield more robust parameter estimates. In this way, a more realistic depiction of superadiabatic convection has important implications for any field that relies on stellar models, including the study of the chemical evolution of the Milky Way Galaxy and exoplanet research.

State-of-the-art one-dimensional (1D) stellar evolution codes rely on simplifying assumptions, such as mixing length theory, in order to describe superadiabatic convection. As a result, 1D stellar structure models do not correctly recover the surface layers of the Sun and other stars with convective envelopes. We present a method that overcomes this structural drawback by employing three-dimensional (3D) hydrodynamic simulations of stellar envelopes: at every time-step of the evolution interpolated 3D envelopes are appended to the 1D structure and are used to supply realistic boundary conditions for the stellar interior. In contrast to previous attempts, our method includes mean 3D turbulent pressure. We apply our method to model the present Sun. The structural shortcomings of standard stellar models lead to systematic errors in the stellar oscillation frequencies inferred from the model. We show that our method fully corrects for this error. Furthermore, we show that our realistic treatment of superadiabatic convection alters the predicted evolution of the Sun. Our results hence have important implications for the characterization of stars. This has ramifications for neighbouring fields, such as exoplanet research and galactic archaeology, for which accurate stellar models play a key role.

1D stellar evolution codes employ rudimentary treatments of turbulent convection. For stars with convective envelopes, this leads to systematic errors in the predicted oscillation frequencies needed for asteroseismology. One way of mending these structural inadequacies is through patching, whereby the outermost layers of 1D models are replaced by the mean stratifications from 3D simulations. In order to viably implement this approach in asteroseismic analysis, interpolation throughout precomputed 3D envelopes is required. We present a method that interpolates throughout precomputed 3D envelopes as a function of effective temperature, surface gravity, and metallicity. We conduct a series of validation tests that demonstrate that the scheme reliably and accurately reproduces the structures of stellar envelopes and apply our method to the Sun as well as two stars observed by Kepler. We parameterize the frequency shift that results from patching and show that the functional forms are evolutionary dependent. In addition we find that neglecting modal effects, such as non-adiabatic energetics, introduces systematic errors in asteroseimically obtained stellar parameters. Both these results suggest that a cautious approach is necessary when utilizing empirical surface corrections in lieu of patching models. Our results have important implications, particularly for characterizing exoplanet systems, where accuracy is of utmost concern.

Current standard solar models neither account properly for the photospheric lithium abundance nor reproduce the inferred solar sound speed profile. Diffusive overshooting at the base of the solar convective envelope has previously been shown to solve either of these model inadequacies. In this paper, we present an analysis of solar models with four different parametrizations of diffusive overshooting. We find that these models are able to recover the correct lithium depletion, regardless of the parametrization, if overshooting is suppressed, during the early evolutionary stages. Further, parametrizations of diffusive overshooting have been shown to improve the inferred sound speed profile. However, none of the presented models are able to simultaneously solve both model inadequacies, showing that diffusive overshooting on its own is deficient to account for observations.

We present a comparison between Monte Carlo simulations and a semi-analytical approach that reproduces the theoretical probability distribution functions of the solar neutrino fluxes, stemming from the $pp$, $pep$, $hep$, $^7\mathrm{Be}$, $^8\mathrm{B}$, $^{13}\mathrm{N}$, $^{15}\mathrm{O}$, and $^{17}\mathrm{F}$ source reactions. We obtain good agreement between the two approaches. Thus, the semi-analytical method yields confidence intervals that closely match those found, based on Monte Carlo simulations, and points towards the same general symmetries of the investigated probability distribution functions. Furthermore, the negligible computational cost of this method is a clear advantage over Monte Carlo simulations, making it trivial to take new observational constraints on the input parameters into account.

Standard 1D stellar evolution models do not correctly reproduce the structure of the outermost layers of stars with convective envelopes. This has been a long-standing problem in stellar modelling affecting both the predicted evolutionary paths and the attributed oscillation frequencies, and indirectly biasing numerous quantities derived from stellar evolution calculations. We present a novel method that mostly eliminates these structural defects by appending mean 3D simulations of stellar envelopes. In contrast to previous attempts we impose the complete structure derived from 3D simulations at each time step during the entire evolution. For this purpose, we interpolate in grids of pre-computed 3D simulations and use the resulting structure as boundary conditions, in order to solve the stellar structure equations for the 1D interior at each time step. Our method provides a continuous transition in many quantities from the interior to the imposed interpolated 3D surface layers. We present a solar calibration model and show that the obtained structure of the surface layers reliably mimics that of the underlying 3D simulations for the present Sun. Moreover, we perform a helioseismic analysis, showing that our method mostly eliminates the structural contribution to the discrepancy between model frequencies and observed p-mode frequencies.

We present a new method for replacing the outermost layers of stellar models with interpolated atmospheres based on results from 3D simulations, in order to correct for structural inadequacies of these layers. This replacement is known as patching. Tests, based on 3D atmospheres from three different codes and interior models with different input physics, are performed. Using solar models, we investigate how different patching criteria affect the eigenfrequencies. These criteria include the depth, at which the replacement is performed, the quantity, on which the replacement is based, and the mismatch in $T_\mathrm{eff}$ and $\log g$ between the un-patched model and patched 3D atmosphere. We find the eigenfrequencies to be unaltered by the patching depth deep within the adiabatic region, while changing the patching quantity or the employed atmosphere grid leads to frequency shifts that may exceed $1\,\mu\mathrm{Hz}$. Likewise, the eigenfrequencies are sensitive to mismatches in $T_\mathrm{eff}$ or $\log g$. A thorough investigation of the accuracy of a new scheme, for interpolating mean 3D stratifications within the atmosphere grids, is furthermore performed. Throughout large parts of the atmosphere grids, our interpolation scheme yields sufficiently accurate results for the purpose of asteroseismology. We apply our procedure in asteroseismic analyses of four \textit{Kepler} stars and draw the same conclusions as in the solar case: Correcting for structural deficiencies lowers the eigenfrequencies, this correction is slightly sensitive to the patching criteria, and the remaining frequency discrepancy between models and observations is less frequency-dependent. Our work shows the applicability and relevance of patching in asteroseismology.