Gauge and constraint degrees of freedom: from analytical to numerical approximations in General Relativity

The harmonic formulation of Einstein's field equations is considered, where the gauge conditions are introduced as dynamical constraints. The difference between the fully constrained approach (used in analytical approximations) and the free evolution one (used in most numerical approximations) is pointed out. As a generalization, quasi-stationary gauge conditions are also discussed, including numerical experiments with the gauge-waves testbed. The complementary 3+1 approach is also considered, where constraints are related instead with energy and momentum first integrals and the gauge must be provided separately. The relationship between the two formalisms is discussed in a more general framework (Z4 formalism). Different strategies in black hole simulations follow when introducing singularity avoidance as a requirement. More flexible quasi-stationary gauge conditions are proposed in this context, which can be seen as generalizations of the current 'freezing shift' prescriptions.Comment: Talk given at the Spanish Relativity Meeting, Tenerife, September 200

A new dissipation term for finite-difference simulations in Relativity

We present a new numerical dissipation algorithm, which can be efficiently used in combination with centered finite-difference methods. We start from a formulation of centered finite-volume methods for Numerical Relativity, in which third-order space accuracy can be obtained by employing just piecewise-linear reconstruction. We obtain a simplified version of the algorithm, which can be viewed as a centered finite-difference method plus some 'adaptive dissipation'. The performance of this algorithm is confirmed by numerical results obtained from 3D black hole simulations.Comment: Talk presented at the Spanish Relativity Meeting (Tenerife 2007

General-relativistic resistive-magnetohydrodynamic simulations of binary neutron stars

We have studied the dynamics of an equal-mass magnetized neutron-star binary within a resistive magnetohydrodynamic (RMHD) approach in which the highly conducting stellar interior is matched to an electrovacuum exterior. Because our analysis is aimed at assessing the modifications introduced by resistive effects on the dynamics of the binary after the merger and through to collapse, we have carried out a close comparison with an equivalent simulation performed within the traditional ideal magnetohydrodynamic approximation. We have found that there are many similarities between the two evolutions but also one important difference: the survival time of the hyper massive neutron star increases in a RMHD simulation. This difference is due to a less efficient magnetic-braking mechanism in the resistive regime, in which matter can move across magnetic-field lines, thus reducing the outward transport of angular momentum. Both the RMHD and the ideal magnetohydrodynamic simulations carried here have been performed at higher resolutions and with a different grid structure than those in previous work of ours [L. Rezzolla, B. Giacomazzo, L. Baiotti, J. Granot, C. Kouveliotou, and M. A. Aloy, Astrophys. J. Letters 732, L6 (2011)], but confirm the formation of a low-density funnel with an ordered magnetic field produced by the black hole--torus system. In both regimes the magnetic field is predominantly toroidal in the highly conducting torus and predominantly poloidal in the nearly evacuated funnel. Reconnection processes or neutrino annihilation occurring in the funnel, none of which we model, could potentially increase the internal energy in the funnel and launch a relativistic outflow, which, however, is not produced in these simulations.Comment: 26 pages, 17 figures; animations available at http://www.southampton.ac.uk/~kd10g13/movies/index.shtm

From Tensor Equations to Numerical Code -- Computer Algebra Tools for Numerical Relativity

In this paper we present our recent work in developing a computer-algebra tool for systems of partial differential equations (PDEs), termed "Kranc". Our work is motivated by the problem of finding solutions of the Einstein equations through numerical simulations. Kranc consists of Mathematica based computer-algebra packages, that facilitate the task of dealing with symbolic tensorial calculations and realize the conversion of systems of partial differential evolution equations into parallelized C or Fortran code.Comment: LaTeX llncs style, 9 pages, 1 figure, to appaer in the proceedings of "SYNASC 2004 - 6th International Symposium on Symbolic and Numeric Algorithms for Scientific Computing", Timisoara, Romania, September 26-30 200

Dynamical damping terms for symmetry-seeking shift conditions

Suitable gauge conditions are fundamental for stable and accurate numerical-relativity simulations of inspiralling compact binaries. A number of well-studied conditions have been developed over the last decade for both the lapse and the shift and these have been successfully used both in vacuum and non-vacuum spacetimes when simulating binaries with comparable masses. At the same time, recent evidence has emerged that the standard "Gamma-driver" shift condition requires a careful and non-trivial tuning of its parameters to ensure long-term stable evolutions of unequal-mass binaries. We present a novel gauge condition in which the damping constant is promoted to be a dynamical variable and the solution of an evolution equation. We show that this choice removes the need for special tuning and provides a shift damping term which is free of instabilities in our simulations and dynamically adapts to the individual positions and masses of the binary black-hole system. Our gauge condition also reduces the variations in the coordinate size of the apparent horizon of the larger black hole and could therefore be useful when simulating binaries with very small mass ratios.Comment: 11 pages, 8 figure

Quantification of tumour heterogenity in MRI

Cancer is the leading cause of death that touches us all, either directly or indirectly. It is estimated that the number of newly diagnosed cases in the Netherlands will increase to 123,000 by the year 2020. General Dutch statistics are similar to those in the UK, i.e. over the last ten years, the age-standardised incidence rate1 has stabilised at around 355 females and 415 males per 100,000. Figure 1 shows the cancer incidence per gender. In the UK, the rise in lifetime risk of cancer is more than one in three and depends on many factors, including age, lifestyle and genetic makeup

General-relativistic resistive magnetohydrodynamics in three dimensions: Formulation and tests

We present a new numerical implementation of the general-relativistic resistive magnetohydrodynamics (MHD) equations within the Whisky code. The numerical method adopted exploits the properties of implicit-explicit Runge-Kutta numerical schemes to treat the stiff terms that appear in the equations for large electrical conductivities. Using tests in one, two, and three dimensions, we show that our implementation is robust and recovers the ideal-MHD limit in regimes of very high conductivity. Moreover, the results illustrate that the code is capable of describing scenarios in a very wide range of conductivities. In addition to tests in flat spacetime, we report simulations of magnetized nonrotating relativistic stars, both in the Cowling approximation and in dynamical spacetimes. Finally, because of its astrophysical relevance and because it provides a severe testbed for general-relativistic codes with dynamical electromagnetic fields, we study the collapse of a nonrotating star to a black hole. We show that also in this case our results on the quasinormal mode frequencies of the excited electromagnetic fields in the Schwarzschild background agree with the perturbative studies within 0.7% and 5.6% for the real and the imaginary part of the l=1 mode eigenfrequency, respectively. Finally we provide an estimate of the electromagnetic efficiency of this process.Comment: 22 pages, 19 figure

On the black hole from merging binary neutron stars: how fast can it spin?

The merger of two neutron stars will in general lead to the formation of a torus surrounding a black hole whose rotational energy can be tapped to potentially power a short gamma-ray burst. We have studied the merger of equal-mass binaries with spins aligned with the orbital angular momentum to determine the maximum spin the black hole can reach. Our initial data consists of irrotational binaries to which we add various amounts of rotation to increase the total angular momentum. Although the initial data violates the constraint equations, the use of the constraint-damping CCZ4 formulation yields evolutions with violations smaller than those with irrotational initial data and standard formulations. Interestingly, we find that a limit of $J/M^2 \simeq 0.89$ exists for the dimensionless spin and that any additional angular momentum given to the binary ends up in the torus rather than in the black hole, thus providing another nontrivial example supporting the cosmic censorship hypothesis.Comment: 4 pages, 2 figures Version to appear in PRD Rapid Communication

Constraint damping of the conformal and covariant formulation of the Z4 system in simulations of binary neutron stars

Following previous work in vacuum spacetimes, we investigate the constraint-damping properties in the presence of matter of the recently developed traceless, conformal and covariant Z4 (CCZ4) formulation of the Einstein equations. First, we evolve an isolated neutron star with an ideal gas equation of state and subject to a constraint-violating perturbation. We compare the evolution of the constraints using the CCZ4 and Baumgarte-Shibata-Shapiro-Nakamura-Oohara-Kojima (BSSNOK) systems. Second, we study the collapse of an unstable spherical star to a black hole. Finally, we evolve binary neutron star systems over several orbits until the merger, the formation of a black hole, and up to the ringdown. We show that the CCZ4 formulation is stable in the presence of matter and that the constraint violations are one or more orders of magnitude smaller than for the BSSNOK formulation. Furthermore, by comparing the CCZ4 and the BSSNOK formulations also for neutron star binaries with large initial constraint violations, we investigate their influence on the errors on physical quantities. We also give a new, simple and robust prescription for the damping parameter that removes the instabilities found when using the fully covariant version of CCZ4 in the evolution of black holes. Overall, we find that at essentially the same computational costs the CCZ4 formulation provides solutions that are stable and with a considerably smaller violation of the Hamiltonian constraint than the BSSNOK formulation. We also find that the performance of the CCZ4 formulation is very similar to another conformal and traceless, but noncovariant formulation of the Z4 system, i.e. the Z4c formulation.Comment: 15 pages, 11 figures; accepted for publication in Phys. Rev.