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

The Globular Cluster System in the Inner Region of M87

1057 globular cluster candidates have been identified in a WFPC2 image of the inner region of M87. The Globular Cluster Luminosity Function (GCLF) can be well fit by a Gaussian profile with a mean value of m_V^0=23.67 +/- 0.07 mag and sigma=1.39 +/- 0.06 mag (compared to m_V^0=23.74 mag and sigma=1.44 mag from an earlier study using the same data by Whitmore it et al. 1995). The GCLF in five radial bins is found to be statistically the same at all points, showing no clear evidence of dynamical destruction processes based on the luminosity function (LF), in contradiction to the claim by Gnedin (1997). Similarly, there is no obvious correlation between the half light radius of the clusters and the galactocentric distance. The core radius of the globular cluster density distribution is R_c=56'', considerably larger than the core of the stellar component (R_c=6.8''). The mean color of the cluster candidates is V-I=1.09 mag which corresponds to an average metallicity of Fe/H = -0.74 dex. The color distribution is bimodal everywhere, with a blue peak at V-I=0.95 mag and a red peak at V-I=1.20 mag. The red population is only 0.1 magnitude bluer than the underlying galaxy, indicating that these clusters formed late in the metal enrichment history of the galaxy and were possibly created in a burst of star/cluster formation 3-6 Gyr after the blue population. We also find that both the red and the blue cluster distributions have a more elliptical shape (Hubble type E3.5) than the nearly spherical galaxy. The average half light radius of the clusters is ~2.5 pc which is comparable to the 3 pc average effective radius of the Milky Way clusters, though the red candidates are ~20% smaller than the blue ones.Comment: 40 pages, 17 figures, 4 tables, latex, accepted for publication in the Ap

Quantum entanglement between electronic and vibrational degrees of freedom in molecules

We consider the quantum entanglement of the electronic and vibrational degrees of freedom in molecules with a tendency towards double welled potentials using model coupled harmonic diabatic potential-energy surfaces. The von Neumann entropy of the reduced density matrix is used to quantify the electron-vibration entanglement for the lowest two vibronic wavefunctions in such a bipartite system. Significant entanglement is found only in the region in which the ground vibronic state contains a density profile that is bimodal (i.e., contains two separate local minima). However, in this region two distinct types of entanglement are found: (1) entanglement that arises purely from the degeneracy of energy levels in the two potential wells and which is destroyed by slight asymmetry, and (2) entanglement that involves strongly interacting states in each well that is relatively insensitive to asymmetry. These two distinct regions are termed fragile degeneracy-induced entanglement and persistent entanglement, respectively. Six classic molecular systems describable by two diabatic states are considered: ammonia, benzene, semibullvalene, pyridine excited triplet states, the Creutz-Taube ion, and the radical cation of the "special pair" of chlorophylls involved in photosynthesis. These chemically diverse systems are all treated using the same general formalism and the nature of the entanglement that they embody is elucidated

Operationally Invariant Measure of the Distance between Quantum States by Complementary Measurements

We propose an operational measure of distance of two quantum states, which conversely tells us their closeness. This is defined as a sum of differences in partial knowledge over a complete set of mutually complementary measurements for the two states. It is shown that the measure is operationally invariant and it is equivalent to the Hilbert-Schmidt distance. The operational measure of distance provides a remarkable interpretation of the information distance between quantum states.Comment: 4 page

A Persistent Disk Wind in GRS 1915+105 with NICER

The bright, erratic black hole X-ray binary GRS 1915+105 has long been a target for studies of disk instabilities, radio/infrared jets, and accretion disk winds, with implications that often apply to sources that do not exhibit its exotic X-ray variability. With the launch of NICER, we have a new opportunity to study the disk wind in GRS 1915+105 and its variability on short and long timescales. Here we present our analysis of 39 NICER observations of GRS 1915+105 collected during five months of the mission data validation and verification phase, focusing on Fe XXV and Fe XXVI absorption. We report the detection of strong Fe XXVI in 32 (>80%) of these observations, with another four marginal detections; Fe XXV is less common, but both likely arise in the well-known disk wind. We explore how the properties of this wind depends on broad characteristics of the X-ray lightcurve: mean count rate, hardness ratio, and fractional RMS variability. The trends with count rate and RMS are consistent with an average wind column density that is fairly steady between observations but varies rapidly with the source on timescales of seconds. The line dependence on spectral hardness echoes known behavior of disk winds in outbursts of Galactic black holes; these results clearly indicate that NICER is a powerful tool for studying black hole winds.Comment: Accepted for publication in ApJL. Comments welcom

Fifteen years of XMM-Newton and Chandra monitoring of Sgr A*: Evidence for a recent increase in the bright flaring rate

We present a study of the X-ray flaring activity of Sgr A* during all the 150 XMM-Newton and Chandra observations pointed at the Milky Way center over the last 15 years. This includes the latest XMM-Newton and Chandra campaigns devoted to monitoring the closest approach of the very red Br-Gamma emitting object called G2. The entire dataset analysed extends from September 1999 through November 2014. We employed a Bayesian block analysis to investigate any possible variations in the characteristics (frequency, energetics, peak intensity, duration) of the flaring events that Sgr A* has exhibited since their discovery in 2001. We observe that the total bright-or-very bright flare luminosity of Sgr A* increased between 2013-2014 by a factor of 2-3 (~3.5 sigma significance). We also observe an increase (~99.9% significance) from 0.27+-0.04 to 2.5+-1.0 day^-1 of the bright-or-very bright flaring rate of Sgr A*, starting in late summer 2014, which happens to be about six months after G2's peri-center passage. This might indicate that clustering is a general property of bright flares and that it is associated with a stationary noise process producing flares not uniformly distributed in time (similar to what is observed in other quiescent black holes). If so, the variation in flaring properties would be revealed only now because of the increased monitoring frequency. Alternatively, this may be the first sign of an excess accretion activity induced by the close passage of G2. More observations are necessary to distinguish between these two hypotheses.Comment: Accepted for publication in MNRA

Criticality and convergence in Newtonian collapse

We study through numerical simulation the spherical collapse of isothermal gas in Newtonian gravity. We observe a critical behavior which occurs at the threshold of gravitational instability leading to core formation. For a given initial density profile, we find a critical temperature, which is of the same order as the virial temperature of the initial configuration. For the exact critical temperature, the collapse converges to a self-similar form, the first member in Hunter's family of self-similar solutions. For a temperature close to the critical value, the collapse first approaches this critical solution. Later on, in the supercritical case, the collapse converges to another self-similar solution, which is called the Larson-Penston solution. In the subcritical case, the gas bounces and disperses to infinity. We find two scaling laws: one for the collapsed mass in the supercritical case and the other for the maximum density reached before dispersal in the subcritical case. The value of the critical exponent is measured to be $\simeq 0.11$ in the supercritical case, which agrees well with the predicted value $\simeq 0.10567$. These critical properties are quite similar to those observed in the collapse of a radiation fluid in general relativity. We study the response of the system to temperature fluctuation and discuss astrophysical implications for the insterstellar medium structure and for the star formation process. Newtonian critical behavior is important not only because it provides a simple model for general relativity but also because it is relevant for astrophysical systems such as molecular clouds.Comment: 15 pages, 8 figures, accepted for publication in PRD, figures 1 and 3 at lower resolution than in journal version, typos correcte

Coherent Ro-vibrational Revivals in a Thermal Molecular Ensemble

We report an experimental and theoretical study of the evolution of vibrational coherence in a thermal ensemble of nitrogen molecules. Rotational dephasing and rephasing of the vibrational coherence is detected by coherent anti-Stokes Raman scattering. The existence of ro-vibrational coupling and the discrete energy spectrum of the rotational bath lead to a whole new class of full and fractional ro-vibrational revivals. Following the rich ro-vibrational dynamics on a nanosecond time scale with sub-picosecond time resolution enables us to determine the second-order ro-vibrational constant $gamma_e$ and assess new possibilities of controlling decoherence.Comment: submitted at Physical Review