Local Invariants and Pairwise Entanglement in Symmetric Multi-qubit System

Pairwise entanglement properties of a symmetric multi-qubit system are analyzed through a complete set of two-qubit local invariants. Collective features of entanglement, such as spin squeezing, are expressed in terms of invariants and a classifcation scheme for pairwise entanglement is proposed. The invariant criteria given here are shown to be related to the recently proposed (Phys. Rev. Lett. 95, 120502 (2005)) generalized spin squeezing inequalities for pairwise entanglement in symmetric multi-qubit states.Comment: 9 pages, 2 figures, REVTEX, Replaced with a published versio

Phase Variation in the Pulse Profile of SMC X-1

We present the results of timing and spectral analysis of X-ray high state observations of the high-mass X-ray pulsar SMC X-1 with Chandra, XMM-Newton, and ROSAT, taken between 1991 and 2001. The source has L_X ~ 3-5 x 10^38 ergs/s, and the spectra can be modeled as a power law plus blackbody with kT_BB \~ 0.18 keV and reprocessed emission radius R_BB ~ 2 x 10^8 cm, assuming a distance of 60 kpc to the source. Energy-resolved pulse profiles show several distinct forms, more than half of which include a second pulse in the soft profile, previously documented only in hard energies. We also detect significant variation in the phase shift between hard and soft pulses, as has recently been reported in Her X-1. We suggest an explanation for the observed characteristics of the soft pulses in terms of precession of the accretion disk.Comment: 4 pages, 4 figures, accepted for publication in ApJL; v2 minor corrections, as will appear in ApJ

Relativistic MHD with Adaptive Mesh Refinement

This paper presents a new computer code to solve the general relativistic magnetohydrodynamics (GRMHD) equations using distributed parallel adaptive mesh refinement (AMR). The fluid equations are solved using a finite difference Convex ENO method (CENO) in 3+1 dimensions, and the AMR is Berger-Oliger. Hyperbolic divergence cleaning is used to control the $\nabla\cdot {\bf B}=0$ constraint. We present results from three flat space tests, and examine the accretion of a fluid onto a Schwarzschild black hole, reproducing the Michel solution. The AMR simulations substantially improve performance while reproducing the resolution equivalent unigrid simulation results. Finally, we discuss strong scaling results for parallel unigrid and AMR runs.Comment: 24 pages, 14 figures, 3 table

Recursive Encoding and Decoding of Noiseless Subsystem and Decoherence Free Subspace

When the environmental disturbace to a quantum system has a wavelength much larger than the system size, all qubits localized within a small area are under action of the same error operators. Noiseless subsystem and decoherence free subspace are known to correct such collective errors. We construct simple quantum circuits, which implement these collective error correction codes, for a small number $n$ of physical qubits. A single logical qubit is encoded with $n=3$ and $n=4$, while two logical qubits are encoded with $n=5$. The recursive relations among the subspaces employed in noiseless subsystem and decoherence free subspace play essential r\^oles in our implementation. The recursive relations also show that the number of gates required to encode $m$ logical qubits increases linearly in $m$.Comment: 9 pages, 3 figure

Suzaku X-ray Spectra and Pulse Profile Variations during the Superorbital Cycle of LMC X-4

We present results from spectral and temporal analyses of Suzaku and RXTE observations of the high mass X-ray binary LMC X-4. Using the full 13 years of available RXTE/ASM data, we apply the ANOVA and Lomb normalized Periodogram methods to obtain an improved superorbital period measurement of 30.32 +/- 0.04 days. The phase-averaged X-ray spectra from Suzaku observations during the high state of the superorbital period can be modeled in the 0.6--50 keV band as the combination of a power-law with Gamma ~ 0.6 and a high-energy cutoff at ~ 25 keV, a blackbody with kT_BB ~ 0.18 keV, and emission lines from Fe K_alpha, O VIII, and Ne IX (X Lyalpha). Assuming a distance of 50 kpc, The source has luminosity L_X ~ 3 x 10^38 ergs s^-1 in the 2--50 keV band, and the luminosity of the soft (blackbody) component is L_BB ~ 1.5 x 10^37 ergs s^-1. The energy resolved pulse profiles show single-peaked soft (0.5-1 keV) and hard (6-10 keV) pulses but a more complex pattern of medium (2-10 keV) pulses; cross-correlation of the hard with the soft pulses shows a phase shift that varies between observations. We interpret these results in terms of a picture in which a precessing disk reprocesses the hard X-rays and produces the observed soft spectral component, as has been suggested for the similar sources Her X-1 and SMC X-1.Comment: 13 emulateapj pages, 11 figures, 4 tables; accepted for publication in Ap

3D simulations of Einstein's equations: symmetric hyperbolicity, live gauges and dynamic control of the constraints

We present three-dimensional simulations of Einstein equations implementing a symmetric hyperbolic system of equations with dynamical lapse. The numerical implementation makes use of techniques that guarantee linear numerical stability for the associated initial-boundary value problem. The code is first tested with a gauge wave solution, where rather larger amplitudes and for significantly longer times are obtained with respect to other state of the art implementations. Additionally, by minimizing a suitably defined energy for the constraints in terms of free constraint-functions in the formulation one can dynamically single out preferred values of these functions for the problem at hand. We apply the technique to fully three-dimensional simulations of a stationary black hole spacetime with excision of the singularity, considerably extending the lifetime of the simulations.Comment: 21 pages. To appear in PR

Critical Phenomena in Neutron Stars I: Linearly Unstable Nonrotating Models

We consider the evolution in full general relativity of a family of linearly unstable isolated spherical neutron stars under the effects of very small, perturbations as induced by the truncation error. Using a simple ideal-fluid equation of state we find that this system exhibits a type-I critical behaviour, thus confirming the conclusions reached by Liebling et al. [1] for rotating magnetized stars. Exploiting the relative simplicity of our system, we are able carry out a more in-depth study providing solid evidences of the criticality of this phenomenon and also to give a simple interpretation of the putative critical solution as a spherical solution with the unstable mode being the fundamental F-mode. Hence for any choice of the polytropic constant, the critical solution will distinguish the set of subcritical models migrating to the stable branch of the models of equilibrium from the set of subcritical models collapsing to a black hole. Finally, we study how the dynamics changes when the numerically perturbation is replaced by a finite-size, resolution independent velocity perturbation and show that in such cases a nearly-critical solution can be changed into either a sub or supercritical. The work reported here also lays the basis for the analysis carried in a companion paper, where the critical behaviour in the the head-on collision of two neutron stars is instead considered [2].Comment: 15 pages, 9 figure

Collapse and black hole formation in magnetized, differentially rotating neutron stars

The capacity to model magnetohydrodynamical (MHD) flows in dynamical, strongly curved spacetimes significantly extends the reach of numerical relativity in addressing many problems at the forefront of theoretical astrophysics. We have developed and tested an evolution code for the coupled Einstein-Maxwell-MHD equations which combines a BSSN solver with a high resolution shock capturing scheme. As one application, we evolve magnetized, differentially rotating neutron stars under the influence of a small seed magnetic field. Of particular significance is the behavior found for hypermassive neutron stars (HMNSs), which have rest masses greater the mass limit allowed by uniform rotation for a given equation of state. The remnant of a binary neutron star merger is likely to be a HMNS. We find that magnetic braking and the magnetorotational instability lead to the collapse of HMNSs and the formation of rotating black holes surrounded by massive, hot accretion tori and collimated magnetic field lines. Such tori radiate strongly in neutrinos, and the resulting neutrino-antineutrino annihilation (possibly in concert with energy extraction by MHD effects) could provide enough energy to power short-hard gamma-ray bursts. To explore the range of outcomes, we also evolve differentially rotating neutron stars with lower masses and angular momenta than the HMNS models. Instead of collapsing, the non-hypermassive models form nearly uniformly rotating central objects which, in cases with significant angular momentum, are surrounded by massive tori.Comment: Submitted to a special issue of Classical and Quantum Gravity based around the New Frontiers in Numerical Relativity meeting at the Albert Einstein Institute, Potsdam, July 17-21, 200

Numerical stability of a new conformal-traceless 3+1 formulation of the Einstein equation

There is strong evidence indicating that the particular form used to recast the Einstein equation as a 3+1 set of evolution equations has a fundamental impact on the stability properties of numerical evolutions involving black holes and/or neutron stars. Presently, the longest lived evolutions have been obtained using a parametrized hyperbolic system developed by Kidder, Scheel and Teukolsky or a conformal-traceless system introduced by Baumgarte, Shapiro, Shibata and Nakamura. We present a new conformal-traceless system. While this new system has some elements in common with the Baumgarte-Shapiro-Shibata-Nakamura system, it differs in both the type of conformal transformations and how the non-linear terms involving the extrinsic curvature are handled. We show results from 3D numerical evolutions of a single, non-rotating black hole in which we demonstrate that this new system yields a significant improvement in the life-time of the simulations.Comment: 7 pages, 2 figure