Search NASA⌕ Search

Engineering topics

Glover, Simon C. O.

Publications and source records attributed to Glover, Simon C. O..

Applying the Tremaine–Weinberg Method to Nearby Galaxies: Stellar-mass-based Pattern Speeds and Comparisons with ISM Kinematics

We apply the Tremaine–Weinberg method to 19 nearby galaxies using stellar mass surface densities and velocities derived from the PHANGS-MUSE survey, to calculate (primarily bar) pattern speeds (Ω{sub P}). After quality checks, we find that around half (10) of these stellar-mass-based measurements are reliable. For those galaxies, we find good agreement between our results and previously published pattern speeds, and we use rotation curves to calculate major resonance locations (corotation radii and Lindblad resonances). We also compare these stellar-mass-derived pattern speeds with Hα (from MUSE) and CO(J = 2 − 1) emission from the PHANGS-ALMA survey. We find that in the case of these clumpy interstellar medium (ISM) tracers, this method erroneously gives a signal that is simply the angular frequency at a representative radius set by the distribution of these clumps (Ω{sub clump}), and that this Ω{sub clump} is significantly different from Ω{sub P} (∼20% in the case of Hα, and ∼50% in the case of CO). Thus, we conclude that it is inadvisable to use “pattern speeds” derived from ISM kinematics. Finally, we compare our derived pattern speeds and corotation radii, along with bar properties, to the global parameters of these galaxies. Consistent with previous studies, we find that galaxies with a later Hubble type have a larger ratio of corotation radius to bar length, more molecular-gas-rich galaxies have higher Ω{sub P}, and more bulge-dominated galaxies have lower Ω{sub P}. Unlike earlier works, however, there are no clear trends between the bar strength and Ω{sub P}, nor between the total stellar mass surface density and the pattern speed.

74 ATOMIC AND MOLECULAR PHYSICS↗

How the First Stars Regulated Star Formation. II. Enrichment by Nearby Supernovae

Metals from Population III (Pop III) supernovae led to the formation of less massive Pop II stars in the early universe, altering the course of evolution of primeval galaxies and cosmological reionization. There are a variety of scenarios in which heavy elements from the first supernovae were taken up into second-generation stars, but cosmological simulations only model them on the largest scales. We present small-scale, high-resolution simulations of the chemical enrichment of a primordial halo by a nearby supernova after partial evaporation by the progenitor star. We find that ejecta from the explosion crash into and mix violently with ablative flows driven off the halo by the star, creating dense, enriched clumps capable of collapsing into Pop II stars. Metals may mix less efficiently with the partially exposed core of the halo, so it might form either Pop III or Pop II stars. Both Pop II and III stars may thus form after the collision if the ejecta do not strip all the gas from the halo. The partial evaporation of the halo prior to the explosion is crucial to its later enrichment by the supernova.

79 ASTRONOMY AND ASTROPHYSICS↗

THE ROLE OF COSMIC-RAY PRESSURE IN ACCELERATING GALACTIC OUTFLOWS

We study the formation of galactic outflows from supernova (SN) explosions with the moving-mesh code AREPO in a stratified column of gas with a surface density similar to the Milky Way disk at the solar circle. We compare different simulation models for SN placement and energy feedback, including cosmic rays (CRs), and find that models that place SNe in dense gas and account for CR diffusion are able to drive outflows with similar mass loading as obtained from a random placement of SNe with no CRs. Despite this similarity, CR-driven outflows differ in several other key properties including their overall clumpiness and velocity. Moreover, the forces driving these outflows originate in different sources of pressure, with the CR diffusion model relying on non-thermal pressure gradients to create an outflow driven by internal pressure and the random-placement model depending on kinetic pressure gradients to propel a ballistic outflow. CRs therefore appear to be non-negligible physics in the formation of outflows from the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗

IMPACT OF SUPERNOVA AND COSMIC-RAY DRIVING ON THE SURFACE BRIGHTNESS OF THE GALACTIC HALO IN SOFT X-RAYS

The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order 10{sup −12} erg cm{sup −2} s{sup −1} deg{sup −2}. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-phase interstellar medium including thermal feedback by supernova explosions as well as cosmic-ray feedback. We include a time-dependent chemical network, self-shielding by gas and dust, and self-gravity. Pure thermal feedback alone is sufficient to produce the observed surface brightness, although it is very sensitive to the supernova rate. Cosmic rays suppress this sensitivity and reduce the surface brightness because they drive cooler outflows. Self-gravity has by far the largest effect because it accumulates the diffuse gas in the disk in dense clumps and filaments, so that supernovae exploding in voids can eject a large amount of hot gas into the halo. This can boost the surface brightness by several orders of magnitude. Although our simulations do not reach a steady state, all simulations produce surface brightness values of the same order of magnitude as the observations, with the exact value depending sensitively on the simulation parameters. We conclude that star formation feedback alone is sufficient to explain the origin of the hot halo gas, but measurements of the surface brightness alone do not provide useful diagnostics for the study of galactic star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Structure analysis of simulated molecular clouds with the Δ-variance

Here, we employ the Δ-variance analysis and study the turbulent gas dynamics of simulated molecular clouds (MCs). Our models account for a simplified treatment of time-dependent chemistry and the non-isothermal nature of the gas. We investigate simulations using three different initial mean number densities of n 0 = 30, 100 and 300 cm -3 that span the range of values typical for MCs in the solar neighbourhood. Furthermore, we model the CO line emission in a post-processing step using a radiative transfer code. We evaluate Δ-variance spectra for centroid velocity (CV) maps as well as for integrated intensity and column density maps for various chemical components: the total, H 2 and 12 CO number density and the integrated intensity of both the 12 CO and 13 CO (J = 1 → 0) lines. The spectral slopes of the Δ-variance computed on the CV maps for the total and H 2 number density are significantly steeper compared to the different CO tracers. We find slopes for the linewidth–size relation ranging from 0.4 to 0.7 for the total and H 2 density models, while the slopes for the various CO tracers range from 0.2 to 0.4 and underestimate the values for the total and H 2 density by a factor of 1.5–3.0. We demonstrate that optical depth effects can significantly alter the Δ-variance spectra. Furthermore, we report a critical density threshold of 100 cm -3 at which the Δ-variance slopes of the various CO tracers change sign. We thus conclude that carbon monoxide traces the total cloud structure well only if the average cloud density lies above this limit.

79 ASTRONOMY AND ASTROPHYSICS↗

A NEW APPROACH TO DETERMINE OPTICALLY THICK H{sub 2} COOLING AND ITS EFFECT ON PRIMORDIAL STAR FORMATION

We present a new method for estimating the H{sub 2} cooling rate in the optically thick regime in simulations of primordial star formation. Our new approach is based on the TreeCol algorithm, which projects matter distributions onto a spherical grid to create maps of column densities for each fluid element in the computational domain. We have improved this algorithm by using the relative gas velocities to weight the individual matter contributions with the relative spectral line overlaps, in order to properly account for the Doppler effect. We compare our new method to the widely used Sobolev approximation, which yields an estimate for the column density based on the local velocity gradient and the thermal velocity. This approach generally underestimates the photon escape probability because it neglects the density gradient and the actual shape of the cloud. We present a correction factor for the true line overlap in the Sobolev approximation and a new method based on local quantities, which fits the exact results reasonably well during the collapse of the cloud, with the error in the cooling rates always being less than 10%. Analytical fitting formulae fail at determining the photon escape probability after formation of the first protostar (error of ∼40%) because they are based on the assumption of spherical symmetry and therefore break down once a protostellar accretion disk has formed. Our method yields lower temperatures and hence promotes fragmentation for densities above ∼10{sup 10} cm{sup –3} at a distance of ∼200 AU from the first protostar. Since the overall accretion rates are hardly affected by the cooling implementation, we expect Pop III stars to have lower masses in our simulations, compared to the results of previous simulations that used the Sobolev approximation.

79 ASTRONOMY AND ASTROPHYSICS↗

The CO-to-H{sub 2} conversion factor across the Perseus molecular cloud

We derive the CO-to-H{sub 2} conversion factor, X {sub CO} = N(H{sub 2})/I {sub CO}, across the Perseus molecular cloud on sub-parsec scales by combining the dust-based N(H{sub 2}) data with the I {sub CO} data from the COMPLETE Survey. We estimate an average X {sub CO} ∼ 3 × 10{sup 19} cm{sup –2} K{sup –1} km{sup –1} s and find a factor of ∼3 variations in X {sub CO} between the five sub-regions in Perseus. Within the individual regions, X {sub CO} varies by a factor of ∼100, suggesting that X {sub CO} strongly depends on local conditions in the interstellar medium. We find that X {sub CO} sharply decreases at A{sub V} ≲ 3 mag but gradually increases at A{sub V} ≳ 3 mag, with the transition occurring at A{sub V} where I {sub CO} becomes optically thick. We compare the N(H I), N(H{sub 2}), I {sub CO}, and X {sub CO} distributions with two models of the formation of molecular gas, a one-dimensional photodissociation region (PDR) model and a three-dimensional magnetohydrodynamic (MHD) model, tracking both the dynamical and chemical evolution of gas. The PDR model based on the steady state and equilibrium chemistry reproduces our data very well but requires a diffuse halo to match the observed N(H I) and I {sub CO} distributions. The MHD model matches our data reasonably well, suggesting that time-dependent effects on H{sub 2} and CO formation are insignificant for an evolved molecular cloud like Perseus. However, we find interesting discrepancies, including a broader range of N(H I), likely underestimated I {sub CO}, and a large scatter of I {sub CO} at small A{sub V} . These discrepancies most likely result from strong compressions and rarefactions and density fluctuations in the MHD model.

79 ASTRONOMY AND ASTROPHYSICS↗

QUANTIFYING OBSERVATIONAL PROJECTION EFFECTS USING MOLECULAR CLOUD SIMULATIONS

The physical properties of molecular clouds are often measured using spectral-line observations, which provide the only probes of the clouds' velocity structure. It is hard, though, to assess whether and to what extent intensity features in position-position-velocity (PPV) space correspond to 'real' density structures in position-position-position (PPP) space. In this paper, we create synthetic molecular cloud spectral-line maps of simulated molecular clouds, and present a new technique for measuring the reality of individual PPV structures. Using a dendrogram algorithm, we identify hierarchical structures in both PPP and PPV space. Our procedure projects density structures identified in PPP space into corresponding intensity structures in PPV space and then measures the geometric overlap of the projected structures with structures identified from the synthetic observation. The fractional overlap between a PPP and PPV structure quantifies how well the synthetic observation recovers information about the three-dimensional structure. Applying this machinery to a set of synthetic observations of CO isotopes, we measure how well spectral-line measurements recover mass, size, velocity dispersion, and virial parameter for a simulated star-forming region. By disabling various steps of our analysis, we investigate how much opacity, chemistry, and gravity affect measurements of physical properties extracted from PPV cubes. For the simulations used here, which offer a decent, but not perfect, match to the properties of a star-forming region like Perseus, our results suggest that superposition induces a ∼40% uncertainty in masses, sizes, and velocity dispersions derived from {sup 13}CO (J = 1-0). As would be expected, superposition and confusion is worst in regions where the filling factor of emitting material is large. The virial parameter is most affected by superposition, such that estimates of the virial parameter derived from PPV and PPP information typically disagree by a factor of ∼2. This uncertainty makes it particularly difficult to judge whether gravitational or kinetic energy dominate a given region, since the majority of virial parameter measurements fall within a factor of two of the equipartition level α ∼ 2.

79 ASTRONOMY AND ASTROPHYSICS↗