1,371 research outputs found
Evolving relativistic fluid spacetimes using pseudospectral methods and finite differencing
We present a new code for solving the coupled Einstein-hydrodynamics
equations to evolve relativistic, self-gravitating fluids. The Einstein field
equations are solved on one grid using pseudospectral methods, while the fluids
are evolved on another grid by finite differencing. We discuss implementation
details, such as the communication between the grids and the treatment of
stellar surfaces, and present code tests.Comment: To appear in the Proceedings of the Eleventh Marcel Grossmann Meetin
Radiation measurements from polar and geosynchronous satellites
The following topics are discussed: (1) cloud effects in climate determination; (2) annual variation in the global heat balance of the earth; (3) the accuracy of precipitation estimates made from passive microwave measurements from satellites; (4) seasonal oceanic precipitation frequencies; (5) determination of mesoscale temperature and moisture fields over land from satellite radiance measurements; and (6) Nimbus 6 scanning microwave spectrometer data evaluation for surface wind and pressure components in tropical storms
Orbiting binary black hole evolutions with a multipatch high order finite-difference approach
We present numerical simulations of orbiting black holes for around twelve
cycles, using a high-order multipatch approach. Unlike some other approaches,
the computational speed scales almost perfectly for thousands of processors.
Multipatch methods are an alternative to AMR (adaptive mesh refinement), with
benefits of simplicity and better scaling for improving the resolution in the
wave zone. The results presented here pave the way for multipatch evolutions of
black hole-neutron star and neutron star-neutron star binaries, where high
resolution grids are needed to resolve details of the matter flow
Black hole-neutron star mergers: effects of the orientation of the black hole spin
The spin of black holes in black hole-neutron star (BHNS) binaries can have a
strong influence on the merger dynamics and the postmerger state; a wide
variety of spin magnitudes and orientations are expected to occur in nature. In
this paper, we report the first simulations in full general relativity of BHNS
mergers with misaligned black hole spin. We vary the spin magnitude from a/m=0
to a/m=0.9 for aligned cases, and we vary the misalignment angle from 0 to 80
degrees for a/m=0.5. We restrict our study to 3:1 mass ratio systems and use a
simple Gamma-law equation of state. We find that the misalignment angle has a
strong effect on the mass of the postmerger accretion disk, but only for angles
greater than ~ 40 degrees. Although the disk mass varies significantly with
spin magnitude and misalignment angle, we find that all disks have very similar
lifetimes ~ 100ms. Their thermal and rotational profiles are also very similar.
For a misaligned merger, the disk is tilted with respect to the final black
hole's spin axis. This will cause the disk to precess, but on a timescale
longer than the accretion time. In all cases, we find promising setups for
gamma-ray burst production: the disks are hot, thick, and hyperaccreting, and a
baryon-clear region exists above the black hole.Comment: 15 pages, 13 figure
Black hole evolution by spectral methods
Current methods of evolving a spacetime containing one or more black holes
are plagued by instabilities that prohibit long-term evolution. Some of these
instabilities may be due to the numerical method used, traditionally finite
differencing. In this paper, we explore the use of a pseudospectral collocation
(PSC) method for the evolution of a spherically symmetric black hole spacetime
in one dimension using a hyperbolic formulation of Einstein's equations. We
demonstrate that our PSC method is able to evolve a spherically symmetric black
hole spacetime forever without enforcing constraints, even if we add dynamics
via a Klein-Gordon scalar field. We find that, in contrast to
finite-differencing methods, black hole excision is a trivial operation using
PSC applied to a hyperbolic formulation of Einstein's equations. We discuss the
extension of this method to three spatial dimensions.Comment: 20 pages, 17 figures, submitted to PR
Black hole-neutron star mergers for 10 solar mass black holes
General relativistic simulations of black hole-neutron star mergers have
currently been limited to low-mass black holes (less than 7 solar mass), even
though population synthesis models indicate that a majority of mergers might
involve more massive black holes (10 solar mass or more). We present the first
general relativistic simulations of black hole-neutron star mergers with 10
solar mass black holes. For massive black holes, the tidal forces acting on the
neutron star are usually too weak to disrupt the star before it reaches the
innermost stable circular orbit of the black hole. Varying the spin of the
black hole in the range a/M = 0.5-0.9, we find that mergers result in the
disruption of the star and the formation of a massive accretion disk only for
large spins a/M>0.7-0.9. From these results, we obtain updated constraints on
the ability of BHNS mergers to be the progenitors of short gamma-ray bursts as
a function of the mass and spin of the black hole. We also discuss the
dependence of the gravitational wave signal on the black hole parameters, and
provide waveforms and spectra from simulations beginning 7-8 orbits before
merger.Comment: 11 pages, 11 figures - Updated to match published versio
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 , 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 : . 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 for our equal mass case
Unequal mass binary neutron star simulations with neutrino transport: Ejecta and neutrino emission
We present 12 new simulations of unequal mass neutron star mergers. The simulations are performed with the SpEC code, and utilize nuclear-theory-based equations of state and a two-moment gray neutrino transport scheme with an improved energy estimate based on evolving the number density. We model the neutron stars with the SFHo, LS220, and DD2 equations of state (EOS) and we study the neutrino and matter emission of all 12 models to search for robust trends between binary parameters and emission characteristics. We find that the total mass of the dynamical ejecta exceeds 0.01  M⊙ only for SFHo with weak dependence on the mass ratio across all models. We find that the ejecta have a broad electron fraction (Y_e) distribution (≈0.06–0.48), with mean 0.2. Y_e increases with neutrino irradiation over time, but decreases with increasing binary asymmetry. We also find that the models have ejecta with a broad asymptotic velocity distribution (≈0.05–0.7c). The average velocity lies in the range 0.2c−0.3c and decreases with binary asymmetry. Furthermore, we find that disk mass increases with binary asymmetry and stiffness of the EOS. The Y_e of the disk increases with softness of the EOS. The strongest neutrino emission occurs for the models with soft EOS. For (anti) electron neutrinos we find no significant dependence of the magnitude or angular distribution or neutrino luminosity with mass ratio. The heavier neutrino species have a luminosity dependence on mass ratio but an angular distribution which does not change with mass ratio
Magnetic effects on the low-T/|W| instability in differentially rotating neutron stars
Dynamical instabilities in protoneutron stars may produce gravitational waves
whose observation could shed light on the physics of core-collapse supernovae.
When born with sufficient differential rotation, these stars are susceptible to
a shear instability (the "low-T/|W| instability"), but such rotation can also
amplify magnetic fields to strengths where they have a considerable impact on
the dynamics of the stellar matter. Using a new magnetohydrodynamics module for
the Spectral Einstein Code, we have simulated a differentially-rotating neutron
star in full 3D to study the effects of magnetic fields on this instability.
Though strong toroidal fields were predicted to suppress the low-T/|W|
instability, we find that they do so only in a small range of field strengths.
Below 4e13 G, poloidal seed fields do not wind up fast enough to have an effect
before the instability saturates, while above 5e14 G, magnetic instabilities
can actually amplify a global quadrupole mode (this threshold may be even lower
in reality, as small-scale magnetic instabilities remain difficult to resolve
numerically). Thus, the prospects for observing gravitational waves from such
systems are not in fact diminished over most of the magnetic parameter space.
Additionally, we report that the detailed development of the low-T/|W|
instability, including its growth rate, depends strongly on the particular
numerical methods used. The high-order methods we employ suggest that growth
might be considerably slower than found in some previous simulations.Comment: REVTeX 4.1, 21 pages, 18 figures, submitting to Physical Review
- …