The Neutron Star Interior Composition Explorer (NICER) recently measured the mass and equatorial radius of the isolated neutron star PSR J0030+0451. We use these measurements to infer the moment of inertia, the quadrupole moment, and the surface eccentricity of an isolated neutron star for the first time, using relations between these quantities that are insensitive to the unknown equation of state of supranuclear matter. We also use these results to forecast the moment of inertia of neutron star A in the double pulsar binary J0737-3039, a quantity anticipated to be directly measured in the coming decade with radio observations. Combining this information with the measurement of the tidal Love number with LIGO/Virgo observations, we propose and implement the first theory-agnostic and equation-of-state-insensitive test of general relativity. Specializing these constraints to a particular modified theory, we find that consistency with general relativity places the most stringent constraint on gravitational parity violation to date, surpassing all other previously reported bounds by seven orders of magnitude.

Gravitational-wave astronomy offers not only new vistas into the realm of astrophysics, but it also opens an avenue for probing, for the first time, general relativity in its strong-field, nonlinear, and dynamical regime, where the theory's predictions manifest themselves in their full glory. We present a study of whether the gravitational wave events detected so far by the LIGO-Virgo scientific collaborations can be used to probe higher-curvature corrections to general relativity. In particular, we focus on two examples: Einstein-dilaton-Gauss-Bonnet and dynamical Chern-Simons gravity. We find that the two events with a low-mass $m \approx 7 M_{\odot}$ BH (GW151226 and GW170608) place stringent constraints on Einstein-dilaton-Gauss-Bonnet gravity, ${\alpha}^{1/2}_{\rm EdGB} \lesssim 5.6$ km, whereas dynamical Chern-Simons gravity remains unconstrained by the gravitational-wave observations analyzed.

The canonical scalar-tensor theory model which exhibits spontaneous scalarization in the strong-gravity regime of neutron stars has long been known to predict a cosmological evolution for the scalar field which generically results in severe violations of present-day Solar System constraints on deviations from general relativity. We study if this tension can be alleviated by generalizing this model to include a disformal coupling between the scalar field $\varphi$ and matter, where the Jordan frame metric ${\tilde g}_{\mu\nu}$ is related to the Einstein frame one $g_{\mu\nu}$ by ${\tilde g}_{\mu\nu}=A(\varphi)^2 (g_{\mu\nu}+\Lambda\, \partial_\mu \varphi \, \partial_\nu\varphi)$. We find that this broader theory admits a late-time attractor mechanism towards general relativity. However, the existence of this attractor requires a value of disformal scale of the order $\Lambda\gtrsim H_0^{-2}$, where $H_0$ is the Hubble parameter of today, which is much larger than the scale relevant for spontaneous scalarization of neutron stars $\Lambda \sim R_s^{2}$ with $R_s (\sim 10^{-22} H_0^{-1})$ being the typical radius of these stars. The large values of $\Lambda$ necessary for the attractor mechanism (i) suppress spontaneous scalarization altogether inside neutron stars and (ii) induce ghost instabilities on scalar field fluctuations, thus preventing a resolution of the tension. We argue that the problem arises because our disformal coupling involves a dimensionful parameter.

We connect two formalisms that describe deformations away from general relativity, one valid in the strong-field regime of neutrons stars and another valid in the radiative regime of gravitational waves: the post-Tolman-Oppenheimer-Volkoff and the parametrized-post-Einsteinian formalisms respectively. We find that post-Tolman-Oppenheimer-Volkoff deformations of the exterior metric of an isolated neutron star induce deformations in the orbital binding energy of a neutron star binary. Such a modification to the binding energy then percolates into the gravitational waves emitted by such a binary, with the leading-order post-Tolman-Oppenheimer-Volkoff modifications introducing a second post-Newtonian order correction to the gravitational wave phase. The lack of support in gravitational wave data for general relativity deformations at this post-Newtonian order can then be used to place constraints post-Tolman-Oppenheimer-Volkoff parameters. As an application, we use the binary neutron star merger event GW170817 to place the constraint $-2.4 \leq \chi \leq 44$ (at 90% credibility) on a combination of post-Tolman-Oppenheimer-Volkoff parameters. We also explore the implications of this result to the possible deformations of the mass-radius relation of neutron stars allowed within this formalism. This work opens the path towards theory-independent tests of gravity, combining astronomical observations of neutron stars and gravitational wave observations.

Much of our physical intuition about black hole quasinormal modes in general relativity comes from the eikonal/geometric optics approximation. According to the well-established eikonal model, the fundamental quasinormal mode represents wavepackets orbiting in the vicinity of the black hole's geodesic photon ring, slowly peeling off towards the event horizon and infinity. Besides its strength as a "visualisation" tool, the eikonal approximation also provides a simple quantitative method for calculating the mode frequency, in close agreement with rigorous numerical results. In this paper we move away from Einstein's theory and its garden-variety black holes and go on to consider spherically symmetric black holes in modified theories of gravity through the lens of the eikonal approximation. The quasinormal modes of such black holes are typically described by a set of coupled wave equations for the various field degrees of freedom. Considering a general, theory-agnostic, system of two equations for two perturbed fields, we derive eikonal formulae for the complex fundamental quasinormal mode frequency. In addition we show that the eikonal modes can be related to the extremum of an effective potential and its associated "photon ring". As an application of our results we consider a specific example of a modified theory of gravity with known black hole quasinormal modes and find that these are well approximated by the eikonal formulae.

An important question when studying the light curves produced by hot spots on the surface of rotating neutron stars is, how might the shape of the light curve be affected by relaxing some of the simplifying assumptions that one would adopt on a first treatment of the problem, such as that of a point-like spot. In this work we explore the dependence of light curves on the size and shape of a single hot spot on the surface of slowly-rotating neutron stars. More specifically, we consider two different shapes for the hot spots (circular and annular) and examine the resulting light curves as functions of the opening angle of the hot spot (for both) and width (for the latter). We find that the point-like approximation can describe light curves reasonably well, if the opening angle of the hot spot is less than $\sim 5^\circ$. Furthermore, we find that light curves from annular spots are almost the same as those of the full circular spot if the opening angle is less than $\sim 35^\circ$ independently of the hot spot's width.

Scalar-tensor theories of gravity where a new scalar degree of freedom couples to the Gauss-Bonnet invariant can exhibit the phenomenon of spontaneous black hole scalarization. These theories admit both the classic black holes predicted by general relativity as well as novel hairy black hole solutions. The stability of the hairy black holes is strongly dependent on the precise form of the scalar-gravity coupling. A radial stability investigation revealed that all scalarized black hole solutions are unstable when the coupling between the scalar field and the Gauss-Bonnet invariant is quadratic in the scalar, whereas stable solutions exist for exponential couplings. Here we elucidate this behavior. We demonstrate that, while the quadratic term controls the onset of the tachyonic instability that gives rise to the black hole hair, the higher-order coupling terms control the nonlinearities that quench that instability, and hence also control the stability of the hairy black hole solutions.

It has recently been shown that nontrivial couplings between a scalar and the Gauss-Bonnet invariant can give rise to black hole spontaneous scalarization. Theories that exhibit this phenomenon are among the leading candidates for testing gravity with upcoming black hole observations. All models considered so far have focused on specific forms for the coupling, neglecting scalar self-interactions. In this work, we take the first steps towards placing this phenomenon on a more robust theoretical footing by considering the leading-order scalar self-interactions as well as the scalar-Gauss-Bonnet coupling. Our approach is consistent with the principles of effective field theory and yields the simplest and most natural model. We find that a mass term for the scalar alters the threshold for the onset of scalarization, and we study the mass range over which scalarized black hole solutions exist. We also demonstrate that the quartic self-coupling is sufficient to produce scalarized solutions that are stable against radial perturbations, without the need to resort to higher-order terms in the Gauss-Bonnet coupling function. Our model therefore represents a canonical model that can be studied further, with the ultimate aim of developing falsifiable tests of black hole scalarization.

The observation of the x-ray pulse profile emitted by hotspots on the surface of neutron stars offers a unique tool to measure the bulk properties of these objects, including their masses and radii. The x-ray emission takes place at the star's surface, where the gravitational field is strong, making these observations an incise probe to examine the curvature of spacetime generated by these stars. Motivated by this and the upcoming data releases by x-ray missions, such as NICER (Neutron star Interior Composition Explorer), we present a complete toolkit to model pulse profiles of rotating neutron stars in scalar-tensor gravity. We find that in this class of theories the presence of the scalar field affects the pulse profile's overall shape, producing strong deviations from the General Relativity expectation. This finding opens the possibility of potentially using x-ray pulse profile data to obtain new constraints on scalar-tensor gravity, if the pulse profile is found to be in agreement with General Relativity.

Physical limits on the equation-of-state (EoS) parameter of a dark energy component non-minimally coupled with the dark matter field are examined in light of the second law of thermodynamics and the positiveness of entropy. Such constraints are combined with observational data sets of type Ia supernovae, baryon acoustic oscillations and the angular acoustic scale of the cosmic microwave background to impose restrictions on the behaviour of the dark matter/dark energy interaction. Considering two EoS parameterisations of the type $w = w_0 + w_a\zeta(z)$, we derive a general expression for the evolution of the dark energy density and show that the combination of thermodynamic limits and observational data provide tight bounds on the $w_0 - w_a$ parameter space.

We infer the mass distribution of neutron stars in binary systems using a flexible Gaussian mixture model and use Bayesian model selection to explore evidence for multi-modality and a sharp cut-off in the mass distribution. We find overwhelming evidence for a bimodal distribution, in agreement with previous literature, and report for the first time positive evidence for a sharp cut-off at a maximum neutron star mass. We measure the maximum mass to be $2.0M_\odot < m_\mathrm{max} < 2.2M_\odot$ (68\%), $2.0M_\odot < m_\mathrm{max}< 2.6M_\odot$ (90\%), and evidence for a cut-off is robust against the choice of model for the mass distribution and to removing the most extreme (highest mass) neutron stars from the dataset. If this sharp cut-off is interpreted as the maximum stable neutron star mass allowed by the equation of state of dense matter, our measurement puts constraints on the equation of state. For a set of realistic equations of state that support $>2M_\odot$ neutron stars, our inference of $m_\mathrm{max}$ is able to distinguish between models at odds ratios of up to $12:1$, whilst under a flexible piecewise polytropic equation of state model our maximum mass measurement improves constraints on the pressure at $3-7\times$ the nuclear saturation density by $\sim 30-50\%$ compared to simply requiring $m_\mathrm{max}> 2M_\odot$. We obtain a lower bound on the maximum sound speed attained inside the neutron star of $c_s^\mathrm{max} > 0.63c$ (99.8\%), ruling out $c_s^\mathrm{max} < c/\sqrt{3}$ at high significance. Our constraints on the maximum neutron star mass strengthen the case for neutron star-neutron star mergers as the primary source of short gamma-ray bursts.

We identify a class of scalar-tensor theories with coupling between the scalar and the Gauss-Bonnet invariant that exhibit spontaneous scalarization for both black holes and compact stars. In particular, these theories formally admit all of the stationary solutions of general relativity, but these are not dynamically preferred if certain conditions are satisfied. Remarkably, black holes exhibit scalarization if their mass lies within one of many narrow bands. We find evidence that scalarization can occur in neutron stars as well.

Certain bulk properties of neutron stars, in particular their moment of inertia, rotational quadrupole moment and tidal Love number, when properly normalized, are related to one another in a nearly equation of state independent way. The goal of this paper is to test these relations with extreme equations of state at supranuclear densities constrained to satisfy only a handful of generic, physically sensible conditions. By requiring that the equation of state be (i) barotropic and (ii) its associated speed of sound be real, we construct a piecewise function that matches a tabulated equation of state at low densities, while matching a stiff equation of state parametrized by its sound speed in the high-density region. We show that the I-Love-Q relations hold to 1 percent with this class of equations of state, even in the extreme case where the speed of sound becomes superluminal and independently of the transition density. We also find further support for the interpretation of the I-Love-Q relations as an emergent symmetry due to the nearly constant eccentricity of isodensity contours inside the star. These results reinforce the robustness of the I-Love-Q relations against our current incomplete picture of physics at supranuclear densities, while strengthening our confidence in the applicability of these relations in neutron star astrophysics.

One of the central goals of the newborn field of gravitational wave astronomy is to test gravity in the highly nonlinear, strong field regime characterizing the spacetime of black holes. In particular, "black hole spectroscopy" (the observation and identification of black hole quasinormal mode frequencies in the gravitational wave signal) is expected to become one of the main tools for probing the structure and dynamics of Kerr black holes. In this paper we take a significant step towards that goal by constructing a "post-Kerr" quasinormal mode formalism. The formalism incorporates a parametrized but general perturbative deviation from the Kerr metric and exploits the well-established connection between the properties of the spacetime's circular null geodesics and the fundamental quasinormal mode to provide approximate, eikonal limit formulae for the modes' complex frequencies. The resulting algebraic toolkit can be used in waveform templates for ringing black holes with the purpose of measuring deviations from the Kerr metric. As a first illustrative application of our framework, we consider the Johannsen-Psaltis deformed Kerr metric and compute the resulting deviation in the quasinormal mode frequency relative to the known Kerr result.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Horndeski gravity holds a special position as the most general extension of Einstein's theory of general relativity with a single scalar degree of freedom and second-order field equations. Because of these features, Horndeski gravity is an attractive phenomenological playground to investigate the consequences of modifications of general relativity in cosmology and astrophysics. We present a review of the progress made so far in the study of compact objects (black holes and neutron stars) within Horndeski gravity. In particular, we review our recent work on slowly rotating black holes and present some new results on slowly rotating neutron stars.

We present a general formulation to analyze the structure of slowly rotating relativistic stars in a broad class of scalar-tensor theories with disformal coupling to matter. Our approach includes theories with generalized kinetic terms, generic scalar field potentials and contains theories with conformal coupling as particular limits. In order to investigate how the disformal coupling affects the structure of relativistic stars, we propose a minimal model of a massless scalar-tensor theory and investigate in detail how the disformal coupling affects the spontaneous scalarization of slowly rotating neutron stars. We show that for negative values of the disformal coupling parameter between the scalar field and matter, scalarization can be suppressed, while for large positive values of the disformal coupling parameter stellar models cannot be obtained. This allows us to put a mild upper bound on this parameter. We also show that these properties can be qualitatively understood by linearizing the scalar field equation of motion in the background of a general-relativistic incompressible star. To address the intrinsic degeneracy between uncertainties in the equation-of-state of neutron stars and gravitational theory, we also show the existence of universal equation-of-state-independent relations between the moment of inertia and compactness of neutron stars in this theory. We show that in a certain range of the theory's parameter space the universal relation largely deviates from that of general relativity, allowing, in principle, to probe the existence of spontaneous scalarization with future observations.

Horndeski's theory of gravity is the most general scalar-tensor theory with a single scalar whose equations of motion contain at most second-order derivatives. A subsector of Horndeski's theory known as "Fab Four" gravity allows for dynamical self-tuning of the quantum vacuum energy, and therefore it has received particular attention in cosmology as a possible alternative to the $\Lambda$CDM model. Here we study compact stars in Fab Four gravity, which includes as special cases general relativity ("George"), Einstein-dilaton-Gauss-Bonnet gravity ("Ringo"), theories with a nonminimal coupling with the Einstein tensor ("John") and theories involving the double-dual of the Riemann tensor ("Paul"). We generalize and extend previous results in theories of the John class and were not able to find realistic compact stars in theories involving the Paul class.

The bulk properties of spherically symmetric stars in general relativity can be obtained by integrating the Tolman-Oppenheimer-Volkoff (TOV) equations. In previous work we developed a "post-TOV" formalism - inspired by parametrized post-Newtonian theory - which allows us to classify in a parametrized, phenomenological form all possible perturbative deviations from the structure of compact stars in general relativity that may be induced by modified gravity at second post-Newtonian order. In this paper we extend the formalism to deal with the stellar exterior, and we compute several potential astrophysical observables within the post-TOV formalism: the surface redshift $z_s$, the apparent radius $R_{\rm app}$, the Eddington luminosity at infinity $L_{\rm E}^\infty$ and the orbital frequencies. We show that, at leading order, all of these quantities depend on just two post-TOV parameters $\mu_1$ and $\chi$, and we discuss the possibility to measure (or set upper bounds on) these parameters.

The lowest neutron star masses currently measured are in the range $1.0-1.1~M_\odot$, but these measurement have either large uncertainties or refer to isolated neutron stars. The recent claim of a precisely measured mass $M/M_{\odot} = 1.174 \pm 0.004$ by Martinez et al [Astrophys. J. 812, 143 (2015)] in a double neutron star system suggests that low-mass neutron stars may be an interesting target for gravitational-wave detectors. Furthermore, Sotani et al [PTEP 2014, 051E01 (2014)] recently found empirical formulas relating the mass and surface redshift of nonrotating neutron stars to the star's central density and to the parameter $\eta\equiv (K_0 L^2)^{1/3}$, where $K_0$ is the incompressibility of symmetric nuclear matter and $L$ is the slope of the symmetry energy at saturation density. Motivated by these considerations, we extend the work by Sotani et al to slowly rotating and tidally deformed neutron stars. We compute the moment of inertia, quadrupole moment, quadrupole ellipticity, tidal and rotational Love number and apsidal constant of slowly rotating neutron stars by integrating the Hartle-Thorne equations at second order in rotation, and we fit all of these quantities as functions of $\eta$ and of the central density. These fits may be used to constrain $\eta$, either via observations of binary pulsars in the electromagnetic spectrum, or via near-future observations of inspiralling compact binaries in the gravitational-wave spectrum.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

We study black hole solutions at first order in the Hartle-Thorne slow-rotation approximation in Horndeski gravity theories. We derive the equations of motion including also cases where the scalar depends linearly on time. In the Hartle-Thorne formalism, all first-order rotational corrections are described by a single frame-dragging function. We show that the frame-dragging function is exactly the same as in general relativity for all known black hole solutions in shift symmetric Horndeski theories, with the exception of theories with a linear coupling to the Gauss-Bonnet invariant. Our results extend previous no-hair theorems for a broad class of Horndeski gravity theories.

Gravitational theories with multiple scalar fields coupled to the metric and each other --- a natural extension of the well studied single-scalar-tensor theories --- are interesting phenomenological frameworks to describe deviations from general relativity in the strong-field regime. In these theories, the $N$-tuple of scalar fields takes values in a coordinate patch of an $N$-dimensional Riemannian target-space manifold whose properties are poorly constrained by weak-field observations. Here we introduce for simplicity a non-trivial model with two scalar fields and a maximally symmetric target-space manifold. Within this model we present a preliminary investigation of spontaneous scalarization for relativistic, perfect fluid stellar models in spherical symmetry. We find that the scalarization threshold is determined by the eigenvalues of a symmetric scalar-matter coupling matrix, and that the properties of strongly scalarized stellar configurations additionally depend on the target-space curvature radius. In preparation for numerical relativity simulations, we also write down the $3+1$ decomposition of the field equations for generic tensor-multi-scalar theories.

Besides their astrophysical interest, compact stars also provide an arena for understanding the properties of theories of gravity that differ from Einstein's general relativity. Numerous studies have shown that different modified theories of gravity can modify the bulk properties (such as mass and radius) of neutron stars for given assumptions on the microphysics. What is not usually stressed though is the strong degeneracy in the predictions of these theories for the stellar mass and radius. Motivated by this observation, in this paper we take an alternative route and construct a stellar structure formalism which, without adhering to any particular theory of gravity, describes in a simple parametrized form the departure from compact stars in general relativity. This "post-Tolman-Oppenheimer-Volkoff (TOV)" formalism for spherical static stars is inspired by the well-known parametrized post-Newtonian theory, extended to second post-Newtonian order by adding suitable correction terms to the fully relativistic TOV equations. We show how neutron star properties are modified within our formalism, paying special attention to the effect of each correction term. We also show that the formalism is equivalent to general relativity with an "effective" (gravity-modified) equation of state.

Some models (such as the Skyrme model, a low-energy effective field theory for QCD) suggest that the high-density matter prevailing in neutron star interiors may be significantly anisotropic. Anisotropy is known to affect the bulk properties of nonrotating neutron stars in General Relativity. In this paper we study the effects of anisotropy on slowly rotating stars in General Relativity. We also consider one of the most popular extensions of Einstein's theory, namely scalar-tensor theories allowing for spontaneous scalarization (a phase transition similar to spontaneous magnetization in ferromagnetic materials). Anisotropy affects the moment of inertia of neutron stars (a quantity that could potentially be measured in binary pulsar systems) in both theories. We find that the effects of scalarization increase (decrease) when the tangential pressure is bigger (smaller) than the radial pressure, and we present a simple criterion to determine the onset of scalarization by linearizing the scalar-field equation. Our calculations suggest that binary pulsar observations may constrain the degree of anisotropy or even, more optimistically, provide evidence for anisotropy in neutron star cores.

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.

We study torsional oscillations of neutron stars in the scalar-tensor theory of gravity using the relativistic Cowling approximation. We compute unperturbed neutron star models adopting realistic equations of state for the neutron star's core and crust. For scalar-tensor theories that allow for spontaneous scalarization, the crust thickness can be significantly smaller than in general relativity. We derive the perturbation equation describing torsional oscillations in scalar-tensor theory, and we solve the corresponding eigenvalue problem to find the oscillation frequencies. The fundamental mode (overtone) frequencies become smaller (larger) than in general relativity for scalarized stellar models. Torsional oscillation frequencies may yield information on the crust microphysics if microphysics effects are not degenerate with strong-gravity effects, such as those due to scalarization. To address this issue, we consider two different models for the equation of state of the crust and we look at the effects of electron screening. The effect of scalarization on torsional oscillation frequencies turns out to be smaller than uncertainties in the microphysics for all spontaneous scalarization models allowed by binary pulsar observations. Our study shows that the observation of quasi-periodic oscillations (QPOs) following giant flares can be used to constrain the microphysics of neutron star crusts, whether spontaneous scalarization occurs or not.

In this work we discuss a general approach for the dark energy thermodynamics considering a varying equation of state (EoS) parameter of the type $\omega(a)=\omega_0+F(a)$ and taking into account the role of a non-zero chemical potential $\mu$. We derive generalized expressions for the entropy density, chemical potential and dark energy temperature $T$ and use the positiveness of the entropy to impose thermodynamic bounds on the EoS parameter $\omega(a)$. In particular, we find that a phantom-like behavior $\omega(a)< -1$ is allowed only when the chemical potential assumes negative values ($\mu<0$).

A significant observational effort has been directed to unveil the nature of the so-called dark energy. However, given the large number of theoretical possibilities, it is possible that such a task cannot be performed on the basis only of the observational data. In this article we discuss some thermodynamic properties of this energy component by assuming that its constituents are massless quanta with a general time-dependent equation-of-state parameter $\omega(z)=\omega_0 + \omega_a f(z)$, where $\omega_0$ and $\omega_a$ are constants and $f(z)$ may assume different forms. We show that very restrictive bounds can be placed on the $w_0$ - $w_a$ space when current observational data are combined with the thermodynamic constraints derived.