Source author record

Christopher L. Fryer

Christopher L. Fryer 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

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

22 published item(s)

preprint2023arXiv

ASCENT - A balloon-borne hard X-ray imaging spectroscopy telescope using transition edge sensor microcalorimeter detectors

Core collapse supernovae are thought to be one of the main sources in the galaxy of elements heavier than iron. Understanding the origin of the elements is thus tightly linked to our understanding of the explosion mechanism of supernovae and supernova nucleosynthesis. X-ray and gamma-ray observations of young supernova remnants, combined with improved theoretical modeling, have resulted in enormous improvements in our knowledge of these events. The isotope ${}^{44}$Ti is one of the most sensitive probes of the innermost regions of the core collapse engine, and its spatial and velocity distribution are key observables. Hard X-ray imaging spectroscopy with the Nuclear Spectroscopic Telescope Array (NuSTAR) has provided new insights into the structure of the supernova remnant Cassiopeia A (Cas A), establishing the convective nature of the supernova engine. However, many questions about the details of this engine remain. We present here the concept for a balloon-borne follow-up mission called ASCENT (A SuperConducting ENergetic x-ray Telescope). ASCENT uses transition edge sensor gamma-ray microcalorimeter detectors with a demonstrated 55 eV Full Width Half Maximum (FWHM) energy resolution at 97 keV. This 8--16-fold improvement in energy resolution over NuSTAR will allow high resolution imaging and spectroscopy of the ${}^{44}$Ti emission. This will allow a detailed reconstruction of gamma-ray line redshifts, widths, and shapes, allowing us to address questions such as: What is the source of the neutron star "kicks"? What is the dominant production pathway for ${}^{44}$Ti? Is the engine of Cas A unique?

preprint2022arXiv

Advancing the Landscape of Multimessenger Science in the Next Decade

The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science.

preprint2021arXiv

Dependence of Dust Formation on the Supernova Explosion

We investigate the properties, composition, and dynamics of dust formation and growth for a diverse set of core-collapse supernovae (CCSNe), varying the progenitor mass, explosion energy, and engine type. These explosions are evolved with a 1-D Lagrangian hydrodynamics code out to several hundred days to model the ejecta as it expands and cools. A multigrain dust nucleation and growth model is applied to these results. We find that higher explosion energies lead to an earlier onset of dust formation, smaller grain sizes, and larger silicate abundances. Further, we see that nuclear burning during the explosion leads to enhanced formation of silicate dust. Finally, we build composite models from our suite to predict the efficiency of CCSNe dust production as a function of metallicity.

preprint2021arXiv

Kilonova Detectability with Wide-Field Instruments

Kilonovae are ultraviolet, optical, and infrared transients powered by the radioactive decay of heavy elements following a neutron star merger. Joint observations of kilonovae and gravitational waves can offer key constraints on the source of Galactic $r$-process enrichment, among other astrophysical topics. However, robust constraints on heavy element production requires rapid kilonova detection (within $\sim 1$ day of merger) as well as multi-wavelength observations across multiple epochs. In this study, we quantify the ability of 13 wide field-of-view instruments to detect kilonovae, leveraging a large grid of over 900 radiative transfer simulations with 54 viewing angles per simulation. We consider both current and upcoming instruments, collectively spanning the full kilonova spectrum. The Roman Space Telescope has the highest redshift reach of any instrument in the study, observing kilonovae out to $z \sim 1$ within the first day post-merger. We demonstrate that BlackGEM, DECam, GOTO, the Vera C. Rubin Observatory's LSST, ULTRASAT, and VISTA can observe some kilonovae out to $z \sim 0.1$ ($\sim$475 Mpc), while DDOTI, MeerLICHT, PRIME, $Swift$/UVOT, and ZTF are confined to more nearby observations. Furthermore, we provide a framework to infer kilonova ejecta properties following non-detections and explore variation in detectability with these ejecta parameters.

preprint2020arXiv

BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code's BSE sector, such as the 'rapid' and the 'delayed' supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster.

preprint2020arXiv

Gamma-rays from kilonova: a potential probe of r-process nucleosynthesis

The mergers of compact binaries with at least one neutron star component are the potential leading sites of the production and ejection of $r$-process elements. Discoveries of galactic binary pulsars, short gamma-ray bursts, and gravitational-wave detections have all been constraining the rate of these events, while the gravitational wave plus broadband electromagnetic coverage of binary neutron star merger (GW170817) has also placed constraints on the properties (mass and composition) of the merger ejecta. But uncertainties and ambiguities in modeling the optical and infrared emission make it difficult to definitively measure the distribution of heavy isotopes in these mergers. In contrast, gamma rays emitted in the decay of these neutron-rich ejecta may provide a more direct measurement of the yields. We calculate the gamma production in remnants of neutron star mergers, considering two epochs: a kilonova epoch, lasting about two weeks, and a much later epoch of tens and hundreds of thousands of years after the merger. For the kilonova epoch, when the expanding ejecta is still only partially transparent to gamma radiation, we use 3D radiative transport simulations to produce the spectra. We show that the gamma-ray spectra associated with beta- and alpha-decay provide a fingerprint of the ejecta properties and, for a sufficiently nearby remnant, may be detectable, even for old remnants. We compare our gamma spectra with the potential detection limits of next generation detectors, including the Lunar Occultation Explorer (LOX), the All-sky Medium Energy Gamma-ray Observatory (AMEGO), and the Compton Spectrometer and Imager (COSI). We show that fission models can be discriminated via the presence of short-lived fission fragments in the remnant spectra.

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

Composition Effects on Kilonova Spectra and Light Curves: I

The merger of neutron star binaries is believed to eject a wide range of heavy elements into the universe. By observing the emission from this ejecta, scientists can probe the ejecta properties (mass, velocity and composition distributions). The emission (a.k.a. kilonova) is powered by the radioactive decay of the heavy isotopes produced in the merger and this emission is reprocessed by atomic opacities to optical and infra-red wavelengths. Understanding the ejecta properties requires calculating the dependence of this emission on these opacities. The strong lines in the optical and infra-red in lanthanide opacities have been shown to significantly alter the light-curves and spectra in these wavelength bands, arguing that the emission in these wavelengths can probe the composition of this ejecta. Here we study variations in the kilonova emission by varying individual lanthanide (and the actinide uranium) concentrations in the ejecta. The broad forest of lanthanide lines makes it difficult to determine the exact fraction of individual lanthanides. Nd is an exception. Its opacities above 1 micron are higher than other lanthanides and observations of kilonovae can potentially probe increased abundances of Nd. Similarly, at early times when the ejecta is still hot (first day), the U opacity is strong in the 0.2-1 micron wavelength range and kilonova observations may also be able to constrain these abundances.

preprint2016arXiv

Particle Acceleration in Relativistic Electron-Ion Outlfows

We use the Los Alamos VPIC code to investigate particle acceleration in relativistic, unmagnetized, collisionless electron-ion plasmas. We run our simulations both with a realistic proton-to-electron mass ratio m_p/m_e = 1836, as well as commonly employed mass ratios of m_p/m_e =100 and 25, and show that results differ among the different cases. In particular, for the physically accurate mass ratio, electron acceleration occurs efficiently in a narrow region of a few hundred inertial lengths near the flow front, producing a power law dN/dgamma ~ gamma^(-p) with p ~ -2 developing over a few decades in energy, while acceleration is weak in the region far downstream. We find 20%, 10%, and 0.2% of the total energy given to the electrons for mass ratios of 25, 100, and 1836 respectively at a time of 2500 (w_p)^-1. Our simulations also show significant magnetic field generation just ahead of and behind the the flow front, with about 1% of the total energy going into the magnetic field for a mass ratio of 25 and 100, and 0.1% for a mass ratio of and 1836. In addition, lower mass ratios show significant fields much further downstream than in the realistic mass ratio case. Our results suggest the region and energetic extent of particle acceleration is directly related to the presence of magnetic field generation. Our work sheds light on the understanding of particle acceleration and emission in gamma-ray bursts, among other relativistic astrophysical outflows, but also underscores the necessity of optimizing numerical and physical parameters, as well as comparing among PIC codes before firm conclusions are drawn from these types of simulations.

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.

preprint2012arXiv

The Growth of the Stellar Seeds of Supermassive Black Holes

One of the most promising explanations for the origin of the billion solar mass black holes (BHs) inferred to power quasars at redshifts z > 6 is that supermassive stars (SMSs) with masses > 10,000 solar masses collapse to form the seed BHs from which they grow. Here we review recent theoretical advances which provide support for this scenario. Firstly, given sufficiently high accretion rates of gas into the cores of primordial protogalaxies, it appears that neither the high energy radiation emitted from the stellar surface nor the limited lifetime of SMSs can prevent their growth to masses of up to > 100,000 solar masses. Secondly, recent cosmological simulations suggest that the high fluxes of molecule-dissociating radiation which may be required in order to achieve such high accretion rates may be more common in the early universe than previously thought. We conclude that the majority of supermassive BHs may originate from SMSs at high redshifts.

preprint2012arXiv

The Growth of the Stellar Seeds of Supermassive Black Holes

The collapse of baryons into extremely massive stars with masses exceeding 10^4 M_Sun in a small fraction of protogalaxies at z > 10 is a promising candidate for the origin of supermassive black holes, some of which grow to a billion solar masses by z ~ 7. We determine the maximum masses such stars can attain by accreting primordial gas. We find that at relatively low accretion rates the strong ionizing radiation of these stars limits their masses to M_* ~ 10^3 M_Sun (dM_acc/dt / 10^-3 M_Sun yr^-1)^8/7, where dM_acc/dt is the rate at which the star gains mass. However, at the higher central infall rates usually found in numerical simulations of protogalactic collapse (>~ 0.1 M_Sun yr^-1), the lifetime of the star instead limits its final mass to >~ 10^6 M_Sun. Furthermore, for the spherical accretion rates at which the star can grow, its ionizing radiation is confined deep within the protogalaxy, so the evolution of the star is decoupled from that of its host galaxy. Lyman alpha emission from the surrounding H II region is trapped in these heavy accretion flows and likely reprocessed into strong Balmer series emission, which may be observable by the James Webb Space Telescope. This, along with strong He II 1640 Angstrom and continuum emission, are likely to be the key observational signatures of the progenitors of supermassive black holes at high redshift.

preprint2011arXiv

Modeling the Multi-Wavelength Emission of Shell-Type Supernova Remnant RX J1713.7-3946

Emission mechanisms of the shell-type supernova remnant (SNR) RX J1713.7-3946 are studied with multi-wavelength observational data from radio, X-ray, GeV $γ$-ray to TeV $γ$-ray band. A Markov Chain Monte Carlo method is employed to explore the high-dimensional model parameter space systematically. Three scenarios for the $γ$-ray emission are investigated: the leptonic, the hadronic and a hybrid one. Thermal emission from the background plasma is also included to constrain the gas density, assuming ionization equilibrium, and a 2$σ$ upper limit of about 0.03 cm$^{-3}$ is obtained as far as thermal energies account for a significant fraction of the dissipated kinetic energy of the SNR shock. Although systematic errors dominate the $χ^2$ of the spectral fit of all models, we find that 1) the leptonic model has the best constrained model parameters, whose values can be easily accommodated with a typical supernova, but gives relatively poor fit to the TeV $γ$-ray data; 2) The hybrid scenario has one more parameter than the leptonic one and improves the overall spectral fit significantly; 3) The hadronic one, which has three more parameters than the leptonic model, gives the best fit to the overall spectrum with relatively not-well-constrained model parameters and very hard spectra of accelerated particles. The uncertainties of the model parameters decrease significantly if the spectral indices of accelerated electrons and protons are the same. The hybrid and hadronic models also require an energy input into high-energy protons, which seems to be too high compared with typical values of a supernova explosion. Further investigations are required to reconcile these observations with SNR theories.

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

Polarized Emission of Sagittarius A*

We explore the parameter space of the two temperature pseudo-Newtonian Keplerian accretion flow model for the millimeter and shorter wavelength emission from Sagittarius A*. A general relativistic ray-tracing code is used to treat the radiative transfer of polarized synchrotron emission from the flow. The synchrotron self-Comptonization and bremsstrahlung emission components are also included. It is shown that the model can readily account for the millimeter to sub-millimeter emission characteristics with an accretion rate of ~6x10^17g.s^-1 and an inclination angle of ~40 deg. However, the corresponding model predicted near-infrared and X-ray fluxes are more than one order of magnitude lower than the observed 'quiescent' state values. While the extended quiescent-state X-ray emission has been attributed to thermal emission from the large-scale accretion flow, the NIR emission and flares are likely dominated by emission regions either within the last stable orbit of a Schwarzschild black hole or associated with outflows. With the viscous parameter derived from numerical simulations, there is still a degeneracy between the electron heating rate and the magnetic parameter. A fully general relativistic treatment with the black hole spin incorporated will resolve these issues.

preprint2008arXiv

Stochastic Electron Acceleration in Shell-Type Supernova Remnants

We study the stochastic electron acceleration by fast mode waves in the turbulent downstream of weakly magnetized collisionless astrophysical shocks. The acceleration is most efficient in a dissipative layer, and the model characteristics are determined by the shock speed, density, magnetic field, and turbulence decay length. The model explains observations of shell-type supernova remnants RX J1713.7-3946 and J0852.0-4622 and can be tested by observations in hard X-rays with the HXMT and NuSTAR or gamma-rays with the GLAST.

preprint2006arXiv

A Radial Velocity Survey of the Cygnus OB2 Association

We conducted a radial velocity survey of the Cygnus OB2 Association over a 6 year (1999 - 2005) time interval to search for massive close binaries. During this time we obtained 1139 spectra on 146 OB stars to measure mean systemic radial velocities and radial velocity variations. We spectroscopically identify 73 new OB stars for the first time, the majority of which are likely to be Association members. Spectroscopic evidence is also presented for a B3Iae classification and temperature class variation (B3 - B8) on the order of 1 year for Cygnus OB2 No. 12. Calculations of the intial mass function with the current spectroscopic sample yield Gamma = -2.2 +/- 0.1. Of the 120 stars with the most reliable data, 36 are probable and 9 are possible single-lined spectroscopic binaries. We also identify 3 new and 8 candidate double-lined spectroscopic binaries. These data imply a lower limit on the massive binary fraction of 30% - 42%. The calculated velocity dispersion for Cygnus OB2 is 2.44 +/- km/s, which is typical of open clusters. No runaway OB stars were found.

preprint2006arXiv

Correlation between Flux and Spectral Index during Flares in Sagittarius A*

Flares in Sagittarius A* are produced by hot plasmas within a few Schwarzschild radii of the supermassive black hole at the Galactic center. The recent detection of a correlation between the spectral index and flux during a near infrared (NIR) flare provides a means to conduct detailed investigations of the plasma heating and radiation processes. We study the evolution of the electron distribution function under the influence of a turbulent magnetic field in a hot collisionless plasma. The magnetic field, presumably generated through instabilities in the accretion flow, can both heat the plasma via resonant wave-particle coupling and cool the electrons via radiation. The electron distribution can generally be approximated as relativistic Maxwellian. To account for the observed correlation, we find that the magnetic field needs to be anti-correlated with the electron ''temperature''. NIR and X-ray light curves are produced for a cooling and a heating phase. The model predicts simultaneous flare activity in the NIR and X-ray bands, which can be compared with observations. These results can be applied to MHD simulations to study the radiative characteristics of collisionless plasmas, especially accretion flows in low-luminosity AGNs.

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.