16,069 research outputs found

    Horizon energy and angular momentum from a Hamiltonian perspective

    Full text link
    Classical black holes and event horizons are highly non-local objects, defined in terms of the causal past of future null infinity. Alternative, (quasi)local definitions are often used in mathematical, quantum, and numerical relativity. These include apparent, trapping, isolated, and dynamical horizons, all of which are closely associated to two-surfaces of zero outward null expansion. In this paper we show that three-surfaces which can be foliated with such two-surfaces are suitable boundaries in both a quasilocal action and a phase space formulation of general relativity. The resulting formalism provides expressions for the quasilocal energy and angular momentum associated with the horizon. The values of the energy and angular momentum are in agreement with those derived from the isolated and dynamical horizon frameworks.Comment: 39 pages, 3 figures, Final Version : content essentially unchanged but many small improvements made in response to referees, a few references adde

    Comparative Monte Carlo Efficiency by Monte Carlo Analysis

    Full text link
    We propose a modified power method for computing the subdominant eigenvalue λ2\lambda_2 of a matrix or continuous operator. Here we focus on defining simple Monte Carlo methods for its application. The methods presented use random walkers of mixed signs to represent the subdominant eigenfuction. Accordingly, the methods must cancel these signs properly in order to sample this eigenfunction faithfully. We present a simple procedure to solve this sign problem and then test our Monte Carlo methods by computing the λ2\lambda_2 of various Markov chain transition matrices. We first computed λ2{\lambda_2} for several one and two dimensional Ising models, which have a discrete phase space, and compared the relative efficiencies of the Metropolis and heat-bath algorithms as a function of temperature and applied magnetic field. Next, we computed λ2\lambda_2 for a model of an interacting gas trapped by a harmonic potential, which has a mutidimensional continuous phase space, and studied the efficiency of the Metropolis algorithm as a function of temperature and the maximum allowable step size Δ\Delta. Based on the λ2\lambda_2 criterion, we found for the Ising models that small lattices appear to give an adequate picture of comparative efficiency and that the heat-bath algorithm is more efficient than the Metropolis algorithm only at low temperatures where both algorithms are inefficient. For the harmonic trap problem, we found that the traditional rule-of-thumb of adjusting Δ\Delta so the Metropolis acceptance rate is around 50% range is often sub-optimal. In general, as a function of temperature or Δ\Delta, λ2\lambda_2 for this model displayed trends defining optimal efficiency that the acceptance ratio does not. The cases studied also suggested that Monte Carlo simulations for a continuum model are likely more efficient than those for a discretized version of the model.Comment: 23 pages, 8 figure

    Modelling the circumstellar medium in RS Ophiuchi and its link to Type Ia supernovae

    Get PDF
    Recent interpretations of narrow, variable absorption lines detected in some Type Ia supernovae suggest that their progenitors are surrounded by dense, circumstellar material. Similar variations detected in the symbiotic recurrent nova system RS Oph, which undergoes thermonuclear outbursts every ~20 years, making it an ideal candidate to investigate the origin of these lines. To this end, we present simulations of multiple mass transfer-nova cycles in RS Oph. We find that the quiescent mass transfer produces a dense, equatorial outflow, i.e., concentrated towards the binary orbital plane, and an accretion disc forms around the white dwarf. The interaction of a spherical nova outburst with these aspherical circumstellar structures produces a bipolar outflow, similar to that seen in HST imaging of the 2006 outburst. In order to produce an ionization structure that is consistent with observations, a mass-loss rate of 5×10−7 M⊙ yr−15 \times 10^{-7}\,\mathrm{M}_{\odot}\,\mathrm{yr}^{-1} from the red giant is required. The simulations also produce a polar accretion flow, which may explain the broad wings of the quiescent H {\alpha} line and hard X-rays. By comparing simulated absorption line profiles to observations of the 2006 outburst, we are able to determine which components arise in the wind and which are due to the novae. We explore the possible behaviour of absorption line profiles as they may appear should a supernova occur in a system like RS Oph. Our models show similarities to supernovae like SN 2006X, but require a high mass-loss rate, M˙∌10−6\dot{M} \sim 10^{-6} to 10−5 M⊙ yr−110^{-5}\,\mathrm{M_\odot}\,\mathrm{yr}^{-1}, to explain the variability in SN 2006X
    • 

    corecore