-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Numerical Parameter Survey of Nonradiative Black Hole Accretion -- Flow Structure and Variability of the Rotation Measure
Authors:
Bijia Pang,
Ue-Li Pen,
Christopher D. Matzner,
Stephen R. Green,
Matthias Liebendörfer
Abstract:
We conduct a survey of numerical simulations to probe the structure and appearance of non-radiative black hole accretion flows like the supermassive black hole at the Galactic centre. We find a generic set of solutions, and make specific predictions for currently feasible rotation measure (RM) observations, which are accessible to current instruments including the EVLA, GMRT and ALMA. The slow tim…
▽ More
We conduct a survey of numerical simulations to probe the structure and appearance of non-radiative black hole accretion flows like the supermassive black hole at the Galactic centre. We find a generic set of solutions, and make specific predictions for currently feasible rotation measure (RM) observations, which are accessible to current instruments including the EVLA, GMRT and ALMA. The slow time variability of the RM is a key quantitative signature of this accretion flow. The time variability of RM can be used to quantitatively measure the nature of the accretion flow, and to differentiate models. Sensitive measurements of RM can be achieved using RM synthesis or using pulsars.
Our energy conserving ideal magneto-hydrodynamical simulations, which achieve high dynamical range by means of a deformed-mesh algorithm, stretch from several Bondi radii to about one thousandth of that radius, and continue for tens of Bondi times. Magnetized flows which lack outward convection possess density slopes around -1, almost independent of physical parameters, and are more consistent with observational constraints than are strongly convective flows We observe no tendency for the flows to become rotationally supported in their centres, or to develop steady outflow.
We support these conclusions with formulae which encapsulate our findings in terms of physical and numerical parameters. We discuss the relation of these solutions to other approaches. The main potential uncertainties are the validity of ideal MHD and the absence of a fully relativistic inner boundary condition. The RM variability predictions are testable with current and future telescopes.
△ Less
Submitted 24 November, 2010; v1 submitted 24 November, 2010;
originally announced November 2010.
-
The Formation of Low-Mass Binary Star Systems Via Turbulent Fragmentation
Authors:
S. S. R. Offner,
K. M. Kratter,
C. D. Matzner,
M. R. Krumholz,
R. I. Klein
Abstract:
We characterize the infall rate onto protostellar systems forming in self-gravitating radiation-hydrodynamic simulations. Using two dimensionless parameters to determine disks' susceptability to gravitational fragmentation, we infer limits on protostellar system multiplicity and the mechanism of binary formation. We show that these parameters give robust predictions even in the case of marginally…
▽ More
We characterize the infall rate onto protostellar systems forming in self-gravitating radiation-hydrodynamic simulations. Using two dimensionless parameters to determine disks' susceptability to gravitational fragmentation, we infer limits on protostellar system multiplicity and the mechanism of binary formation. We show that these parameters give robust predictions even in the case of marginally resolved protostellar disks. We find that protostellar systems with radiation feedback predominately form binaries via turbulent fragmentation, not disk instability, and we predict turbulent fragmentation is the dominant channel for binary formation for low-mass stars. We clearly demonstrate that systems forming in simulations including radiative feedback have fundamentally different parameters than those in purely hydrodynamic simulations.
△ Less
Submitted 18 October, 2010;
originally announced October 2010.
-
The JCMT Legacy Survey of the Gould Belt: a first look at Serpens with HARP
Authors:
S. F. Graves,
J. S. Richer,
J. V. Buckle,
A. Duarte-Cabral,
G. A. Fuller,
M. R. Hogerheijde,
J. E. Owen,
C. Brunt,
H. M. Butner,
B. Cavanagh,
A. Chrysostomou,
E. I. Curtis,
C. J. Davis,
M. Etxaluze,
J. Di Francesco,
P. Friberg,
R. K. Friesen,
J. S. Greaves,
J. Hatchell,
D. Johnstone,
B. Matthews,
H. Matthews,
C. D. Matzner,
D. Nutter,
J. M. C. Rawlings
, et al. (10 additional authors not shown)
Abstract:
The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated,…
▽ More
The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated, with special consideration for high velocity gas. We find the cloud to have a mass of 203 solar masses, and to be gravitationally bound, and that the kinetic energy of the outflowing gas is approximately seventy percent of the turbulent kinetic energy of the cloud. We identify compact outflows towards some of the submillimetre Class 0/I sources in the region
△ Less
Submitted 4 June, 2010;
originally announced June 2010.
-
Stellar Feedback in Molecular Clouds and its Influence on the Mass Function of Young Star Clusters
Authors:
S. Michael Fall,
Mark R. Krumholz,
Christopher D. Matzner
Abstract:
We investigate how the removal of interstellar material by stellar feedback limits the efficiency of star formation in molecular clouds and how this determines the shape of the mass function of young star clusters. In particular, we derive relations between the power-law exponents of the mass functions of the clouds and clusters in the limiting regimes in which the feedback is energy-driven and…
▽ More
We investigate how the removal of interstellar material by stellar feedback limits the efficiency of star formation in molecular clouds and how this determines the shape of the mass function of young star clusters. In particular, we derive relations between the power-law exponents of the mass functions of the clouds and clusters in the limiting regimes in which the feedback is energy-driven and momentum-driven, corresponding to minimum and maximum radiative losses and likely to bracket all realistic cases. We find good agreement between the predicted and observed exponents, especially for momentum-driven feedback, provided the protoclusters have roughly constant mean surface density, as indicated by observations of the star-forming clumps within molecular clouds. We also consider a variety of specific feedback mechanisms, concluding that HII regions inflated by radiation pressure predominate in massive protoclusters, a momentum-limited process when photons can escape after only a few interactions with dust grains. We then present a first estimate of the star formation efficiency in this case, finding that it depends on the masses and sizes of the protoclusters only through their mean surface density, thus ensuring consistency between the observed power-law exponents of the mass functions of the clouds and clusters. The numerical value of this efficiency is also consistent with observations.
△ Less
Submitted 4 February, 2010; v1 submitted 12 October, 2009;
originally announced October 2009.
-
On the Role of Disks in the Formation of Stellar Systems: A Numerical Parameter Study of Rapid Accretion
Authors:
Kaitlin M. Kratter,
Christopher D. Matzner,
Mark R. Krumholz,
Richard I. Klein
Abstract:
We study rapidly accreting, gravitationally unstable disks with a series of global, three dimensional, numerical experiments using the code ORION. In this paper we conduct a numerical parameter study focused on protostellar disks, and show that one can predict disk behavior and the multiplicity of the accreting star system as a function of two dimensionless parameters which compare the disk's ac…
▽ More
We study rapidly accreting, gravitationally unstable disks with a series of global, three dimensional, numerical experiments using the code ORION. In this paper we conduct a numerical parameter study focused on protostellar disks, and show that one can predict disk behavior and the multiplicity of the accreting star system as a function of two dimensionless parameters which compare the disk's accretion rate to its sound speed and orbital period. Although gravitational instabilities become strong, we find that fragmentation into binary or multiple systems occurs only when material falls in several times more rapidly than the canonical isothermal limit. The disk-to-star accretion rate is proportional to the infall rate, and governed by gravitational torques generated by low-m spiral modes. We also confirm the existence of a maximum stable disk mass: disks that exceed ~50% of the total system mass are subject to fragmentation and the subsequent formation of binary companions.
△ Less
Submitted 20 July, 2009;
originally announced July 2009.