Wong et al. (2018) recently performed an encouraging criticism to our paper "Gravitational waves from ultra-short period exoplanets" (Cunha, Silva, Lima 2018) exploring the potentialities of a subset of exoplanets with extremely short periods (less than 80 min) as a possible scientific target to the planned space-based LISA observatory. Here we call attention to some subtleties and limitations underlying the basic criticism which in our view were not properly stressed in their comment. Particularly, simple estimates show that a sphere encircling the Earth with a radius of 250 pc may accommodate a population $ \sim 10^{4}$ ultra-short period exoplanets with characteristic strain of the same order or higher than the ones analyzed in our paper. This means that the question related to the gravitational wave pattern of ultra-short period exoplanets may be surpassed near future by the LISA instrument with new and more definitive data.

In the last two decades, thousands of extrasolar planets were discovered based on different observational techniques, and their number must increase substantially in virtue of the ongoing and near-future approved missions and facilities. It is shown that interesting signatures of binary systems from nearby exoplanets and their parent stars can also be obtained measuring the pattern of gravitational waves that will be made available by the new generation of detectors including the space-based LISA (Laser Interferometer Space Antenna) observatory. As an example, a subset of exoplanets with extremely short periods (less than 80 min) is discussed. All of them have gravitational luminosity, $L_{GW} \sim 10^{30}$ $erg/s$, strain $h \sim 10^{-22}$, frequencies $f_{gw} > 10^{-4}$Hz, and, as such, are within the standard sensitivity curve of LISA. Our analysis suggests that the emitted gravitational wave pattern may also provide an efficient tool to discover ultra short period exoplanets.

The free parameters of a flat accelerating model without dark energy are constrained by using Supernovae type Ia and observational H(z) data. Instead of the vacuum dominance, the present accelerating stage in this modified Einstein-de Sitter cosmology is a consequence of the gravitationally-induced particle production of cold dark matter. The model present a transition from a decelerating to an accelerating regime at low redshifts, and is also able to harmonize a cold dark matter picture with the latest measurements of the Hubble constant H_0, the Supernovae observations (Constitution sample), and the H(z) data.

We present a cosmological model with a solution that self-accelerates at late-times without signs of ghost instabilities on small scales. The model is a natural extension of the Brans-Dicke (BD) theory including a non-linear derivative interaction, which appears in a theory with the Galilean shift symmetry. The existence of the self-accelerating universe requires a negative BD parameter but, thanks to the non-linear term, small fluctuations around the solution are stable on small scales. General relativity is recovered at early times and on small scales by this non-linear interaction via the Vainshtein mechanism. At late time, gravity is strongly modified and the background cosmology shows a phantom-like behaviour and the growth rate of structure formation is enhanced. Thus this model leaves distinct signatures in cosmological observations and it can be distinguished from standard $\Lambda$CDM cosmology.

We study the spectrum of gravitational perturbations around a vacuum de Sitter brane in a 5D asymmetric braneworld model, with induced curvature on the brane. This generalises the stealth acceleration model proposed by Charmousis, Gregory and Padilla (CGP) which realises the Cardassian cosmology in which power law cosmic acceleration can be driven by ordinary matter. Whenever the bulk has infinite volume we find that there is always a perturbative ghost propagating on the de Sitter brane, in contrast to the Minkowski brane case analysed by CGP. We discuss the implication of this ghost for the stealth acceleration model, and identify a limiting case where the ghost decouples as the de Sitter curvature vanishes.

A new kind of accelerating flat model with no dark energy that is fully dominated by cold dark matter (CDM) is investigated. The number of CDM particles is not conserved and the present accelerating stage is a consequence of the negative pressure describing the irreversible process of gravitational particle creation. A related work involving accelerating CDM cosmology has been discussed before the SNe observations [Lima, Abramo & Germano, Phys. Rev. D53, 4287 (1996)]. However, in order to have a transition from a decelerating to an accelerating regime at low redshifts, the matter creation rate proposed here includes a constant term of the order of the Hubble parameter. In this case, $H_0$ does not need to be small in order to solve the age problem and the transition happens even if the matter creation is negligible during the radiation and part of the matter dominated phase. Therefore, instead of the vacuum dominance at redshifts of the order of a few, the present accelerating stage in this sort of Einstein-de Sitter CDM cosmology is a consequence of the gravitational particle creation process. As an extra bonus, in the present scenario does not exist the coincidence problem that plagues models with dominance of dark energy. The model is able to harmonize a CDM picture with the present age of the universe, the latest measurements of the Hubble parameter and the Supernovae observations.

We solve for the behaviour of cosmological perturbations in the Dvali-Gabadadze-Porrati (DGP) braneworld model using a new numerical method. Unlike some other approaches in the literature, our method uses no approximations other than linear theory and is valid on large scales. We examine the behaviour of late-universe density perturbations for both the self-accelerating and normal branches of DGP cosmology. Our numerical results can form the basis of a detailed comparison between the DGP model and cosmological observations.

We study quasi-static perturbations in a cosmological background in the Dvali-Gabadadze-Porrati (DGP) braneworld model. We identify the Vainshtein radius at which the non-linear interactions of the brane bending mode become important in a cosmological background. The Vainshtein radius in the early universe is much smaller than the one in the Minkowski background, but in a self-accelerating universe it is the same as the Minkowski background. Our result shows that the perturbative approach is applicable beyond the Vainshtein radius for weak gravity by taking into account the second order effects of the brane bending mode. The linearised cosmological perturbations are shown to be smoothly matched to the solutions inside the Vainshtein radius. We emphasize the importance of imposing a regularity condition in the bulk by solving the 5D perturbations and we highlight the problem of ad hoc assumptions on the bulk gravity that lead to different conclusions.