Source author record

Steven W. Stahler

Steven W. Stahler 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

12works
3topics
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

12 published item(s)

preprint2019arXiv

The Nature of Class I Sources: Periodic Variables in Orion

We present a quantitative, empirically based argument that at least some Class I sources are low-mass, pre-main-sequence stars surrounded by spatially extended envelopes of dusty gas. The source luminosity arises principally from stellar gravitational contraction, as in optically visible pre-main-sequence stars that lack such envelopes. We base our argument on the fact that some Class I sources in Orion and other star-forming regions have been observed by Spitzer to be periodic variables in the mid-infrared, and with periods consistent with T Tauri rotation rates. Using a radiative transfer code, we construct a variety of dust envelopes surrounding rotating, spotted stars, to see if an envelope that produces a Class I SED at least broadly matches the observed modulations in luminosity. Acceptable envelopes can either be spherical or flattened, and may or may not have polar cavities. The key requirement is that they have a modest equatorial optical depth at the Spitzer waveband of 3.6 $μ$m, typically ${τ_{3.6}}$ ${\approx}$ 0.6. The total envelope mass, based on this limited study, is at most about 0.1 $\text{M}_{\odot}$, less than a typical stellar mass. Future studies should focus on the dynamics of the envelope, to determine whether material is actually falling onto the circumstellar disk.

preprint2014arXiv

Two paths of cluster evolution: global expansion versus core collapse

All gravitationally bound clusters expand, due to both gas loss from their most massive members and binary heating. All are eventually disrupted tidally, either by passing molecular clouds or the gravitational potential of their host galaxies. However, their interior evolution can follow two very different paths. Only clusters of sufficiently large initial population and size undergo the combined interior contraction and exterior expansion that leads eventually to core collapse. In all other systems, core collapse is frustrated by binary heating. These clusters globally expand for their entire lives, up to the point of tidal disruption. Using a suite of direct N-body calculations, we trace the "collapse line" in r_v-N space that separates these two paths. Here, r_v and N are the cluster's initial virial radius and population, respectively. For realistic starting radii, the dividing N-value is from 10^4 to over 10^5. We also show that there exists a minimum population, N_min, for core collapse. Clusters with N < N_min tidally disrupt before core collapse occurs. At the Sun's Galactocentric radius, R_G = 8.5 kpc, we find N_min >~ 300. The minimum population scales with Galactocentric radius as R_G^{-9/8}. The position of an observed cluster relative to the collapse line can be used to predict its future evolution. Using a small sample of open clusters, we find that most lie below the collapse line, and thus will never undergo core collapse. Most globular clusters, on the other hand, lie well above the line. In such a case, the cluster may or may not go through core collapse, depending on its initial size. We show how an accurate age determination can help settle this issue.

preprint2013arXiv

Dynamical Friction in a Gas: The Supersonic Case

Any gravitating mass traversing a relatively sparse gas experiences a retarding force created by its disturbance of the surrounding medium. In a previous contribution (Lee & Stahler 2011), we determined this dynamical friction force when the object's velocity was subsonic. We now extend our analysis to the supersonic regime. As before, we consider small perturbations created in the gas far from the gravitating object, and thereby obtain the net influx of linear momentum over a large, bounding surface. Various terms in the perturbation series formally diverge, necessitating an approximate treatment of the flow streamlines. Nevertheless, we are able to derive exactly the force itself. As in the subsonic case, we find that F=Mdot*V, where Mdot is the rate of mass accretion onto the object and V its instantaneous velocity with respect to distant background gas. Our force law holds even when the object is porous (e.g., a galaxy) or is actually expelling mass in a wind. Quantitatively, the force in the supersonic regime is less than that derived analytically by previous researchers, and is also less than was found in numerical simulations through the mid 1990s. We urge simulators to revisit the problem using modern numerical techniques. Assuming our result to be correct, it is applicable to many fields of astrophysics, ranging from exoplanet studies to galactic dynamics.

preprint2013arXiv

Externally Fed Star Formation: A Numerical Study

We investigate, through a series of numerical calculations, the evolution of dense cores that are accreting external gas up to and beyond the point of star formation. Our model clouds are spherical, unmagnetized configurations with fixed outer boundaries, across which gas enters subsonically. When we start with any near-equilibrium state, we find that the cloud's internal velocity also remains subsonic for an extended period, in agreement with observations. However, the velocity becomes supersonic shortly before the star forms. Consequently, the accretion rate building up the protostar is much greater than the benchmark value c_s^3/G, where c_s is the sound speed in the dense core. This accretion spike would generate a higher luminosity than those seen in even the most embedded young stars. Moreover, we find that the region of supersonic infall surrounding the protostar races out to engulf much of the cloud, again in violation of the observations, which show infall to be spatially confined. Similar problematic results have been obtained by all other hydrodynamic simulations to date, regardless of the specific infall geometry or boundary conditions adopted. Low-mass star formation is evidently a quasi-static process, in which cloud gas moves inward subsonically until the birth of the star itself. We speculate that magnetic tension in the cloud's deep interior helps restrain the infall prior to this event.

preprint2013arXiv

The Dynamics of Ultracompact HII Regions

Many ultracompact HII regions exhibit a cometary morphology in radio continuum emission. In such regions, a young massive star is probably ablating, through its ultraviolet radiation, the molecular cloud clump that spawned it. On one side of the star, the radiation drives an ionization front that stalls in dense molecular gas. On the other side, ionized gas streams outward into the more rarefied environment. This wind is underpressured with respect to the neutral gas. The difference in pressure draws in more cloud material, feeding the wind until the densest molecular gas is dissipated. Recent, time-dependent simulations of massive stars turning on within molecular gas show the system evolving in a direction similar to that just described. Here, we explore a semi-analytic model in which the wind is axisymmetric and has already achieved a steady state. Adoption of this simplified picture allows us to study the dependence of both the wind and its bounding ionization front on the stellar luminosity, the peak molecular density, and the displacement of the star from the center of the clump. For typical parameter values, the wind accelerates transonically to a speed of about 15 km / s, and transports mass outward at a rate of 10^-4 solar masses per year. Stellar radiation pressure acts to steepen the density gradient of the wind.

preprint2012arXiv

Externally Fed Accretion onto Protostars

The asymmetric molecular emission lines from dense cores reveal slow, inward motion in the clouds' outer regions. This motion is present both before and after the formation of a central star. Motivated by these observations, we revisit the classic problem of steady, spherical accretion of gas onto a gravitating point mass, but now include self-gravity of the gas and impose a finite, subsonic velocity as the outer boundary condition. We find that the accretion rate onto the protostar is lower than values obtained for isolated, collapsing clouds, by a factor that is the Mach number of the outer flow. Moreover, the region of infall surrounding the protostar spreads out more slowly, at a speed close to the subsonic, incoming velocity. Our calculation, while highly idealized, provides insight into two longstanding problems -- the surprisingly low accretion luminosities of even the most deeply embedded stellar sources, and the failure so far to detect spatially extended, supersonic infall within their parent dense cores. Indeed, the observed subsonic contraction in the outer regions of dense cores following star formation appears to rule out a purely hydrodynamic origin for these clouds.

preprint2011arXiv

Dynamical Friction in a Gas: The Subsonic Case

We study the force of dynamical friction acting on a gravitating point mass that travels through an extended, isothermal gas. This force is well established in the hypersonic limit, but remains less understood in the subsonic regime. Using perturbation theory, we analyze the changes in gas velocity and density far from the mass. We show analytically that the steady-state friction force is Mdot*V, where Mdot is the mass accretion rate onto an object moving at speed V. It follows that the speed of an object experiencing no other forces declines as the inverse square of its mass. Using a modified version of the classic Bondi-Hoyle interpolation formula for Mdot as a function of V, we derive an analytic expression for the friction force. This expression also holds when mass accretion is thwarted, e.g. by a wind, as long as the wind-cloud interaction is sufficiently confined spatially. Our result should find application in a number of astrophysical settings, such as the motion of galaxies through intracluster gas.

preprint2010arXiv

Molecular Line Profiles from Contracting Dense Cores

We recently proposed that molecular cloud dense cores undergo a prolonged period of quasi-static contraction prior to true collapse. This theory could explain the observation that many starless cores exhibit, through their spectral line profiles, signs of inward motion. We now use our model, together with a publicly available radiative transfer code, to determine the emission from three commonly used species - N2H+, CS, and HCN. A representative dense core of 3 Msun that has been contracting for 1 Myr has line profiles that qualitatively match the observed ones. In particular, optically thick lines have about the right degree of blue-red asymmetry, the empirical hallmark of contraction. The J=2-1 rotational transition of CS only attains the correct type of profile if the species is centrally depleted, as has been suggested by previous studies. These results support the idea that a slow, but accelerating, contraction leads to protostellar collapse. In the future, the kind of analysis presented here can be used to assign ages to individual starless cores.

preprint2010arXiv

Star Clusters Under Stress: Why Small Systems Cannot Dynamically Relax

Utilizing a series of N-body simulations, we argue that gravitationally bound stellar clusters of modest population evolve very differently from the picture presented by classical dynamical relaxation theory. The system's most massive stars rapidly sink towards the center and form binary systems. These binaries efficiently heat the cluster, reversing any incipient core contraction and driving a subsequent phase of global expansion. Most previous theoretical studies demonstrating deep and persistent dynamical relaxation have either conflated the process with mass segregation, ignored three-body interactions, or else adopted the artificial assumption that all cluster members are single stars of identical mass. In such a uniform-mass cluster, binary formation is greatly delayed, as we confirm here both numerically and analytically. The relative duration of core contraction and global expansion is effected by stellar evolution, which causes the most massive stars to die out before they form binaries. In clusters of higher N, the epoch of dynamical relaxation lasts for progressively longer periods. By extrapolating our results to much larger populations, we can understand, at least qualitatively, why some globular clusters reach the point of true core collapse.

preprint2010arXiv

The Dynamical Evolution of the Pleiades

We present the results of a numerical simulation of the history and future development of the Pleiades. This study builds on our previous one that established statistically the present-day structure of this system. Our simulation begins just after molecular cloud gas has been expelled by the embedded stars. We then follow, using an N body code, the stellar dynamical evolution of the cluster to the present and beyond. Our initial state is that which evolves, over the 125 Myr age of the cluster, to a configuration most closely matching the current one. We find that the original cluster, newly stripped of gas, already had a virial radius of 4 pc. This configuration was larger than most observed, embedded clusters. Over time, the cluster expanded further and the central surface density fell by about a factor of two. We attribute both effects to the liberation of energy from tightening binaries of short period. Indeed, the original binary fraction was close to unity. The ancient Pleiades also had significant mass segregation, which persists in the cluster today. In the future, the central density of the Pleiades will continue to fall. For the first few hundred Myr, the cluster as a whole will expand because of dynamical heating by binaries. The expansion process is aided by mass loss through stellar evolution, which weakens the system's gravitational binding. At later times, the Galactic tidal field begins to heavily deplete the cluster mass. It is believed that most open clusters are eventually destroyed by close passage of a giant molecular cloud. Barring that eventuality, the density falloff will continue for as long as 1 Gyr, by which time most of the cluster mass will have been tidally stripped away by the Galactic field.

preprint2009arXiv

The Orbital Decay of Embedded Binary Stars

Young binaries within dense molecular clouds are subject to dynamical friction from ambient gas. Consequently, their orbits decay, with both the separation and period decreasing in time. A simple analytic expression is derived for this braking torque. The derivation utilizes the fact that each binary acts as a quadrupolar source of acoustic waves. The acoustic disturbance has the morphology of a two-armed spiral and carries off angular momentum. From the expression for the braking torque, the binary orbital evolution is also determined analytically. This type of merger may help explain the origin of high-mass stars. If infrared dark clouds, with peak densities up to 10^7 cm^{-3}, contain low-mass binaries, those with separations less than 100 AU merge within about 10^5 yr. During the last few thousand years of the process, the rate of mechanical energy deposition in the gas exceeds the stars' radiative luminosity. Successive mergers may lead to the massive star formation believed to occur in these clouds.

preprint2003arXiv

A Near-Infrared Multiplicity Survey of Class I/Flat-Spectrum Systems in Six Nearby Molecular Clouds

We present new near-IR observations of 76 Class I/flat-spectrum objects in the nearby (d < 320 pc) Perseus, Taurus, Chamaeleon I and II, rho Ophiuchi, and Serpens dark clouds. These observations are part of a larger systematic infrared multiplicity survey of self-embedded objects in the nearest dark clouds. When combined with the results of our previously published near-infrared multiplicity survey, we find a restricted companion star fraction of 14/79 (18% +/- 4%) of the sources surveyed to be binary or higher order multiple systems over a separation range of ~300 - 2000 AU with a magnitude difference $Δ$ K <= 4, and with no correction for background contamination or completeness. This is consistent with the fraction of binary/multiple systems found among older pre-main-sequence T Tauri stars in each of the Taurus, rho Ophiuchi, and Chamaeleon star-forming regions over a similar separation range, as well as the combined companion star fraction for these regions. However, the companion star fraction for solar-type, and lower mass M dwarf, main-sequence stars in the solar neighborhood in this separation range (11% +/- 3%) is approximately one-half that of our sample. Together with multiplicity statistics derived for previously published samples of Class 0 and Class I sources, our study suggests that a significant number of binary/multiple objects may remain to be discovered at smaller separations among our Class I/flat-spectrum sample and/or most of the evolution of binary/multiple systems occurs during the Class 0 phase of early stellar evolution.