Processing of chondritic and planetary material in spiral density waves in the nebula

Abstract— A widely held view of nebular evolution is that during the ～0.5 Ma while interstellar material was collapsing onto the disk, the latter grew in mass to the point of gravitational instability. It responded to this by losing axial symmetry, growing spiral arms that had the capacity to tidally redistribute disk mass (inward) and angular momentum (outward) and prevent further increase in the disk/protosun mass ratio. The spiral arms (density waves) rotated differently than the substance of the nebula, and in some parts of the disk, nebular material may have encountered the arms at supersonic velocities. The disk gas, and solid particles entrained in it, would have been heated to some degree when they passed through shock fronts at the leading edges of the spiral arms. The present paper proposes this was the energetic nebular setting or environment that has long been sought, in which the material now in the planets and chondritic meteorites was thermally processed.     

Low T/|W| dynamical instability in differentially rotating stars: diagnosis with canonical angular momentum

We study the nature of non-axisymmetric dynamical instabilities in differentially rotating stars with both linear eigenmode analysis and hydrodynamic simulations in Newtonian gravity. We especially investigate the following three types of instability; the one-armed spiral instability, the low   T /| W |  bar instability, and the high   T /| W |  bar instability, where T is the rotational kinetic energy and W is the gravitational potential energy. The nature of the dynamical instabilities is clarified by using a canonical angular momentum as a diagnostic. We find that the one-armed spiral and the low   T /| W |  bar instabilities occur around the corotation radius, and they grow through the inflow of canonical angular momentum around the corotation radius. The result is a clear contrast to that of a classical dynamical bar instability in high   T /| W |  . We also discuss the feature of gravitational waves generated from these three types of instability.     

The specific features of the evolution of the viscous protoplanetary circumsolar disk

The problem of angular-momentum and mass transport in the disk is discussed and the disk viscosity is estimated. The evolution of the gas-dust protoplanetary disk at the stage of its formation inside the protostellar (protosolar) accretion envelope is considered. The conditions for the radial growth of the disk are estimated. For the subsequent period, when the central star (young Sun) is in the T Tauri phase, the temporal variations of the radius, mass, and the surface density of the disk, as well as the total mass flux from the disk onto the star (Sun), i.e., the mass accretion rate, are evaluated. The constraints on the initial value of the angular momentum of the protoplanetary circumsolar disk (that is, on the angular momentum of the protosolar cloud) are discussed with due regard for cosmochemical data.Translated from Astronomicheskii Vestnik, Vol. 38, No. 6, 2004, pp. 559–576.Original Russian Text Copyright © 2004 by Makalkin.     

The formation of galaxy discs

Galaxy disc formation must incorporate the multiphase nature of the interstellar medium. The resulting two-phase structure is generated and maintained by gravitational instability and supernova energy input, which yield a source of turbulent viscosity that is able to compete effectively in the protodisc phase with early angular momentum loss of the baryonic component via dynamical friction in the dark halo. Provided that star formation occurs on the viscous drag time-scale, this mechanism provides a means of accounting for disc sizes and radial profiles. The star formation feedback is self-regulated by turbulent gas pressure limited percolation of the supernova remnant heated hot phase, but can run away in gas-rich protodiscs to generate compact starbursts. A simple analytic model is derived for a Schmidt-like global star formation law in terms of the cold gas volume density.     

On the angular momentum in star formation

We discuss the rotation of interstellar clouds which are in a stage immediately before star formation. Cloud collisions seem to be the principal cause of the observed rotation of interstellar clouds. The rotational motion of the clouds is strongly influenced by turbulence.Theories dealing with the resolution of the angular momentum problem in star formation are classified into five major groups. We develop the old idea that the angular momentum of an interstellar cloud passes during star formation into the angular momentum of double star systems and/or circumstellar clouds.It is suggested that a rotating gas cloud contracts into a ring-like structure which fragments into self-gravitating subcondensations. By collisions and gas accretion these subcondensations accrete into binary systems surrounded by circumstellar clouds. Using some rough approximations we find analytical expressions for the semi-major axis of the binary system and for the density of the circumstellar clouds as a function of the initial density and of the initial angular velocity of an interstellar cloud. The obtained values are well within the observational limits.     

The effects of self-gravity on the solar nebula

We study the distribution and transport of angular momentum in a self-gravitating accretion disk formed during the collapse of a rotating gas cloud. Using the surface density for the low-viscosity models and minimum-mass models presented by Cassen and Summers, Poisson's equation is solved explicitly to determine the effects of self-gravitation of the protostellar disk. Analytic expressions for the angular momentum of the central star and other relevant quantities of interest during the formation stage are presented.Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984.On leave from the City College of the City University of New York, U.S.A.     

Triggered planet formation in young stellar clusters

Conventional planet formation models via coagulation of planetesimals require timescales in the range of several 10 or even 100 Myr in the outer regions of a protoplanetary disk. But according to observational data, the lifetime of a protoplanetary disk is limited to about 6 Myr. Therefore the existence of Uranus and Neptune poses a problem. Planet formation via gravitational instability may be a solution for this discrepancy. We present a parameter study of the possibility of gravitationally triggered disk instability. Using a restricted N‐body model which allows for a survey of an extended parameter space, we show that a passing dwarf star with a mass between 0.1 and 1 M_⊙ can probably induce gravitational instabilities in the pre‐planetary solar disk for prograde passages with minimum separations below 80‐170 AU. Inclined and retrograde encounters lead to similar results but require slightly closer passages. Such encounter distances are quite likely in young moderately massive star clusters. The induced gravitational instabilities may lead to enhanced planetesimal formation in the outer regions of the protoplanetary disk, and could therefore be relevant for the formation of Uranus and Neptune. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Does each spiral Galaxy house a slowly growing elliptical one?

N-body simulations performed by us suggest a mechanism for the generation of spiral waves in galaxies in which a mutual quasi-ellipsoidal rotating equilibrium configuration increasing slowly by accretion from the surrounding disk influences the density distribution of stars in the disk such as to give rise to a trailing spiral density wave. Interaction of the spiral wave with the viscous interstellar gas and mutual gravitation between the stars in the disk are believed to influence the form of the spiral. Nevertheless the basic assumption of conventional density wave theory according to which the mutual interaction of stars in the disk is essential for the formation of spirals may not be true.     

Long-term and large-scale viscous evolution of dense planetary rings

Planetary rings are common in the outer Solar System but their origin and long-term evolution is still a matter of debate. It is well known that viscous spreading is a major evolutionary process for rings, as it globally redistributes the disk’s mass and angular momentum, and can lead to the disk’s loosing mass by infall onto the planet or through the Roche limit. However, describing this process is highly dependent on the model used for the viscosity. In this paper we investigate the global and long-term viscous evolution of a circumplanetary disk. We have developed a simple 1D numerical code, but we use a physically realistic viscosity model derived from N-body simulations (Daisaka et al., 2001), and dependent on the disk’s local properties (surface mass density, particle size, distance to the planet). Particularly, we include the effects of gravitational instabilities (wakes) that importantly enhance the disk’s viscosity. This method allows to study the global evolution of the disk over the age of the Solar System.Common estimates of the disk’s spreading time-scales with constant viscosity significantly underestimate the rings’ lifetime. We show that, with a realistic viscosity model, an initially narrow ring undergoes two successive evolutionary stages: (1) a transient rapid spreading when the disk is self-gravitating, with the formation of a density peak inward and an outer region marginally gravitationally stable, and with an emptying time-scale proportional to (where M₀ is the disk’s initial mass), (2) an asymptotic regime where the spreading rate continuously slows down as larger parts of the disk become non-self-gravitating due to the decrease of the surface density, until the disk becomes completely non-self-gravitating. At this point its evolution dramatically slows down, with an emptying time-scale proportional to 1/M₀, which significantly increases the disk’s lifetime compared to the case with constant viscosity. We show also that the disk’s width scales like t^1/4 with the realistic viscosity model, while it scales like t^1/2 in the case of constant viscosity, resulting in much larger evolutionary time-scales in our model. We find however that the present shape of Saturn’s rings looks like a 100 million-years old disk in our simulations. Concerning Jupiter’s, Uranus’ and Neptune’s rings that are faint today, it is not likely that they were much more massive in the past and lost most of their mass due to viscous spreading alone.     

Physics of the primitive solar accretion disk

The theory of viscous accretion disks developed by Lynden-Bell and Pringle has been applied to the evolution of the primitive solar nebula. The additional physical input needed to determine the structure of the disk is described. A series of calculations was carried out using a steady flow approximation to explore the effects on the disk properties of variations in such parameters as the angular momentum and accretion rate of the infalling material from a collapsing interstellar cloud fragment. The more detailed evolutionary calculations involved five cases with various combinations of parameters. It was concluded that the late stages of evolution of the disks would be dominated by the effects of mass loss from the expansion of a hot disk corona into space, and the effects of this were included in the evolutionary calculations. A new theory of comet formation is formulated upon these results. The most important result is the conclusion, which appears to be inescapable, that the primitive solar accretion disk was repeatedly unstable against axisymmetric perturbations, in which rings would form and collapse upon themselves, with the subsequent formation of giant gaseous protoplanets.     

On the formation of protostellar disks

An analysis is presented of the hydrodynamic aspects of the growth of protostellar disks from the accretion (or collapse) of a rotating gas cloud. The size, mass, and radiative properties of protostellar disks are determined by the distribution of mass and angular momentum in the clouds from which they are formed, as well as from the dissipative processes within the disks themselves. The angular momentum of the infalling cloud is redistributed by the action of turbulent viscosity on a shear layer near the surface of the disk (downstream of the accretion shock) and on the radial shear across cylindrical surfaces parallel to the rotation axis. The fraction of gas that is fed into a central core (protostar) during accretion depends on the ratio of the rate of viscous diffusion of angular momentum to the accretion rate; rapid viscous diffusion (or a low accretion rate) promotes a large core-to-disk mass ratio. The continuum radiation spectrum of a highly viscous disk is similar to that of a steady-state accretion disk without mass addition. It is possible to construct models of the primitive solar nebula as an accretion disk, formed by the collapse of a slowly rotating protostellar cloud, and containing the minimum mass required to account for the planets. Other models with more massive disks are also possible.     

The formation of the Hubble sequence of disc galaxies: the effects of early viscous evolution

We investigate a model of disc galaxies whereby viscous evolution of the gaseous disc drives material inwards to form a protobulge. We start from the standard picture of disc formation through the settling of gas into a dark halo potential well, with the disc initially coming into centrifugal equilibrium with detailed conservation of angular momentum. We derive generic analytic solutions for the disc–halo system after adiabatic compression of the dark halo, with free choice of the input virialized dark halo density profile and of the specific angular momentum distribution. We derive limits on the final density profile of the halo in the central regions. Subsequent viscous evolution of the disc is modelled by a variation of the specific angular momentum distribution of the disc, providing analytic solutions to the final disc structure. The assumption that the viscous evolution time-scale and the star formation time-scale are similar leads to predictions of the properties of the stellar components. Focusing on small 'exponential' bulges, i.e., ones that may be formed through a disc instability, we investigate the relationship between the assumed initial conditions, such as halo 'formation', or assembly, redshift z _f, spin parameter λ , baryonic fraction F , and final disc properties such as global star formation time-scale, gas fraction, and bulge-to-disc ratio. We find that the present properties of discs, such as the scalelength, are compatible with a higher initial formation redshift if the redistribution by viscous evolution is included than if it is ignored. We also quantify the dependence of final disc properties on the ratio F λ , thus including the possibility that the baryonic fraction varies from galaxy to galaxy, as perhaps may be inferred from the observations.     

Modeling the luminosity function of galactic low-mass X-ray binaries

The evolution of the family of binaries with a low-mass star and a compact neutron star companion (low-mass X-ray binaries (LMXBs) with neutron stars) ismodeled by the method of population synthesis. Continuous Roche-lobe filling by the optical star in LMXBs is assumed to be maintained by the removal of orbital angular momentum from the binary by a magnetic stellar wind from the optical star and the radiation of gravitational waves by the binary. The developed model of LMXB evolution has the following significant distinctions: (1) allowance for the effect of the rotational evolution of a magnetized compact remnant on themass transfer scenario in the binary, (2) amore accurate allowance for the response of the donor star to mass loss at the Roche-lobe filling stage. The results of theoretical calculations are shown to be in good agreement with the observed orbital period-X-ray luminosity diagrams for persistent Galactic LMXBs and their X-ray luminosity function. This suggests that the main elements of binary evolution, on the whole, are correctly reflected in the developed code. It is shown that most of the Galactic bulge LMXBs at luminosities L _x > 10³⁷ erg s^?1 should have a post-main-sequence Roche-lobe-filling secondary component (low-mass giants). Almost all of the models considered predict a deficit of LMXBs at X-ray luminosities near ～10^36.5 erg s^?1 due to the transition of the binary from the regime of angular momentum removal by a magnetic stellar wind to the regime of gravitational waves (analogous to the widely known period gap in cataclysmic variables, accreting white dwarfs). At low luminosities, the shape of the model luminosity function for LMXBs is affected significantly by their transient behavior-the accretion rate onto the compact companion is not always equal to the mass transfer rate due to instabilities in the accretion disk around the compact object. The best agreement with observed binaries is achieved in the models suggesting that heavy neutron stars with masses 1.4–1.9M _⊙ can be born.     

Energy, angular momentum and wave action associated with density waves in a rotating magnetized gas disc

Both fast and slow magnetohydrodynamic (MHD) density waves propagating in a thin rotating magnetized gas disc are investigated. In the tight-winding or WKBJ regime, the radial variation of MHD density-wave amplitude during wave propagation is governed by the conservation of wave action surface density which travels at a relevant radial group speed C _g. The wave energy surface density and the wave angular momentum surface density are related to by = and = m respectively, where is the angular frequency in an inertial frame of reference and the integer m , proportional to the azimuthal wavenumber, corresponds to the number of spiral arms. Consequently, both wave energy and angular momentum are conserved for spiral MHD density waves. For both fast and slow MHD density waves, net wave energy and angular momentum are carried outward or inward for trailing or leading spirals, respectively. The wave angular momentum flux contains separate contributions from gravity torque, advective transport and magnetic torque. While the gravity torque plays an important role, the latter two can be of comparable magnitudes to the former. Similar to the role of gravity torque, the part of MHD wave angular momentum flux by magnetic torque (in the case of either fast or slow MHD density waves) propagates outward or inward for trailing or leading spirals, respectively. From the perspective of global energetics in a magnetized gas sheet in rotation, trailing spiral structures of MHD density waves are preferred over leading ones. With proper qualifications, the generation and maintenance as well as transport properties of MHD density waves in magnetized spiral galaxies are discussed.     

Fragmentation of a Molecular Cloud Core versus Fragmentation of the Massive Protoplanetary Disk in the Main Accretion Phase

We present three-dimensional numerical simulations on binary formation through fragmentation. The simulations follow gravitational collapse of a molecular cloud core up to growth of the first core by accretion. At the initial stage, the gravity is only slightly dominant over the gas pressure. We made various models by changing initial velocity distribution (rotation speed, rotation law, and bar-mode perturbation). The cloud fragments whenever the cloud rotates sufficiently slowly to allow collapse but faster enough to form a disk before first-core formation. The latter condition is equivalent to Ω₀ t _ff ? 0.05, where Ω₀ and t _ff f denote the initial central angular velocity and the freefall time measured from the central density, and the condition is independent of the initial rotation law and bar-mode perturbation. Fragmentation is classified into six types. When the initial cloud rotates rigidly the cloud collapses to form a adiabatic disk supported by rotation. When the bar-mode perturbation is very minor, the disk deforms to a rotating bar, and the bar fragments. Otherwise, the adiabatic disk evolves into a central core surrounded by a circumstellar disk, and the the circumstellar disk fragments. When the initial cloud rotates differentially, the cloud deforms to a ring or bar in the isothermal collapse phase. The ring fragments into free or more cores, while the bar fragments into only two cores. In the latter case, the core merges due to low orbital angular momentum and new satellite cores form in the later stages.     

Magnetohydrodynamic instabilities and turbulence in accretion disks

In this paper we review the possibilities for magnetohydrodynamic processes to handle the angular momentum transport in accretion disks. Traditionally the angular momentum transport has been considered to be the result of turbulent viscosity in the disk, although the Keplerian flow in accretion disks is linearly stable towards hydrodynamic perturbations. It is on the other hand linearly unstable to some magnetohydrodynamic (MHD) instabilities. The most important instabilities are the Parker and Balbus-Hawley instabilities that are related to the magnetic buoyancy and the shear flow, respectively. We discuss these instabilities not only in the traditional MHD framework, but also in the context of slender flux tubes, that reduce the complexity of the problem while keeping most of the stability properties of the complete problem. In the non-linear regime the instabilities produce turbulence. Recent numerical simulations describe the generation of magnetic fields by a dynamo in the resulting turbulent flow. Eventually such a dynamo may generate a global magnetic field in the disk. The relation of the MHD-turbulence to observations of accretion disks is still obscure. It is commonly believed that magnetic fields can be highly efficient in transporting the angular momentum, but emission lines, short-time scale variability and non-thermal radiation, which a stellar astronomer would take as signs of magnetic variability, are more commonly observed during periods of low accretion rates. Received October 12, 1995 / Accepted November 16, 1995     

Collapsing index:a new method to identify star-forming cores based on ALMA images

Stars form through the gravitational collapse of molecular cloud cores.Before collapsing,the cores are supported by thermal pressure and turbulent motions.A question of critical importance for the understanding of star formation is how to observationally discern whether a core has already initiated gravitational collapse or is still in hydrostatic balance.The canonical method to identify gravitational collapse is based on the observed radial density profile,which would change from Bonnor-Ebert type toward power laws as the core collapses.In practice,due to the projection effect,the resolution limit and other caveats,it has been difficult to directly reveal the dynamical status of cores,particularly in massive star forming regions.We here propose a novel,straightforward diagnostic,namely,the collapsing index(CI),which can be modeled and calculated based on the radial profile of the line width of dense gas.A meaningful measurement of CI requires spatially and spectrally resolved images of optically thin and chemically stable dense gas tracers.ALMA observations are making such data sets increasingly available for massive star forming regions.Applying our method to one of the deepest dense-gas spectral images ever taken toward such a region,namely,the Orion molecular cloud,we detect the dynamical status of selected cores.We observationally distinguished a collapsing core in a massive star forming region from a hydrostatical one.Our approach would help significantly improve our understanding of the interaction between gravity and turbulence within molecular cloud cores in the process of star formation.     

Linear kinetic theory of instabilities of a gravitatingstellar disk

Linear kinetic theory is developed to describe collective oscillations (and their instabilities) propagating in a rapidly rotating disk of stars, representing a highly flattened galaxy. The analysis is carried out for the special case of a self-gravitating, infinitesimally thin, and spatially inhomogeneous system, taking into account the effects both of thermal movements of stars and of gravitational encounters between stars and giant molecular clouds of an interstellar medium. The star–cloud encounters are described with the use of the Landau collision integral. The dynamics of gravity perturbations with rare interparticle encounters is considered. Such a disk is treated by employing the well elaborated mathematical formalisms from plasma perturbation theory using normal-mode analysis. In particular, the method of solving the Boltzmann equation is applied by integration along paths, neglecting the influence of star–cloud encounters on the distribution of stars in the zeroth-order approximation. We are especially interested in important kinetic effects due to wave–star resonances, which we have little knowledge about. The kinetic effects are introduced via a minor drift motion of stars which is computed from the equations of stellar motion in an unperturbed central force field of a galaxy. The dispersion laws for two main branches of disk's oscillations, that is the classical Jeans branch and an additional gradient branch, are deduced. The resonant Landau-type instabilities of hydrodynamically stable Jeans and gradient gravity perturbations is considered to be a long-term generating mechanism for propagating density waves, thereby leading to spiral-like and/or ring-like patterns in the flat galaxies. This revised version was published online in July 2006 with corrections to the Cover Date.     

The role of tidal interactions in star formation

Nearly all of the initial angular momentum of the matter that goes into each forming star must somehow be removed or redistributed during the formation process. The possible transport mechanisms and the possible fates of the excess angular momentum are discussed, and it is argued that transport processes in discs are probably not sufficient by themselves to solve the angular momentum problem, while tidal interactions with other stars in forming binary or multiple systems are likely to be of very general importance in redistributing angular momentum during the star formation process. Most, if not all, stars probably form in binary or multiple systems, and tidal torques in these systems can transfer much of the angular momentum from the gas around each forming star to the orbital motions of the companion stars. Tidally generated waves in circumstellar discs may contribute to the overall redistribution of angular momentum. Stars may gain much of their mass by tidally triggered bursts of rapid accretion, and these bursts could account for some of the most energetic phenomena of the earliest stages of stellar evolution, such as jet-like outflows. If tidal interactions are indeed of general importance, planet-forming discs may often have a more chaotic and violent early evolution than in standard models, and shock heating events may be common. Interactions in a hierarchy of subgroups may play a role in building up massive stars in clusters and in determining the form of the upper initial mass function (IMF) . Many of the processes discussed here have analogues on galactic scales, and there may be similarities between the formation of massive stars by interaction-driven accretion processes in clusters and the buildup of massive black holes in galactic nuclei.