Source author record

Patrick A. Young

Patrick A. Young 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

21works
5topics
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

21 published item(s)

preprint2022arXiv

Mantle Degassing Lifetimes through Galactic Time and the Maximum Age Stagnant-lid Rocky Exoplanets can Support Temperate Climates

The ideal exoplanets to search for life are those within a star's habitable zone. However, even within the habitable zone planets can still develop uninhabitable climate states. Sustaining a temperate climate over geologic ($\sim$Gyr) timescales requires a planet contain sufficient internal energy to power a planetary-scale carbon cycle. A major component of a rocky planet's energy budget is the heat produced by the decay of radioactive elements, especially $^{40}$K, $^{232}$Th, $^{235}$U and $^{238}$U. As the planet ages and these elements decay, this radiogenic energy source dwindles. Here we estimate the probability distribution of the amount of these heat producing elements (HPEs) that enter into rocky exoplanets through Galactic history, by combining the system-to-system variation seen in stellar abundance data with the results from Galactic chemical evolution models. Using these distributions, we perform Monte-Carlo thermal evolution models that maximize the mantle cooling rate. This allows us to create a pessimistic estimate of lifetime a rocky, stagnant-lid exoplanet can support a global carbon cycle and temperate climate as a function of its mass and when it in Galactic history. We apply this framework to a sample of 17 likely rocky exoplanets with measured ages, 7 of which we predict are likely to be actively degassing today despite our pessimistic assumptions. For the remaining planets, including those orbiting TRAPPIST-1, we cannot confidently assume they currently contain sufficient internal heat to support mantle degassing at a rate sufficient to sustain a global carbon cycle or temperate climate without additional tidal heating or undergoing plate tectonics.

preprint2020arXiv

The Influence of Stellar Phosphorus On Our Understanding of Exoplanets and Astrobiology

When searching for exoplanets and ultimately considering their habitability, it is necessary to consider the planet's composition, geophysical processes, and geochemical cycles in order to constrain the bioessential elements available to life. Determining the elemental ratios for exoplanetary ecosystems is not yet possible, but we generally assume that planets have compositions similar to those of their host stars. Therefore, using the Hypatia Catalog of high-resolution stellar abundances for nearby stars, we compare the C, N, Si, and P abundance ratios of main sequence stars with those in average marine plankton, Earth's crust, as well as bulk silicate Earth and Mars. We find that, in general, plankton, Earth, and Mars are N-poor and P-rich compared with nearby stars. However, the dearth of P abundance data, which exists for only ~1% of all stars and 1% of exoplanet hosts, makes it difficult to deduce clear trends in the stellar data, let alone the role of P in the evolution of an exoplanet. Our Sun has relatively high P and Earth biology requires a small, but finite, amount of P. On rocky planets that form around host stars with substantially less P, the strong partitioning of P into the core could rule out the potential for surface P and, consequently, for life on that planet's surface. Therefore, we urge the stellar abundance community to make P observations a priority in future studies and telescope designs.

preprint2020arXiv

Three-dimensional Supernova Models Provide New Insights into the Origins of Stardust

We present the isotope yields of two post-explosion, three-dimensional 15 $M_\odot$ core-collapse supernova models, 15S and 15A, and compare them to the carbon, nitrogen, silicon, aluminum, sulfur, calcium, titanium, iron, and nickel isotopic compositions of SiC stardust. We find that these core-collapse supernova models predict similar carbon and nitrogen compositions to SiC X grains and grains with $^{12}$C/$^{13}$C $<$ 20 and $^{14}$N/$^{15}$N $<$ 60, which we will hereafter refer to as SiC 'D' grains. Material from the interior of a 15 $M_\odot$ explosion reaches high enough temperatures shortly after core collapse to produce the large enrichments of $^{13}$C and $^{15}$N necessary to replicate the compositions of SiC D grains. The innermost ejecta in a core-collapse supernova is operating in the neutrino-driven regime and undergoes fast proton capture after being heated by the supernova shockwave. Both 3-D models predict 0.3 $<$ $^{26}$Al/$^{27}$Al $<$ 1.5, comparable to the ratios seen in SiC X, C, and D grains. Models 15S and 15A, in general, predict very large anomalies in calcium isotopes but do compare qualitatively with the SiC X grain measurements that show $^{44}$Ca and $^{43}$Ca excesses. The titanium isotopic compositions of SiC X grains are well reproduced. The models predict $^{57}$Fe excesses and depletions that are observed in SiC X grains, and in addition predict accurately the $^{60}$Ni/$^{58}$Ni, $^{61}$Ni/$^{58}$Ni, and $^{62}$Ni/$^{58}$Ni ratios in SiC X grains, as a result of fast neutron captures initiated by the propagation of the supernova shockwave. Finally, symmetry has a noticeable effect on the production of silicon, sulfur, and iron isotopes in the SN ejecta.

preprint2020arXiv

Titanium and Iron in the Cassiopeia A Supernova Remnant

Mixing above the proto-neutron star is believed to play an important role in the supernova engine, and this mixing results in a supernova explosion with asymmetries. Elements produced in the innermost ejecta, e.g., ${}^{56}$Ni and ${}^{44}$Ti, provide a clean probe of this engine. The production of ${}^{44}$Ti is particularly sensitive to the exact production pathway and, by understanding the available pathways, we can use ${}^{44}$Ti to probe the supernova engine. Using thermodynamic trajectories from a three-dimensional supernova explosion model, we review the production of these elements and the structures expected to form under the "convective-engine" paradigm behind supernovae. We compare our results to recent X-ray and $γ$-ray observations of the Cassiopeia A supernova remnant.

preprint2019arXiv

A Flexible Bayesian Framework for Assessing Habitability with Joint Observational and Model Constraints

The catalog of stellar evolution tracks discussed in our previous work is meant to help characterize exoplanet host-stars of interest for follow-up observations with future missions like JWST. However, the utility of the catalog has been predicated on the assumption that we would precisely know the age of the particular host-star in question; in reality, it is unlikely that we will be able to accurately estimate the age of a given system. Stellar age is relatively straightforward to calculate for stellar clusters, but it is difficult to accurately measure the age of an individual star to high precision. Unfortunately, this is the kind of information we should consider as we attempt to constrain the long-term habitability potential of a given planetary system of interest. This is ultimately why we must rely on predictions of accurate stellar evolution models, as well a consideration of what we can observably measure (stellar mass, composition, orbital radius of an exoplanet) in order to create a statistical framework wherein we can identify the best candidate systems for follow-up characterization. In this paper we discuss a statistical approach to constrain long-term planetary habitability by evaluating the likelihood that at a given time of observation, a star would have a planet in the 2 Gy continuously habitable zone (CHZ2). Additionally, we will discuss how we can use existing observational data (i.e. data assembled in the Hypatia catalog and the Kepler exoplanet host star database) for a robust comparison to the catalog of theoretical stellar models.

preprint2016arXiv

A Catalog of Stellar Evolution Profiles and the Effects of Variable Composition on Habitable Systems

We present stellar evolution models for 0.5 - 1.2 \Msol at scaled metallicities of 0.1 - 1.5 Z\sol and O/Fe values of 0.44 - 2.28 O/Fe\sol. The time dependent evolution of habitable zone boundaries are calculated for each stellar evolution track based on stellar mass, effective temperature, and luminosity parameterizations. The rate of change of stellar surface quantities and the surrounding habitable zone position are strong functions of all three quantities explored. The range of orbits that remain continuously habitable, or habitable for at least 2 Gyr, are provided. The results show that the detailed chemical characterization of exoplanet host stars and a consideration of their evolutionary history are necessary to assess the likelihood that a planet found in the instantaneous habitable zone has had sufficient time to develop a biosphere capable of producing detectable biosignatures. This model grid is designed for use by the astrobiology and exoplanet communities to efficiently characterize the time evolution of host stars and their habitable zones for planetary candidates of interest.

preprint2016arXiv

A Comparison of Stellar Elemental Abundance Techniques and Measurements

Stellar elemental abundances are important for understanding the fundamental properties of a star or stellar group, such as age and evolutionary history, as well as the composition of an orbiting planet. However, as abundance measurement techniques have progressed, there has been little standardization between individual methods and their comparisons. As a result, different stellar abundance procedures determine measurements that vary beyond quoted error for the same elements within the same stars (Hinkel et al. 2014). The purpose of this paper is to better understand the systematic variations between methods and offer recommendations for producing more accurate results in the future. We have invited a number of participants from around the world (Australia, Portugal, Sweden, Switzerland, and USA) to calculate ten element abundances (C, O, Na, Mg, Al, Si, Fe, Ni, Ba, and Eu) using the same stellar spectra for four stars (HD361, HD10700, HD121504, HD202206). Each group produced measurements for each of the stars using: 1) their own autonomous techniques, 2) standardized stellar parameters, 3) standardized line list, and 4) both standardized parameters and line list. We present the resulting stellar parameters, absolute abundances, and a metric of data similarity that quantifies homogeneity of the data. We conclude that standardization of some kind, particularly stellar parameters, improves the consistency between methods. However, because results did not converge as more free parameters were standardized, it is clear there are inherent issues within the techniques that need to be reconciled. Therefore, we encourage more conversation and transparency within the community such that stellar abundance determinations can be reproducible as well as accurate and precise.

preprint2015arXiv

The Chemical Composition of τ Ceti and Possible Effects on Terrestrial Planets

τ Ceti (HD10700), a G8 dwarf with solar mass of 0.78, is a close (3.65 pc) sun-like star where 5 possibly terrestrial planet candidates (minimum masses of 2, 3.1, 3.5, 4.3, and 6.7 Earth masses) have recently been discovered. We report abundances of 23 elements using spectra from the MIKE spectrograph on Magellan. Using stellar models with the abundances determined here, we calculate the position of the classical habitable zone with time. At the current best fit age, 7.63 Gy, up to two planets (e and f) may be in the habitable zone, depending on atmospheric properties. The Mg/Si ratio of the star is found to be 1.78, which is much greater than for Earth (about 1.2). With a system that has such an excess of Mg to Si ratio it is possible that the mineralogical make-up of planets around τ Ceti could be significantly different from that of Earth, with possible oversaturation of MgO, resulting in an increase in the content of olivine and ferropericlase compared with Earth. The increase in MgO would have a drastic impact on the rheology of the mantles of the planets around τ Ceti.

preprint2015arXiv

The Effects on Supernova Shock Breakout and Swift Light Curves Due to the Mass of the Hydrogen-Rich Envelope

Mass loss remains one of the primary uncertainties in stellar evolution. In the most massive stars, mass loss dictates the circumstellar medium and can significantly alter the fate of the star. Mass loss is caused by a variety of wind mechanisms and also through binary interactions. Supernovae are excellent probes of this mass loss, both the circumstellar material and the reduced mass of the hydrogen-rich envelope. In this paper, we focus on the effects of reducing the hydrogen-envelope mass on the supernova light curve, studying both the shock breakout and peak light curve emission for a wide variety of mass loss scenarios. Even though the trends of this mass loss will be masked somewhat by variations caused by different progenitors, explosion energies, and circumstellar media, these trends have significant effects on the supernova light-curves that should be seen in supernova surveys. We conclude with a comparison of our results to a few key observations.

preprint2014arXiv

Stellar Abundances in the Solar Neighborhood: The Hypatia Catalog

We compile spectroscopic abundance data from 84 literature sources for 50 elements across 3058 stars in the solar neighborhood, within 150 pc of the Sun, to produce the Hypatia Catalog. We evaluate the variability of the spread in abundance measurements reported for the same star by different surveys. We also explore the likely association of the star within the Galactic disk, the corresponding observation and abundance determination methods for all catalogs in Hypatia, the influence of specific catalogs on the overall abundance trends, and the effect of normalizing all abundances to the same solar scale. The resulting large number of stellar abundance determinations in the Hypatia Catalog are analyzed only for thin-disk stars with observations that are consistent between literature sources. As a result of our large dataset, we find that the stars in the solar neighborhood may be reveal an asymmetric abundance distribution, such that a [Fe/H]-rich group near to the mid-plane is deficient in Mg, Si, S, Ca, Sc II, Cr II, and Ni as compared to stars further from the plane. The Hypatia Catalog has a wide number of applications, including exoplanet hosts, thick and thin disk stars, or stars with different kinematic properties.

preprint2013arXiv

A Kine-Chemical Investigation of the AB Dor Moving Group "Stream"

The AB Dor Moving Group consists of a "nucleus" of ~10 stars at d ~ 20 pc, along with dozens of purported "stream" members distributed across the sky. We perform a chemical and kinematic analysis of a subsample of AB Dor stream stars to test whether they constitute a physical stellar group. We use the NEMO Galactic kinematic code to investigate the orbits of the stream members, and perform a chemical abundance analysis using high resolution spectra taken with the Magellan Clay 6.5-m telescope. Using a chi^2 test with the measured abundances for 10 different elements, we find that only half of the purported AB Dor stream members could possibly constitute a statistically chemically homogeneous sample. Some stream members with 3D velocities were hundreds of parsecs from the AB Dor nucleus ~10^8 yr ago, and hence were unlikely to share a common origin. We conclude that the published lists of AB Dor moving group stream members are unlikely to represent the dispersed remnant of a single star formation episode. A subsample of the stream stars appears to be both statistically chemically homogeneous and in the vicinity of the AB Dor nucleus at birth. Their mean metallicity is [Fe/H] = 0.02 +/- 0.02 dex, which we consider representative for the AB Dor group. Finally, we report a strong lower limit on the age of the AB Dor nucleus of >110 Myr based on the pre-MS contraction times for K-type members which have reached the main sequence.

preprint2013arXiv

Can Stellar Mixing Explain the Lack of Type Ib Supernovae in Long-Duration GRBs?

The discovery of supernovae associated with long-duration gamma ray burst observations is primary evidence that the progenitors of these outbursts are massive stars. One of the principle mysteries in understanding these progenitors has been the fact that all of these gamma-ray burst associated supernovae are type Ic supernovae, with no evidence of helium in the stellar atmosphere. Many studies have focused on whether or not this helium is simply hidden from spectral analyses. In this paper, we show results from recent stellar models using new convection algorithms based on our current understanding of stellar mixing. We demonstrate that enhanced convection may lead to severe depletion of stellar helium layers, suggesting that the helium is not observed simply because it is not in the star. We also present light-curves and spectra of these compact helium-depleted stars, compared to models with more conventional helium layers.

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.

preprint2013arXiv

The Age and Stellar Parameters of the Procyon Binary System

The Procyon AB binary system (orbital period 40.838 years, a newly-refined determination), is near and bright enough that the component radii, effective temperatures, and luminosities are very well determined, although more than one possible solution to the masses has limited the claimed accuracy. Preliminary mass determinations for each component are available from HST imaging, supported by ground-based astrometry and an excellent Hipparcos parallax; we use these for our preferred solution for the binary system. Other values for the masses are also considered. We have employed the TYCHO stellar evolution code to match the radius and luminosity of the F5 IV-V primary star to determine the system's most likely age as 1.87 +/- 0.13 Gyr. Since prior studies of Procyon A found its abundance indistinguishable from solar, the solar composition of Asplund, Grevesse & Sauval (Z=0.014) is assumed for the HR Diagram flitting. An unsuccessful attempt to fit using the older solar abundance scale of Grevesse & Sauval (Z=0.019) is also reported. For Procyon B, eleven new sequences for the cooling of non-DA white dwarfs have been calculated, to investigate the dependence of the cooling age on (1) the mass, (2) the core composition, (3) the helium layer mass, and (4) heavy-element opacities in the helium envelope. Our calculations indicate a cooling age of 1.19+/-0.11 Gyr, which implies that the progenitor mass of Procyon B was 2.59(+0.44,-0.26) Msun. In a plot of initial vs final mass of white dwarfs in astrometric binaries or star clusters (all with age determinations), the Procyon B final mass lies several sigma below a straight line fit.

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.

preprint2012arXiv

The Impact of Stellar Abundance Variations on Stellar Habitable Zone Evolution

The high quality spectra required for radial velocity planet searches are well-suited to providing abundances for a wide array of elements in large samples of stars. Abundance ratios of the most common elements relative to Fe are observed to vary by more than a factor of two in planet host candidates. This level of variation has a substantial impact on the evolution of the host star and the extent of its habitable zone. We present stellar models of 1 solar mass stars with custom compositions representing the full range of these non-solar abundance ratios. We find that the effects derived from variation over the observed range of [O/Fe] are particularly dramatic. Habitability lifetimes for some classes of orbits can vary by gigayears for the observed range in [O/Fe].

preprint2010arXiv

Trends in Ti44 and Ni56 from Core-Collapse Supernovae

We compare the yields of Ti44 and Ni56 produced from post-processing the thermodynamic trajectories from three different core-collapse models -- a Cassiopeia A progenitor, a double shock hypernova progenitor, and a rotating 2D explosion -- with the yields from exponential and power-law trajectories. The peak temperatures and densities achieved in these core-collapse models span several of the distinct nucleosynthesis regions we identify, resulting in different trends in the Ti44 and Ni56 yields for different mass elements. The Ti44 and Ni56 mass fraction profiles from the exponential and power-law profiles generally explain the tendencies of the post-processed yields, depending on which regions are traversed by the model. We find integrated yields of Ti44 and Ni56 from the exponential and power-law trajectories are generally within a factor 2 or less of the post-process yields. We also analyze the influence of specific nuclear reactions on the Ti44 and Ni56 abundance evolution. Reactions that affect all yields globally are the 3a, p(e-,nu)n and n(e+,nubar)p. The rest of the reactions are ranked according to their degree of impact on the synthesis of Ti44. The primary ones include Ti44(a,p)V47, Ca40(a,g)Ti44, V45(p,g)Cr46, Ca40(a,p)Sc43, F17(a,p)Ne20, Na21(a,p)Mg24, Sc41(p,g)Ti42, Sc43(p,g)Ti44, Ti44(p,g)V45, and Ni57(p,g)Cu58, along with numerous weak reactions. Our analysis suggests that not all Ti44 need be produced in an a-rich freeze-out in core-collapse events, and that reaction rate equilibria in combination with timescale effects for the expansion profile may account for the paucity of Ti44 observed in supernovae remnants.

preprint2009arXiv

Convection Theory and Sub-photospheric Stratification

As a step toward a complete theoretical integration of 3D compressible hydrodynamic simulations into stellar evolution, convection at the surface and sub-surface layers of the Sun is re-examined, from a restricted point of view, in the language of mixing-length theory (MLT) . Requiring that MLT use a hydrodynamically realistic dissipation length gives a new constraint on solar models. While the stellar structure which results is similar to that obtained by YREC and Garching models, the theoretical picture differs. A new quantitative connection is made between macro-turbulence, micro-turbulence, and the convective velocity scale at the photosphere, which has finite values. The "geometric parameter" in MLT is found to correspond more reasonably with the size of the strong downward plumes which drive convection (Stein and Nordlund 1998), and thus has a physical interpretation even in MLT. Use of 3D simulations of both adiabatic convection and stellar atmospheres will allow the determination of the dissipation length and the geometric parameter (i.e., the entropy jump), with no astronomical calibration. A physically realistic treatment of convection in stellar evolution will require additional modifications beyond MLT, including effects of kinetic energy flux, entrainment (the most dramatic difference from MLT found by Meakin and Arnett 2007), rotation, and magnetic fields (Balbus 2009}.

preprint2009arXiv

Spectra and Light Curves of Failed Supernovae

Astronomers have proposed a number of mechanisms to produce supernova explosions. Although many of these mechanisms are now not considered primary engines behind supernovae, they do produce transients that will be observed by upcoming ground-based surveys and NASA satellites. Here we present the first radiation-hydrodynamics calculations of the spectra and light curves from three of these "failed" supernovae: supernovae with considerable fallback, accretion induced collapse of white dwarfs, and energetic helium flashes (also known as type .Ia supernovae).

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.

preprint2006arXiv

Uncertainties in Supernova Yields I: 1D Explosions

Theoretical nucleosynthetic yields from supernovae are sensitive to both the details of the progenitor star and the explosion calculation. We attempt to comprehensively identify the sources of uncertainties in these yields. In this paper we concentrate on the variations in yields from a single progenitor arising from common 1-dimensional methods of approximating a supernova explosion. Subsequent papers will examine 3-dimensional effects in the explosion and the progenitor, and trends in mass and composition. For the 1-dimensional explosions we find that both elemental and isotopic yields for Si and heavier elements are a sensitive function of explosion energy. Also, piston-driven and thermal bomb type explosions have different yields for the same explosion energy. Yields derived from 1-dimensional explosions are non-unique.