3D simulations of Einstein's equations: symmetric hyperbolicity, live gauges and dynamic control of the constraints

We present three-dimensional simulations of Einstein equations implementing a symmetric hyperbolic system of equations with dynamical lapse. The numerical implementation makes use of techniques that guarantee linear numerical stability for the associated initial-boundary value problem. The code is first tested with a gauge wave solution, where rather larger amplitudes and for significantly longer times are obtained with respect to other state of the art implementations. Additionally, by minimizing a suitably defined energy for the constraints in terms of free constraint-functions in the formulation one can dynamically single out preferred values of these functions for the problem at hand. We apply the technique to fully three-dimensional simulations of a stationary black hole spacetime with excision of the singularity, considerably extending the lifetime of the simulations.Comment: 21 pages. To appear in PR

Improved Laboratory Transition Probabilities for Neutral Chromium and Re-determination of the Chromium Abundance for the Sun and Three Stars

Branching fraction measurements from Fourier transform spectra in conjunction with published radiative lifetimes are used to determine transition probabilities for 263 lines of neutral chromium. These laboratory values are employed to derive a new photospheric abundance for the Sun: log $\epsilon$(Cr I)$_{\odot}$ = 5.64$\pm$0.01 ($\sigma = 0.07$). These Cr I solar abundances do not exhibit any trends with line strength nor with excitation energy and there were no obvious indications of departures from LTE. In addition, oscillator strengths for singly-ionized chromium recently reported by the FERRUM Project are used to determine: log $\epsilon$(Cr II)$_{\odot}$ = 5.77$\pm$0.03 ($\sigma = 0.13$). Transition probability data are also applied to the spectra of three stars: HD 75732 (metal-rich dwarf), HD 140283 (metal-poor subgiant), and CS 22892-052 (metal-poor giant). In all of the selected stars, Cr I is found to be underabundant with respect to Cr II. The possible causes for this abundance discrepancy and apparent ionization imbalance are discussed.Comment: 44 pages, 6 figure

Boundary conditions for hyperbolic formulations of the Einstein equations

In regards to the initial-boundary value problem of the Einstein equations, we argue that the projection of the Einstein equations along the normal to the boundary yields necessary and appropriate boundary conditions for a wide class of equivalent formulations. We explicitly show that this is so for the Einstein-Christoffel formulation of the Einstein equations in the case of spherical symmetry.Comment: 15 pages; text added and typesetting errors corrected; to appear in Classical and Quantum Gravit

High-powered Gravitational News

We describe the computation of the Bondi news for gravitational radiation. We have implemented a computer code for this problem. We discuss the theory behind it as well as the results of validation tests. Our approach uses the compactified null cone formalism, with the computational domain extending to future null infinity and with a worldtube as inner boundary. We calculate the appropriate full Einstein equations in computational eth form in (a) the interior of the computational domain and (b) on the inner boundary. At future null infinity, we transform the computed data into standard Bondi coordinates and so are able to express the news in terms of its standard $N_{+}$ and $N_{\times}$ polarization components. The resulting code is stable and second-order convergent. It runs successfully even in the highly nonlinear case, and has been tested with the news as high as 400, which represents a gravitational radiation power of about $10^{13}M_{\odot}/sec$.Comment: 24 pages, 4 figures. To appear in Phys. Rev.

Measurement of the analyzing power Ay0 for the reaction H(p⃗,d)π+ between 1000 and 1300 MeV

The analyzing power Ay0 of the reaction H(p⃗,d)π+ has been measured at a fixed value of the Mandelstam variable ud=-0.17GeV2 for nine proton energies between 1000 and 1300 MeV. The experiment was performed at SATURNE with the SPES1 spectrometer. The data exhibit structure around √s≃2.37GeV. The origin of this structure could be related to a resonancelike behavior of the 1S0P or 1G4F partial amplitudes

Virtual Compton Scattering and the Generalized Polarizabilities of the Proton at Q^2=0.92 and 1.76 GeV^2

Virtual Compton Scattering (VCS) on the proton has been studied at Jefferson Lab using the exclusive photon electroproduction reaction (e p --> e p gamma). This paper gives a detailed account of the analysis which has led to the determination of the structure functions P_LL-P_TT/epsilon and P_LT, and the electric and magnetic generalized polarizabilities (GPs) alpha_E(Q^2) and beta_M(Q^2) at values of the four-momentum transfer squared Q^2= 0.92 and 1.76 GeV^2. These data, together with the results of VCS experiments at lower momenta, help building a coherent picture of the electric and magnetic GPs of the proton over the full measured Q^2-range, and point to their non-trivial behavior.Comment: version 2: modified according to PRC Editor's and Referee's recommendations. Archival paper for the E93-050 experiment at JLab Hall A. 28 pages, 23 figures, 5 cross-section tables. To be submitted to Phys.Rev.

Numerical Relativity: A review

Computer simulations are enabling researchers to investigate systems which are extremely difficult to handle analytically. In the particular case of General Relativity, numerical models have proved extremely valuable for investigations of strong field scenarios and been crucial to reveal unexpected phenomena. Considerable efforts are being spent to simulate astrophysically relevant simulations, understand different aspects of the theory and even provide insights in the search for a quantum theory of gravity. In the present article I review the present status of the field of Numerical Relativity, describe the techniques most commonly used and discuss open problems and (some) future prospects.Comment: 2 References added; 1 corrected. 67 pages. To appear in Classical and Quantum Gravity. (uses iopart.cls

Backward electroproduction of pi0 mesons on protons in the region of nucleon resonances at four momentum transfer squared Q**2 = 1.0 GeV**2

Exclusive electroproduction of pi0 mesons on protons in the backward hemisphere has been studied at Q**2 = 1.0 GeV**2 by detecting protons in the forward direction in coincidence with scattered electrons from the 4 GeV electron beam in Jefferson Lab's Hall A. The data span the range of the total (gamma* p) center-of-mass energy W from the pion production threshold to W = 2.0 GeV. The differential cross sections sigma_T+epsilon*sigma_L, sigma_TL, and sigma_TT were separated from the azimuthal distribution and are presented together with the MAID and SAID parametrizations.Comment: 17 pages, 11 figures, information can be found at http://hallaweb.jlab.org/experiment/E93-050/vcs.html updated content about SAID analysis updated MAID results following new reference nucl-th/0310041 updated figure