Source author record

Gilberto C. Gómez

Gilberto C. Gómez appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR astro-ph.CO

Catalog footprint

What is connected

10works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A resolution criterion based on characteristic timescales for MHD simulations of molecular clouds

We investigate the effect of numerical magnetic diffusion in magnetohydrodynamic simulations of magnetically supported molecular clouds. To this end, we have performed numerical studies on adaptive mesh isothermal simulations of marginally sub-critical molecular clouds. We find that simulations with low and intermediate resolutions collapse, contrary to what is theoretically expected. However, the simulation with the highest numerical resolution oscillates around an equilibrium state without collapsing. In order to quantify the numerical diffusion of the magnetic field, we ran a second suit of current-sheet simulations in which the numerical magnetic diffusion coefficient can be directly measured, and computed the corresponding diffusion times at various numerical resolutions. On this basis, we propose a criterion for the resolution of magnetic fields in MHD simulations based on requiring that the diffusion time to be larger than the characteristic timescale of the physical process responsible for the dynamic evolution of the structure.

preprint2022arXiv

Different molecular filament widths as tracers of accretion onto filaments

We explore how dense filament widths, when measured using different molecular species, may change as a consequence of gas accretion toward the filament. As a gas parcel falls into the filament, it will experience different density, temperature, and extinction values. The rate at which this environment changes will affect differently the abundance of different molecules. So, a molecule that forms quickly will better reflect the local physical conditions a gas parcel experiences than a slower-forming molecule. Since these differences depend on how the respective timescales compare, the different molecular distributions should reflect how rapidly the environment changes, i.e., the accretion rate toward the filament. We find that the filament widths measured from time-dependent abundances for C2H, CO, CN, CS, and C3H2, are sensitive the most to this effect, being those molecules the ones presenting also the wider filament widths. On the contrary, molecules such as N2H+, NH3, H2CO, HNC or CH3OH are not so sensitive to accretion and present the narrowest filament widths. We propose that ratios of filament widths for different tracers could be a useful tool to estimate the accretion rate onto the filament.

preprint2021arXiv

Density profile evolution during prestellar core collapse: Collapse starts at the large scale

We study the gravitationally-dominated, accretion-driven evolution of a prestellar core. In our model, as the core's density increases, it remains immersed in a constant-density environment and so it accretes from this environment, increasing its mass and reducing its Jeans length. Assuming a power-law density profile $ρ\propto r^{-p}$, we compute the rate of change of the slope $p$, and show that the value $p=2$ is stationary, and furthermore, an attractor. The radial profile of the Jeans length scales as $r^{p/2}$, implying that, for $p<2$, there is a radius below which the region is smaller than its Jeans length, thus appearing gravitationally stable and in need of pressure confinement, while, in reality, it is part of a larger-scale collapse and is undergoing compression by the infalling material. In this region, the infall speed decreases towards the center, eventually becoming subsonic, thus appearing "coherent", without the need for turbulence dissipation. We present a compilation of observational determinations of density profiles in dense cores and show that the distribution of their slopes peaks at $p \sim 1.7$--1.9, supporting the notion that the profile steepens over time. Finally, we discuss the case of magnetic support in a core in which the field scales as $B \propto ρ^β$. For the expected value of $β= 2/3$, this implies that the mass to magnetic flux ratio also decreases towards the central parts of the cores, making them appear magnetically supported, while in reality they may be part of larger collapsing supercritical region. We conclude that local signatures of either thermal or magnetic support are not conclusive evidence of stability, that the gravitational instability of a region must be established at the large scales, and that the prestellar stage of collapse is dynamic rather than quasistatic.

preprint2016arXiv

Energy budget of forming clumps in numerical simulations of collapsing clouds

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 $σ_{v}/R^{1/2} \propto Σ^{1/2}$, where $σ_{v}$ is the velocity dispersion, $R$ is the "radius", and $Σ$ 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 $Σ$ 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 $Σ$ ranges follow Larson-like relations, although characterized by their respective value of $Σ$. 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.

preprint2015arXiv

Inverse Hubble Flows in Molecular Clouds

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 $τ_\mathrm{nl}$ for a density fluctuation to become nonlinear increases without limit as its initial value approaches zero, while in the IHF case $τ_\mathrm{nl} \le τ_\mathrm{ff}$ always, where $τ_\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.

preprint2015arXiv

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.

preprint2013arXiv

Comparing Gaseous and Stellar Orbits in a Spiral Potential

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.

preprint2012arXiv

Aspect Ratio Dependence of the Free-Fall Time for Non-Spherical Symmetries

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.

preprint2011arXiv

The Free-Fall time of finite Sheets and Filaments

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 $ρ$ 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.

preprint2010arXiv

Competitive Accretion in 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^{-Γ}$. 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 $Γ\sim 1$ (where the Salpeter slope is $Γ= 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 $Γ\sim 1$.