We present a unified description of the scenario of Global Hierarchical Collapse and fragmentation (GHC) in molecular clouds (MCs), owing to the continuous decrease of the average Jeans mass in the contracting cloud. GHC constitutes a regime of collapses within collapses, in which small-scale collapses begin at later times, but occur on shorter timescales than large-scale ones. The difference in timescales allows for most of the clouds' mass to be dispersed by feedback from the first massive stars, maintaining the global star formation rate low. All scales accrete from their parent structures. The main features of GHC are: star-forming MCs are in an essentially pressureless regime, which produces filaments that accrete onto clumps and cores ("hubs"). The filaments constitute the collapse flow from cloud to hub scales and may approach a quasi-stationary state; the molecular and dense mass fractions of the clouds increase over time; the first (low-mass) stars appear several Myr after global contraction began; more massive stars appear after a few Myr in massive hubs resulting from the collapse of larger scales; the minimum fragment mass may extend well into the brown-dwarf regime; Bondi-Hoyle-Lyttleton accretion occurs at the protostellar and core scales, accounting for a near-Salpeter IMF; the extreme anisotropy of the filamentary network explains the difficulty in detecting large-scale infall signatures; the balance between inertial and gravitationally-driven motions in clumps evolves during the contraction; prestellar cores adopt Bonnor-Ebert-like profiles, but are contracting ever since early times when they may appear to be unbound and to require pressure confinement; stellar clusters develop radial age and mass segregation gradients. Finally, we discuss the incompatibility between supersonic turbulence and the observed scalings in the molecular hierarchy.

Considering gravitational waves propagating on the most general 4+N-dimensional space-time, we investigate the effects due to the N extra dimensions on the four-dimensional waves. All wave equations are derived in general and discussed. On Minkowski4 times an arbitrary Ricci-flat compact manifold, we find: a massless wave with an additional polarization, the breathing mode, and extra waves with high frequencies fixed by Kaluza-Klein masses. We discuss whether these two effects could be observed.

The removal of the interstellar medium (ISM) of disk galaxies through ram pressure stripping (RPS) has been extensively studied in numerous simulations. Nevertheless, the role of magnetic fields (MF) on the gas dynamics in this process has been hardly studied, although the MF influence on the large-scale disk structure is well established. With this in mind, we present a 3D magnetohydrodynamic (MHD) simulation of face-on RPS of a disk galaxy to study the impact of the galactic MF in the gas stripping. The main effect of including a galactic MF is a flared disk. When the intracluster medium (ICM) wind hits this flared disk, oblique shocks are produced at the interaction interface, where the ISM is compressed, generating a gas inflow from large radii towards the central regions of the galaxy. This inflow is observed for $\sim 150$ Myr and may supply the central parts of the galaxy with material for star formation while the outskirts of the disk are being stripped of gas, thus the oblique shocks can induce and enhance the star formation in the remaining disk. We also observed that the MF alters the shape and structure of the swept gas, giving a smooth appearance in the magnetized case and clumpier and filamentary-like morphology in the hydro case. Finally, we estimated the truncation radius expected for our models using the Gunn-Gott criterion and found that is in agreement with the simulations.

The de Rham-Gabadadze-Tolley massive gravity admits pp-wave backgrounds on which linear fluctuations are shown to undergo time advances for all values of the parameters. The perturbations may propagate in closed time-like curves unless the parameter space is constrained to a line. These classical phenomena take place well within the theory's validity regime.

We analyze the physical properties and energy balance of density enhancements in two SPH simulations of the formation, evolution, and collapse of giant molecular clouds. In the simulations, no feedback is included, so all motions are due either to the initial, decaying turbulence, or to gravitational contraction. We define clumps as connected regions above a series of density thresholds. The resultingfull set of clumps follows the generalized energy-equipartition relation $\sigma_{v}/R^{1/2} \propto \Sigma^{1/2}$, where $\sigma_{v}$ is the velocity dispersion, $R$ is the "radius", and $\Sigma$ is the column density. We interpret this as a natural consequence of gravitational contraction at all scales, rather than virial equilibrium. Nevertheless, clumps with low $\Sigma$ tend to show a large scatter around equipartition. In more than half of the cases, this scatter is dominated by external turbulent compressions that {\it assemble} the clumps, rather than by small-scale random motions that would disperse them. The other half does actually disperse. Moreover, clump sub-samples selected by means of different criteria exhibit different scalings. Sub-samples with narrow $\Sigma$ ranges follow Larson-like relations, although characterized by their respective value of $\Sigma$. Finally, we find that: i) clumps lying in filaments tend to appear sub-virial; ii) high-density cores ($n \ge 10^5$ cm$^3$) that exhibit moderate kinetic energy excesses often contain sink ("stellar") particles, and the excess disappears when the stellar mass is taken into account in the energy balance; iii) cores with kinetic energy excess but no stellar particles are truly in a state of dispersal.

With the use of a background Milky-Way-like potential model, we performed stellar orbital and magnetohydrodynamic (MHD) simulations. As a first experiment, we studied the gaseous response to a bisymmetric spiral arm potential: the widely employed cosine potential model and a self-gravitating tridimensional density distribution based model called PERLAS. Important differences are noticeable in these simulations, while the simplified cosine potential produces two spiral arms for all cases, the more realistic density based model produces a response of four spiral arms on the gaseous disk, except for weak arms -i.e. close to the linear regime- where a two-armed structure is formed. In order to compare the stellar and gas response to the spiral arms, we have also included a detailed periodic orbit study and explored different structural parameters within observational uncertainties. The four armed response has been explained as the result of ultra harmonic resonances, or as shocks with the massive bisymmetric spiral structure, among other. From the results of this work, and comparing the stellar and gaseous responses, we tracked down an alternative explanation to the formation of branches, based only on the orbital response to a self-gravitating spiral arms model. The presence of features such as branches, might be an indication of transiency of the arms.

Motivated by recent numerical simulations of molecular cloud (MC)evolution, in which the clouds engage in global gravitational contraction, and local collapse events culminate significantly earlier than the global collapse, we investigate the growth of density perturbations embedded in a collapsing background, to which we refer as an Inverse Hubble Flow (IHF). We use the standard procedure for the growth of perturbations in a universe that first expands (the usual Hubble Flow) and then recollapses (the IHF). We find that linear density perturbations immersed in an IHF grow faster than perturbations evolving in a static background (the standard Jeans analysis). A fundamental distinction between the two regimes is that, in the Jeans case, the time $\tau_\mathrm{nl}$ for a density fluctuation to become nonlinear increases without limit as its initial value approaches zero, while in the IHF case $\tau_\mathrm{nl} \le \tau_\mathrm{ff}$ always, where $\tau_\mathrm{ff}$ is the free-fall time of the background density. We suggest that this effect, although moderate, implies that small-scale density fluctuations embedded in globally-collapsing clouds must collapse earlier than their parent cloud, regardless of whether the initial amplitude of the fluctuations is moderate or strongly nonlinear, thus allowing the classical mechanism of Hoyle fragmentation to operate in multi-Jeans-mass MCs. More fundamentally, our results show that, contrary to the standard paradigm that fluctuations of all scales grow at the same rate in the linear regime, the hierarchical nesting of the fluctuations of different scales does affect their growth even in the linear stage.

In this work we derive the full 3-D kinematics of the near-infrared outflow HH 223, located in the dark cloud Lynds 723 (L723), where a well-defined quadrupolar CO outflow is found. HH 223 appears projected onto the two lobes of the east-west CO outflow. The radio continuum source VLA 2, towards the centre of the CO outflow, harbours a multiple system of low-mass young stellar objects. One of the components has been proposed to be the exciting source of the east-west CO outflow. From the analisys of the kinematics, we get further evidence on the relationship between the near-infrared and CO outflows and on the location of their exciting source. The proper motions were derived using multi-epoch, narrow-band H$_2$ (2.122 $\mu$m line) images. Radial velocities were derived from the 2.122 $\mu$m line of the spectra. Because of the extended (~5 arcmin), S-shaped morphology of the target, the spectra were obtained with the Multi-Object-Spectroscopy (MOS) observing mode using the instrument LIRIS at the 4.2m William Herschel Telescope. To our knowledge, this work is the first time that MOS observing mode has been successfully used in the near infrared range for an extended target.

We analyze the spiral structure that results in a numerical simulation of a galactic disk with stellar and gaseous components evolving in a potential that includes an axisymmetric halo and bulge. We perform a second simulation without the gas component to observe how it affects the spiral structure in the disk. To quantify this, we use a Fourier analysis and obtain values for the pitch angle and the velocity of the self-excited spiral pattern of the disk. The results show a tighter spiral in the simulation with gaseous component. The spiral structure is consistent with a superposition of waves, each with a constant pattern velocity in given radial ranges.

It is generally assumed that gas in a galactic disk follows closely non self-intersecting periodic stellar orbits. In order to test this common assumption, we have performed MHD simulations of a galactic-like disk under the influence of a spiral galactic potential. We also have calculated the actual orbit of a gas parcel and compared it to stable periodic stellar orbits in the same galactic potential and position. We found that the gaseous orbits approach periodic stellar orbits far from the major orbital resonances only. Gas orbits initialized at a given galactocentric distance but at different azimuths can be different, and scattering is conspicuous at certain galactocentric radii. Also, in contrast to the stellar behaviour, near the 4:1 (or higher order) resonance the gas follows nearly circular orbits, with much shorter radial excursions than the stars. Also, since the gas does not settle into a steady state, the gaseous orbits do not necessarily close on themselves.

We investigate the collapse of non-spherical substructures, such as sheets and filaments, which are ubiquitous in molecular clouds. Such non-spherical substructures collapse homologously in their interiors but are influenced by an edge effect that causes their edges to be preferentially accelerated. We analytically compute the homologous collapse timescales of the interiors of uniform-density, self-gravitating filaments and find that the homologous collapse timescale scales linearly with the aspect ratio. The characteristic timescale for an edge driven collapse mode in a filament, however, is shown to have a square root dependence on the aspect ratio. For both filaments and circular sheets, we find that selective edge acceleration becomes more important with increasing aspect ratio. In general, we find that lower dimensional objects and objects with larger aspect ratios have longer collapse timescales. We show that estimates for star formation rates, based upon gas densities, can be overestimated by an order of magnitude if the geometry of a cloud is not taken into account.

We present new [SII] images of the HH 30 jet and counterjet observed in 2006, 2007, and 2010 that allowed us to measure with improved accuracy the positions and proper motions of the jet and counterjet knots. Our results show that the motion of the knots is essentially ballistic, with the exception of the farthest knots, which trace the large scale C-shape bending of the jet. The observed bending of the jet can be produced by a relative motion of the HH 30 star with respect to its surrounding environment, caused either by a possible proper motion of the HH 30 star, or by the entrainment of environment gas by the red lobe of the nearby L1551-IRS 5 outflow. Alternatively, the bending can be produced by the stellar wind from a nearby CTTS, identified in the 2MASS catalog. The proper motion velocities of the knots of the counterjet show more variations than those of the jet. In particular, we identify two knots of the counterjet that have the same kinematic age but whose velocities differ by almost a factor of two. Thus, it appears that counterjet knots launched simultaneously can be ejected with very different velocities. We confirm that the observed wiggling of the jet and counterjet arises from the orbital motion of the jet source in a binary system. Precession is of secondary importance in shaping the jet. We derive an orbital period $\tau_o=114\pm2$ yr and a mass function $m\mu_c^3=0.014\pm0.006$ $M_\odot$. For a mass of the system of $m=0.45\pm0.04$ $M_\odot$ (the value inferred from the disk kinematics) we obtain a mass $m_j=0.31\pm0.04$ $M_\odot$ for the jet source, a mass $m_c=0.14\pm0.03$ $M_\odot$ for the companion, and a binary separation of $a=18.0\pm0.6$ AU. This binary separation coincides with the value required to account for the size of the inner hole observed in the disk, attributed to tidal truncation in a binary system.

HH 223 is the optical counterpart of a larger scale H2 outflow, driven by the protostellar source VLA 2A, in L723. Its poorly collimated and rather chaotic morphology suggested the Integral Field Spectroscopy (IFS) as an appropriate option to map the emission for deriving the physical conditions and the kinematics. Here we present new results based on the IFS observations made with the INTEGRAL system at the WHT. The brightest knots of HH 223 (\sim16 arcsec, 0.02 pc at a distance of 300 pc) were mapped with a single pointing in the spectral range 6200-7700 A. We obtained the emission-line intensity maps for Halpha, [NII] 6584 A and [SII] 6716, 6731 A, and explored the distribution of the excitation and electron density from [NII]/Halpha, [SII]/Halpha, and [SII] 6716/6731 line-ratio maps. Maps of the radial velocity field were obtained. We analysed the 3D-kinematics by combining the knot radial velocities, derived from IFS data, with the knot proper motions derived from multi-epoch, narrow-band images. The intensity maps built from IFS data reproduced well the morphology found in the narrow-band images. We checked the results obtained from previous long-slit observations with those derived from IFS spectra extracted with a similar spatial sampling. At the positions intersected by the slit, the physical conditions and kinematics derived from IFS are compatible with those derived from long-slit data. In contrast, significant discrepancies were found when the results from long-slit data were compared with the ones derived from IFS spectra extracted at positions shifted a few arcsec from those intersected by the slit. This clearly revealed IFS observations as the best choice to get a reliable picture of the HH emission properties.

Molecular clouds often exhibit filamentary or sheet-like shapes. We compute the free-fall time ($\tff$) for finite, uniform, self-gravitating circular sheets and filamentary clouds of small but finite thickness, so that their volume density $\rho$ can still be defined. We find that, for thin sheets, the free-fall time is larger than that of a uniform sphere with the same volume density by a factor proportional to $\sqrt{A}$, where the aspect ratio $A$ is given by $A=R/h$, $R$ being the sheet's radius and $h$ is its thickness. For filamentary clouds, the aspect ratio is defined as $A=L/\calR$, where $L$ is the filament's half length and $\calR$ is its (small) radius, and the modification factor is a more complicated, although in the limit of large $A$ it again reduces to nearly $\sqrt{A}$. We propose that our result for filamentary shapes naturally explains the ubiquitous configuration of clumps fed by filaments observed in the densest structures of molecular clouds. Also, the longer free-fall times for non-spherical geometries in general may contribute towards partially alleviating the "star-formation conundrum", namely, that the star formation rate in the Galaxy appears to be proceeding in a timescale much larger than the total molecular mass in the Galaxy divided by its typical free-fall time. If molecular clouds are in general formed by thin sheets and long filaments, then their relevant free-fall time may have been systematically underestimated, possibly by factors of up to one order of magnitude.

We report a set of numerical experiments aimed at addressing the applicability of competitive accretion to explain the high-mass end of the stellar initial mass function in a sheet geometry with shallow gravitational potential, in contrast to most previous simulations which have assumed formation in a cluster gravitational potential. Our flat cloud geometry is motivated by models of molecular cloud formation due to large-scale flows in the interstellar medium. The experiments consisted of SPH simulations of gas accretion onto sink particles formed rapidly from Jeans-unstable dense clumps placed randomly in the finite sheet. These simplifications allow us to study accretion with a minimum of free parameters, and to develop better statistics on the resulting mass spectra. We considered both clumps of equal mass and gaussian distributions of masses, and either uniform or spatially-varying gas densities. In all cases, the sink mass function develops a power law tail at high masses, with $dN/dlog M \propto M^{-\Gamma}$. The accretion rates of individual sinks follow $\dot{M} \propto M^2$ at high masses; this results in a continual flattening of the slope of the mass function towards an asymptotic form $\Gamma \sim 1$ (where the Salpeter slope is $\Gamma = 1.35$). The asymptotic limit is most rapidly reached when starting from a relatively broad distribution of initial sink masses. In general the resulting upper mass slope is correlated with the maximum sink mass; higher sink masses are found in simulations with flatter upper mass slopes. Although these simulations are of a highly idealized situation, the results suggest that competitive accretion may be relevant in a wider variety of environments than previously considered, and in particular that the upper mass distribution may generally evolve towards a limiting value of $\Gamma \sim 1$.

It is known that most of the craters on the surface of the Moon were created by the collision of minor bodies of the Solar System. Main Belt Asteroids, which can approach the terrestrial planets as a consequence of different types of resonance, are actually the main responsible for this phenomenon. Our aim is to investigate the impact distributions on the lunar surface that low-energy dynamics can provide. As a first approximation, we exploit the hyberbolic invariant manifolds associated with the central invariant manifold around the equilibrium point L_2 of the Earth - Moon system within the framework of the Circular Restricted Three - Body Problem. Taking transit trajectories at several energy levels, we look for orbits intersecting the surface of the Moon and we attempt to define a relationship between longitude and latitude of arrival and lunar craters density. Then, we add the gravitational effect of the Sun by considering the Bicircular Restricted Four - Body Problem. As further exploration, we assume an uniform density of impact on the lunar surface, looking for the regions in the Earth - Moon neighbourhood these colliding trajectories have to come from. It turns out that low-energy ejecta originated from high-energy impacts are also responsible of the phenomenon we are considering.

We present here new results of two-dimensional hydrodynamical simulations of the eruptive events of the 1840s (the great) and the 1890s (the minor) eruptions suffered by the massive star $\eta$ Car. The two bipolar nebulae commonly known as the Homunculus and the little Homunculus were formed from the interaction of these eruptive events with the underlying stellar wind. As in previous work (Gonzalez et al. 2004a, 2004b), we assume here an interacting, nonspherical multiple-phase wind scenario to explain the shape and the kinematics of both Homunculi, but adopt a more realistic parametrization of the phases of the wind. During the 1890s eruptive event, the outflow speed {\it decreased} for a short period of time. This fact suggests that the little Homunculus is formed when the eruption ends, from the impact of the post-outburst $\eta$ Car wind (that follows the 1890s event) with the eruptive flow (rather than by the collision of the eruptive flow with the pre-outburst wind, as claimed in previous models; Gonzalez et al. 2004a, 2004b). Our simulations reproduce quite well the shape and the observed expansion speed of the large Homunculus. The little Homunculus (which is embedded within the large Homunculus) becomes Rayleigh-Taylor unstable and develop filamentary structures that resembles the spatial features observed in the polar caps. In addition, we find that the interior cavity between the two Homunculi is partially filled by material that is expelled during the decades following the great eruption. This result may be connected with the observed double-shell structure in the polar lobes of the $\eta$ Car nebula. Finally, as in previous work, we find the formation of tenuous, equatorial, high-speed features that seem to be related to the observed equatorial skirt of $\eta$ Car.

(Abridged) The gravitational term for clouds and cores entering in the virial theorem is usually assumed to be equal to the gravitational energy, since the contribution to the gravitational force from the mass distribution outside the volume of integration is assumed to be negligible. Such approximation may not be valid in the presence of an important external net potential. In the present work we analyze the effect of an external gravitational field on the gravitational budget of a density structure. Our cases under analysis are (a) a giant molecular cloud (GMC) with different aspect ratios embedded within a galactic net potential, and (b) a molecular cloud core embedded within the gravitational potential of its parent molecular cloud. We find that for roundish GMCs, the tidal tearing due to the shear in the plane of the galaxy is compensated by the tidal compression in the z direction. The influence of the external effective potential on the total gravitational budget of these clouds is relatively small, although not necessarily negligible. However, for more filamentary GMCs, the external effective potential can be dominant and can even overwhelm self-gravity, regardless of whether its main effect on the cloud is to disrupt it or compress it. This may explain the presence of some GMCs with few or no signs of massive star formation, such as the Taurus or the Maddalena's clouds. In the case of dense cores embedded in their parent molecular cloud, we found that the gravitational content due to the external field may be more important than the gravitational energy of the cores themselves. This effect works in the same direction as the gravitational energy, i.e., favoring the collapse of cores. We speculate on the implications of these results for star formation models.

(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 study the stability of proto-planetary disks with vertical velocity gradients in their equilibrium rotation rates; such gradients are expected to develop when dust settles into the midplane. Using a linear stability analysis of a simple three-layer model, we show that the onset of instability occurs at a larger value of the Richardson number, and therefore for a thicker layer, when the effects of Coriolis forces are included. This analysis also shows that even-symmetry (midplane-crossing) modes develop faster than odd-symmetry ones. These conclusions are corroborated by a large number of nonlinear numerical simulations with two different parameterized prescriptions for the initial (continuous) dust distributions. Based on these numerical experiments, the Richardson number required for marginal stability is more than an order of magnitude larger than the traditional 1/4 value. The dominant modes that grow have horizontal wavelengths of several initial dust scale heights, and in nonlinear stages mix solids fairly homogeneously over a comparable vertical range. We conclude that gravitational instability may be more difficult to achieve than previously thought, and that the vertical distribution of matter within the dust layer is likely globally, rather than locally, determined.