A New Generalized Harmonic Evolution System

A new representation of the Einstein evolution equations is presented that is first order, linearly degenerate, and symmetric hyperbolic. This new system uses the generalized harmonic method to specify the coordinates, and exponentially suppresses all small short-wavelength constraint violations. Physical and constraint-preserving boundary conditions are derived for this system, and numerical tests that demonstrate the effectiveness of the constraint suppression properties and the constraint-preserving boundary conditions are presented.Comment: Updated to agree with published versio

Implementation of higher-order absorbing boundary conditions for the Einstein equations

We present an implementation of absorbing boundary conditions for the Einstein equations based on the recent work of Buchman and Sarbach. In this paper, we assume that spacetime may be linearized about Minkowski space close to the outer boundary, which is taken to be a coordinate sphere. We reformulate the boundary conditions as conditions on the gauge-invariant Regge-Wheeler-Zerilli scalars. Higher-order radial derivatives are eliminated by rewriting the boundary conditions as a system of ODEs for a set of auxiliary variables intrinsic to the boundary. From these we construct boundary data for a set of well-posed constraint-preserving boundary conditions for the Einstein equations in a first-order generalized harmonic formulation. This construction has direct applications to outer boundary conditions in simulations of isolated systems (e.g., binary black holes) as well as to the problem of Cauchy-perturbative matching. As a test problem for our numerical implementation, we consider linearized multipolar gravitational waves in TT gauge, with angular momentum numbers l=2 (Teukolsky waves), 3 and 4. We demonstrate that the perfectly absorbing boundary condition B_L of order L=l yields no spurious reflections to linear order in perturbation theory. This is in contrast to the lower-order absorbing boundary conditions B_L with L<l, which include the widely used freezing-Psi_0 boundary condition that imposes the vanishing of the Newman-Penrose scalar Psi_0.Comment: 25 pages, 9 figures. Minor clarifications. Final version to appear in Class. Quantum Grav

On the Hoop conjecture and the weak cosmic censorship conjecture for the axisymmetric Einstein-Vlasov system

We consider gravitational collapse for the axially symmetric Einstein-Vlasov system. We investigate the weak cosmic censorship conjecture in the case of highly prolate initial data and we investigate the ``only if" part of the Hoop conjecture. Shapiro and Teukolsky initiated a similar study in 1991 \cite{Shapiro1991} where they found support that the weak cosmic censorship conjecture was violated for sufficiently prolate spheroidal initial data. More recently, independent studies of this problem have been carried out by Yoo, Harada and Okawa \cite{Yoo2017} and by East \cite{East2019}. A common feature in these works is that the initial data are dust-like. Dust can be considered as a \textit{singular} case of matter described by the Einstein-Vlasov system. The original motivation by Shapiro and Teukolsky to study this problem is based on the Lin-Mestel-Shu instability for gravitational collapse of uniform spheroids in the case of dust in Newtonian gravity. We argue that the Lin-Mestel-Shu solution is not relevant for studying the weak cosmic censorship of the Einstein-Vlasov system and we argue that dust-like initial data is also not relevant. To investigate collapse of highly prolate spheroidal configurations for the Einstein-Vlasov system is nevertheless interesting in view of the Hoop conjecture. By choosing highly prolate initial data the weak cosmic censorship conjecture is seriously tested. We carry out such a study for initial data which are not dust-like. We find formation of an apparent horizon in all cases we consider, which provides support for the weak cosmic censorship conjecture. In our tests of the Hoop conjecture we compute the polar circumference $\mathcal{C}_{H,p}$ at the time when the apparent horizon forms and find that it is less than 12% above $4\pi M$, where $M$ is the irreducible mass of the apparent horizon, which agrees with the spirit of the Hoop conjecture.Comment: 18 pages, 7 figure

Testing outer boundary treatments for the Einstein equations

Various methods of treating outer boundaries in numerical relativity are compared using a simple test problem: a Schwarzschild black hole with an outgoing gravitational wave perturbation. Numerical solutions computed using different boundary treatments are compared to a `reference' numerical solution obtained by placing the outer boundary at a very large radius. For each boundary treatment, the full solutions including constraint violations and extracted gravitational waves are compared to those of the reference solution, thereby assessing the reflections caused by the artificial boundary. These tests use a first-order generalized harmonic formulation of the Einstein equations. Constraint-preserving boundary conditions for this system are reviewed, and an improved boundary condition on the gauge degrees of freedom is presented. Alternate boundary conditions evaluated here include freezing the incoming characteristic fields, Sommerfeld boundary conditions, and the constraint-preserving boundary conditions of Kreiss and Winicour. Rather different approaches to boundary treatments, such as sponge layers and spatial compactification, are also tested. Overall the best treatment found here combines boundary conditions that preserve the constraints, freeze the Newman-Penrose scalar Psi_0, and control gauge reflections.Comment: Modified to agree with version accepted for publication in Class. Quantum Gra

An axisymmetric evolution code for the Einstein equations on hyperboloidal slices

We present the first stable dynamical numerical evolutions of the Einstein equations in terms of a conformally rescaled metric on hyperboloidal hypersurfaces extending to future null infinity. Axisymmetry is imposed in order to reduce the computational cost. The formulation is based on an earlier axisymmetric evolution scheme, adapted to time slices of constant mean curvature. Ideas from a previous study by Moncrief and the author are applied in order to regularize the formally singular evolution equations at future null infinity. Long-term stable and convergent evolutions of Schwarzschild spacetime are obtained, including a gravitational perturbation. The Bondi news function is evaluated at future null infinity.Comment: 21 pages, 4 figures. Minor additions, updated to agree with journal versio

Solving Einstein's Equations With Dual Coordinate Frames

A method is introduced for solving Einstein's equations using two distinct coordinate systems. The coordinate basis vectors associated with one system are used to project out components of the metric and other fields, in analogy with the way fields are projected onto an orthonormal tetrad basis. These field components are then determined as functions of a second independent coordinate system. The transformation to the second coordinate system can be thought of as a mapping from the original ``inertial'' coordinate system to the computational domain. This dual-coordinate method is used to perform stable numerical evolutions of a black-hole spacetime using the generalized harmonic form of Einstein's equations in coordinates that rotate with respect to the inertial frame at infinity; such evolutions are found to be generically unstable using a single rotating coordinate frame. The dual-coordinate method is also used here to evolve binary black-hole spacetimes for several orbits. The great flexibility of this method allows comoving coordinates to be adjusted with a feedback control system that keeps the excision boundaries of the holes within their respective apparent horizons.Comment: Updated to agree with published versio

Stable radiation-controlling boundary conditions for the generalized harmonic Einstein equations

This paper is concerned with the initial-boundary value problem for the Einstein equations in a first-order generalized harmonic formulation. We impose boundary conditions that preserve the constraints and control the incoming gravitational radiation by prescribing data for the incoming fields of the Weyl tensor. High-frequency perturbations about any given spacetime (including a shift vector with subluminal normal component) are analyzed using the Fourier-Laplace technique. We show that the system is boundary-stable. In addition, we develop a criterion that can be used to detect weak instabilities with polynomial time dependence, and we show that our system does not suffer from such instabilities. A numerical robust stability test supports our claim that the initial-boundary value problem is most likely to be well-posed even if nonzero initial and source data are included.Comment: 27 pages, 4 figures; more numerical results and references added, several minor amendments; version accepted for publication in Class. Quantum Gra

A point process framework for modeling electrical stimulation of the auditory nerve

Model-based studies of auditory nerve responses to electrical stimulation can provide insight into the functioning of cochlear implants. Ideally, these studies can identify limitations in sound processing strategies and lead to improved methods for providing sound information to cochlear implant users. To accomplish this, models must accurately describe auditory nerve spiking while avoiding excessive complexity that would preclude large-scale simulations of populations of auditory nerve fibers and obscure insight into the mechanisms that influence neural encoding of sound information. In this spirit, we develop a point process model of the auditory nerve that provides a compact and accurate description of neural responses to electric stimulation. Inspired by the framework of generalized linear models, the proposed model consists of a cascade of linear and nonlinear stages. We show how each of these stages can be associated with biophysical mechanisms and related to models of neuronal dynamics. Moreover, we derive a semi-analytical procedure that uniquely determines each parameter in the model on the basis of fundamental statistics from recordings of single fiber responses to electric stimulation, including threshold, relative spread, jitter, and chronaxie. The model also accounts for refractory and summation effects that influence the responses of auditory nerve fibers to high pulse rate stimulation. Throughout, we compare model predictions to published physiological data and explain differences in auditory nerve responses to high and low pulse rate stimulation. We close by performing an ideal observer analysis of simulated spike trains in response to sinusoidally amplitude modulated stimuli and find that carrier pulse rate does not affect modulation detection thresholds.Comment: 1 title page, 27 manuscript pages, 14 figures, 1 table, 1 appendi

Reliability and validity of the Dutch dimensional assessment of personality pathology-short form(DAPP-SF), a shortened version of the DAPP-Basic questionnaire

The Dimensional Assessment of Personality Pathology-Basic Questionnaire (DAPP-BQ) appears to be a good choice for the assessment of personality pathology. However, due to its length, administration of the instrument is rather time-consuming, hindering standard inclusion of the DABB-BQ in a battery of assessment instruments at intake. We developed the 136-item DAPP-SF (Short Form), and investigated its psychometric characteristics in various samples, i.e., a community-based sample (n = 487), patients with mood-, anxiety-, and somatoform disorders (n = 1,329), and patients with personality disorders (n = 1,393). Results revealed high internal consistency for almost all dimensions. The factor structure appeared almost identical as compared to the factor structure of the original DAPP-BQ, and was shown to be invariant across the various patient and community samples. Indices for convergent, discriminant and criterion related validity were satisfactory. It is concluded that the good psychometric characteristics of the original DAPP-BQ were preserved in the shortened version of the instrument

Grass silage for biorefinery – Separation efficiency and aerobic stability of silage and solid fraction

201