General relativistic hydrodynamics in curvilinear coordinates

In this paper we report on what we believe is the first successful implementation of relativistic hydrodynamics, coupled to dynamical spacetimes, in spherical polar coordinates without symmetry assumptions. We employ a high-resolution shock-capturing scheme, which requires that the equations be cast in flux-conservative form. One example of such a form is the :Valencia" formulation, which has been adopted in numerous applications, in particular in Cartesian coordinates. Here we generalize this formulation to allow for a reference-metric approach, which provides a natural framework for calculations in curvilinear coordinates. In spherical polar coordinates, for example, it allows for an analytical treatment of the singular r and sin(\theta) terms that appear in the equations. We experiment with different versions of our generalized Valencia formulation in numerical implementations of relativistic hydrodynamics for both fixed and dynamical spacetimes. We consider a number of different tests -- non-rotating and rotating relativistic stars, as well as gravitational collapse to a black hole -- to demonstrate that our formulation provides a promising approach to performing fully relativistic astrophysics simulations in spherical polar coordinates.Comment: 14 pages, 8 figures, version to be published in PR

Numerical Relativity in Spherical Polar Coordinates: Off-center Simulations

We have recently presented a new approach for numerical relativity simulations in spherical polar coordinates, both for vacuum and for relativistic hydrodynamics. Our approach is based on a reference-metric formulation of the BSSN equations, a factoring of all tensor components, as well as a partially implicit Runge-Kutta method, and does not rely on a regularization of the equations, nor does it make any assumptions about the symmetry across the origin. In order to demonstrate this feature we present here several off-centered simulations, including simulations of single black holes and neutron stars whose center is placed away from the origin of the coordinate system, as well as the asymmetric head-on collision of two black holes. We also revisit our implementation of relativistic hydrodynamics and demonstrate that a reference-metric formulation of hydrodynamics together with a factoring of all tensor components avoids problems related to the coordinate singularities at the origin and on the axes. As a particularly demanding test we present results for a shock wave propagating through the origin of the spherical polar coordinate system.Comment: 13 pages, 11 figures; matches version published in PR

S-4B orbital workshop attitude control system study

Saturn S-4B orbital workshop attitude control system analysi

Fully covariant and conformal formulation of the Z4 system in a reference-metric approach: comparison with the BSSN formulation in spherical symmetry

We adopt a reference-metric approach to generalize a covariant and conformal version of the Z4 system of the Einstein equations. We refer to the resulting system as ``fully covariant and conformal", or fCCZ4 for short, since it is well suited for curvilinear as well as Cartesian coordinates. We implement this fCCZ4 formalism in spherical polar coordinates under the assumption of spherical symmetry using a partially-implicit Runge-Kutta (PIRK) method and show that our code can evolve both vacuum and non-vacuum spacetimes without encountering instabilities. Our method does not require regularization of the equations to handle coordinate singularities, nor does it depend on constraint-preserving outer boundary conditions, nor does it need any modifications of the equations for evolutions of black holes. We perform several tests and compare the performance of the fCCZ4 system, for different choices of certain free parameters, with that of BSSN. Confirming earlier results we find that, for an optimal choice of these parameters, and for neutron-star spacetimes, the violations of the Hamiltonian constraint can be between 1 and 3 orders of magnitude smaller in the fCCZ4 system than in the BSSN formulation. For black-hole spacetimes, on the other hand, any advantages of fCCZ4 over BSSN are less evident.Comment: 13 pages, 10 figure

From bcc to fcc: interplay between oscillating long-range and repulsive short-range forces

This paper supplements and partly extends an earlier publication, Phys. Rev. Lett. 95, 265501 (2005). In $d$-dimensional continuous space we describe the infinite volume ground state configurations (GSCs) of pair interactions \vfi and \vfi+\psi, where \vfi is the inverse Fourier transform of a nonnegative function vanishing outside the sphere of radius $K_0$, and $\psi$ is any nonnegative finite-range interaction of range $r_0\leq\gamma_d/K_0$, where $\gamma_3=\sqrt{6}\pi$. In three dimensions the decay of \vfi can be as slow as $\sim r^{-2}$, and an interaction of asymptotic form $\sim\cos(K_0r+\pi/2)/r^3$ is among the examples. At a dimension-dependent density $\rho_d$ the ground state of \vfi is a unique Bravais lattice, and for higher densities it is continuously degenerate: any union of Bravais lattices whose reciprocal lattice vectors are not shorter than $K_0$ is a GSC. Adding $\psi$ decreases the ground state degeneracy which, nonetheless, remains continuous in the open interval $(\rho_d,\rho_d')$, where $\rho_d'$ is the close-packing density of hard balls of diameter $r_0$. The ground state is unique at both ends of the interval. In three dimensions this unique GSC is the bcc lattice at $\rho_3$ and the fcc lattice at $\rho_3'=\sqrt{2}/r_0^3$.Comment: Published versio

Parent-only interventions in the treatment of childhood obesity: a systematic review of randomized controlled trials

Background An effective and cost-effective treatment is required for the treatment of childhood obesity. Comparing parent-only interventions with interventions including the child may help determine this. Methods A systematic review of published and ongoing studies until 2013, using electronic database and manual searches. Inclusion criteria: randomized controlled trials, overweight/obese children aged 5-12 years, parent-only intervention compared with an intervention that included the child, 6 months or more follow-up. Outcomes included measures of overweight. Results Ten papers from 6 completed studies, and 2 protocols for ongoing studies, were identified. Parent-only groups are either more effective than or similarly effective as child-only or parent-child interventions, in the change in degree of overweight. Most studies were at unclear risk of bias for randomization, allocation concealment and blinding of outcome assessors. Two trials were at high risk of bias for incomplete outcome data. Four studies showed higher dropout from parent-only interventions. One study examined programme costs and found parent-only interventions to be cheaper. Conclusions Parent-only interventions appear to be as effective as parent-child interventions in the treatment of childhood overweight/obesity, and may be less expensive. Reasons for higher attrition rates in parent-only interventions need further investigatio

Collapse of Nonlinear Gravitational Waves in Moving-Puncture Coordinates

We study numerical evolutions of nonlinear gravitational waves in moving-puncture coordinates. We adopt two different types of initial data -- Brill and Teukolsky waves -- and evolve them with two independent codes producing consistent results. We find that Brill data fail to produce long-term evolutions for common choices of coordinates and parameters, unless the initial amplitude is small, while Teukolsky wave initial data lead to stable evolutions, at least for amplitudes sufficiently far from criticality. The critical amplitude separates initial data whose evolutions leave behind flat space from those that lead to a black hole. For the latter we follow the interaction of the wave, the formation of a horizon, and the settling down into a time-independent trumpet geometry. We explore the differences between Brill and Teukolsky data and show that for less common choices of the parameters -- in particular negative amplitudes -- Brill data can be evolved with moving-puncture coordinates, and behave similarly to Teukolsky waves

Evolutionary bi-stability in pathogen transmission mode

Many pathogens transmit to new hosts by both infection (horizontal transmission) and transfer to the infected host's offspring (vertical transmission). These two transmission modes require speci®c adap- tations of the pathogen that can be mutually exclusive, resulting in a trade-off between horizontal and vertical transmission. We show that in mathematical models such trade-offs can lead to the simultaneous existence of two evolutionary stable states (evolutionary bi-stability) of allocation of resources to the two modes of transmission. We also show that jumping between evolutionary stable states can be induced by gradual environmental changes. Using quantitative PCR-based estimates of abundance in seed and vege- tative parts, we show that the pathogen of wheat, Phaeosphaeria nodorum, has jumped between two distinct states of transmission mode twice in the past 160 years, which, based on published evidence, we interpret as adaptation to environmental change. The ®nding of evolutionary bi-stability has impli- cations for human, animal and other plant diseases. An ill-judged change in a disease control programme could cause the pathogen to evolve a new, and possibly more damaging, combination of transmission modes. Similarly, environmental changes can shift the balance between transmission modes, with adverse effects on human, animal and plant health

Statistical-mechanical theory of the overall magnetic properties of mesocrystals

The mesocrystal showing both electrorheological and magnetorheological effects is called electro-magnetorheological (EMR) solids. Prediction of the overall magnetic properties of the EMR solids is a challenging task due to the coexistence of the uniaxially anisotropic behavior and structural transition as well as long-range interaction between the suspended particles. To consider the uniaxial anisotropy effect, we present an anisotropic Kirkwood-Fr\"{o}hlich equation for calculating the effective permeabilities by adopting an explicit characteristic spheroid rather than a characteristic sphere used in the derivation of the usual Kirkwood-Fr\"{o}hlich equation. Further, by applying an Ewald-Kornfeld formulation we are able to investigate the effective permeability by including the structural transition and long-range interaction explicitly. Our theory can reduce to the usual Kirkwood-Fr\"{o}hlich equation and Onsager equation naturally. To this end, the numerical simulation shows the validity of monitoring the structure of EMR solids by detecting their effective permeabilities.Comment: 14 pages, 1 figur