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.

The Physical and Chemical Structure of Hot Molecular Cores

We have made self-consistent models of the density and temperature profiles of the gas and dust surrounding embedded luminous objects using a detailed radiative transfer model together with observations of the spectral energy distribution of hot molecular cores. Using these profiles we have investigated the hot core chemistry which results when grain mantles are evaporated, taking into account the different binding energies of the mantle molecules, as well a model in which we assume that all molecules are embedded in water ice and have a common binding energy. We find that most of the resulting column densities are consistent with those observed toward the hot core G34.3+0.15 at a time around 10^4 years after central luminous star formation. We have also investigated the dependence of the chemical structure on the density profile which suggests an observational possibility of constraining density profiles from determination of the source sizes of line emission from desorbed molecules.

The present-day chemical composition of the LMC

High-resolution observations of five OB-type main-sequence stars in the Large Magellanic Cloud (LMC) have been obtained with the UCL echelle spectrograph on the 3.9-m Anglo-Australian Telescope. These spectra have been analysed using LTE model- atmosphere techniques, to derive stellar atmospheric parameters and chemical compositions. As these stars are located within the hydrogen burning main-sequence band, their surface abundances should reflect those of the present-day interstellar medium. Detailed line-by-line differential analyses have been undertaken relative to Galactic comparison stars. We conclude that there exists a general metal deficiency of - 0.31 +/- 0.04 dex within the LMC, and find no significant abundance variations between cluster and field stars. There is also tentative evidence to suggest a lower oxygen to iron abundance ratio, and an over-deficiency of magnesium relative to the other alpha-elements. These are discussed in terms of previous abundance analyses and models of discontinuous (or bursting) star formation within the LMC. Finally, there is some evidence to suggest a greater chemical enrichment of material within the H. region LH104.

The VLT-FLAMES survey of massive stars: stellar parameters and rotational velocities in NGC 3293, NGC 4755 and NGC 6611

An analysis is presented of VLT-FLAMES spectroscopy for three Galactic clusters, NGC3293, NGC4755 and NGC6611. Non-LTE model atmosphere calculations have been used to estimate effective temperatures (from either the helium spectrum or the silicon ionization equilibrium) and gravities (from the hydrogen spectrum). Projected rotational velocities have been deduced from the helium spectrum (for fast and moderate rotators) or the metal line spectrum (for slow rotators). The origin of the low gravity estimates for apparently near main sequence objects is discussed and is related to the stellar rotational velocity. The atmospheric parameters have been used to estimate cluster distances (which are generally in good agreement with previous determinations) and these have been used to estimate stellar luminosities and evolutionary masses. The observed Hertzsprung-Russell diagrams are compared with theoretical predictions and some discrepancies including differences in the main sequence luminosities are discussed. Cluster ages have been deduced and evidence for non-coeval star formation is found for all three of the clusters. Projected rotational velocities for targets in the older clusters, NGC3293 and NGC4755, have been found to be systematically larger than those for the field, confirming recent results in other similar age clusters. The distribution of projected rotational velocities are consistent with a Gaussian distribution of intrinsic rotational velocities. For the relatively unevolved targets in the older clusters, NGC3293 and NGC4755, the peak of the velocity distribution would be 250 km s(-1) with a full-width-half-maximum of approximately 180 km s(-1). For NGC6611, the sample size is relatively small but implies a lower mean rotational velocity. This may be evidence for the spin-down effect due to angular momentum loss through stellar winds, although our results are consistent with those found for very young high mass stars. For all three clusters we deduce present day mass functions with Gamma-values in the range of -1.5 to -1.8, which are similar to other young stellar clusters in the Milky Way.

A survey of [ D2CO] /[ H2CO] and [ N2D^+] /[ N2H^+] ratios towards protostellar cores

Aims.We use observations and models of molecular D/H ratios to probe the physical conditions and chemical history of the gas and to differentiate between gas-phase and grain-surface chemical processing in star forming regions. Methods: As a follow up to previous observations of HDCO/H2CO and DCN/HCN ratios in a selection of low-mass protostellar cores, we have measured D2CO/H2CO and N2D^+/N2H+ ratios in these same sources. For comparison, we have also measured N2D^+/N2H+ ratios towards several starless cores and have searched for N2D+ and deuterated formaldehyde towards hot molecular cores (HMCs) associated with high mass star formation. We compare our results with predictions from detailed chemical models, and to other observations made in these sources. Results: Towards the starless cores and low-mass protostellar sources we have found very high N2D+ fractionation, which suggests that the bulk of the gas in these regions is cold and heavily depleted. The non-detections of N2D+ in the HMCs indicate higher temperatures. We did detect HDCO towards two of the HMCs, with abundances 1-3% of H2CO. These are the first detections of deuterated formaldehyde in high mass sources since Turner (1990) measured HDCO/H2CO and D2CO/H2CO towards the Orion Compact Ridge. Figures 1-5 are only available in electronic form at http://www.aanda.org

The abundance of HNCO and its use as a diagnostic of environment

Aims. We aim to investigate the chemistry and gas phase abundance of HNCO and the variation of the HNCO/CS abundance ratio as a diagnostic of the physics and chemistry in regions of massive star formation. Methods. A numerical-chemical model has been developed which self-consistently follows the chemical evolution of a hot core. The model comprises of two distinct stages. The first stage follows the isothermal, modified free-fall collapse of a molecular dark cloud. This is immediately followed by an increase in temperature which represents the switch on of a central massive star and the subsequent evolution of the chemistry in a hot, dense gas cloud (the hot core). During the collapse phase, gas species are allowed to accrete on to grain surfaces where they can participate in further reactions. During the hot core phase surface species thermally desorb back in to the ambient gas and further chemical evolution takes place. For comparison, the chemical network was also used to model a simple dark cloud and photodissociation regions. Results. Our investigation reveals that HNCO is inefficiently formed when only gas-phase formation pathways are considered in the chemical network with reaction rates consistent with existing laboratory data. This is particularly true at low temperatures but also in regions with temperatures up to ~200 K. Using currently measured gas phase reaction rates, obtaining the observed HNCO abundances requires its formation on grain surfaces – similar to other “hot core” species such as CH3OH. However our model shows that the gas phase HNCO in hot cores is not a simple direct product of the evaporation of grain mantles. We also show that the HNCO/CS abundance ratio varies as a function of time in hot cores and can match the range of values observed. This ratio is not unambiguously related to the ambient UV field as been suggested – our results are inconsistent with the hypothesis of Martín et al. (2008, ApJ, 678, 245). In addition, our results show that this ratio is extremely sensitive to the initial sulphur abundance. We find that the ratio grows monotonically with time with an absolute value which scales approximately linearly with the S abundance at early times.

Discovery of Interstellar Anions in Cepheus and Auriga

We report the detection of microwave emission lines from the hydrocarbon anion C6H- and its parent neutral C6H in the star-forming region L1251A (in Cepheus), and the pre-stellar core L1512 (in Auriga). The carbon-chain-bearing species C4H, HC3N, HC5N, HC7N and C3S are also detected in large abundances. The observations of L1251A constitute the first detections of anions and long- chain polyynes and cyanopolyynes (with more than 5 carbon atoms) in the Cepheus Flare star- forming region, and the first detection of anions in the vicinity of a protostar outside of the Taurus molecular cloud complex, highlighting a wider importance for anions in the chemistry of star formation. Rotational excitation temperatures have been derived from the HC3N hyperfine structure lines, and are found to be 6.2 K for L1251A and 8.7 K for L1512. The anion-to-neutral ratios are 3.6% and 4.1%, respectively, which are within the range of values previously observed in the interstellar medium, and suggest a relative uniformity in the processes governing anion abundances in different dense interstellar clouds. This research contributes towards the growing body of evidence that carbon chain anions are relatively abundant in interstellar clouds throughout the Galaxy, but especially in the regions of relatively high density and high depletion surrounding pre-stellar cores and young, embedded protostars.

A Spitzer space telescope study of Sn 2002hh: An infrared echo from a Type IIP supernova

We present late-time ( 590 - 994 days) mid-IR photometry of the normal but highly reddened Type IIP supernova SN 2002hh. Bright, cool, slowly fading emission is detected from the direction of the supernova. Most of this flux appears not to be driven by the supernova event but instead probably originates in a cool, obscured star formation region or molecular cloud along the line of sight. We also show, however, that the declining component of the flux is consistent with an SN-powered IR echo from a dusty progenitor CSM. Mid-IR emission could also be coming from newly condensed dust and/or an ejecta/CSM impact, but their contributions are likely to be small. For the case of a CSM-IR echo, we infer a dust mass of as little as 0.036 M-. with a corresponding CSM mass of 3.6(0.01/ r(dg)) M-., where rdg is the dust-to-gas mass ratio. Such a CSM would have resulted from episodic mass loss whose rate declined significantly about 28,000 years ago. Alternatively, an IR echo from a surrounding, dense, dusty molecular cloud might also have been responsible for the fading component. Either way, this is the first time that an IR echo has been clearly identified in a Type IIP supernova. We find no evidence for or against the proposal that Type IIP supernovae produce large amounts of dust via grain condensation in the ejecta. However, within the CSM-IR echo scenario, the mass of dust derived implies that the progenitors of the most common of core-collapse supernovae may make an important contribution to the universal dust content.

SN 2009jf: a slow-evolving stripped-envelope core-collapse supernova

We present an extensive set of photometric and spectroscopic data for SN 2009jf, a nearby Type Ib supernova (SN), spanning from ˜20 d before B-band maximum to 1 yr after maximum. We show that SN 2009jf is a slowly evolving and energetic stripped-envelope SN and is likely from a massive progenitor (25-30 Msun). The large progenitor's mass allows us to explain the complete hydrogen plus helium stripping without invoking the presence of a binary companion. The SN occurred close to a young cluster, in a crowded environment with ongoing star formation. The spectroscopic similarity with the He-poor Type Ic SN 2007gr suggests a common progenitor for some SNe Ib and Ic. The nebular spectra of SN 2009jf are consistent with an asymmetric explosion, with an off-centre dense core. We also find evidence that He-rich Ib SNe have a rise time longer than other stripped-envelope SNe, however confirmation of this result and further observations are needed. This paper is based on observations with several telescopes, including NTT(184.D-1151), VLT-UT1(085.D-0750,386.D-0126), NOT, WHT, TNG, PROMPT, Ekar, Calar Alto and Liverpool Telescope.

Chemical evolution in the circumstellar structure of B5 IRS1

The B5 dark cloud has been identified as a site of low-mass star formation. We report a survey of a selection of the molecular species modelled by the B5 dynamical and chemical model of Charnley et al. at the positions of circumstellar HCN clumps in B5 IRS1. All of the key species observed yield either abundances or upper limits to abundances below both the standard and the predicted values, appearing to show evidence of depletion and/or destruction if the transitions observed are thermalized. Our results are not in good agreement with the model, and they bring into question the interpretation of the structure of B5 IRS1 proposed by Fuller et al. It was expected that HCN clump C might exhibit a higher excitation than HCN clump A, since it appeared to be located within the blueshifted molecular outflow. However, there is no significant difference observed between the two clumps, suggesting that the near-infrared and optical nebulosity is evidence of a reflection nebula rather than shocked material in the outflow. Finally, it is observed that our results are more consistent with gas-grain models than with those models excluding gas-grain interaction.

DENSE GAS IN NEARBY GALAXIES .9. A SURVEY FOR HNC

Ten detections and five tentative detections of hydrogen isocyanide (HNC) J=1-0 emission are reported from a survey including sixteen galaxies. Full maps are presented for the nuclear regions of NGC 253 and IC 342, partial maps for Maffei 2, M 82, and M 83. Toward IC 342, the HNC and HCO+ distributions differ from those observed in 12CO, 13CO, HCN, CS, and NH3. This is likely a consequence of the density structure. Relative HNC abundances are with 10(-10)-10(-9) much smaller than those measured in nearby dark clouds and appear to be slightly smaller than those in regions of massive star formation of the Galactic disk. This is consistent with the presence of dense warm gas or a frequent occurrence of shocks in the nuclear regions of the galaxies observed. As in prominent Galactic star forming regions, 3 mm HNC line emission tends to be weaker than the corresponding emission from HCN and HCO+. Toward Arp 220, however, the 3 mm HNC/HCN line intensity ratio is > 1. HNC/HCO+, HNC/CO, and HNC to 20 cm radio continuum luminosity ratios are also particularly large. A possible interpretation is the presence of cool quiescent gas outside the central region which contains the starburst. In the other ultraluminous galaxy observed, NGC 6240, X(HNC) 10 smaller than in Arp 220, demonstrating that the molecular composition in ultraluminous galaxies is far from being uniform.

PHOSPHORUS CHEMISTRY IN DENSE INTERSTELLAR CLOUDS

I have used recent laboratory studies on the reactions of the phosphorus hydride ions, PH(n)+ (n = 0-4) to construct a new model of phosphorus chemistry in interstellar clouds. I find that the non-detection of PN in cold, dark clouds in consistent with the chemical models only if the depletion of phosphorus in large, approximately 10(4) in TMC-1. Although the laboratory studies indicate that organo-phosphorus chains C(n)P can be formed, this large depletion precludes the detection of any phosphorus-bearing moleclues in cold clouds. However, in warm clouds associated with star formation, the depletion of phosphorus may be reduced. In this case one can reproduce the PN abundance toward Orion KL with a depletion factor of about 300. Interestingly, if the organo-phosphorus species are not destroyed by O atoms, I predict fractional abundances in Ori KL of between 10(-11) and 10(-10) for C(n)P (n = 2-4) and HCCP.

Core-collapse supernovae in low-metallicity environments and future all-sky transient surveys

Aims. Massive stars in low-metallicity environments may produce exotic explosions such as long-duration gamma-ray bursts and pair-instability supernovae when they die as core-collapse supernovae (CCSNe). Such events are predicted to be relatively common in the early Universe during the first episodes of star-formation. To understand these distant explosions it is vital to study nearby CCSNe arising in low-metallicity environments to determine if the explosions have different characteristics to those studied locally in high-metallicity galaxies. Many of the nearby supernova searches concentrate their efforts on high star-formation rate galaxies, hence biasing the discoveries to metal rich regimes. Here we determine the feasibility of searching for these CCSNe in metal-poor dwarf galaxies using various survey strategies.

The VLT-FLAMES survey of massive stars: wind properties and evolution of hot massive stars in the Large Magellanic Cloud

We have studied the optical spectra of a sample of 28 O- and early B-type stars in the Large Magellanic Cloud, 22 of which are associated with the young star forming region N11. Our observations sample the central associations of LH9 and LH10, and the surrounding regions. Stellar parameters are determined using an automated fitting method ( Mokiem et al. 2005), which combines the stellar atmosphere code fastwind ( Puls et al. 2005) with the genetic algorithm based optimisation routine PIKAIA ( Charbonneau 1995). We derive an age of 7.0 +/- 1.0 and 3.0 +/- 1.0 Myr for LH9 and LH10, respectively. The age difference and relative distance of the associations are consistent with a sequential star formation scenario in which stellar activity in LH9 triggered the formation of LH10. Our sample contains four stars of spectral type O2. From helium and hydrogen line fitting we find the hottest three of these stars to be similar to 49- 54 kK ( compared to similar to 45- 46 kK for O3 stars). Detailed determination of the helium mass fraction reveals that the masses of helium enriched dwarfs and giants derived in our spectroscopic analysis are systematically lower than those implied by non-rotating evolutionary tracks. We interpret this as evidence for efficient rotationally enhanced mixing leading to the surfacing of primary helium and to an increase of the stellar luminosity. This result is consistent with findings for SMC stars by Mokiem et al. ( 2006). For bright giants and supergiants no such mass discrepancy is found; these stars therefore appear to follow tracks of modestly or non-rotating objects. The set of programme stars was sufficiently large to establish the mass loss rates of OB stars in this Z similar to 1/2 Z(circle dot) environment sufficiently accurate to allow for a quantitative comparison with similar objects in the Galaxy and the SMC. The mass loss properties are found to be intermediate to massive stars in the Galaxy and SMC. Comparing the derived modified wind momenta D-mom as a function of luminosity with predictions for LMC metallicities by Vink et al. ( 2001) yields good agreement in the entire luminosity range that was investigated, i.e. 5.0