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

New Formalism for Numerical Relativity

We present a new formulation of the Einstein equations that casts them in an explicitly first order, flux-conservative, hyperbolic form. We show that this now can be done for a wide class of time slicing conditions, including maximal slicing, making it potentially very useful for numerical relativity. This development permits the application to the Einstein equations of advanced numerical methods developed to solve the fluid dynamic equations, {\em without} overly restricting the time slicing, for the first time. The full set of characteristic fields and speeds is explicitly given.Comment: uucompresed PS file. 4 pages including 1 figure. Revised version adds a figure showing a comparison between the standard ADM approach and the new formulation. Also available at http://jean-luc.ncsa.uiuc.edu/Papers/ Appeared in Physical Review Letters 75, 600 (1995

CP Asymmetries in B to f_0 K_S Decays

We consider the branching ratio and the CP asymmetries in B to f_0(980)K_S decay to the end of determining the deviation of the time-dependent CP asymmetry from sin(2 beta) arising from Standard Model physics. We obtain Delta S_{f_0 K_S} within the context of the QCD factorization framework for the B to f_0(980)K_S decay amplitudes assuming the f_0(980) is a q\bar{q} state and employing a random scan over the theoretical parameter space to assess the possible range in Delta S_{f_0 K_S}. Imposing the value of the experimental branching ratio within 1 sigma and 3 sigma, respectively, of its central value as a constraint, we find the range of Delta S_{f_0 K_S} to be [0.018, 0.033] for a scan in which the parameters are allowed to vary within 1 sigma of their central values and the range [-0.019, 0.064] for a scan in which the parameters vary within 3 sigma of their central values.Comment: 27 pages, 10 figures, references adde

Observation of the semileptonic decays B ---> D* tau- anti-nu(tau) and evidence for B ---> D tau- anti-nu(tau

We present measurements of the semileptonic decays B- --> D0 tau- nubar, B- --> D*0 tau- nubar, B0bar --> D+ tau- nubar, and B0bar --> D*+ tau- nubar, which are potentially sensitive to non--Standard Model amplitudes. The data sample comprises 232x10^6 Upsilon(4S) --> BBbar decays collected with the BaBar detector. From a combined fit to B- and B0bar channels, we obtain the branching fractions B(B --> D tau- nubar) = (0.86 +/- 0.24 +/- 0.11 +/- 0.06)% and B(B --> D* tau- nubar) = (1.62 +/- 0.31 +/- 0.10 +/- 0.05)% (normalized for the B0bar), where the uncertainties are statistical, systematic, and normalization-mode-related

Search for CP violation in the decays D0 ---> K- K+ and D0 ---> pi- pi+

We measure CP-violating asymmetries of neutral charmed mesons in the modes D0 --> K- K+ and D0 --> pi- pi+ with the highest precision to date by using D0 --> K- pi+ decays to correct detector asymmetries. An analysis of 385.8 fb-1 of data collected with the BaBar detector yields values of aCP(KK) = (0.00 +/- 0.34 (stat.) +/- 0.13 (syst.))% and aCP(pipi) = (-0.24 +/- 0.52 (stat.) +/- 0.22 (syst.))%, which agree with Standard Model prediction

Improved estimates of rare K decay matrix-elements from Kl3 decays

The estimation of rare K decay matrix-elements from Kl3 experimental data is extended beyond LO in Chiral Perturbation Theory. Isospin-breaking effects at NLO (and partially NNLO) in the ChPT expansion, as well as QED radiative corrections are now accounted for. The analysis relies mainly on the cleanness of two specific ratios of form-factors, for which the theoretical control is excellent. As a result, the uncertainties on the K+ --> pi+ nu nubar and KL --> pi0 nu nubar matrix-elements are reduced by a factor of about 7 and 4, respectively, and similarly for the direct CP-violating contribution to KL --> pi0 l+ l-. They could be reduced even further with better experimental data for the Kl3 slopes and the K+l3 branching ratios. As a result, the non-parametric errors for B(K --> pi nu nubar) and for the direct CP-violating contributions to B(KL --> pi0 l+ l-) are now completely dominated by those on the short-distance physics.Comment: 16 pages, 1 figure. Numerical analysis updated to include the recent Kl3 data. To appear in Phys. Rev.

Gowdy waves as a test-bed for constraint-preserving boundary conditions

Gowdy waves, one of the standard 'apples with apples' tests, is proposed as a test-bed for constraint-preserving boundary conditions in the non-linear regime. As an illustration, energy-constraint preservation is separately tested in the Z4 framework. Both algebraic conditions, derived from energy estimates, and derivative conditions, deduced from the constraint-propagation system, are considered. The numerical errors at the boundary are of the same order than those at the interior points.Comment: 5 pages, 1 figure. Contribution to the Spanish Relativity Meeting 200

Three dimensional numerical relativity: the evolution of black holes

We report on a new 3D numerical code designed to solve the Einstein equations for general vacuum spacetimes. This code is based on the standard 3+1 approach using cartesian coordinates. We discuss the numerical techniques used in developing this code, and its performance on massively parallel and vector supercomputers. As a test case, we present evolutions for the first 3D black hole spacetimes. We identify a number of difficulties in evolving 3D black holes and suggest approaches to overcome them. We show how special treatment of the conformal factor can lead to more accurate evolution, and discuss techniques we developed to handle black hole spacetimes in the absence of symmetries. Many different slicing conditions are tested, including geodesic, maximal, and various algebraic conditions on the lapse. With current resolutions, limited by computer memory sizes, we show that with certain lapse conditions we can evolve the black hole to about $t=50M$, where $M$ is the black hole mass. Comparisons are made with results obtained by evolving spherical initial black hole data sets with a 1D spherically symmetric code. We also demonstrate that an ``apparent horizon locking shift'' can be used to prevent the development of large gradients in the metric functions that result from singularity avoiding time slicings. We compute the mass of the apparent horizon in these spacetimes, and find that in many cases it can be conserved to within about 5\% throughout the evolution with our techniques and current resolution.Comment: 35 pages, LaTeX with RevTeX 3.0 macros. 27 postscript figures taking 7 MB of space, uuencoded and gz-compressed into a 2MB uufile. Also available at http://jean-luc.ncsa.uiuc.edu/Papers/ and mpeg simulations at http://jean-luc.ncsa.uiuc.edu/Movies/ Submitted to Physical Review

Observation of B+ ---> a(1)+(1260) K0 and B0 ---> a(1)-(1260) K+

We present branching fraction measurements of the decays B^{+} -> a1(1260)^{+} K^{0} and B^{0} to a1(1260)^{-} K^{+} with a1(1260)^{+} -> pi^{-} pi^{+} pi^{+}. The data sample corresponds to 383 million B B-bar pairs produced in e^{+}e^{-} annihilation through the Y(4S) resonance. We measure the products of the branching fractions: B(B^{+}-> a1(1260)^{+} K^{0})B(a1(1260)^{+} -> pi^{-} pi^{+} pi^{+}) = (17.4 +/- 2.5 +/- 2.2) 10^{-6} B(B^{0}-> a1(1260)^{-} K^{+})B(a1(1260)^{-} -> pi^{+} pi^{-} pi^{-}) = (8.2 +/- 1.5 +/- 1.2) 10^{-6}. We also measure the charge asymmetries A_{ch}(B^{+} -> a1(1260)^{+} K^{0})= 0.12 +/- 0.11 +/- 0.02 and A_{ch}(B^{0} -> a1(1260)^{-} K^{+})= -0.16 +/- 0.12 +/- 0.01. The first uncertainty quoted is statistical and the second is systematic

Exploiting gauge and constraint freedom in hyperbolic formulations of Einstein's equations

We present new many-parameter families of strongly and symmetric hyperbolic formulations of Einstein's equations that include quite general algebraic and live gauge conditions for the lapse. The first system that we present has 30 variables and incorporates an algebraic relationship between the lapse and the determinant of the three metric that generalizes the densitized lapse prescription. The second system has 34 variables and uses a family of live gauges that generalizes the Bona-Masso slicing conditions. These systems have free parameters even after imposing hyperbolicity and are expected to be useful in 3D numerical evolutions. We discuss under what conditions there are no superluminal characteristic speeds