Cauchy-perturbative matching and outer boundary conditions I: Methods and tests

We present a new method of extracting gravitational radiation from three-dimensional numerical relativity codes and providing outer boundary conditions. Our approach matches the solution of a Cauchy evolution of Einstein's equations to a set of one-dimensional linear wave equations on a curved background. We illustrate the mathematical properties of our approach and discuss a numerical module we have constructed for this purpose. This module implements the perturbative matching approach in connection with a generic three-dimensional numerical relativity simulation. Tests of its accuracy and second-order convergence are presented with analytic linear wave data.Comment: 13 pages, 6 figures, RevTe

Cauchy-perturbative matching and outer boundary conditions: computational studies

We present results from a new technique which allows extraction of gravitational radiation information from a generic three-dimensional numerical relativity code and provides stable outer boundary conditions. In our approach we match the solution of a Cauchy evolution of the nonlinear Einstein field equations to a set of one-dimensional linear equations obtained through perturbation techniques over a curved background. We discuss the validity of this approach in the case of linear and mildly nonlinear gravitational waves and show how a numerical module developed for this purpose is able to provide an accurate and numerically convergent description of the gravitational wave propagation and a stable numerical evolution.Comment: 20 pages, RevTe

Waveform propagation in black hole spacetimes: evaluating the quality of numerical solutions

We compute the propagation and scattering of linear gravitational waves off a Schwarzschild black hole using a numerical code which solves a generalization of the Zerilli equation to a three dimensional cartesian coordinate system. Since the solution to this problem is well understood it represents a very good testbed for evaluating our ability to perform three dimensional computations of gravitational waves in spacetimes in which a black hole event horizon is present.Comment: 13 pages, RevTeX, to appear in Phys. Rev.

An exact solution for 2+1 dimensional critical collapse

We find an exact solution in closed form for the critical collapse of a scalar field with cosmological constant in 2+1 dimensions. This solution agrees with the numerical simulation done by Pretorius and Choptuik of this system.Comment: 5 pages, 5 figures, Revtex. New comparison of analytic and numerical solutions beyond the past light cone of the singularity added. Two new references added. Error in equation (21) correcte

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

The metal-insulator transition in Si:X: Anomalous response to a magnetic field

The zero-temperature magnetoconductivity of just-metallic Si:P scales with magnetic field, H, and dopant concentration, n, lying on a single universal curve. We note that Si:P, Si:B, and Si:As all have unusually large magnetic field crossover exponents near 2, and suggest that this anomalously weak response to a magnetic field is a common feature of uncompensated doped semiconductors.Comment: 4 pages (including figures

Understanding initial data for black hole collisions

Numerical relativity, applied to collisions of black holes, starts with initial data for black holes already in each other's strong field. The initial hypersurface data typically used for computation is based on mathematical simplifying prescriptions, such as conformal flatness of the 3-geometry and longitudinality of the extrinsic curvature. In the case of head on collisions of equal mass holes, there is evidence that such prescriptions work reasonably well, but it is not clear why, or whether this success is more generally valid. Here we study these questions by considering the ``particle limit'' for head on collisions of nonspinning holes. Einstein's equations are linearized in the mass of the small hole, and described by a single gauge invariant spacetime function psi, for each multipole. The resulting equations have been solved by numerical evolution for collisions starting from various initial separations, and the evolution is studied on a sequence of hypersurfaces. In particular, we extract hypersurface data, that is psi and its time derivative, on surfaces of constant background Schwarzschild time. These evolved data can then be compared with ``prescribed'' data, evolved data can be replaced by prescribed data on any hypersurface, and evolved further forward in time, a gauge invariant measure of deviation from conformal flatness can be evaluated, etc. The main findings of this study are: (i) For holes of unequal mass the use of prescribed data on late hypersurfaces is not successful. (ii) The failure is likely due to the inability of the prescribed data to represent the near field of the smaller hole. (iii) The discrepancy in the extrinsic curvature is more important than in the 3-geometry. (iv) The use of the more general conformally flat longitudinal data does not notably improve this picture.Comment: 20 pages, REVTEX, 26 PS figures include

Statistics of quantum transmission in one dimension with broad disorder

We study the statistics of quantum transmission through a one-dimensional disordered system modelled by a sequence of independent scattering units. Each unit is characterized by its length and by its action, which is proportional to the logarithm of the transmission probability through this unit. Unit actions and lengths are independent random variables, with a common distribution that is either narrow or broad. This investigation is motivated by results on disordered systems with non-stationary random potentials whose fluctuations grow with distance. In the statistical ensemble at fixed total sample length four phases can be distinguished, according to the values of the indices characterizing the distribution of the unit actions and lengths. The sample action, which is proportional to the logarithm of the conductance across the sample, is found to obey a fluctuating scaling law, and therefore to be non-self-averaging, in three of the four phases. According to the values of the two above mentioned indices, the sample action may typically grow less rapidly than linearly with the sample length (underlocalization), more rapidly than linearly (superlocalization), or linearly but with non-trivial sample-to-sample fluctuations (fluctuating localization).Comment: 26 pages, 4 figures, 1 tabl

Negative Magnetoresistance of Granular Metals in a Strong Magnetic Field

The magnetoresistance of a granular superconductor in a strong magnetic field destroying the gap in each grain is considered. It is assumed that the tunneling between grains is sufficiently large such that all conventional effects of localization can be neglected. A non-trivial sensitivity to the magnetic field comes from superconducting fluctuations leading to the formation of virtual Cooper pairs and reducing the density of states. At low temperature, the pairs do not contribute to the macroscopic transport but their existence can drastically reduce the conductivity. Growing the magnetic field one destroys the fluctuations, which improves the metallic properties and leads to the negative magnetoresistance.Comment: 4 pages, 1 figure, RevTe

Positive Magneto-Resistance in Quasi-1D Conductors

We present here a simple qualitative model that interpolates between the high and low temperature properties of quasi-1D conductors. At high temperatures we argue that transport is governed by inelastic scattering whereas at low temperatures the conductance decays exponentially with the electron dephasing length. The crossover between these regimes occurs at the temperature at which the elastic and inelastic scattering times become equal. This model is shown to be in quantitative agreement with the organic conductor $TTT_2I_{3-\delta}$. Within this model, we also show that on the insulating side, the positive magnetoresistance of the form $(H/T)^2$ observed in $TTT_2I_{3-\delta}$ and other quasi-1D conductors can be explained by the role spin-flip scattering plays in the electron dephasing rate.Comment: 4 pages, Latex, no figure