We implement a numerical model reported in the literature to simulate the evolution of a galaxy composed of four matter components, such as: a dark-matter halo; a rotating disk of stars; a spherical bulge of stars and a ring of molecular gas. We show that the evolution of this galaxy model is stable at least for 10 Gyr (Gyr=$10^{9}$ years). We characterize the resulting configuration of this galaxy model by figures of the circular velocity and angular momentum distribution; the tangential and radial components of the velocity; the peak density evolution and the radial density profile. Additionally, we calculate several models of equal-mass galaxy binary collisions, such as: (i) frontal and (ii) oblique (with an impact parameter), (iii) two models with initial conditions taken from a 2-body orbits and (iv) a very close passage. To allow comparison with the galaxy model, we characterize the dynamics of the collision models in an analogous way. Finally, we determine the de Vaucouleurs fitting curves of the radial density profile, on a radial scale of 0-100 kpc, for all the collision models irrespective of the pre-collision trajectory. To study the radial mass density and radial surface density profiles at a smaller radial scale, 0-20 kpc, we use a four-parameters fitting curve.

High-resolution hydrodynamical simulations are presented to follow the gravitational collapse of a uniform turbulent clump, upon which a purely radial compressive velocity pulse is activated in the midst of the evolution of the clump, when its turbulent state has been fully developed. The shape of the velocity pulse is determined basically by two free parameters: the velocity $V_0$ and the initial radial position $r_0$. In the present paper, models are considered in which the velocity $V_0$ takes the values 2, 5, 10, 20, and 50 times the speed of sound of the clump $c_0$, while $r_0$ is fixed for all the models. The collapse of the model with $2 \, c_0$ goes faster as a consequence of the velocity pulse, while the cluster formed in the central region of the isolated clump mainly stays the same. In the models with greater velocity $V_0$, the evolution of the isolated clump is significantly changed, so that a dense shell of gas forms around $r_0$ and moves radially inward. The radial profile of the density and velocity as well as the mass contained in the dense shell of gas are calculated, and it is found that (i) the higher the velocity $V_0$, the less mass is contained in the shell; (ii) there is a critical velocity of the pulse, around $10 \, c_0$, such that for shock models with a lower velocity, there will be a well defined dense central region in the shocked clump surrounded by the shell.

We present high-resolution hydrodynamical simulations aimed at following the gravitational collapse of a gas core, in which a turbulent spectrum of velocity is implemented only initially. We determine the maximal value of the ratio of kinetic energy to gravitational energy, denoted here by $\left(\frac{E_{\rm kin} }{E_{\rm grav}}\right)_{\rm max}$, so that the core (i) will collapse around one free-fall time of time evolution or (ii) will expand unboundedly, because it has a value of $\frac{E_{\rm kin}}{E_{\rm grav}}$ larger than $\left( \frac{E_{\rm kin}}{E_{\rm grav}}\right)_{\rm max}$. We consider core models with a uniform or centrally condensed density profile and with velocity spectra composed of a linear combination of one-half divergence-free turbulence type and the other half of a curl-free turbulence type. We show that the outcome of the core collapse are protostars forming either (i) a multiple system obtained from the fragmentation of filaments and (ii) a single primary system within a long filament. In addition, some properties of these protostars are also determined and compared with those obtained elsewhere.

We consider 28 particle-based simulations aimed at comparing the gravitational collapse of a spherically symmetric, uniform gas core in which two extreme types of turbulent spectra of velocity have been initially induced, so that $\nabla \cdot \vec{v}=0 $ (14 simulations) and $\nabla \times \vec{v}=0 $ (14 simulations). For all the simulations, the ratios of the kinetic energy and thermal energy to the gravitational energy were fixed at $\beta$=0.21 and $\alpha$=0.24, respectively. Most of the simulations finish by forming a single protostar, except for two simulations that form a binary system of protostars. In order to quantify the differences (or similarities) between the two types of simulations, we calculate some integral properties of the resulting protostars, such as the mass M$_f$ and the ratios $\alpha_f$ and $\beta_f$

In this paper we continue the investigation reported by \cite{RMAA} concerning the morphology of binary configurations obtained via the collapse of rotating parent gas structures with total masses in the range of M$_T$= 1 to 5 M$_{\odot}$. Here we extend the mass range and consider the collapse of two uniform gas clumps of M$_T$= 50 and 400 M$_{\odot}$, so that they also rotates rigidly in such a way that its approximate virial parameter takes the values of 0.5, 1.5, and 2.5 and their collapse is induced initially by implementing an azimuthal mass perturbation. To assess the effects of the total mass of the parent gas structure on the nature of the resulting binary configurations, we also consider the collapse of two cores of M$_T$= 1 and 5 M$_{\odot}$. We calculate the collapse of all these parent gas structures using three values of the ratio of thermal energy to potential energy and for two values of the mass perturbation amplitude. We next calculate the binary separations, masses and integral properties of the binary fragments and present them in terms of the total mass of the parent structure. For most of our models, we finally calculate the $\beta$ extreme value, so that a model with a slightly higher $\beta$ value would no longer collapse.

By using the particle-based code Gadget2, we follow the evolution of a gas clump, in which a gravitational collapse is initially induced. The particles representing the gas clump have initially a velocity according to a turbulent spectrum built in a Fourier space of 64$^3$ grid elements. In a very early stage of evolution of the clump, a set of gas particles representing the wind, suddenly move outwards from the clump's center. We consider only two kinds of winds, namely: one with spherical symmetry and a second one being a bipolar collimated jet. In order to assess the dynamical change in the clump due to interaction with the winds, we show iso-velocity and iso-density plots for all our simulations.

Recent VLA and CARMA observations have shown proto-stars in binaries with unprecedented resolution. Specifically, the proto-stellar masses of systems such as CB230 IRS1 and L1165-SMM1 have been detected in the range of $0.1-0.25 \, M_{\odot}$. These are much more massive than the masses generally obtained by numerical simulations of binary formation, around $0.01 \, M_{\odot}$. Motivated by these discrepancies in mass, in this paper we study the formation mass of a binary system as a function of the total mass of its parent core. To achieve this objective, we present a set of numerical simulations of the gravitational collapse of a uniform and rotating core, in which azimuthal symmetric mass seeds are initially implemented in order to favor the formation of a dense filament, out of which a binary system may be formed by direct fragmentation. We first observed that this binary formation process is diminished when the total mass of the parent core $M_0$ is increased; then we increased the level of the ratio of kinetic energy to the gravitational energy, denoted by $\beta$, initially supplied to the rotating core, in order to achieve the desired direct fragmentation of the filament. Next, we found that the mass of the binary fragment increases with the mass of the parent core, as expected. In this paper we confirm this expectation and also measure the binary mass $M_f$ obtained from an initial $M_0$. We then show a schematic diagram $M_0$ vs $\beta$, where the desired binary configurations are located while the ratio of the thermal energy to the gravitational energy $\alpha$ is kept fixed. We also report some basic physical data of the proto-stellar fragments that form the binary system, including the formation mass $M_f$, and its corresponding $\alpha_f$ and $\beta_f$.

We test here the first stage of a route of modifications to be applied to the public GADGET2 code for dynamically identifying accretion centers during the collision process of two adjacent and identical gas cores. Each colliding core has a uniform density profile and rigid body rotation; its mass and size have been chosen to represent the observed core $L1544$; for the thermal and rotational energy ratios with respect to the potential energy, we assume the values $\alpha=0.3$ and $\beta=0.1$, respectively. These values favor the gravitational collapse of the core. We here study cases of both -head-on and off-center collisions, in which the pre-collision velocity increases the initial sound speed of the barotropic gas by up to several times. In a simulation the accretion centers are formed by the highest density particles, so we here report their location and properties in order to realize the collision effects on the collapsing and colliding cores. In one of the models we observe a roughly spherical distribution of accretion centers located at the front wave of the collision. In a forthcoming publication we will apply the full modified GADGET code to study the collision of turbulent cores.

In this paper we study the gravitational collapse of a molecular hydrogen gas cloud composed of a core plus a gas envelope surrounding the core. We numerically simulate the collapse of four cloud models to take a glimpse to the time evolution of several dynamic variables, such as the angular momentum and the $aem$ ratio, as well as the ratios between the thermal and rotational energies with respect to the potential gravitational energy, denoted as $\alpha$ and $\beta$, respectively, among others. We re-take those models introduced by Arreaga et.al (Astronomy and Astrophysics, {\bf Vol. 509}, (2010), pag. A96.) in the present paper in order to produce different outcomes of the collapsing cloud characterized in terms of the aforementioned dynamical variables. Such characterization was missing in the paper by Arreaga et.al (Astronomy and Astrophysics, {\bf Vol. 509}, (2010), pag. A96.), and here we show that the gas envelope extension effects on the collapsing core can be quantitatively compared.

We examine the problem of the collapse and fragmentation of molecular clouds with a Gaussian density distribution with high resolution, double precision numerical simulations using the GADGET-2 code. To describe the thermodynamic properties of the cloud during the collapse -to mimic the rise of temperature predicted by radiative transfer- we use a barotropic equation of state that introduces a critical density to separate the isothermal and adiabatic regimes. We discuss the effects of this critical density in the formation of multiple systems. We confirm the tendency found for Plummer and Gaussian models that if the collapse changes from isothermal to adiabatic at earlier times that occurs for the models with a lower critical density, the collapse is slowed down, and this enhances the fragments' change to survive. However, this effect happens up to a threshold density below which single systems tend to form. On the other hand, by setting a bigger initial perturbation amplitude, the collapse is faster and in some cases a final single object is formed.