We study non-rotating and isotropic strange quark stars in Lorentz-violating theories of gravity, and in particular in Ho\v{r}ava gravity and Einstein-{\ae}ther theory. For quark matter we adopt both linear and non-linear equations of state, corresponding to the MIT bag model and colour flavour locked state, respectively. The new structure equations describing hydrostatic equilibrium generalize the usual Tolman-Oppenheimer-Volkoff (TOV) equations of Einstein's general relativity. A dimensionless parameter $\nu$ measures the deviation from the standard TOV equations, which are recovered in the limit $\nu \rightarrow 0$. We compute the mass, the radius as well as the compactness of the stars, and we show graphically the impact of the parameter $\nu$ on the mass-to-radius profiles for different equations of state describing quark matter. The energy conditions and stability criteria are also considered, and they are all found to be fulfilled.

We report constraints on the nucleon-dark matter particle cross section using the internal luminosity of observed white dwarf stars in the globular cluster Messier 4. Our results cover the parameter space corresponding to relatively light dark matter particles, in the mass range $0.1~GeV-5~GeV$, which is known to be very difficult to be probed by direct dark matter searches. The additional luminosity coming from self-annihilations of dark matter particles captured inside the stars must not exceed the observed luminosity. Imposing that condition, we obtain for the spin independent cross section of light dark matter particles on baryons $\sigma_{N\chi}$ the upper bound: $\sigma_{N\chi} < 4 \times 10^{-41}{\rm cm^2}$.

We study an impact of asymmetric dark matter on properties of the neutron stars and their ability to reach the two solar masses limit, which allows us to present a new range of masses of dark matter particles and their fractions inside the star. Our analysis is based on the observational fact of the existence of three pulsars reaching this limit and on the theoretically predicted reduction of the neutron star maximal mass caused by the accumulation of dark matter in its interior. We also demonstrate that light dark matter particles with masses below 0.2 GeV can create an extended halo around the neutron star leading not to decrease, but to increase of its visible gravitational mass. By using recent results on the spatial distribution of dark matter in the Milky Way, we present an estimate of its fraction inside the neutron stars located in the Galaxy center. We show how the detection of a 2Msun neutron star in the most central region of the Galaxy will impose an upper constraint on the mass of dark matter particles of ~ 60 GeV. Future high precision measurements of the neutron stars maximal mass near the Galactic center, will put a more stringent constraint on the mass of the dark matter particle. This last result is particularly important to prepare ongoing, and future radio and X-ray surveys.

We study radial oscillations of non-rotating neutron stars (NSs) in four-dimensional General Relativity. The interior of the NS was modelled within a recently proposed multicomponent realistic equation of state (EoS) with the induced surface tension (IST). In particular, we considered the IST EoS with two sets of model parameters, that both reproduce all the known properties of normal nuclear matter, give a high quality description of the proton flow constraint, hadron multiplicities created in nuclear-nuclear collisions, consistent with astrophysical observations and the observational data from the NS-NS merger. We computed the 12 lowest radial oscillation modes, their frequencies and corresponding eigenfunctions, as well as the large frequency separation for six selected fiducial NSs (with different radii and masses of 1.2, 1.5 and 1.9 solar masses) of the two distinct model sets. The calculated frequencies show their continuous growth with an increase of the NS central baryon density. Moreover, we found correlations between the behaviour of first eigenfunction calculated for the fundamental mode, the adiabatic index and the speed of sound profile, which could be used to probe the internal structure of NSs with the asteroseismology data.

We investigate the impact of gravity portal to properties of white dwarfs, such as equation-of-state as well as cooling time. We find that the interaction between dark matter spin-zero bosons and electrons in the mean field approximation softens the equation-of-state, and the object coolers down slower compared to the usual case.

Dark photons are particles invoked in some extensions of the Standard Model which could account for at least part of the dark matter content of the Universe. It has been proposed that the production of dark photons in stellar interiors could happen at a rate that depends on both, the dark photon mass and its coupling to Standard Model particles (the kinetic mixing parameter $\chi$). In this work we aim at exploring the impact of dark photon productions in the stellar core of solar mass RGB stars during late evolutionary phases. We demonstrate that near the so-called RGB bump, dark photons production may be an energy sink for the star sufficiently significative to modify the extension of the star convective zones. We show that Asteroseismology is able to detect such variations in the structure, allowing us to predict an upper limit of $\rm 900\ eV$ and $5\times 10^{-15}$ for the mass and kinetic mixing of the dark photons, respectively. We also demonstrate that additional constraints can be derived from the fact that dark photons increase the luminosity of the RGB tip over the current observational uncertainties. This work thus paves the way for an empirical approach to deepen the study of such dark-matter particles.

We propose helioseismology as a new, precision probe of fifth forces at astrophysical scales, and apply it on the most general scalar-tensor theories for dark energy, known as Degenerate Higher-Order Scalar-Tensor theories (DHOST). We explain how the effect of the fifth force on the solar interior leaves an observable imprint on the acoustic oscillations, and under certain assumptions we numerically compute the non-radial pulsation eigenfrequencies within modified gravity. We illustrate its constraining power by showing that helioseismic observations have the potential to improve constraints on the strength of the fifth force by more than $2$ orders of magnitude, as $-1.8 \cdot 10^{-3} \leq Y \leq 1.2 \cdot 10^{-3}$ (at $2\sigma$). This in turn would suggest constraints of similar order for the theory's free functions around a cosmological background ($\alpha_{\text{H}}, \beta_{1}$).

In this work, we study the impact of asymmetric dark matter (ADM) on low-mass main-sequence stars in the Milky Way's nuclear star cluster, where the dark matter (DM) density is expected to be orders of magnitude above what is found near the Sun (${\rho }_{\mathrm{DM}}\gtrsim {10}^{3}\ \mathrm{GeV}\ {\mathrm{cm}}^{-3}$). Using a modified stellar evolution code and considering a DM particle ($m_{\chi} = 4 \text{ GeV}$) with a spin-dependent interaction cross section close to the limits allowed by direct detection, we found that the interactions of ADM with baryons in the star's core can have two separate effects on the evolution of these stars: a decrease in the hydrogen burning rate, extending the duration of the main-sequence of stars with $M ~ 1M_{\odot}$ by a few Gyr; the suppression of the onset of convection in the core of stars with $M \lesssim 1.5M_{\odot}$ and consequent quench of supply for the nuclear reactions. If we consider $\rho_{\text{DM}} > 10^3 \ \text{GeV cm}^{-3}$ (corresponding to the inner 5 pc of the Milky Way), stars lighter than the Sun will have a main-sequence life span comparable to the current age of the universe. Stars heavier than two solar masses are not sensitive to the DM particles considered here.

We obtain an exact analytical solution to Einstein's field equations assuming a non-linear equation-of-state and a particular mass function. Our solution describes the interior of anisotropic color flavor locked strange quark stars. All energy conditions are fulfilled, and therefore the solution obtained here is a realistic solution within General Relativity. The mass-to-radius profile is obtained, and the compactness of the star is computed.

We study the relation between the internal structures of 10 benchmark main-sequence F-stars and their rotational properties. Stellar rotation of main-sequence F-type stars can be characterised by two distinct rotational regimes. Early-type F-stars are usually rapid rotators with periods typically below 10 days, whereas later-type F-stars have longer rotation periods. Specifically, and since the two rotational regimes are tightly connected to the effective temperatures of the stars, we investigate in detail the characteristics of the partial ionisation zones in the outer convective envelopes of these stars, which in turn, depend on the internal temperature profiles. Our study shows that the two rotational regimes might be distinguished by the relative locations of the partial ionisation region of heavy elements and the base of the convective zone. Since in all these stars is expected a dynamo-driven magnetic field where the shear layer between convective and radiative zones (tachocline) plays an important role, this result suggests that the magnetic field may be related to the combined properties of convection and ionisation.

We apply a novel equation of state (EoS) that includes the surface tension contribution induced by interparticle interaction and asymmetry between neutrons and protons, to the study of neutron star (NS) properties. This elaborated EoS is obtained from the virial expansion applied to multicomponent particle mixtures with hard core repulsion. The considered model is in full concordance with all the known properties of normal nuclear matter, provides a high-quality description of the proton flow constraints, hadron multiplicities created during the nuclear-nuclear collision experiments, and equally is consistent with astrophysical data coming from NS observations. The analysis suggests that the best model parameterization gives the incompressibility factor $K_{0}$, symmetry energy $J$, and symmetry energy slope $L$ at normal nuclear density equal to $200$ MeV, $30$ MeV, and $113.28-114.91$ MeV, respectively. The mass-radius relations found for NSs computed with this EoS are consistent with astrophysical observations.

The next generation of solar neutrino detectors will provide a precision measure of the $^8$B electron-neutrino spectrum in the energy range from 1-15 MeV. Although the neutrino spectrum emitted by $^8$B $\beta$-decay reactions in the Sun's core is identical to the neutrino spectrum measured in the laboratory, due to vacuum and matter flavor oscillations, this spectrum will be very different from that measured on Earth by the different solar neutrino experiments. We study how the presence of dark matter (DM) in the Sun's core changes the shape of the $^8$B electron-neutrino spectrum. These modifications are caused by local variations of the electronic density and the $^8$B neutrino source, induced by local changes of the temperature, density and chemical composition. Particularly relevant are the shape changes at low and medium energy range $(E_\nu\le 10 {\; \rm MeV})$, for which the experimental noise level is expected to be quite small. If such a distortion in the $^8$B$\nu_e$ spectrum were to be observed, this would strongly hint in favor of the existence of DM in the Sun's core. The $^8$B electron-neutrino spectrum provides a complementary method to helioseismology and total neutrino fluxes for constraining the DM properties. In particular, we study the impact of light asymmetric DM on solar neutrino spectra. Accurate neutrino spectra measurements could help to determine whether light asymmetric DM exists in the Sun's core, since it has been recently advocated that this type of DM might resolve the solar abundance problem.

Contact interactions between sterile neutrinos and dark matter particles in a hidden sector have been suggested as a good solution to simultaneously resolve the dark matter problem and anomalies in neutrino experiments. In this non-standard particle physics model, sterile and active neutrinos change their through vacuum oscillations and matter (or Mikheyev--Smirnov--Wolfenstein) oscillations, in which the latter mechanism of flavor oscillation depends strongly on the concentration of dark matter in the Sun's core. We found that a large concentration of dark matter in the Sun's interior changes substantially the shape of ${\rm ^8B}$ and $\rm hep$ electron neutrino spectra, but has an insignificant impact on the other neutrino spectra (i.e., $\rm pp$, $\rm pep$, $\rm ^7Be$ and $\rm ^{15}O$, $\rm ^{13}N$ and $\rm ^{17}F$). The strength of the interaction of the dark matter particles with neutrinos depends on an effective coupling constant, $G_\chi$, which is an analog of the Fermi constant for the hidden sector. By using the latest $\rm ^8B$ solar neutrino flux, we found that $G_\chi$ must be smaller than $\rm 0.5\times 10^9$ $G_{\rm F}$ for this particle physics model to be in agreement with the data.

A novel equation of state used for analysis of the heavy ion collision experimental data is generalized to also describe the matter inside neutron stars. This approach differs from others by including an induced surface tension caused by the short range repulsion between particles accounted by presence of their hard cores. This new equation of state respects thermodynamic consistency and is free from the usual causality problems. This makes it applicable over a wide range of temperatures and baryon densities, including the ones inside neutron stars. Despite small number of parameters the present equation of state reproduces properties of the saturated nuclear matter and is consistent with the proton flow constraint as well as with astrophysical data obtained from neutron star observations. Accordingly, we found the sets of parameters that agree with the same ones obtained from the heavy ion collision data analysis.

We predict the sterile neutrino spectrum of some of the key solar nuclear reactions and discuss the possibility of these being observed by the next generation of solar neutrino experiments. By using an up-to-date standard solar model with good agreement with current helioseismology and solar neutrino flux data sets, we found that from solar neutrino fluxes arriving on Earth only 3\%-4\% correspond to the sterile neutrino. The most intense solar sources of sterile neutrinos are the $pp$ and $^7Be$ nuclear reactions with a total flux of $2.2\times 10^{9}\;{\rm cm^2 s^{-1}}$ and $1.8\times 10^{8}\;{\rm cm^2 s^{-1}}$, followed by the $^{13}N$ and $^{15}O$ nuclear reactions with a total flux of $1.9\times 10^{7}\;{\rm cm^2 s^{-1}}$ and $1.7\times 10^{7}\;{\rm cm^2 s^{-1}}$. Moreover, we compute the sterile neutrino spectra of the nuclear proton-proton nuclear reactions -- $pp$, $hep$ and $^8B$ and the carbon-nitrogen-oxygen -- $^{13}N$, $^{15}O$ and $^{17}F$ and the spectral lines of $^7Be$.

Extreme solar activity fluctuations and the occurrence of solar grand minima and maxima episodes, are well established, observed features of the solar cycle. Nevertheless, such extreme activity fluctuations and the dynamics of the solar cycle during Maunder minima-like episodes remain ill-understood. We explore the origin of such extreme solar activity fluctuations and the role of dual poloidal field sources, namely the Babcock-Leighton mechanism and the mean-field alpha effect in the dynamics of the solar cycle. We mainly concentrate on entry and recovery from grand minima episodes such as the Maunder minimum and the dynamics of the solar cycle. We use a kinematic solar dynamo model with a novel set-up in which stochastic perturbations force two distinct poloidal field alpha effects. We explore different regimes of operation of these poloidal sources with distinct operating thresholds, to identify the importance of each. The perturbations are implemented independently in both hemispheres which allows one to study the level of hemispheric coupling and hemispheric asymmetry in the emergence of sunspots. From the simulations performed we identify a few different ways in which the dynamo can enter a grand minima episode. While fluctuations in any of the $\alpha$ effects can trigger intermittency we find that the mean-field alpha effect is crucial for the recovery of the solar cycle from a grand minima episode which a Babcock-Leighton source alone, fails to achieve. Our simulations also demonstrate other cycle dynamics. We conclude that stochastic fluctuations in two interacting poloidal field sources working with distinct operating thresholds is a viable candidate for triggering episodes of extreme solar activity and that the mean-field alpha effect capable of working on weak, sub-equipartition fields is critical to the recovery of the solar cycle following an extended solar minimum.

We use the most recent, complete and independent measurements of masses and radii of white dwarfs in binaries to bound the class of non-trivial modified gravity theories, viable after GW170817/GRB170817, using its effect on the mass-radius relation of the stars. We show that the uncertainty in the latest data is sufficiently small that residual evolutionary effects, most notably the effect of core composition, finite temperature and envelope structure, must now accounted for if correct conclusions about the nature of gravity are to be made. We model corrections resulting from finite temperature and envelopes to a base Hamada-Salpeter cold equation of state and derive consistent bounds on the possible modifications of gravity in the stars' interiors, finding that $Y< 0.14$ at 95\% confidence, an improvement of a factor of three with respect to previous bounds. Finally, our analysis reveals some fundamental degeneracies between the theory of gravity and the precise chemical makeup of white dwarfs.

The high-precision data available from the Kepler satellite allows us to venture in the study of the complex outer convective envelopes of solar-type stars. We use a seismic diagnostic, specialized for investigating the outer layers of solar-type stars, to infer the impact of the ionization processes on the oscillation spectrum, for a sample of Kepler stars. These stars, of spectral type F, cover all of the observational seismic domain of the acoustic oscillation spectrum in solar-type stars. They also cover the range between a cool F-dwarf (\sim 6000 K) and a hotter F-star (\sim 6400 K). Our study reveals the existence of two relevant ionization regions. One of these regions, which is located closer to the surface of the star, is commonly associated with the second ionization of helium, although other chemical species also contribute to ionization. The second region, located deeper in the envelope, is linked with the ionization of heavy elements. Specifically, in this study, we analyze the elements carbon, nitrogen, oxygen, neon, and iron. Both regions can be related to the K electronic shell. We show that, while for cooler stars like the Sun, the influence of this second region on the oscillation frequencies is small; in hotter stars, its influence becomes comparable to the influence of the region of the second ionization of helium. This can guide us in the study of the outer layers of F-stars, specifically with the understanding of phenomena related to rotation and magnetic activity in these stars.

A novel equation of state with the surface tension induced by particles' interactions was generalized to describe the properties of the neutron stars (NSs). In this equation the interaction between particles occurs via the hard core repulsion by taking into account the proper volumes of particles. Recently, this model was successfully applied to the description of the properties of nuclear and hadron matter created in collisions of nucleons. The new approach is free of causality problems and is fully thermodynamically consistent, which enables us to use it for the investigation of the strongly interacting matter phase-diagram properties in a wide range of temperatures and baryon densities, including NSs. Here, we have calculated the mass-radius relations for a compact star using the Tolmann-Oppenheimer-Volkov (TOV) equation for two sets of parameters that satisfy the existing constraints. Accordingly, we found parameter values that are in good agreement with the same ones obtained from the nuclear-nuclear collision data analysis. The astrophysical constraints suggest that the hard core radius of baryons can vary between 0.425 fm and 0.476 fm.

In the present work we study the gravitational effects of condensed dark matter on strange stars. We consider self-interacting dark matter particles with properties consistent with current observational constraints, and dark matter inside the star is modelled as a Bose-Einstein condensate. We integrate numerically the Tolman-Oppenheimer-Volkoff equations in the two-fluid formalism assuming that strange stars are made of up to 4 per cent of dark matter. It is shown that for a mass of the dark matter particles in the range $50 MeV-160 MeV$ strange stars are characterized by a maximum mass and radius similar to the ones found for neutron stars.

We investigate the impact of resonant gravitational waves on quadrupole acoustic modes of Sun-like stars located nearby stellar black hole binary systems (such as GW150914 and GW151226). We find that the stimulation of the low-overtone modes by gravitational radiation can lead to sizeable photometric amplitude variations, much larger than the predictions for amplitudes driven by turbulent convection, which in turn are consistent with the photometric amplitudes observed in most Sun-like stars. For accurate stellar evolution models, using up-to-date stellar physics, we predict photometric amplitude variations of $1$ -- $10^3$ ppm for a solar mass star located at a distance between 1 au and 10 au from the black hole binary, and belonging to the same multi-star system. The observation of such a phenomenon will be within the reach of the Plato mission because telescope will observe several portions of the Milky Way, many of which are regions of high stellar density with a substantial mixed population of Sun-like stars and black hole binaries.

In the present work we study the scalar singlet as well as the Two-Higgs Doublet model inert dark matter particles impact on compact objects, and we provide the first constraints of the parameter space using neutron stars. The models discussed here are characterized by two free parameters, namely the mass $M_\chi$ of the scalar particle that plays the role of the dark matter in the Universe, and a dimensionless coupling constant $\lambda_\chi$ that determines the strength of the interaction of the dark matter particles with the Standard Model Higgs boson. By considering a typical neutron star we were able to obtain constraints on scalar dark matter depending on the DM annihilation cross section and self-interaction coupling constant. Our findings show that i) for heavy DM particles neutron stars can provide us with bounds better that the current limits from direct detection searches only when the self-annihilations of DM particles are negligible and the DM self-interaction coupling constant is very small, while ii) for light DM particles the bounds obtained here are comparable to limits from Higgs invisible decays unless the DM particles are extremely light.

Acoustic quadrupole modes of sunlike stars vibrate when perturbed by a passing gravitational wave generated somewhere in the Universe. Here, we compute the imprint of the gravitational waves on the acoustic spectrum of these stars for gravitational events occurring near the supermassive black hole located at the center of the Milky Way. We found that in most cases the impact of gravitational waves in low-order quadrupole modes is not above the current observational threshold of detectability, although this should be in the reach of the next generation of near infrared observatories and asteroseismology satellite missions. Equally, we found that it is possible to follow the end phase of the coalescence of binaries with large chirp masses, as these phenomena have a unique imprint in the spectra of sunlike stars affecting sequentially several low-order quadrupole modes. Moreover, we discuss the different imprints on the acoustic spectra of the different types of binary systems constituted either by two white dwarfs, two neutron stars, two black holes or a compact star and a massive black hole.

We have analyzed the theoretical model envelopes of eight Kepler F-stars by computing the phase shift of the acoustic waves, $\alpha(\omega)$, and its related function, $\beta(\omega)$. The latter is shown to be a powerful probe of the external stellar layers since it is particularly sensitive to the partial ionization zones located in these upper layers. We found that these theoretical envelopes can be organized into two groups, each of which is characterized by a distinct $\beta(\omega)$ shape that we show to reflect the differences related to the magnitudes of ionization processes. Since $\beta(\omega)$ can also be determined from the experimental frequencies, we compared our theoretical results with the observable $\beta(\omega)$. Using the function $\beta(\omega)$, and with the purpose of quantifying the magnitude of the ionization processes occurring in the outer layers of these stars, we define two indexes, $\Delta \beta_1$ and $\Delta \beta_2$. These indexes allow us to connect the microphysics of the interior of the star with macroscopic observable characteristics. Motivated by the distinct magnetic activity behaviors of F-stars, we studied the relation between the star's rotation period and these indexes. We found a trend, in the form of a power-law dependence, that favors the idea that ionization is acting as an underlying mechanism, which is crucial for understanding the relation between rotation and magnetism and even observational features such as the Kraft break.

Neutrinos coming from the Sun's core are now measured with a high precision, and fundamental neutrino oscillations parameters are determined with a good accuracy. In this work, we estimate the impact that a new neutrino physics model, the so-called generalized Mikheyev-Smirnov-Wolfenstein (MSW) oscillation mechanism, has on the shape of some of leading solar neutrino spectra, some of which will be partially tested by the next generation of solar neutrino experiments. In these calculations, we use a high-precision standard solar model in good agreement with helioseismology data. We found that the neutrino spectra of the different solar nuclear reactions of the proton-proton chains and carbon-nitrogen-oxygen cycle have quite distinct sensitivities to the new neutrino physics. The $HeP$ and $^8B$ neutrino spectra are the ones for which their shapes are more affected when neutrinos interact with quarks in addition to electrons. The shape of the $^{15}O$ and $^{17}F$ neutrino spectra are also modified, although in these cases the impact is much smaller. Finally, the impact in the shape of the $PP$ and $^{13}N$ neutrino spectra is practically negligible.

The high quality data provided by helioseismology, solar neutrino flux measurements, spectral determination of solar abundances, nuclear reactions rates coefficients among other experimental data, leads to the highly accurate prediction of the internal structure of the present Sun - the standard solar model. In this talk, I have discussed how the standard solar model, the best representation of the real Sun, can be used to study the properties of dark matter, for which two complementary approaches have been developed: - to limit the number of theoretical candidates proposed as the dark matter particles, this analysis complements the experimental search of dark matter, and - as a template for the study of the impact of dark matter in the evolution of stars, which possibly occurs for stellar populations formed in regions of high density of dark matter, such as stars formed in the centre of galaxies and the first generations of stars.

So far, direct detection searches have come empty handed in their quest for Dark Matter (DM). Meanwhile, asteroseismology arises as a complementary tool to study DM, as its accumulation in a star can enhance energy transport, by providing a conduction mechanism, producing significant changes in the stellar structure during the course of the star's evolution. The stellar core, particularly affected by the presence of DM, can be investigated through precise asteroseismic diagnostics. We modelled three stars including DM energy transport: the Sun, a slightly less massive and much older star, KIC 7871531 ($0.85 \, M_{\odot}$, $9.41 \, \text{Gyr}$), and a more massive and younger one, KIC 8379927 ($1.12 \, M_{\odot}$, $1.82 \, \text{Gyr}$). We considered both the case of Weakly Interactive Massive Particles, albeit with a low annihilation, and the case of Asymmetric DM for which the number of trapped particles in the star can be much greater. By analysing these models with asteroseismic separation ratios weighted towards the core, we found indications limiting the effective spin-dependent DM-proton coupling for masses of a few GeV. This independent result is very close to the most recent and most stringent direct detection DM constraints.

We use a seismic diagnostic, based on the derivative of the phase shift of the reflected by the surface acoustic waves, to probe the outer layers of the star HD 49933. This diagnostic is particularly sensitive to partial ionization processes occurring above the base of the convective zone. The regions of partial ionization of light elements, hydrogen and helium, have well known seismological signatures. In this work we detect a different seismic signature in the acoustic frequencies, that we showed to correspond to the location where the partial ionization of heavy elements occurs. The location of the corresponding acoustic glitch lies between the region of the second ionization of helium and the base of the convective zone, approximately 5% below the surface of the stars.

We used a consistent and robust solar model to obtain upper limits placed by neutrino telescopes, such as Ice- Cube and Super-Kamiokande, on the Dark Matter-nucleon scattering cross-section, for a general model of Dark Matter with a velocity dependent (p-wave) thermally averaged cross-section. In this picture, the Boltzmann equation for the Dark Matter abundance is numerically solved satisfying the Dark Matter density measured from the Cosmic Microwave Background (CMB). We show that for lower cross-sections and higher masses, the Dark Matter annihilation rate drops sharply, resulting in upper bounds on the scattering cross-section one order of magnitude above those derived from a velocity independent (s-wave) annihilation cross-section. Our results show that upper limits on the scattering cross-section obtained from Dark Matter annihilating in the Sun are sensible to the uncertainty in current standard solar models, fluctuating a maximum of 20 % depending on the annihilation channel.

Stars are massive resonators that may be used as gravitational-wave (GW) detectors with isotropic sensitivity. New insights on stellar physics are being made possible by asteroseismology, the study of stars by the observation of their natural oscillations. The continuous monitoring of oscillation modes in stars of different masses and sizes (e.g., as carried out by NASA's Kepler mission) opens the possibility of surveying the local Universe for GW radiation. Red-giant stars are of particular interest in this regard. Since the mean separation between red giants in open clusters is small (a few light years), this can in principle be used to look for the same GW imprint on the oscillation modes of different stars as a GW propagates across the cluster. Furthermore, the frequency range probed by oscillations in red giants complements the capabilities of the planned eLISA space interferometer. We propose asteroseismology of red giants as a novel approach in the search for gravitational waves.

Sun-like stellar oscillations are excited by turbulent convection and have been discovered in some 500 main sequence and sub-giant stars and in more than 12,000 red giant stars. When such stars are near gravitational wave sources, low-order quadrupole acoustic modes are also excited above the experimental threshold of detectability, and they can be observed, in principle, in the acoustic spectra of these stars. Such stars form a set of natural detectors to search for gravitational waves over a large spectral frequency range, from $10^{-7}$ Hz to $10^{-2}$ Hz. In particular, these stars can probe the $10^{-6}$ Hz -- $10^{-4}$ Hz spectral window which cannot be probed by current conventional gravitational wave detectors, such as SKA and eLISA. The PLATO stellar seismic mission will achieve photospheric velocity amplitude accuracy of $~ {\rm cm/s}$. For a gravitational wave search, we will need to achieve accuracies of the order of $10^{-2}{\rm cm/s}$, i.e., at least one generation beyond PLATO. However, we have found that multi-body stellar systems have the ideal setup for this type of gravitational wave search. This is the case for triple stellar systems formed by a compact binary and an oscillating star. Continuous monitoring of the oscillation spectra of these stars to a distance of up to a kpc could lead to the discovery of gravitational waves originating in our galaxy or even elsewhere in the universe. Moreover, unlike experimental detectors, this observational network of stars will allow us to study the progression of gravitational waves throughout space.

The sunspot activity is the end result of the cyclic destruction and regeneration of magnetic fields by the dynamo action. We propose a new method to analyze the daily sunspot areas data recorded since 1874. By computing the power spectral density of daily data series using the Mexican hat wavelet, we found a power spectrum with a well-defined shape, characterized by three features. The first term is the 22 yr solar magnetic cycle, estimated in our work to be of 18.43 yr. The second term is related to the daily volatility of sunspots. This term is most likely produced by the turbulent motions linked to the solar granulation. The last term corresponds to a periodic source associated with the solar magnetic activity, for which the maximum of power spectral density occurs at 22.67 days. This value is part of the 22-27 day periodicity region that shows an above-average intensity in the power spectra. The origin of this 22.67 day periodic process is not clearly identified, and there is a possibility that it can be produced by convective flows inside the star. The study clearly shows a north-south asymmetry. The 18.43 yr periodical source is correlated between the two hemispheres, but the 22.67 day one is not correlated. It is shown that towards the large timescales an excess occurs in the northern hemisphere, especially near the previous two periodic sources. To further investigate the 22.67 day periodicity we made a Lomb-Scargle spectral analysis. The study suggests that this periodicity is distinct from others found nearby.

During the last century, with the development of modern physics in such diverse fields as thermodynamics, statistical physics, and nuclear and particle physics, the basic principles of the evolution of stars have been successfully well understood. Nowadays, a precise diagnostic of the stellar interiors is possible with the new fields of helioseismology and astroseismology. Even the measurement of solar neutrino fluxes, once a problem in particle physics, is now a powerful probe of the core of the Sun. These tools have allowed the use of stars to test new physics, in particular the properties of the hypothetical particles that constitute the dark matter of the Universe. Here we present recent results obtained using this approach.

Current helioseismology observations allow the determination of the frequencies and surface velocity amplitudes of solar acoustic modes with exceptionally high precision. In some cases, the frequency accuracy is better than one part in a million. We show that there is a distinct possibility that the quadrupole acoustic modes of low order could be excited by gravitational waves (GWs), if the GWs have a strain amplitude in the range $10^{-20}h_{-20}$ with $h_{-20}\sim 1$ or $h_{-20}\sim 10^{3}$, as predicted by several types of GW sources, such as galactic ultracompact binaries or extreme mass ratio inspirals and coalescence of black holes. If the damping rate at low order is $ 10^{-3}\eta_N$ $\mu{\rm Hz}$, with $\eta_N\sim 10^{-3} $ - $1,$ as inferred from the theory of stellar pulsations, then GW radiation will lead to a maximum rms surface velocity amplitude of quadrupole modes of the order of $h_{-20}\eta_N^{-1}\sim$ $10^{-9}$ - $10^{-3}$ ${\rm cm\; s^{-1}}$, on the verge of what is currently detectable via helioseismology. The frequency and sensitivity range probed by helioseismological acoustic modes overlap with, and complement, the capabilities of eLISA for the brightest resolved ultracompact galactic binaries.

The detection of gravity modes produced in the solar radiative zone has been a challenge in modern astrophysics for more than 30 yr and their amplitude in the core is not yet determined. In this Letter, we develop a new strategy to look for standing gravity modes through solar neutrino fluxes. We note that due to a resonance effect, the gravity modes of low degree and low order have the largest impact on the $^{8}B$ neutrino flux. The strongest effect is expected to occur for the dipole mode with radial order $2$, corresponding to periods of about 1.5 hr. These standing gravity waves produce temperature fluctuations that are amplified by a factor of 170 in the boron neutrino flux for the corresponding period, in consonance with the gravity modes. From current neutrino observations, we determine that the maximum temperature variation due to the gravity modes in the Sun's core is smaller than $5.8\times 10^{-4}$. This study clearly shows that due to their high sensitivity to the temperature, the $^8B$ neutrino flux time series is an excellent tool to determine the properties of gravity modes in the solar core. Moreover, if gravity mode footprints are discovered in the $^{8}B$ neutrino flux, this opens a new line of research to probe the physics of the solar core as non-standing gravity waves of higher periods cannot be directly detected by helioseismology but could leave their signature on boron neutrino or on other neutrino fluxes.

Assuming the existence of a primordial asymmetry in the dark sector, we study how DM-baryon long-range interactions, induced by the kinetic mixing of a new $U(1)$ gauge boson and the photon, affects the evolution of the Sun and in turn the sound speed profile obtained from helioseismology. Thanks to the explicit dependence on the exchanged momenta in the differential cross section (Rutherford-like scattering), we find that dark matter particles with a mass of $\sim 10\;{\rm GeV}$, kinetic mixing parameter of the order of $10^{-9}$ and a mediator with a mass smaller than a few MeV improve the agreement between the best solar model and the helioseismic data without being excluded by direct detection experiments. In particular, the \LUX\ detector will soon be able to either constrain or confirm our best fit solar model in the presence of a dark sector with long-range interactions that reconcile helioseismology with thermal neutrino results.

This article reviews some of the leading results obtained in solar dynamo physics by using temporal oscillator models as a tool to interpret observational data and dynamo model predictions. We discuss how solar observational data such as the sunspot number is used to infer the leading quantities responsible for the solar variability during the last few centuries. Moreover, we discuss the advantages and difficulties of using inversion methods (or backward methods) over forward methods to interpret the solar dynamo data. We argue that this approach could help us to have a better insight about the leading physical processes responsible for solar dynamo, in a similar manner as helioseismology has helped to achieve a better insight on the thermodynamic structure and flow dynamics in the Sun's interior.

Dark matter particles with properties identical to dark matter candidates that are hinted at by several international collaborations dedicated to experimental detection of dark matter (DAMA, COGENT, CRESST and CDMS-II, although not, most notably, by LUX), and which also have a dark matter asymmetry identical to the observed baryon asymmetry (Planck and Wilkinson Microwave Anisotropy Probe), may produce a significant impact on the evolution of the first generation of low-metallicity stars. The lifetimes of these stars in different phases of stellar evolution are significantly extended, namely, in the pre-main sequence, main sequence, and red giant phases. In particular, intermediate-mass stars in the red giant phase experience significant changes in their luminosity and chemical composition. The annihilations of dark matter particles affect the interior of the star in such a way that the $3\alpha-$reaction becomes less efficient in the production of carbon and oxygen. This dark matter effect contradicts the excess of carbon and other metals observed today in stars of low mass and low metallicity. Hence, we can impose an upper limit on the dark matter halo density, and therefore on the redshift, at which the first generation of low-metallicity stars formed.

We investigate the effects of a magnetic dipole moment of asymmetric dark matter (DM) in the evolution of the Sun. The dipole interaction can lead to a sizable DM scattering cross section even for light DM, and the asymmetric dark matter can lead to a large DM number density in the Sun. We find that solar model precision tests, using as diagnostic the sound speed profile obtained from helioseismology data, exclude dipolar dark matter particles with a mass larger than 4.3 GeV and magnetic dipole moment larger than 1.6 $\times \;10^{-17}$ e cm.

In the last decade the quality and the amount of observational asteroseismic data that has been made available by space based missions had a tremendous upgrowth. The determination of asteroseismic parameters to estimate the fundamental physical processes occurring in stars' interiors, can be done today in a way that has never been possible before. In this work we propose to compute the seismic observable $\beta$, which is a proxy of the phase shift of the acoustic modes propagating in the envelope of the Sun-like stars. This seismic parameter $\beta$ can be used to identify rapid variation regions usually known as glitches. We show that a small variation in the structure, smaller than 1% in the sound speed, produces a glitch in the acoustic potential that could explain the oscillatory character of $\beta$. This method allows us to determine the location and the thickness of the glitch with precision. We applied this idea to the Sun-like star $\alpha$ Centauri A and found a glitch located at approximately $0.94\,R$ (1400 s) and with a thickness of 0.2% of the stars' radius. This is fully consistent with the data and also validates other seismic tests.

The impact of global acoustic modes on the 8B neutrino flux time series is computed for the first time. It is shown that the time fluctuations of the 8B neutrino flux depend on the amplitude of acoustic eigenfunctions in the region where the 8B neutrino flux is produced: modes with low n (or order) that have eigenfunctions with a relatively large amplitude in the Sun's core, strongly affect the neutrino flux; conversely, modes with high n that have eigenfunctions with a minimal amplitude in the Sun's core have a very small impact on the neutrino flux. It was found that the global modes with a larger impact on the 8B neutrino flux have a frequency of oscillation in the interval 250 \mu Hz to 500 \mu Hz (or a period in the interval 30 minutes to 70 minutes), such as the f-modes (n=0) for the low degrees, radial modes of order n smaller or equal to 3, and the dipole mode of order n=1. Their corresponding neutrino eigenfunctions are very sensitive to the solar inner core and are unaffected by the variability of the external layers of the solar surface. If time variability of neutrinos is observed for these modes, it will lead to new ways of improving the sound speed profile inversion in the central region of the Sun.

Neutrinos are produced in several neutrino nuclear reactions of the proton-proton chain and carbon-nitrogen-oxygen cycle that take place at different radius of the Sun's core. Hence, measurements of solar neutrino fluxes provide a precise determination of the local temperature. The accumulation of non-annihilating light dark matter particles (with masses between 5 GeV and 16 GeV in the Sun produces a change in the local solar structure, namely, a decrease in the central temperature of a few percent. This variation depends on the properties of the dark matter particles, such as the mass of the particle and its spin-independent scattering cross-section on baryon-nuclei, specifically, the scattering with helium, oxygen, and nitrogen among other heavy elements. This temperature effect can be measured in almost all solar neutrino fluxes. In particular, by comparing the neutrino fluxes generated by stellar models with current observations, namely 8B neutrino fluxes, we find that non-annihilating dark matter particles with a mass smaller than 10 GeV and a spin-independent scattering cross-section with heavy baryon-nuclei larger than 3 x 10^{-37} cm^-2 produce a variation in the 8B neutrino fluxes that would be in conflict with current measurements.

Recent observations of solar twin stars with planetary systems like the Sun, have uncovered that these present a peculiar surface chemical composition. This is believed to be related to the formation of earth-like planets. This suggests that twin stars have a radiative interior that is richer in heavy elements than their envelopes. Moreover, the current standard solar model does not fully agree with the helioseismology data and solar neutrino flux measurements. In this work, we find that this agreement can improve if the Sun has mass loss during the pre-main sequence, as was previously shown by other groups. Despite this better agreement, the internal composition of the Sun is still uncertain, especially for elements heavier than helium. With the goal of inferring the chemical abundance of the solar interior, we tested several chemical compositions. We found that heavy element abundances influence the sound speed and solar neutrinos equally. Nevertheless, the carbon-nitrogen-oxygen (CNO;13N, 15O and 17F) neutrino fluxes are the most affected; this is due to the fact that contrarily to proton-proton (pp, pep, 8B and 7Be) neutrino fluxes, the CNO neutrino fluxes are less dependent on the total luminosity of the star. Furthermore, if the central solar metallicity increases by 30%, as hinted by the solar twin stars observations, this new solar model predicts that 13N, 15O and 17F neutrino fluxes increase by 25%-80% relative to the standard solar model. Finally, we highlight that the next generation of solar neutrino experiments will not only put constraints on the abundances of carbon, oxygen and nitrogen, but will also give some information about their radial distribution.

The electronic density in the Sun's inner core is inferred from the 8B, 7Be and pep neutrino flux measurements of the Super-Kamiokande, SNO and Borexino experiments. We have developed a new method in which we use the KamLAND detector determinations of the neutrino fundamental oscillation parameters: the mass difference and the vacuum oscillation angle. Our results suggest that the solar electronic density in the Sun's inner core (for a radius smaller than 10% of the solar radius) is well above the current prediction of the standard solar model, and by as much as 25%. A potential confirmation of these preliminary findings can be achieved when neutrino detectors are able to reduce the error of the electron-neutrino survival probability by a factor of 15.

We report recent results on the impact of asymmetric dark matter (DM) particles on low-mass stars. First, we found that the small convective core expected in stars with masses between 1.1 and 1.3 Msun is suppressed due to DM cooling. Moreover, stars with masses below 1 Msun have their central temperatures and densities more strongly influenced by DM than in the solar case. We were able to put limits to the DM mass and spin-dependent DM-proton scattering cross section by comparing the modelling of the nearby star Alpha Cen B with photometric, spectroscopic and asteroseismic observations.

Solar neutrinos coming from different nuclear reactions are now detected with a high statistics. Consequently, an accurate spectroscopic analysis of the neutrino fluxes arriving on the Earth's detectors become available, in the context of neutrino oscillations. In this work, we explore the possibility of using this information to infer the radial profile of the electronic density in the solar core. So, we discuss the constraints on the Sun's density and chemical composition that can be determined from solar neutrino observations. This approach constitutes an independent and alternative diagnostic to the helioseismic investigations already done. The direct inversion method, that we propose to get the radial solar electronic density profile, is almost independent of the solar model.

We report the first constraints on the properties of weakly interacting low-mass dark matter (DM) particles using asteroseismology. The additional energy transport mechanism due to accumulated asymmetric DM particles modifies the central temperature and density of low-mass stars and suppresses the convective core expected in 1.1-1.3 Ms stars even for an environmental DM density as low as the expected in the solar neighborhood. An asteroseismic modeling of the stars KIC 8006161, HD 52265 and Alpha Cen B revealed small frequency separations significantly deviated from the observations, leading to the exclusion of a region of the DM parameter space mass versus spin-dependent DM-proton scattering cross section comparable with present experimental constraints.

The dark matter content of the Universe is likely to be a mixture of matter and antimatter, perhaps comparable to the measured asymmetric mixture of baryons and antibaryons. During the early stages of the Universe, the dark matter particles are produced in a process similar to baryogenesis, and dark matter freeze-out depends on the dark matter asymmetry and the annihilation cross section (s-wave and p-wave annihilation channels). In these \eta-parametrised asymmetric dark matter models (\eta-ADM), the dark matter particles have an annihilation cross section close to the weak interaction cross section, and a value of \eta-dark matter asymmetry close to the baryon asymmetry \eta_B. Furthermore, we assume that dark matter scattering of baryons, namely, the spin-independent scattering cross section, is of the same order as the range of values suggested by several theoretical particle physics models used to explain the current unexplained events reported in the DAMA/LIBRA, CoGeNT and CRESST experiments. Here, we constrain \eta-ADM models by investigating the impact of such a type of dark matter on the evolution of the Sun, namely, the flux of solar neutrinos and helioseismology. We find that dark matter particles with a mass smaller than 15 GeV, a spin-independent scattering cross section on baryons of the order of a picobarn, and an \eta-dark matter asymmetry with a value in the interval 10^{-12}-10^{-10}, would induce a change in solar neutrino fluxes in disagreement with current neutrino flux measurements. A natural consequence of this model is suppressed annihilation, thereby reducing the tension between indirect and direct dark matter detection experiments, but the model also allows a greatly enhanced annihilation cross section. All the cosmological \eta-asymmetric dark matter scenarios that we discuss are in agreement with the current WMAP measured values.

We develop and subsequently explore the solution space of a simple flux transport dynamo model that incorporates a time dependent large scale meridional circulation. Based on recent observations we prescribed an analytical form for the amplitude of this circulation and study its impact in the evolution of the magnetic field. We find that cyclic variations in the amplitude and frequency of the meridional flow affect the strength of the solar cycle. Variations in the amplitude of the fluctuations influence the shape of the solar cycle but are only relevant to the cycle's strength variations when they occur at a frequency different from or out of phase of the solar cycle's.

We propose a new approach to test possible corrections to Newtonian gravity using solar physics. The high accuracy of current solar models and new precise observations allow one to constrain corrections to standard gravity at unprecedented levels. Our case study is Eddington-inspired gravity, an attractive modified theory of gravity which results in non-singular cosmology and collapse. The theory is equivalent to standard gravity in vacuum, but it sensibly differs from it within matter, for instance it affects the evolution and the equilibrium structure of the Sun, giving different core temperature profiles, deviations in the observed acoustic modes and in solar neutrino fluxes. Comparing the predictions from a modified solar model with observations, we constrain the coupling parameter of the theory, |kappa_g| < 3x10^5 m^5 s^2 / kg. Our results show that the Sun can be used to efficiently constraint alternative theories of gravity.

In order to characterize how dark matter (DM) annihilation inside stars changes the aspect of a stellar cluster we computed the evolution until the ignition of the He burning of stars from 0.7 to 3.5 M_sun within halos of DM with different characteristics. We found that, when a cluster is surrounded by a dense DM halo, the positions of the cluster' stars in the H-R diagram have a brighter and hotter turn-off point than in the classical scenario without DM, therefore giving the cluster a younger appearance. The high DM densities required to produce these effects are expected only in very specific locations, such as near the center of our Galaxy. In particular, if DM is formed by the 8 GeV WIMPs recently invoked to reconcile the results from direct detection experiments, then this signature is predicted for halos of DM with a density rho_DM = 3 \cdot 10^5 GeV cm^3. A DM density gradient inside the stellar cluster would result in a broader main sequence, turn-off and red giant branch regions. Moreover, we found that for very high DM halo densities the bottom of the isochrones in the H-R diagram rises to higher luminosities, leading to a characteristic signature on the stellar cluster. We argue that this signature could be used to indirectly probe the presence of DM particles in the location of a cluster.

We studied the rate at which stars capture dark matter (DM) particles, considering different assumptions regarding the DM characteristics and in particular investigating how the stellar physics influences the capture rate. Two scenarios were considered: first, we assumed the maximal values for the spin-dependent and spin-independent DM particle-nucleon scattering cross sections allowed by the limits from direct detection experiments. Second, we considered that both scattering cross sections are of the same order, with the aim of studying the dependencies of the capture rate on stellar elements other than hydrogen. We found that the characteristics of the capture rate are very different in the two scenarios. Furthermore, we quantified the uncertainties on the computed capture rate (C_x) and on the ratio between the luminosities from DM annihilations and thermonuclear reactions (L_x/L_nuc) derived from an imprecise knowledge of the stellar structure and DM parameters. For instance, while an uncertainty of 10% on the typical DM velocity leads to similar errors on the computed C_x and L_x/L_nuc, the same uncertainty on the stellar mass becomes more relevant and duplicates the errors. Our results may be used to evaluate the reliability of the computed capture rate for the hypothetical use of stars other than the Sun as DM probes.

We studied the formation and evolution of low-mass stars within halos with high concentration of dark matter (DM) particles, using a highly sophisticated expression to calculate the rate at which DM particles are captured inside the star. For very high DM densities in the host halo (\rho_{\chi}>10^10 GeV cm^-3 for a 1 M_{\odot} star), we found that young stars stop sooner their gravitational collapse in the pre-Main Sequence phase, reaching states of equilibrium in which DM annihilation is their only source of energy. The lower effective temperature of these stars, which depends on the properties of the DM particles and DM halo, may be used as an alternative method to investigate the nature of DM.

The annihilation of huge quantities of captured dark matter (DM) particles inside low-mass stars has been shown to change some of the stellar properties, such as the star's effective temperature or the way the energy is transported throughout the star. While in the classical picture, without DM, a star of 1 M_sun is expected to have a radiative interior during the main sequence, the same star evolving in a halo of DM with a density rho_x > 10^8 GeV cm^-3 will develop a convective core in order to evacuate the energy from DM annihilation in a more efficient way. This convective core leaves a discontinuity in the density and sound-speed profiles that can be detected by the analysis of the stellar oscillations. In this paper we present an approach towards the use of asteroseismology to detect the signature produced by the presence of DM inside a star, and we propose a new methodology to infer the properties of a DM halo from the stellar oscillations (such as the product of the DM density and the DM particle-nucleon scattering cross-section).

Although helioseismology has been used as an effective tool for studying the physical mechanisms acting in most of the solar interior, the microscopic and dynamics of the deep core is still not well understood. Helioseismological anomalies may be partially resolved if the Sun captures light, non-annihilating dark matter particles, a currently discussed dark matter candidate that is motivated by recent direct detection limits. Once trapped, such particles (4-10 GeV) naturally fill the solar core. With the use of a well-defined stellar evolution code that takes into account an accurate description of the capture of dark matter particles by the Sun, we investigate the impact of such particles in its inner core. Even a relatively small amount of dark matter particles in the solar core will leave an imprint on the absolute frequency values of gravity modes, as well as the equidistant spacing between modes of the same degree. The period separation for gravity modes could reveal changes of up to 3% for annihilating dark matter and of up to 20% for non-annihilating dark matter. This effect is most pronounced in the case of the gravity dipole (l=1) modes.

In this work we use a low order dynamo model and study under which conditions can it reproduce solar grand minima. We begin by building the phase space of a proxy for the toroidal component of the solar magnetic field and we develop a model, derived from mean field dynamo theory, that gives the time evolution of the toroidal field. This model is characterized by a non-linear oscillator whose coefficients retain most of the physics behind dynamo theory. In the derivation of the model we also include stochastic oscillations in the $\alpha$ effect. We found no evidences that stochastic fluctuations in a linear $\alpha$ effect can trigger grand minima episodes in this model. In contrast, the model used points out that possible mechanism that can trigger grand minima should involve the meridional circulation, magnetic diffusivity or field intensification by buoyancy driven instabilities.

The formation and evolution of low-mass stars within dense halos of dark matter (DM) leads to evolution scenarios quite different from the classical stellar evolution. As a result of our detailed numerical work, we describe these new scenarios for a range of DM densities on the host halo, a range of scattering cross sections of the DM particles considered, and for stellar masses from 0.7 to 3 M_{\odot}. For the first time, we also computed the evolution of young low-mass stars in their Hayashi track in the pre-main sequence phase and found that, for high DM densities, these stars stop their gravitational collapse before reaching the main sequence, in agreement with similar studies on first stars. Such stars remain indefinitely in an equilibrium state with lower effective temperatures (|\Delta T_eff| > 10^3 K for a star of one solar mass), the annihilation of captured DM particles in their core being the only source of energy. In the case of lower DM densities, these proto-stars continue their collapse and progress through the main sequence burning hydrogen at a lower rate. A star of 1 M_{\odot} will spend a time greater than the current age of the universe consuming all the hydrogen in its core if it evolves in a halo with DM density \rho_{\chi}=10^9 GeV cm^-3. We also show the strong dependence of the effective temperature and luminosity of these stars on the characteristics of the DM particles and how this can be used as an alternative method for DM research.

The solar magnetism is no more considered as a purely superficial phenomenon. The SoHO community has shown that the length of the solar cycle depends on the transition region between radiation and convection. Nevertheless, the internal solar (stellar) magnetism stays poorly known. Starting in 2008, the American instrument HMI/SDO and the European microsatellite PICARD will enrich our view of the Sun-Earth relationship. Thus obtaining a complete MHD solar picture is a clear objective for the next decades and it requires complementary observations of the dynamics of the radiative zone. For that ambitious goal, space prototypes are being developed to improve gravity mode detection. The Sun is unique to progress on the topology of deep internal magnetic fields and to understand the complex mechanisms which provoke photospheric and coronal magnetic changes and possible longer cycles important for human life. We propose the following roadmap in Europe to contribute to this "impressive" revolution in Astronomy and in our Sun-Earth relationship: SoHO (1995-2007), PICARD (2008-2010), DynaMICS (2009-2017) in parallel to SDO (2008-2017) then a world-class mission located at the L1 orbit or above the solar poles.

The experimental uncertainties in G between different experiments have important implications for helioseismology. We show that these uncertainties for the standard solar model lead to a range in the value of the square of the sound speed in the nuclear region that is as much as 0.15% higher than the inverted helioseismic sound speed. While a lower value of G is preferred for the standard model, any definite prediction is masked by the uncertainties in the solar models available in the literature. However future refinements of helioseismology with an accuracy of the order of 1E-3 to 1E-4 in the square of the sound speed, especially in combination with precision measurements of the 8B solar neutrino flux should be capable of independently testing these experimental values of G.

If the Milky Way is populated by WIMPs as predicted by cosmological models of the large-scale structure of the universe and as motivated by SUSY, the capture of high-mass WIMPs by the Sun would affect the temperature, density and chemical composition of the solar core. We use the sound speed and the density profiles inferred from the helioseismic instruments on the Solar and Heliospheric Observatory (SOHO) to discuss the effect of WIMP accretion and annihilation on the evolution of the Sun. The WIMP transport of energy inside the Sun is not negligible for WIMPs with a mass smaller than 60 GeV and annihilating WIMPs with <sigma_a v>~ 10^{-27}cm^3/sec. WIMP-accreting models with WIMP masses smaller than 30 GeV are in conflict with the most recent seismic data. We combine our new constraints with the analysis of predicted neutrino fluxes from annihilating WIMPs in the solar core. Working in the framework of the Minimal Supersymmetric Standard Model and considering the neutralino as the best dark matter particle candidate, we find that supersymmetric models, consistent with solar seismic data and with recent measurements of dark matter relic density, lead to a measured muon flux on Earth in the range of 1 to 10^4 km^{-2} yr^{-1}, for neutralino masses between 30 and 400 GeV.

SNO measurements strongly constrain the central temperature of the Sun within a precision of much less than 1%. This result can be used to constrain the parameter space of SUSY dark matter particle candidates. In this first analysis we find a lower limit for the WIMP mass of 60 GeV, well above the WIMP evaporation limit of 10 GeV. Furthermore, in the event that WIMPs create a quasi-isothermal core within the Sun, they will produce a peculiar distribution of the solar neutrino fluxes measured on Earth. Typically, a WIMP with a mass of 100 GeV and annihilation cross-section of 1E-34 cm^3/sec, will decrease the standard solar model neutrino predictions, by up to 4% for the Cl, by 3% for the heavy water, and by 1% for the Ga detectors.

The presence of dark matter in the solar neighbourhood can be tested in the framework of stellar evolution theory by using the new results of helioseismology. If weakly interacting massive particles accumulate in the center of the Sun, they can provide an additional mechanism for transferring energy from the solar core. The presence of these particles produces a change in the local luminosity of the Sun of the order of 0.1%, which is now within the reach of seismic solar experiments. We find that typical WIMPs with a mass of 60 GeV, annihilation cross-section of 10^{-32} cm^3/sec, and scalar scattering interaction between 10^{-36} cm^2 and 10^{-41} cm^2, the lower end of this range being that suggested by the DAMA direct detection experiment, modify the radial profile of the square of the sound speed in the solar core by \~ 0.1 %. This level of change is strongly constrained by the most recent helioseismological results. Significantly more massive WIMPs have a much smaller effect on the solar core temperature profile. However future helioseismological experiments have the potential exploring most, if not all, of the WIMP parameter space accessible to current and planned direct detection experiments.