Observational evidence and theoretical arguments postulate that outflows may play a significant role in the advection-dominated accretion discs (ADAFs). While the azimuthal viscosity is the main focus of most previous studies in this context, recent studies indicated that disc structure can also be affected by the radial viscosity. In this work, we incorporate these physical ingredients and the toroidal component of the magnetic field to explore their roles in the steady-state structure of ADAFs. We thereby present a set of similarity solutions where outflows contribute to the mass loss, angular momentum removal, and the energy extraction. Our solutions indicate that the radial viscosity causes the disc to rotate with a slower rate, whereas the radial gas velocity increases. For strong winds, the infall velocity may be of order the Keplerian speed if the radial viscosity is considered and the saturated conduction parameter is high enough. We show that the strength of magnetic field and of wind can affect the effectiveness of radial viscosity.

We study the evolution of the protoplanetary discs (PPDs) in the presence of magnetically driven winds with the stress relations motivated by the non-ideal MHD disc simulations. Contribution of the magnetic winds in the angular momentum removal and mass loss are described using these relations which are quantified in terms of the plasma parameter. Evolution of the essential disc quantities including the surface density, accretion rate and wind mass loss rate are studied for a wide range of the model parameters. Two distinct phases of the disc evolution are found irrespective of the adopted input parameters. While at the early phase of the disc evolution, global disc quantities such as its total mass and magnetic flux undergo non-significant reductions, their rapid declines are found in the second phase of evolution. Duration of each phase, however, depends upon the model parameters including magnetic wind strength. Our model predicts that contributions of the magnetic winds in the disc evolution are significant during the second phase. We then calculated locus of points in the plane of the accretion rate and total disc mass corresponding to an ensemble of evolving PPDs. Our theoretical isochrone tracks exhibit reasonable fits to the observed PPDs in star forming regions Lupus and $\sigma$-Orion.

In protoplanetary discs (PPDs) consisting of gas and dust particles, fluid instabilities induced by the drag force, including secular gravitational instability (SGI) can facilitate planet formation. Although SGI subject to the axisymmetric perturbations was originally studied in the absence of gas feedback and it then generalized using a two-fluid approach, the fate of the nonaxisymmetric SGI, in either case, is an unexplored problem. We present a linear perturbation analysis of the nonaxisymmetric SGI in a PPD by implementing a two-fluid model. We explore the growth of the local, nonaxisymmetric perturbations using a set of linearized perturbation equations in a sheared frame. The nonaxisymmetric perturbations display significant growth during a finite time interval even when the system is stable against the axisymmetric perturbations. Furthermore, the surface density perturbations do not show continuous growth but are temporally amplified. We also study cases where the dust component undergoes amplification whereas the gas component remains stable. The amplitude amplification, however, strongly depends on the model parameters. In the minimum mass solar nebula (MMSN), for instance, the dust fluid amplification at the radial distance 100 au occurs when the Stokes number is about unity. But the amplification factor reduces as the dust and gas coupling becomes weaker. Furthermore, perturbations with a larger azimuthal wavelength exhibit a larger amplification factor.

Recent numerical simulations of magnetized accretion discs show that the radial-azimuthal component of the stress tensor due to the magnetorotational instability (MRI) is well represented by a power-law function of the gas pressure rather than a linear relation which has been used in most of the accretion disc studies. The exponent of this power-law function which depends on the net flux of the imposed magnetic field is reported in the range between zero and unity. However, the physical consequences of this power-law stress-pressure relation within the framework of the standard disc model have not been explored so far. In this study, the structure of an accretion disc with a power-law stress-pressure relation is studied using analytical solutions in the steady-state and time-dependent cases. The derived solutions are applicable to different accreting systems, and as an illustrative example, we explore structure of protoplanetary discs using these solutions. We show that the slopes of the radial surface density and temperature distributions become steeper with decreasing the stress exponent. However, if the disc opacity is dominated by icy grains and value of the stress exponent is less than about $0.5$, the surface density and temperature profiles become so steep that make them unreliable. We also obtain analytical solutions for the protoplanetary discs which are irradiated by the host star. Using these solutions, we find that the effect of the irradiation becomes more significant with decreasing the stress exponent.

Lifetimes of protoplanetary disks (PPDs) are believed to be severely constrained by material depleting mechanisms, including photoevaporative winds due to the host star radiation or external radiation sources. Most previous studies focused on exploring the role of the winds in the exposed PPDs with a single star; however, exploring the evolution of the circumbinary disks with the photoevaporative winds driven by the host star radiation and external radiation sources deserves further investigation. In this study, we investigate the evolution of the circumbinary PPDs with the photoevaporative winds induced by external far ultraviolet (FUV) radiation field. We show that this mass-loss process can significantly constrain properties of a circumbinary PPD, including its lifetime, mass and radius. The lifetime of a circumbinary PPD, for instance, is found by a factor of about two longer than a similar circumstellar disk and this enhancement strongly depends on the viscosity parameter. But our model shows that viscosity dependence of the disk lifetime in the circumbinary case is more pronounced compared to the circumstellar case. We also show that dispersal of a circumbinary PPD occurs over a longer time as the disk temperature distribution becomes steeper. Our results also imply that dead zone in a photoevaporative circumbinary PPD extends over a larger radial range in comparison to a circumstellar disk counterpart. We also show that our calculations are in agreement with the observed circumbinary PPDs orbiting equal-mass binaries.

We present an analytical model to investigate the production of pebbles and their radial transport through a protoplanetary disk (PPD) with magnetically driven winds. While most of the previous analytical studies in this context assume that the radial turbulent coefficient is equal to the vertical dust diffusion coefficient, in the light of the results of recent numerical simulations, we relax this assumption by adopting effective parametrisations of the turbulent coefficients involved in terms of the strength of the magnetic fields driving the wind. Theoretical studies have already pointed out that even in the absence of winds, these coefficients are not necessarily equal, though its consequences regarding pebble production have not been explored. In this paper, we investigate the evolution of the pebble production line, the radial mass flux of the pebbles and their corresponding surface density as a function of the plasma parameter at the disk midplane. Our analysis explicitly demonstrates that the presence of magnetically-driven winds in a PPD leads to considerable reduction of the rate and duration of the pebble delivery. We show that when the wind is strong, the core growth in mass due to the pebble accretion is so slow that it is unlikely that a core could reach a pebble isolation mass during a PPD lifetime. When the mass of a core reaches this critical value, pebble accretion is halted due to core-driven perturbations in the gas. With decreasing wind strength, however, pebble accretion may, in a shorter time, increase the mass of a core to the pebble isolation mass.

We implement a new semi-analytical approach to investigate radially self-similar solutions for the steady-state advection-dominated accretion flows (ADAFs). We employ the usual $\alpha$-prescription for the viscosity and all the components of the energy-momentum tensor are considered. In this case, in the spherical coordinate, the problem reduces to a set of eighth-order, nonlinear differential equations with respect to the latitudinal angle $\theta$. Using the Fourier expansions for all the flow quantities, we convert the governing differential equations to a large set of nonlinear algebraic equations for the Fourier coefficients. Using the Newton-Raphson method we solve the algebraic equations and ADAF properties are investigated over a wide range of the model parameters. We also show that the implemented series are truly convergent. The main advantage of our numerical method is that it does not suffer from the usual technical restrictions which may arise for solving ADAF differential equations near to the polar axis. In order to check the reliability of the approach, we recover some widely studied solutions. We also introduce a new varying $\alpha$ viscosity model and new outflow and inflow solutions for ADAFs are presented using Fourier expansion series.

We present a linear analysis of the radiative Rayleigh-Taylor (RT) instability in the presence of magnetic field for both optically thin and thick regimes. When the flow is optically thin, magnetic field not only stabilizes perturbations with short wavelengths, but also growth rate of the instability at long wavelengths is reduced compared to a nonmagnetized case. Then, we extend our analysis to the optically thick flows with a conserved total specific entropy and properties of the unstable perturbations are investigated in detail. Growth rate of the instability at short wavelengths is suppressed due to the presence of the magnetic field, however, growth rate is nearly constant at long wavelengths because of the radiation field. Since the radiative bubbles around massive protostars are subject to the RT instability, we also explore implications of our results in this context. In the nonmagnetized case, the growth time-scale of the instability for a typical bubble is found less than one thousand years which is very short compared to the typical star formation time-scale. Magnetic field with a reasonable strength significantly increases the growth time-scale to more than hundreds of thousands years. The instability, furthermore, is more efficient at large wavelengths, whereas in the non-magnetized case, growth rate at short wavelengths is more significant.

We study local stability of the advection-dominated optically thick (slim) and optically thin discs with purely toroidal magnetic field and the radial viscous force using a linear perturbation analysis. Our dispersion relation indicates that the presence of magnetic fields and radial viscous force cannot give rise to any new mode of the instability. We find, however, that growth rate of the thermal mode in the slim discs and that of the acoustic modes in the slim and optically thin discs are dramatically affected by the radial viscous force. This force tends to strongly decrease the growth rate of the outward-propagating acoustic mode (O-mode) at the short-wavelength limit, but it causes a slim disc to become thermally more unstable. We find that growth rate of the thermal mode increases in the presence of radial viscous force. This enhancement is more significant when the viscosity parameter is large. We also show that growth rate of the O-mode reduces when radial viscous force is considered. The growth rates of the thermal and acoustic modes depend on the magnetic field. Although the instability of O-mode for a stronger magnetic field case has a higher growth rate, the thermal mode of the slim discs can be suppressed when the magnetic field is strong. The inertial-acoustic instability of a magnetized disc may explain the quasi-periodic oscillations (QPOs) from the black holes.

We present a new set of analytical solutions to model the steady state structure of a protoplanetary disc with a magnetically-driven wind. Our model implements a parametrization of the stresses involved and the wind launching mechanism in terms of the plasma parameter at the disc midplane, as suggested by the results of recent, local MHD simulations. When wind mass-loss is accounted for, we find that its rate significantly reduces the disc surface density, particularly in the inner disc region. We also find that models that include wind mass-loss lead to thinner dust layers. As an astrophysical application of our models, we address the case of HL Tau, whose disc exhibits a high accretion rate and efficient dust settling at its midplane. These two observational features are not easy to reconcile with conventional accretion disc theory, where the level of turbulence needed to explain the high accretion rate would prevent a thin dust layer. Our disc model that incorporates both mass-loss and angular momentum removal by a wind is able to account for HL Tau observational constraints concerning its high accretion rate and dust layer thinness.

We obtain a new set of analytical solutions for the evolution of a self-gravitating accretion disc by holding the Toomre parameter close to its threshold, and obtaining the stress parameter from the cooling rate. In agreement with the previous numerical solutions, furthermore, the accretion rate is assumed to be independent of the disc radius. Extreme situations where the entire disc is either optically thick or optically thin are studied independently and the obtained solutions can be used for exploring the early or the final phases of a protoplanetary disc evolution. Our solutions exhibit decay of the accretion rate as a power-law function of the age of the system with the exponent -0.75 and -1.04 for optically thick and thin cases, respectively. Our calculations permit us to explore evolution of the snow line analytically. Location of the snow line in the optically thick regime evolves as a power-law function of time with the exponent -0.16, however, when the disc is optically thin, location of the snow line as a function of time with the exponent -0.7 has a stronger dependence on time. It means that in an optically thin disc inward migration of the snow line is faster than an optically thick disc.

We obtain two-dimensional exact analytic solutions for the structure of the hot accretion flows without wind. We assume that the only non-zero component of the stress tensor is $T_{r\varphi}$. Furthermore we assume that the value of viscosity coefficient $\alpha$ varies with $\theta$. We find radially self-similar solutions and compare them with the numerical and the analytical solutions already studied in the literature. The no-wind solution obtained in this paper may be applied to the nuclei of some cool-core clusters.

Self-gravitating accretion discs in a gravitoturbulent state, including radiation and gas pressures, are studied using a set of new analytical solutions. While the Toomre parameter of the disc remains close to its critical value for the onset of gravitational instability, the dimensionless stress parameter is uniquely determined from the thermal energy reservoir of the disc and its cooling rate. Our solutions are applicable to the accretion discs with dynamically important radiation pressure like in the quasars discs. We show that physical quantities of a gravitoturbulent disc in the presence of radiation are significantly modified compared to solutions with only gas pressure. We show that the dimensionless stress parameter is an increasing function of the radial distance so that its steepness strongly depends on the accretion rate. In a disc without radiation its slope is 4.5, however, we show that in the presence of radiation, it varies between 2 and 4.5 depending on the accretion rate and the central mass. As for the surface density, we find a shallower profile with an exponent -2 in a disc with sub-Eddington accretion rate compared to a similar disc, but without radiation, where its surface density slope is -3 independent of the accretion rate. We then investigate gravitational stability of the disc when the stress parameter reaches to its critical value. In order to self-consistently determine the fragmentation boundary, however, it is shown that the critical value of the stress parameter is a power-law function of the ratio of gas pressure and the total pressure and its exponent is around 1.7. We also estimate the maximum mass of the central black hole using our analytical solutions.

We study drag-driven instability in a protoplanetary disc consisting of a layer of single-sized dust particles which are coupled to the magnetized gas aerodynamically and the particle-to-gas feedback is included. We find a dispersion relation for axisymmetric linear disturbances and growth rate of the unstable modes are calculated numerically. While the secular gravitational instability in the absence of particle-to-gas feedback predicts the dust layer is unstable, magnetic fields significantly amplifies the instability if the Toomre parameter for the gas component is fixed. We also show that even a weak magnetic field is able to amplify the instability more or less irrespective of the dust-gas coupling.

We study stability of a dust layer in a gaseous disc subject to the linear axisymmetric perturbations. Instead of considering single-size particles, however, the population of dust particles is assumed to consist of two grain species. Dust grains exchange momentum with the gas via the drag force and their self-gravity is also considered. We show that the presence of two grain sizes can increase the efficiency of the linear growth of drag-driven instability in the protoplanetary discs. A second dust phase with a small mass, comparing to the first dust phase, would reduce the growth timescale even by a factor of two or more especially when its coupling to the gas is weak. It means that once a certain amount of large dust particles form, even though it is much smaller than that of small dust particles, the dust layer becomes more unstable and dust clumping are accelerated. Thus, presence of dust particles with various sizes must be considered in studies of dust clumping in protoplanetary discs where both large and small dust grains are present.

We study orbit of a pressure-confined cloud in the broad-line region (BLR) of active galactic nuclei (AGNs) when the combined effects of the central gravity and anisotropic radiation pressure and the drag force are considered. Physical properties of the intercloud gas such as its pressure and dynamic viscosity are defined as power-law functions of the radial distance. For a drag force proportional to the relative velocity of a cloud and the background gas, a detailed analysis of the orbits is performed for different values of the input parameters. We also present analytical solutions for a situation where the intercloud pressure is uniform and the viscosity is proportional to the inverse square of the radial distance. Our analytical and numerical solutions demonstrate decay of the orbits because of considering the drag force so that a cloud will eventually fall onto the central region after so called time-of-flight. We found that time-of-flight of a BLR cloud is proportional to the inverse of the dimensionless drag coefficient. We discuss if time-of-flight becomes shorter than the life time of the whole system, then existence of mechanisms for continually forming BLR clouds is needed.

Position of the ice line in a circumbinary disc is determined using a simplified and illustrative model. Main sources of the heat in the energy balance of the disc, i.e. heating by the turbulence, irradiation by the components of the binary and the tidal heating are considered. Our goal is to clarify role of the tidal heating in the position of the ice line. When viscous heating and irradiation of the binary are considered, ice line lies interior to the inner radius of the disc in most of the binaries represented by our parameter survey. But tidal heating significantly extends position of the ice line to a larger radius, so that a smaller fraction of the circumbinaries' population may have ice lines interior to the inner radius of the disc.

The standard Advection-Dominated Accretion Flow (ADAF) is studied using a set of self-similar analytical solutions in the spherical coordinates. Our new solutions are useful for studying ADAFs without dealing with the usual mathematical complexity. We assume the $r\varphi$ component of the stress tensor dominates and the latitudinal component of the velocity is negligible. Moreover, the fluid is incompressible and the solutions are radially self-similar. We show that our analytical solutions display most of the important properties of ADAFs which have already been obtained by the detailed numerical solutions. According to our solutions, the density and the pressure of the flow decreases from the equator to the polar regions and this reduction depends on the amount of the advected energy. We also show analytically that an ADAF tends to a quasi-spherical configuration as more energy is advected with the radial flow.

Using a multi-fluid approach, we investigate streaming and thermal instabilities of the electron-ion plasma with homogeneous cold cosmic rays propagating perpendicular to the background magnetic field. Perturbations are considered to be also across the magnetic field. The back-reaction of cosmic rays resulting in strong streaming instabilities is taken into account. It is shown that for sufficiently short wavelength perturbations, the growth rates can exceed the growth rate of cosmic-ray streaming instability along the magnetic field found by Nekrasov & Shadmehri (2012), which is in its turn considerably larger than the growth rate of the Bell instability (2004). The thermal instability is shown not to be subject to the action of cosmic rays in the model under consideration. The dispersion relation for the thermal instability has been derived which includes sound velocities of plasma and cosmic rays, Alfv\'{e}n and cosmic-ray streaming velocities. The relation between these parameters determines the kind of thermal instability ranging from the Parker (1953) to the Field (1965) instabilities. The results obtained can be useful for a more detailed investigation of electron-ion astrophysical objects such as supernova remnant shocks, galaxy clusters and others including the dynamics of streaming cosmic rays.

Steady-state advective accretion flows in the presence of thermal conduction are studied. All three components of velocity in a spherical coordinates are considered and the flow displays both inflowing and outflowing regions according to our similarity solutions. Thermal conductivity provides latitudinal energy transport and so, the flow rotates more slowly and becomes hotter with increasing thermal conductivity coefficient. We also show that opening angle of the outflow region decreases as thermal conduction becomes stronger.

We study linear theory of the magnetized Rayleigh-Taylor instability in a system consisting of ions and neutrals. Both components are affected by a uniform vertical gravitational field. We consider ions and neutrals as two separate fluid systems where they can exchange momentum through collisions. However, ions have direct interaction with the magnetic field lines but neutrals are not affected by the field directly. The equations of our two-fluid model are linearized and by applying a set of proper boundary conditions, a general dispersion relation is derived for our two superposed fluids separated by a horizontal boundary. We found two unstable modes for a range of the wavenumbers. It seems that one of the unstable modes corresponds to the ions and the other one is for the neutrals. Both modes are reduced with increasing the collision rate of the particles and the ionization fraction. We show that if the two-fluid nature is considered, RT instability would not be suppressed and also show that the growth time of the perturbations increases. As an example, we apply our analysis to the Local clouds which seems to have arisen because of the RT instability. Assuming that the clouds are partially ionized, we find that the growth rate of these clouds increases in comparison to a fully ionized case.

Self-similar and semi-analytical solutions are found for the height-averaged equations govern the dynamical behavior of a polytropic, self-gravitating disk under the effects of winds, around the nascent object. In order to describe time evolution of the system, we adopt a radius dependent mass loss rate, then highlight its importance on both the traditional $\alpha$ and innovative $\beta$ models of viscosity prescription. In agreement with some other studies, our solutions represent that Toomre parameter is less than one in most regions on the $\beta$-disk which indicate that in such disks gravitational instabilities can occur in various distances from the central accretor and so the $\beta$-disk model might provide a good explanation of how the planetary systems form. The purpose of the present work is twofold. First, examining the structure of disk with wind in comparison to no-wind solution; and second, to see if the adopted viscosity prescription affects significantly the dynamical behavior of the disk-wind system. We also considered the temperature distribution in our disk by a polytropic condition. The solutions imply that, under our boundary conditions, the radial velocity is larger for $\alpha$-disks and increases as wind becomes stronger in both viscosity models. Also, we noticed that the disk thickness increases by amplifying the wind or adopting larger values for polytropic exponent $\gamma$. It also may globally decrease if one prescribe $\beta$-model for the viscosity. Moreover, in both viscosity models, surface density and mass accretion rate reduce as wind gets stronger or $\gamma$ increases.

Given the fact that accretion discs are associated with their parent molecular cloud, we studied its effects as a constraint on the outer boundary of the viscous-resistive polytropic self-gravitating accretion flows subject to the mass and angular momentum loss.

Using the multi-fluid approach, we investigate streaming and thermal instabilities of the electron-ion plasma with homogeneous cold cosmic rays drifting perpendicular to the background magnetic field. Perturbations across the magnetic field are considered. The back-reaction of cosmic rays resulting in the streaming instability is taken into account. The thermal instability is shown not to be subject to the action of cosmic rays in the model under consideration. The dispersion relation for the thermal instability has been derived which includes sound velocities of plasma and cosmic rays, Alfv\'{e}n and cosmic ray drift velocities. The relation between these parameters determines the kind of thermal instability from Parker's to Field's type instability. The results obtained can be useful for a more detailed the investigation of electron-ion astrophysical objects such as galaxy clusters including the dynamics of streaming cosmic rays.

When a gravitating point mass moves subsonically through a magnetized and isothermal medium, the dynamical structure of the flow is studied far from the mass using a perturbation analysis. Analytical solutions for the first-order density and the velocity perturbations are presented. Validity of our solutions is restricted to the cases where the Alfven velocity in the ambient medium is less than the accretor's velocity. The density field is less dense because of the magnetic effects according to the solutions and the dynamical friction force becomes lower as the strength of the magnetic field increases.

The purpose of this study is to analyze the dynamical role of a radiation field on the growth rate of the unstable Kelvin - Helmholtz (KH) perturbations. As a first step toward this purpose, the analyze is done in a general way, irrespective of applying the model to a specific astronomical system. The transition zone between the two layers of the fluid is ignored. Then, we perform a linear analysis and by imposing suitable boundary conditions and considering a radiation field, we obtain appropriate dispersion relation. Unstable modes are studied by solving the dispersion equation numerically, and then growth rates of them are obtained. By analyzing our dispersion relation, we show that for a wide range of the input parameters, the radiation field has a destabilizing effect on KH instability. In eruptions of the galaxies or supermassive stars, the radiation field is dynamically important and because of the enhanced KH growth rates in the presence of the radiation; these eruptions can inject more momentum and energy into their environment and excite more turbulent motions.

We investigate the streaming and thermal instabilities of the electron-ion plasma with homogeneous cold cosmic rays drifting perpendicular to the background magnetic field in the multi-fluid approach. One-dimensional perturbations along the magnetic field are considered. The induced return current of the background plasma and back-reaction of cosmic rays are taken into account. It is shown that the cosmic ray back-reaction results in the streaming instability having considerably larger growth rates than that due to the return current of the background plasma. This increase is by a factor of the square root of the ratio of the background plasma mass density to the cosmic ray mass density. The maximal growth rates and corresponding wave numbers are found. The thermal instability is shown to be not subject to the action of cosmic rays in the model under consideration. The dispersion relation for the thermal instability includes ion inertia. In the limit of fast thermal energy exchange between electrons and ions, the isobaric and isochoric growth rates are derived. The results obtained can be useful for the investigation of the electron-ion astrophysical objects such as galaxy clusters including the dynamics of streaming cosmic rays.

Gravitational stability of a disc consisting of the gaseous and the stellar components are studied in the linear regime when the gaseous component is turbulent. A phenomenological approach is adopted to describe the turbulence, in which not only the effective surface density but the velocity dispersion of the gaseous component are both scale-dependent as power-law functions of the wavenumber of the perturbations. Also, the stellar component which has gravitational interaction with the gas is considered as a fluid. We calculate growth rate of the perturbations and in the most of the studied cases, the stability of the disc highly depends on the existence of the stars and the exponents of the functions for describing the turbulence. Our analysis suggests that the conventional gas and star threshold is not adequate for analyzing stability of the two-component discs when turbulence is considered.

We study role of the Hall diffusion in the magnetic star-disc interaction. In a simplified steady state configuration, the total torque is calculated in terms of the fastness parameter and a new term because of the Hall diffusion. We show the total torque reduces as the Hall term becomes more significant. Also, the critical fastness parameter (at which the total torque is zero) reduces because of the Hall diffusion.

Force due to the self-gravity of the disc in the vertical direction is considered to study its possible effects on the structure of a magnetized advection-dominated accretion disc. We present steady-sate self similar solutions for the dynamical structure of such a type of the accretion flows. Our solutions imply reduced thickness of the disc because of the self-gravity. It also imply that the thickness of the disc will increase by adding the magnetic field strength.

We investigate electromagnetic buoyancy instabilities of the electron-ion plasma with the heat flux based on not the magnetohydrodynamic (MHD) equations, but using the multicomponent plasma approach when the momentum equations are solved for each species. We consider a geometry in which the background magnetic field, gravity, and stratification are directed along one axis. The nonzero background electron thermal flux is taken into account. Collisions between electrons and ions are included in the momentum equations. No simplifications usual for the one-fluid MHD-approach in studying these instabilities are used. We derive a simple dispersion relation, which shows that the thermal flux perturbation generally stabilizes an instability for the geometry under consideration. This result contradicts to conclusion obtained in the MHD-approach. We show that the reason of this contradiction is the simplified assumptions used in the MHD analysis of buoyancy instabilities and the role of the longitudinal electric field perturbation which is not captured by the ideal MHD equations. Our dispersion relation also shows that the medium with the electron thermal flux can be unstable, if the temperature gradients of ions and electrons have the opposite signs. The results obtained can be applied to the weakly collisional magnetized plasma objects in laboratory and astrophysics.

We develop a general theory of buoyancy instabilities in the electron-ion plasma with the electron heat flux based not upon MHD equations, but using a multicomponent plasma approach in which the momentum equation is solved for each species. We investigate the geometry in which the background magnetic field is perpendicular to the gravity and stratification. General expressions for the perturbed velocities are given without any simplifications. Collisions between electrons and ions are taken into account in the momentum equations in a general form, permitting us to consider both weakly and strongly collisional objects. However, the electron heat flux is assumed to be directed along the magnetic field that implies a weakly collisional case. Using simplifications justified for an investigation of buoyancy instabilities with the electron thermal flux, we derive simple dispersion relations both for collisionless and collisional cases for arbitrary directions of the wave vector. The collisionless dispersion relation considerably differs from that obtained in the MHD framework and is similar to the Schwarzschild's criterion. This difference is connected with simplified assumptions used in the MHD analysis of buoyancy instabilities and with the role of the longitudinal electric field perturbation which is not captured by the ideal MHD equations. The results obtained can be applied to clusters of galaxies and other astrophysical objects.

The possibility that the stellar initial mass function (IMF) arises mostly from cloud structure is investigated with fractal Brownian motion (fBm) clouds that have power-law power spectra. An fBm cloud with a realistic projected power spectrum slope of $\beta=2.8$ is found to have a mass function for clumps exceeding a threshold density that is a power-law with a slope of $\alpha=2.35$, the same as in the Salpeter IMF. Any hierarchically structured cloud has a clump mass function with about the same slope. This result implies that turbulent interstellar clouds produce dense substructure with the observed pre-stellar core mass function built in from the start. Details of the clump formation processes are not critical. The conversion of clumps into stars involves a second step. A one-to-one correspondence between clump mass and star mass is not necessary to convert the clump mass spectrum into an IMF with the same power-law slope. As long as clumps have an internal stellar IMF from sub-fragmentation, protostellar accretion, coalescence and other processes, and the characteristic mass for this internal IMF scales with the clump mass, then the IMF slope above the minimum characteristic mass will equal the clump mass slope. A detailed review of IMF models illustrates the prominence of cloud structure as a major component in a wide class of theories. Tests are proposed to determine the relative importance of cloud structure and competitive accretion in the IMF.

We study Kelvin-Helmholtz (KH) instability at the interface of a disc and corona system by doing a linear perturbation analysis. The disc is assumed to be thin, however, the corona is considered to be nearly quasispherical because of its high temperature. Under these circumstances, the interface is subject to the KH instability for a given set of the input parameters. Growth rates of the KH unstable modes are calculated for a wide range of the input parameters. We show that for a certain range of the perturbations, the unstable KH perturbations are growing with time scales comparable to the inverse of the angular velocity of the accretion disc (dynamical time scale). Thus, KH instability at the interface of a disc-corona may have enough time to affect the dynamical structure of its underlying accretion disc by possible exchange of the mass, angular momentum or even energy. Our linear analysis shows that KH instability may provide a mechanism for such exchanges between a disc and its corona.

We develop a theory of buoyancy instabilities of the electron-ion plasma with the heat flux based on not the MHD equations, but using the multicomponent plasma approach. We investigate a geometry in which the background magnetic field, gravity, and stratification are directed along one axis. No simplifications usual for the MHD-approach in studying these instabilities are used. The background electron thermal flux and collisions between electrons and ions are included. We derive the simple dispersion relation, which shows that the thermal flux perturbation generally stabilizes an instability. There is a narrow region of the temperature gradient, where an instability is possible. This result contradicts to a conclusion obtained in the MHD-approach. We show that the reason of this contradiction is the simplified assumptions used in the MHD analysis of buoyancy instabilities and the role of the longitudinal electric field perturbation, which is not captured by the MHD equations. Our dispersion relation also shows that a medium with the electron thermal flux can be unstable, if the temperature gradient of ions and electrons have the opposite signs. The results obtained can be applied to ICM and clusters of galaxies.

(abridged) We develop a model which describes the coevolution of the mass function of dense cores and of the IMF in a protocluster clump. In the model, cores injected in the clump evolve under the effect of gas accretion. Accretion onto the cores follows a time-dependent accretion rate that describes accretion in a turbulent medium. Once the accretion timescales of cores exceed their contraction timescales, they are turned into stars. We include the effect of feedback by the newly formed massive stars through their stellar winds. A fraction of the wind's energy is assumed to counter gravity and disperse the gas from the protocluster and as a consequence, quench further star formation. The latter effect sets the final IMF of the cluster. We apply our model to a clump that is expected to resemble the progenitor clump of the Orion Nebula Cluster (ONC). Our model is able to reproduce both the shape and normalization of the ONC's IMF and the mass function of dense cores in Orion. The complex features of the ONC's IMF,i.e., a shallow slope in the mass range ~0.3-2.5 Msol,a steeper slope in the mass range ~2.5-12 Msol, and a nearly flat tail at the high mass end are reproduced. The model predicts a 'rapid' star formation process with an age spread for the stars of 2.3 10^5 yr which is consistent with the fact that 80% of the ONC's stars have ages of <=0.3 Myr. The model predicts a primordial mass segregation with the most massive stars being born in the region between 2-4 times the core radius of the cluster. In parallel, the model also reproduces, simultaneously, the mass function of dense cores in Orion. We study the effects of varying the model parameters on the resulting IMF and show that the IMF of stellar clusters is expected to show significant variations, provided variations in the clumps and cores properties exist.

We study magnetothermal instability in the ionized plasmas including the effects of Ohmic, ambipolar and Hall diffusion. Magnetic field in the single fluid approximation does not allow transverse thermal condensations, however, non-ideal effects highly diminish the stabilizing role of the magnetic field in thermally unstable plasmas. Therefore, enhanced growth rate of thermal condensation modes in the presence of the diffusion mechanisms speed up the rate of structure formation.

We study dynamical effects of cosmics rays (CRs) on the thermal instability in the linear regime. CRs and the thermal plasma are treated as two different interacting fluids, in which CRs can diffuse along the magnetic field lines. We show that growth rate of the magnetothermal condensation mode is reduced because of the existence of CRs and this stabilizing effect depends on the diffusion coefficient and the ratio of the CRs pressure to gas pressure. Thus, a slower rate of structure formation via thermal instability is predicted when CRs are considered.

We study thermal instability in a magnetized and partially ionized plasma with charged dust particles. Our linear analysis shows that the growth rate of the unstable modes in the presence of dust particles strongly depends on the ratio of the cooling rate and the modified dust-cyclotron frequency. If the cooling rate is less than the modified dust-cyclotron frequency, then growth rate of the condensation modes does not modify due to the existence of the charged dust particles. But when the cooling rate is greater than (or comparable to) the modified dust-cyclotron frequency, the growth rate of unstable modes increases because of the dust particles. Also, wavenumber of the perturbations corresponding to the maximum growth rate shifts to the smaller values (larger wavelengths) as the cooling rate becomes larger than the modified dust-cyclotron frequency. We show that growth rate of the condensation modes increases with the electrical charge of the dust particles.

We study effects of winds on the time evolution of isothermal, self-gravitating accretion discs by adopting a radius dependent mass loss rate because of the existence of the wind. Our similarity and semi-analytical solution describes time evolution of the system in the slow accretion limit. The disc structure is distinct in the inner and outer parts, irrespective of the existence of the wind. We show that existence of wind will lead to a reduction of the surface density in the inner and outer parts of the disc in comparison to a no-wind solution. Also, the radial velocity significantly increases in the outer part of the disc, however, the accretion rate decreases due to the reduced surface density in comparison to the no-wind solution. In the inner part of the disc, mass loss due to the wind is negligible according to our solution. But the radial size of this no-wind inner region becomes smaller for stronger winds.

We present self-similar solutions for advection -dominated accretion flows with thermal conduction in the presence of outflows. Possible effects of outflows on the accretion flow are parametrized and a saturated form of thermal conduction, as is appropriate for the weakly-collisional regime of interest, is included in our model. While the cooling effect of outflows is noticeable, thermal conduction provides an extra heating source. In comparison to accretion flows without winds, we show that the disc rotates faster and becomes cooler because of the angular momentum and energy flux which are taking away by the winds. But thermal conduction opposes the effects of winds and not only decreases the rotational velocity, but increases the temperature. However, reduction of the surface density and the enhanced accretion velocity are amplified by both of the winds and the thermal conduction. We find that for stronger outflows, a higher level of saturated thermal conduction is needed to significantly modify the physical profiles of the accretion flow.

We investigate the linear theory of Kelvin-Helmholtz instability at the interface between a partially ionized dusty outflow and the ambient material analytically. We model the interaction as a multifluid system in a planar geometry. The unstable modes are independent from the charge polarity of the dust particles. Although our results show a stabilizing effect for charged dust particles, the growth time scale of the growing modes gradually becomes independent of the mass or charge of the dust particles when the magnetic field strength increases. We show that growth time scale decreases with increasing the magnetic field. Also, as the mass of the dust particles increases, the growth time scale of the unstable mode increases.

We study the steady-state structure of an accretion disc with a corona surrounding a central, rotating, magnetized star. We assume that the magneto-rotational instability is the dominant mechanism of angular momentum transport inside the disc and is responsible for producing magnetic tubes above the disc. In our model, a fraction of the dissipated energy inside the disc is transported to the corona via these magnetic tubes. This energy exchange from the disc to the corona which depends on the disc physical properties is modified because of the magnetic interaction between the stellar magnetic field and the accretion disc. According to our fully analytical solutions for such a system, the existence of a corona not only increases the surface density but reduces the temperature of the accretion disc. Also, the presence of a corona enhances the ratio of gas pressure to the total pressure. Our solutions show that when the strength of the magnetic field of the central neutron star is large or the star is rotating fast enough, profiles of the physical variables of the disc significantly modify due to the existence of a corona.

We study the effect of an imposed magnetic field on the motion of charged dust particles in magnetically active regions of a protoplanetary disc. Assuming a power law structure for the vertical and the toroidal components of the magnetic field for the regions beyond magnetically dead region of the disc, the radial and the vertical velocities of the charged particles, in the asymptotic case of small particles, are calculated analytically. While grains with radii smaller than a critical radius significantly are affected by the magnetic force, motion of the particles with larger radii is independent of the magnetic field. The critical radius depends on the magnetic geometry and the charge of the grains. Assuming that a grain particle has one elementary charge and the physical properties of the disc correspond to a minimum-mass solar nebula, we show that only micron-sized grains are affected by the magnetic force. Also, charge polarity determines direction of the radial velocity. For such small particles, both the radial and the vertical velocities increase due to the magnetic force.

We present a model to explain the mass segregation and shallow mass functions observed in the central parts of dense and young starburst stellar clusters. The model assumes that the initial pre-stellar cores mass function resulting from the turbulent fragmentation of the proto-cluster cloud is significantly altered by the cores coalescence before they collapse to form stars. With appropriate, yet realistic parameters, this model based on the competition between cores coalescence and collapse reproduces the mass spectra of the well studied Arches cluster. Namely, the slopes at the intermediate and high mass ends are reproduced, as well as the peculiar bump observed at 6 M_sol. This coalescence-collapse process occurs on short timescale of the order of one fourth the free fall time of the proto-cluster cloud (i.e., a few 10^{4} years), suggesting that mass segregation in Arches and similar clusters is primordial. The best fitting model implies the total mass of the Arches cluster is 1.45 10^{5} M_sol, which is slightly higher than the often quoted, but completeness affected, observational value of a few 10^{4} M_sol. The derived star formation efficiency is ~30 percent which implies that the Arches cluster is likely to be gravitationally bound.

We study the linear theory of Kelvin-Helmholtz instability in a layer of ions and neutrals with finite thickness. In the short wavelength limit the thickness of the layer has a negligible effect on the growing modes. However, perturbations with wavelength comparable to layer's thickness are significantly affected by the thickness of the layer. We show that the thickness of the layer has a stabilizing effect on the two dominant growing modes. Transition between the modes not only depends on the magnetic strength, but also on the thickness of the layer.

(Abridged) We present numerical hydrodynamical simulations of the formation, evolution and gravitational collapse of isothermal molecular cloud cores. A compressive wave is set up in a constant sub-Jeans density distribution of radius r = 1 pc. As the wave travels through the simulation grid, a shock-bounded spherical shell is formed. The inner shock of this shell reaches and bounces off the center, leaving behind a central core with an initially almost uniform density distribution, surrounded by an envelope consisting of the material in the shock-bounded shell, with a power-law density profile that at late times approaches a logarithmic slope of -2 even in non-collapsing cases. The resulting density structure resembles a quiescent core of radius < 0.1 pc, with a Bonnor-Ebert-like (BE-like) profile, although it has significant dynamical differences: it is initially non-self-gravitating and confined by the ram pressure of the infalling material, and consequently, growing continuously in mass and size. With the appropriate parameters, the core mass eventually reaches an effective Jeans mass, at which time the core begins to collapse. Thus, there is necessarily a time delay between the appearance of the core and the onset of its collapse, but this is not due to the dissipation of its internal turbulence as it is often believed. These results suggest that pre-stellar cores may approximate Bonnor-Ebert structures which are however of variable mass and may or may not experience gravitational collapse, in qualitative agreement with the large observed frequency of cores with BE-like profiles.

We investigate the time evolution of the mass distribution of pre-stellar cores (PSCs) and their transition to the initial stellar mass function (IMF) in the central parts of a molecular cloud (MC) under the assumption that the coalescence of cores is important. Our aim is to explain the observed shallow IMF in dense stellar clusters such as the Arches cluster. The initial distributions of PSCs at various distances from the MC center are those of gravitationally unstable cores resulting from the gravo-turbulent fragmentation of the MC. As time evolves, there is a competition between the PSCs rates of coalescence and collapse. Whenever the local rate of collapse is larger than the rate of coalescence in a given mass bin, cores are collapsed into stars. With appropriate parameters, we find that the coalescence-collapse model reproduces very well all the observed characteristics of the Arches stellar cluster IMF; Namely, the slopes at high and low mass ends and the peculiar bump observed at ~5-6 M_sol. Our results suggest that today's IMF of the Arches cluster is very similar to the primordial one and is prior to the dynamical effects of mass segregation becoming important

We study the dynamical structure of a self-gravitating disc with corona around a super-massive black hole. Assuming that the magneto-rotational-instability (MRI) responsible for generating the turbulent stresses inside the disc is also the source for a magnetically dominated corona, a fraction of the power released when the disc matter accretes is transported to and dissipated in the corona. This has major effect on the structure of the disc and its gravitational (in)stability according to our analytical and self-consistent solutions. We determine the radius where the disc crosses the inner radius of gravitational instability and forms the first stars. Not only the location of this radius which may extend to very large distances from the central black hole, but also the mass of the first stars highly depends on the input parameters, notably the viscous coefficient, mass of the central object and the accretion rate. For accretion discs around quasi-stellar objects (QSOs) and the Galactic center, we determine the self-gravitating radius and the mass of the first clumps. Comparing the cases with a corona and without a corona for typical discs around QSOs or the Galactic center, when the viscosity coefficient is around 0.3, we show that the self-gravitating radius decreases by a factor of approximately 2, but the mass of the fragments increases with more or less the same factor. Existence of a corona implies a more gravitationally unstable disc according to our results. The effect of a corona on the instability of the disc is more effective when the viscosity coefficient increases.

We study the instantaneous virial balance of clumps and cores (CCs) in 3D simulations of driven, MHD, isothermal molecular clouds (MCs). The models represent a range of magnetic field strengths in MCs from subcritical to non-magnetic regimes. We identify CCs at different density thresholds, and for each object, we calculate all the terms that enter the Eulerian form of the virial theorem (EVT). A CC is considered gravitationally bound when the gravitational term in the EVT is larger than the amount for the system to be virialized, which is more stringent than the condition that it be large enough to make the total volume energy negative. We also calculate, quantities commonly used in the observations to indicate the state of gravitational boundedness of CCs such as the Jeans number J_c, the mass-to magnetic flux ratio mu_c, and the virial parameter alpha_vir. Our results show that: a) CCs are dynamical out-of-equilibrium structures. b) The surface energies are of the same order than their volume counterparts c) CCs are either in the process of being compressed or dispersed by the velocity field. Yet, not all CCs that have a compressive net kinetic energy are gravitationally bound. d) There is no 1-to-1 correspondence between the state of gravitational boundedness of a CC as described by the virial analysis or as implied by the classical indicators. In general, in the virial analysis, we observe that only the inner regions of the objects are gravitationally bound, whereas J_c, alpha_vir, and mu_c estimates tend to show that they are more gravitationally bound at the lowest threshold levels and more magnetically supercritical. g) We observe, in the non-magnetic simulation, the existence of a bound core with structural and dynamical properties that resemble those of the Bok globule Barnard 68 (B68).

We study the nonlinear dynamics of a warped or twisted accretion disc, in which the viscosity coefficients are assumed to be locally proportional to the rotational velocity (beta- prescription). Using asymptotic methods for thin discs, dynamical equations of the disc are obtained in a warped spherical polar coordinates. These equations are solved by the method of the separation of the variables. This analysis constitutes an analogous study of the nonlinear theory of an alpha model warped disc which has been studied by Ogilvie (1999). We have compared our results with Ogilvie's analysis. The dynamical behaviors of these models have also been discussed. Our results show that different viscosity prescriptions and magnitudes (alpha and beta prescriptions) affect dynamics of a warped accretion disc. Therefore it can be important in determining the viscosity law even for a warped disc.

The steady-state structure of self-gravitating, magnetized accretion disks is studied using a set of self-similar solutions which are appropriate in the outer regions. The disk is assumed to be isothermal and the magnetic field outside of the disk is treated in a phenomenological way. However, the internal field is determined self-consistently. The behaviour of the solutions are investigated by changing the input parameters of the model, i.e. mass accretion rate, coefficients of viscosity and resistivity, and the magnetic field configuration.

Thermal conduction has been suggested as a possible mechanism by which sufficient extra heating is provided in radiation-dominated accretion flows. We consider the extreme case in which the generated energy due to the viscosity and the energy transported by a saturated form of thermal conduction are balanced by the advection cooling. For the steady-state structure of such accretion flows a set of self-similar solution are presented. Based on these solutions while the radial and the rotational velocities are both sub-Keplerian, increasing the level of thermal conduction has the effects of decreasing the rotational velocity, but increasing the radial velocity. Conduction provides extra heating and the temperature of the gas increases with thermal conduction.

We investigate the effects of a large-scale magnetic field with open field lines on the steady-state structure of a radiation-dominated accretion disk, using self-similarity technique. The disk is supposed to be turbulent and possesses an effective viscosity and an effective magnetic diffusivity. We consider the extreme case in which the generated energy due to viscous and magnetic dissipation is balanced by the advection cooling. While the magnetic field outside of the disk is treated in a phenomenological way, the internal field is determined self-consistently. Magnetized and nonmagnetized solutions have the same radial dependence, irrespective of the values of the input parameters. Generally, our self-similar solutions are very sensitive to the viscosity or diffusivity coefficients. For example, the density and the rotation velocity increase when the viscosity coefficient decreases. The gas rotates with sub-Keplerian angular velocity with a factor less than unity which depends on the magnetic field configuration. Magnetic field significantly reduce disk thickness, however, tends to increase the radial velocity comparing to the nonmagnetic self-similar solutions.

When the gas of a magnetized filamentary cloud obeys a polytropic equation of state, gravitational collapse of the cloud is studied using a simplified model. We concentrate on the radial distribution and restrict ourselves to the purely toroidal magnetic field. If the axial motions and poloidal magnetic fields are sufficiently weak, we could reasonably expect our solutions to be a good approximation. We show that while the filament experiences gravitational condensation and the density at the center increases, the toroidal flux-to-mass ratio remains constant. A series of spatial profiles of density, velocity and magnetic field for several values of the toroidal flux-to-mass ratio and the polytropic index, is obtained numerically and discussed.

We discuss the lifetimes and evolution of clumps and cores formed as turbulent density fluctuations in nearly isothermal molecular clouds. In the non-magnetic case, clumps are unlikely to reach a hydrostatic state, and instead are expected to either proceed directly to collapse, or else ``rebound'' towards the mean pressure and density of the parent cloud. Rebounding clumps are delayed in their re-expansion by their self-gravity. From a simple virial calculation, we find re-expansion times of a few free-fall times. In the magnetic case, we present a series of driven-turbulence, ideal-MHD isothermal numerical simulations in which we follow the evolution of clumps and cores in relation to the magnetic criticality of their ``parent clouds'' (the numerical boxes). In subcritical boxes, magnetostatic clumps do not form. A few moderately-gravitationally bound clumps form which however are dispersed by the turbulence in < 1.3 Myr. An estimate of the ambipolar diffusion (AD) time scale t_AD in these cores gives t_AD > 1.3 Myr, only slightly longer than the dynamical times. In supercritical boxes, some cores become locally supercritical and collapse in typical times ~ 1 Myr. We also observe longer-lived supercritical cores that however do not collapse because they are smaller than the local Jeans length. Fewer clumps and cores form in these simulations than in their non-magnetic counterpart. Our results suggest that a) A fraction of the cores may not form stars, and may correspond to some of the observed starless cores. b) Cores may be out-of-equilibrium structures, rather than quasi-magnetostatic ones. c) The magnetic field may help reduce the star formation efficiency by reducing the probability of core formation, rather than by significantly delaying the collapse of individual cores.

A coalescence model using the observed properties of pre-stellar condensations (PSCs) shows how an initially steep IMF that might be characteristic of primordial cloud fragmentation can change into a Salpeter IMF or shallower IMF in a cluster of normal density after one dynamical time, even if the PSCs are collapsing on their own dynamical time. The model suggests that top-heavy IMFs in some starburst clusters originate with PSC coalescence.

A self-similar solution for time evolution of quasi-spherical, self-gravitating accretion flows is obtained under the assumption that the generated heat by viscosity is retained in the flow. The solutions are parameterized by the ratio of the mass of the accreting gas to the central object mass and the viscosity coefficient. While the density and the pressure are obtained simply by solving a set of ordinary differential equations, the radial and the rotational velocities are presented analytically. Profiles of the density and the rotational velocities show two distinct features. Low density outer accreting flow with relatively flat rotation velocity, surrounds an inner high density region. In the inner part, the rotational velocity increases from the center to a transition radius where separates the inner and outer portions. We show that the behaviour of the solutions in the inner region depends on the ratio of the heat capacities, $\gamma$, and the viscosity coefficient, $\alpha$.

A model for axisymmetric magnetized accretion flow is proposed. The dominant mechanism of energy dissipation is assumed to be the magnetic diffusivity due to turbulence in the accretion flow. In analogy to the advection-dominated accretion flow (ADAF) solutions, a constant fraction of the resistively dissipated energy is stored in the accreting gas and the rest is radiated. The general self-similar solutions which describe a resistive and {\it nonrotating} flow with purely poloidal magnetic field is first introduced. The radial dependence of physical quantities is identical to viscous ADAF solutions. Although the main focus of this study is nonrotating magnetized accretion flow, for rotating flow with both poloidal and toroidal components of magnetic field we find the radial scaling of solutions similar to the nonrotating case. We show that the accretion and the rotation velocities are both below the Keplerian rate, irrespective of the amount of cooling. We show that the set of the equations is reduced to a second order differential equation for a nonrotating flow. The geometrical shape of the disk changes depending on the fraction of the resistively dissipated energy which is stored in the accreting gas. However, there is a hot low-density gas above the disk in almost all cases. The net accretion rate is calculated for a set of illustrative parameters.

We discuss the lifetimes and evolution of dense cores formed as turbulent density fluctuations in magnetized, isothermal molecular clouds. We consider numerical simulations in which we measure the cores' magnetic criticality and Jeans stability in relation to the magnetic criticality of their ``parent clouds'' (the numerical boxes). In subcritical boxes, dense cores do not form, and collapse does not occur. In supercritical boxes, some cores collapse, being part of larger clumps that are supercritical from the start, and whose central, densest regions (the cores) are initially subcritical, but rapidly become supercritical, presumably by accretion along field lines. Numerical artifacts are ruled out. The time scales for cores to go from subcritical to supercritical and then collapse are a few times the core free-fall time, $\tfc$. Our results suggest that cores are out-of-equilibrium, transient structures, rather than quasi-magnetostatic configurations.

Under the assumptions that molecular clouds are nearly spatially and temporally isothermal and that the density peaks (``cores'') within them are formed by turbulent fluctuations, we argue that cores cannot reach a hydrostatic (or magneto-static) state as a consequence of their formation process. In the non-magnetic case, stabilization requires a Bonnor-Ebert truncation at a finite radius, which is not feasible for a single-temperature flow, unless it amounts to a shock, which is clearly a dynamical feature. Instead, in this case, cores must be dynamical entities that can either be pushed into collapse, or else ``rebound'' towards the mean pressure and density as the parent cloud. Nevertheless, rebounding cores are delayed in their re-expansion by their own self-gravity. We give a crude estimate for the re-expansion time as a function of the closeness of the final compression state to the threshold of instability, finding typical values $\sim 1$ Myr, i.e., of the order of a few free-fall times. Our results support the notion that not all cores observed in molecular clouds need to be on route to forming stars, but that instead a class of ``failed cores'' should exist, which must eventually re-expand and disperse, and which can be identified with observed starless cores. In the magnetic case, recent observational and theoretical work suggests that all cores are critical or supercritical, and are thus qualitatively equivalent to the non-magnetic case. Our results support the notion that the entire star formation process is dynamical, with no intermediate hydrostatic stages.

The mm-wave continuum sources (MCS) in Ophiuchus have mutual collision rates less than their collapse rates by a factor of 10 to 100, suggesting most will form stars without further interactions. However, this ratio of rates would have exceeded unity in the past if the MCS were 2.5 times larger than they are now. Such a high previous ratio was not likely, suggesting that the MCS grew from smaller sizes at a constant high density, possible behind turbulent shocks. Proximity constraints are even more important in massive clusters, where massive stars with the same or greater space density are more strongly interacting than the Ophiuchus MCS. As a result, the density contrast for MCS must be larger in massive clusters than it is in Ophiuchus, or significant coalescence will occur in the protostellar phase, possible forming a massive black hole. A proportionality to the second power of the Mach number allows the MCS cores to collapse independently. These results suggest that stars in dense clusters generally form on a dynamical time by the continuous collection and rapid collapse of turbulence-shocked gas. Implications of proximity constraints on the initial stellar mass function are also discussed.

We obtain self-similar solutions that describe the dynamics of a self-gravitating, rotating, viscous system. We use simplifying assumptions; but explicitly include viscosity and the cooling due to the dissipation of energy. By assuming that the turbulent dissipation of energy is as power law of the density and the speed v_{rms} and for a power-law dependence of viscosity on the density, pressure, and rotational velocity, we investigate turbulent cooling flows. It has been shown that for the cylindrically and the spherically cooling flows the similarity indices are the same, and they depend only on the exponents of the dissipation rate and the viscosity model. Depending on the values of the exponents, which the mechanisms of the dissipation and viscosity determine them, we may have solutions with different general physical properties. The conservation of the total mass and the angular momentum of the system strongly depends on the mechanisms of energy dissipation and the viscosity model.

We discuss the virial balance of all members of a cloud ensemble in numerical simulations of self-gravitating MHD turbulence. We first discuss the choice of reference frame for evaluating the terms entering the virial theorem (VT), concluding that the balance of each cloud should be measured in its own reference frame. We then report preliminary results suggesting that a) the clouds are far from virial equilibrium, with the ``geometric'' (time derivative) terms dominating the VT. b) The surface terms in the VT are as important as the volume ones, and tend to decrease the action of the latter. c) This implies that gravitational binding should be considered including the surface terms in the overall balance.

The self-similar equilibrium models of self-gravitating, rotating, isothermal systems are investigated analytically. In these models the rotation velocity is constant and the density varies as $\frac{f(\theta, \phi)}{r^2}$, where $r$ and $\theta$ are the spherical radius and the co-latitude, respectively. The nonaxisymmetric solutions contain three free parameters, one of the parameters depends on the rotation velocity. These parameters determine the overall shape of the density distribution. By assuming that the dominant component of the magnetic field is purely toroidal and the ratio of the purely toroidal magnetic pressure to the gas pressure, $\alpha$, is spatially constant, the axisymmetric solutions generalized so as the effect of magnetic field could be studied. We find that the equilibria of axially symmetric systems yield ellipsoids or spheres only when the ratio of rotation velocity to the sound speed is taken to be $\sqrt{2\alpha}$.

We study the dynamics of a self-gravitating cooling filamentary cloud using a simplified model. We concentrate on the radial distribution and restrict ourselves to quasi-hydrostatic, cylindrically symmetric cooling flows. For a power-law dependence of cooling function on the temperature, self-similar solutions which describe quasi-hydrostatic cooling flows are derived. We consider obtically thin filaments with a constant mass per unit length and the solutions are parameterized by their line masses. There is no polytropic relation between the density and the pressure. The filament experiences radiative condensation, irrespective of the \gamma, the gas specific heat ratio. So, the filament becomes denser due to the quasi-hydrostatic flows and the density at the center increases in proportion to (t_{0}-t)^{-1}, where t denotes the time. The term, t_{0}, denotes an epoch at which the central density increases infinitely. We also found that the radius of the filament decreases in proportion to (t_{0}-t)^{0.5}.