HCN production in Titan's atmosphere: Coupling quantum chemistry and disequilibrium atmospheric modeling

This is the author accepted manuscript. the final version is available from IOP Publishing via the DOI in this recordHydrogen cyanide (HCN) is a critical reactive source of nitrogen for building key biomolecules relevant for the origin of life. Still, many HCN reactions remain uncharacterized by experiments and theory, and the complete picture of HCN production in planetary atmospheres is not fully understood. To improve this situation, we develop a novel technique making use of computational quantum chemistry, experimental data, and atmospheric numerical simulations. First, we use quantum chemistry simulations to explore the entire field of possible reactions for a list of primary species in N2-, CH4-, and H2-dominated atmospheres. In this process, we discover 33 new reactions with no previously known rate coefficients. From here, we develop a consistent reduced atmospheric hybrid chemical network (CRAHCN) containing experimental values when available and our calculated rate coefficients otherwise. Next, we couple CRAHCN to a 1D chemical kinetic model (ChemKM) to compute the HCN abundance as a function of atmospheric depth on Titan. Our simulated atmospheric HCN profile agrees very well with the Cassini observations. CRAHCN contains 104 reactions; however, nearly all of the simulated atmospheric HCN profile can be obtained using a scaled-down network of only 19 dominant reactions. From here, we form a complete picture of HCN chemistry in Titan's atmosphere, from the dissociation of the main atmospheric species, down to the direct production of HCN along four major channels. One of these channels was first discovered and characterized in Pearce et al. and this work.NSERCEuropean Union Horizon 202

The lower mass function of the young open cluster Blanco 1: from 30 Mjup to 3 Mo

We performed a deep wide field optical survey of the young (~100-150 Myr) open cluster Blanco1 to study its low mass population well down into the brown dwarf regime and estimate its mass function over the whole cluster mass range.The survey covers 2.3 square degrees in the I and z-bands down to I ~ z ~ 24 with the CFH12K camera. Considering two different cluster ages (100 and 150 Myr), we selected cluster member candidates on the basis of their location in the (I,I-z) CMD relative to the isochrones, and estimated the contamination by foreground late-type field dwarfs using statistical arguments, infrared photometry and low-resolution optical spectroscopy. We find that our survey should contain about 57% of the cluster members in the 0.03-0.6 Mo mass range, including 30-40 brown dwarfs. The candidate's radial distribution presents evidence that mass segregation has already occured in the cluster. We took it into account to estimate the cluster mass function across the stellar/substellar boundary. We find that, between 0.03Mo and 0.6Mo, the cluster mass distribution does not depend much on its exact age, and is well represented by a single power-law, with an index alpha=0.69 +/- 0.15. Over the whole mass domain, from 0.03Mo to 3Mo, the mass function is better fitted by a log-normal function with m0=0.36 +/- 0.07Mo and sigma=0.58 +/- 0.06. Comparison between the Blanco1 mass function, other young open clusters' MF, and the galactic disc MF suggests that the IMF, from the substellar domain to the higher mass part, does not depend much on initial conditions. We discuss the implications of this result on theories developed to date to explain the origin of the mass distribution.Comment: 18 pages, 15 figures and 5 tables accepted in A&

Black Hole Models of Quasars

Observations of active galactic nuclei are interpreted in terms of a theoretical model involving accretion onto a massive black hole. Optical quasars and Seyfert galaxies are associated with holes accreting near the Eddington rate and radio galaxies with sub-critical accretion. It is argued that magnetic fields are largely responsible for extracting energy and angular momentum from black holes and disks. Recent studies of electron-positron pair plasmas and their possible role in establishing the emergent X-ray spectrum are reviewed. The main evolutionary properties of active galactic nuclei can be interpreted in terms of a simple model in which black holes accrete gas at a rate dictated by the rate of gas supply which decreases with cosmic time. It may be worth searching for eclipsing binary black holes in lower power Seyferts

Overview of the Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) project

The Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) seeks to characterize 20 OB-dominated young clusters and their environs at distances d ≤ 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria and maps of nebular gas (particularly of hot X-ray-emitting plasma) and dust. A catalog of MYStIX Probable Complex Members with several hundred OB stars and 31,784 low-mass pre-main sequence stars is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of H II regions, and the survivability of protoplanetary disks in H II regions. This paper gives an introduction and overview of the project, covering the data analysis methodology and application to two star-forming regions: NGC 2264 and the Trifid Nebula. © 2013. The American Astronomical Society. All rights reserved.We thank J. Forbrich and P. Teixeira (Univ. Vienna) for useful discussion about NGC 2264. The MYStIX project is supported at Penn State by NASA grant NNX09AC74G, NSF grant AST-0908038, and theChandra ACIS Team contract SV4- 74018 (PIs: G. Garmire & L. Townsley), issued by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of NASA under contract NAS8-03060. M. S. Povich was supported by an NSF Astronomy and Astrophysics Postdoctoral Fellowship under award AST-0901646. This research made use of data products from the Chandra Data Archive and the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory (California Institute of Technology) under a contract with NASA. The United Kingdom Infrared Telescope is operated by the Joint Astronomy Centre on behalf of the Science and Technology Facilities Council of the U.K. This work is based in part on data obtained as part of the UKIRT Infrared Deep Sky Survey and in part on data obtained in UKIRT Director’s Discretionary Time. This research used data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. The HAWK-I near-infrared observations were collected with the High Acuity Wide-field K-band Imager instrument on the ESO 8 m Very Large Telescope at Paranal Observatory, Chile, under ESO programme 60.A-9284(K). This research has also made use of NASA’s Astrophysics Data System Bibliographic Services, the SIMBAD database operated at the Centre de Donnees ´ Astronomique de Strasbourg, and SAOImage DS9 software developed by Smithsonian Astrophysical Observatory

Radiatively heated, protoplanetary discs with dead zones. I. Dust settling and thermal structure of discs around M stars

The irradiation of protoplanetary discs by central stars is the main heating mechanism for discs, resulting in their flared geometric structure. In a series of papers, we investigate the deep links between 2D self-consistent disc structure and planetary migration in irradiated discs, focusing particularly on those around M stars. In this first paper, we analyse the thermal structure of discs that are irradiated by an M star by solving the radiative transfer equation by means of a Monte Carlo code. Our simulations of irradiated hydrostatic discs are realistic and self-consistent in that they include dust settling with multiple grain sizes (N=15), the gravitational force of an embedded planet on the disc, and the presence of a dead zone (a region with very low levels of turbulence) within it. We show that dust settling drives the temperature of the mid-plane from an $r^{-3/5}$ distribution (well mixed dust models) toward an $r^{-3/4}$. The dead zone, meanwhile, leaves a dusty wall at its outer edge because dust settling in this region is enhanced compared to the active turbulent disc at larger disc radii. The disc heating produced by this irradiated wall provides a positive gradient region of the temperature in the dead zone in front of the wall. This is crucially important for slowing planetary migration because Lindblad torques are inversely proportional to the disc temperature. Furthermore, we show that low turbulence of the dead zone is self-consistently induced by dust settling, resulting in the Kelvin-Helmholtz instability (KHI). We show that the strength of turbulence arising from the KHI in the dead zone is $\alpha=10^{-5}$.Comment: 19 pages, 20 figures, 3 tables, accepted for publication in MNRA

From clump to disc scales in W3 IRS4 A case study of the IRAM NOEMA large programme CORE

Context. High-mass star formation typically takes place in a crowded environment, with a higher likelihood of young forming stars affecting and being affected by their surroundings and neighbours, as well as links between different physical scales affecting the outcome. However, observational studies are often focused on either clump or disc scales exclusively. Aims. We explore the physical and chemical links between clump and disc scales in the high-mass star formation region W3 IRS4, a region that contains a number of different evolutionary phases in the high-mass star formation process, as a case-study for what can be achieved as part of the IRAM NOrthern Extended Millimeter Array (NOEMA) large programme named CORE: “Fragmentation and disc formation in high-mass star formation”. Methods. We present 1.4 mm continuum and molecular line observations with the IRAM NOEMA interferometer and 30 m telescope, which together probe spatial scales from ~0.3−20′′ (600−40 000 AU or 0.003−0.2 pc at 2 kpc, the distance to W3). As part of our analysis, we used XCLASS to constrain the temperature, column density, velocity, and line-width of the molecular emission lines. Results. The W3 IRS4 region includes a cold filament and cold cores, a massive young stellar object (MYSO) embedded in a hot core, and a more evolved ultra-compact (UC)H II region, with some degree of interaction between all components of the region that affects their evolution. A large velocity gradient is seen in the filament, suggesting infall of material towards the hot core at a rate of 10−3−10−4 M⊙ yr−1, while the swept up gas ring in the photodissociation region around the UCH II region may be squeezing the hot core from the other side. There are no clear indications of a disc around the MYSO down to the resolution of the observations (600 AU). A total of 21 molecules are detected, with the abundances and abundance ratios indicating that many molecules were formed in the ice mantles of dust grains at cooler temperatures, below the freeze-out temperature of CO (≲35 K). This contrasts with the current bulk temperature of ~50 K, which was obtained from H2CO. Conclusions. CORE observations allow us to comprehensively link the different structures in the W3 IRS4 region for the first time. Our results argue that the dynamics and environment around the MYSO W3 IRS4 have a significant impact on its evolution. This context would be missing if only high resolution or continuum observations were available

Obtaining Statistics of Turbulent Velocity from Astrophysical Spectral Line Data

Turbulence is a crucial component of dynamics of astrophysical fluids dynamics, including those of ISM, clusters of galaxies and circumstellar regions. Doppler shifted spectral lines provide a unique source of information on turbulent velocities. We discuss Velocity-Channel Analysis (VCA) and its offspring Velocity Coordinate Spectrum (VCS) that are based on the analytical description of the spectral line statistics. Those techniques are well suited for studies of supersonic turbulence. We stress that a great advantage of VCS is that it does not necessary require good spatial resolution. Addressing the studies of mildly supersonic and subsonic turbulence we discuss the criterion that allows to determine whether a traditional tool for such a research, namely, Velocity Centroids are dominated by density or velocity. We briefly discuss the use of higher order correlations as the means to study intermittency of turbulence. We discuss observational data available and prospects of the field.Comment: 12 pages, Invited Talk, Penetrating Bars Through Masks of Cosmic Dust, Pilanesberg National Park, South Afrika, 7 June-12 June 200

Multi-scale view of star formation in IRAS 21078+5211: From clump fragmentation to disk wind

In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively

IRAS 23385+6053: An embedded massive cluster in the making

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars

Disk fragmentation in high-mass star formation. High-resolution observations towards AFGL 2591-VLA 3

Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources