Source author record

Christopher D. Matzner

Christopher D. Matzner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2022arXiv

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

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.

preprint2022arXiv

How to search for Mirror Stars with Gaia

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.

preprint2022arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2019arXiv

Astronomy in a Low-Carbon Future

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.

preprint2019arXiv

Negative Dynamical Friction on compact objects moving through dense gas

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.

preprint2016arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

Star cluster formation with stellar feedback and large-scale inflow

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.

preprint2014arXiv

Evolution of Accretion Disks in Tidal Disruption Events

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.

preprint2014arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

Shock Emergence in Supernovae: Limiting Cases and Accurate Approximations

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2011arXiv

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

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.

preprint2011arXiv

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

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.

preprint2011arXiv

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

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.

preprint2010arXiv

Numerical Parameter Survey of Nonradiative Black Hole Accretion -- Flow Structure and Variability of the Rotation Measure

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.

preprint2010arXiv

Stellar Feedback in Molecular Clouds and its Influence on the Mass Function of Young Star Clusters

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.

preprint2009arXiv

On the Role of Disks in the Formation of Stellar Systems: A Numerical Parameter Study of Rapid Accretion

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.

preprint2007arXiv

Protostellar outflow-driven turbulence

Protostellar outflows crisscross the regions of star cluster formation, stirring turbulence and altering the evolution of the forming cluster. We model the stirring of turbulent motions by protostellar outflows, building on an observation that the scaling law of supersonic turbulence implies a momentum cascade analogous to the energy cascade in Kolmogorov turbulence. We then generalize this model to account for a diversity of outflow strengths, and for outflow collimation, both of which enhance turbulence. For a single value of its coupling coefficient the model is consistent with turbulence simulations by Li & Nakamura and, plausibly, with observations of the NGC 1333 cluster-forming region. Outflow-driven turbulence is strong enough to stall collapse in cluster-forming regions for several crossing times, relieving the mismatch between star formation and turbulent decay rates. The predicted line-width-size scaling implies radial density indices between -1 and -2 for regions supported by outflow-driven turbulence, with a tendency for steeper profiles in regions that are more massive or have higher column densities.