Tracking Black Holes in Numerical Relativity

This work addresses and solves the problem of generically tracking black hole event horizons in computational simulation of black hole interactions. Solutions of the hyperbolic eikonal equation, solved on a curved spacetime manifold containing black hole sources, are employed in development of a robust tracking method capable of continuously monitoring arbitrary changes of topology in the event horizon, as well as arbitrary numbers of gravitational sources. The method makes use of continuous families of level set viscosity solutions of the eikonal equation with identification of the black hole event horizon obtained by the signature feature of discontinuity formation in the eikonal's solution. The method is employed in the analysis of the event horizon for the asymmetric merger in a binary black hole system. In this first such three dimensional analysis, we establish both qualitative and quantitative physics for the asymmetric collision; including: 1. Bounds on the topology of the throat connecting the holes following merger, 2. Time of merger, and 3. Continuous accounting for the surface of section areas of the black hole sources.Comment: 14 pages, 16 figure

Phylogenetic patterns recover known HIV epidemiological relationships and reveal common transmission of multiple variants.

The growth of human immunodeficiency virus (HIV) sequence databases resulting from drug resistance testing has motivated efforts using phylogenetic methods to assess how HIV spreads1-4. Such inference is potentially both powerful and useful for tracking the epidemiology of HIV and the allocation of resources to prevention campaigns. We recently used simulation and a small number of illustrative cases to show that certain phylogenetic patterns are associated with different types of epidemiological linkage5. Our original approach was later generalized for large next-generation sequencing datasets and implemented as a free computational pipeline6. Previous work has claimed that direction and directness of transmission could not be established from phylogeny because one could not be sure that there were no intervening or missing links involved7-9. Here, we address this issue by investigating phylogenetic patterns from 272 previously identified HIV transmission chains with 955 transmission pairs representing diverse geography, risk groups, subtypes, and genomic regions. These HIV transmissions had known linkage based on epidemiological information such as partner studies, mother-to-child transmission, pairs identified by contact tracing, and criminal cases. We show that the resulting phylogeny inferred from real HIV genetic sequences indeed reveals distinct patterns associated with direct transmission contra transmissions from a common source. Thus, our results establish how to interpret phylogenetic trees based on HIV sequences when tracking who-infected-whom, when and how genetic information can be used for improved tracking of HIV spread. We also investigate limitations that stem from limited sampling and genetic time-trends in the donor and recipient HIV populations

A robust method for calculating interface curvature and normal vectors using an extracted local level set

The level-set method is a popular interface tracking method in two-phase flow simulations. An often-cited reason for using it is that the method naturally handles topological changes in the interface, e.g. merging drops, due to the implicit formulation. It is also said that the interface curvature and normal vectors are easily calculated. This last point is not, however, the case in the moments during a topological change, as several authors have already pointed out. Various methods have been employed to circumvent the problem. In this paper, we present a new such method which retains the implicit level-set representation of the surface and handles general interface configurations. It is demonstrated that the method extends easily to 3D. The method is validated on static interface configurations, and then applied to two-phase flow simulations where the method outperforms the standard method and the results agree well with experiments.Comment: 31 pages, 18 figure

Experimental determination of the complete spin structure for anti-proton + proton -> anti-\Lambda + \Lambda at anti-proton beam momentum of 1.637 GeV/c

The reaction anti-proton + proton -> anti-\Lambda + \Lambda -> anti-proton + \pi^+ + proton + \pi^- has been measured with high statistics at anti-proton beam momentum of 1.637 GeV/c. The use of a transversely-polarized frozen-spin target combined with the self-analyzing property of \Lambda/anti-\Lambda decay allows access to unprecedented information on the spin structure of the interaction. The most general spin-scattering matrix can be written in terms of eleven real parameters for each bin of scattering angle, each of these parameters is determined with reasonable precision. From these results all conceivable spin-correlations are determined with inherent self-consistency. Good agreement is found with the few previously existing measurements of spin observables in anti-proton + proton -> anti-\Lambda + \Lambda near this energy. Existing theoretical models do not give good predictions for those spin-observables that had not been previously measured.Comment: To be published in Phys. Rev. C. Tables of results (i.e. Ref. 24) are available at http://www-meg.phys.cmu.edu/~bquinn/ps185_pub/results.tab 24 pages, 16 figure

Solar Magnetic Tracking. I. Software Comparison and Recommended Practices

Feature tracking and recognition are increasingly common tools for data analysis, but are typically implemented on an ad-hoc basis by individual research groups, limiting the usefulness of derived results when selection effects and algorithmic differences are not controlled. Specific results that are affected include the solar magnetic turnover time, the distributions of sizes, strengths, and lifetimes of magnetic features, and the physics of both small scale flux emergence and the small-scale dynamo. In this paper, we present the results of a detailed comparison between four tracking codes applied to a single set of data from SOHO/MDI, describe the interplay between desired tracking behavior and parameterization of tracking algorithms, and make recommendations for feature selection and tracking practice in future work.Comment: In press for Astrophys. J. 200

Generic Tracking of Multiple Apparent Horizons with Level Flow

We report the development of the first apparent horizon locator capable of finding multiple apparent horizons in a ``generic'' numerical black hole spacetime. We use a level-flow method which, starting from a single arbitrary initial trial surface, can undergo topology changes as it flows towards disjoint apparent horizons if they are present. The level flow method has two advantages: 1) The solution is independent of changes in the initial guess and 2) The solution can have multiple components. We illustrate our method of locating apparent horizons by tracking horizon components in a short Kerr-Schild binary black hole grazing collision.Comment: 13 pages including figures, submitted to Phys Rev

Searches for Gauge Mediated SUSY Breaking at LEP

Searches for neutralinos and sleptons with arbitrary lifetimes, as predicted in the framework of the Gauge Mediated Supersymmetry Breaking (GMSB) model, have been performed by the four LEP collaborations. No evidence for these particles has been found in the data recorded at center-of-mass energies up to sqrt(s)=209GeV. Therefore constraints on the production cross-sections and particle masses as well as interpretations in the framework of the GMSB model are presented.Comment: 6 pages, 4 figures, To appear in the proceedings of the 37th Rencontres de Moriond, Electroweak Interactions and Unified Theories, Les Arcs, France, 9-16 March 200