The Galactic Branches as a Possible Evidence for Transient Spiral Arms

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.Comment: 17 pages, 9 figures. Accepted for publication in MNRA

Multi-scale accretion in dense cloud cores and the delayed formation of massive stars

The formation mechanism of massive stars remains one of the main open problems in astrophysics, in particular the relationship between the mass of the most massive stars, and that of the cores in which they form. Numerical simulations of the formation and evolution of large molecular clouds, within which dense cores and stars form self-consistently, show in general that the cores' masses increase in time, and also that the most massive stars tend to appear later (by a few to several Myr) than lower-mass stars. Here we present an idealized model that incorporates accretion onto the cores as well as onto the stars, in which the core's mass growth is regulated by a ``gravitational choking'' mechanism that does not involve any form of support. This process is of purely gravitational origin, and causes some of the mass accreted onto the core to stagnate there, rather than being transferred to the central stars. Thus, the simultaneous mass growth of the core and of the stellar mass can be computed. In addition, we estimate the mass of the most massive allowed star before its photoionizing radiation is capable of overcoming the accretion flow onto the core. This model constitutes a proof-of-concept for the simultaneous growth of the gas reservoir and the stellar mass, the delay in the formation of massive stars observed in cloud-scale numerical simulations, the need for massive, dense cores in order to form massive stars, and the observed correlation between the mass of the most massive star and the mass of the cluster it resides in. Also, our model implies that by the time massive stars begin to form in a core, a number of low-mass stars are expected to have already formed.Comment: Submitted to MNRAS. Originally submitted to Nature Astronomy, but withdrawn from that journal after not having received a reviewer's report for over four months. Comments welcom

Effects of Non-Circular Motions on Azimuthal Color Gradients

Assuming that density waves trigger star formation, and that young stars preserve the velocity components of the molecular gas where they are born, we analyze the effects that non-circular gas orbits have on color gradients across spiral arms. We try two approaches, one involving semi-analytical solutions for spiral shocks, and another with magnetohydrodynamic (MHD) numerical simulation data. We find that, if non-circular motions are ignored, the comparison between observed color gradients and stellar population synthesis models would in principle yield pattern speed values that are systematically too high for regions inside corotation, with the difference between the real and the measured pattern speeds increasing with decreasing radius. On the other hand, image processing and pixel averaging result in systematically lower measured spiral pattern speed values, regardless of the kinematics of stellar orbits. The net effect is that roughly the correct pattern speeds are recovered, although the trend of higher measured $\Omega_p$ at lower radii (as expected when non-circular motions exist but are neglected) should still be observed. We examine the Martinez-Garcia et al. (2009) photometric data and confirm that this is indeed the case. The comparison of the size of the systematic pattern speed offset in the data with the predictions of the semi-analytical and MHD models corroborates that spirals are more likely to end at Outer Lindblad Resonance, as these authors had already found.Comment: 32 pages, 15 figures, accepted to Ap

COMPETITIVE ACCRETION IN A SHEET GEOMETRY AND THE STELLAR IMF

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$. Comment: 20 pages, 12 figure

Errors in kinematic distances and our image of the Milky Way Galaxy

Errors in the kinematic distances, under the assumption of circular gas orbits, were estimated by performing synthetic observations of a model disk galaxy. It was found that the error is < 0.5 kpc for most of the disk when the measured rotation curve was used, but larger if the real rotation curve is applied. In both cases, the error is significantly larger at the positions of the spiral arms. The error structure is such that, when kinematic distances are used to develope a picture of the large scale density distribution, the most significant features of the numerical model are significantly distorted or absent, while spurious structure appears. By considering the full velocity field in the calculation of the kinematic distances, most of the original density structures can be recovered.Comment: Accepted for publication in A

Una solución vía el método de yamabe de un problema en el exponente critico de sobolev

En este artículo se probara que el siguiente problema semilineal elíptico con valores de frontera