SEARCH · Engineering Papers

Results for “stochastic processes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,513 records · Page 84

All sources NASA NTRS DOE OSTI

DOE OSTI2015

SOLAR INTERACTING PROTONS VERSUS INTERPLANETARY PROTONS IN THE CORE PLUS HALO MODEL OF DIFFUSIVE SHOCK ACCELERATION AND STOCHASTIC RE-ACCELERATION

With the first observations of solar γ-rays from the decay of pions, the relationship of protons producing ground level enhancements (GLEs) on the Earth to those of similar energies producing the γ-rays on the Sun has been debated. These two populations may be either independent and simply coincident in large flares, or they may be, in fact, the same population stemming from a single accelerating agent and jointly distributed at the Sun and also in space. Assuming the latter, we model a scenario in which particles are accelerated near the Sun in a shock wave with a fraction transported back to the solar surface to radiate, while the remainder is detected at Earth in the form of a GLE. Interplanetary ions versus ions interacting at the Sun are studied for a spherical shock wave propagating in a radial magnetic field through a highly turbulent radial ray (the acceleration core) and surrounding weakly turbulent sector in which the accelerated particles can propagate toward or away from the Sun. The model presented here accounts for both the first-order Fermi acceleration at the shock front and the second-order, stochastic re-acceleration by the turbulence enhanced behind the shock. We find that the re-acceleration is important in generating the γ-radiation and we also find that up to 10% of the particle population can find its way to the Sun as compared to particles escaping to the interplanetary space.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Neutrino and cosmic-ray emission and cumulative background from radiatively inefficient accretion flows in low-luminosity active galactic nuclei

We study high-energy neutrino and cosmic-ray (CR) emission from the cores of low-luminosity active galactic nuclei (LLAGN). In LLAGN, the thermalization of particles is expected to be incomplete in radiatively inefficient accretion flows (RIAF), allowing the existence of non-thermal particles. In this work, assuming stochastic particle acceleration due to turbulence in RIAFs, we solve the Fokker–Planck equation and calculate spectra of escaping neutrinos and CRs. The RIAF in LLAGN can emit CR protons with ≳10 PeV energies and TeV–PeV neutrinos generated via pp and/or pγ reactions. We find that, if ∼1% of the accretion luminosity is carried away by non-thermal ions, the diffuse neutrino intensity from the cores of LLAGN may be as high as E{sub ν}{sup 2}Φ{sub ν}∼3×10{sup −8} GeV cm{sup −2} s{sup −1} sr{sup −1}, which can be compatible with the observed IceCube data. This result does not contradict either of the diffuse gamma-ray background observed by Fermi or observed diffuse CR flux. Our model suggests that, although very-high-energy gamma-rays may not escape, radio-quiet active galactic nuclei with RIAFs can emit GeV gamma-rays, which could be used for testing the model. We also calculate the neutron luminosity from RIAFs of LLAGN, and discuss a strong constraint on the model of jet mass loading mediated by neutrons from the diffuse neutrino observation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Simultaneous multi-band radio and X-ray observations of the galactic center magnetar SGR 1745–2900

We report on multi-frequency, wideband radio observations of the Galactic Center magnetar (SGR 1745–2900) with the Green Bank Telescope for ∼100 days immediately following its initial X-ray outburst in 2013 April. We made multiple simultaneous observations at 1.5, 2.0, and 8.9 GHz, allowing us to examine the magnetar's flux evolution, radio spectrum, and interstellar medium parameters (such as the dispersion measure (DM), the scattering timescale, and its index). During two epochs, we have simultaneous observations from the Chandra X-ray Observatory, which permitted the absolute alignment of the radio and X-ray profiles. As with the two other radio magnetars with published alignments, the radio profile lies within the broad peak of the X-ray profile, preceding the X-ray profile maximum by ∼0.2 rotations. We also find that the radio spectral index γ is significantly negative between ∼2 and 9 GHz; during the final ∼30 days of our observations γ∼−1.4, which is typical of canonical pulsars. The radio flux has not decreased during this outburst, whereas the long-term trends in the other radio magnetars show concomitant fading of the radio and X-ray fluxes. Finally, our wideband measurements of the DMs taken in adjacent frequency bands in tandem are stochastically inconsistent with one another. Based on recent theoretical predictions, we consider the possibility that the DM is frequency-dependent. Despite having several properties in common with the other radio magnetars, such as L{sub X,qui}/L{sub rot}≲1, an increase in the radio flux during the X-ray flux decay has not been observed thus far in other systems.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

CALIBRATING UV STAR FORMATION RATES FOR DWARF GALAXIES FROM STARBIRDS

Integrating our knowledge of star formation (SF) traced by observations at different wavelengths is essential for correctly interpreting and comparing SF activity in a variety of systems and environments. This study compares extinction corrected integrated ultraviolet (UV) emission from resolved galaxies with color–magnitude diagram (CMD) based star formation rates (SFRs) derived from resolved stellar populations and CMD fitting techniques in 19 nearby starburst and post-starburst dwarf galaxies. The data sets are from the panchromatic Starburst Irregular Dwarf Survey and include deep legacy GALEX UV imaging, Hubble Space Telescope optical imaging, and Spitzer MIPS imaging. For the majority of the sample, the integrated near-UV fluxes predicted from the CMD-based SFRs—using four different models—agree with the measured, extinction corrected, integrated near-UV fluxes from GALEX images, but the far-UV (FUV) predicted fluxes do not. Furthermore, we find a systematic deviation between the SFRs based on integrated FUV luminosities and existing scaling relations, and the SFRs based on the resolved stellar populations. This offset is not driven by different SF timescales, variations in SFRs, UV attenuation, nor stochastic effects. This first comparison between CMD-based SFRs and an integrated FUV emission SFR indicator suggests that the most likely cause of the discrepancy is the theoretical FUV–SFR calibration from stellar evolutionary libraries and/or stellar atmospheric models. We present an empirical calibration of the FUV-based SFR relation for dwarf galaxies, with uncertainties, which is ∼53% larger than previous relations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

A PILOT DEEP SURVEY FOR X-RAY EMISSION FROM fuvAGB STARS

We report the results of a pilot survey for X-ray emission from a newly discovered class of AGB stars with far-ultraviolet excesses (fuvAGB stars) using XMM-Newton and Chandra. We detected X-ray emission in three of six fuvAGB stars observed—the X-ray fluxes are found to vary in a stochastic or quasi-periodic manner on roughly hour-long timescales, and simultaneous UV observations using the Optical Monitor on XMM for these sources show similar variations in the UV flux. These data, together with previous studies, show that X-ray emission is found only in fuvAGB stars. From modeling the spectra, we find that the observed X-ray luminosities are ∼(0.002–0.2) L{sub ⊙} and the X-ray-emitting plasma temperatures are ∼(35–160) × 10{sup 6} K. The high X-ray temperatures argue against the emission arising in stellar coronae, or directly in an accretion shock, unless it occurs on a WD companion. However, none of the detected objects is a known WD-symbiotic star, suggesting that if WD companions are present, they are relatively cool (<20,000 K). In addition, the high X-ray luminosities specifically argue against emission originating in the coronae of main-sequence companions. We discuss several models for the X-ray emission and its variability and find that the most likely scenario for the origin of the X-ray (and FUV) emission involves accretion activity around a companion star, with confinement by strong magnetic fields associated with the companion and/or an accretion disk around it.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

STAR CLUSTER PROPERTIES IN TWO LEGUS GALAXIES COMPUTED WITH STOCHASTIC STELLAR POPULATION SYNTHESIS MODELS

We investigate a novel Bayesian analysis method, based on the Stochastically Lighting Up Galaxies (slug) code, to derive the masses, ages, and extinctions of star clusters from integrated light photometry. Unlike many analysis methods, slug correctly accounts for incomplete initial mass function (IMF) sampling, and returns full posterior probability distributions rather than simply probability maxima. We apply our technique to 621 visually confirmed clusters in two nearby galaxies, NGC 628 and NGC 7793, that are part of the Legacy Extragalactic UV Survey (LEGUS). LEGUS provides Hubble Space Telescope photometry in the NUV, U, B, V, and I bands. We analyze the sensitivity of the derived cluster properties to choices of prior probability distribution, evolutionary tracks, IMF, metallicity, treatment of nebular emission, and extinction curve. We find that slug's results for individual clusters are insensitive to most of these choices, but that the posterior probability distributions we derive are often quite broad, and sometimes multi-peaked and quite sensitive to the choice of priors. In contrast, the properties of the cluster population as a whole are relatively robust against all of these choices. We also compare our results from slug to those derived with a conventional non-stochastic fitting code, Yggdrasil. We show that slug's stochastic models are generally a better fit to the observations than the deterministic ones used by Yggdrasil. However, the overall properties of the cluster populations recovered by both codes are qualitatively similar.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

THE IMPACT OF THE IONOSPHERE ON GROUND-BASED DETECTION OF THE GLOBAL EPOCH OF REIONIZATION SIGNAL

The redshifted 21 cm line of neutral hydrogen (H i), potentially observable at low radio frequencies (∼50–200 MHz), is a promising probe of the physical conditions of the intergalactic medium during Cosmic Dawn and the Epoch of Reionization (EoR). The sky-averaged H i signal is expected to be extremely weak (∼100 mK) in comparison to the Galactic foreground emission (∼10{sup 4} K). Moreover, the sky-averaged spectra measured by ground-based instruments are affected by chromatic propagation effects (∼tens of kelvin) originating in the ionosphere. We analyze data collected with the upgraded Broadband Instrument for Global Hydrogen Reionization Signal system deployed at the Murchison Radio-astronomy Observatory to assess the significance of ionospheric effects on the detection of the global EoR signal. The ionospheric effects identified in these data are, particularly during nighttime, dominated by absorption and emission. We measure some properties of the ionosphere, such as the electron temperature (T{sub e} ≈ 470 K at nighttime), magnitude, and variability of optical depth (τ{sub 100} {sub MHz} ≈ 0.01 and δτ ≈ 0.005 at nighttime). According to the results of a statistical test applied on a large data sample, very long integrations (∼100 hr collected over approximately 2 months) lead to increased signal-to-noise ratio even in the presence of ionospheric variability. This is further supported by the structure of the power spectrum of the sky temperature fluctuations, which has flicker noise characteristics at frequencies ≳10{sup −5} Hz, but becomes flat below ≈10{sup −5} Hz. Hence, we conclude that the stochastic error introduced by the chromatic ionospheric effects tends to zero in an average. Therefore, the ionospheric effects and fluctuations are not fundamental impediments preventing ground-based instruments from integrating down to the precision required by global EoR experiments, provided that the ionospheric contribution is properly accounted for in the data analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

SENSITIVITY OF COSMIC-RAY PROTON SPECTRA TO THE LOW-WAVENUMBER BEHAVIOR OF THE 2D TURBULENCE POWER SPECTRUM

In this study, a novel ab initio cosmic ray (CR) modulation code that solves a set of stochastic transport equations equivalent to the Parker transport equation, and that uses output from a turbulence transport code as input for the diffusion tensor, is introduced. This code is benchmarked with a previous approach to ab initio modulation. The sensitivity of computed galactic CR proton spectra at Earth to assumptions made as to the low-wavenumber behavior of the two-dimensional (2D) turbulence power spectrum is investigated using perpendicular mean free path expressions derived from two different scattering theories. Constraints on the low-wavenumber behavior of the 2D power spectrum are inferred from the qualitative comparison of computed CR spectra with spacecraft observations at Earth. Another key difference from previous studies is that observed and inferred CR intensity spectra at 73 AU are used as boundary spectra instead of the usual local interstellar spectrum. Furthermore, the results presented here provide a tentative explanation as to the reason behind the unusually high galactic proton intensity spectra observed in 2009 during the recent unusual solar minimum.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

TIME-DEPENDENT STOCHASTIC ACCELERATION MODEL FOR FERMI BUBBLES

We study stochastic acceleration models for the Fermi bubbles. Turbulence is excited just behind the shock front via Kelvin–Helmholtz, Rayleigh–Taylor, or Richtmyer–Meshkov instabilities, and plasma particles are continuously accelerated by the interaction with the turbulence. The turbulence gradually decays as it goes away from the shock fronts. Adopting a phenomenological model for the stochastic acceleration, we explicitly solve the temporal evolution of the particle energy distribution in the turbulence. Our results show that the spatial distribution of high-energy particles is different from those for a steady solution. We also show that the contribution of electrons that escaped from the acceleration regions significantly softens the photon spectrum. The photon spectrum and surface brightness profile are reproduced by our models. If the escape efficiency is very high, the radio flux from the escaped low-energy electrons can be comparable to that of the WMAP haze. We also demonstrate hadronic models with the stochastic acceleration, but they are unlikely in the viewpoint of the energy budget.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

NONLINEAR EFFECTS IN PARTICLE TRANSPORT IN STOCHASTIC MAGNETIC FIELDS

Collisional particle transport in stochastic magnetic fields is studied using a semi-analytical method. The aim is to determine the influence of the nonlinear effects that occur in the magnetic field line random walk on particle transport. We show that particle transport coefficients can be strongly influenced by the magnetic line trapping. The conditions that correspond to these nonlinear regimes are determined. We also analyze the effects produced by the space variation of the large-scale magnetic field. We show that an average drift is generated by the gradient of the magnetic field, which strongly increases and reverses its orientation in the nonlinear regime.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

FLUCTUATIONS AND FLARES IN THE ULTRAVIOLET LINE EMISSION OF COOL STARS: IMPLICATIONS FOR EXOPLANET TRANSIT OBSERVATIONS

Variations in stellar flux can potentially overwhelm the photometric signal of a transiting planet. Such variability has not previously been well-characterized in the ultraviolet lines used to probe the inflated atmospheres surrounding hot Jupiters. Therefore, we surveyed 38 F-M stars for intensity variations in four narrow spectroscopic bands: two enclosing strong lines from species known to inhabit hot Jupiter atmospheres, C II λλ1334, 1335 and Si III λ1206; one enclosing Si IV λλ1393, 1402; and 36.5 Å of interspersed continuum. For each star/band combination, we generated 60 s cadence lightcurves from archival Hubble Space Telescope Cosmic Origins Spectrograph and Space Telescope Imaging Spectrograph time-tagged photon data. Within these lightcurves, we characterized flares and stochastic fluctuations as separate forms of variability. Flares: we used a cross-correlation approach to detect 116 flares. These events occur in the time-series an average of once per 2.5 hr, over 50% last 4 minutes or less, and most produce the strongest response in Si IV. If the flare occurred during a transit measurement integrated for 60 minutes, 90/116 would destroy the signal of an Earth, 27/116 Neptune, and 7/116 Jupiter, with the upward bias in flux ranging from 1% to 109% of quiescent levels. Fluctuations: photon noise and underlying stellar fluctuations produce scatter in the quiescent data. We model the stellar fluctuations as Gaussian white noise with standard deviation σ {sub x}. Maximum likelihood values of σ {sub x} range from 1% to 41% for 60 s measurements. These values suggest that many cool stars will only permit a transit detection to high confidence in ultraviolet resonance lines if the radius of the occulting disk is ≳1 R{sub J} . However, for some M dwarfs this limit can be as low as several R {sub ⊕}.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

STOCHASTIC MODEL OF THE SPIN DISTRIBUTION OF DARK MATTER HALOS

We employ a stochastic approach to probing the origin of the log-normal distributions of halo spin in N-body simulations. After analyzing spin evolution in halo merging trees, it was found that a spin change can be characterized by a stochastic random walk of angular momentum. Also, spin distributions generated by random walks are fairly consistent with those directly obtained from N-body simulations. We derived a stochastic differential equation from a widely used spin definition and measured the probability distributions of the derived angular momentum change from a massive set of halo merging trees. The roles of major merging and accretion are also statistically analyzed in evolving spin distributions. Several factors (local environment, halo mass, merging mass ratio, and redshift) are found to influence the angular momentum change. The spin distributions generated in the mean-field or void regions tend to shift slightly to a higher spin value compared with simulated spin distributions, which seems to be caused by the correlated random walks. We verified the assumption of randomness in the angular momentum change observed in the N-body simulation and detected several degrees of correlation between walks, which may provide a clue for the discrepancies between the simulated and generated spin distributions in the voids. However, the generated spin distributions in the group and cluster regions successfully match the simulated spin distribution. We also demonstrated that the log-normality of the spin distribution is a natural consequence of the stochastic differential equation of the halo spin, which is well described by the Geometric Brownian Motion model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

The probability equation for the cosmological comoving curvature perturbation

Fluctuations of the comoving curvature perturbation with wavelengths larger than the horizon length are governed by a Langevin equation whose stochastic noise arise from the quantum fluctuations that are assumed to become classical at horizon crossing. The infrared part of the curvature perturbation performs a random walk under the action of the stochastic noise and, at the same time, it suffers a classical force caused by its self-interaction. By a path-interal approach and, alternatively, by the standard procedure in random walk analysis of adiabatic elimination of fast variables, we derive the corresponding Kramers-Moyal equation which describes how the probability distribution of the comoving curvature perturbation at a given spatial point evolves in time and is a generalization of the Fokker-Planck equation. This approach offers an alternative way to study the late time behaviour of the correlators of the curvature perturbation from infrared effects.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

Local stochastic non-Gaussianity and N-body simulations

Large-scale clustering of highly biased tracers of large-scale structure has emerged as one of the best observational probes of primordial non-Gaussianity of the local type (i.e. f{sub NL}{sup local}). This type of non-Gaussianity can be generated in multifield models of inflation such as the curvaton model. Recently, Tseliakhovich, Hirata, and Slosar showed that the clustering statistics depend qualitatively on the ratio of inflaton to curvaton power ξ after reheating, a free parameter of the model. If ξ is significantly different from zero, so that the inflaton makes a non-negligible contribution to the primordial adiabatic curvature, then the peak-background split ansatz predicts that the halo bias will be stochastic on large scales. In this paper, we test this prediction in N-body simulations. We find that large-scale stochasticity is generated, in qualitative agreement with the prediction, but that the level of stochasticity is overpredicted by ≈ 30%. Other predictions, such as ξ independence of the halo bias, are confirmed by the simulations. Surprisingly, even in the Gaussian case we do not find that halo model predictions for stochasticity agree consistently with simulations, suggesting that semi-analytic modeling of stochasticity is generally more difficult than modeling halo bias.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

A new analytic solution for 2nd-order Fermi acceleration

A new analytic solution for 2nd-order Fermi acceleration is presented. In particular, we consider time-dependent rates for stochastic acceleration, diffusive and convective escape as well as adiabatic losses. The power law index q of the turbulence spectrum is unconstrained and can therefore account for Kolmogorov (q = 5/3) and Kraichnan (q = 3/2) turbulence, Bohm diffusion (q = 1) as well as the hard-sphere approximation (q = 2). This considerably improves beyond solutions known to date and will prove a useful tool for more realistic modelling of 2nd-order Fermi acceleration in a variety of astrophysical environments.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

Diffusive propagation of cosmic rays from supernova remnants in the Galaxy. I: spectrum and chemical composition

In this paper we investigate the effect of stochasticity in the spatial and temporal distribution of supernova remnants on the spectrum and chemical composition of cosmic rays observed at Earth. The calculations are carried out for different choices of the diffusion coefficient D(E) experienced by cosmic rays during propagation in the Galaxy. In particular, at high energies we assume that D(E)∝E{sup δ}, with δ = 1/3 and δ = 0.6 being the reference scenarios. The large scale distribution of supernova remnants in the Galaxy is modeled following the distribution of pulsars, with and without accounting for the spiral structure of the Galaxy. We find that the stochastic fluctuations induced by the spatial and temporal distribution of supernovae, together with the effect of spallation of nuclei, lead to mild but sensible violations of the simple, leaky-box-inspired rule that the spectrum observed at Earth is N(E)∝E{sup −α} with α = γ+δ, where γ is the slope of the cosmic ray injection spectrum at the sources. Spallation of nuclei, even with the small rates appropriate for He, may account for small differences in spectral slopes between different nuclei, possibly providing an explanation for the recent CREAM observations. For δ = 1/3 we find that the slope of the proton and helium spectra are ∼ 2.67 and ∼ 2.6 respectively (with fluctuations depending on the realization of source distribution) at energies around ∼ 1 TeV (to be compared with the measured values of 2.66±0.02 and 2.58±0.02). For δ = 0.6 the hardening of the He spectra is not observed. The stochastic effects discussed above cannot be found in ordinary propagation calculations, such as GALPROP, where these effects and the point like nature of the sources are not taken into account. We also comment on the effect of time dependence of the escape of cosmic rays from supernova remnants, and of a possible clustering of the sources in superbubbles. In a second paper we will discuss the implications of these different scenarios for the anisotropy of cosmic rays.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

Diffusive propagation of cosmic rays from supernova remnants in the Galaxy. II: anisotropy

In this paper we investigate the effects of stochasticity in the spatial and temporal distribution of supernova remnants on the anisotropy of cosmic rays observed at Earth. The calculations are carried out for different choices of the diffusion coefficient D(E) experienced by cosmic rays during propagation in the Galaxy. The propagation and spallation of nuclei (with charge 1 ≤ Z ≤ 26) are taken into account. At high energies (E > 1 TeV) we assume that D(E)∝(E/Z){sup δ}, with δ = 1/3 and δ = 0.6 being the reference scenarios. The large scale distribution of supernova remnants in the Galaxy is modeled following the distribution of pulsars with and without accounting for the spiral structure of the Galaxy. Our calculations allow us to determine the contribution to anisotropy resulting from both the large scale distribution of SNRs in the Galaxy and the random distribution of the nearest remnants. The naive expectation that the anisotropy amplitude scales as δ{sub A}∝D(E) is shown to be a wild oversimplification of reality which does not reflect in the predicted anisotropy for any realistic distribution of the sources. The fluctuations in the anisotropy pattern are dominated by nearby sources, so that predicting or explaining the observed anisotropy amplitude and phase becomes close to impossible. Nevertheless, the results of our calculations, when compared to the data, allow us to draw interesting conclusions in terms of the propagation scenario to be preferred both in terms of the energy dependence of the diffusion coefficient and of the size of the halo. We find that the very weak energy dependence of the anisotropy amplitude below 10{sup 5} GeV, as observed by numerous experiments, as well as the rise at higher energies, can best be explained if the diffusion coefficient is D(E)∝E{sup 1/3}. Faster diffusion, for instance with δ = 0.6, leads in general to an exceedingly large anisotropy amplitude. The spiral structure introduces interesting trends in the energy dependence of the anisotropy pattern, which qualitatively reflect the trend seen in the data. The inhomogeneous spatial distribution of the sources in the Galactic disc induces a large scale anisotropy which is not sensitive to the stochastic nature of nearby SNRs: we find that this additional contribution to δ{sub A} becomes more important for large values of the size of the halo, H. The two terms are comparable in size for H ∼ 2 kpc which corresponds to the scale height of the gradient of the spatial distribution of SNRs in the Galaxy. The dependence on energy of δ{sub A}(E) is close to monotonic when the large-scale, regular term dominates, and does not seem to reflect the observed anisotropy amplitude. Both contributions to the total anisotropy are illustrated and discussed with the help of semi-analytical results.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

Cosmological perturbations of self-accelerating universe in nonlinear massive gravity

We study cosmological perturbations of self-accelerating universe solutions in the recently proposed nonlinear theory of massive gravity, with general matter content. While the broken diffeomorphism invariance implies that there generically are 2 tensor, 2 vector and 2 scalar degrees of freedom in the gravity sector, we find that the scalar and vector degrees have vanishing kinetic terms and nonzero mass terms. Depending on their nonlinear behavior, this indicates either nondynamical nature of these degrees or strong couplings. Assuming the former, we integrate out the 2 vector and 2 scalar degrees of freedom. We then find that in the scalar and vector sectors, gauge-invariant variables constructed from metric and matter perturbations have exactly the same quadratic action as in general relativity. The difference from general relativity arises only in the tensor sector, where the graviton mass modifies the dispersion relation of gravitational waves, with a time-dependent effective mass. This may lead to modification of stochastic gravitational wave spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗