The Electron Temperature Gradient in the Galactic Disk

We derive the electron temperature gradient in the Galactic disk using a sample of HII regions that spans Galactocentric distances 0--17 kpc. The electron temperature was calculated using high precision radio recombination line and continuum observations for more than 100 HII regions. Nebular Galactocentric distances were calculated in a consistent manner using the radial velocities measured by our radio recombination line survey. The large number of nebulae widely distributed over the Galactic disk together with the uniformity of our data provide a secure estimate of the present electron temperature gradient in the Milky Way. Because metals are the main coolants in the photoionized gas, the electron temperature along the Galactic disk should be directly related to the distribution of heavy elements in the Milky Way. Our best estimate of the electron temperature gradient is derived from a sample of 76 sources for which we have the highest quality data. The present gradient in electron temperature has a minimum at the Galactic Center and rises at a rate of 287 +/- 46 K/kpc. There are no significant variations in the value of the gradient as a function of Galactocentric radius or azimuth. The scatter we find in the HII region electron temperatures at a given Galactocentric radius is not due to observational error, but rather to intrinsic fluctuations in these temperatures which are almost certainly due to fluctuations in the nebular heavy element abundances. Comparing the HII region gradient with the much steeper gradient found for planetary nebulae suggests that the electron temperature gradient evolves with time, becoming flatter as a consequence of the chemical evolution of the Milky Way's disk.Comment: 43 pages, 9 figures (accepted for publication in the ApJ

An open extensible tool environment for Event-B

Abstract. We consider modelling indispensable for the development of complex systems. Modelling must be carried out in a formal notation to reason and make meaningful conjectures about a model. But formal modelling of complex systems is a difficult task. Even when theorem provers improve further and get more powerful, modelling will remain difficult. The reason for this that modelling is an exploratory activity that requires ingenuity in order to arrive at a meaningful model. We are aware that automated theorem provers can discharge most of the onerous trivial proof obligations that appear when modelling systems. In this article we present a modelling tool that seamlessly integrates modelling and proving similar to what is offered today in modern integrated development environments for programming. The tool is extensible and configurable so that it can be adapted more easily to different application domains and development methods.

Stability conditions and Stokes factors

Let A be the category of modules over a complex, finite-dimensional algebra. We show that the space of stability conditions on A parametrises an isomonodromic family of irregular connections on P^1 with values in the Hall algebra of A. The residues of these connections are given by the holomorphic generating function for counting invariants in A constructed by D. Joyce.Comment: Very minor changes. Final version. To appear in Inventione

An 8.5 GHz Arecibo survey of Carbon Recombination Lines toward Ultra-compact \HII regions: Physical properties of dense molecular material

We report here on a survey of carbon recombination lines (RLs) near 8.5 GHz toward 17 ultra-compact \HII regions (\UCHII s). Carbon RLs are detected in 11 directions, indicating the presence of dense photodissociation regions (PDRs) associated with the \UCHII s. In this paper, we show that the carbon RLs provide important, complementary information on the kinematics and physical properties of the ambient medium near \UCHII s. Non-LTE models for the carbon line forming region are developed, assuming that the PDRs surround the \UCHII s, and we constrained the model parameters by multi-frequency RL data. Modeling shows that carbon RL emission near 8.5 GHz is dominated by stimulated emission and hence we preferentially observe the PDR material that is in front of the \UCHII continuum. We find that the relative motion between ionized gas and the associated PDR is about half that estimated earlier, and has an RMS velocity difference of 3.3 \kms. Our models also give estimates for the PDR density and pressure. We found that the neutral density of PDRs is typically $>$ 5 $\times$ 10$^5$ \cmthree and \UCHII s can be embedded in regions with high ambient pressure. Our results are consistent with a pressure confined \HII region model where the stars are moving relative to the cloud core. Other models cannot be ruled out, however. Interestingly, in most cases, the PDR pressure is an order of magnitude larger than the pressure of the ionized gas. Further investigation is needed to understand this large pressure difference.Comment: 28 pages, 7 figures, 5 tables (accepted for publication in ApJ

Integrating microbial ecology into ecosystem models: challenges and priorities

Microbial communities can potentially mediate feedbacks between global change and ecosystem function, owing to their sensitivity to environmental change and their control over critical biogeochemical processes. Numerous ecosystem models have been developed to predict global change effects, but most do not consider microbial mechanisms in detail. In this idea paper, we examine the extent to which incorporation of microbial ecology into ecosystem models improves predictions of carbon (C) dynamics under warming, changes in precipitation regime, and anthropogenic nitrogen (N) enrichment. We focus on three cases in which this approach might be especially valuable: temporal dynamics in microbial responses to environmental change, variation in ecological function within microbial communities, and N effects on microbial activity. Four microbially-based models have addressed these scenarios. In each case, predictions of the microbial-based models differ—sometimes substantially—from comparable conventional models. However, validation and parameterization of model performance is challenging. We recommend that the development of microbial-based models must occur in conjunction with the development of theoretical frameworks that predict the temporal responses of microbial communities, the phylogenetic distribution of microbial functions, and the response of microbes to N enrichment

On the S-matrix renormalization in effective theories

This is the 5-th paper in the series devoted to explicit formulating of the rules needed to manage an effective field theory of strong interactions in S-matrix sector. We discuss the principles of constructing the meaningful perturbation series and formulate two basic ones: uniformity and summability. Relying on these principles one obtains the bootstrap conditions which restrict the allowed values of the physical (observable) parameters appearing in the extended perturbation scheme built for a given localizable effective theory. The renormalization prescriptions needed to fix the finite parts of counterterms in such a scheme can be divided into two subsets: minimal -- needed to fix the S-matrix, and non-minimal -- for eventual calculation of Green functions; in this paper we consider only the minimal one. In particular, it is shown that in theories with the amplitudes which asymptotic behavior is governed by known Regge intercepts, the system of independent renormalization conditions only contains those fixing the counterterm vertices with $n \leq 3$ lines, while other prescriptions are determined by self-consistency requirements. Moreover, the prescriptions for $n \leq 3$ cannot be taken arbitrary: an infinite number of bootstrap conditions should be respected. The concept of localizability, introduced and explained in this article, is closely connected with the notion of resonance in the framework of perturbative QFT. We discuss this point and, finally, compare the corner stones of our approach with the philosophy known as ``analytic S-matrix''.Comment: 28 pages, 10 Postscript figures, REVTeX4, submitted to Phys. Rev.

3-He in the Milky Way Interstellar Medium: Ionization Structure

The cosmic abundance of the 3-He isotope has important implications for many fields of astrophysics. We are using the 8.665 GHz hyperfine transition of 3-He+ to determine the 3-He/H abundance in Milky Way HII regions and planetary nebulae. This is one in a series of papers in which we discuss issues involved in deriving accurate 3-He/H abundance ratios from the available measurements. Here we describe the ionization correction we use to convert the 3-He+/H+ abundance, y3+, to the 3-He/H abundance, y3. In principle the nebular ionization structure can significantly influence the y3 derived for individual sources. We find that in general there is insufficient information available to make a detailed ionization correction. Here we make a simple correction and assess its validity. The correction is based on radio recombination line measurements of H+ and 4-He+, together with simple core-halo source models. We use these models to establish criteria that allow us to identify sources that can be accurately corrected for ionization and those that cannot. We argue that this effect cannot be very large for most of the sources in our observational sample. For a wide range of models of nebular ionization structure we find that the ionization correction factor varies from 1 to 1.8. Although large corrections are possible, there would have to be a conspiracy between the density and ionization structure for us to underestimate the ionization correction by a substantial amount.Comment: 36 pages, 4 figures To appear Astrophysical Journal, 20 August 2007, vol 665, no

The Chemical Evolution Carousel of Spiral Galaxies : Azimuthal Variations of Oxygen Abundance in NGC1365

19 pages, 13 figures. Accepted to ApJThe spatial distribution of oxygen in the interstellar medium of galaxies is the key to understanding how efficiently metals that are synthesized in massive stars can be redistributed across a galaxy. We present here a case study in the nearby spiral galaxy NGC1365 using 3D optical data obtained in the TYPHOON Program. We find systematic azimuthal variations of the HII region oxygen abundance imprinted on a negative radial gradient. The 0.2 dex azimuthal variations occur over a wide radial range of 0.3 to 0.7 R25 and peak at the two spiral arms in NGC1365. We show that the azimuthal variations can be explained by two physical processes: gas undergoes localized, sub-kpc scale self-enrichment when orbiting in the inter-arm region, and experiences efficient, kpc scale mixing-induced dilution when spiral density waves pass through. We construct a simple chemical evolution model to quantitatively test this picture and find that our toy model can reproduce the observations. This result suggests that the observed abundance variations in NGC1365 are a snapshot of the dynamical local enrichment of oxygen modulated by spiral-driven, periodic mixing and dilution.Peer reviewedFinal Published versio

A 31 GHz Survey of Low-Frequency Selected Radio Sources

The 100-m Robert C. Byrd Green Bank Telescope (GBT) and the Owens Valley Radio Observatory (OVRO) 40-m radio telescope have been used to conduct a survey of 3165 known extragalactic radio sources over 143 square degrees of the sky. Target sources were selected from the NRAO VLA Sky Survey in fields observed by the Cosmic Background Imager (CBI); most are extragalactic active galactic nuclei (AGN) with 1.4 GHz flux densities of 3 to 10 mJy. The resulting 31 GHz catalogs are presented in full online. Using a Maximum-Likelihood analysis to obtain an unbiased estimate of the distribution of the 1.4 to 31 GHz spectral indices of these sources, we find a mean 31 to 1.4 GHz flux ratio of 0.110 +/- 0.003 corresponding to a spectral index of alpha=-0.71 +/- 0.01 (S ~ nu^alpha); 9.0 +/- 0.8 % of sources have alpha > -0.5 and 1.2 +/- 0.2 % have alpha > 0. By combining this spectral index distribution with 1.4 GHz source counts we predict 31 GHz source counts in the range 1 mJy < S_31 < 4 mJy, N(>S_31) = (16.7 +/- 1.7) deg^2 (S_31/1 mJy)^(-0.80 +/- 0.07). We also assess the contribution of mJy-level (S_1.4 < 3.4 mJy) radio sources to the 31 GHz CMB power spectrum, finding a mean power of ell (ell+1) C^src_ell/(2 pi) = 44 +/- 14 micro-Kelvin^2 and a 95% upper limit of 80 micro-Kelvin^2 at ell = 2500. Including an estimated contribution of 12 micro-Kelvin^2 from the population of sources responsible for the turn-up in counts below S_1.4 = 1 mJy this amounts to 21 +/- 7 % of what is needed to explain the CBI high-ell excess signal, 275 +/- 63 micro-Kelvin^2. These results are consistent with other measurements of the 31 GHz point source foreground.Comment: Replace with accepted version, incorporating referee comments. Significant revisions for clarity, added comparison to recent results. Quoted error bars on source counts are larger and the range of validity more limited; results on CBI excess unchange