Arxiv:2007.13451V2 [Gr-Qc] 24 Nov 2020 on GR and Its Modiﬁcations [31–33], but Also Delivers Model Particles As Well As Dark Matter Candidates [56– Drawbacks

Early evolutionary tracks of low-mass stellar objects in modiﬁed gravity

Aneta Wojnar1, ∗ 1Laboratory of Theoretical Physics, Institute of Physics, University of Tartu, W. Ostwaldi 1, 50411 Tartu, Estonia Using a simple model of low-mass stellar objects we have shown modified gravity impact on their early evolution, such as Hayashi tracks, radiative core development, effective temperature, masses, and luminosities. We have also suggested that the upper mass’ limit of fully convective stars on the Main Sequence might be different than commonly adopted.

I. INTRODUCTION is, in the so-called mass gap [38]), which merged with a black hole of 23M , provided even more questions for Working on modified gravity does not make one to for- theoretical physics of compact objects. get the elegance and success of Einstein’s theory of grav- However, it turns out that there is a class of stellar ity, being already confirmed by many observations [1]; objects, with the internal structure much better under- even more, General Relativity (GR) still delights when stood than that of neutron stars, which might be used to one of its mysterious predictions, such as the existence of constrain theories of gravity. It is a family of low-mass black holes, is directly affirmed by the finding of gravi- stars (LMS) [39–41] which includes such ordinary objects tational waves as a result of black holes’ binary mergers as M dwarfs (also called red dwarfs), which are cool Main [2] as well as soon after the imaging of the shadow of the Sequence stars with masses in the range [0.09 − 0.6]M , supermassive black hole of M87 [3] (see [4] for a review). brown dwarfs (stars which failed to join the Main Se- Despite the successes, there are still shortcomings quence, with masses below 0.09M ) [42, 43] or infant which, among many others, the modified gravity com- stars, that is, pre-main-sequence (PMS) stars which have munity is trying to solve by proposing extensions or/and not yet started hydrogen burning [44]. alternatives to GR. The lack of detection of dark mat- LMS are the most common stellar objects - around 70 ter and dark energy sources in order to be in agreement percent of stars in the Milky Way are red dwarfs, evolv- with the cosmological concordance model [5–10], issues ing slowly because of their small size and low masses. with unification with the high energy physics [11, 12] as Their importance becomes clear [45] when one wants to well as existence of space-time singularities [13] are just understand the properties of distant galaxies - a signif- sequent motivations for modified theories of gravity. icant part of baryonic mass is contained in LMS [46]; In order to answer the mentioned ambiguities, some globular clusters, being the oldest objects in the Universe of the extensions of GR propose to consider more gen- with reliable age determination, mainly consists of such eral than Einstein-Hilbert gravitational action (with non- objects [47]. Another, maybe the most important, argu- linear terms of curvature scalars) [18], or to include min- ment for studying LMS is our relationship with the Sun: imally or non-minimally coupled scalar fields [19, 20], the knowledge on its past and future, that is, the evo- to incorporate additional geometric ingredients [21], and lution of a star with dependent planets like the Earth, to treat the physical constants as dynamical quantities crucially determines our fate and survival possibilities. [22, 23], as a few examples. Some of those theories can be There are already many discoveries [48–55] of exoplanets already constrained by gravitational wave observations orbiting low-mass stars whose habitable properties es- [24–29]. sentially depend on star’s characteristics, from which the Astrophysics also provides additional constraints (one most important one is its mass. Moreover, these objects of the most exciting is the neutron stars’ merger [30]) are also studied in order to test properties of Standard

arXiv:2007.13451v2 [gr-qc] 24 Nov 2020 on GR and its modifications [31–33], but also delivers Model particles as well as dark matter candidates [56– drawbacks. The observations of neutron stars (NS) with 58].There are also works constraining the gravitational masses of two solar ones [34–36] clearly demonstrate how theories using data obtained from the Sun and white much is not understood yet when we try to construct a dwarfs observations [59, 60], such as, for instance, an- NS model with our current knowledge on gravitational alyzing the effect of an additional Yukawa component for interactions and nuclear physics at the extreme densi- gravity in the stellar equilibrium [61], the role of a non- ties. The lack of information on the conditions at the minimal coupling between matter and curvature [62] or NS center, where densities are above the nuclear satura- the effect of an ungravity component [63, 64]. 14 3 tion density ρ ≈ 2.8 × 10 gr/cm , forces to extrapolate The early evolution of low-mass stars has not been yet forms of equation of state (EoS). The recent discovery of examined in the context of modified gravity according to an object of 2.6 solar masses [37] (being above the heav- our knowledge (for a recent review on stellar structure in iest known NS and below the lightest black hole, that modified gravity see [65] while on astrophysical test see [66]. Main Sequence and red giant stars in modified gravity were studied in [67–69].). Due to that fact, we would ∗ [email protected] like to answer the questions if modified gravity may have 2 any significant effect on early life of such objects as well It can be shown [18] that one may rewrite the field as what kind of consequences might be carried by the equations as a dynamical equation for the conformal met- positive answer to that issue. Accordingly, we will fo- ric hµν [71, 72] and the undynamic scalar field defined as cus on the simplest example of metric-affine theories of Φ = f 0(R): gravity, that is, Palatini f(R) gravity; nevertheless, the 1 1 analogous analysis can be made for any other theory of R¯µν − hµν R¯ = κT¯µν − hµν U¯(Φ) (6a) gravity which modifies the Newtonian hydrostatic equi- 2 2 librium and related to that stellar equations. The action of Palatini f(R) gravity, which is the sim- ΦR¯ − (Φ2 U¯(Φ))0 = 0 (6b) plest generalization of GR, has the following form where we have introduced U¯(Φ) = RΦ−f(R) and appro- Z Φ2 1 √ 4 ¯ −1 S = Sg + Sm = −gf(R)d x + Sm[gµν , ψm], (1) priate energy momentum tensor Tµν = Φ Tµν . It has 2κ been already shown [73–76] that this representation of µν the Palatini gravity significantly simplifies considerations where R = R gµν is the Ricci scalar constructed with the metric g and Ricci tensor built of the independent on particular physical problems. connection Γ.ˆ Thus the common assumption on g- metricity of the connection is abandoned. Let us notice II. PALATINI STARS that we use the (− + ++) metric signature convention 8πG while κ = − 4 [70]. The variation of 1 with respect to c The stellar structure in the metric-affine theory (for the the metric gµν gives detailed review on that topic see [65]) was studied mainly 1 in the context of spherical-symmetric solutions and mass- f 0(R)R − f(R)g = κT , (2) µν 2 µν µν radius relation [77–86], the last one given by the modified Tolman-Oppenheimer-Volkoff equations. Possible issues where Tµν is the energy momentum tensor of the matter and their solutions were discussed in [87–96]. Works on 2 δSm field, obtained in the standard way Tµν = − √ . −g δgµν stability problems can be found in [76, 97–101]. Non- Later on it will be assumed to be a perfect fluid. Here, the relativistic stars which are our concern were considered primes denote derivatives with respect to the function’s in [102–107]. In what follows, we will use some results 0 df(R) derived in [76, 102–104]. argument: f (R) = dR . The variation with respect to the independent connection Γˆ provides A. Non-relativistic Palatini stars √ 0 µν ∇ˆ β( −gf (R)g ) = 0. (3) It was demonstrated [102, 103] for the Starobinski ˆ We immediately notice that ∇β is the covariant deriva- model tive calculated with respect to Γ,ˆ that is, it is the Levi- 2 Civita connection of the conformal metric f(R) = R + βR (7)

0 that non-relativistic Palatini stars can be described by hµν = f (R)gµν . (4) the following equations A very helpful equation, called the structural equation, dp Gm(˜r)ρ(˜r) is obtained from the trace of 2 taken with respect to gµν = − , (8) dr˜ Φ(˜r)˜r2 f 0(R)R − 2f(R) = κT, (5) Z r˜ m = 4πx2ρ(x)dx , (9) where T is the trace of the energy-momentum tensor. 0 For some chosen functional f(R) it is possible to solve wherer ˜2 = Φ(˜r)r2 and Φ(˜r) ≡ f 0(R) = 1 + 2κc2βρ(˜r). the structural equation 5 in order to obtain the relation After transforming back to the Jordan frame, taking the R = R(T ). Taylor expansion around β = 0 we may write down the It is a well-known fact, derived easily from 5, that in modified hydrostatic equilibrium equation as the vacuum the Palatini gravity provides Einstein vacuum solution with the cosmological constant indepen- p0 = −gρ(1 + κc2β[rρ0 − 3ρ]) , (10) dently of the f(R) form. Moreover, in the case of analytic f(R) functions it was shown [14] that the center-of-mass where g = const is the surface gravity, which can be orbits are the same as in GR while the modifications of approximated on the star’s surface (that is, on the pho- energy and momentum which appear in Euler equation tosphere, which is often taken as the surface of a star) are not sensitive to the experiments performed for the so- as lar system orbits so far. The situation may change when Gm(r) GM g ≡ ∼ , (11) atomic level experiments will be available [15–17]. r2 R2 3 where M = m(R). Let us notice that the transformation Main Sequence. Thus, the modified Lane-Emden equa- of the mass function m(˜r) to m(r) depends on the energy tion has the following form density which on the non-relativistic star’s surface will drop to zero. Due to that fact, we approximate the mass 2 2 2 1 dΦ 2 1 d √ 2κ c ρcα (Φ + 2 ξ dξ ) function to the one of the very familiar form m0(r) = Φξ θ − θn+1 = − √ θn, ξ dξ2 n + 1 4πr2ρ(r) in the Jordan frame, such that one has Φ (21) n 2 m00 = 8πrρ + 4πr2ρ0. (12) where Φ = 1 + 2αθ with α defined as α = κc βρc, while the dimensionless variables are given by We use it in 10, together with 11 written after differenta- tion with respect to r as d2m/dr2 = 2m/r2, in order to n n+1 r = rcξ, ρ = ρcθ , p = pcθ , (22) find the following form (n + 1)p r2 = c , (23) 0 g c 2 p = −gρ 1 + 8β . (13) 4πGρc c2r The heat transport with respect to radiative and con- where pc and ρc are the central pressures and densities 1 ductive processes is given by [108] and n = γ−1 is the polytropic index of the polytropic equation of state 20. Let us notice that in the standard ∂T 3 κ l = − rc , (14) version of the Lane-Emden equation (when α = 0) one ∂m 64π2ac r4T 3 deals with only one parameter n whose value is related to where l is the local luminosity, the radiation density con- a type of the star. It also indicates the stable stars’ con- −15 erg figurations (see for example [109]). The most important stant is a = 7.57 × 10 cm3K4 and the opacity κrc is given by values are n = 3/2, modelling cores of fully convective stars and small white dwarfs, and n = 3, used for high 1 1 1 = + (15) mass white dwarfs and approximated analysis of main- κrc κrad κcd sequence stars. In a modified version we have the extra parameter α, coming from the modification of GR, used with κ being the radiative opacity while κ the con- rad cd for constraining the theory against observational data (or ductive one. Writing 13 as in some cases, by theoretical analysis, as for instance the ∂p Gm Gm mentioned stability problem, discussed in [107]). = − 1 + 8β (16) ∂m 4πr4 c2r3 The equation 21 has two exact solutions [104] for n = {0, 1} and using it together with the heat transport 14 one has −1 ξ2 ξ2 − 15 ∂T 3κrcl Gm θ = − + 1, θ = , (24) = 1 + 8β . (17) n=0 n=1 2 2 2 ∂p 16πacGmT 3 c2r3 6 2ακ c ρc(10 + ξ ) Similarly as in the standard case, we define a gradient thus for the other values of the index n one needs to solve describing the temperature variation with depth the equation numerically [102, 103]. Let us notice that the equation depends on the central energy density which d ln T ∇rad := (18) is a common feature of Palatini theories of gravity. That d ln p rad is, the theory introduces new energy-density dependent contributions which distinguishes it from other proposals which in Palatini case takes a form extending GR [93]. 3κ lp Gm −1 Using the solutions of the modified Lane-Emden equa- ∇ = rc 1 + 8β . (19) rad 16πacGmT 4 c2r3 tion 21 one may obtain the star’s mass, radius, central density, and temperature via the well-known expressions

B. Polytropic Palatini stars 3 M = 4πrc ρcωn, (25) n K 3−n n−1 Since in the further part we will consider polytropic R = γn M n−3 , (26) stars whose equation of state is given by the simple G power-law relation 3M ρ = δ , (27) c n 4πR3 p = Kργ , (20) 1 Kµ n T = ρc θn, (28) it is convenient to recall now the Palatini Lane-Emden kB equation [102]. Its solutions will be needed to the discussion on the pre-main-sequence phase of the stellar evolu- where kB is the Boltzmann constant and µ the mean tion, as well as to describe fully convective stars on the molecular weight. It should be commented that the con- 4 stants 29 and 31 we may rewrite the stability condition 33 for the equation of state ρ = ρ(p, T, µ) with the homogeneous chemi- 2 3 ξ Φ 2 dθ cal composition µ (which results as dµ = 0 for both, the ωn = − |ξ=ξ , (29) 1 Φξ dξ R element and surrounding) in the following form 1 + 2 ξ Φ n−1 1 n n−3 3−n 3−n 1 dT 1 dT γn = (4π) (n + 1) ωn ξR, (30) − + > 0. (34) ξ T dr e T dr s δ = − R . (31) n − 1 Φ 2 dθ dr 3 Φ |ξ=ξR Multiplying it by the term −p one obtains the stability 1+ 1 ξ ξ dξ dp 2 Φ criterion differ from their GR forms because of the new Φ- d ln T d ln T dependent terms [104]. It is so since in order to obtain < , or ∇ < ∇ . (35) d ln p d ln p e the equation 21, the Einstein frame was used and finally s e one has to come back to the Jordan one by performing If the element changes adiabatically, we may write ∇ = the conformal transformation. e ∇ad while if the energy is transported by the radia- Let us notice that in the case of polytropies the equation (and conduction), then ∇ = ∇rad. This is the tion 19 can be written in terms of solutions of the modi- Schwarzschild criterion for the stable star’s layer: fied Lane-Emden equation [102, 103]

∇rad < ∇ad. (36) 3κ lp 4α −1 ∇ = rc 1 − , (32) rad 4 When perfect, monatomic gas is considered, then the adi- 16πacGmT 3δn abatic temperature gradient can be shown to be ∇ad=0.4 2 (see e.g. [110]). with α = κc βρc. However, if ∇rad is too high, that is, we are dealing with large flux F = l/(4πr2) or very opaque matter, or ∇ has a depression, the LHS of 36 will be bigger than III. A TOY MODEL FOR FULLY CONVECTIVE rad ∇ and a part of the flux will be carried by the convec- STARS IN PALATINI GRAVITY ad tion, so ∇= 6 ∇rad. Thus, the condition for the convective energy transport in some region of the star is ∇rad > ∇ad. A. A brief comment on dynamical instability Let us just comment that the chemical composition gradient ∇µ, which will appear on the RHS when Apart from the radiative and conductive energy trans- non-homogeneous chemical composition is considered port briefly mentioned in II A, convection is another phe- (Ledoux criterion), has a stabilizing effect. nomenon which may have an important role in the heat Since it was shown that in Palatini gravity the radia- transport in some regions of the star. It appears when tive gradient is modified, it will also have an effect on the small fluctuations of functions and variables describing Schwarzschild criterion, that is, the convection appears a spherical symmetric star, which are always present in when the star’s interior, grow: that causes mixing of the stellar −1 material as well as it may be an agent of energy transport 3κrcpl Gm 1 + 8β > ∇ . (37) through one region to another. 16πacGmT 4 c2r3 ad We will focus on an ideal gas, therefore ρ ∼ p/T . Let us consider an element ρe which remains always in the Therefore, depending on the Starobinsky parameter β, pressure balance with the surrounding ρs, so Dp := pe − the modification can have a stabilizing or destabilizing ps = 0. When the element is slightly hotter with DT > 0, effect. from the ideal gas relation we have Dρ < 0 - the element is lighter than the surrounding material and will be lifted upwards by the buoyancy forces from r to r + ∆r. The B. Convective stars change of the element density risen by dr is written as 1. Hayashi tracks dρ dρ Dρ = − ∆r. (33) dr e dr s We will consider a fully ionized monatomic gas with the temperature T and mean molecular weight µ fulling If on the considered layer we deal with Dρ > 0, the ele- the interior of a convective star up to the photosphere. ment is heavier than the surrounding and will be drawn We assume that the photosphere lies in r ∼ R, where back to its original position, so the perturbation is re- R is the star’s radius, and thus, as already mentioned, moved and we deal with a stable configuration. the stratification ∇e = d ln T/d ln p = ∇ad is adiabatic Assuming that the energy is not being exchanged dur- and equaled to 2/5. In such a case it turns out that ing that process (that is, the element rises adiabatically), the equation of state can be written as the polytropic 5 equation of state 20 with the index n = 3/2. Using the applying the absorption law 43 one gets ideal gas relation in the polytropic EoS 20 2 ! 3 M 1 − 4α µp p = 8.12 × 1014 3δ , (45) ρ = , (38) ph 4.5 1 NAkBT LTph Zµ 2 where NA and kB are the Avogardo and Boltzmann con- where we have already used the Stefan-Boltzmann law 2 4 stants, respectively, one may write L = 4πσR Tph with Teff|r=R ≡ Tph. On the other hand, from 39 taken on the photosphere 1+n NAkB with n = 3/2 we have p = KT˜ 1+n, K˜ = K−n. (39) µ 3 −2/3 2/5 5 µ 4π 2G 1 3 2 5 5 5 Let us notice that however K is a constant, it depends Teff|r=R = 1 M R pph. NAkB ξ(−θ0) 2 5 on modified gravity, since the formula (46) 1 Applying again the Stefan-Boltzmann law to the above n 4π G 1− 1 3 −1 K = M n R n (40) expression, one gets ξn+1(−θ0 (ξ ))n−1 n + 1 R n R 1 3 2 5 ! 11 L 2 Mp µ includes the solutions of the modified Lane-Emden equa- −6 ph Tph = 9.196 × 10 0 5 . (47) tion 21. So K and K˜ vary not only from star to star, −θ ξR but it also depends on the modified gravity model (here via the solutions of the modified Lane-Emden equation The pressure appearing in the above relation must be for given Starobinsky parameter β). equaled to the gravitational pressure taken on the pho- The EoS 39 holds as far as the photosphere. The photosphere, that is, the equation 45. After some algebraic manipulation and rescaling the result to the solar values tosphere is a visible surface with the temperature Teff which satisfies the Stefan-Boltzmann equation (that is, L → L/L and M → M/M , where L and M are photosphere is a surface from which the radiation is emit- solar luminosity and solar mass, respectively, one finds 1 ted into space while Teff is the temperature of a black 4 17  4α 3  1 7 1− 3δ body). For the photosphere, the optical depth τ takes 102 51 13 L M Z the value τ = 2/3. Above the photosphere one deals T = 2487.77µ 51  √  K. ph L M  5 0  with an atmosphere that we assume to be radiative for ξR −θ which the absorption law is given by a simple relation (48) κ = κ piT j. (41) Treating the star’s mass and mean molecular weight abs 0 as parameters (so each star with a mass M and uniform Since we will consider rather cool stars, for temperatures composition µ will have its own evolutionary track drawn by 48), and for given metallicity Z, we got a familiar form in the range 3000 . T . 6000K, its surface layer is domi- nated by H− opacity [108], which for hydrogen mass frac- relating the effective temperature and luminosity of the tion X ≈ 0.7 is pre-main-sequence star. The tracks given by the equation 48 are Hayashi tracks [111], that is, almost vertical 1 9 2 −1 κH− = κ0ρ 2 T cm g , (42) lines on the right-hand side of the H-R diagram. They are followed by the infant stars with masses not exceeding −31 Z where κ0 ≈ 2.5 × 10 0.02 . A metal mass fraction three solar masses, having nearly constant low effective Z is taken from the range 0.001 . Z . 0.03, with 0.02 temperature. The stars on the Hayashi tracks are fully being the solar metallicity. For the ideal gas, 42 can be convective apart from the radiative photosphere. The rewritten as relation 48 shows that their effective temperatures are

1 8.5 2 −1 almost independent of luminosity which means that Teﬀ κH− = κgp 2 T cm g , (43) is nearly independent of the way how the energy is gen-

1 erated. µ 2 −33 1 where κg = κ0 ≈ 1.371 × 10 Zµ 2 . Let us notice that the temperature coefficient is too NAkB As already mentioned, the photosphere can be de- low (it should be around 4000K) but this is the result of fined at the radius for which the optical depth with a a very simplified calculation. However, as already men- mean opacity κ (averaged over the stellar atmosphere) is tioned, the considered toy-model is good enough to ex- equaled to 2/3: amine the modified gravity effects. A few Hayashi tracks with respect to the parameter α are drawn in the Figure ∞ Z 2 1 for a star with mass M = 0.25M , mean molecular τ(r) = κ ρdr = . (44) 3 weight µ = 0.618, and solar metallicity Z = 0.02. Bigger r masses and smaller amounts of metals (Z < 0.02) give Using this relation in order to integrate the hydrostatic higher effective temperatures with a similar pattern as in equilibrium equation 13 with r = R and M = m(R), and Figure 1 for different parameters α. 6

on the H-R diagram (their luminosities remain almost constant while their eﬀective temperatures increase). When a star on the Hayashi track contracts, its luminosity decreases and thus it may happen that the α=-0.4 star’s interior becomes radiative (the turning point from Hayashi to Henyey track). However, the case when fully α=-0.1 convective star reaches the Main Sequence before devel- oping a radiative core also occurs, and this is the situation which we would like to study in our toy model. α=0 As discussed in III A, a chemically homogeneous layer is convective if α=0.1 ∇rad > ∇ad. (49)

α=0.4 A radiative region in the star’s center starts to develop when ∇rad drops below ∇ad; thus let us examine the limiting case, that is, when ∇rad = ∇ad. Therefore, for 3 the polytropic model 39 with the polytropic index n = 2 and the assumption that the deep interior opacity can be described by the Kramers’ absorption law 41 with i = 1 and j = −4.5, the temperature gradient ∇rad after entering the numerical values for the constants is Figure 1. [color online] The Hayashi tracks of a star with mass 5 0 M = 0.25M , metallicity Z = 0.02, and chemical composi- 70 δ3/2ξ θ κ0L ∇rad = 1.564 × 10 5 2 3 3.5 , (50) tion µ = 0.618 with respect to a few values of the parameter (4α − 3δ3/2) µ M R T α, given by the equation 48. where L is the total luminosity (after the homology law of contracting stars [115]). Applying the central temper- Coming back to the relation 48, let us emphasize again ature solution 28 together with the homology contraction that the difference between the GR relation and ours is argument not only given by the presence of the parameter α but also by the values of ξ and θ0, which are different in 2/3 R µδ3/2 M modified gravity - they are obtained for a given value of T = 6.679 × 10−16 , (51) c ξ5/3(−θ0)1/3 R α by solving the modified Lane-Emden equation 21. Even a small change in the constants, which will be our case, and the Stefan-Boltzmann law one gets changes the position of the Hayashi track on the H-R di- 1.25 agram. Therefore, for a non-zero Starobinsky parameter L 10.83 0 2.167 κ0 we will deal with a shift of such a path which immedi- −12 ξ (−θ ) L ∇rad = 2.71×10 1.33 8.5 5.5 , (52) ately leads to the possibility to constrain the theory by δ3/2 (3δ3/2 − 4α) µ M−1 Teff the observations of fully convective pre-main sequence stars following Hayashi tracks. Especially useful for such where we have defined M−1 = M/(0.1M ). For the ideal analysis could be T Tauri stars (for the review, see for ex- gas model that we are using the adiabatic gradient is ample [44]) which are newly formed low-mass stars, very ∇ad = 0.4, thus the minimum luminosity active and variable, in the process of contracting to the 1.064 3 0.8 Main Sequence, and which just started to be visible in δ ( δ3/2 − α) T L = 9.89 × 107L 3/2 4 eff M 4.4, the optical range. min 8.67 0 1.73 −1 ξ (−θ ) κ0 (53) which is, for the GR values in the case of a star with 2. Fully convective stars on Main Sequence Teff = 4000 K, mass M = 0.1M and opacity 62, around GR −4 33 Lmin = 0.66×10 L where L ≈ 4×10 ergs/s. Then, Let us now consider a fully convective star (apart from we may write for the arbitrary low-mass star the radiative region on the surface) at the end of its jour- 1.064 3 ney on the Hayashi track, that is, a PMS star approach- δ ( δ3/2 − α) Lmod ≈ 181.8LGR 3/2 4 . (54) ing the Main Sequence. Such an object, depending on its min min ξ8.67(−θ0)1.73 mass, may become a Brown Dwarf [43], convective star on the Main Sequence or, before reaching the Main Se- It may however happen that a low-mass star moving quence, it may start following the Henyey track as soon as along the Hayashi track can cross the Main Sequence the radiative core appears [112–114]. The last evolution- without reaching the minimum luminosity 53 needed for ary track is represented by the almost horizontal lines developping a radiative region in its center. It means 7 that such a star may reach the Main Sequence being fully The obtained expression, as we will see below, is sensi- convective (apart from the photosphere), or, if its mass tive to opacity (here to Kramers’ opacity with κ0). More- does not exceed the value ∼ 0.09M , it will fail to be a over, one usually solves much more complicated equa- Main Sequence star, that is, the star will become a Brown tions numerically, having tabulated values of opacities Dwarf ([42], for the review see [43], in modified gravity and taking into account nuclear reaction rates, more ac- [103, 116–118]). curate EoS, non-grey atmosphere models, ...[41, 119]). Let us assume that such a fully convective star has From such a theoretical analysis it is well-known that the smaller luminosity but very close to the one obtained fully convective stars on the Main Sequence have masses above (so it almost has the luminosity Lmin), and that from the range (∼ 0.09 − 0.35)M . The early analysis its mass is large enough to burn hydrogen. It means that gave the upper bound 0.26M [120] while 0.35M was we have a fully convective star on the Main Sequence, obtained by [41]. In our very simplified model one gets while another star with the luminosity (and mass) a bit masses a bit above the minimum one needed for hydro- bigger than Lmin will have a radiative center. For both gen burning. However, as argumented above, we may use stars the energy generation per unit mass in the process this theoretical crude result to see the modified gravity of hydrogen burning can be obtained by the following influence on the stars’ masses. power-law form [43] Using the GR values for the Lane-Emden solutions (α = 0), as well as αd = 4.82 and the degree of the T s ρ u−1 degeneracy electron pressure as η = 9.4 [43] one gets ˙ = ˙ , (55) pp c T ρ c c GR 0.9 −0.11 M−1 = 31.17µ Teffκ0 , (60) with s ≈ 6.31 and u ≈ 2.28, while Tc and ρc are as usual the central temperature and density, respectively, while taking the mean molecular weight µ = 0.618 and obtained from the near center solution of 21. As before, the effective temperature as Teff = 4000K: the hydrogen fraction is taken as X ≈ 0.75, while ˙ ≈ c −0.11 −9 s u−1 −1 −1 M = 4.86M κ . (61) 3.4 × 10 Tc ρc ergs g s [43]. Integrating it over 0 the stellar volume one gets the luminosity from hydrogen burning Let us consider two estimated opacities of the Kramers’ form 41 with i = 1, j = −4.5: the total bound-free and

ξR Z 2 free-free opacities [108]: 3 2 n(u+ 3 s) LHB = 4πrc ρc˙c ξ θ dξ, (56) 0 bf 25 Z(1 + X) 2 −1 κ0 ≈ 4 × 10 µ cm g , (62) where θ is a solution of the generalized Lane-Emden equa- NAkB tion. Following the result found in [103], we have for the ff 22 (X + Y )(1 + X) 2 −1 κ0 ≈ 4 × 10 µ cm g , (63) quadratic Palatini case that the luminosity from hydro- NAkB gen burning is given by the following expression for which the masses are, taking X = 0.75 and Z = 0.02, δ5.487 η10.15 respectively L = 1.53 × 107L 3/2 M 11.977 , HB ω γ16.46 −1 (η + α )16.46 3/2 3/2 d M = 0.099M ,M = 0.135M . (64) (57) bf ff where η measures the degree of the degeneracy electron As already mentioned, from our very simplified model 5µe pressure of the star and αd ≡ 2µ ≈ 4.82. Let us notice we got the mass values too small. Even so, we use these again that the above relation is modified by the values values as the reference ones to see how much they are of δ, ω, and γ. Let us just mention that the above lumi- affected by modification of our gravitational model. The nosity for GR values for a star with mass M = 0.1M results are presented in the Table I for which we have −4 is around LHB = 0.61 × 10 L [43, 103], being thus a calculated a few masses for both opacity models with re- very faint red dwarf star which belongs to the class M. It spect to the parameter α (and the corresponding quan- can be also written for a general low-mass star as tities depending on it, such as ω, γ, δ 29). δ5.487 Lmod ≈ 198LGR 3/2 . (58) HB HB 16.46 IV. CONCLUSIONS ω3/2γ3/2

Equaling the two luminosities 53 and 57 allows to ﬁnd We would like to stress again that the presented model the mass of the biggest fully convective star on the Main is too simple to describe real PMS stars and fully con- Sequence (a star onset of radiative core development) vective M dwarfs on the Main Sequence. We have not

0.9 0.11 2.173 2.173 0.132 taken into account the magnetic ﬁeld [119, 121], ther- µ Teﬀ (αd + η) γ ω monuclear depletion of the light elements [122, 123, 125], M−1 = 1.7 1.34 0.11 0.58 1.14 0 0.23 . η κ0 δ3/2 ξ (−θ ) rotation [124], protostellar initial conditions, more accu- (59) rate atmosphere and opacities models. However, our aim 8

was to demonstrate how modified gravity can affect the α Mbf/M Mff/M macroscopic values of such objects, such as mass (which -0.4 0.047 0.065 is the most crucial quantity in order to study stars’ evo- -0.1 0.083 0.114 lution), radius, effective temperature and luminosity as well as how it may change their early evolutionary tracks. 0 (GR) 0.099 0.135 Due to that fact, the considered model is close enough to 0.1 0.12 0.159 young low-mass stellar objects in order to show such an 0.4 0.18 0.24 effect. The discussion undertaken in the subsection III B 1 Table I. Numerical values of maximal masses (in solar masses) of fully convective stars on the Main Sequence for different proved that the Hayashi tracks in modified gravity are values of α = κc2βρ . significantly shifted to right or left in the H-R diagram c relative to their GR counterpart. In the GR case (when α = 0), the equation 48 describes a path of effective temperature which is taken by a PMS star with the mass criterion is also altered by the theory of gravity in such M and mean molecular weight µ - at any stage the tem- a way that the extra term appearing in the criterion has perature cannot fall below that value because the star stabilizing or destabilizing effect. Depending on the sign could enter into the forbidden region with lower tem- of the parameter, the radiative core can develop quicker peratures on the H-R diagram [111]. Since in Palatini or slower (it means, for smaller or bigger masses) than gravity this value is modified, the PMS stars with known for the masses predicted by GR [41]. In consequence, mass/luminosity moving along the Hayashi tracks could improving our toy model and then confronting it against be used to constraint the theory [126, 127]. By con- accurate empirical masses could be a powerful tool to sidering more realistic models which take into account constraint the modified gravity. Such empirical mass- missing physics (for example proper treating of the at- luminosity relations are already available (see e.g. [133– mosphere, opacity, degeneracy, rotation,...), analysis of 136]) as well as more direct and model-independent ap- Hayashi tracks of the PMS, such as T-Tauri stars [44], proaches for mass determination using detached eclipsing will be a powerful tool to constrain gravitational theories binaries (for a review see [137, 138]). which modify the stellar equations. The accurate mass determination is crucial when one It should be also noticed that for a star with mass M, shapes the evolutionary track of a star - we have shown metallicity Z, and uniform composition µ the effective that the minimum mass a star needs to have in order to temperature can be higher (lower) because of the extra develop a radiative core can be different 59 than the one term coming from the modification of the gravitational given by theoretical models using Newtonian hydrostatic theory (in other words, the star would follow a neighbor’s equilibrium. When low-mass stars considered, fully con- track of the one given by GR). Having a slightly different vective stars and stars with radiative core are modelled evolutionary track means that the star may stay longer in different ways, thus knowing that theoretical range of or shorter on the pre-main sequence phase which has an fully convective stars can have other upper limit (from effect on the total stars’ luminosity contributing to the 0.09 to 0.35M in GR based theoretical models [41]) galaxy brightness (a very nice discussion on that topic might improve existing numerical models. Moreover, can be found in [68]). that could also shed light on a discrepancy between pre- Another important point which also appeared in our dicted and dynamical masses of the M dwarfs and PMS consideration is the direct effect of modified gravity on stars with masses below 0.5M , discussed in details in, the equation of state. It was already pointed out in [94] for example, [126, 139]. (therein EiBI theory [128–131] studied) and later in [132] The general conclusion, together with the previous (general metric-affine theories) that the modified gravity works on low-mass stars in modified gravity [103, 116– introduces gravitational backreaction on the fluid parti- 118], is that not only extreme environments such as black cles, resulting as an additional pressure which should be holes, neutron stars and early/late cosmology are a back- taken into account. Commonly used EoSs, such as poly- ground for testing theories of gravity. Low-mass stellar tropic ones, are obtained by considering motion of parti- objects give additional, if not simpler (taking into ac- cles in flat spacetime so it may happen that the pressure count still unknown internal features of neutron stars or may not be local in high curvature regime which means conditions of early stage of the universe), possibility to that the covariant assumption for the pressure might not have a closer look at a bunch of gravitational theories, for hold [65]. In our simple model we have used a modified one deals with better understood density regimes allow- polytropic EoS: the polytropic constant K 40 turn out ing to examine eventual effects caused by such proposals. to depend on solutions of generalized Lane-Emden equation and thus is sensitive to modifications of gravitational equations. ACKNOWLEDGMENTS A similar situation happens with the temperature gradient which is used for examinations of dynamical sta- The author would like to thank Jakub Ostrowski for bility. As demonstrated, the Schwarzschild (or Ledoux) enlightening discussions. 9

The work is supported by the European Union through the ERDF CoE grant TK133.

[1] C. M. Will, Living Rev. Rel. 17 (2014) 4. [30] B. P. Abbott et al. [LIGO Scientific and Virgo Collab- [2] B. P. Abbott et al. [LIGO Scientific and Virgo Collab- orations], Phys. Rev. Lett. 116 (2016) 061102. orations], Phys. Rev. Lett. 119 (2017) 161101. [31] E. Berti et al., Class. Quant. Grav. 32 (2015) 243001. [3] K. Akiyama et al. [Event Horizon Telescope Collabora- [32] A. Casalino, M. Rinaldi, L. Sebastiani, S. Vagnozzi, tion], Astrophys. J. 875 (2019) L1. Phys. Dark Univ. 22 (2018) 108 [4] L. Barack et al., arXiv:1806.05195 [gr-qc]. [33] A. Casalino, M. Rinaldi, L. Sebastiani, S. Vagnozzi, [5] E. J. Copeland, M. Sami, and S. Tsujikawa, Int. J. Mod. Class. Quant. Grav. 36 (2019) 017001 Phys. D 15 (2006) 1753. [34] M. Linares, T. Shahbaz, and J. Casares, The Astrophys- [6] S. Nojiri and S. D. Odintsov, Int. J. Geom. Meth. Mod. ical Journal 859 (2018) 54. Phys. 4 (2007) 115. [35] J. Antoniadis et al., Science 340 (2012) 6131. [7] S. Nojiri, S.D. Odintsov, V.K. Oikonomou, Modified [36] F. Crawford, M. S. E. Roberts, J. W. T. Hessels, S. M. Gravity Theories on a Nutshell: Inflation, Bounce and Ransom, M. Livingstone, C. R. Tam and V. M. Kaspi, Late-time Evolution, Physics Reports 692 (2017). Astrophys. J. 652 (2006) 1499. [8] S. Nojiri, S.D. Odintsov, Unified cosmic history in mod- [37] R. Abbott et al, The Astrophysical Journal 896 L44 ified gravity: from F (R) theory to Lorentz non-invariant (2020). models, Physics Reports 505 (2011). [38] M.C. Straight, J. Sakstein, E.J. Baxter, [9] S. Capozziello and M. Francaviglia, Gen. Rel. Grav. 40 arXiv:2009.10716 . (2008) 357. [39] G. Chabrier. I. Baraffe, Annu. Rev. Astron. Astrophys., [10] S. M. Carroll, A. De Felice, V. Duvvuri, D. A. Easson, 38: 337-77 (2000). M. Trodden, and M. S. Turner, Phys. Rev. D 71 (2005) [40] G. Laughlin, P. Bodenheimer, F.C. Adams, The Astro- 063513. physical Journal, 482: 420-432, 1997. [11] L. Parker and D. J. Toms, “Quantum Field The- [41] G. Chabrier, I. Baraffe, Astronomy and Astrophysics, ory in Curved Spacetime: Quantized Fields and Grav- v.327, p.1039-1053 (1997). ity” (Cambridge University Press, Cambridge, England, [42] S.S. Kumar, The Astrophysical Journal 137 (1963): 2009). 1121. [12] N. D. Birrel and P. C. W. Davies, “Quantum Fields [43] A. Burrows and J. Liebert, Rev. Mod. Phys. 65 (1993) in Curved Space” (Cambridge University Press, Cam- 301. bridge, England, 1982). [44] C. Bertout, T Tauri stars: Wild as dust, Anny. Rev. [13] J. M. M. Senovilla and D. Garfinkle, Class. Quant. Astron. Astrophys. 27: 351-95, 1989. Grav. 32 (2015) 124008. [45] M. Catelan, AIP Conf.Proc. 930: 39-90, 2007. [14] J.D. Toniato, R.C. Rodrigues, A. Wojnar, Phys.Rev.D [46] P. Kroupa, Science, 295, 82 (2002). 101 (2020) 6, 064050. [47] A.R. Sandage, The Astronomical Journal 58 (1953): 61- [15] P.K. Schwartz, D. Giulini, Phys. Rev. A 100, 052116 75. (2019). [48] A. Cassan et al, Nature 481.7380 (2012): 167-169. [16] G.J. Olmo, Phys. Rev.D77, 084021 (2008). [49] Y. Alibert, W. Benz, Astronomy and Astrophysics, 598, [17] G.J. Olmo, Phys. Rev. Lett. 98, 061101 (2007). L5 (2017). [18] A. De Felice and S. Tsujikawa, Living Rev. Rel. 13 [50] M. Mayor, et al, Astronomy and Astrophysics 507.1 (2010) 3. (2009): 487-494. [19] C. H. Brans and R. H. Dicke, Phys. Rev. 124 (1961) [51] X. Bonfils et al, Astronomy and Astrophysics 549 925. (2013): A109. [20] P. G. Bergmann, Int. J. Theor. Phys. 1 (1968) 25. [52] J. Yang, Y. Liu, Y. Hu, D.S. Abbot, The Astrophysical [21] J. BeltránJiménez, L. Heisenberg and T. S. Koivisto, Journal, 796: L22, 2014. arXiv:1903.06830 [hep-th]. [53] B. Hansen, International Journal of Astrobiology, 14(2), [22] M. P. Dabrowski and K. Marosek, JCAP 1302 (2013) 267-278 (2015). 012. [54] R.M. Ramirez, L. Kaltenegger, The Astronomical Jour- [23] K. Leszczynska, A. Balcerzak and M. P. Dabrowski, nal 797: L25 (2014). JCAP 1502 (2015) 012. [55] A. Bohn et al, The Astrophysical Journal 898: L16 [24] J. M. Ezquiaga and M. Zumalacárregui, Phys. Rev. (2020). Lett. 119 (2017) 251304. [56] G.G. Raffelt, Astrophysics probes of particle physics, [25] T. Baker, E. Bellini, P. G. Ferreira, M. Lagos, J. Noller Physics Reports, Vol 333, 2000. and I. Sawicki, Phys. Rev. Lett. 119 (2017) 251301. [57] G.G. Raffelt, Stars as laboratories for fundamental [26] P. Creminelli and F. Vernizzi, Phys. Rev. Lett. 119 physics: The astrophysics of neutrinos, axions, and (2017) 251302. other weakly interacting particles, University of Chicago [27] D. Langlois, R. Saito, D. Yamauchi and K. Noui, Phys. press, 1996. Rev. D 97 (2018) 061501. [58] J.P.P. Vieira, C.J.A.P. Martins, M.J.P.F.G. Monteiro, [28] J. Sakstein and B. Jain, Phys. Rev. Lett. 119 (2017) Phys.Rev.D 86 (2012) 043003. 251303. [59] I.D. Saltas, I. Lopes, Phys. Rev. Lett. 123 (2019) 9, [29] L. Lombriser and A. Taylor, JCAP 1603 (2016) 031. 091103. 10

[60] R.K. Jain, C. Kouvaris, N. G. Nielsen, Phys. Rev. Lett. 040. 116, 151103. [96] G.J. Olmo, D. Rubiera-Garcia, arXiv:2007.04065. [61] O. Bertolami and J. Páramos,Phys. Rev. D 77, 084018. [97] T-H. Sham, L-M. Lin, P. T. Leung, Phys. Rev. D 86 [62] O. Bertolami and J. Páramos,Phys. Rev. D 71, 023521. (2012) 064015. [63] O. Bertolami, J. PPáramos,and P. Santos, Phys. Rev. [98] H. Sotani, Phys. Rev. D 89 (2014) 124037. D 80, 022001. [99] M. Z. Bhatti, Z. Yousaf, Zarnoor, Gen. Rel. Grav. 51 [64] O. Bertolami, H. Mariji, Phys. Rev. D 93, 104046. (2019) 144. [65] G.J. Olmo, D. Rubiera-Gracia, A. Wojnar, Stellar struc- [100] P. Pani, T. Delsate, and V. Cardoso, Phys. Rev. D 85 ture models in modified theories of gravity: lessons and (2012) 084020. challenges, Physics Reports 876 (2020) 1-75. [101] W-X. Feng, C-Q. Geng, L-W Luo, Chin. Phys. C 43 [66] T. Baker et al The Novel Probes Project–Tests of Grav- (2019) 083107. ity on Astrophysical Scales, arXiv:1908.03430 (2019). [102] A. Wojnar, Eur. Phys. J. C79 (2019) no.1, 51. [67] P. Chang, L. Hui, The Astrophysical Journal 732.1 [103] G.J. Olmo, D. Rubiera-Garcia, A. Wojnar, Phys.Rev.D (2011): 25. 100 (2019) 4, 044020. [68] A-Ch. Davis et al, Physical Review D 85.12 (2012): [104] A. Sergyeyev, A. Wojnar, Eur.Phys.J.C 80 (2020) 4, 123006. 313. [69] S. Chowdhury, T. Sarkar, arXiv:2008.12264. [105] S. Banerjee, S. Shankar, T.P. Singh, JCAP 2017 (2017) [70] S. Weinberg, Gravitation and Cosmology: principles 004. and Applications of the General Theory of Relativity, [106] C. Wibisono, A. Sulaksono, Int. J. Mod. Phys. D, 27 John Wiley and Sons (1972). (2018) 1850051. [71] A. Stachowski, M. Szydlowski, A. Borowiec, Eur. Phys. [107] A. Wojnar, Acta Phys.Polon.Supp. 13 (2020) 249. J. C77, 406 (2017). [108] C. J. Hansen, S. D. Kawaler, V. Trimble, Stellar Interi- [72] M. Szydlowski, A. Stachowski, A. Borowiec, Eur. Phys. ors. Physical Principles, Structure, and Evolution, 2nd J. C77, 603 (2017). edition, Springer 2004. [73] V.I. Afonso, G.J. Olmo, D. Rubiera-Garcia, Phys.Rev.D [109] N. K. Glendenning, Compact Stars, Springer 1996 97 (2018) 2, 021503. [110] R. Kippenhahn, A. Weigert, A, Weiss, Stellar structure [74] V.I. Afonso, G.J. Olmo, E. Orazi, D. Rubiera-Garcia, and evolution, Springer, 2012. Eur.Phys.J.C 78 (2018) 10, 866. [111] C. Hayashi, Publication of the Astronomical Society of [75] V.I. Afonso, G.J. Olmo, E. Orazi, D. Rubiera-Garcia, Japan, Vol.13, 450-452 (1961). Phys.Rev.D 99 (2019) 4, 044040. [112] L. Henyey, R. Lelevier, and R.D. Levee, Publications of [76] A. Wojnar, Eur. Phys. J. C78 (2018) no.5, 421. the Astronomical Society of the Pacific 67.396 (1955): [77] K. Kainulainen, V. Reijonen, D. Sunhede, Phys. Rev. 154-160. D. 76 (2007) 043503. [113] L. Henyey, J.E. Forbes, N.L. Gould, The Astrophysical [78] F. A. T. Pannia, F. Garcia, S. E. P. Bergliaffa, M. Orel- Journal 139 (1964): 306. lana, M., G. E. Romero, Gen. Rel. Grav. 49 (2017) 25. [114] L. Henyey, M.S. Vardya, P. Bodenheimer, The Astro- [79] D. E. Barraco, V. H. Hamity, Phys. Rev. D 62 (2000) physical Journal 142 (1965): 841. 044027. [115] M. Szydlowski, Marek, A.J. Maciejewski, Journal of [80] V. Reijonen, [arXiv:0912.0825 [gr-qc]]. Physics A: Mathematical and General 37.10 (2004): [81] G. Panotopoulos, Gen. Rel. Grav. 49 (2017) 69. 3501. [82] T. Harko, FSN Lobo, MK Mak, SV Sushkov, Phys. Rev. [116] J. Sakstein, Phys. Rev. Lett. 115 (2015) 201101. D 88 (2013) 044032. [117] J. Sakstein, Phys. Rev. D 92 (2015) 124045. [83] T. Harko, FSN Lobo, MK Mak, SV Sushkov, Phys. Rev. [118] M. Crisostomi, M. Lewandowski and F. Vernizzi, Phys. D 88 (2013) 044032. Rev. D 100 no.2, 024025 (2019). [84] A. I. Qauli, M. Iqbal, A. Sulaksono, H.S. Ramadhan, [119] M. Cantiello, J. Braithwaite, The Astrophysical Jour- Phys. Rev. D 93 (2016) 104056. nal, 883: 106, 2019. [85] A. I. Qauli, A. Sulaksono, H. S. Ramadhan, I. Husin, [120] C. Hayashi, T. Nakano, Progress of Theoretical Physics, arXiv:1710.03988. 30(4), 460-474, (1963). [86] B. Danila, T. Harko, F. S. N. Lobo, M. K. Mak, Phys. [121] J. Morin et al., Monthly Notices of the Royal Astronom- Rev. D 95 (2017) 044031. ical Society 384.1 (2008): 77-86. [87] E. Barausse, T. P. Sotiriou, J. C. Miller, Class. Quant. [122] L. Bildsten, E.F. Brown, C.D. Matzner, G. Ushomirsky, Grav. 25 (2008) 062001. The Astrophysical Journal 482.1 (1997): 442. [88] E. Barausse, T. P. Sotiriou, J C. Miller, Class. Quant. [123] G. Ushomirsky, C.D. Matzner, E.F. Brown, L. Bildsten, Grav. 25 (2008) 105008. V.G. Hilliard, P.C Schroeder, The Astrophysical Jour- [89] E. Barausse, T. P. Sotiriou, J C. Miller, EAS Publica- nal 497 : 253-266, 1998. tions Series 30 (2008) 189. [124] N.J. Wright, et al, Monthly Notices of the Royal Astro- [90] P. Pani and T. P. Sotiriou, Phys. Rev. Lett. 109 (2012) nomical Society 479.2 (2018): 2351-2360. 251102 [125] A. Wojnar, arXiv:2009.10983 [91] Y-H. Sham, P. T. Leung, L-M. Lin, Phys. Rev. D 87 [126] L.A. Hillenbrand, R.J. White, The Astrophysical Jour- (2013) 061503. nal 604.2 (2004): 741. [92] J. Barrientos O. and G. F. Rubilar, Phys. Rev. D 93 [127] C.A. Tout, M. Livio, I.A. Bonnell, Mon. Not. R. Astron. (2016) 024021. Soc. 310, 360-376 (1999). [93] G.J. Olmo, Phys.Rev.D 78 (2008) 104026. [128] M. Banados and P. G. Ferreira, Phys.Rev.Lett. 105, [94] H-Ch. Kim, Phys. Rev. D 89 (2014) 064001. 011101 (2010). [95] A. Mana, L. Fatibene, M. Ferraris, JCAP 2015 (2015) [129] T. Delsate and J. Steinhoff, Phys. Rev. Lett. 109, 021101 11

(2012). tronomical Society, Vol. 481, Issue 1, 2018. [130] J. Beltran Jimenez, L. Heisenberg, G. J. Olmo and [136] A.W. Mann et al., The Astrophysical Journal, 871: 63, D. Rubiera-Garcia, Phys. Rept. 727 (2018) 1 2019. [131] . N. Vollick, Phys. Rev. D 69, 064030 (2004). [137] G. Chabrier et al., The mass-radius relationship from [132] A. Delhom-Latorre, G.J. Olmo, M. Ronco, Physics Let- solar-type stars to terrestrial planets: a review, AIP ters B 780 (2018): 294-299. Conference Proceedings. Vol. 1094. No. 1. American In- [133] X. Delfosse, T. Forveille, D. SA˜©gransan, J.L. Beuzit, stitute of Physics, 2009. S. Udry, C. Perrier, M. Mayor, Astronomy and Astro- [138] A. Serenelli et al., arXiv:2006.10868. physics, 364, 217-224 (2000). [139] L. Siess, E. Dufour, M. Forestini, Astron. Astrophys. [134] A.W. Mann et al., The Astrophysical Journal, 804: 64, 358: 593-599, 2000. 2015. [135] S.G. Parsons et al., Monthly Notices of the Royal As-