Search

What if the dark matter Sector is truly dark, self-interacting and unreachable by terrestrial experiments? How could we find hints of such dark sector if it is experimentally unreachable by any terrestrial experiment? We study a low-scale mirror Standard Model which can act as a model for dark matter, which interacts only gravitationally with the Standard Model particles. The mirror Standard Model sector particles are stable particles can comprise a measurable part of the dark matter of the Universe. These mirror Standard Model particles acquire mass through a low-temperature dark first order phase transition. We examine in detail this dark phase transition and we indicate how stochastic gravitational waves can be generated through this transition. For the model we use, the produced energy spectrum of the gravitational waves can be detected by the Square Kilometer Array. Moreover, we propose a possible way to detect effects of the particle nature of dark matter, using observational data coming from the guitar nebula, which can work if dark matter is collisional, so interacting dark matter. Without specifying a model for interacting dark matter, thus choosing an agnostic approach for interacting dark matter, we assume that the guitar nebula bow shock is generated by the interaction of the high speed neutron star that passes through the interstellar medium, which is assumed to be comprised by interacting dark matter and hydrogen. Our main proposal is that the opening angle of the bow shock can be directly related to the speed of sound of the dark matter particles, and a large angle of the bow shock could be a strong indicator that the interstellar medium is comprised by collisional dark matter and hydrogen gas. This is motivated by the Bosma effect., Comment: Revised version, PRD Accepted

In this work we investigate the phenomenological implications of several non-trivial axion-Higgs couplings, which cover most of the possible non-perturbative scenarios. Specifically we consider the combination of having higher order non-renormalizable Higgs-axion couplings originating from a weakly coupled effective theory combined with non-perturbative couplings of the form $\sim \epsilon \Lambda_c^2|H|^2\cos (\frac{\phi}{f_a})$. Since we consider the misalignment axion, the non-perturbative couplings can be expanded in the form of a perturbation expansion in powers of $\phi/f_a$, thus after the electroweak symmetry breaking, the effective potential of the axion is drastically affected by these terms. We investigate the phenomenological implications of these terms for various values of the mass scale $\Lambda_c$, and some scenarios are theoretically disfavored, while other scenarios with non-perturbative Higgs-axion couplings of the form $\sim \epsilon \Lambda_c^2|H|^2\cos (\frac{\phi}{f_a})$ with $\Lambda_c\sim 10^{-10}\times m_a$ and $m_a\sim 10^{-10}\,$eV, lead to a characteristic pattern in the stochastic gravitational wave background via the deformation of the background equation of state parameter occurring at frequencies probed by the Einstein Telescope. We also consider loop effects from the Higgs sector caused by the term $\sim \epsilon \Lambda_c^2|H|^2\cos (\frac{\phi}{f_a})$ if we close the Higgs in one loop., Comment: PRD accepted

In this work we study the inflationary framework of Einstein-Gauss-Bonnet theories which produce a primordial gravitational wave speed that respects the constraint imposed by the GW170817 event, namely $\left| c_T^2 - 1 \right| < 6 \times 10^{-15}$ in natural units. For these theories in general, the scalar Gauss-Bonnet coupling function $\xi (\phi)$ is arbitrary and unrelated to the scalar potential. We develop the inflationary formalism of this theory, and present analytic forms for the slow-roll indices, the observational indices and the $e$-foldings number, assuming solely a slow-roll era and also that the constraints imposed by the GW170817 event hold true. We present in detail an interesting class of models that produces a viable inflationary era. Finally, we investigate the behavior of the Einstein-Gauss-Bonnet coupling function during the reheating era, in which case the evolution of the Hubble rate satisfies a constant-roll-like relation $\dot{H}=\delta H^2$. As we show, the behavior of the scalar Gauss-Bonnet coupling during the reheating is different from that of the inflationary era. This is due to the fact that the evolution governing the reheating is different compared to the inflationary era., Comment: PLB Accepted

In this work we highlight an important perspective for the complete understanding of the stochastic gravitational background structure. The stochastic gravitational wave background is perhaps the most important current and future tool towards pinpointing the early Universe phenomenology related with the inflationary era and the subsequent reheating era. Many mysteries are inherent to the stochastic spectrum so in this work we highlight the fact that the complete understanding of early Universe physics and of astrophysical processes requires data from many distinct frequency band ranges. The combination of these data will provide a deeper and better understanding of the physics that forms the stochastic gravitational wave background, in both cases that it is of cosmological or astrophysical origin. We also discuss how the reheating temperature may be determined by combining multi-band frequency data from gravitational wave experiments and we also discuss how the shape of the gravitational wave energy spectrum can help us better understand the physical processes that formed it., Comment: Physics of the Dark Universe accepted, contribution to the special issue ''Proceedings of Cosmology and Astrophysics in the Gravitational Wave Astronomy Era''

We consider the implications of a scalar field interacting with the dark matter fluid on the energy spectrum of the primordial gravitational waves. We choose an interaction type which before a critical matter density $\rho_m^c$ during the reheating era or early radiation domination era, the scalar field loses energy transferring it to the dark matter fluid, while after the critical matter density $\rho_m^c$ the dark matter fluid loses energy transferring it to the scalar field. The scalar field is assumed to have an exponential potential and at the critical matter density with $\rho_m\sim \rho_m^c$, at which point, the interaction between the scalar and the dark matter fluid is switched off, we demand that the effective equation of state of scalar field is described by a matter dominated era. This is crucial since it affects the behavior of the trajectories and the fixed points of the two dimensional dynamical system composed by the dark matter fluid and the scalar field. Specifically the phase space contains two stable dark matter dominated final attractors and two unstable stiff era dominated fixed points. Thus, there exists the remarkable possibility that the Universe might feel the passing of the scalar field through the unstable kinetic dominated fixed points, during the reheating era, with the total equation of state parameter of the Universe being deformed to be larger than $w=1/3$. This deformation of the total equation of state parameter during the reheating era can potentially have significant effects on the energy spectrum of the primordial gravitational waves. Also the model we use also contains an $F(R)$ gravity which controls in a dominant way inflation and the late-time acceleration, in a phenomenologically viable way., Comment: PRD Accepted

In this work we consider an axionic scalar-tensor theory of gravity and its effects on static neutron stars. The axionic theory is considered in the regime in which the axion oscillates around its potential minimum, which cosmologically occurs post-inflationary, when the Hubble rate is of the same order as the axion mass. We construct the TOV equations for this axionic theory and for a spherically symmetric static spacetime and we solve these numerically using a quite robust double shooting LSODA based python integration method. Regarding the equations of state, we used nine mainstream and quite popular ones, namely, the WFF1, the SLy, the APR, the MS1, the AP3, the AP4, the ENG, the MPA1 and the MS1b, using the piecewise polytropic description for each. From the extracted data we calculate the Jordan frame masses and radii, and we confront the resulting phenomenology with five well-known neutron star constraints. As we demonstrate, the AP3, the ENG and the MPA1 equations of state yield phenomenologically viable results which are compatible with the constraints, with the MPA1 equation of state enjoying an elevated role among the three. The reason is that the MPA1 fits well the phenomenological constraints. A mentionable feature is the fact that all the viable phenomenologically equations of state produce maximum masses which are in the mass-gap region with $M_{max}>2.5M_{\odot}$, but lower that the causal 3 solar masses limit. We also compare the neutron star phenomenology produced by the axionic scalar-tensor theory with the phenomenology produced by inflationary attractors scalar-tensor theories., Comment: CQG Accepted

The significance of the electroweak phase transition is undeniable, and although initially it was believed that it was second-order, it is now believed that it is a first-order transition. However, it is not a strong first-order phase transition in the context of the Standard Model and the remedy to this issue is to use the Higgs portal and directly couple the Higgs to a hidden scalar sector. This can result in a strong electroweak phase transition, while the couplings to a hidden scalar are constrained by several phenomenological constraints, such as the sphaleron rate criterion and the branching ratio of the Higgs to invisible channels. In this work, we consider the standard singlet extensions of the Standard Model, including dimension-six non-renormalizable operators that couple a real singlet scalar field with the Higgs doublet. As a result, we examine the effects of those Higgs-singlet couplings on the electroweak phase transition. The effective theory, where the non-renormalizable couplings originate from, is considered to be active beyond 15$\,$TeV. As we show, the Universe experiences a two-step electroweak phase transition, a primary phase transition in the singlet sector at a high temperature, and then a subsequent first-order phase transition from the singlet vacuum to the electroweak vacuum. The singlet's phase transition can either be second-order or first-order, depending on the singlet mass and its couplings to the Higgs. In particular, we show that the dimension-six operator assists in generating a strong electroweak phase transition in regions of the parameter space that were excluded in the previous singlet extensions of the Standard Model. This is further apparent for low singlet masses \(m_S < m_H/2\) which are rarely taken into account in the literature due to the invisible branching ratio of the Higgs boson., Comment: PRD Accepted, Abstract reduced due to arXiv restrictions, minor typos fixed

In this work we aim to revive the interest for non-minimal derivative coupling theories of gravity, in light of the GW170817 event. These theories include a string motivated non-minimal kinetic term for the scalar field of the form $\sim \xi(\phi) G^{\mu\nu}\partial_\mu\phi\partial_\nu\phi$ and predict that the primordial tensor perturbations have a speed that it is distinct from the speed of light. Due to the fact that the Universe is classical during and in the post-inflationary epoch, there is no fundamental reason for the graviton to change its mass, the GW170817 event severely constrained these theories. We analyze and formalize the inflationary phenomenology of these theories, using both the latest Planck data and the GW170817 event data to constrain these theories. Due to the fact that there is no constraint in the choice of the scalar potential and the non-minimal coupling function, we provide several classes of viable models, using convenient forms of the minimal coupling function in terms of the scalar potential, aiming for analyticity, and we discuss the advantages and disadvantages of each viable model., Comment: NPB Accepted

In this work we shall prove that the tensor spectral index of the primordial tensor perturbations for GW170817-compatible Einstein-Gauss-Bonnet theories, takes the approximate simplified form $n_{\mathcal{T}}\simeq 2\left(-1+\frac{1}{\lambda(\phi)} \right)\epsilon_1$ at leading order, with $\lambda (\phi)$ being a function of the scalar field which depends on the scalar field potential and the second derivative of the scalar-Gauss-Bonnet coupling $\xi''(\phi)$. With our analysis we aim to provide a definitive criterion for selecting Einstein-Gauss-Bonnet models that can provide a blue-tilted inflationary phenomenology, by simply looking at the scalar potential and the scalar-Gauss-Bonnet coupling. We shall prove this using two distinct approaches and as we show the tilt of the tensor spectral index is determined by the values of the potential $V(\phi)$ and of scalar-Gauss-Bonnet coupling at first horizon crossing. Specifically the blue-tilted tensor spectral index can occur when $\xi''(\phi_*)V(\phi_*)>0$ at first horizon crossing., Comment: EPL Accepted

In this work we study the GW170817-compatible Einstein-Gauss-Bonnet theories during the reheating and the end of inflationary era. Given the scalar field potential $V(\phi)$ which can have some intrinsic importance for the theory, determining the scalar coupling function $\xi(\phi)$ can be cumbersome due to lack of analyticity. The GW170817 observation constrains the scalar coupling function and the scalar field potential to have some interdependence, thus during the slow-roll era one can calculate the scalar coupling function. However, when the slow-roll era ends, it is expected that the scalar coupling function should have a different form and the same applies for the reheating era, assuming that the scalar potential of the theory does not change. In this work we exactly aim to highlight this feature of Einstein-Gauss-Bonnet theories, as the Universe evolves through distinct sequential evolution eras, and we focus on how to determine the scalar coupling function during the various evolutionary eras, from inflation to the reheating era. Regarding both the end of the inflationary era and the reheating era, it is found that the Hubble rate obeys a constant-roll-like condition of the form $\dot{H}=\delta H^2$, thus the determination of the scalar Gauss-Bonnet function $\xi(\phi)$ is reduced to solving a differential equation. A mentionable feature of the era exactly at the end of inflation is that the Klein-Gordon equation is decoupled from the field equations, because the Gauss-Bonnet invariant is zero. We provide several examples of interest to support our arguments., Comment: AoP accepted

In this work we consider the realization of a variant emergent Universe scenario in the context of $f(R)$ gravity. We use well-known reconstruction techniques existing in the literature, and we find the approximate form of the vacuum $f(R)$ that reproduces the specific variant emergent Universe scale factor in the large curvature approximation. As we show, in the variant emergent Universe scenario, the Hubble horizon shrinks primordially and the Universe undergoes in an accelerated expansion era. In a perturbation theory approach, the scalar and tensor curvature perturbations can be expressed in terms of the phenomenological indices $\epsilon_i , i=1,3,4$, usually used for inflationary phenomenology, and we extract the spectral indices for the scalar and tensor perturbations, along with the tensor-to-scalar ratio expressed in terms of the perturbation indices. Using a powerful genetic algorithm we investigate in depth the parameter space of the model, quantified by several free parameters, in order to study the viability of the model when this is compared with the Planck 2018 data. We also investigate the implications of the vacuum $f(R)$ gravity we found on the Big Bang nucleosynthesis., Comment: IJMPD Accepted. arXiv admin note: text overlap with arXiv:2309.04850

In this work, we propose a model of Einstein--Gauss-Bonnet gravity coupled with two scalar fields. The two scalar fields are considered to be ``frozen'' or they become non-dynamical by employing appropriate constraints in terms of Lagrange multiplier fields. We show that, even in the case that the arbitrary spherically symmetric spacetime is dynamical, we can construct a model where the wormhole spacetime is a stable solution in this framework. We especially concentrate on the model reproducing the dynamical wormhole, where the wormhole appears in a finite-time interval. We investigate the propagation of the gravitational wave in the wormhole spacetime background and we show that the propagation speed is different from that of light $\to$ light in general, and there is a difference in the speeds between the incoming propagating wave and the outgoing propagating gravitational wave., Comment: arXiv admin note: text overlap with arXiv:2311.06932

In this work, we examine the propagation of gravitational waves in cosmological and astrophysical spacetimes in the context of Einstein--Gauss-Bonnet gravity, in view of the GW170817 event. The perspective we approach the problem is to obtain a theory which can produce a gravitational wave speed that is equal to that of light in the vacuum, or at least the speed can be compatible with the constraints imposed by the GW170817 event. As we show, in the context of Einstein--Gauss-Bonnet gravity, the propagation speed of gravity waves in cosmological spacetimes can be compatible with the GW170817 event, and we reconstruct some viable models. However, the propagation of gravity waves in spherically symmetric spacetimes violates the GW170817 constraints, thus it is impossible for the gravitational wave that propagates in a spherically symmetric spacetime to have a propagating speed which is equal to that of light in the vacuum. The same conclusion applies to the Einstein--Gauss-Bonnet theory with two scalars. We discuss the possible implications of our results on spherically symmetric spacetimes., Comment: PRD Accepted

The far infrared physics is a fascinating topic for theoretical physics, since the foundation of quantum field theory and neutrinos seem to be strongly related with the far infrared physics of our Universe. In this work we shall explore the possibility of a late-time thermal phase transition caused by the axion-neutrino interactions. The axion is assumed to be the misalignment axion which is coupled primordially to a chiral symmetric neutrino. The chiral symmetry is supposed to be broken either spontaneously or explicitly, and two distinct phenomenological models of axion-neutrinos are constructed. The axion behaves as cold dark matter during all its evolution eras, however if we assume that the axion and the neutrino fields interact coherently in a classical way as fields, or as ensembles, then we consider thermal effects in the axion sector, due to the values of operators $\phi$ for the axion and $\bar{\nu}\nu$ due to the neutrinos. The thermal equilibrium between the two has no effect to the axion effective potential for a wide temperature range. As we show, contrary to the existing literature, the axion never becomes destabilized due to the finite temperature effects, however if axion-Higgs higher order non-renormalizable operators are present in the Lagrangian, the axion potential is destabilized in the temperature range $T\sim 0.1\,$MeV down to $T\sim 0.01\,$eV and a first order phase transition takes place. The initial axion vacuum decays to the energetically more favorable axion vacuum, and the latter decays to the Higgs vacuum which is more preferable. This late-time phase transition might take place in the redshift range $z\sim 385-37$ and thus it may cause density fluctuations in the post-recombination era., Comment: Accepted in Physics of the Dark Universe, abstract reduced due to arXiv limitations

Singularities in any physical theory are either remarkable indicators of the unknown underlying fundamental theory, or indicate a change in the description of the physical reality. In General Relativity there are three fundamental kinds of singularities that might occur, firstly the black hole spacelike crushing singularities, e.g. in the Schwarzschild case and two cosmological spacelike singularities appearing in finite-time, namely, the Big Bang singularity and the Big Rip singularity. In the case of black hole and Big Bang singularity, the singularity indicates that the physics is no longer described by the classical gravity theory but some quantum version of gravity is probably needed. The Big Rip is a future singularity which appears in the context of General Relativity due to a phantom scalar field needed to describe the dark energy era. Apart from the Big Rip singularity, a variety of finite-time future singularities, such as, sudden singularity, Big Freeze singularity, generalized sudden singularity, $w$-singularity and so on, are allowed in various class of cosmological models irrespective of their origin. The occurrence of these finite-time singularities has been intensively investigated in the context of a variety of dark energy, modified gravity, and other alternative cosmological theories. These singularities suggest that the current cosmological scenario is probably an approximate version of a fundamental theory yet to be discovered. In this review we provide a concrete overview of the cosmological theories constructed in the context of Einstein's General Relativity and modified gravity theories that may lead to finite-time cosmological singularities. We also discuss various approaches suggested in the literature that could potentially prevent or mitigate finite-time singularities within the cosmological scenarios., Comment: Invited Review Article from Physics Reports; Published version; 138 pages including 7 tables, 1 figure and the references

We examine the effects of a stiff pre-recombination era on the present day's energy spectrum of the primordial gravitational waves. If the background total equation of state parameter at the pre-recombination era is described by a kination era one, this directly affects the modes with characteristic wavenumbers which reenter the Hubble horizon during this stiff era. The stiff era causes a broken-power-law effect on the energy spectrum of the gravitational waves. We use two approaches, one model agnostic and a specific model that can realize this scenario. In all cases, the inflationary era can be realized either by some theory leading to a standard red-tilted tensor spectral index or by some theory which has a mild tensor spectral index $n_{\mathcal{T}}=0.17-0.37$ like an Einstein-Gauss-Bonnet theory. For the model agnostic scenario case, the NANOGrav signal can be explained by the stiff pre-recombination era combined with an inflationary era with a mild blue-tilted tensor spectral index $n_{\mathcal{T}}=0.37$ and a low-reheating temperature $T_R\sim 0.1$GeV. In the same case, the red-tilted inflationary theory signal can be detectable by the future LISA, BBO and DECIGO experiments. The model dependent approach is based on a Higgs-axion model which can yield multiple deformations of the background total equation of state parameter, causing multiple broken-power-law behaviors occurring in various eras before and after the recombination era. In this case, the NANOGrav signal is explained by this model in conjunction with an inflationary era with a really mild blue-tilted tensor spectral index $n_{\mathcal{T}}=0.17$ and a low-reheating temperature $T_R\sim 20\,$GeV. In this case, the signal can be detectable by the future Litebird experiment, which is a very characteristic pattern in the tail of the primordial gravitational wave energy spectrum., Comment: Refs added, clarifications added

We review recent trends on inflationary dynamics in the context of viable modified gravity theories. After providing a general overview of the inflationary paradigm emphasizing on what problems of hot Big Bang theory inflation solves, and a somewhat introductory presentation of single field inflationary theories with minimal and non-minimal couplings, we review how inflation can be realized in terms of several string motivated models of inflation, which involve Gauss-Bonnet couplings of the scalar field, higher order derivatives of the scalar field, and some subclasses of viable Horndeski theories. We also present and analyze inflation in the context of Chern-Simons theories of gravity, including various subcases and generalizations of string corrected modified gravities which also contain Chern-Simons correction terms, with the scalar field being identified with the invisible axion, which is the most viable to date dark matter candidate. We also provide a detailed account of vacuum $f(R)$ gravity inflation, and also inflation in $f(R,\phi)$ and kinetic-corrected $f(R,\phi)$ theories of gravity. In the end of the review we discuss the technique for calculating the overall effect of modified gravity on the waveform of the standard general relativistic gravitational wave form., Comment: Prepared review for the Feature Articles 2023 of the journal Symmetry

In this work we present several characteristic examples of theories of gravity and particle physics scenarios that may yield an observable energy spectrum of stochastic primordial gravitational waves, compatible with the 2023 NANOGrav observations. The resulting theories yield a flat or a peak-like energy spectrum, and we further seek the conditions which if hold true, the energy spectrum can be compatible with the recent NANOGrav stochastic gravitational wave detection. As we show, in most cases a blue tilted spectrum combined with a relatively low reheating temperature is needed, the scale of which is determined by whether the radiation domination era is ordinary or it is an abnormal radiation domination era. One intriguing Higgs-axion model, which predicts short slow-roll eras for the axion field at the post-electroweak breaking epoch, which eventually change the total equation of state parameter at the reheating era, can explain the NANOGrav signal, if a blue tilted tensor spectral index inflationary era precedes the reheating era, and a reheating temperature of the order $\mathcal{O}(400)\,$GeV. This specific model produces an energy spectrum of primordial gravitational waves with a characteristic peak that is detectable from both the NANOGrav and future LISA experiment, but not from the future Einstein telescope., Comment: PRD Accepted

We study a power-law $F(R)$ gravity with an early dark energy term, that can describe both the early-time and the late-time acceleration of the Universe. We confront this scenario with recent observational data including the Pantheon Type Ia supernovae, measurements of the Hubble parameter $H(z)$ (Cosmic Chronometers), data from Baryon Acoustic Oscillations and standard rulers data from the Cosmic Microwave Background (CMB) radiation. The model demonstrates some achievements in confronting with these observations and can be compared with the $\Lambda$-Cold-Dark-Matter model. In particular, in both models we obtain very close estimates for the Hubble constant $H_0$, but it is not true for $\Omega_m^0$. The early dark energy term supports viability of the considered $F(R)$ gravity model., Comment: PLB Accepted

In this work we study static neutron stars in the context of several inflationary models which are popular in cosmology. These inflationary models are non-minimally coupled scalar theories which yield a viable inflationary phenomenology in both Jordan and Einstein frames. By considering the constraints from inflationary theories, which basically determine the values of the potential strength, usually considered as a free parameter in astrophysical neutron star works, we construct and solve the Tolman-Oppenheimer-Volkoff equations using a solid python-3 LSODA integrator. For our study we consider several popular inflationary models, such as the universal attractors, the $R^p$ attractors (three distinct model values), the induced inflation, the quadratic inflation, the Higgs inflation and the $a$-attractors (two distinct model values) and for the following popular equations of state the WFF1, the SLy, the APR, the MS1, the AP3, the AP4, the ENG, the MPA1 and the MS1b. We construct the $M-R$ diagram and we confront the resulting theory with theoretical and observational constraints. As we demonstrate, remarkably, all the neutron stars produced by all the inflationary models we considered are compatible with all the constraints for the MPA1 equation of state. It is notable that for this particular equation of state, the maximum masses of the neutron stars are in the mass-gap region with $M>2.5M_{\odot}$, but lower than the 3 solar masses causal limit. We also make the observation that as the NICER constraints are pushed towards larger radii, as for example in the case of the black widow pulsar PSR J0952-0607, it seems that equations of state that produce neutron stars with maximum masses in the mass gap region, with $M>2.5M_{\odot}$, but lower than the 3 solar masses causal limit, are favored and are compatible with the modified NICER constraints., Comment: PRD Accepted, abstract reduced due to arXiv restrictions

In this paper, we worked in the framework of an inflationary $f(R,T)$ theory, in the presence of a canonical scalar field. More specifically, the $f(R,T)=\gamma R+2\kappa\alpha T$ gravity. The values of the dimensionless parameters $\alpha$ and $\gamma$ are taken to be $\alpha \geq 0$ and $0 < \gamma \leq 1$. The motivation for that study was the striking similarities between the slow-roll parameters of the inflationary model used in this work and the ones obtained by the rescaled Einstein-Hilbert gravity inflation $f(R)=\alpha R$. We examined a variety of potentials to determine if they agree with the current Planck Constraints. In addition, we checked whether these models satisfy the Swampland Criteria and we specified the exact region of the parameter space that produces viable results for each model. As we mention in Section IV the inflationary $f(R,T)$ theory used in this work can not produce a positive $n_T$ which implies that the stochastic gravitational wave background will not be detectable., Comment: IJMPD Accepted

In this article we present an alternative formalism for the inflationary phenomenology of rescaled Einstein-Gauss-Bonnet models which are in agreement with the GW170817 event. By constraining the propagation velocity of primordial tensor perturbations, an approximate form for the time derivative of the scalar field coupled to the Gauss-Bonnet density is extracted. In turn, the overall degrees of freedom decrease and similar to the case of the canonical scalar field, only one scalar function needs to be designated, while the other is extracted from the continuity equation of the scalar field. We showcase explicitly that the slow-roll indices can be written in a closed form as functions of three dimensionless parameters, namely $x=\frac{1}{2\alpha}\bigg(\frac{\kappa\xi'}{\xi''}\bigg)^2$, $\beta=8H^2\xi''$ and $\gamma=\frac{\xi'\xi'''}{\xi''^2}$ and in turn, we prove that the Einstein-Gauss-Bonnet model can in fact produce a blue-tilted tensor spectral index if the condition $\beta\geq1$ is satisfied, which is possible only for Einstein-Gauss-Bonnet models with $\xi''(\phi_k)>0$. Afterwards, a brief comment on the running of the spectral indices is made where it is shown that $a_{\mathcal{S}}(k_*)$ and $a_{\mathcal{T}}(k_*)$ in the constrained case are approximately of the order $\mathcal{O}(10^{-3})$, if not smaller. Last but not least, we examine the conditions under which the Swampland criteria are satisfied. We connect the tracking condition related to scalar field theories with the present models, and we highlight the important feature of the models we propose that the tracking condition can be satisfied only if the Swampland criteria are simultaneously satisfied, however the cases with $\xi\sim1/V$ and $\xi\sim V$ are excluded, as they cannot describe the inflationary era properly., Comment: PRD Accepted, abstract reduced due to arXiv limitations

This study focuses on the static neutron star perspective for two types of cosmological inflationary attractor theories, namely the induced inflationary attractors and the quadratic inflationary attractors. The two cosmological models can be discriminated cosmologically, since one of the two does not provide a viable inflationary phenomenology, thus in this paper we investigate the predictions of these theories for static neutron stars, mainly focusing on the mass and radii of neutron stars. We aim to demonstrate that although the models have different inflationary phenomenology, the neutron star phenomenology predictions of the two models are quite similar. We solve numerically the Tolman-Oppenheimer-Volkoff equations in the Einstein frame using a powerful double shooting numerical technique, and after deriving the mass-radius graphs for three types of polytropic equations of state, we derive the Jordan frame mass and radii. With regard the equations of state we use polytropic equation of state with the small density part being either the WFF1, the APR or the intermediate stiffness equation of state SLy. The results of our models will be confronted with quite stringent recently developed constraints on the radius of neutron stars with specific mass. As we show, the only equation of state which provides results compatible with the constraints is the SLy, for both the quadratic and induced inflation attractors. Thus nowadays, scalar-tensor descriptions of neutron stars are quite scrutinized due to the growing number of constraining observations, which eventually may also constrain theories of inflation., Comment: CQG Accepted

We investigate the effects of short axion kination eras on the energy spectrum of the primordial gravitational waves corresponding to modes that re-enter the Hubble horizon at the post-electroweak symmetry breaking epoch, well within the radiation domination era. Our main assumption is the existence of an extremely weakly coupled hidden sector between the Higgs and the axion, materialized by higher order non-renormalizable dimension six and dimension eight operators, active at a scale M of the order 20-100TeV. This new physics scale M which is way higher than the electroweak scale, is motivated by the lack of new particle observations in the large hadron collider to date, beyond the electroweak scale. Once the electroweak symmetry breaking occurs at T\sim 100 GeV, the axion potential acquires a new minimum due to the new terms generated by the electroweak breaking, and the axion oscillations at the origin are destabilized. In effect after some considerable amount of time, the axion rolls swiftly to its new minimum, experiencing a short kination epoch, where its energy density redshifts as $\rho_a\sim a^{-6}$. After it reaches the new minimum, since the latter is energetically less favorable that the Higgs minimum, it decays to the Higgs minimum and the Universe is described again by the Higgs minimum. The axion returns to the origin and commences again oscillations initiated by quantum fluctuations, redshifts as dark matter, and the same procedure is repeated perpetually. These short axion kination eras may disturb the background total equation of state parameter during the radiation domination era. As we show, the energy spectrum of the gravitational waves mainly depends on how many times the short axion kination epochs occur., Comment: PRD Accepted, Abstract reduced due to arXiv limitations

Under the main assumption that the axion scalar field mainly composes the dark matter in the Universe, in this paper we shall extend the formalism of kinetic axion $R^2$ gravity to include Gauss-Bonnet terms non-minimally coupled to the axion field. As we demonstrate, this non-trivial Gauss-Bonnet term has dramatic effects on the inflationary phenomenology and on the kinetic axion scenario. Specifically, in the context of our formalism, the kinetic axion ceases to be kinetically dominated at the end of the inflationary era, since the condition $\dot{\phi}\simeq 0$ naturally emerges in the theory. Thus, unlike the case of kinetic axion $R^2$ gravity, the Gauss-Bonnet corrected kinetic axion $R^2$ gravity leads to an inflationary era which is not further extended and the reheating era commences right after the inflationary era, driven by the $R^2$ fluctuations., Comment: Published in the Physics of the Dark Universe

In this work we study the neutron star phenomenology of $R^p$ attractor theories in the Einstein frame. The Einstein frame $R^p$ attractor theories have the attractor property that they originate from a large class of Jordan frame scalar theories with arbitrary non-minimal coupling. These theories in the Einstein frame provide a viable class of inflationary models, and in this work we investigate their implications on static neutron stars. We numerically solve the Tolman-Oppenheimer-Volkoff equations in the Einstein frame, for three distinct equations of state, and we provide the mass-radius diagrams for several cases of interest of the $R^p$ attractor theories. We confront the results with several timely constraints on the radii of specific mass neutron stars, and as we show, only a few cases corresponding to specific equations of state pass the stringent tests on neutron stars phenomenology., Comment: MNRAS Accepted

Motivated by the fact that the pre-inflationary era may evolve in an exotic way, in this work we formalize anisotropic evolution in the context of modified gravity, focusing on pre-inflationary and near the vicinity of the inflationary epochs. We specialize on specific metrics like Bianchi and Taub and we formalize the inflationary theory in vacuum $F(R)$ gravity, in $F(R)$ gravity with an extra scalar field and in Gauss-Bonnet gravity. We discuss the qualitative effects of the anisotropies on the evolution of the Universe and also we consider several specific solutions, like the de Sitter solution, in both the isotropic and anisotropic contexts. Furthermore, several exotic modified gravity cosmological solutions, like the ones which contain finite time singularities, are also discussed in brief., Comment: NPB Accepted

In this work we shall study a new class of attractor models which we shall call generalized $R^p$-attractor models. This class of models is based on a generalization of the Einstein frame potential of $R^p$ $f(R)$ gravity models in the Jordan frame. We present the attractor properties of the corresponding non-minimally coupled Jordan frame theory, and we calculate the observational indices of inflation in the Einstein frame. As we show, there is a large class of non-minimally coupled scalar theories, with an arbitrary non-minimal coupling which satisfies certain conditions, that yield the same Einstein frame potential, this is why these models are characterized attractors. As we demonstrate, the generalized $R^p$-attractor models are viable and well fitted within the Planck constraints. This includes the subclass of the generalized $R^p$-attractor models, namely the Einstein frame potential of $R^p$ inflation in the Jordan frame, a feature also known in the literature. We also highlight an important issue related to the $R^p$ inflation in the Jordan frame, which is known to be non-viable. By conformal invariance, the $R^p$ inflation model should also be viable in the Jordan frame, which is not. We pinpoint the source of the problem using two different approaches in the $f(R)$ gravity Jordan frame, and as we demonstrate, the problem arises in the literature due to some standard simplifications made for the sake of analyticity. We demonstrate the correct way to analyze $R^p$ inflation in the Jordan frame, using solely the slow-roll conditions., Comment: IJMPD Accepted

In this work we examine the effects of a pre-inflationary de Sitter bounce on the energy spectrum of the primordial gravitational waves. Specifically we assume that the Universe is described by several evolution patches, starting with a de Sitter pre-inflationary bounce which is followed by an quasi-de Sitter slow-roll inflationary era, followed by a constant equation of state parameter abnormal reheating era, which is followed by the radiation and matter domination eras and the late-time acceleration eras. The bounce and the inflationary era can be realized by vacuum $f(R)$ gravity and the abnormal reheating and the late-time acceleration eras by the synergy of $f(R)$ gravity and the prefect matter fluids present. Using well-known reconstruction techniques we find which $f(R)$ gravity can realize each evolution patch, except from the matter and radiation domination eras which are realized by the corresponding matter fluids. Accordingly, we calculate the damping factor of the primordial de Sitter bounce, and as we show, the signal can be detected by only one gravitational wave future experiment, in contrast to the case in which the bounce is absent. We discuss in detail the consequences of our results and the future perspectives., Comment: Nuclear Physics B in press. arXiv admin note: text overlap with arXiv:2209.09781

In this work we consider a combined theoretical framework comprised by $F(R)$ gravity and a kinetic scalar field. The kinetic energy of the scalar field dominates over its potential for all cosmic times, and the kinetic scalar potential is chosen to be small and non-trivial. In this case, we show that the primordial gravitational wave energy spectrum of vacuum $F(R)$ gravity is significantly enhanced and can be detectable in future interferometers. The kinetic scalar thus affects significantly the inflationary era, since it extends its duration, but also has an overall amplifying effect on the energy spectrum of pure $F(R)$ gravity primordial gravitational waves. The form of the signal is characteristic for all these theories, since it is basically flat and should be detectable from all future gravitational wave experiments for a wide range of frequencies, unless some unknown damping factor occurs due to some unknown physical process., Comment: Astroparticle Physics in press

In this work we consider the effect of an $R^2$ term on the kinetic misalignment axion theory. By using the slow-roll assumptions during inflation and the field equations, we construct an autonomous dynamical system for the kinetic axion, including the effects of the $R^2$ term and we solve numerically the dynamical system. As we demonstrate, the pure kinetic axion attractor is transposed to the right in the field phase space, and it is no longer $(\phi,\dot{\phi})=(\langle \phi \rangle,0)$, but it is $(\phi,\dot{\phi})=(\langle \phi '\rangle,0)$, with $\langle \phi '\rangle\neq 0$ some non-zero value of the scalar field with $\langle \phi '\rangle> \langle \phi \rangle$. This feature indicates that the kinetic axion mechanism is enhanced, and the axion oscillations are further delayed, compared with the pure kinetic axion case. The phenomenological implications on the duration of the inflationary era, on the commencing of the reheating era and the reheating temperature, are also discussed., Comment: EPL in press

In this short review, we discuss how Earth's climatological and geological history and also how the shadows of galactic black holes might reveal our Universe's past evolution. Specifically we point out that a pressure singularity that occurred in our Universe's past might have left its imprint on Earth's geological and climatological history and on the shadows of cosmological black holes. Our approach is based on the fact that the $H_0$ tension problem may be resolved if some sort of abrupt physics change occurred in our Universe $70-150\,$Myrs ago, an abrupt change that deeply affected the Cepheid parameters. We review how such an abrupt physics change might have been caused in our Universe by a smooth passage of it through a pressure finite-time singularity. Such finite-time singularities might occur in modified gravity and specifically in $F(R)$ gravity, so we show how modified gravity might drive this type of evolution, without resorting to peculiar cosmic fluids or scalar fields. The presence of such a pressure singularity can distort the elliptic trajectories of bound objects in the Universe, causing possible geological and climatological changes on Earth, if its elliptic trajectory around the Sun might have changed. Also, such a pressure singularity affects directly the circular photon orbits around supermassive galactic black holes existing at cosmological redshift distances, thus the shadows of some cosmological black holes at redshifts $z\leq 0.01$, might look different in shape, compared with the SgrA* and M87* supermassive black holes. This feature however can be checked experimentally in the very far future., Comment: Invited article accepted in the journal Universe, special issue ''Modified Gravity Approaches to the Tensions of LCDM''

The Laser Interferometer Space Antenna (LISA) mission, scheduled for launch in the mid-2030s, is a gravitational wave observatory in space designed to detect sources emitting in the millihertz band. LISA is an ESA flagship mission, currently entering the Phase B development phase. It is expected to help us improve our understanding about our Universe by measuring gravitational wave sources of different types, with some of the sources being at very high redshifts $z\sim 20$. On the 23rd of February 2022 we organized the 1$^\mathrm{st}$ {\it LISA in Greece Workshop}. This workshop aimed to inform the Greek scientific and tech industry community about the possibilities of participating in LISA science and LISA mission, with the support of the Hellenic Space Center (HSC). In this white paper, we summarize the outcome of the workshop, the most important aspect of it being the inclusion of $15$ Greek researchers to the LISA Consortium, raising our total number to $22$. At the same time, we present a road-map with the future steps and actions of the Greek Gravitational Wave community with respect to the future LISA mission., Comment: 11 pages, 2 figures

The black hole observations obtained so far indicate one thing: similar "donuts" exist in the sky. But what if some of the observed black hole shadows that will obtained in the future are different from the others? In this work the aim is to show that a difference in the shadow of some observed black holes in the future, might explain the $H_0$-tension problem. In this letter we investigate the possible effects of a pressure cosmological singularity on the circular photon orbits and the shadow of galactic supermassive black holes at cosmological redshifts. Since the pressure singularity is a global event in the Universe, the effects of the pressure singularity will be imposed on supermassive black holes at a specific redshift. As we show, the pressure singularity affects the circular photon orbits around cosmological black holes described by the McVittie metric, and specifically, for some time before the time instance that the singularity occurs, the photon orbits do not exist. We discuss the possible effects of the absence of circular photon orbits on the shadow of these black holes. Our idea indicates that if a pressure singularity occurred in the near past, then this could have a direct imprint on the shadow of supermassive galactic black holes at the redshift corresponding to the time instance that the singularity occurred in the past. Thus, if a sample of shadows is observed in the future for redshifts $z\leq 0.01$, and for a specific redshift differences are found in the shadows, this could be an indication that a pressure singularity occurred, and this global event might resolve the $H_0$-tension as discussed in previous work. However, the observation of several shadows at redshifts $z\leq 0.01$ is rather a far future task., Comment: EPL Accepted

In this work we investigate the quantitative effects of the misalignment kinetic axion on $R^2$ inflation. Due to the fact that the kinetic axion possesses a large kinetic energy which dominates its potential energy, during inflation its energy density redshifts as stiff matter fluid and evolves in a constant-roll way, making the second slow-roll index to be non-trivial. At the equations of motion level, the $R^2$ term dominates the evolution, thus the next possible effect of the axion could be found at the cosmological perturbations level, via the second slow-roll index which is non-trivial. As we show, the latter elegantly cancels from the observational indices, however, the kinetic axion extends the duration of the inflationary era to an extent that it may cause a 15$\%$ decrease in the tensor-to-scalar ratio of the vacuum $R^2$ model. This occurs because as the $R^2$ model approaches its unstable quasi-de Sitter attractor in the phase space of $F(R)$ gravity due to the $\langle R^2 \rangle $ fluctuations, the kinetic axion dominates over the $R^2$ inflation and in effect the background equation of state is described by a stiff era, or equivalently a kination era, different from the ordinary radiation domination era. This in turn affects the duration of the inflationary era, increasing the $e$-foldings number up to $5$ $e$-foldings in some cases, depending on the reheating temperature, which in turn has a significant quantitative effect on the observational indices of inflation and especially on the tensor-to-scalar ratio., Comment: PRD Accepted

In this work we shall provide a model-independent general calculation of the running of the spectral index for vacuum $F(R)$ gravities. We shall exploit the functional form of the spectral index and of the tensor-to-scalar ratio in order to present a general $n_s-r$ relation for vacuum $F(R)$ gravity theories. As we show, viable $F(R)$ gravity theories can be classified to two classes of models regarding their prediction for the running spectral index. The $R^2$-attractor models predict a running of the spectral index in the range $-10^{-3}

In this work we investigate the effects of a geometrically generated early dark energy era on the energy spectrum of the primordial gravitational waves. The early dark energy era, which we choose it to have a constant equation of state parameter $w$, is synergistically generated by an appropriate $f(R)$ gravity in the presence of matter and radiation perfect fluids. As we demonstrate, the predicted signal for the energy spectrum of the $f(R)$ primordial gravitational waves is amplified and can be detectable, for various reheating temperatures, especially for large reheating temperatures. The signal amplitude depends on the duration of the early dark energy era and on the value of the dark energy equation of state parameter, with the most latter affecting more crucially the amplification. Specifically the amplification occurs when the equation of state parameter approaches the de Sitter value $w=-1$. Regarding the duration of the early dark energy era, we find that the largest amplification occurs when the early dark energy era commences at a temperature $T=0.85\,$eV until $T=7.8\,$eV. Moreover we study a similar scenario in which amplification occurs, where the early dark energy era commences at $T=0.29\,$eV and lasts until the temperature is increased by $\Delta T\sim 1.7\,$eV., Comment: Accepted in Symmetry, invited article

In this work we shall consider the effects of a geometrically generated post-inflationary era on the energy spectrum of the primordial gravitational waves. Specifically, we shall consider a post-inflationary constant equation of state era, generated by the synergistic effect of $f(R)$ gravity and of radiation and matter perfect fluids. Two cases of interest shall be studied, one with equation of state parameter $w=-1/3$, in which case the Universe neither accelerates nor decelerates, and one with $w=0$ so an early matter domination era. For the evaluation of the inflationary observational indices which is relevant for the calculation of the gravitational waves energy spectrum, we also took into account the effects of the constant equation of state parameter era, on the $e$-foldings number. In both the $w=-1/3$ and $w=0$ cases, the energy spectrum of the primordial gravitational waves is amplified, but for the $w=0$ case, the effect is stronger., Comment: Accepted in Annalen de Physik

In this work, we shall provide an $F(R)$ gravity theoretical framework for solving the $H_0$-tension. Specifically, by exploiting the $F(R)$ gravity correspondence with a scalar-tensor theory, we shall provide a condition in which when it is satisfied, the $H_0$-tension is alleviated. The condition that remedies the $H_0$-tension restricts the corresponding $F(R)$ gravity, and we present in brief the theoretical features of the constrained $F(R)$ gravity theory in both the Jordan and Einstein frames. The condition that may remedy the $H_0$-tension is based on the existence of a metastable de Sitter point that occurs for redshifts near the recombination. This metastable de Sitter vacuum restricts the functional form of the $F(R)$ gravity in the Jordan frame. We also show that by appropriately choosing the $F(R)$ gravity, along with the theoretical solution offered for the $H_0$-tension problem, one may also provide a unified description of the inflationary era with the late-time accelerating era, in terms of two extra de Sitter vacua. We propose a new approach to $F(R)$ gravity by introducing a new class of integral $F(R)$ gravity functions, which may be wider than the usual class expressed in terms of elementary $F(R)$ gravity functions. Finally, the Einstein frame inflationary dynamics formalism is briefly discussed., Comment: NPB Accepted

In this work we shall study the late-time dynamics of several $F(R)$ gravity models. By appropriately expressing the field equations in terms of the redshift and of a statefinder function, we shall solve numerically the field equations using appropriate physical motivated initial conditions. We consider models which, by construction, are described by a nearly $R^2$-model at early epochs and we fine tune the parameters to achieve viability and compatibility with the latest Planck constraints at late times. Such models provide a unified description of inflation and dark energy era and notably a common feature of all the models is the presence of dark energy oscillations. Furthermore, we show that, in contrast to general relativistic fluids and scalar field descriptions, a large spectrum of different dark energy physics is generated by simple $F(R)$ gravity models, varying from phantom, to nearly de Sitter and to quintessential dark energy eras., Comment: IJMPD Accepted

In this paper we shall consider an axionic Chern-Simons corrected $f(R)$ gravity theoretical framework, and we shall study the chirality of the generated primordial gravitational waves. Particularly, we shall consider two main axion models, the canonical misalignment axion model and the kinetic axion model, both of which provide an interesting particle phenomenology, in the presence of $R^2$ terms in the inflationary Lagrangian. The axion does not affect significantly the background evolution during the inflationary era, which is solely controlled by $R^2$ gravity. However, the due to the presence of the Chern-Simons term, the tensor perturbations are directly affected, and our aim is to reveal the extent of effects of the Chern-Simons term on the gravitational waves modes, for both the axion models. We aim to produce analytical descriptions of the primordial tensor modes behavior, and thus we solve analytically the evolution equations of the tensor modes, for a nearly de Sitter primordial evolution controlled by the $R^2$ gravity. We focus the analytical study on superhorizon and subhorizon modes. For the misalignment model, we were able to produce analytic solutions for both the subhorizon and superhorizon modes, in which case we found the behavior of the circular polarization function. Our results indicate that the produced tensor spectrum is strongly chiral. For the kinetic axion model though, analytic solutions are obtained only for the superhorizon modes. In order to have a grasp of the behavior of the chirality of the tensor modes, we studied the chirality of the superhorizon modes, however a more complete study is needed, which is impossible to do analytically though., Comment: PRD Accepted

In this work we shall calculate in detail the effect of an GW170817-compatible Einstein-Gauss-Bonnet theory on the energy spectrum of the primordial gravitational waves. The spectrum is affected by two characteristics, the overall amplification/damping factor caused by the GW170817-compatible Einstein-Gauss-Bonnet theory and by the tensor spectral index and the tensor-to-scalar ratio. We shall present the formalism for studying the inflationary dynamics and post-inflationary dynamics of GW170817-compatible Einstein-Gauss-Bonnet theories for all redshifts starting from the radiation era up to the dark energy era. We exemplify our formalism by using two characteristic models, which produce viable inflationary and dark energy eras. As we demonstrate, remarkably the overall damping/amplification factor is of the order of unity, thus the GW170817-compatible Einstein-Gauss-Bonnet models affect the primordial gravitational waves energy spectrum only via their tensor spectral index and the tensor-to-scalar ratio. Both models have a blue tilted tensor spectrum, and therefore the predicted energy spectrum of the primordial gravity waves can be detectable by most of the future gravitational waves experiments, for various reheating temperatures., Comment: Accepted in Astroparticle Physics

String theory enjoys an elevated role among quantum gravity theories, since it seems to be the most consistent UV completion of general relativity and the Standard Model. However, it is hard to verify the existence of this underlying theory on terrestrial accelerators. One way to probe string theory is to study its imprints on the low-energy effective inflationary Lagrangian, which are quantified in terms of high energy correction terms. It is highly likely, thus, to find higher order curvature terms combined with string moduli, that is scalar fields, since both these types of interactions and matter fields appear in string theory. In this work we aim to stress the probability that the inflationary dynamics are controlled by the synergy of scalar fields and higher order curvature terms. Specifically, we shall consider a well motivated quantum corrected canonical scalar field theory, with the quantum corrections being of the $\mathcal{R}^2$ type. The reason for choosing minimally coupled scalar theory is basically because if scalar fields are evaluated in their vacuum configuration, they will either be minimally coupled or conformally coupled. Here we choose the former case, and the whole study shall be performed in the string frame (Jordan frame), in contrast to similar studies in the literature where the Einstein frame two scalar theory is considered. We derive the field equations of the quantum-corrected theory at leading order and we present the form the slow-roll indices obtain for the quantum corrected theory. We exemplify our theoretical framework by using the quadratic inflation model, and as we show, the $\mathcal{R}^2$ quantum corrected quadratic inflation model produces a viable inflationary phenomenology, in contrast with the simple quadratic inflation model., Comment: NPB Accepted

In this work we shall exhaustively study the effects of modified gravity on the energy spectrum of the primordial gravitational waves background. S. Weinberg has also produced significant works related to the primordial gravitational waves with the most important one being the effects of neutrinos on primordial gravitational waves. With this sort review, our main aim is to gather all the necessary information for studying the effects of modified gravity on primordial gravitational waves in a concrete and quantitative way and in a single paper. After reviewing all the necessary techniques for extracting the general relativistic energy spectrum, and how to obtain in a WKB way the modified gravity damping or amplifying factor, we concentrate on specific forms of modified gravity of interest. The most important parameter involved for the calculation of the effects of modified gravity on the energy spectrum is the parameter $a_M$ which we calculate for the cases of $f(R,\phi)$ gravity, Chern-Simons-corrected $f(R,\phi)$ gravity, Einstein-Gauss-Bonnet-corrected $f(R,\phi)$ gravity, and higher derivative extended Einstein-Gauss-Bonnet-corrected $f(R,\phi)$ gravity. The exact forms of $a_M$ is presented explicitly for the first time in the literature. With regard to Einstein-Gauss-Bonnet-corrected $f(R,\phi)$ gravity, and higher derivative extended Einstein-Gauss-Bonnet-corrected $f(R,\phi)$ gravity theories, we focus on the case that the gravitational wave propagating speed is equal to that of light's in vacuum. We provide expressions for $a_M$ expressed in terms of the cosmic time and in terms of the redshift, which can be used directly for the numerical calculation of the effect of modified gravity on the primordial gravitational wave energy spectrum., Comment: Invited review, Symmetry accepted

In order to describe inflation in general relativity, scalar fields must inevitably be used, with all the setbacks of that description. On the other hand, $f(R)$ gravity and other modified gravity theories seem to provide a unified description of early and late-time dynamics without resorting to scalar or phantom theories. The question is, can modified gravity affect directly the mysterious radiation domination era? Addressing this question is the focus in this work, and we shall consider the case for which in the early stages of the radiation domination era, namely during the reheating era, the background equation of state parameter is different from $w=1/3$. As we show, in the context of $f(R)$ gravity, an abnormal reheating era can affect the primordial gravitational wave energy spectrum today. Since future interferometers will exactly probe this era, which consists of subhorizon modes that reentered the horizon during the early stages of the radiation domination era, the focus in this work is how a short abnormal reheating era that deviates from the standard perfect fluid pattern with $w\neq 1/3$, and generated by higher order curvature terms, can affect the primordial gravitational wave energy spectrum. Using a WKB approach, we calculate the effect of an $f(R)$ gravity generated abnormal reheating era, and as we show the primordial gravitational wave spectrum is significantly amplified, a result which is in contrast to the general relativistic case, where the effect is minor., Comment: Fortschritte der Physik Accepted

In this work we shall use a bottom-up approach for obtaining viable inflationary Einstein-Gauss-Bonnet models which are also compatible with the GW170817 event. Specifically, we shall use a recently developed theoretical framework in which we shall specify only the tensor-to-scalar ratio, in terms of the $e$-foldings number. Starting from the tensor-to-scalar ratio, we shall reconstruct from it the Einstein-Gauss-Bonnet theory which can yield such a tensor-to-scalar ratio, finding the scalar potential and the Gauss-Bonnet coupling scalar function as functions of the $e$-foldings number. Accordingly, the calculation of the spectral index of the primordial scalar perturbations, and of the tensor spectral index easily is greatly simplified and these observational indices can easily be found. After presenting the general formalism for the bottom-up reconstruction, we exemplify our findings by presenting several Einstein-Gauss-Bonnet models of interest which yield a viable inflationary phenomenology. These models have also an interesting common characteristic, which is a blue tilted tensor spectral index. We also investigate the predicted energy spectrum of the primordial gravitational waves for these Einstein-Gauss-Bonnet models, and as we show, all the models yield a detectable primordial wave energy power spectrum., Comment: CQG accepted, some typos corrected, one author name modified, doi added

In the frame of Gauss-Bonnet gravity and in the limit of $D\to 4$, based on the fact that spherically symmetric solution derived using any of regularization schemes will be the same form as the original theory \cite{Banerjee:2020yhu,Hennigar:2020lsl,Casalino:2020kbt,Aoki:2020lig}, we derive a new interior spherically symmetric solution assuming specific forms of the metric potentials that have two constants . Using the junction condition we determine these two constants. By using the data of the star EXO $1785-248 $, whose mass is $ { M=1.3\pm0.2 \,{\textrm M}_\odot}$ and radius $ { l=8.849\pm0.4\,{\textrm km}}$, we calculate the numerical values of these constants, in terms of the dimensionful coupling parameter of the Gauss-Bonnet term, and eventually, we get real values for these constants. In this regard, we show that the components of the energy-momentum tensor has a finite value at the center of the star as well as a smaller value to the surface of the star. Moreover, we show that the equations of the state behave in a non-linear way due to the impact of the Gauss-Bonnet term. Using the Tolman-Oppenheimer-Volkoff equation, the adiabatic index, and stability in the static state we show that the model under consideration is always stable. Finally, the solution of this study is matched with observational data of other pulsars showing satisfactory results, Comment: 24 pages, Latex, 5 figures, Accepted in Journal of Symmetry

In the recent literature it has been shown that the $H_0$ tension may be eliminated if an abrupt physics transition changed the Cepheid parameters in the near past of the Universe, nearly $70-150\,$Myrs ago. In this letter we stress the possibility that this abrupt transition was caused by the smooth passage of our Universe through a pressure finite-time cosmological singularity. Being a non-crushing type singularity the pressure singularity can leave its imprints in the Universe, since it occurs globally and literally everywhere. We discuss how this scenario could easily be realized by $F(R)$ gravity, with the strong energy conditions being satisfied without the need for a scalar field or specific matter fluids. We also stress the fact that the pressure singularity can affect the effective gravitational constant of $F(R)$ gravity. Moreover, we stress the fact that pressure singularities can disrupt the trajectories of bound objects in the Universe, which is also pointed out in the literature, even in the context of general relativity. We also show numerically in a general relativistic framework that elliptic trajectories are distorted and changed to different elliptic trajectories when the Universe passes through the pressure singularity. Such a disruption of the trajectories could have tidal effects on the surface of the earth, for example on sea waters and oceans, regarding the distortion of moon's elliptic trajectory. Accordingly, the distortion of earth's trajectory around the sun could have affected climatologically earth $70-150\,$Myrs ago., Comment: EPL Accepted

Neutron stars are perfect candidates to investigate the effects of a modified gravity theory, since the curvature effects are significant and more importantly, potentially testable. In most cases studied in the literature in the context of massive scalar-tensor theories, inflationary models were examined. The most important of scalar-tensor models is the Higgs model, which, depending on the values of the scalar field, can be approximated by different scalar potentials, one of which is the inflationary. Since it is not certain how large the values of the scalar field will be at the near vicinity and inside a neutron star, in this work we will answer the question, which potential form of the Higgs model is more appropriate in order for it to describe consistently a static neutron star. As we will show numerically, the non-inflationary Higgs potential, which is valid for certain values of the scalar field in the Jordan frame, leads to extremely large maximum neutron star masses, however the model is not self-consistent, because the scalar field approximation used for the derivation of the potential, is violated both at the center and at the surface of the star. These results shows the uniqueness of the inflationary Higgs potential, since it is the only approximation for the Higgs model, that provides self-consistent results., Comment: Invited article to appear in the Symmetry MDPI special issue "The Nuclear Physics of Neutron Stars"

The vacuum $R^2$ model is known to generate a quasi-de Sitter evolution for inflation, using solely the slow-roll assumptions. Using standard reconstruction techniques, we demonstrate that the $f(R)$ gravity which actually realizes the quasi-de Sitter evolution is not simply the $R^2$ model but a deformed $R^2$ model which contains extra terms in addition to the $R^2$ model. We analyze in detail the inflationary dynamics of the deformed $R^2$ model and we demonstrate that the predictions are quite close to the ones of the pure $R^2$ model, regardless the values of the free parameters. Basically the deformed $R^2$ model is also a single parameter inflationary model, exactly like the ordinary $R^2$ model. In contrast to the early-time era, where the deformed $R^2$ model is quite similar to the $R^2$ model, at late times, the phenomenological picture is different. The deformed $R^2$ model describes stronger gravity compared to the ordinary $R^2$ with an effectively smaller effective Planck mass. We propose the addition of an early dark energy term which does not affect at all the inflationary era, but takes over the control of the late-time dynamics at late times. We study in some detail the predicted dark energy era evolution, and we demonstrate that the dark energy corrected deformed $R^2$ model can describe a viable dark energy era, compatible with the Planck constraints on cosmological parameters. Furthermore the model is distinct from the $\Lambda$-Cold-Dark-Matter model, but shows a qualitatively similar behavior., Comment: PRD Accepted