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

Strongly hyperbolic second order Einstein's evolution equations

BSSN-type evolution equations are discussed. The name refers to the Baumgarte, Shapiro, Shibata, and Nakamura version of the Einstein evolution equations, without introducing the conformal-traceless decomposition but keeping the three connection functions and including a densitized lapse. It is proved that a pseudo-differential first order reduction of these equations is strongly hyperbolic. In the same way, densitized Arnowitt-Deser-Misner evolution equations are found to be weakly hyperbolic. In both cases, the positive densitized lapse function and the spacelike shift vector are arbitrary given fields. This first order pseudodifferential reduction adds no extra equations to the system and so no extra constraints.Comment: LaTeX, 16 pages, uses revtex4. Referee corections and new appendix added. English grammar improved; typos correcte

Calculating Phases Between B => K* pi Amplitudes

A phase $\Delta \Phi$ between amplitudes for $B^0 \to K^{*0} \pi^0$ and $B^0 \to K^{*+} \pi^-$ plays a crucial role in a method for constraining Cabibbo-Kobayashi-Maskawa (CKM) parameters. We present a general argument for destructive interference between amplitudes for $B^0 \to K^{*+} \pi^-$ and $B^0 \to K^{*0} \pi^0$ forming together a smaller $I(K^* \pi) = 3/2$ amplitude. Applying flavor SU(3) and allowing for conservative theoretical uncertainties, we obtain lower limits on $|\Delta \Phi|$ and its charge-conjugate. Values of these two phases favored by the BaBar Collaboration are in good agreement with our bounds.Comment: 12 pages, 6 figures; slight revisions and clarification

Minimal Flavor Constraints for Technicolor

We analyze the constraints on the the vacuum polarization of the standard model gauge bosons from a minimal set of flavor observables valid for a general class of models of dynamical electroweak symmetry breaking. We will show that the constraints have a strong impact on the self-coupling and masses of the lightest spin-one resonances. Our analysis is applicable to any four and higher dimensional extension of the standard model reducing to models of dynamical electroweak symmetry breaking.Comment: 26 pages, we have added appendix C and some references and corrected some typo

Restriction on the energy and luminosity of e+e- storage rings due to beamstrahlung

The role of beamstrahlung in high-energy e+e- storage-ring colliders (SRCs) is examined. Particle loss due to the emission of single energetic beamstrahlung photons is shown to impose a fundamental limit on SRC luminosities at energies 2E_0 >~ 140 GeV for head-on collisions and 2E_0 >~ 40 GeV for crab-waist collisions. With beamstrahlung taken into account, we explore the viability of SRCs in the E_0=240-500 GeV range, which is of interest in the precision study of the Higgs boson. At 2E_0=240 GeV, SRCs are found to be competitive with linear colliders; however, at 2E_0=400-500 GeV, the attainable SRC luminosity would be a factor 15-25 smaller than desired.Comment: Latex, 5 pages. v2 differs only by minor changes is abstract and introduction, one reference is added. v3 corresponds to the paper published in PR

Improved Determination of the CKM Angle alpha from B to pi pi decays

Motivated by a recent paper that compares the results of the analysis of the CKM angle alpha in the frequentist and in the Bayesian approaches, we have reconsidered the information on the hadronic amplitudes, which helps constraining the value of alpha in the Standard Model. We find that the Bayesian method gives consistent results irrespective of the parametrisation of the hadronic amplitudes and that the results of the frequentist and Bayesian approaches are equivalent when comparing meaningful probability ranges or confidence levels. We also find that from B to pi pi decays alone the 95% probability region for alpha is the interval [80^o,170^o], well consistent with recent analyses of the unitarity triangle where, by using all the available experimental and theoretical information, one gets alpha = (93 +- 4)^o. Last but not least, by using simple arguments on the hadronic matrix elements, we show that the unphysical region alpha ~ 0, present in several experimental analyses, can be eliminated.Comment: 16 pages, 7 figure

Electric dipole moments from spontaneous CP violation in SU(3)-flavoured SUSY

The SUSY flavour problem is deeply related to the origin of flavour and hence to the origin of the SM Yukawa couplings themselves. Since all CP-violation in the SM is restricted to the flavour sector, it is possible that the SUSY CP problem is related to the origin of flavour as well. In this work, we present three variations of an SU(3) flavour model with spontaneous CP violation. Such models explain the hierarchy in the fermion masses and mixings, and predict the structure of the flavoured soft SUSY breaking terms. In such a situation, both SUSY flavour and CP problems do not exist. We use electric dipole moments and lepton flavour violation processes to distinguish between these models, and place constraints on the SUSY parameter space.Comment: 8 pages, 2 Figures, to appear in the DISCRETE'08 proceeding

First order hyperbolic formalism for Numerical Relativity

The causal structure of Einstein's evolution equations is considered. We show that in general they can be written as a first order system of balance laws for any choice of slicing or shift. We also show how certain terms in the evolution equations, that can lead to numerical inaccuracies, can be eliminated by using the Hamiltonian constraint. Furthermore, we show that the entire system is hyperbolic when the time coordinate is chosen in an invariant algebraic way, and for any fixed choice of the shift. This is achieved by using the momentum constraints in such as way that no additional space or time derivatives of the equations need to be computed. The slicings that allow hyperbolicity in this formulation belong to a large class, including harmonic, maximal, and many others that have been commonly used in numerical relativity. We provide details of some of the advanced numerical methods that this formulation of the equations allows, and we also discuss certain advantages that a hyperbolic formulation provides when treating boundary conditions.Comment: To appear in Phys. Rev.

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