Axisymmetric core collapse simulations using characteristic numerical relativity

We present results from axisymmetric stellar core collapse simulations in general relativity. Our hydrodynamics code has proved robust and accurate enough to allow for a detailed analysis of the global dynamics of the collapse. Contrary to traditional approaches based on the 3+1 formulation of the gravitational field equations, our framework uses a foliation based on a family of outgoing light cones, emanating from a regular center, and terminating at future null infinity. Such a coordinate system is well adapted to the study of interesting dynamical spacetimes in relativistic astrophysics such as stellar core collapse and neutron star formation. Perhaps most importantly this procedure allows for the unambiguous extraction of gravitational waves at future null infinity without any approximation, along with the commonly used quadrupole formalism for the gravitational wave extraction. Our results concerning the gravitational wave signals show noticeable disagreement when those are extracted by computing the Bondi news at future null infinity on the one hand and by using the quadrupole formula on the other hand. We have strong indication that for our setup the quadrupole formula on the null cone does not lead to physical gravitational wave signals. The Bondi gravitational wave signals extracted at infinity show typical oscillation frequencies of about 0.5 kHz.Comment: 17 pages, 18 figures, submitted to Phys. Rev.

First-order quasilinear canonical representation of the characteristic formulation of the Einstein equations

We prescribe a choice of 18 variables in all that casts the equations of the fully nonlinear characteristic formulation of general relativity in first--order quasi-linear canonical form. At the analytical level, a formulation of this type allows us to make concrete statements about existence of solutions. In addition, it offers concrete advantages for numerical applications as it now becomes possible to incorporate advanced numerical techniques for first order systems, which had thus far not been applicable to the characteristic problem of the Einstein equations, as well as in providing a framework for a unified treatment of the vacuum and matter problems. This is of relevance to the accurate simulation of gravitational waves emitted in astrophysical scenarios such as stellar core collapse.Comment: revtex4, 7 pages, text and references added, typos corrected, to appear in Phys. Rev.

Tethered subsatellite study

The results are presented of studies performed relating to the feasibility of deploying a subsatellite from the shuttle by means of a tether. The dynamics, the control laws, the aerodynamics, the heating, and some communication considerations of the tethered subsatellite system are considered. Nothing was found that prohibits the use of a subsatellite joined to the shuttle by a long (100 km) tether. More detailed studies directed at specific applications are recommended

Bondian frames to couple matter with radiation

A study is presented for the non linear evolution of a self gravitating distribution of matter coupled to a massless scalar field. The characteristic formulation for numerical relativity is used to follow the evolution by a sequence of light cones open to the future. Bondian frames are used to endow physical meaning to the matter variables and to the massless scalar field. Asymptotic approaches to the origin and to infinity are achieved; at the boundary surface interior and exterior solutions are matched guaranteeing the Darmois--Lichnerowicz conditions. To show how the scheme works some numerical models are discussed. We exemplify evolving scalar waves on the following fixed backgrounds: A) an atmosphere between the boundary surface of an incompressible mixtured fluid and infinity; B) a polytropic distribution matched to a Schwarzschild exterior; C) a Schwarzschild- Schwarzschild spacetime. The conservation of energy, the Newman--Penrose constant preservation and other expected features are observed.Comment: 20 pages, 6 figures; to appear in General Relativity and Gravitatio

Numerical Evolution of axisymmetric vacuum spacetimes: a code based on the Galerkin method

We present the first numerical code based on the Galerkin and Collocation methods to integrate the field equations of the Bondi problem. The Galerkin method like all spectral methods provide high accuracy with moderate computational effort. Several numerical tests were performed to verify the issues of convergence, stability and accuracy with promising results. This code opens up several possibilities of applications in more general scenarios for studying the evolution of spacetimes with gravitational waves.Comment: 11 pages, 6 figures. To appear in Classical and Quantum Gravit

Scalar field induced oscillations of neutron stars and gravitational collapse

We study the interaction of massless scalar fields with self-gravitating neutron stars by means of fully dynamic numerical simulations of the Einstein-Klein-Gordon perfect fluid system. Our investigation is restricted to spherical symmetry and the neutron stars are approximated by relativistic polytropes. Studying the nonlinear dynamics of isolated neutron stars is very effectively performed within the characteristic formulation of general relativity, in which the spacetime is foliated by a family of outgoing light cones. We are able to compactify the entire spacetime on a computational grid and simultaneously impose natural radiative boundary conditions and extract accurate radiative signals. We study the transfer of energy from the scalar field to the fluid star. We find, in particular, that depending on the compactness of the neutron star model, the scalar wave forces the neutron star either to oscillate in its radial modes of pulsation or to undergo gravitational collapse to a black hole on a dynamical timescale. The radiative signal, read off at future null infinity, shows quasi-normal oscillations before the setting of a late time power-law tail.Comment: 12 pages, 13 figures, submitted to Phys. Rev.

Gravitational waves from axisymmetrically oscillating neutron stars in general relativistic simulations

Gravitational waves from oscillating neutron stars in axial symmetry are studied performing numerical simulations in full general relativity. Neutron stars are modeled by a polytropic equation of state for simplicity. A gauge-invariant wave extraction method as well as a quadrupole formula are adopted for computation of gravitational waves. It is found that the gauge-invariant variables systematically contain numerical errors generated near the outer boundaries in the present axisymmetric computation. We clarify their origin, and illustrate it possible to eliminate the dominant part of the systematic errors. The best corrected waveforms for oscillating and rotating stars currently contain errors of magnitude $\sim 10^{-3}$ in the local wave zone. Comparing the waveforms obtained by the gauge-invariant technique with those by the quadrupole formula, it is shown that the quadrupole formula yields approximate gravitational waveforms besides a systematic underestimation of the amplitude of $O(M/R)$ where $M$ and $R$ denote the mass and the radius of neutron stars. However, the wave phase and modulation of the amplitude can be computed accurately. This indicates that the quadrupole formula is a useful tool for studying gravitational waves from rotating stellar core collapse to a neutron star in fully general relativistic simulations. Properties of the gravitational waveforms from the oscillating and rigidly rotating neutron stars are also addressed paying attention to the oscillation associated with fundamental modes

Characteristic Evolution and Matching

I review the development of numerical evolution codes for general relativity based upon the characteristic initial value problem. Progress in characteristic evolution is traced from the early stage of 1D feasibility studies to 2D axisymmetric codes that accurately simulate the oscillations and gravitational collapse of relativistic stars and to current 3D codes that provide pieces of a binary black hole spacetime. Cauchy codes have now been successful at simulating all aspects of the binary black hole problem inside an artificially constructed outer boundary. A prime application of characteristic evolution is to extend such simulations to null infinity where the waveform from the binary inspiral and merger can be unambiguously computed. This has now been accomplished by Cauchy-characteristic extraction, where data for the characteristic evolution is supplied by Cauchy data on an extraction worldtube inside the artificial outer boundary. The ultimate application of characteristic evolution is to eliminate the role of this outer boundary by constructing a global solution via Cauchy-characteristic matching. Progress in this direction is discussed.Comment: New version to appear in Living Reviews 2012. arXiv admin note: updated version of arXiv:gr-qc/050809

A remarkable record of the genus Pseudolucia from Bolivia (Lepidoptera: Lycaenidae)

The occurrence of a taxon morphologically close to Pseudolucia jujuyensis Bálint, Eisele & Johnson, 2000 is recorded in dry habitats of Torotoro Dinosaurs National Park, Potosí, Bolivia. This record remarkably extends the range of Pseudolucia by almost 800 km northwards in austral South America. Five specimens were available for examinations, hence wing-pattern, genitalia and mitochondrial DNA were analysed. However, the taxonomy of the specimens could not be satisfactorily resolved in relation to P. jujuyensis, for which only the holotype exists. The females use Cuscuta for ovipositing, what is supposedly the larval host – a remarkable character of the chilensis species group of Pseudolucia, which includes P. jujuyensis. On the basis of molecular markers it was revealed that the Torotoro population is the sister to the rest of the chilensis species group, which together are the clade sister to the rest of the genus

Understanding Hematopoietic Stem Cell Development through Functional Correlation of Their Proliferative Status with the Intra-aortic Cluster Architecture

During development, hematopoietic stem cells (HSCs) emerge in the aorta-gonad-mesonephros (AGM) region through a process of multi-step maturation and expansion. While proliferation of adult HSCs is implicated in the balance between self-renewal and differentiation, very little is known about the proliferation status of nascent HSCs in the AGM region. Using Fucci reporter mice that enable in vivo visualization of cell-cycle status, we detect increased proliferation during pre-HSC expansion followed by a slowing down of cycling once cells start to acquire a definitive HSC state, similar to fetal liver HSCs. We observe time-specific changes in intra-aortic hematopoietic clusters corresponding to HSC maturation stages. The proliferative architecture of the clusters is maintained in an orderly anatomical manner with slowly cycling cells at the base and more actively proliferating cells at the more apical part of the cluster, which correlates with c-KIT expression levels, thus providing an anatomical basis for the role of SCF in HSC maturation