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.
We provide additional information on our recent study of the electromagnetic emission produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. As anticipa
ted in a recent letter, our results show that although a dual-jet structure is present, the associated luminosity is ~ 100 times smaller than the total one, which is predominantly quadrupolar. We here discuss the details of our implementation of the equations in which the force-free condition is not implemented at a discrete level, but rather obtained via a damping scheme which drives the solution to satisfy the correct condition. We show that this is important for a correct and accurate description of the current sheets that can develop in the course of the simulation. We also study in greater detail the three-dimensional charge distribution produced as a consequence of the inspiral and show that during the inspiral it possesses a complex but ordered structure which traces the motion of the two black holes. Finally, we provide quantitative estimates of the scaling of the electromagnetic emission with frequency, with the diffused part having a dependence that is the same as the gravitational-wave one and that scales as L ~ Omega^{10/3-8/3}, while the collimated one scales as L ~ Omega^{5/3-6/3}, thus with a steeper dependence than previously estimated. We discuss the impact of these results on the potential detectability of dual jets from supermassive black holes and the steps necessary for more accurate estimates.
We present a new formulation of the Einstein equations based on a conformal and traceless decomposition of the covariant form of the Z4 system. This formulation combines the advantages of a conformal decomposition, such as the one used in the BSSNOK
formulation (i.e. well-tested hyperbolic gauges, no need for excision, robustness to imperfect boundary conditions) with the advantages of a constraint-damped formulation, such as the generalized harmonic one (i.e. exponential decay of constraint violations when these are produced). We validate the new set of equations through standard tests and by evolving binary black hole systems. Overall, the new conformal formulation leads to a better behavior of the constraint equations and a rapid suppression of the violations when they occur. The changes necessary to implement the new conformal formulation in standard BSSNOK codes are very small as are the additional computational costs.
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 the
se 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.
