Ninja data analysis with a detection pipeline based on the Hilbert-Huang Transform

The Ninja data analysis challenge allowed the study of the sensitivity of data analysis pipelines to binary black hole numerical relativity waveforms in simulated Gaussian noise at the design level of the LIGO observatory and the VIRGO observatory. We analyzed NINJA data with a pipeline based on the Hilbert Huang Transform, utilizing a detection stage and a characterization stage: detection is performed by triggering on excess instantaneous power, characterization is performed by displaying the kernel density enhanced (KD) time-frequency trace of the signal. Using the simulated data based on the two LIGO detectors, we were able to detect 77 signals out of 126 above SNR 5 in coincidence, with 43 missed events characterized by signal to noise ratio SNR less than 10. Characterization of the detected signals revealed the merger part of the waveform in high time and frequency resolution, free from time-frequency uncertainty. We estimated the timelag of the signals between the detectors based on the optimal overlap of the individual KD time-frequency maps, yielding estimates accurate within a fraction of a millisecond for half of the events. A coherent addition of the data sets according to the estimated timelag eventually was used in a characterization of the event.Comment: Accepted for publication in CQG, special issue NRDA proceedings 200

Two-dimensional quantum black holes: Numerical methods

We present details of a new numerical code designed to study the formation and evaporation of 2-dimensional black holes within the CGHS model. We explain several elements of the scheme that are crucial to resolve the late-time behavior of the spacetime, including regularization of the field variables, compactification of the coordinates, the algebraic form of the discretized equations of motion, and the use of a modified Richardson extrapolation scheme to achieve high-order convergence. Physical interpretation of our results will be discussed in detail elsewhere

Requirements elicitation interviews and applications for an underground remote-piloted aerial system

Abstract: Being underground is dangerous, it really is. People shouldn’t go there if they don’t need to, especially in an age where driverless cars by Google and shopping delivery by Amazon are on the brink of implementation in unstructured public environments; but held up, not by technology, but by wrangling about legislation covering liability determination, and fears of hacking drive-by-wire systems. Underground mining is a secure (from an IT perspective), structured, and predictable environment with limited and controlled human access – an ideal place for automated or remotely piloted systems. This paper outlines the semi-structured interview process that was executed to determine the requirements for an underground remotely piloted aerial system (URPAS). The potential applications are explored; including search and rescue, the business-orientated activity of regularly scanning rockpasses to predict and prevent blockages, and scanning a blocked rockpass to determine the blockage location and structure. An even more ’out-of-the-box’ application, using the vehicle to deliver explosives to the underside of a blockage, was included in the brainstorming discussions. Interviews with ten mining and unmanned aerial systems experts were conducted with a questionnaire as the primary data collection tool. The questionnaires were analysed to determine the representatively of the sample set, and therefore the validity of the gathered data, based upon expertise ratings in each of the relevant areas of knowledge: mining, surveying and mapping, and remote-piloted aerial systems. The goal was to identify the key performance requirements of a U-RPAS, and determine the feasibility of such a system being developed. A specialized company providing a scanning service emerged the preferred implementation method, and the rationale for this choice is presented. As context, the sub-system prototypes used in the brainstorming section of the interviews are presented, as are the implementation scenarios discussed in the interviews. It is thus shown that this method of requirements elicitation is suitable for this type of technology implementation project

Comparisons of binary black hole merger waveforms

This a particularly exciting time for gravitational wave physics. Ground-based gravitational wave detectors are now operating at a sensitivity such that gravitational radiation may soon be directly detected, and recently several groups have independently made significant breakthroughs that have finally enabled numerical relativists to solve the Einstein field equations for coalescing black-hole binaries, a key source of gravitational radiation. The numerical relativity community is now in the position to begin providing simulated merger waveforms for use by the data analysis community, and it is therefore very important that we provide ways to validate the results produced by various numerical approaches. Here, we present a simple comparison of the waveforms produced by two very different, but equally successful approaches--the generalized harmonic gauge and the moving puncture methods. We compare waveforms of equal-mass black hole mergers with minimal or vanishing spins. The results show exceptional agreement for the final burst of radiation, with some differences attributable to small spins on the black holes in one case.Comment: Revtex 4, 5 pages. Published versio

A Dynamical Systems Approach to Schwarzschild Null Geodesics

The null geodesics of a Schwarzschild black hole are studied from a dynamical systems perspective. Written in terms of Kerr-Schild coordinates, the null geodesic equation takes on the simple form of a particle moving under the influence of a Newtonian central force with an inverse-cubic potential. We apply a McGehee transformation to these equations, which clearly elucidates the full phase space of solutions. All the null geodesics belong to one of four families of invariant manifolds and their limiting cases, further characterized by the angular momentum L of the orbit: for |L|>|L_c|, (1) the set that flow outward from the white hole, turn around, then fall into the black hole, (2) the set that fall inward from past null infinity, turn around outside the black hole to continue to future null infinity, and for |L|<|L_c|, (3) the set that flow outward from the white hole and continue to future null infinity, (4) the set that flow inward from past null infinity and into the black hole. The critical angular momentum Lc corresponds to the unstable circular orbit at r=3M, and the homoclinic orbits associated with it. There are two additional critical points of the flow at the singularity at r=0. Though the solutions of geodesic motion and Hamiltonian flow we describe here are well known, what we believe is a novel aspect of this work is the mapping between the two equivalent descriptions, and the different insights each approach can give to the problem. For example, the McGehee picture points to a particularly interesting limiting case of the class (1) that move from the white to black hole: in the limit as L goes to infinity, as described in Schwarzschild coordinates, these geodesics begin at r=0, flow along t=constant lines, turn around at r=2M, then continue to r=0. During this motion they circle in azimuth exactly once, and complete the journey in zero affine time.Comment: 14 pages, 3 Figure

Black Hole Mergers and Unstable Circular Orbits

We describe recent numerical simulations of the merger of a class of equal mass, non-spinning, eccentric binary black hole systems in general relativity. We show that with appropriate fine-tuning of the initial conditions to a region of parameter space we denote the threshold of immediate merger, the binary enters a phase of close interaction in a near-circular orbit, stays there for an amount of time proportional to logarithmic distance from the threshold in parameter space, then either separates or merges to form a single Kerr black hole. To gain a better understanding of this phenomena we study an analogous problem in the evolution of equatorial geodesics about a central Kerr black hole. A similar threshold of capture exists for appropriate classes of initial conditions, and tuning to threshold the geodesics approach one of the unstable circular geodesics of the Kerr spacetime. Remarkably, with a natural mapping of the parameters of the geodesic to that of the equal mass system, the scaling exponent describing the whirl phase of each system turns out to be quite similar. Armed with this lone piece of evidence that an approximate correspondence might exist between near-threshold evolution of geodesics and generic binary mergers, we illustrate how this information can be used to estimate the cross section and energy emitted in the ultra relativistic black hole scattering problem. This could eventually be of use in providing estimates for the related problem of parton collisions at the Large Hadron Collider in extra dimension scenarios where black holes are produced.Comment: 16 pages, 12 figures; updated to coincide with journal versio

Requirements degradation for the creation of a first prototype

The requirements engineering process is typically executed, irrespective of the process model chosen, for the final commercially viable system. The system requirements generated are for a system deployed and used in its final form and function. However, the first prototype that is generated is typically representative of a minimum viable technology, and represents a degraded set of the initial system requirements specification. Typically, a first prototype is used as a technology demonstrator, and its failure or success will determine the continuation of the project, with success triggering the allocation of additional financial and personal resources. This paper explores techniques for requirements degradation that can be used to form the system requirements specification for the first prototype. A requirements Engineering methodology is proposed based upon a survey of literature. It takes into consideration the characteristics of the project, i.e. a market driven, technology implementation research project with limited budget and a flexible timeline executed in an academic environment. The techniques must take into cognizance the main risk items, and core requirements, that need to be demonstrated in the minimum viable technology to secure the future of the project. The degradation cannot undermine or jeopardize the future success of the commercially viable system in determining the subset of requirements for the minimum viable technology

Malaria intervention scale-up in Africa : effectiveness predictions for health programme planning tools, based on dynamic transmission modelling

Scale-up of malaria prevention and treatment needs to continue to further important gains made in the past decade, but national strategies and budget allocations are not always evidence-based. Statistical models were developed summarizing dynamically simulated relations between increases in coverage and intervention impact, to inform a malaria module in the Spectrum health programme planning tool.; The dynamic Plasmodium falciparum transmission model OpenMalaria was used to simulate health effects of scale-up of insecticide-treated net (ITN) usage, indoor residual spraying (IRS), management of uncomplicated malaria cases (CM) and seasonal malaria chemoprophylaxis (SMC) over a 10-year horizon, over a range of settings with stable endemic malaria. Generalized linear regression models (GLMs) were used to summarize determinants of impact across a range of sub-Sahara African settings.; Selected (best) GLMs explained 94-97 % of variation in simulated post-intervention parasite infection prevalence, 86-97 % of variation in case incidence (three age groups, three 3-year horizons), and 74-95 % of variation in malaria mortality. For any given effective population coverage, CM and ITNs were predicted to avert most prevalent infections, cases and deaths, with lower impacts for IRS, and impacts of SMC limited to young children reached. Proportional impacts were larger at lower endemicity, and (except for SMC) largest in low-endemic settings with little seasonality. Incremental health impacts for a given coverage increase started to diminish noticeably at above ~40 % coverage, while in high-endemic settings, CM and ITNs acted in synergy by lowering endemicity. Vector control and CM, by reducing endemicity and acquired immunity, entail a partial rebound in malaria mortality among people above 5 years of age from around 5-7 years following scale-up. SMC does not reduce endemicity, but slightly shifts malaria to older ages by reducing immunity in child cohorts reached.; Health improvements following malaria intervention scale-up vary with endemicity, seasonality, age and time. Statistical models can emulate epidemiological dynamics and inform strategic planning and target setting for malaria control

Simulation of Binary Black Hole Spacetimes with a Harmonic Evolution Scheme

A numerical solution scheme for the Einstein field equations based on generalized harmonic coordinates is described, focusing on details not provided before in the literature and that are of particular relevance to the binary black hole problem. This includes demonstrations of the effectiveness of constraint damping, and how the time slicing can be controlled through the use of a source function evolution equation. In addition, some results from an ongoing study of binary black hole coalescence, where the black holes are formed via scalar field collapse, are shown. Scalar fields offer a convenient route to exploring certain aspects of black hole interactions, and one interesting, though tentative suggestion from this early study is that behavior reminiscent of "zoom-whirl" orbits in particle trajectories is also present in the merger of equal mass, non-spinning binaries, with appropriately fine-tuned initial conditions.Comment: 16 pages, 14 figures; replaced with published versio