In this work we use Bayesian inference to quantitatively reconstruct the solar properties most relevant to the solar composition problem using as inputs the information provided by helioseismic and solar neutrino data. In particular, we use a Gaussian process to model the functional shape of the opacity uncertainty to gain flexibility and become as free as possible from prejudice in this regard. With these tools we first readdress the statistical significance of the solar composition problem. Furthermore, starting from a composition unbiased set of standard solar models we are able to statistically select those with solar chemical composition and other solar inputs which better describe the helioseismic and neutrino observations. In particular, we are able to reconstruct the solar opacity profile in a data driven fashion, independently of any reference opacity tables, obtaining a 4% uncertainty at the base of the convective envelope and 0.8% at the solar core. When systematic uncertainties are included, results are 7.5% and 2% respectively. In addition we find that the values of most of the other inputs of the standard solar models required to better describe the helioseismic and neutrino data are in good agreement with those adopted as the standard priors, with the exception of the astrophysical factor $S_{11}$ and the microscopic diffusion rates, for which data suggests a 1% and 30% reduction respectively. As an output of the study we derive the corresponding data driven predictions for the solar neutrino fluxes.

We compute a new generation of standard solar models (SSMs) that includes recent updates on some important nuclear reaction rates and a more consistent treatment of the equation of state. Models also include a novel and flexible treatment of opacity uncertainties based on opacity kernels, required in the light of recent theoretical and experimental works on radiative opacity. Two large sets of SSMs, each based on a different canonical set of solar abundances with high and low metallicity (Z), are computed to determine model uncertainties and correlations among different observables. We present detailed comparisons of high- and low-Z models against different ensembles of solar observables including solar neutrinos, surface helium abundance, depth of convective envelope and sound speed profile. A global comparison, including all observables, yields a p-value of 2.7$\sigma$ for the high-Z model and 4.7$\sigma$ for the low-Z one. When the sound-speed differences in the narrow region of $0.65 < r/R_{sun} < 0.70$ are excluded from the analysis, results are 0.9$\sigma$ and 3.0$\sigma$ for high- and low-Z models respectively. These results show that: high-Z models agree well with solar data but have a systematic problem right below the bottom of the convective envelope linked to steepness of molecular weight and temperature gradients, and that low-Z models lead to a much more general disagreement with solar data. We also show that, while simple parametrizations of opacity uncertainties can strongly alleviate the solar abundance problem, they are insufficient to substantially improve the agreement of SSMs with helioseismic data beyond that obtained for high-Z models due to the intrinsic correlations of theoretical predictions. All the detailed information of the new solar models, can be found online: http://www.ice.cat/personal/aldos/Solar_Data.html

Because of the large neutron excess of $^{22}$Ne, this isotope rapidly sediments in the interior of the white dwarfs. This process releases an additional amount of energy, thus delaying the cooling times of the white dwarf. This influences the ages of different stellar populations derived using white dwarf cosmochronology. Furthermore, the overabundance of $^{22}$Ne in the inner regions of the star, modifies the Brunt-V\"ais\"al\"a frequency, thus altering the pulsational properties of these stars. In this work, we discuss the impact of $^{22}$Ne sedimentation in white dwarfs resulting from Solar metallicity progenitors ($Z=0.02$). We performed evolutionary calculations of white dwarfs of masses $0.528$, $0.576$, $0.657$ and $0.833$ M$_{\sun}$, derived from full evolutionary computations of their progenitor stars, starting at the Zero Age Main Sequence all the way through central hydrogen and helium burning, thermally-pulsing AGB and post-AGB phases. Our computations show that at low luminosities ($\log(L/L_{\sun})\la -4.25$), $^{22}$Ne sedimentation delays the cooling of white dwarfs with Solar metallicity progenitors by about 1~Gyr. Additionally, we studied the consequences of $^{22}$Ne sedimentation on the pulsational properties of ZZ~Ceti white dwarfs. We find that $^{22}$Ne sedimentation induces differences in the periods of these stars larger than the present observational uncertainties, particularly in more massive white dwarfs.

We study the impact on the Sun of an exotic energy-loss channel caused by plasmon decay into fermionic minicharged particles with charge $\epsilon e$ and mass $m_f$. We compare solar models with this extra emission to helioseismological and neutrino data, obtaining a bound $\epsilon <2.2\times 10^{-14}\ ({\rm 95\%\, CL})$ for $m_f\lesssim 25 \, \text{eV}$. Our result is comparable to previous limits from the cooling of globular cluster stars, while at the same time it is better understood and takes theoretical and observational errors into account.

We present a new statistical analysis that combines helioseismology (sound speed, surface helium and convective radius) and solar neutrino observations (the $^8$B and $^7$Be fluxes) to place upper limits to the properties of non standard weakly interacting particles. Our analysis includes theoretical and observational errors, accounts for tensions between input parameters of solar models and can be easily extended to include other observational constraints. We present two applications to test the method: the well studied case of axions and axion-like particles and the more novel case of low mass hidden photons. For axions we obtain an upper limit at $3\sigma$ for the axion-photon coupling constant of $g_{a\gamma}\,<\,4.1 \cdot 10^{-10} \rm{GeV^{-1}}$. For hidden photons we obtain the most restrictive upper limit available accross a wide range of masses for the product of the kinetic mixing and mass of $\chi m < 1.8 \cdot 10^{-12} \rm{eV}$ at $3\sigma$. Both cases improve the previous solar constraints based on the Standard Solar Models showing the power of using a global statistical approach.