Gravitational waves from black hole collisions via an eclectic approach

We present the first results in a new program intended to make the best use of all available technologies to provide an effective understanding of waves from inspiralling black hole binaries in time for imminent observations. In particular, we address the problem of combining the close-limit approximation describing ringing black holes and full numerical relativity, required for essentially nonlinear interactions. We demonstrate the effectiveness of our approach using general methods for a model problem, the head-on collision of black holes. Our method allows a more direct physical understanding of these collisions indicating clearly when non-linear methods are important. The success of this method supports our expectation that this unified approach will be able to provide astrophysically relevant results for black hole binaries in time to assist gravitational wave observations.Comment: 4 pages, 3 eps figures, Revte

Improved Implementation of Point Location in General Two-Dimensional Subdivisions

We present a major revamp of the point-location data structure for general two-dimensional subdivisions via randomized incremental construction, implemented in CGAL, the Computational Geometry Algorithms Library. We can now guarantee that the constructed directed acyclic graph G is of linear size and provides logarithmic query time. Via the construction of the Voronoi diagram for a given point set S of size n, this also enables nearest-neighbor queries in guaranteed O(log n) time. Another major innovation is the support of general unbounded subdivisions as well as subdivisions of two-dimensional parametric surfaces such as spheres, tori, cylinders. The implementation is exact, complete, and general, i.e., it can also handle non-linear subdivisions. Like the previous version, the data structure supports modifications of the subdivision, such as insertions and deletions of edges, after the initial preprocessing. A major challenge is to retain the expected O(n log n) preprocessing time while providing the above (deterministic) space and query-time guarantees. We describe an efficient preprocessing algorithm, which explicitly verifies the length L of the longest query path in O(n log n) time. However, instead of using L, our implementation is based on the depth D of G. Although we prove that the worst case ratio of D and L is Theta(n/log n), we conjecture, based on our experimental results, that this solution achieves expected O(n log n) preprocessing time.Comment: 21 page

Topological quantum many-body scars in quantum dimer models on the kagome lattice

We present a class of quantum dimer models on the kagome lattice with full translational invariance that feature a quantum many-body scar state of analytically known entanglement properties within their spectra. Using exact diagonalization on lattices of up to 60 sites, we show that nonscar states conform to the eigenstate thermalization hypothesis. Specifically, we show that energies are distributed according to the Gaussian ensemble expected of their respective symmetry sector, illustrate the existence of the scar from bipartite entanglement properties, and demonstrate revival phenomena in studies of fidelity dynamics

Nonlinear and Perturbative Evolution of Distorted Black Holes; 2, Odd-parity Modes

We compare the fully nonlinear and perturbative evolution of nonrotating black holes with odd-parity distortions utilizing the perturbative results to interpret the nonlinear results. This introduction of the second polarization (odd-parity) mode of the system, and the systematic use of combined techniques brings us closer to the goal of studying more complicated systems like distorted, rotating black holes, such as those formed in the final inspiral stage of two black holes. The nonlinear evolutions are performed with the 3D parallel code for Numerical Relativity, {Cactus}, and an independent axisymmetric code, {Magor}. The linearized calculation is performed in two ways: (a) We treat the system as a metric perturbation on Schwarzschild, using the Regge-Wheeler equation to obtain the waveforms produced. (b) We treat the system as a curvature perturbation of a Kerr black hole (but here restricted to the case of vanishing rotation parameter a) and evolve it with the Teukolsky equation The comparisons of the waveforms obtained show an excellent agreement in all cases

MUBs inequivalence and affine planes

There are fairly large families of unitarily inequivalent complete sets of N+1 mutually unbiased bases (MUBs) in C^N for various prime powers N. The number of such sets is not bounded above by any polynomial as a function of N. While it is standard that there is a superficial similarity between complete sets of MUBs and finite affine planes, there is an intimate relationship between these large families and affine planes. This note briefly summarizes "old" results that do not appear to be well-known concerning known families of complete sets of MUBs and their associated planes.Comment: This is the version of this paper appearing in J. Mathematical Physics 53, 032204 (2012) except for format changes due to the journal's style policie

Unstable fields in Kerr spacetimes

We show that both the interior region $r<M-\sqrt{M^2-a^2}$ of a Kerr black hole and the $a^2>M^2$ Kerr naked singularity admit unstable solutions of the Teukolsky equation for any value of the spin weight. For every harmonic number there is at least one axially symmetric mode that grows exponentially in time and decays properly in the radial directions. These can be used as Debye potentials to generate solutions for the scalar, Weyl spinor, Maxwell and linearized gravity field equations on these backgrounds, satisfying appropriate spatial boundary conditions and growing exponentially in time, as shown in detail for the Maxwell case. It is suggested that the existence of the unstable modes is related to the so called "time machine" region, where the axial Killing vector field is time-like, and the Teukolsky equation, restricted to axially symmetric fields, changes its character from hyperbolic to elliptic

Measurement of adsorption of a single component from the liquid phase : modelling investigation and sensitivity analysis

In this work, we consider an alternative approach for the measurement of adsorption from the liquid phase. Consider a mixture consisting of a non-adsorbed component (B) and an adsorbed component (A) present at some low concentration. Initially, a feed of component B only flows through a column packed with an adsorbent. Then, the feed is switched to the mixture of A and B. As soon as the mixture enters the column, there will be a reduction in the outlet flow rate as component A leaves the liquid phase and passes into the adsorbed phase. There are three stages to this work. The first is to develop overall and component balances to show how the amount adsorbed of component A can be determined from the variation in the column outlet flow rate. The second is to determine the actual variation in the column outlet flow rate for both plug flow and axial-dispersed plug flow. The final stage is to consider the suitability of a gravity-fed system to deliver the feed to the column. An analysis of the results shows that the experimental arrangement should be able to accurately monitor adsorption from the liquid phase where the mass fraction of the solute is of the order of 1%: the limiting experimental factor is how constant the volumetric flow rate of the liquid feed can be maintained

Conductance of Tomonaga-Luttinger liquid wires and junctions with resistances

We study the effect that resistive regions have on the conductance of a quantum wire with interacting electrons which is connected to Fermi liquid leads. Using the bosonization formalism and a Rayleigh dissipation function to model the power dissipation, we use both scattering theory and Green's function techniques to derive the DC conductance. The resistive regions are generally found to lead to incoherent transport. For a single wire, we find that the resistance adds in series to the contact resistance of h/e^2 for spinless electrons, and the total resistance is independent of the Luttinger parameter K_W of the wire. We numerically solve the bosonic equations to illustrate what happens when a charge density pulse is incident on the wire; the results depend on the parameters of the resistive and interacting regions in interesting ways. For a junction of Tomonaga-Luttinger liquid wires, we use a dissipationless current splitting matrix to model the junction. For a junction of three wires connected to Fermi liquid leads, there are two families of such matrices; we find that the conductance matrix generally depends on K_W for one family but is independent of K_W for the other family, regardless of the resistances present in the system.Comment: 6 pages, 3 figures; added a discussion of time reversal invariance; this is the published versio

Black Hole--Scalar Field Interactions in Spherical Symmetry

We examine the interactions of a black hole with a massless scalar field using a coordinate system which extends ingoing Eddington-Finkelstein coordinates to dynamic spherically symmetric-spacetimes. We avoid problems with the singularity by excising the region of the black hole interior to the apparent horizon. We use a second-order finite difference scheme to solve the equations. The resulting program is stable and convergent and will run forever without problems. We are able to observe quasi-normal ringing and power-law tails as well an interesting nonlinear feature.Comment: 16 pages, 26 figures, RevTex, to appear in Phys. Rev.

S=1/2 chains and spin-Peierls transition in TiOCl

We study TiOCl as an example of an S=1/2 layered Mott insulator. From our analysis of new susceptibility data, combined with LDA and LDA+U band structure calculations, we conclude that orbital ordering produces quasi-one-dimensional spin chains and that TiOCl is a new example of Heisenberg-chains which undergo a spin-Peierls transition. The energy scale is an order of magnitude larger than that of previously known examples. The effects of non-magnetic Sc impurities are explained using a model of broken finite chains.Comment: 5 pages, 5 figures (color); details on crystal growth added; to be published in Phys. Rev.