BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Peter Anninos

Peter Anninos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.IMastro-phastro-ph.GAastro-ph.SR

Catalog footprint

What is connected

17works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

General Relativistic Implicit Monte Carlo Radiation-Hydrodynamics

We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++: an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of the radiation-hydrodynamics equations, but generalized for both Newtonian and relativistic regimes. A multiple reference frame approach is used to geodesically transport photon packets (and solve the hydrodynamics equations) in the coordinate frame, while radiation-matter interactions are handled either in the fluid or electron frames then communicated via Lorentz boosts and orthonormal tetrad bases attached to the fluid. We describe a method for constructing estimators of radiation moments using path-weighting that generalizes to arbitrary coordinate systems in flat or curved spacetime. Absorption, emission, scattering, and relativistic Comptonization are among the matter interactions considered in this report. We discuss our formulations and numerical methods, and validate our models against a suite of radiation and coupled radiation-hydrodynamics test problems in both flat and curved spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Resolution Study of Thermonuclear Initiation in White Dwarf Tidal Disruption Events

We study the initiation of thermonuclear detonations in tidally disrupted white dwarf stars by intermediate-mass ($10^3 M_\odot$) black holes. The length scales required to resolve the initiation mechanism are not easily reached in three-dimensions, so instead we have devised two-dimensional proxy models which, together with a logarithmic gridding strategy, can adequately capture detonation wave fronts as material undergoes simultaneous compression and stretching from tidal forces. We consider 0.15 and 0.6 solar mass white dwarf stars parameterized by tidal strengths in the range $β=4~\text{to}~23$. High spatial resolution elucidates the manner and conditions leading to thermonuclear detonation, linking the initiation sequence to stellar composition and tidal strength. All of our models suffer sustained detonations triggered by a combination of adiabatic compression, mild thermonuclear preconditioning, and collisional heating, in degrees depending primarily on tidal strength. We find many diagnostics, such as temperature, total released energy, and iron group products, are fairly well-converged (better than 10%) at resolutions below 10 km along the scale height of the orbital plane. The exceptions are intermediate mass transients like calcium, which remain uncertain up to factors of two even at 1 km resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-Frequency General Relativistic Radiation-Hydrodynamics with $\bf{M}_1$ Closure

We report on recent upgrades to our general relativistic radiation-magnetohydrodynamics code, Cosmos++, which expands the two-moment, $\bf{M}_1$, radiation treatment from grey to multi-frequency transport, including Doppler and gravitational frequency shifts. The solver accommodates either photon (Bose-Einstein) or neutrino (Fermi-Dirac) statistical distribution functions with absorption, emission, and elastic scattering processes. An implicit scheme is implemented to simultaneously solve the primitive inversion problem together with the radiation-matter coupling source terms, providing stability over a broad range of opacities and optical depths where the interactions terms can be stiff. We discuss our formulations and numerical methods, and validate our methods against a wide variety of test problems spanning optically thin to thick regimes in flat, weakly curved, and strongly curved spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical Simulations of Optically Thick Accretion onto a Black Hole - II. Rotating Flow

In this paper we report on recent upgrades to our general relativistic radiation magnetohydrodynamics code, Cosmos++, including the development of a new primitive inversion scheme and a hybrid implicit-explicit solver with a more general closure relation for the radiation equations. The new hybrid solver helps stabilize the treatment of the radiation source terms, while the new closure allows for a much broader range of optical depths to be considered. These changes allow us to expand by orders of magnitude the range of temperatures, opacities, and mass accretion rates, and move a step closer toward our goal of performing global simulations of radiation-pressure-dominated black hole accretion disks. In this work we test and validate the new method against an array of problems. We also demonstrate its ability to handle super-Eddington, quasi-spherical accretion. Even with just a single proof-of-principle simulation, we already see tantalizing hints of the interesting phenomenology associated with the coupling of radiation and gas in super-Eddington accretion flows.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Three-Dimensional Moving-Mesh Simulations of Galactic Center Cloud G2

Using three-dimensional, moving-mesh simulations, we investigate the future evolution of the recently discovered gas cloud G2 traveling through the galactic center. We consider the case of a spherical cloud initially in pressure equilibrium with the background. Our suite of simulations explores the following parameters: the equation of state, radial profiles of the background gas, and start times for the evolution. Our primary focus is on how the fate of this cloud will affect the future activity of Sgr A*. From our simulations we expect an average feeding rate in the range of $5-19 \times 10^{-8}$ solar masses per year beginning in 2013 and lasting for at least 7 years (our simulations stop in year 2020). The accretion varies by less than a factor of three on timescales <1 month, and shows no more than a factor of 10 difference between the maximum and minimum. These rates are comparable to the current estimated accretion rate in the immediate vicinity of Sgr A*, although they represent only a small (<5%) increase over the current expected feeding rate at the effective inner boundary of our simulations. Therefore, the break up of cloud G2 may have only a minimal effect on the brightness and variability of Sgr A* over the next decade. This is because current models of the galactic center predict that most of the gas will be caught up in outflows. However, if the accreted G2 material can remain cold, it may not mix well with the hot, diffuse background gas, and instead accrete efficiently onto Sgr A*. Further observations of G2 will give us an unprecedented opportunity to test this idea. The break up of the cloud itself may also be observable. By tracking the amount of cloud energy that is dissipated during our simulations, we are able to get a rough estimate of the luminosity associated with its tidal disruption; we find values of a few $10^{36}$ erg/s.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Update on 3D Moving Mesh Simulations of Galactic Center Cloud G2

Using three-dimensional, moving-mesh simulations, we investigate the future evolution of the recently discovered gas cloud G2 traveling through the galactic center. From our simulations we expect an average feeding rate onto Sgr A* in the range of $(5-19) \times 10^{-8} M_\odot\mathrm{~yr}^{-1}$ beginning in 2014. This accretion varies by less than a factor of three on timescales of about 1 month, and shows no more than a factor of 10 difference between the maximum and minimum observed rates within any given model. These rates are comparable to the current estimated accretion rate in the immediate vicinity of Sgr A*, although they represent only a small (<10%) increase over the current expected feeding rate at the effective inner boundary of our simulations $(r_\mathrm{acc} = 750 R_S \sim 10^{15} \mathrm{cm})$. We also explore multiple possible equations of state to describe the gas. In examining the Br-$γ$ light curves produced from our simulations, we find that all of our isothermal models predict significant (factor of 10) enhancements in the luminosity of G2 as it approaches pericenter, in conflict with observations. Models that instead allow the cloud to heat as it is compressed do better at matching observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical Simulations of Optically Thick Accretion onto a Black Hole - I. Spherical Case

Modeling the radiation generated by accreting matter is an important step towards realistic simulations of black hole accretion disks, especially at high accretion rates. To this end, we have recently added radiation transport to the existing general relativistic magnetohydrodynamic code, Cosmos++. However, before attempting to model radiative accretion disks, we have tested the new code using a series of shock tube and Bondi (spherical inflow) problems. The four radiative shock tube tests, first presented by Farris et al. (2008), have known analytic solutions, allowing us to calculate errors and convergence rates for our code. The Bondi problem only has an analytic solution when radiative processes are ignored, but is pertinent because it is closer to the physics we ultimately want to study. In our simulations, we include Thomson scattering and thermal bremsstrahlung in the opacity, focusing exclusively on the super-Eddington regime. Unlike accretion onto bodies with solid surfaces, super-Eddington accretion onto black holes does not produce super-Eddington luminosity. In our examples, despite accreting at up to 300 times the Eddington rate, our measured luminosity is always several orders of magnitude below Eddington.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamical Bar-Mode Instability in Differentially Rotating Magnetized Neutron Stars

This paper presents a numerical study over a wide parameter space of the likelihood of the dynamical bar-mode instability in differentially rotating magnetized neutron stars. The innovative aspect of this study is the incorporation of magnetic fields in such a context, which have thus far been neglected in the purely hydrodynamical simulations available in the literature. The investigation uses the Cosmos++ code which allows us to perform three dimensional simulations on a cylindrical grid at high resolution. A sample of Newtonian magneto-hydrodynamical simulations starting from a set of models previously analyzed by other authors without magnetic fields has been performed, providing estimates of the effects of magnetic fields on the dynamical bar-mode deformation of rotating neutron stars. Overall, our results suggest that the effect of magnetic fields are not likely to be very significant in realistic configurations. Only in the most extreme cases are the magnetic fields able to suppress growth of the bar mode.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Application of the Cubed-Sphere Grid to Tilted Black-Hole Accretion Disks

In recent work we presented the first results of global general relativistic magnetohydrodynamic (GRMHD) simulations of tilted (or misaligned) accretion disks around rotating black holes. The simulated tilted disks showed dramatic differences from comparable untilted disks, such as asymmetrical accretion onto the hole through opposing "plunging streams" and global precession of the disk powered by a torque provided by the black hole. However, those simulations used a traditional spherical-polar grid that was purposefully underresolved along the pole, which prevented us from assessing the behavior of any jets that may have been associated with the tilted disks. To address this shortcoming we have added a block-structured "cubed-sphere" grid option to the Cosmos++ GRMHD code, which will allow us to simultaneously resolve the disk and polar regions. Here we present our implementation of this grid and the results of a small suite of validation tests intended to demonstrate that the new grid performs as expected. The most important test in this work is a comparison of identical tilted disks, one evolved using our spherical-polar grid and the other with the cubed-sphere grid. We also demonstrate an interesting dependence of the early-time evolution of our disks on their orientation with respect to the grid alignment. This dependence arises from the differing treatment of current sheets within the disks, especially whether they are aligned with symmetry planes of the grid or not.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

3D Relativistic MHD Simulation of a Tilted Accretion Disk Around a Rapidly Rotating Black Hole

We posit that accreting compact objects, including stellar mass black holes and neutron stars as well as supermassive black holes, may undergo extended periods of accretion during which the angular momentum of the disk at large scales is misaligned with that of the compact object. In such a scenario, Lense-Thirring precession caused by the rotating compact object can dramatically affect the disk. In this presentation we describe results from a three-dimensional relativistic magnetohydrodynamic simulation of an MRI turbulent disk accreting onto a tilted rapidly rotating black hole. For this case, the disk does not achieve the commonly described Bardeen-Petterson configuration; rather, it remains nearly planar, undergoing a slow global precession. Accretion from the disk onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is a consequence of the non-sphericity of the gravitational spacetime of the black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Computational Cosmology: from the Early Universe to the Large Scale Structure

In order to account for the observable Universe, any comprehensive theory or model of cosmology must draw from many disciplines of physics, including gauge theories of strong and weak interactions, the hydrodynamics and microphysics of baryonic matter, electromagnetic fields, and spacetime curvature, for example. Although it is difficult to incorporate all these physical elements into a single complete model of our Universe, advances in computing methods and technologies have contributed significantly towards our understanding of cosmological models, the Universe, and astrophysical processes within them. A sample of numerical calculations (and numerical methods) applied to specific issues in cosmology are reviewed in this article: from the Big Bang singularity dynamics to the fundamental interactions of gravitational waves; from the quark-hadron phase transition to the large scale structure of the Universe. The emphasis, although not exclusively, is on those calculations designed to test different models of cosmology against the observed Universe.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

A Phase Space Approach to the Gravitational Arrow of Time

We attempt to find a function that characterizes gravitational clumping and that increases monotonically as inhomogeneity increases. We choose $S = lnΩ$ as the candidate ``gravitational entropy'' function, where $Ω$ is the phase-space volume below the Hamiltonian H of the system under consideration. We compute $Ω$ for transverse electromagnetic waves and for gravitational wave, radiation and density perturbations in an expanding FLRW universe. These calculations are carried out in the linear regime under the assumption that the phases of the oscillators comprising the system are random. Entropy is thus attributed to the lack of knowledge of the exact field configuration. We find that $Ω$, and hence $lnΩ$ behaves as required. We also carry out calculations for Bianchi IX cosmological models and find that, even in this homogeneous case, the function can be interpreted sensibly. We compare our results with Penrose's C^2 hypothesis. Because S is defined to resemble the fundamental statistical mechanics definition of entropy, we are able to recover the entropy in a variety of familiar circumstances including, evidently, black-hole entropy. The results point to the utility of the relativistic ADM Hamiltonian formalism in establishing a connection between general relativity and statistical mechanics, although fully nonlinear calculations will need to be performed to remove any doubt.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Dynamics of Gravitational Waves in 3D: Formulations, Methods, and Tests

The dynamics of gravitational waves is investigated in full 3+1 dimensional numerical relativity, emphasizing the difficulties that one might encounter in numerical evolutions, particularly those arising from non-linearities and gauge degrees of freedom. Using gravitational waves with amplitudes low enough that one has a good understanding of the physics involved, but large enough to enable non-linear effects to emerge, we study the coupling between numerical errors, coordinate effects, and the nonlinearities of the theory. We discuss the various strategies used in identifying specific features of the evolution. We show the importance of the flexibility of being able to use different numerical schemes, different slicing conditions, different formulations of the Einstein equations (standard ADM vs. first order hyperbolic), and different sets of equations (linearized vs. full Einstein equations). A non-linear scalar field equation is presented which captures some properties of the full Einstein equations, and has been useful in our understanding of the coupling between finite differencing errors and non-linearites. We present a set of monitoring devices which have been crucial in our studying of the waves, including Riemann invariants, pseudo-energy momentum tensor, hamiltonian constraint violation, and fourier spectrum analysis.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

The Near-Linear Regime of Gravitational Waves in Numerical Relativity

We report on a systematic study of the dynamics of gravitational waves in full 3D numerical relativity. We find that there exists an interesting regime in the parameter space of the wave configurations: a near-linear regime in which the amplitude of the wave is low enough that one expects the geometric deviation from flat spacetime to be negligible, but nevertheless where nonlinearities can excite unstable modes of the Einstein evolution equations causing the metric functions to evolve out of control. The implications of this for numerical relativity are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Three dimensional numerical relativity: the evolution of black holes

We report on a new 3D numerical code designed to solve the Einstein equations for general vacuum spacetimes. This code is based on the standard 3+1 approach using cartesian coordinates. We discuss the numerical techniques used in developing this code, and its performance on massively parallel and vector supercomputers. As a test case, we present evolutions for the first 3D black hole spacetimes. We identify a number of difficulties in evolving 3D black holes and suggest approaches to overcome them. We show how special treatment of the conformal factor can lead to more accurate evolution, and discuss techniques we developed to handle black hole spacetimes in the absence of symmetries. Many different slicing conditions are tested, including geodesic, maximal, and various algebraic conditions on the lapse. With current resolutions, limited by computer memory sizes, we show that with certain lapse conditions we can evolve the black hole to about $t=50M$, where $M$ is the black hole mass. Comparisons are made with results obtained by evolving spherical initial black hole data sets with a 1D spherically symmetric code. We also demonstrate that an ``apparent horizon locking shift'' can be used to prevent the development of large gradients in the metric functions that result from singularity avoiding time slicings. We compute the mass of the apparent horizon in these spacetimes, and find that in many cases it can be conserved to within about 5\% throughout the evolution with our techniques and current resolution.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

The Head-On Collision of Two Equal Mass Black Holes: Numerical Methods

The head-on collision of two nonrotating axisymmetric equal mass black holes is treated numerically. We take as initial data the single parameter family of time-symmetric solutions discovered by Misner which consists of two Einstein-Rosen bridges that can be placed arbitrarily distant from one another. A number of problems associated with previous attempts to evolve these data sets have been overcome. In this article, we discuss our choices for coordinate systems, gauges, and the numerical algorithms that we have developed to evolve this system.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

The Collision of Two Black Holes

We study the head-on collision of two equal mass, nonrotating black holes. We consider a range of cases from holes surrounded by a common horizon to holes initially separated by about $20M$, where $M$ is the mass of each hole. We determine the waveforms and energies radiated for both the $\ell = 2$ and $\ell=4$ waves resulting from the collision. In all cases studied the normal modes of the final black hole dominate the spectrum. We also estimate analytically the total gravitational radiation emitted, taking into account the tidal heating of horizons using the membrane paradigm, and other effects. For the first time we are able to compare analytic calculations, black hole perturbation theory, and strong field, nonlinear numerical calculations for this problem, and we find excellent agreement.

0 citations0 reviewsNo code linkedOpen access