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

Simulating binary neutron stars: dynamics and gravitational waves

We model two mergers of orbiting binary neutron stars, the first forming a black hole and the second a differentially rotating neutron star. We extract gravitational waveforms in the wave zone. Comparisons to a post-Newtonian analysis allow us to compute the orbital kinematics, including trajectories and orbital eccentricities. We verify our code by evolving single stars and extracting radial perturbative modes, which compare very well to results from perturbation theory. The Einstein equations are solved in a first order reduction of the generalized harmonic formulation, and the fluid equations are solved using a modified convex essentially non-oscillatory method. All calculations are done in three spatial dimensions without symmetry assumptions. We use the \had computational infrastructure for distributed adaptive mesh refinement.Comment: 14 pages, 16 figures. Added one figure from previous version; corrected typo

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

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

Thermal phenomenology of hadrons from 200 AGeV S+S collisions

We develop a complete and consistent description for the hadron spectra from heavy ion collisions in terms of a few collective variables, in particular temperature, longitudinal and transverse flow. To achieve a meaningful comparison with presently available data, we also include the resonance decays into our picture. To disentangle the influences of transverse flow and resonance decays in the $m_T$-spectra, we analyse in detail the shape of the $m_T$-spectra.Comment: 31 pages, 13 figs in seperate uuencoded file, for LaTeX, epsf.sty and dvips, TPR-93-16 and BNL-(no number yet

Hydrodynamical assessment of 200 AGeV collisions

We are analyzing the hydrodynamics of 200 A GeV S+S collisions using a new approach which tries to quantify the uncertainties arising from the specific implementation of the hydrodynamical model. Based on a previous phenomenological analysis we use the global hydrodynamics model to show that the amount of initial flow, or initial energy density, cannot be determined from the hadronic momentum spectra. We additionally find that almost always a sizeable transverse flow deve- lops, which causes the system to freeze out, thereby limiting the flow velocity in itself. This freeze-out dominance in turn makes a distinction between a plasma and a hadron resonance gas equation of state very difficult, whereas a pure pion gas can easily be ruled out from present data. To complete the picture we also analyze particle multiplicity data, which suggest that chemical equilibrium is not reached with respect to the strange particles. However, the over- population of pions seems to be at most moderate, with a pion chemical potential far away from the Bose divergence.Comment: 19 pages, 11 figs in separate uuencoded file, for LateX, epsf.tex, dvips, TPR-94-5 and BNL-(no number yet

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

Spin observables of the reactions NN -> DeltaN and pd -> Delta (pp)(1S0) in collinear kinematics

A general formalism for double and triple spin-correlations of the reaction NN -> DeltaN is developed for the case of collinear kinematics. A complete polarization experiment allowing to reconstruct all of the four amplitudes describing this process is suggested. Furthermore, the spin observables of the inelastic charge-exchange reaction pd -> Delta^0(pp)(1S0) are analyzed in collinear kinematics within the single pN scattering mechanism involving the subprocess pn -> Delta^0p. The full set of spin observables related to the polarization of one or two initial particles and one final particle is obtained in terms of three invariant amplitudes of the reaction pd -> Delta (pp)(1S0) and the transition form factor d->(pp)(1S0). A complete polarization experiment for the reaction pd -> Delta^0(pp)(1S0) is suggested which allows one to determine three independent combinations of the four amplitudes of the elementary subprocess NN -> DeltaN.Comment: 12 pages, 1 figur

Excitons in quasi-one dimensional organics: Strong correlation approximation

An exciton theory for quasi-one dimensional organic materials is developed in the framework of the Su-Schrieffer-Heeger Hamiltonian augmented by short range extended Hubbard interactions. Within a strong electron-electron correlation approximation, the exciton properties are extensively studied. Using scattering theory, we analytically obtain the exciton energy and wavefunction and derive a criterion for the existence of a $B_u$ exciton. We also systematically investigate the effect of impurities on the coherent motion of an exciton. The coherence is measured by a suitably defined electron-hole correlation function. It is shown that, for impurities with an on-site potential, a crossover behavior will occur if the impurity strength is comparable to the bandwidth of the exciton, corresponding to exciton localization. For a charged impurity with a spatially extended potential, in addition to localization the exciton will dissociate into an uncorrelated electron-hole pair when the impurity is sufficiently strong to overcome the Coulomb interaction which binds the electron-hole pair. Interchain coupling effects are also discussed by considering two polymer chains coupled through nearest-neighbor interchain hopping $t_{\perp}$ and interchain Coulomb interaction $V_{\perp}$. Within the $t$ matrix scattering formalism, for every center-of-mass momentum, we find two poles determined only by $V_{\perp}$, which correspond to the interchain excitons. Finally, the exciton state is used to study the charge transfer from a polymer chain to an adjacent dopant molecule.Comment: 24 pages, 23 eps figures, pdf file of the paper availabl