Source author record

Gabriel Rockefeller

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

11 published item(s)

preprint2015arXiv

Generating Optimal Initial Conditions for Smoothed Particle Hydrodynamics Simulations

We review existing SPH setup methods and outline their advantages, limitations and drawbacks. We present a new method for constructing initial conditions for smoothed particle hydrodynamics (SPH) simulations, which may also be of interest for N-body simulations, and demonstrate this method on a number of applications. This new method is inspired by adaptive binning techniques using weighted Voronoi tesselations. Particles are placed and iteratively moved based on their proximity to neighboring particles and the desired spatial resolution. This new method can satisfy arbitrarily complex spatial resolution requirements.

preprint2014arXiv

The Fallback Mechanisms in Core-Collapse Supernovae

Although the details of the core-collapse supernova mechanism are not fully understood, it is generally accepted that the energy released in the collapse produces a shock that disrupts the star and produces the explosion. Some of the stellar material does not receive enough energy to escape the potential well of the newly formed neutron star and it falls back on to the core. This fallback plays an important role in setting compact remnant masses and has been invoked to explain a number of observed phenomena: late-time neutrino emission, r-process element production, peculiar supernovae, and long-duration gamma-ray bursts. To examine the link between fallback and these observed phenomena, we must understand the nature of fallback and currently, multiple mechanisms have been proposed that lead to different characteristics of the fallback. In this paper, we model a series of three-dimensional explosions by following the fallback for two different massive progenitors to study both the nature of the fallback and the role of asymmetric explosions on this fallback. One of the most-cited mechanisms behind fallback argues that it is caused by the reverse shock produced when the outward moving supernova shock decelerates in the relatively flat density profile of the stellar envelope. We show that this mechanism plays a minor, if any, role in driving fallback and most supernova fallback is produced by prompt mechanisms.

preprint2013arXiv

First Simulations of Core- Collapse Supernovae to Supernova Remnants with SNSPH

We present the first 3-dimensional simulations following the evolution of supernova shocks from their inception in the stellar core through the development of a supernova remnant into the Sedov phase. Our set of simulations use two different progenitors and two different conditions for the structure of the circumstellar environment. These calculations demonstrate the role that supernova instabilities (the instabilities that develop as the shock drive through the star) play in defining the structure and long-term development of instabilities in supernova remnants. We also present a first investigation of the mixing between stellar and interstellar matter as the supernova evolves into a young supernova remnant.

preprint2012arXiv

A Case Study of Small Scale Structure Formation in 3D Supernova Simulations

It is suggested in observations of supernova remnants that a number of large- and small-scale structures form at various points in the explosion. Multidimensional modeling of core-collapse supernovae has been undertaken since SN1987A, and both simulations and observations suggest/show that Rayleigh-Taylor instabilities during the explosion is a main driver for the formation of structure in the remnants. We present a case study of structure formation in 3D in a \msol{15} supernova for different parameters. We investigate the effect of moderate asymmetries and different resolutions of the formation and morphology of the RT unstable region, and take first steps at determining typical physical quantities (size, composition) of arising clumps. We find that in this progenitor the major RT unstable region develops at the He/OC interface for all cases considered. The RT instabilities result in clumps that are overdense by 1-2 orders of magnitude with respect to the ambient gas, have size scales on the level of a few % of the remnant diameter, and are not diffused after the first $\sim30$ yrs of the remnant evolution, in the absence of a surrounding medium.

preprint2011arXiv

Convective-reactive proton-C12 combustion in Sakurai's object (V4334 Sagittarii) and implications for the evolution and yields from the first generations of stars

Depending on mass and metallicity as well as evolutionary phase, stars occasionally experience convective-reactive nucleosynthesis episodes. We specifically investigate the situation when nucleosynthetically unprocessed, H-rich material is convectively mixed with a He-burning zone, for example in convectively unstable shell on top of electron-degenerate cores in AGB stars, young white dwarfs or X-ray bursting neutron stars. Such episodes are frequently encountered in stellar evolution models of stars of extremely low or zero metal content [...] We focus on the convective-reactive episode in the very-late thermal pulse star Sakurai's object (V4334 Sagittarii). Asplund etal. (1999) determined the abundances of 28 elements, many of which are highly non-solar, ranging from H, He and Li all the way to Ba and La, plus the C isotopic ratio. Our simulations show that the mixing evolution according to standard, one-dimensional stellar evolution models implies neutron densities in the He that are too low to obtain a significant neutron capture nucleosynthesis on the heavy elements. We have carried out 3D hydrodynamic He-shell flash convection [...] we assume that the ingestion process of H into the He-shell convection zone leads only after some delay time to a sufficient entropy barrier that splits the convection zone [...] we obtain significantly higher neutron densities (~few 10^15 1/cm^3) and reproduce the key observed abundance trends found in Sakurai's object. These include an overproduction of Rb, Sr and Y by about 2 orders of magnitude higher than the overproduction of Ba and La. Such a peculiar nucleosynthesis signature is impossible to obtain with the mixing predictions in our one-dimensional stellar evolution models. [...] We determine how our results depend on uncertainties of nuclear reaction rates, for example for the C13(α, n)O16 reaction.

preprint2011arXiv

Remnants of Binary White Dwarf Mergers

We carry out a comprehensive smooth particle hydrodynamics simulation survey of double-degenerate white dwarf binary mergers of varying mass combinations in order to establish correspondence between initial conditions and remnant configurations. We find that all but one of our simulation remnants share general properties such as a cold, degenerate core surrounded by a hot disk, while our least massive pair of stars forms only a hot disk. We characterize our remnant configurations by the core mass, the rotational velocity of the core, and the half-mass radius of the disk. We also find that some of our simulations with very massive constituent stars exhibit helium detonations on the surface of the primary star before complete disruption of the secondary. However, these helium detonations are insufficiently energetic to ignite carbon, and so do not lead to prompt carbon detonations.

preprint2011arXiv

Simulating the Common Envelope Phase of a Red Giant Using SPH and Uniform Grid Codes

We use three-dimensional hydrodynamical simulations to study the rapid infall phase of the common envelope interaction of a red giant branch star of mass equal to 0.88 \msun and a companion star of mass ranging from 0.9 down to 0.1 \msun. We first compare the results obtained using two different numerical techniques with different resolutions, and find overall very good agreement. We then compare the outcomes of those simulations with observed systems thought to have gone through a common envelope. The simulations fail to reproduce those systems in the sense that most of the envelope of the donor remains bound at the end of the simulations and the final orbital separations between the donor's remnant and the companion, ranging from 26.8 down to 5.9 \rsun, are larger than the ones observed. We suggest that this discrepancy vouches for recombination playing an essential role in the ejection of the envelope and/or significant shrinkage of the orbit happening in the subsequent phase.

preprint2010arXiv

56Ni Production in Double Degenerate White Dwarf Collisions

We present a comprehensive study of white dwarf collisions as an avenue for creating type Ia supernovae. Using a smooth particle hydrodynamics code with a 13-isotope, α-chain nuclear network, we examine the resulting 56Ni yield as a function of total mass, mass ratio, and impact parameter. We show that several combinations of white dwarf masses and impact parameters are able to produce sufficient quantities of 56Ni to be observable at cosmological distances. We find the 56Ni production in double-degenerate white dwarf collisions ranges from sub-luminous to the super-luminous, depending on the parameters of the collision. For all mass pairs, collisions with small impact parameters have the highest likelihood of detonating, but 56Ni production is insensitive to this parameter in high-mass combinations, which significantly increases their likelihood of detection. We also find that the 56Ni dependence on total mass and mass ratio is not linear, with larger mass primaries producing disproportionately more 56Ni than their lower mass secondary counterparts, and symmetric pairs of masses producing more 56Ni than asymmetric pairs.

preprint2010arXiv

The effect of 12C + 12C rate uncertainties on the weak s-process component

The contribution by massive stars (M > 9 solar masses) to the weak s-process component of the solar system abundances is primarily due to the 22Ne neutron source, which is activated near the end of helium-core burning. The residual 22Ne left over from helium-core burning is then reignited during carbon burning, initiating further s-processing that modifies the isotopic distribution. This modification is sensitive to the stellar structure and the carbon burning reaction rate. Recent work on the 12C + 12C reaction suggests that resonances located within the Gamow peak may exist, causing a strong increase in the astrophysical S-factor and consequently the reaction rate. To investigate the effect of an increased rate, 25 solar mass stellar models with three different carbon burning rates, at solar metallicity, were generated using the Geneva Stellar Evolution Code (GENEC) with nucleosynthesis post-processing calculated using the NuGrid Multi-zone Post-Processing Network code (MPPNP). The strongest rate caused carbon burning to occur in a large convective core rather than a radiative one. The presence of this large convective core leads to an overlap with the subsequent convective carbon-shell, significantly altering the initial composition of the carbon-shell. In addition, an enhanced rate causes carbon-shell burning episodes to ignite earlier in the evolution of the star, igniting the 22Ne source at lower temperatures and reducing the neutron density.

preprint2006arXiv

Probing the Density in the Galactic Center Region: Wind-Blown Bubbles and High-Energy Proton Constraints

Recent observations of the Galactic center in high-energy gamma-rays (above 0.1TeV) have opened up new ways to study this region, from understanding the emission source of these high-energy photons to constraining the environment in which they are formed. We present a revised theoretical density model of the inner 5pc surrounding Sgr A* based on the fact that the underlying structure of this region is dominated by the winds from the Wolf-Rayet stars orbiting Sgr A*. An ideal probe and application of this density structure is this high energy gamma-ray emission. We assume a proton-scattering model for the production of these gamma-rays and then determine first whether such a model is consistent with the observations and second whether we can use these observations to further constrain the density distribution in the Galactic center.

preprint2006arXiv

The Environments around Long-Duration Gamma-Ray Burst Progenitors

Gamma-ray burst (GRB) afterglow observations have allowed us to significantly constrain the engines producing these energetic explosions. Te redshift and position information provided by these afterglows have already allowed us to limit the progenitors of GRBs to only a few models. The afterglows may also provide another observation that can place further constraints on the GRB progenitor: measurements telling us about the environments surrounding GRBs. Current analyses of GRB afterglows suggest that roughly half of long-duration gamma-ray bursts occur in surroundings with density profiles that are uniform. We study the constraints placed by this observation on both the classic ``collapsar'' massive star progenitor and its relative, the ``helium-merger'' progenitor. We study several aspects of wind mass-loss and find that our modifications to the standard Wolf-Rayet mass-loss paradigm are not sufficient to produce constant density profiles. Although this does not rule out the standard ``collapsar'' progenitor, it does suggest a deficiency with this model. We then focus on the He-merger models and find that such progenitors can fit this particular constraint well. We show how detailed observations can not only determine the correct progenitor for GRBs, but also allow us to study binary evolution physics.