Hamiltonian Relaxation

Due to the complexity of the required numerical codes, many of the new formulations for the evolution of the gravitational fields in numerical relativity are not tested on binary evolutions. We introduce in this paper a new testing ground for numerical methods based on the simulation of binary neutron stars. This numerical setup is used to develop a new technique, the Hamiltonian relaxation (HR), that is benchmarked against the currently most stable simulations based on the BSSN method. We show that, while the length of the HR run is somewhat shorter than the equivalent BSSN simulation, the HR technique improves the overall quality of the simulation, not only regarding the satisfaction of the Hamiltonian constraint, but also the behavior of the total angular momentum of the binary. The latest quantity agrees well with post-Newtonian estimations for point-mass binaries in circular orbits.Comment: More detailed description of the numerical implementation added and some typos corrected. Version accepted for publication in Class. and Quantum Gravit

The constraints as evolution equations for numerical relativity

The Einstein equations have proven surprisingly difficult to solve numerically. A standard diagnostic of the problems which plague the field is the failure of computational schemes to satisfy the constraints, which are known to be mathematically conserved by the evolution equations. We describe a new approach to rewriting the constraints as first-order evolution equations, thereby guaranteeing that they are satisfied to a chosen accuracy by any discretization scheme. This introduces a set of four subsidiary constraints which are far simpler than the standard constraint equations, and which should be more easily conserved in computational applications. We explore the manner in which the momentum constraints are already incorporated in several existing formulations of the Einstein equations, and demonstrate the ease with which our new constraint-conserving approach can be incorporated into these schemes.Comment: 10 pages, updated to match published versio

The Milky Way Project: A statistical study of massive star formation associated with infrared bubbles

The Milky Way Project citizen science initiative recently increased the number of known infrared bubbles in the inner Galactic plane by an order of magnitude compared to previous studies. We present a detailed statistical analysis of this dataset with the Red MSX Source catalog of massive young stellar sources to investigate the association of these bubbles with massive star formation. We particularly address the question of massive triggered star formation near infrared bubbles. We find a strong positional correlation of massive young stellar objects (MYSOs) and H II regions with Milky Way Project bubbles at separations of < 2 bubble radii. As bubble sizes increase, a statistically significant overdensity of massive young sources emerges in the region of the bubble rims, possibly indicating the occurrence of triggered star formation. Based on numbers of bubble-associated RMS sources we find that 67+/-3% of MYSOs and (ultra)compact H II regions appear associated with a bubble. We estimate that approximately 22+/-2% of massive young stars may have formed as a result of feedback from expanding H II regions. Using MYSO-bubble correlations, we serendipitously recovered the location of the recently discovered massive cluster Mercer 81, suggesting the potential of such analyses for discovery of heavily extincted distant clusters.Comment: 16 pages, 17 figures. Accepted for publication in ApJ, comments welcome. Milky Way Project public data release available at http://www.milkywayproject.org/dat

Renormalization of the charged scalar field in curved space

The DeWitt-Schwinger proper time point-splitting procedure is applied to a massive complex scalar field with arbitrary curvature coupling interacting with a classical electromagnetic field in a general curved spacetime. The scalar field current is found to have a linear divergence. The presence of the external background gauge field is found to modify the stress-energy tensor results of Christensen for the neutral scalar field by adding terms of the form $(eF)^2$ to the logarithmic counterterms. These results are shown to be expected from an analysis of the degree of divergence of scalar quantum electrodynamics.Comment: 24 pages REVTe

Tips for implementing multigrid methods on domains containing holes

As part of our development of a computer code to perform 3D `constrained evolution' of Einstein's equations in 3+1 form, we discuss issues regarding the efficient solution of elliptic equations on domains containing holes (i.e., excised regions), via the multigrid method. We consider as a test case the Poisson equation with a nonlinear term added, as a means of illustrating the principles involved, and move to a "real world" 3-dimensional problem which is the solution of the conformally flat Hamiltonian constraint with Dirichlet and Robin boundary conditions. Using our vertex-centered multigrid code, we demonstrate globally second-order-accurate solutions of elliptic equations over domains containing holes, in two and three spatial dimensions. Keys to the success of this method are the choice of the restriction operator near the holes and definition of the location of the inner boundary. In some cases (e.g. two holes in two dimensions), more and more smoothing may be required as the mesh spacing decreases to zero; however for the resolutions currently of interest to many numerical relativists, it is feasible to maintain second order convergence by concentrating smoothing (spatially) where it is needed most. This paper, and our publicly available source code, are intended to serve as semi-pedagogical guides for those who may wish to implement similar schemes.Comment: 18 pages, 11 figures, LaTeX. Added clarifications and references re. scope of paper, mathematical foundations, relevance of work. Accepted for publication in Classical & Quantum Gravit

The discrete energy method in numerical relativity: Towards long-term stability

The energy method can be used to identify well-posed initial boundary value problems for quasi-linear, symmetric hyperbolic partial differential equations with maximally dissipative boundary conditions. A similar analysis of the discrete system can be used to construct stable finite difference equations for these problems at the linear level. In this paper we apply these techniques to some test problems commonly used in numerical relativity and observe that while we obtain convergent schemes, fast growing modes, or ``artificial instabilities,'' contaminate the solution. We find that these growing modes can partially arise from the lack of a Leibnitz rule for discrete derivatives and discuss ways to limit this spurious growth.Comment: 18 pages, 22 figure

Boosted three-dimensional black-hole evolutions with singularity excision

Binary black hole interactions provide potentially the strongest source of gravitational radiation for detectors currently under development. We present some results from the Binary Black Hole Grand Challenge Alliance three- dimensional Cauchy evolution module. These constitute essential steps towards modeling such interactions and predicting gravitational radiation waveforms. We report on single black hole evolutions and the first successful demonstration of a black hole moving freely through a three-dimensional computational grid via a Cauchy evolution: a hole moving ~6M at 0.1c during a total evolution of duration ~60M

Extending the lifetime of 3D black hole computations with a new hyperbolic system of evolution equations

We present a new many-parameter family of hyperbolic representations of Einstein's equations, which we obtain by a straightforward generalization of previously known systems. We solve the resulting evolution equations numerically for a Schwarzschild black hole in three spatial dimensions, and find that the stability of the simulation is strongly dependent on the form of the equations (i.e. the choice of parameters of the hyperbolic system), independent of the numerics. For an appropriate range of parameters we can evolve a single 3D black hole to $t \simeq 600 M$ -- $1300 M$, and are apparently limited by constraint-violating solutions of the evolution equations. We expect that our method should result in comparable times for evolutions of a binary black hole system.Comment: 11 pages, 2 figures, submitted to PR

Stable characteristic evolution of generic 3-dimensional single-black-hole spacetimes

We report new results which establish that the accurate 3-dimensional numerical simulation of generic single-black-hole spacetimes has been achieved by characteristic evolution with unlimited long term stability. Our results cover a selection of distorted, moving and spinning single black holes, with evolution times up to 60,000M.Comment: 4 pages, 3 figure

Gravitational wave extraction and outer boundary conditions by perturbative matching

We present a method for extracting gravitational radiation from a three-dimensional numerical relativity simulation and, using the extracted data, to provide outer boundary conditions. The method treats dynamical gravitational variables as nonspherical perturbations of Schwarzschild geometry. We discuss a code which implements this method and present results of tests which have been performed with a three dimensional numerical relativity code