Evolution of Magnetic Fields in Freely Decaying Magnetohydrodynamic Turbulence

We study the evolution of magnetic fields in freely decaying magnetohydrodynamic turbulence. By quasi-linearizing the Navier-Stokes equation, we solve analytically the induction equation in quasi-normal approximation. We find that, if the magnetic field is not helical, the magnetic energy and correlation length evolve in time respectively as E_B \propto t^{-2(1+p)/(3+p)} and \xi_B \propto t^{2/(3+p)}, where p is the index of initial power-law spectrum. In the helical case, the magnetic helicity is an almost conserved quantity and forces the magnetic energy and correlation length to scale as E_B \propto (log t)^{1/3} t^{-2/3} and \xi_B \propto (log t)^{-1/3} t^{2/3}.Comment: 4 pages, 2 figures; accepted for publication in PR

Making use of geometrical invariants in black hole collisions

We consider curvature invariants in the context of black hole collision simulations. In particular, we propose a simple and elegant combination of the Weyl invariants I and J, the {\sl speciality index} ${\cal S}$. In the context of black hole perturbations $\cal S$ provides a measure of the size of the distortions from an ideal Kerr black hole spacetime. Explicit calculations in well-known examples of axisymmetric black hole collisions demonstrate that this quantity may serve as a useful tool for predicting in which cases perturbative dynamics provide an accurate estimate of the radiation waveform and energy. This makes ${\cal S}$ particularly suited to studying the transition from nonlinear to linear dynamics and for invariant interpretation of numerical results.Comment: 4 pages, 3 eps figures, Revte

The Lazarus project: A pragmatic approach to binary black hole evolutions

We present a detailed description of techniques developed to combine 3D numerical simulations and, subsequently, a single black hole close-limit approximation. This method has made it possible to compute the first complete waveforms covering the post-orbital dynamics of a binary black hole system with the numerical simulation covering the essential non-linear interaction before the close limit becomes applicable for the late time dynamics. To determine when close-limit perturbation theory is applicable we apply a combination of invariant a priori estimates and a posteriori consistency checks of the robustness of our results against exchange of linear and non-linear treatments near the interface. Once the numerically modeled binary system reaches a regime that can be treated as perturbations of the Kerr spacetime, we must approximately relate the numerical coordinates to the perturbative background coordinates. We also perform a rotation of a numerically defined tetrad to asymptotically reproduce the tetrad required in the perturbative treatment. We can then produce numerical Cauchy data for the close-limit evolution in the form of the Weyl scalar $\psi_4$ and its time derivative $\partial_t\psi_4$ with both objects being first order coordinate and tetrad invariant. The Teukolsky equation in Boyer-Lindquist coordinates is adopted to further continue the evolution. To illustrate the application of these techniques we evolve a single Kerr hole and compute the spurious radiation as a measure of the error of the whole procedure. We also briefly discuss the extension of the project to make use of improved full numerical evolutions and outline the approach to a full understanding of astrophysical black hole binary systems which we can now pursue.Comment: New typos found in the version appeared in PRD. (Mostly found and collected by Bernard Kelly

Scattering of Dirac and Majorana Fermions off Domain Walls

We investigate the interaction of fermions having both Dirac and left-handed and right-handed Majorana mass terms with vacuum domain walls. By solving the equations of motion in thin-wall approximation, we calculate the reflection and transmission coefficients for the scattering of fermions off walls.Comment: 6 pages, 1 figure, some typos corrected, one reference added, major revisions, title changed, version accepted for publication in Phys. Rev.

Probing the QCD vacuum with an abelian chromomagnetic field: A study within an effective model

We study the response of the QCD vacuum to an external abelian chromomagnetic field in the framework of a non local Nambu-Jona Lasinio model with the Polyakov loop. We use the Lattice results on the deconfinement temperature of the pure gauge theory to compute the same quantity in the presence of dynamical quarks. We find a linear relationship between the deconfinement temperature with quarks and the squared root of the applied field strength, $gH$, in qualitative (and to some extent also quantitative) agreement with existing Lattice calculations. On the other hand, we find a discrepancy on the approximate chiral symmetry restoration: while Lattice results suggest the deconfinement and the chiral restoration remain linked even at non-zero value of $gH$, our results are consistent with a scenario in which the two transitions are separated as $gH$ is increased.Comment: 14 pages, 7 figures, RevTeX4. Published version, with enlarged abstract and minor changes in the main tex

Testing the Isotropy of the Universe with Type Ia Supernovae

We analyze the magnitude-redshift data of type Ia supernovae included in the Union and Union2 compilations in the framework of an anisotropic Bianchi type I cosmological model and in the presence of a dark energy fluid with anisotropic equation of state. We find that the amount of deviation from isotropy of the equation of state of dark energy, the skewness \delta, and the present level of anisotropy of the large-scale geometry of the Universe, the actual shear \Sigma_0, are constrained in the ranges -0.16 < \delta < 0.12 and -0.012 < \Sigma_0 < 0.012 (1\sigma C.L.) by Union2 data. Supernova data are then compatible with a standard isotropic universe (\delta = \Sigma_0 = 0), but a large level of anisotropy, both in the geometry of the Universe and in the equation of state of dark energy, is allowed.Comment: 12 pages, 7 figures, 2 tables. Union2 analysis added. New references added. To appear in Phys. Rev.

A perturbative solution for gravitational waves in quadratic gravity

We find a gravitational wave solution to the linearized version of quadratic gravity by adding successive perturbations to the Einstein's linearized field equations. We show that only the Ricci squared quadratic invariant contributes to give a different solution of those found in Einstein's general relativity. The perturbative solution is written as a power series in the $\beta$ parameter, the coefficient of the Ricci squared term in the quadratic gravitational action. We also show that, for monochromatic waves of a given angular frequency $\omega$, the perturbative solution can be summed out to give an exact solution to linearized version of quadratic gravity, for $0<\omega<c/\mid\beta\mid^{1/2}$. This result may lead to implications to the predictions for gravitational wave backgrounds of cosmological origin.Comment: 9 pages, to appear in CQ

Constraints on the anisotropy of dark energy

If the equation of state of dark energy is anisotropic there will be additional quadrupole anisotropy in the cosmic microwave background induced by the time dependent anisotropic stress quantified in terms of $\Delta w$. Assuming that the entire amplitude of the observed quadrupole is due to this anisotropy, we conservatively impose a limit of $|\Delta w| < 2.1\times 10^{-4}$ for any value of $w\ge -1$ assuming that $\Omega_{\rm m}<0.5$. This is considerably tighter than that which comes from SNe. Stronger limits, upto a factor of 10, are possible for specific values of $\Omega_{\rm m}$ and $w$. Since we assume this component is uncorrelated with the stochastic component from inflation, we find that both the expectation value and the sample variance are increased. There no improvement in the likelihood of an anomalously low quadrupole as suggested by previous work on an elliptical universe

Gravitational-Wave Recoil from the Ringdown Phase of Coalescing Black Hole Binaries

The gravitational recoil or "kick" of a black hole formed from the merger of two orbiting black holes, and caused by the anisotropic emission of gravitational radiation, is an astrophysically important phenomenon. We combine (i) an earlier calculation, using post-Newtonian theory, of the kick velocity accumulated up to the merger of two non-spinning black holes, (ii) a "close-limit approximation" calculation of the radiation emitted during the ringdown phase, and based on a solution of the Regge-Wheeler and Zerilli equations using initial data accurate to second post-Newtonian order. We prove that ringdown radiation produces a significant "anti-kick". Adding the contributions due to inspiral, merger and ringdown phases, our results for the net kick velocity agree with those from numerical relativity to 10-15 percent over a wide range of mass ratios, with a maximum velocity of 180 km/s at a mass ratio of 0.38.Comment: 9 pages, 5 figures; to appear in Class. Quant. Gra