The observation of GW170817 binary neutron star (BNS) merger event has imposed strong bounds on the speed of gravitational waves (GWs) locally, inferring that the speed of GWs propagation is equal to the speed of light. Current GW detectors in operation will not be able to observe BNS merger to long cosmological distance, where possible cosmological corrections on the cosmic expansion history are expected to play an important role, specially for investigating possible deviations from general relativity. Future GW detectors designer projects will be able to detect many coalescences of BNS at high $z$, such as the third generation of the ground GW detector called Einstein Telescope (ET) and the space-based detector deci-hertz interferometer gravitational wave observatory (DECIGO). In this paper, we relax the condition $c_T/c = 1$ to investigate modified GW propagation where the speed of GWs propagation is not necessarily equal to the speed of light. Also, we consider the possibility for the running of the Planck mass corrections on modified GW propagation. We parametrize both corrections in terms of an effective GW luminosity distance and we perform a forecast analysis using standard siren events from BNS mergers, within the sensitivity predicted for the ET and DECIGO. We find at high $z$ very strong forecast bounds on the running of the Planck mass, namely $\mathcal{O}(10^{-1})$ and $\mathcal{O}(10^{-2})$ from ET and DECIGO, respectively. Possible anomalies on GW propagation are bound to $|c_T/c - 1| \leq 10^{-2} \,\,\, (10^{-2})$ from ET (DECIGO), respectively. We finally discuss the consequences of our results on modified gravity phenomenology.

Two mechanisms of gravitational waves (GWs) emission in fast-spinning white dwarfs (WDs) are investigated: accretion of matter and magnetic deformation. In both cases, the GW emission is generated by an asymmetry around the rotation axis of the star. However, in the first case, the asymmetry is due to the amount of accreted matter on the magnetic poles, while in the second case it is due to the intense magnetic field. We have estimated the GW amplitude and luminosity for three binary systems that have a fast-spinning magnetized WD, namely, AE Aquarii, AR Scorpii and RX J0648.0-4418. We find that, for the first mechanism, the systems AE Aquarii and RX J0648.0-4418 can be observed by the space detectors BBO and DECIGO if they have an amount of accreted mass of $\delta m \geq 10^{-5}M_{\odot }$. For the second mechanism, the three systems studied require that the WD have a magnetic field above $\sim 10^{9}$ G to emit GWs that can be detected by BBO. We also verified that, in both mechanisms, the gravitational luminosity has an irrelevant contribution to the spindown luminosity of these three systems. Therefore, other mechanisms of energy emission are needed to explain the spindown of these objects.

Gravitational waves (GWs) emission due to magnetic deformation mechanism is applied for Soft Gamma Repeaters (SGRs) and Anomalous X-Ray Pulsars(AXPs), described as fast-spinning and magnetized white dwarfs (WDs). The emission is caused by the asymmetry around the rotation axis of the star generated by its own intense magnetic field. Thus, for the first time in the literature, it is estimated the GWs counterpart for SGRs/AXPs described as WD pulsars. We find that some SGRs/AXPs can be observed by the space detectors BBO and DECIGO. In particular, 1E 1547.0-5408 and SGR 1806-20 could be detected in 1 year of observation, whereas SGR 1900+14, CXOU J171405.7-381031, Swift J1834.9-0846 and SGR 1627-41 could be observed with a 5-year observation time. We also found that SGRs/AXPs as highly magnetized neutron stars are far below the sensitivity curves of BBO and DECIGO. This result indicates that a possible detection of continuous GWs originated from these objects would corroborate the WD pulsar model.

The direct detection of gravitational waves (GWs) opened a new chapter in the modern cosmology to probe possible deviations from the general relativity (GR) theory. In the present work, we investigate for the first time the modified GW form propagation from the inspiraling of compact binary systems within the context of $f(T)$ gravity in order to obtain new forecasts/constraints on the free parameter of the theory. First, we show that the modified waveform differs from the GR waveform essentially due to induced corrections on the GWs amplitude. Then, we discuss the forecasts on the $f(T)$ gravity assuming simulated sources of GWs as black hole binaries, neutron star binaries and black hole - neutron star binary systems, which emit GWs in the frequency band of the Advanced LIGO (aLIGO) interferometer and of the third generation Einstein Telescope (ET). We show that GWs sources detected within the aLIGO sensitivity can return estimates of the same order of magnitude of the current cosmological observations. On the other hand, detection within the ET sensitivity can improve by up to 2 orders of magnitude the current bound on the $f(T)$ gravity. Therefore, the statistical accuracy that can be achieved by future ground based GW observations, mainly with the ET detector (and planed detectors with a similar sensitivity), can allow strong bounds on the free parameter of the theory, and can be decisive to test the theory of gravitation.

We discuss some aspects of de Araujo, Coelho & Costa (2016, 2017) concerning the role of a time dependent magnetic ellipticity on the pulsars' braking indices and on the putative gravitational waves (GWs) these objects can emit. Since only nine of more than 2000 known pulsars have accurately measured braking indices, it is of interest to extend this study to all known pulsars, in particular as regards GW generation. In summary, our results show a pessimistic prospect for the detection of GWs generated by these pulsars, even for advanced detectors such as aLIGO and AdVirgo, and the planned Einstein Telescope, would not be able to detect these pulsar, if the ellipticity has magnetic origin.

We investigate the propagation of primordial gravitational waves within the context of the Horndeski theories, for this, we present a generalized transfer function quantifying the sub-horizon evolution of gravitational waves modes after they enter the horizon. We compare the theoretical prediction of the modified primordial gravitational waves spectral density with the aLIGO, Einstein telescope, LISA, gLISA and DECIGO sensitivity curves. Assuming reasonable and different values for the free parameters of the theory (in agreement with the event GW170817 and stability conditions of the theory), we note that the gravitational waves amplitude can vary significantly in comparison with general relativity. We find that in some cases the gravitational primordial spectrum can cross the sensitivity curves for DECIGO detector with the maximum frequency sensitivity to the theoretical predictions around 0.05 - 0.30 Hz. From our results, it is clear that the future generations of space based interferometers can bring new perspectives to probing modifications in general relativity.

The attempt to describe the recent accelerated expansion of the Universe includes different propositions for dark energy models and modified gravity theories. Establish their features in order to discriminate and even rule out part of these models using observational data is a fundamental issue of cosmology. In the present work we consider a class of extended theories of gravity (ETG) that are minimally coupled to the ordinary matter fields. In this context, and assuming a homogeneous and isotropic spacetime, we derive the energy conditions for this ETG class. We then put constraints on these conditions using a model-independent approach to reconstruct the deceleration function along with the Joint Light-curve Analysis (JLA) supernova sample, 11 baryon acoustic oscillation and 22 cosmic-chronometer measurements. We also consider an additional condition imposing the strong energy condition only on the ordinary matter. This is to guarantee the presence of only attractive matter in the energy-momentum tensor, at least in the redshift range of the observations, i.e., the recent accelerated expansion of the Universe is due solely to the modifications in the gravity theory. The main result of this work is a general reconstruction of the energy conditions valid for every considered ETG.

The presence of multifractality in a time series shows different correlations for different time scales as well as intermittent behaviour that cannot be captured by a single scaling exponent. The identification of a multifractal nature allows for a characterization of the dynamics and of the intermittency of the fluctuations in non-linear and complex systems. In this study, we search for a possible multifractal structure (multifractality signature) of the flux variability in the quasar 3C 273 time series for all electromagnetic wavebands at different observation points, and the origins for the observed multifractality. This study is intended to highlight how the scaling behaves across the different bands of the selected candidate which can be used as an additional new technique to group quasars based on the fractal signature observed in their time series and determine whether quasars are non-linear physical systems or not. The Multifractal Detrended Moving Average algorithm (MFDMA) has been used to study the scaling in non-linear, complex and dynamic systems. To achieve this goal, we applied the backward ({\theta} = 0) MFDMA method for one-dimensional signals. We observe weak multifractal (close to monofractal) behaviour in some of the time series of our candidate except in the mm, UV and X-ray bands. The non-linear temporal correlation is the main source of the observed multifractality in the time series whereas the heaviness of the distribution contributes less.

Among the various possibilities to probe the theory behind the recent accelerated expansion of the universe, the energy conditions (ECs) are of particular interest, since it is possible to confront and constrain the many models, including different theories of gravity, with observational data. In this context, we use the ECs to probe any alternative theory whose extra term acts as a cosmological constant. For this purpose, we apply a model-independent approach to reconstruct the recent expansion of the universe. Using Type Ia supernova, baryon acoustic oscillations and cosmic-chronometer data, we perform a Markov Chain Monte Carlo analysis to put constraints on the effective cosmological constant $\Omega^0_{\rm eff}$. By imposing that the cosmological constant is the only component that possibly violates the ECs, we derive lower and upper bounds for its value. For instance, we obtain that $0.59 < \Omega^0_{\rm eff} < 0.91$ and $0.40 < \Omega^0_{\rm eff} < 0.93$ within, respectively, $1\sigma$ and $3\sigma$ confidence levels. In addition, about 30\% of the posterior distribution is incompatible with a cosmological constant, showing that this method can potentially rule it out as a mechanism for the accelerated expansion. We also study the consequence of these constraints for two particular formulations of the bimetric massive gravity. Namely, we consider the Visser's theory and the Hassan and Roses's massive gravity by choosing a background metric such that both theories mimic General Relativity with a cosmological constant. Using the $\Omega^0_{\rm eff}$ observational bounds along with the upper bounds on the graviton mass we obtain constraints on the parameter spaces of both theories.

In one of our previous articles we have considered the role of a time dependent magnetic ellipticity on pulsars' braking indices and on the putative gravitational waves these objects can emit. Since only nine of more than 2000 known pulsars have accurately measured braking indices, it is of interest to extend this study to all known pulsars, in particular to what concerns the gravitational waves generation. To do so, as shown in our previous article, we need to know some pulsars' observable quantities such as: periods and their time derivatives, and estimated distances to the Earth. Moreover, we also need to know the pulsars' masses and radii, for which, in here we are adopting current fiducial values. Our results show that the gravitational wave amplitude is at best $h \sim 10^{-28}$. This leads to a pessimistic prospect for the detection of gravitational waves generated by these pulsars, even for Advanced LIGO and Advanced Virgo, and the planned Einstein Telescope, whether the ellipticity has magnetic origin.

We investigate the cosmological production of gravitational waves in a nonsingular flat cosmology powered by a "running vacuum" energy density described by $\rho_{\Lambda}\equiv\rho_{\Lambda}(H)$, a phenomenological expression potentially linked with the renormalization group approach in quantum field theory in curved spacetimes. The model can be interpreted as a particular case of the class recently discussed by Perico et al. (Phys. Rev. D {\bf 88}, 063531, 2013) which is termed complete in the sense that the cosmic evolution occurs between two extreme de Sitter stages (early and late time de Sitter phases). {The gravitational wave equation is derived and its time-dependent part numerically integrated since the primordial de Sitter stage. The generated spectrum of gravitons is also compared with the standard calculations where an abrupt transition, from the early de Sitter to the radiation phase, is usually assumed.} It is found that the stochastic background of gravitons is very similar to the one predicted by the cosmic concordance model plus inflation except at higher frequencies ($\nu \gtrsim 100$ kHz). This remarkable signature of a "running vacuum" cosmology combined with the proposed high frequency gravitational wave detectors and measurements of the CMB polarization (B-modes) may provide a new window to confront more conventional models of inflation.

We study the role of time dependent magnetic ellipticities ($\epsilon_{B}$) on the calculation of the braking index of pulsars. Moreover, we study the consequences of such a $\epsilon_{B}$ on the amplitude of gravitational waves (GWs) generated by pulsars with measured braking indices. We show that, since the ellipticity generated by the magnetic dipole is extremely small, the corresponding amplitude of GWs is much smaller than the amplitude obtained via the spindown limit.

We study the putative emission of gravitational waves (GWs) in particular for pulsars with measured braking index. We show that the appropriate combination of both GW emission and magnetic dipole brakes can naturally explain the measured braking index, when the surface magnetic field and the angle between the magnetic dipole and rotation axes are time dependent. Then we discuss the detectability of these very pulsars by aLIGO and the Einstein Telescope. We call attention to the realistic possibility that aLIGO can detect the GWs generated by at least some of these pulsars, such as Vela, for example.

The geosynchronous Laser Interferometer Space Antenna (gLISA) is a space-based gravitational wave (GW) mission that, for the past five years, has been under joint study at the Jet Propulsion Laboratory, Stanford University, the National Institute for Space Research (I.N.P.E., Brazil), and Space Systems Loral. If flown at the same time as the LISA mission, the two arrays will deliver a joint sensitivity that accounts for the best performance of both missions in their respective parts of the mHz band. This simultaneous operation will result in an optimally combined sensitivity curve that is "white" from a few mHz to 1 Hz, making the two antennas capable of detecting, with high signal-to-noise ratios (SNRs), coalescing black-hole binaries (BHBs) with masses in the range (10 - 100 million) solar masses. Their ability of jointly tracking, with enhanced SNR, signals similar to that observed by the Advanced Laser Interferometer Gravitational Wave Observatory (aLIGO) on September 14, 2015 (the GW150914 event) will result in a larger number of observable small-mass binary black-holes and an improved precision of the parameters characterizing these sources. Together, LISA, gLISA and aLIGO will cover, with good sensitivity, a frequency band from a tenth of a mHz to a kHz frequency band.

Recently, a braking index for the pulsar PSR J1640-4631 has been measured. With a braking index of $n = 3.15 \pm 0.03$, this pulsar has the highest braking index ever measured. As it is well known, a pure magnetic dipole brake yields $n = 3$, whereas a pure gravitational wave brake yields $n = 5$. Therefore, each of these mechanisms alone can not account for the braking index found for PSR J1640-4631. Here we show that such a braking index can be accounted for if the spindown model combines magnetic dipole and gravitational wave brakes. Then, we briefly discuss the detectability of this pulsar by aLIGO and the planned Einstein Telescope. In particular, we show that the amplitude of the gravitational wave that comes from our model is around a factor four lower than the amplitude modeled exclusively by gravitational wave energy loss. Another interesting outcome of our modeling is that we are able to obtain the ellipticity from the braking index and other pulsar parameters.

A study of binary systems composed of two point particles with different masses in the linear regime of the characteristic formulation of general relativity in a Minkowski background is provided. The present paper generalizes a previous study find in the literature. The boundary conditions at the world tubes generated by the particle's orbits are explored, when the metric variables are decomposed in spin-weighted spherical harmonics. The power lost by the emission of gravitational waves is computed using the News Bondi's functions. The power found is nothing but the well-known result for the power obtained by Peters and Mathews in a different approach. This agreement shows the consistency of the approach considered here. The contribution to the gravitational radiation of several multipole terms is also shown.

We firstly revisit the energy loss mechanism known as quantum vacuum friction (QVF), clarifying some of its subtleties. Then we investigate the observables that could easily differentiate QVF from the classical magnetic dipole radiation for pulsars with braking indices (n) measured accurately. We show this is specially the case for the time evolution of a pulsar's magnetic dipole direction ($\dot{\phi}$) and surface magnetic field ($\dot{B}_0$). As it is well known in the context of the classic magnetic dipole radiation, $n<3$ would only be possible for positive $(\dot{B}_0/B_0 + \dot{\phi}/\tan\phi)$, which, for instance, leads to $\dot{B}_0>0$ ($\dot{\phi}>0$) when $\phi$ ($B_0$) is constant. On the other hand, we show that QVF can result in very contrasting predictions with respect to the above ones. Finally, even in the case $\dot{B}_0$ in both aforesaid models for a pulsar has the same sign, for a given $\phi$, we show that they give rise to different associated timescales, which could be another way to falsify QVF.

Gravitational waves are perturbations in the spacetime that propagate at the speed of light. The study of such phenomenon is interesting because many cosmological processes and astrophysical objects, such as binary systems, are potential sources of gravitational radiation and can have their emissions detected in the near future by the next generation of interferometric detectors. Concerning the astrophysical objects, an interesting case is when there are several sources emitting in such a way that there is a superposition of signals, resulting in a smooth spectrum which spans a wide range of frequencies, the so-called stochastic background. In this paper, we are concerned with the stochastic backgrounds generated by compact binaries (i.e. binary systems formed by neutron stars and black holes) in the coalescing phase. In particular, we obtain such backgrounds by employing a new method developed in our previous studies.

Binary Systems are the most studied sources of gravitational waves. The mechanisms of emission and the behavior of the orbital parameters are well known and can be written in analytic form in several cases. Besides, the strongest indication of the existence of gravitational waves has arisen from the observation of binary systems. On the other hand, when the detection of gravitational radiation becomes a reality, one of the observed pattern of the signals will be probably of stochastic background nature, which are characterized by a superposition of signals emitted by many sources around the universe. Our aim here is to develop an alternative method of calculating such backgrounds emitted by cosmological compact binary systems during their periodic or quasiperiodic phases. We use an analogy with a problem of Statistical Mechanics in order to perform this sum as well as taking into account the temporal variation of the orbital parameters of the systems. Such a kind of background is of particular importance since it could well form an important foreground for the planned gravitational wave interferometers DECI-Hertz Interferometer Gravitational wave Observatory (DECIGO), Big Bang Observer (BBO), Laser Interferometer Space Antenna (LISA) or Evolved LISA (eLISA), Advanced Laser Interferometer Gravitational-Wave Observatory (ALIGO) and Einstein Telescope (ET).

In the study of gravitational waves (GWs), the stochastic background generated by compact binary systems are among the most important kinds of signals. The reason for such an importance has to do with their probable detection by the interferometric detectors [such as the Advanced LIGO (ALIGO) and Einstein Telescope (ET)] in the near future. In this paper we are concerned with, in particular, the stochastic background of GWs generated by double neutron star (DNS) systems in circular orbits during their periodic and quasi--periodic phases. Our aim here is to describe a new method to calculate such spectra, which is based on an analogy with a problem of Statistical Mechanics. Besides, an important characteristic of our method is to consider the time evolution of the orbital parameters.

Binary systems emit gravitational waves in a well-known pattern; for binaries in circular orbits, the emitted radiation has a frequency that is twice the orbital frequency. Systems in eccentric orbits, however, emit gravitational radiation in the higher harmonics too. In this paper, we are concerned with the stochastic background of gravitational waves generated by double neutron star systems of cosmological origin in eccentric orbits. We consider in particular the long-lived systems, that is, those binaries for which the time to coalescence is longer than the Hubble time ($\sim 10$Gyr). Thus, we consider double neutron stars with orbital frequencies ranging from $10^{-8}$ to $2\times 10^{-6}$Hz. Although in the literature some papers consider the spectra generated by eccentric binaries, there is still space for alternative approaches for the calculation of the backgrounds. In this paper, we use a method that consists in summing the spectra that would be generated by each harmonic separately in order to obtain the total background. This method allows us to clearly obtain the influence of each harmonic on the spectra. In addition, we consider different distribution functions for the eccentricities in order to investigate their effects on the background of gravitational waves generated. At last, we briefly discuss the detectability of this background by space-based gravitational wave antennas and pulsar timing arrays.

We analyze the trajectories of three geostationary satellites forming the GEOstationary GRAvitational Wave Interferometer (GEOGRAWI)~\cite{tinto}, a space-based laser interferometer mission aiming to detect and study gravitational radiation in the ($10^{-4} - 10$) Hz band. The combined effects of the gravity fields of the Earth, the Sun and the Moon make the three satellites deviate from their nominally stationary, equatorial and equilateral configuration. Since changes in the satellites relative distances and orientations could negatively affect the precision of the laser heterodyne measurements, we have derived the time-dependence of the inter-satellite distances and velocities, the variations of the polar angles made by the constellation's three arms with respect to a chosen reference frame, and the time changes of the triangle's enclosed angles. We find that, during the time between two consecutive station-keeping maneuvers (about two weeks), the relative variations of the inter-satellite distances do not exceed a value of $0.05$ percent, while the relative velocities between pairs of satellites remain smaller than about $0.7 \ {\rm m/s}$. In addition, we find the angles made by the arms of the triangle with the equatorial plane to be periodic functions of time whose amplitudes grow linearly with time; the maximum variations experienced by these angles as well as by those within the triangle remain smaller than $3$ arc-minutes, while the East-West angular variations of the three arms remain smaller than about $15$ arc-minutes during the two-weeks period. The relatively small variations of these orbit parameters result into a set of system functional and performance requirements that are less stringent than those characterizing an interplanetary mission.

We analyze the sensitivities of a geostationary gravitational wave interferometer mission operating in the sub-Hertz band. Because of its smaller armlength, in the lower part of its accessible frequency band ($10^{-4} - 2 \times 10^{-2}$ Hz) our proposed Earth-orbiting detector will be less sensitive, by a factor of about seventy, than the Laser Interferometer Space Antenna (LISA) mission. In the higher part of its band instead ($2 \times 10^{-2} - 10$ Hz), our proposed interferometer will have the capability of observing super-massive black holes (SMBHs) with masses smaller than $\sim 10^{6}$ M$_{\odot}$. With good event rates for these systems, a geostationary interferometer will be able to accurately probe the astrophysical scenarios that account for their formation.

It is believed that soon after the Planck era, spacetime should have a semi-classical nature. Therefore, it is unavoidable to modify the theory of General Relativity or look for alternative theories of gravitation. An interesting possibility found in the literature considers two geometric counter-terms to regularize the divergences of the effective action. These counter-terms are responsible for a higher order derivative metric theory of gravitation. In the present letter we investigate how isotropisation occurs. For this reason a single solution is chosen throughout this article. We obtain perturbatively, by two different methods, that the tensor and scalar components emerge naturally during the isotropisation process. In this sense our result provides a numerical example to Stelle's well known result on classical gravity with higher derivates. Our entire analysis is restricted to the particular Bianchi type $I$ case.

It is presented a brief review on gravitational waves (GWs). It is shown how the wave equation is obtained from Einstein's equations and how many and how are the polarization modes of these waves. It is discussed the reasons why GWs sources should be of astrophysical or cosmological origin. Thus, it is discussed what would be the most likely sources of GWs to be detected by the detectors of GWs currently in operation and those that should be operational in the future, emphasizing in particular the sources involving compact objects. The compact objects such as neutron stars, black holes and binary systems involving compact stars can be important sources of GWs. Last but not least, it is discussed the GWs astrophysics that is already possible to do, in particular involving the compact objects.

A way to probe alternative theories of gravitation is to study if they could account for the structures of the universe. We then modified the well-known Gadget-2 code to probe alternative theories of gravitation through galactic dynamics. As an application, we simulate the evolution of spiral galaxies to probe alternative theories of gravitation whose weak field limits have a Yukawa-like gravitational potential. These simulations show that galactic dynamics can be used to constrain the parameters associated with alternative theories of gravitation. It is worth stressing that the recipe given in the present study can be applied to any other alternative theory of gravitation in which the superposition principle is valid.

We propose a Geostationary Gravitational Wave Interferometer (GEOGRAWI) mission concept for making observations in the sub-Hertz band. GEOGRAWI is expected to meet some of LISA's science goals in the lower part of its accessible frequency band ($10^{-4} - 2 \times 10^{-2}$ Hz), and to outperform them by a large margin in the higher-part of it ($2 \times 10^{-2} - 10$ Hz). As a consequence of its Earth-bound orbit, GEOGRAWI is significantly less expensive than the interplanetary LISA mission and could be either an entirely US mission or managed and operated by NASA in partnership with the Brazilian Space Agency.

The coalescence history of massive black holes has been derived from cosmological simulations, in which the evolution of those objects and that of the host galaxies are followed in a consistent way. The present study indicates that supermassive black holes having masses greater than $\sim 10^{9} M_{\odot}$ underwent up to 500 merger events along their history. The derived coalescence rate per comoving volume and per mass interval permitted to obtain an estimate of the expected detection rate distribution of gravitational wave signals ("ring-down") along frequencies accessible by the planned interferometers either in space (LISA) or in the ground (Einstein). For LISA, in its original configuration, a total detection rate of about $15 yr^{-1}$ is predicted for events having a signal-to-noise ratio equal to 10, expected to occur mainly in the frequency range $4-9 mHz$. For the Einstein gravitational wave telescope, one event each 14 months down to one event each 4 years is expected with a signal-to-noise ratio of 5, occurring mainly in the frequency interval $10-20 Hz$. The detection of these gravitational signals and their distribution in frequency would be in the future an important tool able to discriminate among different scenarios explaining the origin of supermassive black holes.

In this work, we explore some cosmological implications of the model proposed by M. Visser in 1998. In his approach, Visser intends to take in account mass for the graviton by means of an additional bimetric tensor in the Einstein's field equations. Our study has shown that a consistent cosmological model arises from Visser's approach. The most interesting feature is that an accelerated expansion phase naturally emerges from the cosmological model, and we do not need to postulate any kind of dark energy to explain the current observational data for distant type Ia supernovae (SNIa).

Theories of gravity for which gravitons can be treated as massive particles have presently been studied as realistic modifications of General Relativity, and can be tested with cosmological observations. In this work, we study the ability of a recently proposed theory with massive gravitons, the so-called Visser theory, to explain the measurements of luminosity distance from the Union2 compilation, the most recent Type-Ia Supernovae (SNe Ia) dataset, adopting the current ratio of the total density of non-relativistic matter to the critical density ($\Omega_m$) as a free parameter. We also combine the SNe Ia data with constraints from Baryon Acoustic Oscillations (BAO) and CMB measurements. We find that, for the allowed interval of values for $\Omega_m$, a model based on Visser's theory can produce an accelerated expansion period without any dark energy component, but the combined analysis (SNe Ia + BAO + CMB) shows that the model is disfavored when compared with $\Lambda$CDM model.

In the standard scenario of the Newtonian gravity, a late-type galaxy (i.e., a spiral galaxy) is well described by a disk and a bulge embedded in a halo mainly composed by dark matter. In Brownstein-Moffat gravity, there is a claim that late-type galaxy systems would not need to have halos, avoiding as a result the dark matter problem, i.e., a modified gravity (non-Newtonian) would account for the galactic structure with no need of dark matter. In the present paper, we probe this claim via numerical simulations. Instead of using a "static galaxy," where the centrifugal equilibrium is usually adopted, we probe the Brownstein-Moffat gravity dynamically via numerical $N$-body simulations.

Cosmological Gravitational Waves (GWs) are usually associated with the transverse-traceless part of the metric perturbations in the context of the theory of cosmological perturbations. These modes are just the usual polarizations `+' and `x' which appear in the general relativity theory. However, in the majority of the alternative theories of gravity, GWs can present more than these two polarization states. In this context, the Newman-Penrose formalism is particularly suitable for evaluating the number of non-null GW modes. In the present work we intend to take into account these extra polarization states for cosmological GWs in alternative theories of gravity. As an application, we derive the dynamical equations for cosmological GWs for two specific theories, namely, a general scalar-tensor theory which presents four polarization states and a massive bimetric theory which is in the most general case with six polarization states for GWs. The mathematical tool presented here is quite general, so it can be used to study cosmological perturbations in all metric theories of gravity.

We study the generation of a stochastic gravitational wave (GW) background produced by a population of neutron stars (NSs) which go over a hadron-quark phase transition in its inner shells. We obtain, for example, that the NS phase transition, in cold dark matter scenarios, could generate a stochastic GW background with a maximum amplitude of $h_{\rm BG} \sim 10^{-24}$, in the frequency band $\simeq 20-2000 {\rm Hz}$ for stars forming at redshifts of up to $z\simeq 20.$ We study the possibility of detection of this isotropic GW background by correlating signals of a pair of `advanced' LIGO observatories.

Since the Newtonian gravitation is largely used to model with success the structures of the universe, such as galaxies and clusters of galaxies, for example, a way to probe and constrain alternative theories, in the weak field limit, is to apply them to model the structures of the universe. We then modified the well known Gadget-2 code to probe alternative theories of gravitation through galactic dynamics. In particular, we modified the Gadget-2 code to probe alternatives theories whose weak field limits have a Yukawa-like gravitational potential. As a first application of this modified Gadget-2 code we simulate the evolution of elliptical galaxies. These simulations show that galactic dynamics can be used to constrain the parameters associated with alternative theories of gravitation.

We present studies of the collapse of neutron stars that undergo a hadron-quark phase transition. A spherical Lagrangian hydrodynamic code has been written. As initial condition we take different neutron star configurations taking into account its density, energy density and pressure distribution. The phase transition is imposed at different evolution times. We have found that a significant amount of matter on the surface can be ejected while the remaining star rings in the fundamental and first pressure modes.

In the weak field limit general relativity reduces, as is well known, to the Newtonian gravitation. Alternative theories of gravity, however, do not necessarily reduce to Newtonian gravitation; some of them, for example, reduce to Yukawa-like potentials instead of the Newtonian potential. Since the Newtonian gravitation is largely used to model with success the structures of the universe, such as for example galaxies and clusters of galaxies, a way to probe and constrain alternative theories, in the weak field limit, is to apply them to model the structures of the universe. In the present study, we consider how to probe Yukawa-like potentials using N-body numerical simulations.

We study the generation of a stochastic gravitational wave (GW) background produced by a population of neutron stars (NSs) which go over a hadron-quark phase transition in its inner shells. We obtain, for example, that the NS phase transition, in cold dark matter scenarios, could generate a stochastic GW background with a maximum amplitude of $h_{\rm BG} \sim 10^{-24}$, in the frequency band $\nu_{\rm{obs}} \simeq 20-2000 {\rm Hz}$ for stars forming at redshifts of up to $z\simeq 20.$ We study the possibility of detection of this isotropic GW background by correlating signals of a pair of Advanced LIGO observatories.

The direct observation of gravitational waves (GWs) in the near future, and the corresponding determination of the number of independent polarizations, is a powerful tool to test general relativity and alternative theories of gravity. In the present work we use the Newman-Penrose formalism to characterize GWs in quadratic gravity and in a particular class of f(R) Lagrangians. We find that both quadratic gravity and the f(R) theory belong to the most general invariant class of GWs, i.e., they can present up to six independent polarizations of GWs. For a particular combination of the parameters, we find that quadratic gravity can present up to five polarizations states. On the other hand, if we use the Palatini approach for f(R) theories, GWs present only the usual two transverse-traceless polarizations such as in general relativity. Thus, we conclude that the observation of GWs can strongly constrain the suitable formalism for these theories.

In this work we give special attention to the bimetric theory of gravitation with massive gravitons proposed by Visser in 1998. In his theory, a prior background metric is necessary to take in account the massive term. Although in the great part of the astrophysical studies the Minkowski metric is the best choice to the background metric, it is not possible to consider this metric in cosmology. In order to keep the Minkowski metric as background in this case, we suggest an interpretation of the energy-momentum conservation in Visser's theory, which is in accordance with the equivalence principle and recovers naturally the special relativity in the absence of gravitational sources. Although we do not present a general proof of our hypothesis we show its validity in the simple case of a plane and dust-dominated universe, in which the `massive term' appears like an extra contribution for the energy density.

We present a new solution for the rotation curves of galactic disks with gravitational potential of the Yukawa type. We follow the technique employed by Toomre in 1963 in the study of galactic disks in the Newtonian theory. This new solution allows an easy comparison between the Newtonian solution and the Yukawian one. Therefore, constraints on the parameters of theories of gravitation can be imposed, which in the weak field limit reduce to Yukawian potentials. We then apply our formulae to the study of rotation curves for a zero-thickness exponential disk and compare it with the Newtonian case studied by Freeman in 1970. As an application of the mathematical tool developed here, we show that in any theory of gravity with a massive graviton (this means a gravitational potential of the Yukawa type), a strong limit can be imposed on the mass (m_g) of this particle. For example, in order to obtain a galactic disk with a scale length of b ~ 10 kpc, we should have a massive graviton of m_g << 10^{-59} g. This result is much more restrictive than those inferred from solar system observations.

We briefly review the properties of quasi-normal modes of neutron stars and black holes. We analyse the consequences of a possible detection of such modes via the gravitational waves associated with them, especially addressing our study to the Brazilian spherical antenna, on which a possible detection would occur at 3.0-3.4 kHz. A question related to any putative gravitational wave detection concerns the source that produces it. We argue that, since the characteristic damping times for the gravitational waves of neutron stars and black holes are different, a detection can distinguish between them, and also distinguish the neutron stars oscillating modes. Moreover, since the source can be identified by its characteristic damping time, we are able to extract information about the neutron star or black hole. This information would lead, for example, to a strong constraint in the nuclear matter equation of state, namely the compression modulus should be K=220 MeV.

There is in the literature a number of papers addressing the stochastic background of gravitational waves (GWs) generated by an ensemble of astrophysical sources. The main ingredient in such studies is the so called star formation rate density (SFRD), which gives the density of stars formed per unit time. Some authors argue, however, that there is, in the equation that determines the amplitude of the stochastic background of GWs, an additional (1+z) term dividing the SFRD, which would account for the effect of cosmic expansion onto the time variable. We argue here that the inclusion of this additional term is wrong. In order to clarify where the inclusion of the (1+z) term is really necessary, we briefly discuss the calculation of event rates in the study of GRBs (gamma ray bursts) from cosmological origin.

In this paper we investigate the epochs in which the Universe started accelerating and when it began to become dark energy-dominated (i.e., the dynamics of the expansion of the Universe dominated by the dark energy). We provide analytic expressions to calculate the redshifts of these epochs as a function of density parameters. Moreover, we review and discuss cosmological models with a dark energy component, which can have an interesting characteristic, namely, they never stop accelerating. This holds even if the Universe is at present time either flat, open, or closed. If the dark energy is the cosmological constant the Universe will eventually end up undergoing an exponentially expansion phase, and the total density parameter converging to $\Omega=1$. This is exactly what is considered in inflationary scenario to generate the initial conditions for the big bang. One can then argue that the Universe begun with an inflationary phase and will end up with another inflationary phase. Thus, it follows that in both the early and the late Universe $\Omega \to 1$. We also discuss the above issues in the context of the XCDM parametrization.

Different studies show that dark matter of non-baryonic origin might exist. There have been experimental evidences that at least one form of dark matter has been detected through microlensing effects. This form of dark matter is named MACHOs (Massive Astrophysical Compact Halo Objects). The MACHO collaboration estimated that the masses of these objects are to be in the range 0.15-0.95 solar masses, where the most probable mass is of 0.5 solar masses. Some authors argue that MACHOs could be black holes, and that they could form binary systems, BHMACHO binaries. As is well known binary systems are sources of gravitational waves. The Brazilian spherical antenna will operate in the frequency band of 3.0-3.4 kHz, sensitive to binaries of a pair of 0.5 solar mass black holes just before coalescing. In the present work we study the detectability of these putative BHMACHO binaries by the Brazilian spherical antenna Mario Schenberg.

The Brazilian spherical antenna (Schenberg) is planned to detect high frequency gravitational waves (GWs) ranging from 3.0 kHz to 3.4 kHz. There is a host of astrophysical sources capable of being detected by the Brazilian antenna, namely: core collapse in supernova events; (proto)neutron stars undergoing hydrodynamical instability; f-mode unstable neutron stars, caused by quakes and oscillations; excitation of the first quadrupole normal mode of 4-9 solar mass black holes; coalescence of neutron stars and/or black holes; exotic sources such as bosonic or strange matter stars rotating at 1.6 kHz; and inspiralling of mini black hole binaries. We here address our study in particular to the neutron stars, which could well become f-mode unstable producing therefore GWs. We estimate, for this particular source of GWs, the event rates that in principle can be detected by Schenberg and by the Dutch Mini-Grail antenna.

We study the generation of a gravitational wave (GW) background produced from a population of core-collapse supernovae, which form black holes in scenarios of structure formation of the Universe. We obtain, for example, that a pre-galactic population of black holes, formed at redshifts $z\simeq 30-10$, could generate a stochastic GW background with a maximum amplitude of $h_{\rm BG}\simeq 10^{-24}$ in the frequency band $\nu_{\rm obs}\simeq 30-470 {\rm Hz}$ (considering a maximum efficiency of generation of GWs, namely, $\epsilon_{\rm GW}=7\times 10^{-4}$). In particular, we discuss what astrophysical information could be obtained from a positive, or even a negative, detection of such a GW background produced in scenarios such as those studied here. One of them is the possibility of obtaining the initial and final redshifts of the emission period from the observed spectrum of GWs.

We study the generation of a stochastic gravitational wave (GW) background produced from a population of core-collapse supernovae, which form black holes in scenarios of structure formation. We obtain, for example, that the formation of a population (Population III) of black holes, in cold dark matter scenarios, could generate a stochastic GW background with a maximum amplitude of $h_{\rm BG} \simeq 10^{-24}$ and corresponding closure energy density of $\Omega_{\rm{GW}}\sim 10^{-7}$, in the frequency band $\nu_{\rm{obs}} \simeq 30-470 {\rm Hz}$ (assuming a maximum efficiency of generation of GWs, namely, $\epsilon_{\rm GW_{\rm max}} = 7\times 10^{-4}$) for stars forming at redshifts $z\simeq 30-10.$ We show that it will be possible in the future to detect this isotropic GW background by correlating signals of a pair of `advanced' LIGO observatories (LIGO III) at a signal-to-noise ratio of $\simeq 40$. We discuss what astrophysical information could be obtained from a positive (or even a negative) detection of such a GW background generated in scenarios such as those studied here. One of them is the possibility of obtaining the initial and final redshifts of the emission period from the observed spectrum of GWs.

A spherical Lagrangian hydrodynamical code has been written to study the formation of cosmological structures in the early Universe. In this code we take into account the presence of collisionless non-baryonic cold dark matter (CDM), the cosmological constant and a series of physical processes present during and after the recombination era, such as photon drag resulting from the cosmic background radiation and hydrogen molecular production. We follow the evolution of the structure since the recombination era until the present epoch. As an application of this code we study the formation of voids starting from negative density perturbations which evolved during and after the recombination era. We analyse a set of COBE-normalized models, using different spectra to see their influence on the formation of voids. Our results show that large voids with diameters ranging from 10h^{-1} Mpc up to 50h^{-1} Mpc can be formed in a universe model dominated by the cosmological constant (\Omega_\Lambda ~ 0.8). This particular scenario is capable of forming large and deep empty regions (with density contrasts \delta < -0.6). Our results also show that the physical processes acting on the baryonic matter produce a transition region where the radius of the dark matter component is greater than the baryonic void radius. The thickness of this transition region ranges from about tens of kiloparsecs up to a few megaparsecs, depending on the spectrum considered. Putative objects formed near voids and within the transition region would have a different amount of baryonic/dark matter when compared with \Omega_b/\Omega_d. If one were to use these galaxies to determine, by dynamical effects or other techniques, the quantity of dark matter present in the Universe, the result obtained would be only local and not representative of the Universe as a whole.

Recently Fuller & Shi proposed that the gravitational collapse of supermassive objects ($M \gtrsim 10^4M_\odot$) could be a cosmological source of $\gamma$-ray bursts (GRBs). The major advantage of their model is that supermassive object collapses are far more energetic than solar mass-scale compact mergers. Also, in their proposal the seeds of supermassive black holes (SMBHs) thus formed could give rise to the SMBHs observed at the center of many galaxies. We argue here that, besides the generation of GRBs, there could well occur a strong generation of gravitational waves (GWs) during the formation of SMBHs. As a result, the rate of such GW bursts could be as high as the rate of GRBs in the model by Fuller & Shi. In this case, the detection of GRBs and bursts of GWs should occur with a small time difference. We also argue that the GWs produced by the SMBHs studied here could be detected when the Laser Interferometric Space Antenna (LISA) becomes operative.

A solution to the Einstein field equations that represents a rigidly rotating dust accompanied by a thin matter shell of the same type is found.

General Relativity predicts that binary systems of stars produce gravitational waves of significant intensity. Here we are particularly interested in the cataclysmic variable binaries (CVs). These systems emit low frequency gravitational waves, $ f < 10^{-3} Hz$. We present here a catalog of CVs and argue that part of them are capable of being detected by the Laser Interferometer Space Antenna (LISA).

A continuous stochastic background of gravitational waves (GWs) for burst sources is produced if the mean time interval between the occurrence of bursts is smaller than the average time duration of a single burst at the emission, i.e., the so called duty cycle must be greater than one. To evaluate the background of GWs produced by an ensemble of sources, during their formation, for example, one needs to know the average energy flux emitted during the formation of a single object and the formation rate of such objects as well. In many cases the energy flux emitted during an event of production of GWs is not known in detail, only characteristic values for the dimensionless amplitude and frequencies are known. Here we present a shortcut to calculate stochastic backgrounds of GWs produced from cosmological sources. For this approach it is not necessary to know in detail the energy flux emitted at each frequency. Knowing the characteristic values for the ``lumped'' dimensionless amplitude and frequency we show that it is possible to calculate the stochastic background of GWs produced by an ensemble of sources.

In this article we extend the study performed in our previous article on the collapse of primordial objects. We here analyze the behavior of the physical parameters for clouds ranging from $10^7M_\odot$ to $10^{15}M_\odot$. We studied the dynamical evolution of these clouds in two ways: purely baryonic clouds and clouds with non-baryonic dark matter included. We start the calculations at the beginning of the recombination era, following the evolution of the structure until the collapse (that we defined as the time when the density contrast of the baryonic matter is greater than $10^4$). We analyze the behavior of the several physical parameters of the clouds (as, e.g., the density contrast and the velocities of the baryonic matter and the dark matter) as a function of time and radial position in the cloud. In this study all physical processes that are relevant to the dynamical evolution of the primordial clouds, as for example photon-drag (due to the cosmic background radiation), hydrogen molecular production, besides the expansion of the Universe, are included in the calculations. In particular we find that the clouds, with dark matter, collapse at higher redshift when we compare the results with the purely baryonic models. As a general result we find that the distribution of the non-baryonic dark matter is more concentrated than the baryonic one. It is important to stress that we do not take into account the putative virialization of the non-baryonic dark matter, we just follow the time and spatial evolution of the cloud solving its hydrodynamical equations. We studied also the role of the cooling-heating processes in the purely baryonic clouds.

We present here studies of collapse of purely baryonic Population III objects with masses ranging from $10M_\odot$ to $10^6M_\odot$. A spherical Lagrangian hydrodynamic code has been written to study the formation and evolution of the primordial clouds, from the beginning of the recombination era ($z_{rec} \sim 1500$) until the redshift when the collapse occurs. All the relevant processes are included in the calculations, as well as, the expansion of the Universe. As initial condition we take different values for the Hubble constant and for the baryonic density parameter (considering however a purely baryonic Universe), as well as different density perturbation spectra, in order to see their influence on the behavior of the Population III objects evolution. We find, for example, that the first mass that collapses is $8.5\times10^4M_\odot$ for $h=1$, $\Omega=0.1$ and $\delta_i={\delta\rho / \rho}=(M / M_o)^{-1/3}(1+z_{rec})^{-1}$ with the mass scale $M_o=10^{15}M_\odot$. For $M_o=4\times10^{17}M_\odot$ we obtain $4.4\times10^{4}M_\odot$ for the first mass that collapses. The cooling-heating and photon drag processes have a key role in the collapse of the clouds and in their thermal history. Our results show, for example, that when we disregard the Compton cooling-heating, the collapse of the objects with masses $>8.5\times10^4M_\odot$ occurs earlier. On the other hand, disregarding the photon drag process, the collapse occurs at a higher redshift.

In this letter we suggest a scenario for simultaneous emission of gravitational-wave and $\gamma$-ray bursts (GRBs) from soft gamma-ray repeaters (SGRs). we argue that both of the radiations can be generated by a super-Eddington accreting neutron stars in X-ray binaries. In this model a supercritical accretion transient takes back onto the remnant star the disk leftover by the hydrodynamic instability phase of a low magnetized, rapidly rotating neutron star in a X-ray binary system. We estimate the rise timescale $\Delta t_c = 0.21 ms$, minimum mass accretion rate needed to trigger the $\gamma$-ray emission, $\dot{M}_\lambda = 4.5 \times 10^{28} g$, and its effective associated temperature $T_{eff} = 740 keV$, and the timescale for repeating a burst of $\gamma$-rays $\Delta \tau_R = 11.3 yr$. Altogether, we find the associated GW amplitude and frequency to be $h_c = 2.7 \times 10^{-23}/{(Hz)}^{1/2}$ and $f_{gw} = 966 Hz$, for a source distance $\sim 55 kpc$. Detectability of the pulses by t he forthcoming GW anntenas is discussed and found likely.

The possible existence of a primordial magnetic field in the universe has been previously investigated in many articles. Studies involving the influence of a magnetic field in the nucleosyntesis era, studies considering the effects in the formation of structures during the radiation era and the matter era have been considered. We here assume the existence of a primordial magnetic field and study its effect, in particular, in the formation of voids. The study is twofold: to put constraints on the strength of the magnetic field during the recombination era and to preview its effects on the formation of voids.

We have developed a numerical code to study the deformation ($\varepsilon = (I_{zz}-I_{xx})/I_{zz}$, where $I_{ii}$ are the moments of inertia) of neutron stars in rapidly rotation in a fully general relativistic calculation. We have found that the deformation is larger, depending on the angular velocity, than is generally assumed for gravitational wave estimations. Calculations were performed by employing the Bethe-Johnson I EOS (equation of state) and a new set of models by the Frankfurt group including $\Lambda$ hyperons for several choices of their coupling constants to ordinary nucleons. Possible implications for gravitational wave searches are briefly discussed.

In this work, we analyse the temperature fluctuations of the cosmic microwave background radiation observed by COBE and show that the distribution can be fitted by a fractal distribution with a fractal dimension $ D= 1.43 \pm 0.07 $. This value is in close agreement with the fractal dimension obtained by Coleman and Pietronero (1992) and Luo and Schramm (1992) from galaxy-galaxy and cluster-cluster correlations up to $ \sim 100 h^{-1} Mpc$. The fact that the observed temperature fluctuations correspond to scales much larger than $ 100 h^{-1} Mpc$ and are signatures of the primordial density fluctuations at the recombination layer suggests that the structure of the matter at the early universe was already fractal and thus non-homogeneous on those scales. This result may have important consequences for the theoretical framework that describes the universe.