No Arabic abstract
The formation of stars occurs in the dense molecular cloud phase of the interstellar medium. Observations and numerical simulations of molecular clouds have shown that supersonic magnetised turbulence plays a key role for the formation of stars. Simulations have also shown that a large fraction of the turbulent energy dissipates in shock waves. The three families of MHD shocks --- fast, intermediate and slow --- distinctly compress and heat up the molecular gas, and so provide an important probe of the physical conditions within a turbulent cloud. Here we introduce the publicly available algorithm, SHOCKFIND, to extract and characterise the mixture of shock families in MHD turbulence. The algorithm is applied to a 3-dimensional simulation of a magnetised turbulent molecular cloud, and we find that both fast and slow MHD shocks are present in the simulation. We give the first prediction of the mixture of turbulence-driven MHD shock families in this molecular cloud, and present their distinct distributions of sonic and Alfvenic Mach numbers. Using subgrid one-dimensional models of MHD shocks we estimate that ~0.03 % of the volume of a typical molecular cloud in the Milky Way will be shock heated above 50 K, at any time during the lifetime of the cloud. We discuss the impact of this shock heating on the dynamical evolution of molecular clouds.
We have performed Smoothed Particle Magneto-Hydrodynamics (SPMHD) calculations of colliding clouds to investigate the formation of massive stellar clusters, adopting a timestep criterion to prevent large divergence errors. We find that magnetic fields do not impede the formation of young massive clusters (YMCs), and the development of high star formation rates, although we do see a strong dependence of our results on the direction of the magnetic field. If the field is initially perpendicular to the collision, and sufficiently strong, we find that star formation is delayed, and the morphology of the resulting clusters is significantly altered. We relate this to the large amplification of the field with this initial orientation. We also see that filaments formed with this configuration are less dense. When the field is parallel to the collision, there is much less amplification of the field, dense filaments form, and the formation of clusters is similar to the purely hydrodynamical case. Our simulations reproduce the observed tendency for magnetic fields to be aligned perpendicularly to dense filaments, and parallel to low density filaments. Overall our results are in broad agreement with past work in this area using grid codes.
We present 3D hydrodynamic simulations of the adiabatic interaction of a shock with a dense, spherical cloud. We compare how the nature of the interaction changes with the Mach number of the shock, $M$, and the density contrast of the cloud, $chi$. We examine the differences with 2D axisymmetric calculations, perform detailed resolution tests, and compare inviscid results to those obtained with the inclusion of a $k$-$epsilon$ subgrid turbulence model. Resolutions of 32-64 cells per cloud radius are the minimum necessary to capture the dominant dynamical processes in 3D simulations, while the 3D inviscid and $k$-$epsilon$ simulations typically show very good agreement. Clouds accelerate and mix up to 5 times faster when they are poorly resolved. The interaction proceeds very similarly in 2D and 3D - although non-azimuthal modes lead to different behaviour, there is very little effect on key global quantities such as the lifetime of the cloud and its acceleration. In particular, we do not find significant differences in the hollowing or voiding of the cloud between 2D and 3D simulations with $M=10$ and $chi=10$, which contradicts previous work in the literature.
The turbulent destruction of a cloud subject to the passage of an adiabatic shock is studied. We find large discrepancies between the lifetime of the cloud and the analytical result of Hartquist et al. (1986). These differences appear to be due to the assumption in Hartquist et al. that mass-loss occurs largely as a result of lower pressure regions on the surface of the cloud away from the stagnation point, whereas in reality Kelvin-Helmholtz (KH) instabilities play a dominant role in the cloud destruction. We find that the true lifetime of the cloud (defined as when all of the material from the core of the cloud is well mixed with the intercloud material in the hydrodynamic cells) is about 6 times t_KHD, where t_KHD is the growth timescale for the most disruptive, long-wavelength, KH instabilities. These findings have wide implications for diffuse sources where there is transfer of material between hot and cool phases. The properties of the interaction as a function of Mach number and cloud density contrast are also studied. The interaction is milder at lower Mach numbers with the most marked differences occuring at low shock Mach numbers when the postshock gas is subsonic with respect to the cloud (i.e. M < 2.76). Material stripped off the cloud only forms a long tail-like feature if the density contrast of the cloud to the ambient medium, chi > 1e3.
We present an analytical model to identify thin discs in galaxies, and apply this model to a sample of SDSS MaNGA galaxies. This model fits the velocity and velocity dispersion fields of galaxies with regular kinematics. By introducing two parameters $zeta$ related to the comparison of the models asymmetric drift correction to the observed gas kinematics and $eta$ related to the dominant component of a galaxy, we classify the galaxies in the sample as disc-dominated, non-disc-dominated, or disc-free indicating galaxies with a dominating thin disc, a non-dominating thin disc, or no thin disc detection with our method, respectively. The dynamical mass resulting from our model correlates with stellar mass, and we investigate discrepancies by including gas mass and variation of the initial mass function. As expected, most spiral galaxies in the sample are disc-dominated, while ellipticals are predominantly disc-free. Lenticular galaxies show a dichotomy in their kinematic classification, which is related to their different star formation rates and gas fractions. We propose two possible scenarios to explain these results. In the first scenario, disc-free lenticulars formed in more violent processes than disc-dominated ones, while in the second scenario, the quenching processes in lenticulars lead to a change in their kinematic structures as disc-dominated lenticulars evolve to disc-free ones.
Winds from massive stars have velocities of 1000 km/s or more, and produce hot, high pressure gas when they shock. We develop a theory for the evolution of bubbles driven by the collective winds from star clusters early in their lifetimes, which involves interaction with the turbulent, dense interstellar medium of the surrounding natal molecular cloud. A key feature is the fractal nature of the hot bubbles surface. The large area of this interface with surrounding denser gas strongly enhances energy losses from the hot interior, enabled by turbulent mixing and subsequent cooling at temperatures T = 10^4-10^5 K where radiation is maximally efficient. Due to the extreme cooling, the bubble radius scales differently (R ~ t^1/2) from the classical Weaver77 solution, and has expansion velocity and momentum lower by factors of 10-10^2 at given R, with pressure lower by factors of 10^2 - 10^3. Our theory explains the weak X-ray emission and low shell expansion velocities of observed sources. We discuss further implications of our theory for observations of the hot bubbles and cooled expanding shells created by stellar winds, and for predictions of feedback-regulated star formation in a range of environments. In a companion paper, we validate our theory with a suite of hydrodynamic simulations.