Source author record

C. E. Fields

C. E. Fields 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

7works
3topics
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

7 published item(s)

preprint2022arXiv

Testing White Dwarf Age Estimates using Wide Double White Dwarf Binaries from Gaia EDR3

White dwarf (WD) stars evolve simply and predictably, making them reliable age indicators. However, self-consistent validation of the methods for determining WD total ages has yet to be widely performed. This work uses 1565 wide ( > 100 au) WD+WD binaries and 24 new triples containing at least two WDs to test the accuracy and validity of WD total age determinations. For these 1589 wide double-WD binaries and triples, we derive total ages of each WD using photometric data from all-sky surveys, in conjunction with Gaia parallaxes and current hydrogen-atmosphere WD models. Ignoring initial-to-final-mass relations and considering only WD cooling ages, we find that roughly 21-36% of the more massive WDs in a system have a shorter cooling age. Since more massive WDs should be born as more massive main-sequence stars, we attribute this unphysical disagreement as evidence of prior mergers or the presence of an unresolved companion, suggesting that roughly 21-36% of wide WD+WD binaries were once triples. Among the 423 wide WD+WD pairs that pass high-fidelity cuts, we find that 25% total age uncertainties are generally appropriate for WDs with masses > 0.63 Msun and temperatures < 12,000 K, and provide suggested inflation factors for age uncertainties for higher-mass WDs. Overall, WDs return reliable stellar ages, but we detail cases where total ages are least reliable, especially for WDs < 0.63 Msun.

preprint2021arXiv

The Three-Dimensional Collapse of a Rapidly Rotating 16 $M_{\odot}$ Star

We report on the three-dimensional (3D) hydrodynamic evolution to iron core-collapse of a rapidly rotating 16 $M_{\odot}$ star. For the first time, we follow the 3D evolution of the angular momentum (AM) distribution in the iron core and convective shell burning regions for the final 10 minutes up to and including gravitational instability and core-collapse. In 3D, we find that convective regions show efficient AM transport that leads to an AM profile that differs in shape and magnitude from $\texttt{MESA}$ within a few shell convective turnover timescales. For different progenitor models, such as those with tightly coupled Si/O convective shells, efficient AM transport in 3D simulations could lead to a significantly different AM distribution in the stellar interior affecting estimates of the natal neutron star or black hole spin. Our results suggest that 3D AM transport in convective and rotating shell burning regions are critical components in models of massive stars and could qualitatively alter the explosion outcome and inferred compact remnant properties.

preprint2020arXiv

On The Development of Multidimensional Progenitor Models For Core-collapse Supernovae

Multidimensional hydrodynamic simulations of shell convection in massive stars suggest the development of aspherical perturbations that may be amplified during iron core-collapse. These perturbations have a crucial and qualitative impact on the delayed neutrino-driven core-collapse supernova explosion mechanism by increasing the total stress behind the stalled shock. In this paper, we investigate the properties of a 15 \msun model evolved in 1-,2-, and 3-dimensions (3D) for the final $\sim$424 seconds before gravitational instability and iron core-collapse using MESA and the FLASH simulation framework. We find that just before collapse, our initially perturbed fully 3D model reaches angle-averaged convective velocity magnitudes of $\approx$ 240-260 km s$^{-1}$ in the Si- and O-shell regions with a Mach number $\approx$ 0.06. We find the bulk of the power in the O-shell resides at large scales, characterized by spherical harmonic orders ($\ell$) of 2-4, while the Si-shell shows broad spectra on smaller scales of $\ell\approx30-40$. Both convective regions show an increase in power at $\ell=5$ near collapse. We show that the 1D \texttt{MESA} model agrees with the convective velocity profile and speeds of the Si-shell when compared to our highest resolution 3D model. However, in the O-shell region, we find that \texttt{MESA} predicts speeds approximately \emph{four} times slower than all of our 3D models suggest. All eight of the multi-dimensional stellar models considered in this work are publicly available.

preprint2020arXiv

The evolutionary roads leading to low effective spins, high black hole masses, and O1/O2 rates of LIGO/Virgo binary black holes

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

preprint2016arXiv

On Variations Of Pre-Supernova Model Properties

We explore the variation in single star 15-30 $\rm{M}_{\odot}$, non-rotating, solar metallicity, pre-supernova MESA models due to changes in the number of isotopes in a fully-coupled nuclear reaction network and adjustments in the mass resolution. Within this two-dimensional plane we quantitatively detail the range of core masses at various stages of evolution, mass locations of the main nuclear burning shells, electron fraction profiles, mass fraction profiles, burning lifetimes, stellar lifetimes, and compactness parameter at core-collapse for models with and without mass loss. Up to carbon burning we generally find mass resolution has a larger impact on the variations than the number of isotopes, while the number of isotopes plays a more significant role in determining the span of the variations for neon, oxygen and silicon burning. Choice of mass resolution dominates the variations in the structure of the intermediate convection zone and secondary convection zone during core and shell hydrogen burning respectively, where we find a minimum mass resolution of $\approx$ 0.01 $\rm{M}_{\odot}$ is necessary to achieve convergence in the helium core mass at the $\approx$5% level. On the other hand, at the onset of core-collapse we find $\approx$30% variations in the central electron fraction and mass locations of the main nuclear burning shells, a minimum of $\approx$127 isotopes is needed to attain convergence of these values at the $\approx$10% level.

preprint2016arXiv

Properties of Carbon-Oxygen White Dwarfs From Monte Carlo Stellar Models

We investigate properties of carbon-oxygen white dwarfs with respect to the composite uncertainties in the reaction rates using the stellar evolution toolkit, Modules for Experiments in Stellar Astrophysics (MESA) and the probability density functions in the reaction rate library STARLIB. These are the first Monte Carlo stellar evolution studies that use complete stellar models. Focusing on 3 M$_{\odot}$ models evolved from the pre main-sequence to the first thermal pulse, we survey the remnant core mass, composition, and structure properties as a function of 26 STARLIB reaction rates covering hydrogen and helium burning using a Principal Component Analysis and Spearman Rank-Order Correlation. Relative to the arithmetic mean value, we find the width of the 95\% confidence interval to be $ΔM_{\rm 1TP}$ $\approx$ 0.019 M$_{\odot}$ for the core mass at the first thermal pulse, $Δ$$t_{\rm{1TP}}$ $\approx$ 12.50 Myr for the age, $Δ\log(T_{\rm c}/{\rm K}) \approx$ 0.013 for the central temperature, $Δ\log(ρ_{\rm c}/{\rm g \ cm}^{-3}) \approx$ 0.060 for the central density, $ΔY_{\rm{e,c}} \approx$ 2.6$\times$10$^{-5}$ for the central electron fraction, $ΔX_{\rm c}(^{22}\rm{Ne}) \approx$ 5.8$\times$10$^{-4}$, $ΔX_{\rm c}(^{12}\rm{C}) \approx$ 0.392, and $ΔX_{\rm c}(^{16}\rm{O}) \approx$ 0.392. Uncertainties in the experimental $^{12}$C($α,γ)^{16}\rm{O}$, triple-$α$, and $^{14}$N($p,γ)^{15}\rm{O}$ reaction rates dominate these variations. We also consider a grid of 1 to 6 M$_{\odot}$ models evolved from the pre main-sequence to the final white dwarf to probe the sensitivity of the initial-final mass relation to experimental uncertainties in the hydrogen and helium reaction rates.

preprint2015arXiv

On Carbon Burning in Super Asymptotic Giant Branch Stars

We explore the detailed and broad properties of carbon burning in Super Asymptotic Giant Branch (SAGB) stars with 2755 MESA stellar evolution models. The location of first carbon ignition, quenching location of the carbon burning flames and flashes, angular frequency of the carbon core, and carbon core mass are studied as a function of the ZAMS mass, initial rotation rate, and mixing parameters such as convective overshoot, semiconvection, thermohaline and angular momentum transport. In general terms, we find these properties of carbon burning in SAGB models are not a strong function of the initial rotation profile, but are a sensitive function of the overshoot parameter. We quasi-analytically derive an approximate ignition density, $ρ_{ign} \approx 2.1 \times 10^6$ g cm$^{-3}$, to predict the location of first carbon ignition in models that ignite carbon off-center. We also find that overshoot moves the ZAMS mass boundaries where off-center carbon ignition occurs at a nearly uniform rate of $ΔM_{\rm ZAMS}$/$Δf_{\rm{ov}}\approx$ 1.6 $M_{\odot}$. For zero overshoot, $f_{\rm{ov}}$=0.0, our models in the ZAMS mass range $\approx$ 8.9 to 11 $M_{\odot}$ show off-center carbon ignition. For canonical amounts of overshooting, $f_{\rm{ov}}$=0.016, the off-center carbon ignition range shifts to $\approx$ 7.2 to 8.8 $M_{\odot}$. Only systems with $f_{\rm{ov}}$ $\geq 0.01$ and ZAMS mass $\approx$ 7.2-8.0 $M_{\odot}$ show carbon burning is quenched a significant distance from the center. These results suggest a careful assessment of overshoot modeling approximations on claims that carbon burning quenches an appreciable distance from the center of the carbon core.