Resonance bifurcations of robust heteroclinic networks

Robust heteroclinic cycles are known to change stability in resonance bifurcations, which occur when an algebraic condition on the eigenvalues of the system is satisfied and which typically result in the creation or destruction of a long-period periodic orbit. Resonance bifurcations for heteroclinic networks are more complicated because different subcycles in the network can undergo resonance at different parameter values, but have, until now, not been systematically studied. In this article we present the first investigation of resonance bifurcations in heteroclinic networks. Specifically, we study two heteroclinic networks in $\R^4$ and consider the dynamics that occurs as various subcycles in each network change stability. The two cases are distinguished by whether or not one of the equilibria in the network has real or complex contracting eigenvalues. We construct two-dimensional Poincare return maps and use these to investigate the dynamics of trajectories near the network. At least one equilibrium solution in each network has a two-dimensional unstable manifold, and we use the technique developed in [18] to keep track of all trajectories within these manifolds. In the case with real eigenvalues, we show that the asymptotically stable network loses stability first when one of two distinguished cycles in the network goes through resonance and two or six periodic orbits appear. In the complex case, we show that an infinite number of stable and unstable periodic orbits are created at resonance, and these may coexist with a chaotic attractor. There is a further resonance, for which the eigenvalue combination is a property of the entire network, after which the periodic orbits which originated from the individual resonances may interact. We illustrate some of our results with a numerical example.Comment: 46 pages, 20 figures. Supplementary material (two animated gifs) can be found on http://www.maths.leeds.ac.uk/~alastair/papers/KPR_res_net_abs.htm

Testing the no-hair theorem with GW150914

We analyze gravitational-wave data from the first LIGO detection of a binary black-hole merger (GW150914) in search of the ringdown of the remnant black hole. Using observations beginning at the peak of the signal, we find evidence of the fundamental quasinormal mode and at least one overtone, both associated with the dominant angular mode ($\ell=m=2$), with $3.6\sigma$ confidence. A ringdown model including overtones allows us to measure the final mass and spin magnitude of the remnant exclusively from postinspiral data, obtaining an estimate in agreement with the values inferred from the full signal. The mass and spin values we measure from the ringdown agree with those obtained using solely the fundamental mode at a later time, but have smaller uncertainties. Agreement between the postinspiral measurements of mass and spin and those using the full waveform supports the hypothesis that the GW150914 merger produced a Kerr black hole, as predicted by general relativity, and provides a test of the no-hair theorem at the ${\sim}10\%$ level. An independent measurement of the frequency of the first overtone yields agreement with the no-hair hypothesis at the ${\sim 20}\%$ level. As the detector sensitivity improves and the detected population of black hole mergers grows, we can expect that using overtones will provide even stronger tests.Comment: v2: journal versio

Upper bounds on success probabilities in linear optics

We develop an abstract way of defining linear-optics networks designed to perform quantum information tasks such as quantum gates. We will be mainly concerned with the nonlinear sign shift gate, but it will become obvious that all other gates can be treated in a similar manner. The abstract scheme is extremely well suited for analytical as well as numerical investigations since it reduces the number of parameters for a general setting. With that we show numerically and partially analytically for a wide class of states that the success probability of generating a nonlinear sign shift gate does not exceed 1/4 which to our knowledge is the strongest bound to date.Comment: 8 pages, typeset using RevTex4, 5 EPS figure

The effect of symmetry breaking on the dynamics near a structurally stable heteroclinic cycle between equilibria and a periodic orbit

The effect of small forced symmetry breaking on the dynamics near a structurally stable heteroclinic cycle connecting two equilibria and a periodic orbit is investigated. This type of system is known to exhibit complicated, possibly chaotic dynamics including irregular switching of sign of various phase space variables, but details of the mechanisms underlying the complicated dynamics have not previously been investigated. We identify global bifurcations that induce the onset of chaotic dynamics and switching near a heteroclinic cycle of this type, and by construction and analysis of approximate return maps, locate the global bifurcations in parameter space. We find there is a threshold in the size of certain symmetry-breaking terms below which there can be no persistent switching. Our results are illustrated by a numerical example

On the equivalence of the Langevin and auxiliary field quantization methods for absorbing dielectrics

Recently two methods have been developed for the quantization of the electromagnetic field in general dispersing and absorbing linear dielectrics. The first is based upon the introduction of a quantum Langevin current in Maxwell's equations [T. Gruner and D.-G. Welsch, Phys. Rev. A 53, 1818 (1996); Ho Trung Dung, L. Kn\"{o}ll, and D.-G. Welsch, Phys. Rev. A 57, 3931 (1998); S. Scheel, L. Kn\"{o}ll, and D.-G. Welsch, Phys. Rev. A 58, 700 (1998)], whereas the second makes use of a set of auxiliary fields, followed by a canonical quantization procedure [A. Tip, Phys. Rev. A 57, 4818 (1998)]. We show that both approaches are equivalent.Comment: 7 pages, RevTeX, no figure

Treating instabilities in a hyperbolic formulation of Einstein's equations

We have recently constructed a numerical code that evolves a spherically symmetric spacetime using a hyperbolic formulation of Einstein's equations. For the case of a Schwarzschild black hole, this code works well at early times, but quickly becomes inaccurate on a time scale of 10-100 M, where M is the mass of the hole. We present an analytic method that facilitates the detection of instabilities. Using this method, we identify a term in the evolution equations that leads to a rapidly-growing mode in the solution. After eliminating this term from the evolution equations by means of algebraic constraints, we can achieve free evolution for times exceeding 10000M. We discuss the implications for three-dimensional simulations.Comment: 13 pages, 9 figures. To appear in Phys. Rev.

Black Hole Area in Brans-Dicke Theory

We have shown that the dynamics of the scalar field $\phi (x)= ``G^{-1}(x)"$ in Brans-Dicke theories of gravity makes the surface area of the black hole horizon {\it oscillatory} during its dynamical evolution. It explicitly explains why the area theorem does not hold in Brans-Dicke theory. However, we show that there exists a certain non-decreasing quantity defined on the event horizon which is proportional to the black hole entropy for the case of stationary solutions in Brans-Dicke theory. Some numerical simulations have been demonstrated for Oppenheimer-Snyder collapse in Brans-Dicke theory.Comment: 12 pages, latex, 5 figures, epsfig.sty, some statements clarified and two references added, to appear in Phys. Rev.

Black Hole--Scalar Field Interactions in Spherical Symmetry

We examine the interactions of a black hole with a massless scalar field using a coordinate system which extends ingoing Eddington-Finkelstein coordinates to dynamic spherically symmetric-spacetimes. We avoid problems with the singularity by excising the region of the black hole interior to the apparent horizon. We use a second-order finite difference scheme to solve the equations. The resulting program is stable and convergent and will run forever without problems. We are able to observe quasi-normal ringing and power-law tails as well an interesting nonlinear feature.Comment: 16 pages, 26 figures, RevTex, to appear in Phys. Rev.

Numerical simulations of neutron star-black hole binaries in the near-equal-mass regime

Simulations of neutron star-black hole (NSBH) binaries generally consider black holes with masses in the range $(5-10)M_\odot$, where we expect to find most stellar mass black holes. The existence of lower mass black holes, however, cannot be theoretically ruled out. Low-mass black holes in binary systems with a neutron star companion could mimic neutron star-neutron (NSNS) binaries, as they power similar gravitational wave (GW) and electromagnetic (EM) signals. To understand the differences and similarities between NSNS mergers and low-mass NSBH mergers, numerical simulations are required. Here, we perform a set of simulations of low-mass NSBH mergers, including systems compatible with GW170817. Our simulations use a composition and temperature dependent equation of state (DD2) and approximate neutrino transport, but no magnetic fields. We find that low-mass NSBH mergers produce remnant disks significantly less massive than previously expected, and consistent with the post-merger outflow mass inferred from GW170817 for moderately asymmetric mass ratio. The dynamical ejecta produced by systems compatible with GW170817 is negligible except if the mass ratio and black hole spin are at the edge of the allowed parameter space. That dynamical ejecta is cold, neutron-rich, and surprisingly slow for ejecta produced during the tidal disruption of a neutron star : $v\sim (0.1-0.15)c$. We also find that the final mass of the remnant black hole is consistent with existing analytical predictions, while the final spin of that black hole is noticeably larger than expected -- up to $\chi_{\rm BH}=0.84$ for our equal mass case