The gravitational collapse of Plummer protostellar clouds

¹ Centro de Investigación en Física, Universidad de Sonora, Blvd. Rosales y Blvd. Encinas s/n Col. Centro, Hermosillo 83000, Sonora, Mexico e-mail: garreaga@cajeme.cifus.uson.mx
² Instituto Nacional de Investigaciones Nucleares, Carretera Mexico-Toluca s/n, La Marquesa, Ocoyoacac 52750, Estado de Mexico, Mexico e-mail: jaime.klapp@inin.gob.mx
³ Facultad de Ciencias, Universidad Autónoma del Estado de Mexico, El Cerrillo Piedras Blancas 50200, Estado de Mexico, Mexico e-mail: fer_7mx@yahoo.com.mx

Received: 2 February 2009
Accepted: 22 October 2009

Abstract

The growing observational evidence that main-sequence and pre-main-sequence protostars are found in binary and multiple systems suggests that they are formed by a fragmentation of collapsing molecular cloud cores. In this paper we present the results of a set of numerical simulations aimed to study the gravitational collapse and fragmentation of a centrally condensed cloud with a nearly flat central density region and a surrounding gas envelope. In order to describe this cloud structure, we use the Plummer radial density profile which satisfies the observed fact that protostellar clouds have a flat central density profile in their innermost region. We consider the cloud to be made of molecular hydrogen and describe the thermodynamics with a single barotropic equation of state, which includes a critical density as a unique free parameter that determines a thermodynamical change on the collapsing gas: from an isothermal to an adiabatic regime. In this paper we consider four different values for the initial radius of the cloud, ranging from 2675 to . In the models, for each ratio of the cloud to core radius, we use two critical density values: and . When the adiabatic change regime starts earlier, we find interesting gas structures as a result of the collapse, although these structures are different according to the initial mass content of the envelope and the initial angular velocity of the cloud. When the thermodynamical change occurs later, i.e., for , we observe that the previously found structure is almost erased to give place instead to a single clump of gas without any adorning spiral arms. In general, we find that as the extension of the envelope mass increases, the possibility of a model to produce a multiple system decreases. This is a result of the initial configuration of our models, namely that with bigger envelopes their cores have a lower ratio of rotational to gravitational energy , a lower ratio of thermal plus rotational to gravitational energy , and a lower angular velocity , which induces a stronger collapse which in turn contributes to the destruction of the structure that is formed during the initial phases of the collapse. Thus in a sufficient quantity rotational energy is crucial for the fragmentation to occur and survive.