Monte carlo simulation for molecular gas dynamics

The dynamics of low-density flows is governed by the Boltzmann equation of the kinetic theory of gases. This is a nonlinear integro-differential equation and, in general, numerical methods must be used to obtain its solution. The present paper, after a brief review of Direct Simulation Monte Carlo (DSMC) methods due to Bird, and Belotserkovskii and Yanitskii, studies the details of theDSMC method of Deshpande for mono as well as multicomponent gases. The present method is a statistical particle-in-cell method and is based upon the Kac-Prigogine master equation which reduces to the Boltzmann equation under the hypothesis of molecular chaos. The proposed Markoff model simulating the collisions uses a Poisson distribution for the number of collisions allowed in cells into which the physical space is divided. The model is then extended to a binary mixture of gases and it is shown that it is necessary to perform the collisions in a certain sequence to obtain unbiased simulation.

Lurking in ULIRGs: Molecular Gas in local Merging Galaxies

Mergers and interacting galaxies are pivotal to the evolution of galaxies in the universe. They are the sites of prodigious star formation and key to understanding the starburst processes: the physical and chemical properties and the dynamics of the molecular gas. ULIRGs or Ultraluminous Infrared Galaxies are a result of many of these mergers. They host extreme starbursts, AGNs, and mergers. They are the perfect laboratory to probe the connection between starbursts, black hole accretion and mergers and to further our understanding of star formation and merging.

NGC 6240 and Arp 220 can be considered the founding members of this very active class of objects. They are in different stages of merging and hence are excellent case studies to further our understanding about the merging process. We have imaged the dense star-forming regions of these galaxies at sub-arcsec resolution with CARMA C and B Configurations (2" and 0.5 - 0.8"). Multi-band imaging allows excitation analysis of HCN, HCO⁺, HNC, and CS along with CO transitions to constrain the properties of the gas. Our dataset is unique in that we have observed these lines at similar resolutions and high sensitivity which can be used to derive line ratios of faint high excitation lines.

Arp 220 has not had confirmed X-ray AGN detections for either nuclei. However, our observations indicate HCN/HNC ratios consistent with the chemistry of X-ray Dominated Regions (XDRs) -- a likely symptom of AGN. We calculated the molecular Hydrogen densities using each of the molecular species and conclude that assuming abundances of HNC and HCO⁺ similar to those in galactic sources are incorrect in the case of ULIRGs. The physical conditions in the dense molecular gas in ULIRGs alter these abundances. The derived H2 volume densities are ~ 5 x 10⁴ cm^-3 in both Arp 220 nuclei and ~ 10⁴ cm^-3 in NGC 6240.

Molecular gas in the galaxy M83 : Its distribution, kinematics, and relation to star formation

The barred spiral galaxy M83 (NGC5236) has been observed in the 12CO J=1–0 and J=2–1 millimetre lines with the Swedish-ESO Submillimetre Telescope (SEST). The sizes of the CO maps are 100×100, and they cover the entire optical disk. The CO emission is strongly peaked toward the nucleus. The molecular spiral arms are clearly resolved and can be traced for about 360º. The total molecular gas mass is comparable to the total Hi mass, but H2 dominates in the optical disk. Iso-velocity maps show the signature of an inclined, rotating disk, but also the effects of streaming motions along the spiral arms. The dynamical mass is determined and compared to the gas mass. The pattern speed is determined from the residual velocity pattern, and the locations of various resonances are discussed. The molecular gas velocity dispersion is determined, and a trend of decreasing dispersion with increasing galactocentric radius is found. A total gas (H2+Hi+He) mass surface density map is presented, and compared to the critical density for star formation of an isothermal gaseous disk. The star formation rate (SFR) in the disk is estimated using data from various star formation tracers. The different SFR estimates agree well when corrections for extinctions, based on the total gas mass map, are made. The radial SFR distribution shows features that can be associated with kinematic resonances. We also find an increased star formation efficiency in the spiral arms. Different Schmidt laws are fitted to the data. The star formation properties of the nuclear region, based on high angular resolution HST data, are also discussed.

Tidally compressed gas in centers of early-type and ultraluminous galaxies

In this paper we propose that the compressive tidal held in the centers of flat-core early-type galaxies and ultraluminous galaxies compresses molecular clouds producing dense gas observed in the centers of these galaxies. The effect of galactic tidal fields is usually considered disruptive in the literature. However, for some galaxies, the mass profile flattens toward the center and the resulting galactic tidal field is not disruptive, but instead it is compressive within the flat-core region. We have used the virial theorem to determine the minimum density of a molecular cloud to be stable and gravitationally bound within the tidally compressive region of a galaxy. We have applied the mechanism to determine the mean molecular cloud densities in the centers of a sample of flat-core, early-type galaxies and ultraluminous galaxies. For early-type galaxies with a core-type luminosity profile, the tidal held of the galaxy is compressive within half the core radius. We have calculated the mean gas densities for molecular gas in a sample of early-type galaxies which have already been detected in CO emission, and we obtain mean densities of [n] similar to 10(3)-10(6) cm(-3) within the central 100 pc radius. We also use our model to calculate the molecular cloud densities in the inner few hundred parsecs of a sample of ultraluminous galaxies. From the observed rotation curves of these galaxies we show that they have a compressive core within their nuclear region. Our model predicts minimum molecular gas densities in the range 10(2)-10(4) cm(-3) in the nuclear gas disks; the smaller values are applicable typically for galaxies with larger core radii. The resulting density values agree well with the observed range. Also, for large core radii, even fairly low-density gas (similar to 10(2) cm(-3)) can remain bound and stable close to the galactic center.

Drag on an ellipsoid particle in free-molecular plasma flow

Based on the analysis of molecular gas dynamics, the drag and moment acting on an ellipsoid particle of revolution X-2/a(2) + Y-2/a(2) + Z(2)/c(2) = 1, as an example of nonspherical particles, are studied under the condition of free-molecular plasma flow with thin plasma sheaths. A nonzero moment which causes nonspherical particle self-oscillation and self-rotation around its own axis in the plasma flow-similar to the pitching moment in aerodynamics-is discovered for the first time. When the ratio of axis length c/a is unity, the moment is zero and the drag formula are reduced to the well-known results of spherical particles. The effects of the particle-plasma relative velocity, the plasma temperature, and the particle materials on the drag and moment are also investigated.

Rarefied gas dynamics: fundamentals, simulations and micro flows

A survey of molecular line emission towards ultracompact HII regions

We have used the JCMT to survey molecular line emission towards 14 ultracompact HII regions (G5.89, G9.62, G10.30, G10.47, G12.21, G13.87, G29.96, G31.41, G34.26, G43.89, G45.12, G45.45, G45.47, and G75.78). For each source, we observed up to ten 1 GHz bands between 200 and 350 GHz, covering lines of more than 30 species including multiple transitions of CO isotopes, CH3OH, CH3CCH, CH3CN and HCOOCH3, and sulphuretted molecules. The number of transitions detected varied by a factor of 20 between sources; which were chosen following observations of high-excitation ammonia (Cesaroni et al. 1994a) and methyl cyanide (Olmi et al. 1993). In half our sample (the line-poor sources), only (CO)-O-17: (CO)-O-18, SO, (CS)-S-34 and CH3OH were detected. In the line-rich sources, we detected over 150 lines, including high excitation lines of CH3CN, HCOOCH3; C2H5CN, CH3OH, and CH3CCH. We have calculated the physical conditions of the molecular gas. To reproduce the emission from the line-rich sources requires both a hot, dense compact core and an ambient cloud consisting of less dense, cooler gas. The hot cores, which are less than 0.1 pc in size; reach densities of at least 10(8) cm(-3) and temperatures of more than 80 K. The line-poor sources can be modelled without a hot core by a 20-30 K, 10(5) cm(-3) cloud. We find no correlation between the size of the HII region and the current physical conditions in the molecular environment. A comparison with chemical models (Millar et al. 1997) confirms that grain surface chemistry is important in hot cores.

Gas streaming motions towards the nucleus of M81

We present two-dimensional stellar and gaseous kinematics of the inner 120 x 250 pc2 of the LINER/Seyfert 1 galaxy M81, from optical spectra obtained with the Gemini Multi-Object Spectrograph (GMOS) integral field spectrograph on the Gemini-North telescope at a spatial resolution of approximate to 10 pc. The stellar velocity field shows circular rotation and, overall, is very similar to the published large-scale velocity field, but deviations are observed close to the minor axis which can be attributed to stellar motions possibly associated with a nuclear bar. The stellar velocity dispersion of the bulge is 162 +/- 15 km s-1, in good agreement with previous measurements and leading to a black hole mass of M(BH) = 5.5+3.6(-2.0) x 107 M(circle dot) based on the M(BH)-Sigma relationship. The gas kinematics is dominated by non-circular motions and the subtraction of the stellar velocity field reveals blueshifts of approximate to-100 km s-1 on the far side of the galaxy and a few redshifts on the near side. These characteristics can be interpreted in terms of streaming towards the centre if the gas is in the plane. On the basis of the observed velocities and geometry of the flow, we estimate a mass inflow rate in ionized gas of approximate to 4.0 x 10-3 M(circle dot) yr-1, which is of the order of the accretion rate necessary to power the LINER nucleus of M81. We have also applied the technique of principal component analysis (PCA) to our data, which reveals the presence of a rotating nuclear gas disc within approximate to 50 pc from the nucleus and a compact outflow, approximately perpendicular to the disc. The PCA combined with the observed gas velocity field shows that the nuclear disc is being fed by gas circulating in the galaxy plane. The presence of the outflow is supported by a compact jet seen in radio observations at a similar orientation, as well as by an enhancement of the [O i]/H alpha line ratio, probably resulting from shock excitation of the circumnuclear gas by the radio jet. With these observations we are thus resolving both the feeding - via the nuclear disc and observed gas inflow, and the feedback - via the outflow, around the low-luminosity active nucleus of M81.

Ammonia observations of the nearby molecular cloud MBM 12

We present NH_3(1,1) and (2,2) observations of MBM 12, the closest known molecular cloud (65-pc distance), aimed at finding evidence for on-going star formation processes, No local temperature (with a T_rot upper limit of 12 K) or linewidth enhancement is found, which suggests that the area of the cloud that we have mapped (∼ 15-arcmin size) is not currently forming stars. Therefore this nearby 'starless' molecular gas region is an ideal laboratory to study the physical conditions preceding new star formation. A radio continuum source has been found in Very Large Array archive data, close to but outside the NH_3 emission. This source is likely to be a background object.

Kinetic theory for sheared granular flows

Rapid granular flows are far-from-equilibrium-driven dissipative systems where the interaction between the particles dissipates energy, and so a continuous supply of energy is required to agitate the particles and facilitate the rearrangement required for the flow. This is in contrast to flows of molecular fluids, which are usually close to equilibrium, where the molecules are agitated by thermal fluctuations. Sheared granular flows form a class of flows where the energy required for agitating the particles in the flowing state is provided by the mean shear. These flows have been studied using the methods of kinetic theory of gases, where the particles are treated in a manner similar to molecules in a molecular gas, and the interactions between particles are treated as instantaneous energy-dissipating binary collisions. The validity of the assumptions underlying kinetic theory, and their applicability to the idealistic case of dilute sheared granular flows are first discussed. The successes and challenges for applying kinetic theory for realistic dense sheared granular flows are then summarised. (C) 2014 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Diagnostics on an atmospheric pressure plasma jet

The atmospheric pressure plasma jet (APPJ) is a homogeneous non-equilibrium discharge at ambient pressure. It operates with a noble base gas and a percentage-volume admixture of a molecular gas. Applications of the discharge are mainly based on reactive species in the effluent. The effluent region of a discharge operated in helium with an oxygen admixture has been investigated. The optical emission from atomic oxygen decreases with distance from the discharge but can still be observed several centimetres in the effluent. Ground state atomic oxygen, measured using absolutely calibrated two-photon laser induced fluorescence spectroscopy, shows a similar behaviour. Detailed understanding of energy transport mechanisms requires investigations of the discharge volume and the effluent region. An atmospheric pressure plasma jet has been designed providing excellent diagnostics access and a simple geometry ideally suited for modelling and simulation. Laser spectroscopy and optical emission spectroscopy can be applied in the discharge volume and the effluent region.

Sulphur chemistry and evolution in hot cores

We compare the results of our JCMT spectral line survey of molecular gas towards ultracompact HII regions with the predictions of models of sulphur chemistry in hot cores. We investigate the range of evolutionary models that are consistent with the observed physical conditions and chemical abundances, and see to what extent it is possible to constrain core ages by comparing abundances with the predictions of chemical models. The observed abundance ratios vary little from source to source, suggesting that all the sources are at a similar evolutionary stage. The models are capable of predicting the observed abundances of H2S, SO, SO2, and CS. The models fail to predict the amount of OCS observed, suggesting that an alternative formation route is required. An initial H2S abundance from grain mantle evaporation of similar to 10(-7) is preferred.

Large- and small-scale structure of the intermediate- and high-velocity clouds towards the LMC and SMC.

We employ Ca II K and Na I D interstellar absorption-line spectroscopy of early-type stars in the Large and Small Magellanic Clouds (LMC, SMC) to investigate the large- and small-scale structure in foreground intermediate- and high-velocity clouds (I/HVCs). Data include FLAMES-GIRAFFE Ca II K observations of 403 stars in four open clusters, plus FEROS or UVES spectra of 156 stars in the LMC and SMC. The FLAMES observations are amongst the most extensive probes to date of Ca II structures on ∼20 arcsec scales in Magellanic I/HVCs. From the FLAMES data within a 0 ∘.∘.∘.5 field of view, the Ca II K equivalent width in the I/HVC components towards three clusters varies by factors of ≥10. There are no detections of molecular gas in absorption at intermediate or high velocities, although molecular absorption is present at LMC and Galactic velocities towards some sightlines. The FEROS/UVES data show Ca II K I/HVC absorption in ∼60 per cent of sightlines. The range in the Ca II/Na I ratio in I/HVCs is from –0.45 to +1.5 dex, similar to previous measurements for I/HVCs. In 10 sightlines we find Ca II/O I ratios in I/HVC gas ranging from 0.2 to 1.5 dex below the solar value, indicating either dust or ionization effects. In nine sightlines I/HVC gas is detected in both H I and Ca II at similar velocities, implying that the two elements form part of the same structure.

Recherche d’étoiles jeunes de faible masse dans le voisinage solaire

Formées lors de l’effondrement gravitationnel d’un nuage de gaz moléculaire, les étoiles naissantes auront différentes masses variant entre 0.08 et environ 100M . La majorité de la population stellaire de la Galaxie est constituée d’étoiles dont la masse est inférieure à environ 0.6 M . Le dernier évènement de formation stellaire dans le voisinage solaire s’est produit dans la bulle locale il y a au plus 100 millions d’années, vraisemblablement provoqué par le passage d’une onde de choc dans le bras local de la Galaxie. C’est ainsi que se formèrent de jeunes associations d’étoiles dont les membres se caractérisent en particulier par une vitesse spatiale et une position commune dans la Galaxie. Les associations jeunes étant peu densément peuplées et relativement proches du Soleil, leurs membres se font plutôt rares et dispersés sur toute la voûte céleste. Jusqu’à présent, surtout les étoiles les plus massives (brillantes) ont été répertoriées. Les étoiles jeunes de faible masse, constituant la majorité de la population, restent pour la plupart à être identifiées. Les étoiles jeunes de faible masse représentent une population clef pour contraindre les modèles évolutifs des étoiles M et des naines brunes. Elles sont également d’excellentes candidates pour chercher des exoplanètes via les techniques d’imagerie directe. Ce mémoire présente une nouvelle méthode utilisant un modèle cinématique enrichi d’une analyse statistique Bayesienne pour identifier des étoiles jeunes de faible masse dans les associations beta Pictoris, Tucana-Horologium et AB Doradus. À partir d’un échantillon de 1080 étoiles K et M, toutes comportant des indicateurs de jeunesse tels l’émission Halpha et une forte luminosité dans les rayons X, leurs propriétés cinématiques (mouvement propre) et photométriques sont analysées pour en extraire 98 candidates hautement probables membres d’une des trois associations. Une confirmation de leur statut comme membre nécessitera en particulier une mesure de leur vitesse radiale (prédit par notre analyse) et une mesure de la largeur équivalente du lithium à 6708 Å pour mieux contraindre leur âge.

Étude de variabilité photométrique infrarouge de naines brunes

Les naines brunes sont des objets astronomiques de faible masse ( 0.012 - 0.075 M_Sun ) et de basse température ( T < 3,500 K ). Bien qu’elles se forment comme des étoiles, c’est-à-dire par l’effondrement d’un nuage de gaz moléculaire, les naines brunes n’ont pas une masse suffisante pour entretenir des réactions de fusion nucléaire en leur coeur. Les naines brunes relativement chaudes (type L) sont recouvertes de nuages de poussière mais ces derniers disparaissent progressivement de l’atmosphère lorsque la température chute sous les 1,500 K (type T). Les naines brunes près de la transition L/T devraient donc être partiellement recouvertes de nuages. De par leur rotation relativement rapide (2 h - 12 h), le couvert nuageux inhomogène des naines brunes devrait produire une variabilité photométrique observable en bande J (1.2 um), la longueur d’onde à laquelle les nuages ont la plus forte opacité. Ce mémoire présente les résultats d’une recherche de variabilité photométrique infrarouge pour une dizaine de naines brunes de type spectral près de la transition L/T. Les observations, obtenues à l’Observatoire du Mont-Mégantic, ont permis le suivi photométrique en bande J de neuf cibles. Une seule d’entre elles, SDSS J105213.51+442255.7 (T0.5), montre des variations périodiques sur une période d’environ 3 heures avec une amplitude pic-à-pic variant entre 40 et 80 mmag. Pour les huit autres cibles, on peut imposer des limites (3 sigma) de variabilité périodique à moins de 15 mmag pour des périodes entre 1 et 6 heures. Ces résultats supportent l’hypothèse qu’un couvert nuageux partiel existe pour des naines brunes près de la transition L/T mais ce phénomène demeure relativement peu fréquent.