Showing 1–50 of 78 results for author: Matzner, C

Playing with FIRE: A Galactic Feedback-Halting Experiment Challenges Star Formation Rate Theories

Authors: Shivan Khullar, Christopher D. Matzner, Norman Murray, Michael Y. Grudić, Dávid Guszejnov, Andrew Wetzel, Philip F. Hopkins

Abstract: Stellar feedback influences the star formation rate (SFR) and the interstellar medium of galaxies in ways that are difficult to quantify numerically, because feedback is an essential ingredient of realistic simulations. To overcome this, we conduct a feedback-halting experiment starting with a Milky Way-mass galaxy in the FIRE-2 simulation framework. Terminating feedback, and comparing to a simula… ▽ More Stellar feedback influences the star formation rate (SFR) and the interstellar medium of galaxies in ways that are difficult to quantify numerically, because feedback is an essential ingredient of realistic simulations. To overcome this, we conduct a feedback-halting experiment starting with a Milky Way-mass galaxy in the FIRE-2 simulation framework. Terminating feedback, and comparing to a simulation in which feedback is maintained, we monitor how the runs diverge. We find that without feedback, interstellar turbulent velocities decay. There is a marked increase of dense material, while the SFR increases by over an order of magnitude. Importantly, this SFR boost is a factor of $\sim$15-20 larger than is accounted for by the increased free fall rate caused by higher densities. This implies that feedback moderates the star formation efficiency per free-fall time more directly than simply through the density distribution. To probe changes at the scale of giant molecular clouds (GMCs), we identify GMCs using density and virial parameter thresholds, tracking clouds as the galaxy evolves. Halting feedback stimulates rapid changes, including a proliferation of new bound clouds, a decrease of turbulent support in loosely-bound clouds, an overall increase in cloud densities, and a surge of internal star formation. Computing the cloud-integrated SFR using several theories of turbulence regulation, we show that these theories underpredict the surge in SFR by at least a factor of three. We conclude that galactic star formation is essentially feedback-regulated on scales that include GMCs, and that stellar feedback affects GMCs in multiple ways. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 20 pages, 10 figures, under ongoing review at ApJ

Electromagnetic Signatures of Mirror Stars

Authors: Isabella Armstrong, Berkin Gurbuz, David Curtin, Christopher Matzner

Abstract: Mirror Stars are a generic prediction of dissipative dark matter models, including minimal atomic dark matter and twin baryons in the Mirror Twin Higgs. Mirror Stars can capture regular matter from the interstellar medium through extremely suppressed kinetic mixing interactions between the regular and the dark photon. This accumulated "nugget" will draw heat from the mirror star core and emit high… ▽ More Mirror Stars are a generic prediction of dissipative dark matter models, including minimal atomic dark matter and twin baryons in the Mirror Twin Higgs. Mirror Stars can capture regular matter from the interstellar medium through extremely suppressed kinetic mixing interactions between the regular and the dark photon. This accumulated "nugget" will draw heat from the mirror star core and emit highly characteristic X-ray and optical signals. In this work, we devise a general parameterization of mirror star nugget properties that is independent of the unknown details of mirror star stellar physics, and use the Cloudy spectral synthesis code to obtain realistic and comprehensive predictions for the thermal emissions from optically thin mirror star nuggets. We find that mirror star nuggets populate an extremely well-defined and narrow region of the HR diagram that only partially overlaps with the white dwarf population. Our detailed spectral predictions, which we make publicly available, allow us to demonstrate that optically thin nuggets can be clearly distinguished from white dwarf stars by their continuum spectrum shape, and from planetary nebulae and other optically thin standard sources by their highly exotic emission line ratios. Our work will enable realistic mirror star telescope searches, which may reveal the detailed nature of dark matter. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 12 pages, 6 figures + references

Misaligned precessing jets are choked by the accretion disk wind

Authors: Wenbin Lu, Tatsuya Matsumoto, Christopher D. Matzner

Abstract: We analytically and numerically study the hydrodynamic propagation of a precessing jet in the context of tidal disruption events (TDEs) where the star's angular momentum is misaligned with the black hole spin. We assume that a geometrically thick accretion disk undergoes Lense-Thirring precession around the black hole spin axis and that the jet is aligned with the instantaneous disk angular moment… ▽ More We analytically and numerically study the hydrodynamic propagation of a precessing jet in the context of tidal disruption events (TDEs) where the star's angular momentum is misaligned with the black hole spin. We assume that a geometrically thick accretion disk undergoes Lense-Thirring precession around the black hole spin axis and that the jet is aligned with the instantaneous disk angular momentum. At large spin-orbit misalignment angles, the duty cycle along a given direction that the jet sweeps across is much smaller than unity. The faster jet and slower disk wind alternately fill a given angular region, which leads to strong shock interactions between the two. We show that precessing jets can only break out of the wind confinement if the misalignment angle is less than a few times the jet opening angle. The very small event rate of observed jetted TDEs is then explained by the condition of double alignment: both the stellar angular momentum and the observer's line of sight are nearly aligned with the black hole spin. In most TDEs, the jets are initially choked by the disk wind and may only break out later when the disk eventually aligns itself with the spin axis due to the viscous damping of the precession. Such late-time jets may produce delayed radio rebrightening as seen in many optically selected TDEs. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 11 pages, 6 figures, plus appendix. MNRAS submitted. Comments welcome!

Nonlinear acoustics and shock dynamics in isentropic atmospheres

Authors: Tamar Faran, Christopher D. Matzner, Eliot Quataert

Abstract: Nonlinear acoustic evolution is often discussed in the context of wave-steepening that leads to shock formation, and is of special interest in applications where the shock continues to strengthen due to a narrowing of its channel or the stratification of the medium. Accurate scalings govern low amplitude waves and strong shocks, but connecting these phases, or describing waves that are nonlinear f… ▽ More Nonlinear acoustic evolution is often discussed in the context of wave-steepening that leads to shock formation, and is of special interest in applications where the shock continues to strengthen due to a narrowing of its channel or the stratification of the medium. Accurate scalings govern low amplitude waves and strong shocks, but connecting these phases, or describing waves that are nonlinear from the outset, generally requires simulation. We address this problem using the fact that waves within a plane-parallel, isentropic and gravitationally stratified atmosphere are described by exact simple-wave solutions, thanks to the conservation of Riemann invariants in a freely falling reference frame. Our solutions enable us to discriminate waves that reflect from those that form shocks, and to capture wave and shock evolution using an ordinary differential equation. For several relevant values of the adiabatic index $γ$ the solutions are explicit; furthermore, nonlinear wave reflection from a free surface can be described analytically for $γ=3$. Comparison to hydrodynamic simulations shows that our analytic shock approximation is accurate up to moderate ($\sim$ few--15) Mach numbers, where the accuracy increases with the adiabatic index. Our solutions also imply that an initially subsonic pulse is unable to unbind mass from the atmosphere without significantly increasing its entropy. △ Less

Submitted 23 December, 2023; v1 submitted 2 October, 2023; originally announced October 2023.

Comments: Paper withdrawn due to a conceptual error regarding backward travelling perturbations (section 3). We are working on an update

Observability of flashes from ejecta crashes in aspherical supernovae, with application to SN 2008D

Authors: Benjamin Scully, Christopher D. Matzner, Almog Yalinewich

Abstract: A new class of transient, which has been hypothesized to accompany the explosion of an aspherical compact supernova, would arise when streams of ejecta collide outside the star. However, conditions that favour the prompt release of radiation from the collision, such as a diffuse stellar envelope, disfavour the creation of non-radial ejecta in the first place. To determine whether the collision can… ▽ More A new class of transient, which has been hypothesized to accompany the explosion of an aspherical compact supernova, would arise when streams of ejecta collide outside the star. However, conditions that favour the prompt release of radiation from the collision, such as a diffuse stellar envelope, disfavour the creation of non-radial ejecta in the first place. To determine whether the collision can both occur and be visible, we simulate aspherical explosions using the HUJI-RICH moving-mesh hydrodynamics code and analyze them in terms of diffusion measures defined for individual fluid elements. While our simulations are highly idealized, they connect to realistic explosions via a single dimensionless parameter. Defining two measures of the importance of diffusivity (two versions of the inverse Péclet number), we find that one varies in a way that indicates colliding ejecta can release a photon flash, while the other does not. Examining the x-ray transient XT 080109 associated with supernova SN 2008D, we find that its fluence and duration are consistent with the properties of an ejecta collision in the aspherical model that is most likely to emit a flash. Our results give tentative evidence for the possibility of collision-induced flashes for a narrow and radius-dependent range of asphericity, and motivate future radiation hydrodynamics simulations. △ Less

Submitted 17 August, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 8 pages, 3 figures, accepted for publication in MNRAS, references updated

Origin of high-velocity ejecta and early red excess emission in the infant Type Ia supernova 2021aefx

Authors: Yuan Qi Ni, Dae-Sik Moon, Maria R. Drout, Christopher D. Matzner, Kelvin C. C. Leong, Sang Chul Kim, Hong Soo Park, Youngdae Lee

Abstract: \object{SN 2021aefx} is a normal Type Ia Supernova (SN) with red excess emission over the first $\sim$ 2 days. We present detailed analysis of this SN using our high-cadence KMTNet multi-band photometry, spectroscopy, and publicly available data. We provide the first measurements of its epochs of explosion (MJD 59529.32 $\pm$ 0.16) as well as ``first light'' (MJD 59529.85 $\pm$ 0.55) associated wi… ▽ More \object{SN 2021aefx} is a normal Type Ia Supernova (SN) with red excess emission over the first $\sim$ 2 days. We present detailed analysis of this SN using our high-cadence KMTNet multi-band photometry, spectroscopy, and publicly available data. We provide the first measurements of its epochs of explosion (MJD 59529.32 $\pm$ 0.16) as well as ``first light'' (MJD 59529.85 $\pm$ 0.55) associated with the main ejecta ${\rm{^{56}Ni}}$ distribution. This places our first detection of SN 2021aefx at $\sim -$0.5 hours since ``first light'', indicating the presence of additional power sources. Our peak-spectrum confirms its Type Ia sub-classification as intermediate between Core-Normal and Broad-Line, and we estimate the ejecta mass to be $\sim$ 1.34 $M_{\odot}$. The pre-peak spectral evolution identifies fast-expanding material reaching $>$ 40,000 km s$^{-1}$ (the fastest ever observed in Type Ia SNe) and at least two distinct homologously-expanding ejecta components: (1) a normal-velocity (12,400 km s$^{-1}$) component consistent with the typical photospheric evolution of Chandrasekhar-mass ejecta; and (2) a high-velocity (23,500 km s$^{-1}$) component visible during the first $\sim$ 3.6 days post-explosion, which locates the component within the outer $<$ 16\% of the ejecta mass. Asymmetric, subsonic explosion processes producing a non-spherical photosphere provide an explanation for the simultaneous presence of the two components, as well as the red excess emission via a slight ${\rm{^{56}Ni}}$ enrichment in the outer $\sim$ 0.5\% of the ejecta mass. Our spectrum from 300 days post-peak advances the constraint against non-degenerate companions and further supports a near-Chandrasekhar-mass explosion origin. Off-center ignited delayed-detonations of Chandrasekhar-mass white dwarfs may be responsible for the observed features of SN 2021aefx in some normal Type Ia SNe. △ Less

Submitted 5 April, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: Submitted for publication in ApJ. 29 pages, 14 figures, 4 tables

Velocity-Coherent Substructure in TMC-1: Inflow and Fragmentation

Authors: Simon E. T. Smith, Rachel Friesen, Antoine Marchal, Jaime E. Pineda, Paola Caselli, Michael Chun-Yuan Chen, Spandan Choudhury, James Di Francesco, Adam Ginsburg, Helen Kirk, Chris Matzner, Anna Punanova, Samantha Scibelli, Yancy Shirley

Abstract: Filamentary structures have been found nearly ubiquitously in molecular clouds and yet their formation and evolution is still poorly understood. We examine a segment of Taurus Molecular Cloud 1 (TMC-1) that appears as a single, narrow filament in continuum emission from dust. We use the Regularized Optimization for Hyper-Spectral Analysis (ROHSA), a Gaussian decomposition algorithm which enforces… ▽ More Filamentary structures have been found nearly ubiquitously in molecular clouds and yet their formation and evolution is still poorly understood. We examine a segment of Taurus Molecular Cloud 1 (TMC-1) that appears as a single, narrow filament in continuum emission from dust. We use the Regularized Optimization for Hyper-Spectral Analysis (ROHSA), a Gaussian decomposition algorithm which enforces spatial coherence when fitting multiple velocity components simultaneously over a data cube. We analyze HC$_5$N (9-8) line emission as part of the Green Bank Ammonia Survey (GAS) and identify three velocity-coherent components with ROHSA. The two brightest components extend the length of the filament, while the third component is fainter and clumpier. The brightest component has a prominent transverse velocity gradient of $2.7 \pm 0.1$ km s$^{-1}$ pc$^{-1}$ that we show to be indicative of gravitationally induced inflow. In the second component, we identify regularly spaced emission peaks along its length. We show that the local minima between pairs of adjacent HC$_5$N peaks line up closely with submillimetre continuum emission peaks, which we argue is evidence for fragmentation along the spine of TMC-1. While coherent velocity components have been described as separate physical structures in other star-forming filaments, we argue that the two bright components identified in HC$_5$N emission in TMC-1 are tracing two layers in one filament: a lower density outer layer whose material is flowing under gravity towards the higher density inner layer of the filament. △ Less

Submitted 6 February, 2023; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: 17 pages, 8 figures; Accepted for publication to MNRAS

A low-cost ultraviolet-to-infrared absolute quantum efficiency characterization system of detectors

Authors: Ajay S. Gill, Mohamed M. Shaaban, Aaron Tohuvavohu, Suresh Sivanandam, Roberto G. Abraham, Seery Chen, Maria R. Drout, Deborah Lokhorst, Christopher D. Matzner, Stefan W. Mochnacki, Calvin B. Netterfield

Abstract: We present a low-cost ultraviolet to infrared absolute quantum efficiency detector characterization system developed using commercial off-the-shelf components. The key components of the experiment include a light source,a regulated power supply, a monochromator, an integrating sphere, and a calibrated photodiode. We provide a step-by-step procedure to construct the photon and quantum efficiency tr… ▽ More We present a low-cost ultraviolet to infrared absolute quantum efficiency detector characterization system developed using commercial off-the-shelf components. The key components of the experiment include a light source,a regulated power supply, a monochromator, an integrating sphere, and a calibrated photodiode. We provide a step-by-step procedure to construct the photon and quantum efficiency transfer curves of imaging sensors. We present results for the GSENSE 2020 BSI CMOS sensor and the Sony IMX 455 BSI CMOS sensor. As a reference for similar characterizations, we provide a list of parts and associated costs along with images of our setup. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: 15 pages, 14 figures, SPIE Astronomical Telescopes and Instrumentation, Conference 12191: X-ray, Optical, and Infrared Detectors for Astronomy X, Paper Number: 12191-39 (Montreal, July 2022)

The origin and evolution of the normal Type Ia SN 2018aoz with infant-phase reddening and excess emission

Authors: Yuan Qi Ni, Dae-Sik Moon, Maria R. Drout, Abigail Polin, David J. Sand, Santiago GonzÁlez-GaitÁn, Sang Chul Kim, Youngdae Lee, Hong Soo Park, D. Andrew Howell, Peter E. Nugent, Anthony L. Piro, Peter J. Brown, LluÍs Galbany, Jamison Burke, Daichi Hiramatsu, Griffin Hosseinzadeh, Stefano Valenti, Niloufar Afsariardchi, Jennifer E. Andrews, John Antoniadis, Rachael L. Beaton, K. Azalee Bostroem, Raymond G. Carlberg, S. Bradley Cenko , et al. (18 additional authors not shown)

Abstract: SN~2018aoz is a Type Ia SN with a $B$-band plateau and excess emission in the infant-phase light curves $\lesssim$ 1 day after first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN~2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show the SN is int… ▽ More SN~2018aoz is a Type Ia SN with a $B$-band plateau and excess emission in the infant-phase light curves $\lesssim$ 1 day after first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN~2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show the SN is intermediate between two subtypes of normal Type Ia: Core-Normal and Broad-Line. The excess emission could have contributions from the radioactive decay of surface iron-peak elements as well as ejecta interaction with either the binary companion or a small torus of circumstellar material. Nebular-phase limits on H$α$ and He~I favour a white dwarf companion, consistent with the small companion size constrained by the low early SN luminosity, while the absence of [O~I] and He~I disfavours a violent merger of the progenitor. Of the two main explosion mechanisms proposed to explain the distribution of surface iron-peak elements in SN~2018aoz, the asymmetric Chandrasekhar-mass explosion is less consistent with the progenitor constraints and the observed blueshifts of nebular-phase [Fe~II] and [Ni~II]. The helium-shell double-detonation explosion is compatible with the observed lack of C spectral features, but current 1-D models are incompatible with the infant-phase excess emission, $B_{\rm max}-V_{\rm max}$ color, and absence of nebular-phase [Ca~II]. Although the explosion processes of SN~2018aoz still need to be more precisely understood, the same processes could produce a significant fraction of Type Ia SNe that appear normal after $\sim$ 1 day. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: Submitted for publication in ApJ. 35 pages, 16 figures, 7 tables

Infant-phase reddening by surface Fe-peak elements in a normal Type Ia Supernova

Authors: Yuan Qi Ni, Dae-Sik Moon, Maria R. Drout, Abigail Polin, David J. Sand, Santiago Gonzalez-Gaitan, Sang Chul Kim, Youngdae Lee, Hong Soo Park, D. Andrew Howell, Peter E. Nugent, Anthony L. Piro, Peter J. Brown, Lluis Galbany, Jamison Burke, Daichi Hiramatsu, Griffin Hosseinzadeh, Stefano Valenti, Niloufar Afsariardchi, Jennifer E. Andrews, John Antoniadis, Iair Arcavi, Rachael L. Beaton, K. Azalee Bostroem, Raymond G. Carlberg , et al. (19 additional authors not shown)

Abstract: Type Ia Supernovae are thermonuclear explosions of white dwarf stars. They play a central role in the chemical evolution of the Universe and are an important measure of cosmological distances. However, outstanding questions remain about their origins. Despite extensive efforts to obtain natal information from their earliest signals, observations have thus far failed to identify how the majority of… ▽ More Type Ia Supernovae are thermonuclear explosions of white dwarf stars. They play a central role in the chemical evolution of the Universe and are an important measure of cosmological distances. However, outstanding questions remain about their origins. Despite extensive efforts to obtain natal information from their earliest signals, observations have thus far failed to identify how the majority of them explode. Here, we present infant-phase detections of SN 2018aoz from a brightness of -10.5 absolute AB magnitudes -- the lowest luminosity early Type Ia signals ever detected -- revealing a hitherto unseen plateau in the $B$-band that results in a rapid redward color evolution between 1.0 and 12.4 hours after the estimated epoch of first light. The missing $B$-band flux is best-explained by line-blanket absorption from Fe-peak elements in the outer 1% of the ejected mass. The observed $B-V$ color evolution of the SN also matches the prediction from an over-density of Fe-peak elements in the same outer 1% of the ejected mass, whereas bluer colors are expected from a purely monotonic distribution of Fe-peak elements. The presence of excess nucleosynthetic material in the extreme outer layers of the ejecta points to enhanced surface nuclear burning or extended sub-sonic mixing processes in some normal Type Ia Supernova explosions. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: Accepted for publication in Nature Astronomy. Main text = 8 pages, 4 figures, 1 table; Full document = 46 pages, with Methods, Supplementary Information, 7 Supplementary figures, 2 Supplementary tables and references. Nat Astron (2022)

How to search for Mirror Stars with Gaia

Authors: Aaron Howe, Jack Setford, David Curtin, Christopher D. Matzner

Abstract: We show for the first time how to conduct a direct search for dark matter using Gaia observations. Its public astrometric data may contain the signals of mirror stars, exotic compact objects made of atomic dark matter with a tiny kinetic mixing between the dark and SM photon. Mirror stars capture small amounts of interstellar material in their cores, leading to characteristic optical/IR and X-ray… ▽ More We show for the first time how to conduct a direct search for dark matter using Gaia observations. Its public astrometric data may contain the signals of mirror stars, exotic compact objects made of atomic dark matter with a tiny kinetic mixing between the dark and SM photon. Mirror stars capture small amounts of interstellar material in their cores, leading to characteristic optical/IR and X-ray emissions. We develop the detailed pipeline for conducting a mirror star search using data from Gaia and other stellar catalogues, and demonstrate our methodology by conducting a search for toy mirror stars with a simplified calculation of their optical/IR emissions over a wide range of mirror star and hidden sector parameters. We also obtain projected exclusion bounds on the abundance and properties of mirror stars if no candidates are found, demonstrating that Gaia is a new and uniquely powerful probe of atomic dark matter. Our study provides the blueprint for a realistic mirror star search that includes a more complete treatment of the captured interstellar gas in the future. △ Less

Submitted 28 June, 2022; v1 submitted 10 December, 2021; originally announced December 2021.

Comments: 19 pages, 13 figures. Accepted for publication in JHEP; added various clarifications, references, typo fixes; no change in results or conclusions

Are massive dense clumps truly sub-virial? A new analysis using Gould Belt ammonia data

Authors: Ayushi Singh, Christopher D. Matzner, Rachel K. Friesen, Peter G. Martin, Jaime E. Pineda, Erik W. Rosolowsky, Felipe Alves, Ana Chacón-Tanarro, Hope How-Huan Chen, Michael Chun-Yuan Chen, Spandan Choudhury, James Di Francesco, Jared Keown, Helen Kirk, Anna Punanova, Youngmin Seo, Yancy Shirley, Adam Ginsburg, Stella S. R. Offner, Héctor G. Arce, Paola Caselli, Alyssa A. Goodman, Philip C. Myers, Elena Redaelli

Abstract: Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of… ▽ More Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of this conclusion, we use ammonia line data from the Green Bank Ammonia Survey and Planck-calibrated dust emission data from Herschel to estimate the masses and kinetic and gravitational energies for dense clumps in the Gould Belt clouds. We show that several types of systematic error can enhance the appearance of low kinetic-to-gravitational energy ratios: insufficient removal of foreground and background material; ignoring the kinetic energy associated with velocity differences across a resolved cloud; and over-correcting for stratification when evaluating the gravitational energy. Using an analysis designed to avoid these errors, we find that the most massive Gould Belt clumps harbor virial motions, rather than sub-virial ones. As a byproduct, we present a catalog of masses, energies, and virial energy ratios for 85 Gould Belt clumps. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: Submitted to ApJ

The density structure of supersonic self-gravitating turbulence

Authors: Shivan Khullar, Christoph Federrath, Mark R. Krumholz, Christopher D. Matzner

Abstract: We conduct numerical experiments to determine the density probability distribution function (PDF) produced in supersonic, isothermal, self-gravitating turbulence of the sort that is ubiquitous in star-forming molecular clouds. Our experiments cover a wide range of turbulent Mach number and virial parameter, allowing us for the first time to determine how the PDF responds as these parameters vary,… ▽ More We conduct numerical experiments to determine the density probability distribution function (PDF) produced in supersonic, isothermal, self-gravitating turbulence of the sort that is ubiquitous in star-forming molecular clouds. Our experiments cover a wide range of turbulent Mach number and virial parameter, allowing us for the first time to determine how the PDF responds as these parameters vary, and we introduce a new diagnostic, the dimensionless star formation efficiency versus density ($ε_{\rm ff}(s)$) curve, which provides a sensitive diagnostic of the PDF shape and dynamics. We show that the PDF follows a universal functional form consisting of a log-normal at low density with two distinct power law tails at higher density; the first of these represents the onset of self-gravitation, and the second reflects the onset of rotational support. Once the star formation efficiency reaches a few percent, the PDF becomes statistically steady, with no evidence for secular time-evolution at star formation efficiencies from about five to 20 percent. We show that both the Mach number and the virial parameter influence the characteristic densities at which the log-normal gives way to the first power-law, and the first to the second, and we extend (for the former) and develop (for the latter) simple theoretical models for the relationship between these density thresholds and the global properties of the turbulent medium. △ Less

Submitted 7 August, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: 18 pages, 10 figures, accepted in MNRAS, revised version

Wave-Driven Shocks in Stellar Outbursts: Dynamics, Envelope Heating, and Nascent Blastwaves

Authors: Christopher D. Matzner, Stephen Ro

Abstract: We address the shocks from acoustic pulses and wave trains in general one-dimensional flows, with an emphasis on the application to super-Eddington outbursts in massive stars. Using approximate adiabatic invariants, we generalize the classical equal-area technique in its integral and differential forms. We predict shock evolution for the case of an initially sinusoidal but finite wave train, with… ▽ More We address the shocks from acoustic pulses and wave trains in general one-dimensional flows, with an emphasis on the application to super-Eddington outbursts in massive stars. Using approximate adiabatic invariants, we generalize the classical equal-area technique in its integral and differential forms. We predict shock evolution for the case of an initially sinusoidal but finite wave train, with separate solutions for internal shocks and head or tail shocks, and demonstrate detailed agreement with numerical simulations. Our internal shock solution motivates improved expressions for the shock heating rate. Our solution for head and tail shocks demonstrates that these preserve dramatically more wave energy to large radii and have a greater potential for the direct ejection of matter. This difference highlights the importance of the waveform for shock dynamics. Our weak-shock analysis predicts when shocks will become strong and provides a basis from which this transition can be addressed. We use it to estimate the mass ejected by sudden sound pulses and weak central explosions. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: 11 pages, 5 figures. ApJ, accepted

On Linear and Nonlinear Acoustics in Stratified, Variable-Area Ducts and Atmospheres, and Lighthill's Proposition

Authors: Christopher D. Matzner, Stephen Ro

Abstract: We consider linear and nonlinear waves in a stratified hydrostatic fluid within a channel of variable area, under the restriction of one-dimensional flow. We derive a modified version of Riemann's invariant that is related to the wave luminosity. This quantity obeys a simple dynamical equation in linear theory, from which the rules of wave reflection are easily discerned; and it is adiabatically c… ▽ More We consider linear and nonlinear waves in a stratified hydrostatic fluid within a channel of variable area, under the restriction of one-dimensional flow. We derive a modified version of Riemann's invariant that is related to the wave luminosity. This quantity obeys a simple dynamical equation in linear theory, from which the rules of wave reflection are easily discerned; and it is adiabatically conserved in the high-frequency limit. Following a suggestion by Lighthill, we apply the linear adiabatic invariant to predict mildly nonlinear waves. This incurs only moderate error. We find that Lighthill's criterion for shock formation is essentially exact for leading shocks, and for shocks within high-frequency waves. We conclude that approximate invariants can be used to accurately predict the self-distortion of low-amplitude acoustic pulses, as well as the dissipation patterns of weak shocks, in complicated environments such as stellar envelopes. We also identify fully nonlinear solutions for a restricted class of problems. △ Less

Submitted 13 November, 2020; originally announced November 2020.

Comments: 15 pages, 4 figures, submitted to J. Fluid Mech

The Nickel Mass Distribution of Stripped-Envelope Supernovae: Implications for Additional Power Sources

Authors: Niloufar Afsariardchi, Maria R. Drout, David Khatami, Christopher D. Matzner, Dae-Sik Moon, Yuan Qi Ni

Abstract: We perform a systematic study of the $^{56}$Ni mass ($M_{\rm Ni}$) of 27 stripped envelope supernovae (SESNe) by modeling their light-curve tails, highlighting that use of ``Arnett's rule'' overestimates $M_{\rm Ni}$ for SESN by a factor of $\sim$2. Recently, \citet{Khatami2019} presented a new model relating the peak time ($t_{\rm p}$) and luminosity ($L_{\rm p}$) of a radioactive-powered SN to i… ▽ More We perform a systematic study of the $^{56}$Ni mass ($M_{\rm Ni}$) of 27 stripped envelope supernovae (SESNe) by modeling their light-curve tails, highlighting that use of ``Arnett's rule'' overestimates $M_{\rm Ni}$ for SESN by a factor of $\sim$2. Recently, \citet{Khatami2019} presented a new model relating the peak time ($t_{\rm p}$) and luminosity ($L_{\rm p}$) of a radioactive-powered SN to its $M_{\rm Ni}$ that addresses several limitations of Arnett-like models, but depends on a dimensionless parameter, $β$. Using observed $t_{\rm p}$, $L_{\rm p}$, and tail-measured $M_{\rm Ni}$ values for 27 SESN, we observationally calibrate $β$ for the first time. Despite scatter, we demonstrate that the model of \citet{Khatami2019} with empirically-calibrated $β$ values provides significantly improved measurements of $M_{\rm Ni}$ when only photospheric data is available. However, these observationally-constrained $β$ values are systematically lower than those inferred from numerical simulations, primarily because the observed sample has significantly higher (0.2-0.4 dex) $L_{\rm p}$ for a given $M_{\rm Ni}$. While effects due to composition, mixing, and asymmetry can increase $L_{\rm p}$ current models cannot explain the systematically low $β$ values. However, the discrepancy can be alleviated if $\sim$7--50\% of $L_{\rm p}$ for the observed sample originates from sources other than $^{56}$Ni. Either shock cooling or magnetar spin-down could provide the requisite luminosity. Finally, we find that even with our improved measurements, the $M_{\rm Ni}$ values of SESN are still a factor of $\sim$3 larger than those of hydrogen-rich Type II SN, indicating that these supernovae are inherently different in terms of their progenitor initial mass distributions or explosion mechanisms. △ Less

Submitted 15 June, 2021; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: 27 pages, 12 figures, accepted for publication in ApJ

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

Authors: Che-Yu Chen, Erica A. Behrens, Jasmin E. Washington, Laura M. Fissel, Rachel K. Friesen, Zhi-Yun Li, Jaime E. Pineda, Adam Ginsburg, Helen Kirk, Samantha Scibelli, Felipe Alves, Elena Redaelli, Paola Caselli, Anna Punanova, James Di Francesco, Erik Rosolowsky, Stella S. R. Offner, Peter G. Martin, Ana Chacón-Tanarro, Hope H. -H. Chen, Michael C. -Y. Chen, Jared Keown, Youngmin Seo, Yancy Shirley, Hector G. Arce , et al. (4 additional authors not shown)

Abstract: The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular cl… ▽ More The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud. △ Less

Submitted 24 March, 2020; originally announced March 2020.

Comments: 18 pages, 11 figures, accepted for publication in MNRAS

Velocity-coherent Filaments in NGC 1333: Evidence for Accretion Flow?

Authors: Michael Chun-Yuan Chen, James Di Francesco, Erik Rosolowsky, Jared Keown, Jaime E. Pineda, Rachel K. Friesen, Paola Caselli, How-Huan Chen, Christopher D. Matzner, Stella S. Offner, Anna Punanova, Elena Redaelli, Samantha Scibelli, Yancy Shirley

Abstract: Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficie… ▽ More Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficient, robust, and automatic method to fit ammonia lines with multiple velocity components, generalizable to other molecular species. We also present CRISPy, a Python package to identify filament spines in 3D images (e.g., position-position-velocity cubes), along with a complementary technique to sort fitted velocity components into velocity-coherent filaments. In NGC 1333, we find a wealth of velocity gradient structures on a beam-resolved scale of ~0.05 pc. Interestingly, these local velocity gradients are not randomly oriented with respect to filament spines and their perpendicular, i.e., radial, component decreases in magnitude towards the spine for many filaments. Together with remarkably constant velocity gradients on larger scales along many filaments, these results suggest a scenario in which gas falling onto filaments is progressively damped and redirected to flow along these filaments. △ Less

Submitted 26 February, 2020; originally announced February 2020.

Comments: Accepted to ApJ

Journal ref: ApJ 891 (2020) 84

Simulation of a Compact Object with Outflows Moving Through a Gaseous Background

Authors: Xinyu Li, Philip Chang, Yuri Levin, Christopher D. Matzner, Philip J. Armitage

Abstract: A compact object moving relative to surrounding gas accretes material and perturbs the density of gas in its vicinity. In the classical picture of Bondi-Hoyle-Lyttleton accretion, the perturbation takes the form of an overdense wake behind the object, which exerts a dynamical friction drag. We use hydrodynamic simulations to investigate how the accretion rate and strength of dynamical friction are… ▽ More A compact object moving relative to surrounding gas accretes material and perturbs the density of gas in its vicinity. In the classical picture of Bondi-Hoyle-Lyttleton accretion, the perturbation takes the form of an overdense wake behind the object, which exerts a dynamical friction drag. We use hydrodynamic simulations to investigate how the accretion rate and strength of dynamical friction are modified by the presence of outflow from the compact object. We show that the destruction of the wake by an outflow reduces dynamical friction, and reverses its sign when the outflow is strong enough, in good quantitative agreement with analytic calculations. For a strong isotropic outflow, the outcome on scales that we have simulated is a negative dynamical friction, i.e., net acceleration. For jet-like outflows driven by reprocessed accretion, both the rate of accretion and the magnitude of dynamical friction drop for more powerful jets. The accretion rate is strongly intermittent when the jet points to the same direction as the motion of the compact object. The dynamical effects of outflows may be important for the evolution of compact objects during the common envelope phase of binary systems, and for accreting compact objects and massive stars encountering AGN discs. △ Less

Submitted 30 March, 2020; v1 submitted 14 December, 2019; originally announced December 2019.

Comments: Accepted by MNRAS

Constraints on the binary black hole hypothesis for system LB-1

Authors: R. -F. Shen, C. D. Matzner, A. W. Howard, W. Zhang

Abstract: At about 70 solar masses, the recently-discovered dark object orbited by a B-type star in the system LB-1 is difficult to understand as the end point of standard stellar evolution, except as a binary black hole (BBH). LB-1 shows a strong, broad H-alpha emission line that is best attributed to a gaseous disk surrounding the dark mass. We use the observed H-alpha line shape, particularly its wing ex… ▽ More At about 70 solar masses, the recently-discovered dark object orbited by a B-type star in the system LB-1 is difficult to understand as the end point of standard stellar evolution, except as a binary black hole (BBH). LB-1 shows a strong, broad H-alpha emission line that is best attributed to a gaseous disk surrounding the dark mass. We use the observed H-alpha line shape, particularly its wing extension, to constrain the inner radius of the disk and thereby the separation of a putative BBH. The hypothesis of a current BBH is effectively ruled out on the grounds that its merger time must be a small fraction of the current age of the B star. The hypothesis of a previous BBH that merged to create the current dark mass is also effectively ruled out by the low orbital eccentricity, due to the combination of mass loss and kick resulted from gravitational wave emission in any past merger. We conclude that the current dark mass is a single black hole produced by the highly mass-conserving, monolithic collapse of a massive star. △ Less

Submitted 28 November, 2019; originally announced November 2019.

Comments: 5 pages, submitted to ApJL

The Formation of Stars -- From Filaments to Cores to Protostars and Protoplanetry Disks

Authors: James Di Francesco, Helen Kirk, Doug Johnstone, Ralph Pudritz, Shantanu Basu, Sarah Sadavoy, Laura Fissel, Lewis Knee, Mehrnoosh Tahani, Rachel Friesen, Simon Coudé, Erik Rosolowsky, Nienke van der Marel, Michel Fich, Christine Wilson, Chris Matzner, Ruobing Dong, Brenda Matthews, Gerald Schieven

Abstract: Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length sc… ▽ More Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length scales leads to the formation and evolution of substructures within clouds, including filaments, clumps, cores, disks, outflows, the protostars/YSOs themselves, and planets. In this white paper, prepared for the 2020 Long Range Plan panel which will recommend Canada's future directions for astronomy, we describe the observational and theoretical leadership in the star formation field that Canada's vibrant community has demonstrated over the past decade. Drawing from this extensive background, we identify five key questions that must be addressed for further progress to be made in understanding star formation in the next decade. Addressing these questions will improve our understanding of the dynamics of the dense gas and the role of the magnetic field in star formation, the optical properties of the dust used to trace mass and magnetic fields, the sources of variability in star-forming objects on short timescales, and the physical processes that specifically promote the clustering of stars. We further highlight key facilities in which Canada should become involved to continue making progress in this field. Single-dish facilities we recommend include LSST, trans-atmospheric far-infrared telescopes like BLAST-TNG and SPICA, and ground-based telescopes like JCMT, GBT, and CCAT-p. Interferometric facilities we recommend include ALMA, ngVLA, and SKA1. △ Less

Submitted 9 November, 2019; v1 submitted 5 November, 2019; originally announced November 2019.

Comments: 11 pages, a contributed white paper prepared for Canada's 2020 Long Range Plan decadal process

Astronomy in a Low-Carbon Future

Authors: Christopher D. Matzner, Nicolas B. Cowan, René Doyon, Vincent Hénault-Brunet, David Lafrenère, Martine Lokken, Peter G. Martin, Sharon Morsink, Magdalen Normandeau, Nathalie Ouellette, Mubdi Rahman, Joel Roediger, James Taylor, Rob Thacker, Marten van Kerkwijk

Abstract: The global climate crisis poses new risks to humanity, and with them, new challenges to the practices of professional astronomy. Avoiding the more catastrophic consequences of global warming by more than 1.5 degrees requires an immediate reduction of greenhouse gas emissions. According to the 2018 United Nations Intergovernmental Panel report, this will necessitate a 45% reduction of emissions by… ▽ More The global climate crisis poses new risks to humanity, and with them, new challenges to the practices of professional astronomy. Avoiding the more catastrophic consequences of global warming by more than 1.5 degrees requires an immediate reduction of greenhouse gas emissions. According to the 2018 United Nations Intergovernmental Panel report, this will necessitate a 45% reduction of emissions by 2030 and net-zero emissions by 2050. Efforts are required at all levels, from the individual to the governmental, and every discipline must find ways to achieve these goals. This will be especially difficult for astronomy with its significant reliance on conference and research travel, among other impacts. However, our long-range planning exercises provide the means to coordinate our response on a variety of levels. We have the opportunity to lead by example, rising to the challenge rather than reacting to external constraints. We explore how astronomy can meet the challenge of a changing climate in clear and responsible ways, such as how we set expectations (for ourselves, our institutions, and our granting agencies) around scientific travel, the organization of conferences, and the design of our infrastructure. We also emphasize our role as reliable communicators of scientific information on a problem that is both human and planetary in scale. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: 16 pages. A White Paper prepared for the Canadian Long Range Plan 2020

Optical Transient from an Explosion Close to the Stellar Surface

Authors: Almog Yalinewich, Christopher D. Matzner

Abstract: We study the hydrodynamic evolution of an explosion close to the stellar surface, and give predictions for the radiation from such an event. We show that such an event will give rise to a multi-wavelength transient. We apply this model to describe a precursor burst to the peculiar supernova iPTF14hls, which occurred in 1954, sixty year before the supernova. We propose that the new generation of op… ▽ More We study the hydrodynamic evolution of an explosion close to the stellar surface, and give predictions for the radiation from such an event. We show that such an event will give rise to a multi-wavelength transient. We apply this model to describe a precursor burst to the peculiar supernova iPTF14hls, which occurred in 1954, sixty year before the supernova. We propose that the new generation of optical surveys might detect similar transients, and they can be used to identify supernova progenitors well before the explosion. △ Less

Submitted 17 July, 2019; originally announced July 2019.

KSP-SN-2016kf: a long-rising H-rich Type II Supernova with unusually high $^{56}$Ni mass discovered in the KMTNet Supernova Program

Authors: Niloufar Afsariardchi, Dae-Sik Moon, Maria R. Drout, Santiago González-Gaitán, Yuan Qi Ni, Christopher D. Matzner, Sang Chul Kim, Youngdae Lee, Hong Soo Park, Avishay Gal-Yam, Giuliano Pignata, Bon-Chul Koo, Stuart Ryder, Sang-Mok Cha, Yongseok Lee

Abstract: We present the discovery and the photometric and spectroscopic study of H-rich Type II supernova (SN) KSP-SN-2016kf (SN2017it) observed in the KMTNet Supernova Program in the outskirts of a small irregular galaxy at $z\simeq0.043$ within a day from the explosion. Our high-cadence, multi-color ($BVI$) light curves of the SN show that it has a very long rise time ($t_\text{rise}\simeq 20$ days in… ▽ More We present the discovery and the photometric and spectroscopic study of H-rich Type II supernova (SN) KSP-SN-2016kf (SN2017it) observed in the KMTNet Supernova Program in the outskirts of a small irregular galaxy at $z\simeq0.043$ within a day from the explosion. Our high-cadence, multi-color ($BVI$) light curves of the SN show that it has a very long rise time ($t_\text{rise}\simeq 20$ days in $V$ band), a moderately luminous peak ($M_V\simeq -$17.6 mag), a notably luminous and flat plateau ($M_V\simeq -$17.4 mag and decay slope $s\simeq0.53$ mag per 100 days), and an exceptionally bright radioactive tail. Using the color-dependent bolometric correction to the light curves, we estimate the $^{56}$Ni mass powering the observed radioactive tail to be $0.10\pm0.01$ M$_\odot$, making it a H-rich Type II SN with one of the largest $^{56}$Ni masses observed to date. The results of our hydrodynamic simulations of the light curves constrain the mass and radius of the progenitor at the explosion to be $\sim$15 M$_\odot$ (evolved from a star with an initial mass of $\sim$ 18.8 M$_\odot$) and $\sim1040$ R$_\odot$, respectively, with the SN explosion energy of $\sim 1.3\times 10^{51}$ erg s$^{-1}$. The above-average mass of the KSP-SN-2016kf progenitor, together with its low metallicity $ Z/Z_\odot \simeq0.1-0.4$ obtained from spectroscopic analysis, is indicative of a link between the explosion of high-mass red supergiants and their low-metallicity environment. The early part of the observed light curves shows the presence of excess emission above what is predicted in model calculations, suggesting there is interaction between the ejecta and circumstellar material. We further discuss the implications of the high progenitor initial mass and low-metallicity environment of KSP-SN-2016kf on our understanding of the origin of Type II SNe. △ Less

Submitted 21 June, 2019; originally announced June 2019.

Comments: 18 pages, 12 figures, accepted for publication in ApJ

Negative Dynamical Friction on compact objects moving through dense gas

Authors: Andrei Gruzinov, Yuri Levin, Christopher D. Matzner

Abstract: An over-dense wake is created by a gravitating object moving through a gaseous medium, and this wake pulls back on the object and slows it down. This is conventional dynamical friction in a gaseous medium. We argue that if the object drives a sufficiently powerful outflow, the wake is destroyed and instead an extended under-dense region is created behind the object. In this case the overall gravit… ▽ More An over-dense wake is created by a gravitating object moving through a gaseous medium, and this wake pulls back on the object and slows it down. This is conventional dynamical friction in a gaseous medium. We argue that if the object drives a sufficiently powerful outflow, the wake is destroyed and instead an extended under-dense region is created behind the object. In this case the overall gravitational force is applied in the direction of the object's motion, producing a negative dynamical friction (NDF). Black holes in dense gas drive powerful outflows and may experience the NDF, although extensive numerical work is probably needed to demonstrate or refute this conclusively. NDF may be important for stellar-mass black holes and neutron stars inside `common envelopes' in binary systems, for stellar mass black holes inside AGN discs, or for massive black holes growing through super-Eddington accretion in early Universe. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: 7 pages, 3 figures, submitted to MNRAS

Virial Ratio: Direct Evaluation from Molecular Cloud Data and the Challenges of Improving Accuracy

Authors: Ayushi Singh, Christopher D. Matzner, Peter H. Jumper

Abstract: The virial ratio between kinetic and gravitational terms provides key insight into the balance of forces that confine a molecular cloud, but the clumpy and filamentary structures of resolved clouds make it difficult to evaluate this ratio in a consistent way. For clouds with resolved maps of column density as well as a line tracer, we demonstrate that the gravitational energy can be estimated dire… ▽ More The virial ratio between kinetic and gravitational terms provides key insight into the balance of forces that confine a molecular cloud, but the clumpy and filamentary structures of resolved clouds make it difficult to evaluate this ratio in a consistent way. For clouds with resolved maps of column density as well as a line tracer, we demonstrate that the gravitational energy can be estimated directly from observations in a manner similar to the kinetic energy. This offers improved diagnostic power and consistency. Disentangling a cloud from foreground and background materials is a persistent challenge, for which we introduce a strategy based on Abel's transform. We provide proofs of principle using simulated clouds. △ Less

Submitted 13 June, 2019; v1 submitted 30 April, 2019; originally announced May 2019.

Comments: Published to ApJ, 9 pages, 6 figures, 1 table

Droplets I: Pressure-Dominated Sub-0.1 pc Coherent Structures in L1688 and B18

Authors: Hope How-Huan Chen, Jaime E. Pineda, Alyssa A. Goodman, Andreas Burkert, Stella S. R. Offner, Rachel K. Friesen, Philip C. Myers, Felipe Alves, Hector G. Arce, Paola Caselli, Ana Chacon-Tanarro, Michael Chun-Yuan Chen, James Di Francesco, Adam Ginsburg, Jared Keown, Helen Kirk, Peter G. Martin, Christopher Matzner, Anna Punanova, Elena Redaelli, Erik Rosolowsky, Samantha Scibelli, Young Min Seo, Yancy Shirley, Ayushi Singh

Abstract: We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transiti… ▽ More We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process. △ Less

Submitted 15 May, 2019; v1 submitted 26 September, 2018; originally announced September 2018.

Comments: Accepted by ApJ in April, 2019

Journal ref: 2019ApJ...877...93C

Capture and Escape: The Dependence of Radiation Forces on Clumping in Dusty Envelopes

Authors: Peter H. Jumper, Christopher D. Matzner

Abstract: Dust barriers effectively capture the photon momentum of a central light source, but low-density channels, along with re-emission at longer wavelengths, enhance its escape. We use Monte Carlo simulations to study the effects of inhomogeneity on radiation forces imparted to a dust envelope around a central star. We survey the strength and scale of an inhomogeneous perturbation field, as well as the… ▽ More Dust barriers effectively capture the photon momentum of a central light source, but low-density channels, along with re-emission at longer wavelengths, enhance its escape. We use Monte Carlo simulations to study the effects of inhomogeneity on radiation forces imparted to a dust envelope around a central star. We survey the strength and scale of an inhomogeneous perturbation field, as well as the optical depth of its spherical reference state. We run at moderate numerical resolution, relying on our previous resolution study for calibration of the associated error. We find that inhomogeneities matter most when their scale exceeds the characteristic mean free path. As expected, they tend to reduce the net radiation force and extend its range; however, there is significant variance among realizations. Within our models, force integrals correlate with the emergent spectral energy distribution, given a specified set of dust properties. A critical issue is the choice of integral measures of the radiation force: for strong deviations from spherical symmetry the relevant measures assess radial forces relative to the cloud centre of mass. Of these, we find the virial term due to radiation to be the least stochastic of several integral measures in the presence of inhomogeneities. △ Less

Submitted 6 September, 2018; v1 submitted 14 June, 2018; originally announced June 2018.

Comments: 11 pages, 15 figures, 1 table; Update to MNRAS version with an additional bibliographic correction; comments welcome

Aspherical Supernovae: Effects on Early Light Curves

Authors: Niloufar Afsariardchi, Christopher D. Matzner

Abstract: Early light from core-collapse supernovae, now detectable in high-cadence surveys, holds clues to a star and its environment just before it explodes. However, effects that alter the early light have not been fully explored. We highlight the possibility of non-radial flows at the time of shock breakout. These develop in sufficiently non-spherical explosions if the progenitor is not too diffuse. Whe… ▽ More Early light from core-collapse supernovae, now detectable in high-cadence surveys, holds clues to a star and its environment just before it explodes. However, effects that alter the early light have not been fully explored. We highlight the possibility of non-radial flows at the time of shock breakout. These develop in sufficiently non-spherical explosions if the progenitor is not too diffuse. When they do develop, non-radial flows limit ejecta speeds and cause ejecta-ejecta collisions. We explore these phenomena and their observational implications, using global, axisymmetric, non-relativistic FLASH simulations of simplified polytropic progenitors, which we scale to representative stars. We develop a method to track photon production within the ejecta, enabling us to estimate band-dependent light curves from adiabatic simulations. Immediate breakout emission becomes hidden as an oblique flow develops. Non-spherical effects lead the shock-heated ejecta to release a more constant luminosity at a higher, evolving color temperature at early times, effectively mixing breakout light with the early light curve. Collisions between non-radial ejecta thermalize a small fraction of the explosion energy; we address emission from these collisions in a subsequent paper. △ Less

Submitted 9 March, 2018; originally announced March 2018.

Comments: 23 pages, 20 figures, Accepted for publication in ApJ

Accretion and Feedback in Star Cluster Formation

Authors: Christopher D. Matzner

Abstract: Star cluster formation is unlikely to be a sudden event: instead, matter will flow to a cluster's formation site over an extended period, even as stars form and inject energy to the region. A cluster's gaseous precursor must persist under the competing influences of accretion and feedback for several crossing times, insofar as star formation is a slow process. The new-born stellar cluster should t… ▽ More Star cluster formation is unlikely to be a sudden event: instead, matter will flow to a cluster's formation site over an extended period, even as stars form and inject energy to the region. A cluster's gaseous precursor must persist under the competing influences of accretion and feedback for several crossing times, insofar as star formation is a slow process. The new-born stellar cluster should therefore preserve a memory of this competition. Using analytical approximations we assess the dynamical state of the gas, mapping regimes in which various types of feedback are weak or strong. Protostellar outflows, radiation pressure, and ionized gas pressure are accounted for. Comparison to observations shows that feedback is often incapable of expelling gas in the more massive, rapidly-accreting clusters, but feedback may nevertheless starve accretion by acting on larger scales. △ Less

Submitted 4 December, 2017; originally announced December 2017.

Comments: Proceedings of the Star Formation in Different Environments, ICISE, Quy Nhon, Vietnam, 2016, (eds. D. Johnstone, T. Hoang, F. Nakamura, Q. Nguyen-Luong, and J. T. Tranh Van)

Report number: ICISE/SFDE16/28

The Green Bank Ammonia Survey: Unveiling the Dynamics of the Barnard 59 star-forming Clump

Authors: E. Redaelli, F. O. Alves, P. Caselli, J. E. Pineda, R. K. Friesen, A. Chacón-Tanarro, C. D. Matzner, A. Ginsburg, E. Rosolowsky, J. Keown, S. S. R. Offner, J. Di Francesco, H. Kirk, P. C. Myers, A. Hacar, A. Cimatti, H. H. Chen, M. C. Chen, Y. M. Seo, K. I. Lee

Abstract: Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe N… ▽ More Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency. △ Less

Submitted 6 November, 2017; originally announced November 2017.

Comments: 18 pages, 18 figures

Radiation Forces on Dust Envelopes

Authors: Peter H. Jumper, Christopher D. Matzner

Abstract: We address in detail the radiation forces on spherical dust envelopes around luminous stars, and numerical solutions for these forces, as a first step toward more general dust geometries. Two physical quantities, a normalized force and a force-averaged radius, suffice to capture the overall effects of radiation forces; these combine to yield the radiation term in the virial theorem. In addition to… ▽ More We address in detail the radiation forces on spherical dust envelopes around luminous stars, and numerical solutions for these forces, as a first step toward more general dust geometries. Two physical quantities, a normalized force and a force-averaged radius, suffice to capture the overall effects of radiation forces; these combine to yield the radiation term in the virial theorem. In addition to the optically thin and thick regimes, the wavelength dependence of dust opacity allows for an intermediate case in which starlight is easily trapped but infrared radiation readily escapes. Scattering adds a non-negligible force in this intermediate regime. We address all three regimes analytically and provide approximate formulae for the force parameters, for arbitrary optical depth and inner dust temperature. Turning to numerical codes, we examine the convergence properties of the Monte Carlo code Hyperion run in Cartesian geometry. We calibrate both Hyperion and our analytical estimate using the DUSTY code, run in spherical geometry. We find that Monte Carlo codes tend to underestimate the radiation force when the mean free path of starlight is not well resolved, as this causes the inner dust temperature, and therefore the inner Rosseland opacity, to be too low. We briefly discuss implications for more complicated radiation transfer problems. △ Less

Submitted 29 June, 2018; v1 submitted 25 October, 2017; originally announced October 2017.

Comments: 10 pages, 5 figures, 6 tables, accepted to MNRAS, comments welcome; updates posting from the previous "submitted to MNRAS" version, incorporating revisions suggested by the reviewer

The Green Bank Ammonia Survey: Observations of Hierarchical Dense Gas Structures in Cepheus-L1251

Authors: Jared Keown, James Di Francesco, Helen Kirk, Rachel K. Friesen, Jaime E. Pineda, Erik Rosolowsky, Adam Ginsburg, Stella S. R. Offner, Paola Caselli, Felipe Alves, Ana Chacón-Tanarro, Anna Punanova, Elena Redaelli, Young Min Seo, Christopher D. Matzner, Michael Chun-Yuan Chen, Alyssa A. Goodman, How-Huan Chen, Yancy Shirley, Ayushi Singh, Hector G. Arce, Peter Martin, Philip C. Myers

Abstract: We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc… ▽ More We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ

The Green Bank Ammonia Survey: Dense Cores Under Pressure in Orion A

Authors: Helen Kirk, Rachel K. Friesen, Jaime E. Pineda, Erik Rosolowsky, Stella S. R. Offner, Christopher D. Matzner, Philip C. Myers, James Di Francesco, Paola Caselli, Felipe O. Alves, Ana Chacón-Tanarro, How-Huan Chen, Michael Chun-Yuan Chen, Jared Keown, Anna Punanova, Young Min Seo, Yancy Shirley, Adam Ginsburg, Christine Hall, Ayushi Singh, Héctor G. Arce, Alyssa A. Goodman, Peter Martin, Elena Redaelli

Abstract: We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity… ▽ More We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined. △ Less

Submitted 17 August, 2017; originally announced August 2017.

Comments: Accepted for publication in ApJ Aug 10, 2017

The Green Bank Ammonia Survey (GAS): First Results of NH3 mapping the Gould Belt

Authors: Rachel K. Friesen, Jaime E. Pineda, Erik Rosolowsky, Felipe Alves, Ana Chacón-Tanarro, Hope How-Huan Chen, Michael Chun-Yuan Chen, James Di Francesco, Jared Keown, Helen Kirk, Anna Punanova, Youngmin Seo, Yancy Shirley, Adam Ginsburg, Christine Hall, Stella S. R. Offner, Ayushi Singh, Héctor G. Arce, Paola Caselli, Alyssa A. Goodman, Peter G. Martin, Christopher Matzner, Philip C. Myers, Elena Redaelli

Abstract: We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt… ▽ More We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, $σ_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of $σ_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available. △ Less

Submitted 20 April, 2017; originally announced April 2017.

Comments: 33 pages, 27 figures, accepted to ApJS. Datasets are publicly available: https://dataverse.harvard.edu/dataverse/GAS_DR1

Shock Dynamics in Stellar Outbursts: I. Shock formation

Authors: Stephen Ro, Christopher D. Matzner

Abstract: Wave-driven outflows and non-disruptive explosions have been implicated in pre-supernova outbursts, supernova impostors, LBV eruptions, and some narrow-line and superluminous supernovae. To model these events, we investigate the dynamics of stars set in motion by strong acoustic pulses and wave trains, focusing here on nonlinear wave propagation, shock formation, and an early phase of the developm… ▽ More Wave-driven outflows and non-disruptive explosions have been implicated in pre-supernova outbursts, supernova impostors, LBV eruptions, and some narrow-line and superluminous supernovae. To model these events, we investigate the dynamics of stars set in motion by strong acoustic pulses and wave trains, focusing here on nonlinear wave propagation, shock formation, and an early phase of the development of a weak shock. We identify the shock formation radius, showing that a heuristic estimate based on crossing characteristics matches an exact expansion around the wave front and verifying both with numerical experiments. Our general analytical condition for shock formation applies to one-dimensional motions within any static environment, including both eruptions and implosions, and can easily be extended to non-stationary flows. We also consider the early phase of shock energy dissipation. We find that waves of super-Eddington acoustic luminosity always create shocks, rather than damping by radiative diffusion. Therefore, shock formation is integral to super-Eddington outbursts. △ Less

Submitted 28 December, 2016; originally announced December 2016.

Comments: 9 pages, 6 figures, submitted to Astrophysical Journal, comments welcome

On the Launching and Structure of Radiatively Driven Winds in Wolf-Rayet Stars

Authors: Stephen Ro, Christopher D. Matzner

Abstract: Hydrostatic models of Wolf-Rayet stars typically contain low-density outer envelopes that inflate the stellar radii by a factor of several and are capped by a denser shell of gas. Inflated envelopes and density inversions are hallmarks of envelopes that become super-Eddington as they cross the iron-group opacity peak, but these features disappear when mass loss is sufficiently rapid. We re-examine… ▽ More Hydrostatic models of Wolf-Rayet stars typically contain low-density outer envelopes that inflate the stellar radii by a factor of several and are capped by a denser shell of gas. Inflated envelopes and density inversions are hallmarks of envelopes that become super-Eddington as they cross the iron-group opacity peak, but these features disappear when mass loss is sufficiently rapid. We re-examine the structures of steady, spherically symmetric wind solutions that cross a sonic point at high optical depth, identifying the physical mechanism by which outflow affects the stellar structure, and provide an improved analytical estimate for the critical mass loss rate above which extended structures are erased. Weak-flow solutions below this limit resemble hydrostatic stars even in supersonic zones; however, we infer that these fail to successfully launch optically thick winds. Wolf-Rayet envelopes will therefore likely correspond to the strong, compact solutions. We also find that wind solutions with negligible gas pressure are stably stratified at and below the sonic point. This implies that convection is not the source of variability in Wolf-Rayet stars, as has been suggested; but, acoustic instabilities provide an alternative explanation. Our solutions are limited to high optical depths by our neglect of Doppler enhancements to the opacity, and do not account for acoustic instabilities at high Eddington factors; yet they provide useful insights into Wolf-Rayet stellar structures. △ Less

Submitted 23 February, 2016; originally announced February 2016.

Comments: 14 pages, 8 figures, Accepted by ApJ

Star cluster formation with stellar feedback and large-scale inflow

Authors: Christopher D. Matzner, Peter H. Jumper

Abstract: During star cluster formation, ongoing mass accretion is resisted by stellar feedback in the form of protostellar outflows from the low-mass stars and photo-ionization and radiation pressure feedback from the massive stars. We model the evolution of cluster-forming regions during a phase in which both accretion and feedback are present, and use these models to investigate how star cluster formatio… ▽ More During star cluster formation, ongoing mass accretion is resisted by stellar feedback in the form of protostellar outflows from the low-mass stars and photo-ionization and radiation pressure feedback from the massive stars. We model the evolution of cluster-forming regions during a phase in which both accretion and feedback are present, and use these models to investigate how star cluster formation might terminate. Protostellar outflows are the strongest form of feedback in low-mass regions, but these cannot stop cluster formation if matter continues to flow in. In more massive clusters, radiation pressure and photo-ionization rapidly clear the cluster-forming gas when its column density is too small. We assess the rates of dynamical mass ejection and of evaporation, while accounting for the important effect of dust opacity on photo-ionization. Our models are consistent with the census of protostellar outflows in NGC 1333 and Serpens South, and with the dust temperatures observed in regions of massive star formation. Comparing observations of massive cluster-forming regions against our model parameter space, and against our expectations for accretion-driven evolution, we infer that massive-star feedback is a likely cause of gas disruption in regions with velocity dispersions less than a few kilometers per second, but that more massive and more turbulent regions are too strongly bound for stellar feedback to be disruptive. △ Less

Submitted 10 November, 2015; originally announced November 2015.

Comments: 19 pages, 10 figures, accepted by ApJ

Molecular and Ionized Hydrogen in 30 Doradus. I. Imaging Observations

Authors: Sherry C. C. Yeh, Ernest R. Seaquist, Christopher D. Matzner, Eric W. Pellegrini

Abstract: We present the first fully calibrated H$_2$, 1-0 S(1) image of the entire 30 Doradus nebula. The observations were conducted using the NOAO Extremely Wide-Field Infrared Imager on the CTIO 4-meter Blanco Telescope. Together with a NEWFIRM Br$γ$ image of 30 Doradus, our data reveal the morphologies of the warm molecular gas and ionized gas in 30 Doradus. The brightest H$_2$-emitting area, which ext… ▽ More We present the first fully calibrated H$_2$, 1-0 S(1) image of the entire 30 Doradus nebula. The observations were conducted using the NOAO Extremely Wide-Field Infrared Imager on the CTIO 4-meter Blanco Telescope. Together with a NEWFIRM Br$γ$ image of 30 Doradus, our data reveal the morphologies of the warm molecular gas and ionized gas in 30 Doradus. The brightest H$_2$-emitting area, which extends from the northeast to the southwest of R136, is a photodissociation region viewed face-on, while many clumps and pillar features located at the outer shells of 30 Doradus are photodissociation regions viewed edge-on. Based on the morphologies of H$_2$, Br$γ$, $^{12}$CO, and 8$μ$m emission, the H$_2$ to Br$γ$ line ratio and Cloudy models, we find that the H$_2$ emission is formed inside the photodissociation regions of 30 Doradus, 2 - 3 pc to the ionization front of the HII region, in a relatively low-density environment $<$ 10$^4$ cm$^{-3}$. Comparisons with Br$γ$, 8$μ$m, and CO emission indicate that H$_2$ emission is due to fluorescence, and provide no evidence for shock excited emission of this line. △ Less

Submitted 20 April, 2015; originally announced April 2015.

Comments: 10 pages, 6 figures, 2 tables; accepted for publication in ApJ

Oblique Shock Breakout in Supernovae and Gamma-Ray Bursts: II. Numerical Solutions For Non-Relativistic Pattern Speeds

Authors: Pegah Salbi, Christopher D. Matzner, Stephen Ro, Yuri Levin

Abstract: Non-spherical explosions develop non-radial flows as the pattern of shock emergence progresses across the stellar surface. In supernovae these flows can limit ejecta speeds, stifle shock breakout emission, and cause collisions outside the star. Similar phenomena occur in stellar and planetary collisions, tidal disruption events, accretion-induced collapses, and propagating detonations. We present… ▽ More Non-spherical explosions develop non-radial flows as the pattern of shock emergence progresses across the stellar surface. In supernovae these flows can limit ejecta speeds, stifle shock breakout emission, and cause collisions outside the star. Similar phenomena occur in stellar and planetary collisions, tidal disruption events, accretion-induced collapses, and propagating detonations. We present two-dimensional, nested-grid Athena simulations of non-radial shock emergence in a frame comoving with the breakout pattern, focusing on the adiabatic, non-relativistic limit in a plane stratified envelope. We set boundary conditions using a known self-similar solution and explore the role of box size and resolution on the result. The shock front curves toward the stellar surface, and exhibits a kink from which weak discontinuities originate. Flow around the point of shock emergence is neither perfectly steady nor self-similar. Waves and vortices, which are not predominantly due to grid effects, emanate from this region. The post-shock flow is deflected along the stellar surface, and its pressure disturbs the stellar atmosphere upstream of the emerging shock. We use the numerical results and their analytical limits to predict the effects of radiation transfer and gravity, which are not included in our simulations. △ Less

Submitted 27 May, 2014; originally announced May 2014.

Comments: 15 pages, 12 figures, submitted to ApJ

Oblique Shock Breakout in Supernovae and Gamma-Ray Bursts: I. Dynamics and Observational Implications

Authors: Christopher D. Matzner, Yuri Levin, Stephen Ro

Abstract: In a non-spherical stellar explosion, non-radial motions become important near the stellar surface. For realistic deviations from spherical symmetry, non-radial flow dramatically alters the dynamics and emission of shock emergence on a significant fraction of the surface. The breakout flash is stifled, ejecta speeds are limited, and matter is cast sideways. Non-radial ejection allows for collision… ▽ More In a non-spherical stellar explosion, non-radial motions become important near the stellar surface. For realistic deviations from spherical symmetry, non-radial flow dramatically alters the dynamics and emission of shock emergence on a significant fraction of the surface. The breakout flash is stifled, ejecta speeds are limited, and matter is cast sideways. Non-radial ejection allows for collisions outside the star, which may engender a new type of transient. Strongly oblique breakouts are most easily produced in compact stellar progenitors, such as white dwarfs and stripped-envelope core collapse supernovae. We study the shock structure and post-shock acceleration using conservation laws, a similarity analysis, and an approximate theory for oblique shocks. The shock is likely to extend vertically from the stellar surface, then kink before joining a deep asymptotic solution. Outflow from the region crossed by an oblique shock is probably unsteady and may affect the surface ahead of the main shock. We comment on the implications for several notable explosions in which the non-spherical dynamics described in this paper are likely to play an important role. We also briefly consider relativistic and superluminal pattern speeds. △ Less

Submitted 28 October, 2013; originally announced October 2013.

Comments: 12 pages, 1 figure, accepted in Ap.J

Shock Emergence in Supernovae: Limiting Cases and Accurate Approximations

Authors: Stephen Ro, Christopher D. Matzner

Abstract: We examine the dynamics of accelerating normal shocks in stratified planar atmospheres, providing accurate fitting formulae for the scaling index relating shock velocity to the initial density and for the post-shock acceleration factor as functions of the polytropic and adiabatic indices which parameterize the problem. In the limit of a uniform initial atmosphere there are analytical formulae for… ▽ More We examine the dynamics of accelerating normal shocks in stratified planar atmospheres, providing accurate fitting formulae for the scaling index relating shock velocity to the initial density and for the post-shock acceleration factor as functions of the polytropic and adiabatic indices which parameterize the problem. In the limit of a uniform initial atmosphere there are analytical formulae for these quantities. In the opposite limit of a very steep density gradient the solutions match the outcome of shock acceleration in exponential atmospheres. △ Less

Submitted 25 June, 2013; originally announced June 2013.

Evolution of Accretion Disks in Tidal Disruption Events

Authors: Rong-Feng Shen, Christopher D. Matzner

Abstract: During a stellar tidal disruption event (TDE), an accretion disk forms as stellar debris returns to the disruption site and circularizes. Rather than being confined within the circularizing radius, the disk can spread to larger radii to conserve angular momentum. A spreading disk is a source of matter for re-accretion at rates which can exceed the later stellar fallback rate, although a disk wind… ▽ More During a stellar tidal disruption event (TDE), an accretion disk forms as stellar debris returns to the disruption site and circularizes. Rather than being confined within the circularizing radius, the disk can spread to larger radii to conserve angular momentum. A spreading disk is a source of matter for re-accretion at rates which can exceed the later stellar fallback rate, although a disk wind can suppress its contribution to the central black hole accretion rate. A spreading disk is detectible through a break in the central accretion rate history, or, at longer wavelengths, by its own emission. We model the evolution of TDE disk size and accretion rate, by accounting for the time-dependent fallback rate, for the influence of wind losses in the early, advective stage, and for the possibility of thermal instability for accretion rates intermediate between the advection-dominated and gas-pressure dominated states. The model provides a dynamic basis for modeling TDE light curves. All or part of a young TDE disk will precess as a solid body due to Lense-Thirring effect, and precession may manifest itself as quasi-periodic modulation of light curve. The precession period increases with time. Applying our results to the jetted TDE candidate Swift J1644+57, whose X-ray light curve shows numerous quasi-periodic dips, we argue that the data best fit a scenario in which a main-sequence star was fully disrupted by an intermediate mass black hole on an orbit significantly inclined from the black hole equator, with the apparent jet shutoff at t= 500 d corresponding to a disk transition from the advective state to the gas-pressure dominated state. △ Less

Submitted 28 January, 2014; v1 submitted 23 May, 2013; originally announced May 2013.

Comments: 23 pages, 11 figures, 2 appendices, accepted for ApJ

Journal ref: 2014, ApJ, 784, 87

Line Emission from Radiation-Pressurized HII Regions I: Internal Structure and Line Ratios

Authors: Sherry C. C. Yeh, Silvia Verdolini, Mark R. Krumholz, Christopher D. Matzner, Alexander G. G. M. Tielens

Abstract: The emission line ratios [OIII]5007/H-beta and [NII]6584/H-alpha have been adopted as an empirical way to distinguish between the fundamentally different mechanisms of ionization in emission-line galaxies. However, detailed interpretation of these diagnostics requires calculations of the internal structure of the emitting HII regions, and these calculations depend on the assumptions one makes abou… ▽ More The emission line ratios [OIII]5007/H-beta and [NII]6584/H-alpha have been adopted as an empirical way to distinguish between the fundamentally different mechanisms of ionization in emission-line galaxies. However, detailed interpretation of these diagnostics requires calculations of the internal structure of the emitting HII regions, and these calculations depend on the assumptions one makes about the relative importance of radiation pressure and stellar winds. In this paper we construct a grid of quasi-static HII region models to explore how choices about these parameters alter HII regions' emission line ratios. We find that, when radiation pressure is included in our models, HII regions reach a saturation point beyond which further increases in the luminosity of the driving stars does not produce any further increase in effective ionization parameter, and thus does not yield any further alteration in an HII region's line ratio. We also show that, if stellar winds are assumed to be strong, the maximum possible ionization parameter is quite low. As a result of this effect, it is inconsistent to simultaneously assume that HII regions are wind-blown bubbles and that they have high ionization parameters; some popular HII region models suffer from this inconsistency. Our work in this paper provides a foundation for a companion paper in which we embed the model grids we compute here within a population synthesis code that enables us to compute the integrated line emission from galactic populations of HII regions. △ Less

Submitted 19 March, 2013; originally announced March 2013.

Comments: 8 pages, 4 figures, accepted for publication in ApJ

Line Emission from Radiation-Pressurized HII Region II: Dynamics and Population Synthesis

Authors: Silvia Verdolini, Sherry C. C. Yeh, Mark R. Krumholz, Christopher D. Matzner, Alexander G. G. M. Tielens

Abstract: Optical and infrared emission lines from HII regions are an important diagnostic used to study galaxies, but interpretation of these lines requires significant modeling of both the internal structure and dynamical evolution of the emitting regions. Most of the models in common use today assume that HII region dynamics are dominated by the expansion of stellar wind bubbles, and have neglected the c… ▽ More Optical and infrared emission lines from HII regions are an important diagnostic used to study galaxies, but interpretation of these lines requires significant modeling of both the internal structure and dynamical evolution of the emitting regions. Most of the models in common use today assume that HII region dynamics are dominated by the expansion of stellar wind bubbles, and have neglected the contribution of radiation pressure to the dynamics, and in some cases also to the internal structure. However, recent observations of nearby galaxies suggest that neither assumption is justified, motivating us to revisit the question of how HII region line emission depends on the physics of winds and radiation pressure. In a companion paper we construct models of single HII regions including and excluding radiation pressure and winds, and in this paper we describe a population synthesis code that uses these models to simulate galactic collections of HII regions with varying physical parameters. We show that the choice of physical parameters has significant effects on galactic emission line ratios, and that in some cases the line ratios can exceed previously claimed theoretical limits. Our results suggest that the recently-reported offset in line ratio values between high-redshift star-forming galaxies and those in the local universe may be partially explained by the presence of large numbers of radiation pressured-dominated HII regions within them. △ Less

Submitted 19 March, 2013; originally announced March 2013.

Comments: 12 pages, 11 figures, accepted for publication in ApJ

OB Associations at the Upper End of the Milky Way Luminosity Function

Authors: Mubdi Rahman, Christopher D. Matzner, Dae-Sik Moon

Abstract: The Milky Way's most luminous, young and massive (M > 10^4 Msun) star clusters and OB associations have largely evaded detection despite knowledge of their surrounding H II regions. We search for these clusters and associations within the 40 star forming complexes from Rahman & Murray in the 13 most luminous WMAP free-free emission sources of the Galaxy. Selecting for objects with the dust-reddene… ▽ More The Milky Way's most luminous, young and massive (M > 10^4 Msun) star clusters and OB associations have largely evaded detection despite knowledge of their surrounding H II regions. We search for these clusters and associations within the 40 star forming complexes from Rahman & Murray in the 13 most luminous WMAP free-free emission sources of the Galaxy. Selecting for objects with the dust-reddened colors of OB stars, we identify new candidate associations using the 2MASS point source catalog. In 40 star forming complexes searched, 22 contain cluster/association candidates with sizes and masses in the range of 3' - 26' and 10^{2.3} - 10^{5} Msun. Of the 22 candidates, at least 7 have estimated masses > 10^4 Msun, doubling the number of such massive clusters known in the Galaxy. Applying our method to a statistically similar set of test locations, we estimate that 3.0 +/- 0.6 of our 22 candidate associations are unrelated to the star forming complexes. In addition, the apparent extinctions of our candidate associations correlate well with the predictions from a Galactic model. These facts, along with the clear detection of a known OB association and the previous spectral verification of one cluster found by this method, validate our method. Only one of the searched WMAP sources remains without a candidate. In 8 of the most luminous WMAP sources, the candidate associations can account for the observed free-free flux. With our new compilation, the Galactic census of young, massive stellar associations may now be about half complete. △ Less

Submitted 14 November, 2012; originally announced November 2012.

Comments: 25 pages, 14 figures. To appear in ApJ

Journal ref: ApJ, 2013, 766, 135

Ionization Parameter as a Diagnostic of Radiation and Wind Pressures in H II Regions and Starburst Galaxies

Authors: Sherry C. C. Yeh, Christopher D. Matzner

Abstract: The ionization parameter U is potentially useful for measuring radiation pressure feedback from massive star clusters, as it reflects the radiation-to-gas-pressure ratio and is readily derived from mid-infrared line ratios. We consider several effects which determine the apparent value of U in HII regions and galaxies. An upper limit is set by the compression of gas by radiation pressure. The pres… ▽ More The ionization parameter U is potentially useful for measuring radiation pressure feedback from massive star clusters, as it reflects the radiation-to-gas-pressure ratio and is readily derived from mid-infrared line ratios. We consider several effects which determine the apparent value of U in HII regions and galaxies. An upper limit is set by the compression of gas by radiation pressure. The pressure from stellar winds and the presence of neutral clumps both reduce U for a given radiation intensity. The most intensely irradiated regions are selectively dimmed by internal dust absorption of ionizing photons, inducing observational bias on galactic scales. We explore these effects analytically and numerically, and use them to interpret previous observational results. We find that radiation confinement sets the upper limit log_10 U = -1 seen in individual regions. Unresolved starbursts display a maximum value of ~ -2.3. While lower, this is also consistent with a large portion of their HII regions being radiation dominated, given the different technique used to interpret unresolved regions, and given the bias caused by dust absorption. We infer that many individual, strongly illuminated regions cannot be dominated by stellar winds, and that even when averaged on galactic scales, shocked wind pressures cannot be large compared to radiation pressure. Therefore, most HII regions cannot be adiabatic wind bubbles. Our models imply a metallicity dependence in the physical structure and dust attenuation of radiation-dominated regions, both of which should vary strongly across a critical metallicity of about one-twentieth solar. △ Less

Submitted 27 June, 2012; originally announced June 2012.

Comments: 19 pages, 14 figures, accepted for publication in ApJ

Journal ref: 2012ApJ...757..108Y

Spectroscopic Confirmation of the Dragonfish Association: The Galaxy's Most Luminosity OB Association

Authors: Mubdi Rahman, Dae-Sik Moon, Christopher D. Matzner

Abstract: Young OB associations with masses greater than 10$^4$ M$_{\sun}$ have been inferred to exist in the Galaxy but have largely evaded detection. Recently, a candidate OB association has been identified within the most luminous star forming complex in the Galaxy, the Dragonfish Nebula. We identify 18 young, massive stars with near-infrared spectroscopy from a sample of 50 members within the candidate… ▽ More Young OB associations with masses greater than 10$^4$ M$_{\sun}$ have been inferred to exist in the Galaxy but have largely evaded detection. Recently, a candidate OB association has been identified within the most luminous star forming complex in the Galaxy, the Dragonfish Nebula. We identify 18 young, massive stars with near-infrared spectroscopy from a sample of 50 members within the candidate OB association, including 15 O-type, and three Luminous Blue Variables or Wolf-Rayet stars. This number matches the expected yield of massive stars from the candidate association, confirming its existence and ability to power the parent star forming complex. These results demonstrate the existence of a 10$^5$ M$_{\sun}$ OB association, more powerful than any previously known in the Galaxy, comparable in mass only to Westerlund 1. Further, the results also validate the color selection method used to identify the association, adding credence to others discovered in the same way. △ Less

Submitted 14 November, 2011; originally announced November 2011.

Comments: 13 pages, 3 figures, 1 table. Accepted to ApJL

Journal ref: 2011, ApJL, 743, L28

Coasting external shock in wind medium: an origin for the X-ray plateau decay component in Swift GRB afterglows

Authors: Rongfeng Shen, Christopher D. Matzner

Abstract: The plateaus observed in about one half of the early X-ray afterglows are the most puzzling feature in gamma-ray bursts (GRBs) detected by Swift. By analyzing the temporal and spectral indices of a large X-ray plateau sample, we find that 55% can be explained by external, forward shock synchrotron emission produced by a relativistic ejecta coasting in a ρ~ r^{-2}, wind-like medium; no energy injec… ▽ More The plateaus observed in about one half of the early X-ray afterglows are the most puzzling feature in gamma-ray bursts (GRBs) detected by Swift. By analyzing the temporal and spectral indices of a large X-ray plateau sample, we find that 55% can be explained by external, forward shock synchrotron emission produced by a relativistic ejecta coasting in a ρ~ r^{-2}, wind-like medium; no energy injection into the shock is needed. After the ejecta collects enough medium and transitions to the adiabatic, decelerating blastwave phase, it produces the post-plateau decay. For those bursts consistent with this model, we find an upper limit for the initial Lorentz factor of the ejecta, Γ_0 \leq 46 (ε_e/0.1)^{-0.24} (ε_B/0.01)^{0.17}; the isotropic equivalent total ejecta energy is E_{iso} ~ 10^{53} (ε_e/0.1)^{-1.3} (ε_B/0.01)^{-0.09} (t_b/10^4 s) erg, where ε_e and ε_B are the fractions of the total energy at the shock downstream that are carried by electrons and the magnetic field, respectively, and t_b is the end of the plateau. Our finding supports Wolf-Rayet stars as the progenitor stars of some GRBs. It raises intriguing questions about the origin of an intermediate-Γ_0 ejecta, which we speculate is connected to the GRB jet emergence from its host star. For the remaining 45% of the sample, the post-plateau decline is too rapid to be explained in the coasting-in-wind model, and energy injection appears to be required. △ Less

Submitted 29 September, 2011; v1 submitted 15 September, 2011; originally announced September 2011.

Comments: 11 pages, 5 figures, to appear in ApJ, proof-corrected version, added more references

Journal ref: 2012, ApJ, 744, 36

The Global Evolution of Giant Molecular Clouds II: The Role of Accretion

Authors: Nathan J. Goldbaum, Mark R. Krumholz, Christopher D. Matzner, Christopher F. McKee

Abstract: We present virial models for the global evolution of giant molecular clouds. Focusing on the presence of an accretion flow, and accounting for the amount of mass, momentum, and energy supplied by accretion and star formation feedback, we are able to follow the growth, evolution, and dispersal of individual giant molecular clouds. Our model clouds reproduce the scaling relations observed in both ga… ▽ More We present virial models for the global evolution of giant molecular clouds. Focusing on the presence of an accretion flow, and accounting for the amount of mass, momentum, and energy supplied by accretion and star formation feedback, we are able to follow the growth, evolution, and dispersal of individual giant molecular clouds. Our model clouds reproduce the scaling relations observed in both galactic and extragalactic clouds. We find that accretion and star formation contribute contribute roughly equal amounts of turbulent kinetic energy over the lifetime of the cloud. Clouds attain virial equilibrium and grow in such a way as to maintain roughly constant surface densities, with typical surface densities of order 50 - 200 Msun pc^-2, in good agreement with observations of giant molecular clouds in the Milky Way and nearby external galaxies. We find that as clouds grow, their velocity dispersion and radius must also increase, implying that the linewidth-size relation constitutes an age sequence. Lastly, we compare our models to observations of giant molecular clouds and associated young star clusters in the LMC and find good agreement between our model clouds and the observed relationship between H ii regions, young star clusters, and giant molecular clouds. △ Less

Submitted 15 June, 2011; v1 submitted 30 May, 2011; originally announced May 2011.

Comments: 23 Pages, 9 Figures. Accepted to ApJ