The reader can find in the literature a lot of different techniques to study the dynamics of a given system and also, many suitable numerical integrators to compute them. Notwithstanding the recent work of Maffione et al. (2011a) for mappings, a detailed comparison among the widespread indicators of chaos in a general system is still lacking. Such a comparison could lead to select the most efficient algorithms given a certain dynamical problem. Furthermore, in order to choose the appropriate numerical integrators to compute them, more comparative studies among numerical integrators are also needed. This work deals with both problems. We first extend the work of Maffione et al. (2011) for mappings to the 2D H\’enon & Heiles (1964) potential, and compare several variational indicators of chaos: the Lyapunov Indicator (LI); the Mean Exponential Growth Factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI) and its generalized version, the Generalized Alignment Index (GALI); the Fast Lyapunov Indicator (FLI) and its variant, the Orthogonal Fast Lyapunov Indicator (OFLI); the Spectral Distance (D) and the Dynamical Spectras of Stretching Numbers (SSNs). We also include in the record the Relative Lyapunov Indicator (RLI), which is not a variational indicator as the others. Then, we test a numerical technique to integrate Ordinary Differential Equations (ODEs) based on the Taylor method implemented by Jorba & Zou (2005) (called taylor), and we compare its performance with other two well-known efficient integrators: the Prince & Dormand (1981) implementation of a Runge-Kutta of order 7-8 (DOPRI8) and a Bulirsch-St\”oer implementation. These tests are run under two very different systems from the complexity of their equations point of view: a triaxial galactic potential model and a perturbed 3D quartic oscillator.

The reader can find in the literature a lot of different techniques to study the dynamics of a given system and also, many suitable numerical integrators to compute them. Notwithstanding the recent work of Maffione et al. (2011a) for mappings, a detailed comparison among the widespread indicators of chaos in a general system is still lacking. Such a comparison could lead to select the most efficient algorithms given a certain dynamical problem. Furthermore, in order to choose the appropriate numerical integrators to compute them, more comparative studies among numerical integrators are also needed. This work deals with both problems. We first extend the work of Maffione et al. (2011) for mappings to the 2D H\’enon & Heiles (1964) potential, and compare several variational indicators of chaos: the Lyapunov Indicator (LI); the Mean Exponential Growth Factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI) and its generalized version, the Generalized Alignment Index (GALI); the Fast Lyapunov Indicator (FLI) and its variant, the Orthogonal Fast Lyapunov Indicator (OFLI); the Spectral Distance (D) and the Dynamical Spectras of Stretching Numbers (SSNs). We also include in the record the Relative Lyapunov Indicator (RLI), which is not a variational indicator as the others. Then, we test a numerical technique to integrate Ordinary Differential Equations (ODEs) based on the Taylor method implemented by Jorba & Zou (2005) (called taylor), and we compare its performance with other two well-known efficient integrators: the Prince & Dormand (1981) implementation of a Runge-Kutta of order 7-8 (DOPRI8) and a Bulirsch-St\”oer implementation. These tests are run under two very different systems from the complexity of their equations point of view: a triaxial galactic potential model and a perturbed 3D quartic oscillator.

The reader can find in the literature a lot of different techniques to study the dynamics of a given system and also, many suitable numerical integrators to compute them. Notwithstanding the recent work of Maffione et al. (2011a) for mappings, a detailed comparison among the widespread indicators of chaos in a general system is still lacking. Such a comparison could lead to select the most efficient algorithms given a certain dynamical problem. Furthermore, in order to choose the appropriate numerical integrators to compute them, more comparative studies among numerical integrators are also needed. This work deals with both problems. We first extend the work of Maffione et al. (2011) for mappings to the 2D H\’enon & Heiles (1964) potential, and compare several variational indicators of chaos: the Lyapunov Indicator (LI); the Mean Exponential Growth Factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI) and its generalized version, the Generalized Alignment Index (GALI); the Fast Lyapunov Indicator (FLI) and its variant, the Orthogonal Fast Lyapunov Indicator (OFLI); the Spectral Distance (D) and the Dynamical Spectras of Stretching Numbers (SSNs). We also include in the record the Relative Lyapunov Indicator (RLI), which is not a variational indicator as the others. Then, we test a numerical technique to integrate Ordinary Differential Equations (ODEs) based on the Taylor method implemented by Jorba & Zou (2005) (called taylor), and we compare its performance with other two well-known efficient integrators: the Prince & Dormand (1981) implementation of a Runge-Kutta of order 7-8 (DOPRI8) and a Bulirsch-St\”oer implementation. These tests are run under two very different systems from the complexity of their equations point of view: a triaxial galactic potential model and a perturbed 3D quartic oscillator.

We present the discovery and spectroscopic follow-up of a nearby late-type L dwarf (2M0614+3950), and two extremely wide very-low-mass binary systems (2M0525-7425AB and 2M1348-1344AB), resulting from our search for common proper motion pairs containing ultracool components in the Two Micron All Sky Survey (2MASS) and the Wide-field Infrared Survey Explorer (WISE) catalogs. The near-infrared spectrum of 2M0614+3950 indicates a spectral type L9 \pm 1 object residing at a distance of 26.1 \pm 1.3 pc. The optical spectrum of the 2M0525-7425 primary reveals an M3.0 \pm 0.5 dwarf, accompanied by a secondary previously classified as L2. The system has an angular separation of ~44″, equivalent to ~2000 AU at the 45.7 \pm 2.5 pc distance. Using optical and infrared spectra, respectively, we classify the components of 2M1348-1344AB as M4.5 \pm 0.5 and T6 \pm 1. The angular separation of ~68″ is equivalent to ~1300 AU at the distance of 19.2 \pm 0.9 pc. 2M1348-1344AB is one of only five very wide (separation > 1000 AU) systems containing late T dwarfs known to date.

We present low-resolution MMT Hectospec spectroscopy of 594 candidate Monoceros stream member stars. Based on strong color-magnitude diagram overdensities, we targeted three fields within the stream’s footprint, with 178 deg < l < 203 deg and -25 deg < b < 25 deg. By comparing the measured iron abundances with those expected from smooth Galactic components alone, we measure, for the first time, the spectroscopic metallicity distribution function for Monoceros. We find the stream to be chemically distinct from both the thick disk and halo, with [Fe/H] = -1, and do not detect a trend in the stream's metallicity with Galactic longitude. Passing from b = +25 deg to b = -25 deg the median Monoceros metallicity trends upward by 0.1 dex, though uncertainties in modeling sample contamination by the disk and halo make this a marginal detection. In each field, we find Monoceros to have an intrinsic [Fe/H] dispersion of 0.10-0.22 dex. From the CaII K line, we measure [Ca/Fe] for a subsample of metal poor program stars with -1.1 < [Fe/H] < -0.5. In two of three fields, we find calcium deficiencies qualitatively similar to previously reported [Ti/Fe] underabundances in Monoceros and the Sagittarius tidal stream. Further, using 90 spectra of thick disk stars in the Monoceros pointings with b ~ +/-25 deg, we detect a 0.22 dex north/south metallicity asymmetry coincident with known stellar density asymmetry at R_GC ~ 12 kpc and |Z| ~ 1.7 kpc. Our median Monoceros [Fe/H] = -1.0 and its relatively low dispersion naturally fit the expectation for an appropriately luminous M_V ~ -13 dwarf galaxy progenitor.

When a gravitating point mass moves subsonically through a magnetized and isothermal medium, the dynamical structure of the flow is studied far from the mass using a perturbation analysis. Analytical solutions for the first-order density and the velocity perturbations are presented. Validity of our solutions is restricted to the cases where the Alfven velocity in the ambient medium is less than the accretor’s velocity. The density field is less dense because of the magnetic effects according to the solutions and the dynamical friction force becomes lower as the strength of the magnetic field increases.

Ultraviolet observations of the QSO 3C 263 (zem = 0.652) with COS and FUSE reveal O VI absorption systems at z = 0.06342 and 0.14072 . WIYN multi-object spectrograph observations provide information about the galaxies associated with the absorbers. The multi-phase system at z = 0.06342 traces cool photoionized gas and warm collisionally ionized gas associated with a L ~ 0.31L* compact spiral emission line galaxy with an impact parameter of 63 kpc. The cool photoionized gas in the absorber is well modeled with log U ~ -2.6, log N(H) ~17.8, log n(H) ~ -3.3 and [Si/H] = -0.14\pm0.23. The collisionally ionized gas containing C IV and O VI probably arises in cooling shock heated transition temperature gas with log T ~ 5.5. The absorber is likely tracing circumgalactic gas enriched by gas ejected from the spiral emission line galaxy. The simple system at z = 0.14072 only contains O VI and broad and narrow H I. The O VI with b = 33.4\pm11.9 km s-1 is likely associated with the broad H I {\lambda}1215 absorption with b = 86.7\pm15.4 km s-1. The difference in Doppler parameters implies the detection of a very large column of warm gas with log T = 5.61(+0.16, -0.25), log N(H) = 19.54(+0.26, -0.44) and [O/H] = -1.48 (+0.46, -0.26). This absorber is possibly associated with a 1.6L* absorption line galaxy with an impact parameter of 617 kpc although an origin in warm filament gas or in the halo of a fainter galaxy is more likely.

We report the results of a near-infrared imaging study of a $7.8 \times 7.8$ arcmin$^2$ region centered on the 6.7 GHz methanol maser associated with the RCW 34 star forming region using the 1.4m IRSF telescope at Sutherland. A total of 1283 objects were detected simultaneously in J, H, and K for an exposure time of 10800 seconds. The J-H, H-K two-colour diagram revealed a strong concentration of more than 700 objects with colours similar to what is expected of reddened classical T Tauri stars. The distribution of the objects on the K {\it vs} J-K colour-magnitude diagram is also suggestive that a significant fraction of the 1283 objects is lower mass pre-main sequence stars. We also present the luminosity function for the subset of about 700 pre-main sequence stars and show that it suggests ongoing star formation activity for about $10^7$ years. An examination of the spatial distribution of the pre-main sequence stars shows that the fainter (older) part of the population is more dispersed over the observed region and the brighter (younger) subset is more concentrated around the position of the O8.5V star. This suggests that the physical effects of the O8.5V star and the two early B-type stars on the remainder of the cloud out of which they formed, could have played a role in the onset of the more recent episode of star formation in RCW 34.

We derive the magnitude of fluctuations in total synchrotron intensity in the Milky Way, both from observations and from theory under various assumption about the relation between cosmic rays and interstellar magnetic fields. Given the relative magnitude of the fluctuations in the Galactic magnetic field suggested by the Faraday rotation and polarization data, the observations are inconsistent with local energy equipartition between cosmic rays and magnetic fields. Our analysis of synchrotron fluctuations suggests that the distribution of cosmic rays is nearly uniform at scales of order $100\p$, in contrast to that of the interstellar magnetic field. A conservative upper limit on the relative variations in the cosmic ray number density is 0.2–0.4 at the scales of order 100\,pc. Our results are consistent with a mild anticorrelation between cosmic-ray and magnetic energy densities at these scales. Energy equipartition between cosmic rays and magnetic fields still may hold, but at scales exceeding 1\,kpc.

HIFI GOT C+ Galactic plane [CII] spectral survey has detected strong emission at the spiral arm tangencies. We use the unique viewing geometry of the Scutum-Crux (S-C) tangency near i = 30degs to detect the warm ionized medium (WIM) component traced by [CII] and to study the effects of spiral density waves on Interstellar Medium (ISM) gas. We compare [CII] velocity features with ancillary HI, 12CO and 13CO data near tangent velocities at each longitude to separate the cold neutral medium and the warm neutral + ionized components in the S-C tangency, then we identify [CII] emission at the highest velocities without any contribution from 12CO clouds, as WIM. We present the GOT C+ results for the S-C tangency. We interpret the diffuse and extended excess [CII] emission at and above the tangent velocities as arising in the electron-dominated warm ionized gas in the WIM. We derive an electron density in the range of 0.2 – 0.9 cm^-3 at each longitude, a factor of several higher than the average value from Halpha and pulsar dispersion. We interpret the excess [CII] in S-C tangency as shock compression of the WIM induced by the spiral density waves.

We study the star formation rate (SFR) – stellar mass (M*) relation in a self-consistent manner from 0 < z 10 at 1 < z < 1.5), 2) red star-forming galaxies with low levels of dust obscuration and low specific SFRs (11%), and 3) dusty, blue star-forming galaxies with high specific SFRs (7%). The remaining 28% comprises quiescent galaxies. Galaxies on the "normal" star formation sequence show strong trends of increasing dust attenuation with stellar mass and a decreasing specific SFR, with an observed scatter of 0.25 dex (0.17 dex intrinsic scatter). The dusty, blue galaxies reside in the upper envelope of the star formation sequence with remarkably similar spectral shapes at all masses, suggesting that the same physical process is dominating the stellar light. The red, low-dust star-forming galaxies may be in the process of shutting off and migrating to the quiescent population.

We study the star formation rate (SFR) – stellar mass (M*) relation in a self-consistent manner from 0 < z 10 at 1 < z < 1.5), 2) red star-forming galaxies with low levels of dust obscuration and low specific SFRs (12%), and 3) dusty, blue star-forming galaxies with high specific SFRs (6%). The remaining 29% comprises quiescent galaxies. Galaxies on the "normal" star formation sequence show strong trends of increasing dust attenuation and decreasing sSFR with stellar mass, with an observed scatter of 0.26 dex (0.19 dex intrinsic scatter). The dusty, blue galaxies reside in the upper envelope of the star formation sequence with remarkably similar spectral shapes at all masses, suggesting that the same physical process is dominating the stellar light. The red, low-dust star-forming galaxies may be in the process of shutting off and migrating to the quiescent population.

We analyse the kinematics of disc stars observed by the RAVE survey in and beyond the Solar neighbourhood.We detect significant overdensities in the velocity distributions using a technique based on the wavelet transform.We find that the main local kinematic groups are large scale features, surviving at least up to ~1 kpc from the Sun in the direction of anti-rotation, and also at ~700 pc below the Galactic plane.We also find that for regions located at different radii than the Sun, the known groups appear shifted in the velocity plane. For example, the Hercules group has a larger azimuthal velocity for regions inside the Solar circle and a lower value outside. We have also discovered a new group at (U, V) = (92,-22) km/s in the Solar neighbourhood and confirmed the significance of other previously found groups. Some of these trends detected for the first time are consistent with dynamical models of the effects of the bar and the spiral arms. More modelling is required to definitively characterise the non-axisymmetric components of our Galaxy using these groups.

We present results from integral field spectroscopy with PMAS. The observed field contains: five protoplanetary discs (also known as proplyds), the high-velocity jet HH 514 and a bowshock. Spatial distribution maps are obtained for different emission line fluxes, the c(H{\beta}) coefficient, electron densities and temperatures, ionic abundances of different ions from collisionally excited lines (CELs), C2+ and O2+ abundances from recombination lines (RLs) and the abundance discrepancy factor of O2+, ADF(O2+). We find that collisional de-excitation has a major influence on the line fluxes in the proplyds. If this is not properly accounted for then physical conditions deduced from commonly used line ratios will be in error, leading to unreliable chemical abundances for these objects. We obtain the intrinsic emission of the proplyds 177-341, 170-337 and 170-334 by a direct subtraction of the background emission, though the last two present some background contamination due to their small sizes. A detailed analysis of 177-341 spectra reveals the presence of high-density gas (3.8\times10^5 cm^-3) in contrast to the typical values observed in the background gas of the nebula (3800 cm^-3). We also explore how the background subtraction could be affected by the possible opacity of the proplyd. We construct a physical model for the proplyd 177-341 finding a good agreement between the predicted and observed line ratios. Finally, we find that the use of reliable physical conditions returns an ADF(O2+) about zero for the intrinsic spectra of 177-341, while the background emission presents the typical ADF(O2+) observed in the Orion Nebula. We conclude that the presence of high-density ionized gas is severely affecting the abundances determined from CELs and, therefore, those from RLs should be considered as a better approximation to the true abundances.

I report on progress in my ongoing work with Professor Jayant Murthy concerning the origin and nature of the diffuse ultraviolet background radiation over the sky. We have obtained and are reducing a vast trove of Voyager ultraviolet spectrometer observations of the diffuse background shortward of Lyman alpha, including for the first time measurements made from the outermost regions of the solar system, where noise from solar-system scattered (and then grating-scattered) solar Lyman alpha is lowest. Also, we have obtained and are investigating the complete set of GALEX observations of the diffuse ultraviolet background longward of Lyman alpha. Preliminary investigation appears to confirm that longward of Lyman alpha there exists a component of the diffuse ultraviolet background that is not dust-scattered starlight.

We present the images and kinematics of circumnuclear molecular gas from 100 pc scale down to 10 pc scale in nearby active galactic nuclei (AGNs) using the Submillimeter Array (SMA) and the Plateau de Bure Interferometer (PdBI). We have observed several nearby galaxies that host AGNs, such as the nearest radio galaxy Centaurus A (NGC 5128), the Seyfert 2 galaxy M51 (NGC 5194), the Seyfert 2 galaxy NGC 1068, the Seyfert 1 galaxy NGC 1097, and the Seyfert 2 / starburst composite galaxy NGC 4945, in CO lines to see whether the molecular gas distribution, kinematics, and physical conditions at 10 – 100 pc scale follows the AGN unified model or not. In 100 pc scale, most of the circumnuclear molecular gas shows smooth velocity gradient, suggesting a regular rotating feature, and also shows abnormal line ratios, suggesting the existence of active sources to make the circumnuclear molecular gas dense and/or warm conditions or abnormal chemical compositions. In 10 pc scale, on the other hand, the molecular gas kinematics shows various characteristics, some shows very disturbed kinematics such as a jet-entrained feature in the galaxies that have jets, but some still shows regular rotation feature in a galaxy that does not have obvious jets. These results indicate that the kinematics and physical/chemical conditions of the circumnuclear molecular gas at the scale less than 100 pc is highly affected by the AGN activities, and at this scale, there is no clear evidence of any unified feature seen in the circumnuclear molecular gas.

A debate has arisen regarding the importance of stationary versus eruptive mass loss for massive star evolution. The reason is that stellar winds have been found to be clumped, which results in the reduction of unclumped empirical mass-loss rates. Most stellar evolution models employ theoretical mass-loss rates which are already reduced by a moderate factor of ~2-3 compared to non-corrected empirical rates. A key question is whether these reduced rates are of the correct order of magnitude, or if they should be reduced even further, which would mean that the alternative of eruptive mass loss becomes necessary. Here we introduce the transition mass-loss rate (dM/dt)_trans between O and Wolf-Rayet (WR) stars. Its novelty is that it is model independent. All that is required is postulating the spectroscopic transition point in a given data-set, and determining the stellar luminosity, which is far less model dependent than the mass-loss rate. The transition mass-loss rate is subsequently used to calibrate stellar wind strength by its application to the Of/WNh stars in the Arches cluster. Good agreement is found with two alternative modelling/theoretical results, suggesting that the rates provided by current theoretical models are of the right order of magnitude in the ~50Msun mass range. Our results do not confirm the specific need for eruptive mass loss as Luminous Blue Variables, and current stellar evolution modelling for Galactic massive stars seems sound. Mass loss through alternative mechanisms might still become necessary at lower masses, and/or metallicities, and the quantification of alternative mass loss is desirable.

Context: Hydrogen peroxide (HOOH) was recently detected toward \rho Oph A. Subsequent astrochemical modeling that included reactions in the gas phase and on the surface of dust grains was able to explain the observed abundance, and highlighted the importance of grain chemistry in the formation of HOOH as an intermediate product in water formation. This study also predicted that the hydroperoxyl radical HO2, the precursor of HOOH, should be detectable. Aims: We aim at detecting the hydroperoxyl radical HO2 in \rho Oph A. Methods: We used the IRAM 30m and the APEX telescopes to target the brightest HO2 lines at about 130 and 260 GHz. Results: We detect five lines of HO2 (comprising seven individual molecular transitions). The fractional abundance of HO2 is found to be about 1e-10, a value similar to the abundance of HOOH. This observational result is consistent with the prediction of the above mentioned astrochemical model, and thereby validates our current understanding of the water formation on dust grains. Conclusions: This detection, anticipated by a sophisticated gas-grain chemical model, demonstrates that models of grain chemistry have improved tremendously and that grain surface reactions now form a crucial part of the overall astrochemical network.

We present an SPH parameter study of the dynamical effect of photoionization from O–type stars on star–forming clouds of a range of masses and sizes during the time window before supernovae explode. Our model clouds all have the same degree of turbulent support initially, the ratio of turbulent kinetic energy to gravitational potential energy being set to $E_{\rm kin}/|E_{\rm pot}|$=0.7. We allow the clouds to form stars and study the dynamical effects of the ionizing radiation from the massive stars or clusters born within them. We find that dense filamentary structures and accretion flows limit the quantities of gas that can be ionized, particularly in the higher density clusters. More importantly, the higher escape velocities in our more massive (10$^{6}$M$_{\odot}$) clouds prevent the HII regions from sweeping up and expelling significant quantities of gas, so that the most massive clouds are largely dynamically unaffected by ionizing feedback. However, feedback has a profound effect on the lower–density 10$^{4}$ and 10$^{5}$M$_{\odot}$ clouds in our study, creating vast evacuated bubbles and expelling tens of percent of the neutral gas in the 3Myr timescale before the first supernovae are expected to detonate, resulting in clouds highly porous to both photons and supernova ejecta.

The Milky Way’s dwarf spheroidal satellites include the nearest, smallest and least luminous galaxies known. They also exhibit the largest discrepancies between dynamical and luminous masses. This article reviews the development of empirical constraints on the structure and kinematics of dSph stellar populations and discusses how this phenomenology translates into constraints on the amount and distribution of dark matter within dSphs. Some implications for cosmology and the particle nature of dark matter are discussed, and some topics/questions for future study are identified.

RR Lyrae variables and the stellar constituents of globular clusters are employed to establish the cosmic distance scale and age of the universe. However, photometry for RR Lyrae variables in the globular clusters M3, M15, M54, M92, NGC2419, and NGC6441 exhibit a dependence on the clustercentric distance. For example, variables and stars positioned near the crowded high-surface brightness cores of the clusters may suffer from photometric contamination, which invariably affects a suite of inferred parameters (e.g., distance, color excess, absolute magnitude, etc.). The impetus for this study is to mitigate the propagation of systematic uncertainties by increasing awareness of the pernicious impact of contaminated and radial-dependent photometry.

Interactions of cosmic rays with the interstellar gas and radiation fields of the Milky Way provide the majority of the gamma rays observed by the Fermi Gamma Ray Space Telescope. In addition to the gas which is densely concentrated along the Galactic Disk, hydrodynamical simulations and observational evidence favor the presence of a halo of hot (T~10^6 K) ionized hydrogen (H_II), extending with non-negligible densities out to the virial radius of the Milky Way. We show that cosmic ray collisions with this circum-galactic gas should be expected to provide a significant flux of gamma rays, on the order of 10% of the observed isotopic gamma ray background at energies above 1 GeV. In addition, gamma rays originating from the extended H_II halos of other galaxies along a given line-of-sight should contribute to this background at a similar level.

We report a bimodality in the azimuthal angle distribution of gas around galaxies as traced by MgII absorption: Halo gas prefers to exist near the projected galaxy major and minor axes. The bimodality is demonstrated by computing the mean azimuthal angle probability distribution function using 88 spectroscopically confirmed MgII absorption-selected galaxies [W_r(2796)> 0.1A] and 35 spectroscopically confirmed non-absorbing galaxies [W_r(2796)<0.1A] imaged with HST and SDSS. The azimuthal angle distribution for non-absorbers is flat, indicating no azimuthal preference for gas characterized by W_r(2796)<0.1A. We find that blue star-forming galaxies clearly drive the bimodality. We compute an azimuthal angle dependent MgII absorption covering fraction and find that it is enhanced by as much as 20-30% along the major and minor axes. The equivalent width distribution for gas along the major axis is likely skewed toward weaker MgII absorption than for gas along the projected minor axis. These combined results are highly suggestive that the bimodality is driven by gas accreted along the galaxy major axis and outflowing along the galaxy minor axis. The opening angle of outflows is 2.5 times larger than for accreting gas. We find the probably of detecting outflows is 60%, implying that winds are more commonly observed. This scenario is consistent with ideas of galaxy evolution were star-forming galaxies accrete new gas reservoirs, forming new stars and producing winds, while red early-type galaxies exist passively due to a lack of new gas reservoirs to form new stars.

Using high resolution Chandra data, we report the presence of a weak X-ray point source coincident with the nucleus of NGC 4178, a late-type bulgeless disk galaxy known to have high ionization mid-infrared (mid-IR) lines typically associated with active galactic nuclei (AGNs). Although the faintness of this source precludes a direct spectral analysis, we are able to infer its basic spectral properties using hardness ratios. X-ray modeling, combined with the nuclear mid-IR characteristics, suggests that NGC 4178 may host a highly absorbed AGN accreting at a high rate with a bolometric luminosity on order of 10^43 ergs/s. The black hole mass estimate, based on our Chandra data and archival VLA data using the most recent fundamental plane relations is \sim 10^4 – 10^5 M\odot, possibly the lowest mass nuclear black hole currently known. There are also three off-nuclear sources, two with a similar brightness to the nuclear source at 36″ and 32″ from the center. As with the nuclear source, hardness ratios are used to estimate spectra for these two sources, and both are consistent with a simple power- law model with absorption. These two sources have X-ray luminosities of the order of \sim 10^38 ergs/s, which place them at the threshold between X-ray binaries and ultra-luminous X-ray sources (ULXs). The third off-nuclear source, located 49″ from the center, is the brightest source detected, with an X-ray luminosity of \sim 10^40 ergs/s. Its spectrum is well-fit with an absorbed power law model, suggesting that it is a ULX. We also fit its spectrum with the Bulk Motion Comptonization (BMC) model and suggest that this source is consistent with an intermediate-mass black hole (IMBH) of mass (6\times2)\times10^3 M\odot.

We present analysis of the UV-spectrum of the low-z AGN IRAS-F22456-5125 obtained with the Cosmic Origins Spectrograph on board the Hubble Space Telescope. The spectrum reveals six main kinematic components, spanning a range of velocities of up to 800 km s-1, which for the first time are observed in troughs associated with CII, CIV, NV, SiII, SiIII, SiIV and SIV. We also obtain data on the OVI troughs, which we compare to those available from an earlier FUSE epoch. Column densities measured from these ions allow us to derive a well-constrained photoionization solution for each outflow component. Two of these kinematic components show troughs associated with transitions from excited states of SiII\ and CII. The number density inferred from these troughs, in combination with the deduced ioinization parameter, allows us to determine the distance to these outflow components from the central source. We find these components to be at a distance of ~ 10 kpc. The distances and the number densities derived are consistent with the outflow being part of a galactic wind.

We present a series of new high-sensitivity and high-resolution radio-continuum images of M31 at \lambda=20 cm (\nu=1.4 GHz). These new images were produced by merging archived 20 cm radio-continuum observations from the Very Large Array (VLA) telescope. Images presented here are sensitive to rms=60 \mu Jy and feature high angular resolution (<10"). A complete sample of discrete radio sources have been catalogued and analysed across 17 individual VLA projects. We identified a total of 864 unique discrete radio sources across the field of M31. One of the most prominent regions in M31 is the ring feature for which we estimated total integrated flux of 706 mJy at \lambda=20 cm. We compare here, detected sources to those listed in Gelfand et al. (2004) at \lambda=92 cm and find 118 sources in common to both surveys. The majority (61%) of these sources exhibit a spectral index of \alpha <-0.6 indicating that their emission is predominantly non-thermal in nature. That is more typical for background objects.

New X-ray (XMM-Newton) and JHKs (OMM) observations for members of the star cluster Alessi 95, which Turner et al.(2012) discovered hosts the classical Cepheid SU Cas, were used in tandem with UCAC3 (proper motion) and 2MASS observations to determine precise cluster parameters: E(J-H)=0.08+-0.02 and d=405+-15 pc. The ensuing consensus among cluster, pulsation, IUE, and trigonometric distances (d=414+-5(se)+-10(sd) pc) places SU Cas in a select group of nearby fundamental Cepheid calibrators (Delta Cep, Zeta Gem). High-resolution X-ray observations may be employed to expand that sample as the data proved pertinent for identifying numerous stars associated with SU Cas. Acquiring X-ray observations of additional fields may foster efforts to refine Cepheid calibrations used to constrain H_0.

In order to determine if the material ablated from high-velocity clouds (HVCs) is a significant source of low-velocity high ions (C IV, N V, and O VI) such as those found in the Galactic halo, we simulate the hydrodynamics of the gas and the time-dependent ionization evolution of its carbon, nitrogen, and oxygen ions. Our suite of simulations examines the ablation of warm material from clouds of various sizes, densities, and velocities as they pass through the hot Galactic halo. The ablated material mixes with the environmental gas, producing an intermediate-temperature mixture that is rich in high ions and that slows to the speed of the surrounding gas. We find that the slow mixed material is a significant source of the low-velocity O VI that is observed in the halo, as it can account for at least ~1/3 of the observed O VI column density. Hence, any complete model of the high ions in the halo should include the contribution to the O VI from ablated HVC material. However, such material is unlikely to be a major source of the observed C IV, presumably because the observed C IV is affected by photoionization, which our models do not include. We discuss a composite model that includes contributions from HVCs, supernova remnants, a cooling Galactic fountain, and photoionization by an external radiation field. By design, this model matches the observed O VI column density. This model can also account for most or all of the observed C IV, but only half of the observed N V.

We present new Herschel-SPIRE imaging spectroscopy (194-671 microns) of the bright starburst galaxy M82. Covering the CO ladder from J=4-3 to J=13-12, spectra were obtained at multiple positions for a fully sampled ~ 3 x 3 arcminute map, including a longer exposure at the central position. We present measurements of 12CO, 13CO, [CI], [NII], HCN, and HCO+ in emission, along with OH+, H2O+ and HF in absorption and H2O in both emission and absorption, with discussion. We use a radiative transfer code and Bayesian likelihood analysis to model the temperature, density, column density, and filling factor of multiple components of molecular gas traced by 12CO and 13CO, adding further evidence to the high-J lines tracing a much warmer (~ 500 K), less massive component than the low-J lines. The addition of 13CO (and [CI]) is new and indicates that [CI] may be tracing different gas than 12CO. No temperature/density gradients can be inferred from the map, indicating that the single-pointing spectrum is descriptive of the bulk properties of the galaxy. At such a high temperature, cooling is dominated by molecular hydrogen. Photon-dominated region (PDR) models require higher densities than those indicated by our Bayesian likelihood analysis in order to explain the high-J CO line ratios, though cosmic-ray enhanced PDR models can do a better job reproducing the emission at lower densities. Shocks and turbulent heating are likely required to explain the bright high-J emission.

We use a composite gravitational galactic model consisting of a disk, a halo, a massive nucleus and a strong nuclear bar, in order to study the connections between global and local parameters in a realistic dynamical system. The local model is constructed from a two-dimensional perturbed harmonic oscillator and can be derived by expanding the global model in the vicinity of the central stable Lagrange equilibrium point. The frequencies of oscillations are not arbitrary, but they are connected with all the parameters involved with the global model. Moreover, the value of the local energy is also connected with the value of the global energy. Low and high energy stars in the global model display chaotic motion. Comparison with previous research reveals that the presence of the massive nucleus is responsible for the chaotic motion of the low energy stars. In the local motion, the low energy stars show interesting resonance phenomena, but the chaotic motion, if any, is negligible. On the contrary, the high energy stars do not show bounded motion in the local model. This is an indication of particular activity near the center of galaxies possessing massive nuclei.

We present a new publicly available tool (DustPol) aimed to model the polarised thermal dust emission. The module DustPol, which is publicly available, is part of the ARTIST (Adaptable Radiative Transfer Innovations for Submillimetre Telescopes) package, which also offers tools for modelling the polarisation of line emission together with a model library and a Python-based user interface. DustPol can easily manage analytical as well as pre-gridded models to generate synthetic maps of the Stokes I, Q, and U parameters. These maps are stored in FITS format which is straightforwardly read by the data reduction software used, e.g., by the Atacama Large Millimeter Array (ALMA). This turns DustPol into a powerful engine for the prediction of the expected polarisation features of a source observed with ALMA or the Planck satellite as well as for the interpretation of existing submillimetre observations obtained with other telescopes. DustPol allows the parameterisation of the maximum degree of polarisation and we find that, in a prestellar core, if there is depolarisation, this effect should happen at densities of 10^6 cm-3 or larger. We compare a model generated by DustPol with the observational polarisation data of the low-mass Class 0 object NGC 1333 IRAS 4A, finding that the total and the polarised emission are consistent.

We have used high resolution (~2.3″) observations of the local (D = 46 Mpc) luminous infrared galaxy Arp 299 to map out the physical properties of the molecular gas which provides the fuel for its extreme star formation activity. The 12CO J=3-2, 12CO J=2-1 and 13CO J=2-1 lines were observed with the Submillimeter Array and the short spacings of the 12CO J=2-1 and J=3-2 observations have been recovered using James Clerk Maxwell Telescope single dish observations. We use the radiative transfer code RADEX to estimate the physical properties (density, column density and temperature) of the different regions in this system. The RADEX solutions of the two galaxy nuclei, IC 694 and NGC 3690, are consistent with a wide range of gas components, from warm moderately dense gas with T_{kin} > 30 K and n(H_{2}) ~ 0.3 – 3 x 10^{3} cm^{-3} to cold dense gas with T_{kin} ~ 10-30 K and n(H_{2}) > 3 x 10^{3} cm^{-3}. The overlap region is shown to have a better constrained solution with T_{\rm{kin}}$ ~ 10-50 K and n(H_{2}) ~ 1-30 x 10^{3} cm^{-3}. We estimate the gas masses and star formation rates of each region in order to derive molecular gas depletion times. The depletion times of all regions (20-60 Myr) are found to be about 2 orders of magnitude lower than those of normal spiral galaxies. This rapid depletion time can probably be explained by a high fraction of dense gas on kiloparsec scales in Arp 299. We estimate the CO-to-H_{2} factor, \alpha_{co} to be 0.4 \pm 0.3 (3 x 10^{-4}/ x_{CO}) M_{sol} (K km s^{-1} pc^{2})^{-1} for the overlap region. This value agrees well with values determined previously for more advanced merger systems.

The gravitational properties of torus have been investigated. It is shown that the torus can be formed from the test particles orbiting in the gravitational field of the central mass. In this case the toroidal distribution is achieved due to significant spread of inclinations and eccentricities of their orbits. To investigate a self-gravity of the torus we considered the N-body problem for the torus located in the gravitational field of the central mass. It is shown that in the equilibrium state the cross-section of the torus has an oval shape with Gaussian density distribution. The dependence of obscuring efficiency as a function of torus inclination is found.

We present high angular resolution observations (0.5″x0.3″) carried out with the Submillimeter Array (SMA) toward the AFGL2591 high-mass star forming region. Our SMA images reveal a clear chemical segregation within the AFGL2591 VLA 3 hot core, where different molecular species (Type I, II and III) appear distributed in three concentric shells. This is the first time that such a chemical segregation is ever reported at linear scales <3000 AU within a hot core. While Type I species (H2S and 13CS) peak at the AFGL2591 VLA 3 protostar, Type II molecules (HC3N, OCS, SO and SO2) show a double-peaked structure circumventing the continuum peak. Type III species, represented by CH3OH, form a ring-like structure surrounding the continuum emission. The excitation temperatures of SO2, HC3N and CH3OH (185+-11 K, 150+-20 K and 124+-12 K, respectively) show a temperature gradient within the AFGL2591 VLA 3 envelope, consistent with previous observations and modeling of the source. By combining the H2S, SO2 and CH3OH images, representative of the three concentric shells, we find that the global kinematics of the molecular gas follow Keplerian-like rotation around a 40 Mo-star. The chemical segregation observed toward AFGL2591 VLA 3 is explained by the combination of molecular UV photo-dissociation and a high-temperature (~1000 K) gas-phase chemistry within the low extinction innermost region in the AFGL2591 VLA 3 hot core.

We present a comparative study of recent works on merger-timescales with dynamical friction and find a strong contrast between idealized/isolated mergers (Boylan-Kolchin et al. 2008) and mergers from a cosmological volume (Jiang et al. 2008). Our study measures the duration of mergers in a cosmological N-body simulation of dark matter, with emphasis on higher redshifts (z < 10) and a lower mass range. In our analysis we consider and compare two merger definitions; tidal disruption and coalescence. We find that the merger-time formula proposed by Jiang et al. (2008) describes our results well and conclude that cosmologically motivated merger-time formulae provide a more versatile and statistically robust approximation for practical applications such as semi-analytic/hybrid models.

We present a comparative study of recent works on merger-timescales with dynamical friction and find a strong contrast between idealized/isolated mergers (Boylan-Kolchin et al. 2008) and mergers from a cosmological volume (Jiang et al. 2008). Our study measures the duration of mergers in a cosmological N-body simulation of dark matter, with emphasis on higher redshifts (z < 10) and a lower mass range. In our analysis we consider and compare two merger definitions; tidal disruption and coalescence. We find that the merger-time formula proposed by Jiang et al. (2008) describes our results well and conclude that cosmologically motivated merger-time formulae provide a more versatile and statistically robust approximation for practical applications such as semi-analytic/hybrid models.

We present Herschel images from the HOBYS key program of the Eagle Nebula (M16) in the far-infrared and sub-millimetre, using the PACS and SPIRE cameras at 70{\mu}m, 160{\mu}m, 250{\mu}m, 350{\mu}m, 500{\mu}m. M16, home to the Pillars of Creation, is largely under the influence of the nearby NGC6611 high-mass star cluster. The Herschel images reveal a clear dust temperature gradient running away from the centre of the cavity carved by the OB cluster. We investigate the heating effect of NGC6611 on the entire M16 star-forming complex seen by Herschel including the diffuse cloud environment and the dense filamentary structures identified in this region. In addition, we interpret the three-dimensional geometry of M16 with respect to the nebula, its surrounding environment, and the NGC6611 cavity. The dust temperature and column density maps reveal a prominent eastern filament running north-south and away from the high-mass star-forming central region and the NGC6611 cluster, as well as a northern filament which extends around and away from the cluster. The dust temperature in each of these filaments decreases with increasing distance from the NGC6611 cluster, indicating a heating penetration depth of \sim 10 pc in each direction in 3 – 6 \times 10^{22} cm-2 column density filaments. We show that in high-mass star-forming regions OB clusters impact the temperature of future star-forming sites, modifying the initial conditions for collapse and effecting the evolutionary criteria of protostars developed from spectral energy distributions. Possible scenarios for the origin of the morphology seen in this region are discussed, including a western equivalent to the eastern filament, which was destroyed by the creation of the OB cluster and its subsequent winds and radiation.

The central region of the Milky Way provides a unique laboratory for a systematic, spatially-resolved population study of evolved massive stars of various types in a relatively high metallicity environment. We have conducted a multi-wavelength data analysis of 180 such stars or candidates, most of which were drawn from a recent large-scale HST/NICMOS narrow-band Pa-a survey, plus additional 14 Wolf-Rayet stars identified in earlier ground-based spectroscopic observations of the same field. The multi-wavelength data include broad-band IR photometry measurements from HST/NICMOS, SIRIUS, 2MASS, Spitzer/IRAC, and Chandra X-ray observations. We correct for extinctions toward individual stars, improve the Pa-a line equivalent width measurements, quantify the substantial mid-IR dust emission associated with WC stars, and find X-ray counterparts. In the process, we identify 10 foreground sources, some of which may be nearby cataclysmic variables. The WN stars in the Arches and Central clusters show correlations between the Pa-a equivalent width and the adjacent continuum emission. However, the WN stars in the latter cluster are systematically dimmer than those in the Arches cluster, presumably due to the different ages of the two clusters. In the EW-magnitude plot, WNL stars, WC stars and OB supergiants roughly fall into three distinct regions. We estimate that the dust mass associated with individual WC stars in the Quintuplet cluster can reach 1e-5 M, or more than one order of magnitude larger than previous estimates. Thus WC stars could be a significant source of dust in the galaxies of the early universe. Nearly half of the evolved massive stars in the GC are located outside the three known massive stellar clusters. The ionization of several compact HII regions can be accounted for by their enclosed individual evolved massive stars, which thus likely formed in isolation or in small groups.

We have measured electron-ion recombination for Fe XII forming Fe XI using a merged beams configuration at the heavy-ion storage ring TSR located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. The measured merged beams recombination rate coefficient (MBRRC) for collision energies from 0 to 1500 eV is presented. This work uses a new method for determining the absolute MBRRC based on a comparison of the ion beam decay rate with and without the electron beam on. For energies below 75 eV, the spectrum is dominated by dielectronic recombination (DR) resonances associated with 3s-3p and 3p-3d core excitations. At higher energies we observe contributions from 3-N’ and 2-N’ core excitations DR. We compare our experimental results to state-of-the-art multi-configuration Breit-Pauli (MCBP) calculations and find significant differences, both in resonance energies and strengths. We have extracted the DR contributions from the measured MBRRC data and transformed them into a plasma recombination rate coefficient (PRRC) for temperatures in the range of 10^3 to 10^7 K. We show that the previously recommended DR data for Fe XII significantly underestimate the PRRC at temperatures relevant for both photoionized plasmas (PPs) and collisionaly ionized plasmas (CPs). This is to be contrasted with our MCBP PRRC results which agree with the experiment to within 30% at PP temperatures and even better at CP temperatures. We find this agreement despite the disagreement shown by the detailed comparison between our MCBP and experimental MBRRC results. Lastly, we present a simple parameterized form of the experimentally derived PRRC for easy use in astrophysical modelling codes.

Observations of the dust emission using millimetre/submillimetre bolometer arrays can be contaminated by molecular line flux, such as flux from 12CO. As the brightest molecular line in the submillimetre, it is important to quantify the contribution of CO flux to the dust continuum bands. Conversion factors were used to convert molecular line integrated intensities to flux detected by bolometer arrays in mJy per beam. These factors were calculated for 12CO line integrated intensities to the SCUBA-2 850 {\mu}m and 450 {\mu}m bands. The conversion factors were then applied to HARP 12CO 3-2 maps of NGC 1333 in the Perseus complex and NGC 2071 and NGC 2024 in the Orion B molecular cloud complex to quantify the respective 12CO flux contribution to the 850 {\mu}m dust continuum emission. Sources with high molecular line contamination were analysed in further detail for molecular outflows and heating by nearby stars to determine the cause of the 12CO contribution. The majority of sources had a 12CO 3-2 flux contribution under 20 per cent. However, in regions of molecular outflows, the 12CO can dominate the source dust continuum (up to 79 per cent contamination) with 12CO fluxes reaching \sim 68 mJy per beam.

Context. The putative tori surrounding the accretion disks of active galactic nuclei (AGNs) play a fundamental role in the unification scheme of AGNs. Infrared long-baseline interferometry allows us to study the inner dust distribution in AGNs with unprecedented spatial resolution over a wide infrared wavelength range. Aims. Near- and mid-infrared interferometry is used to investigate the milli-arcsecond-scale dust distribution in the type 1.5 Seyfert nucleus of NGC 3783. Methods. We observed NGC 3783 with the VLTI/AMBER instrument in the K-band and compared our observations with models. Results. From the K-band observations, we derive a ring-fit torus radius of 0.74 +/- 0.23 mas or 0.16 +/- 0.05 pc. We compare this size with infrared interferometric observations of other AGNs and UV/optical-infrared reverberation measurements. For the interpretation of our observations, we simultaneously model our near- and mid-infrared visibilities and the SED with a temperature/density-gradient model including an additional inner hot 1400 K ring component.

We are carrying out multi-frequency radio continuum observations, using the Australia Telescope Compact Array, to systematically search for collimated ionized jets towards high-mass young stellar objects (HMYSOs). Here we report observations at 1.4, 2.4, 4.8 and 8.6 GHz, made with angular resolutions of about 7, 4, 2, and 1 arcsec, respectively, towards six objects of a sample of 33 southern HMYSOs thought to be in very early stages of evolution. The objects in the sample were selected from radio and infrared catalogs by having positive radio spectral indices and being luminous (L_bol > 20,000 L_sun), but underluminous in radio emission compared to that expected from its bolometric luminosity. This criteria makes the radio sources good candidates for being ionized jets. As part of this systematic search, two ionized jets have been discovered: one previously published and the other reported here. The rest of the observed candidates correspond to three hypercompact hii regions and two ultracompact hii regions. The two jets discovered are associated with two of the most luminous (70,000 and 100,000 Lsun) HMYSOs known to harbor this type of objects, showing that the phenomena of collimated ionized winds appears in the formation process of stars at least up to masses of ~ 20 M_sun and provides strong evidence for a disk-mediated accretion scenario for the formation of high-mass stars. From the incidence of jets in our sample, we estimate that the jet phase in high-mass protostars lasts for 40,000 yr.

We present a multi-wavelength photometric study of ~15,000 resolved stars in the nearby spiral galaxy M83 (NGC5236, D=4.61Mpc) based on Hubble Space Telescope Wide Field Camera 3 observations using four filters: F336W, F438W, F555W, and F814W. We select 50 regions (an average size of 260 pc by 280 pc) in the spiral arm and inter-arm areas of M83, and determine the age distribution of the luminous stellar populations in each region. This is accomplished by correcting for extinction towards each individual star by comparing its colors with predictions from stellar isochrones. We compare the resulting luminosity weighted mean ages of the luminous stars in the 50 regions with those determined from several independent methods, including the number ratio of red-to-blue supergiants, morphological appearance of the regions, surface brightness fluctuations, and the ages of clusters in the regions. We find reasonably good agreement between these methods. We also find that young stars are much more likely to be found in concentrated aggregates along spiral arms, while older stars are more dispersed. These results are consistent with the scenario that star formation is associated with the spiral arms, and stars form primarily in star clusters and then disperse on short timescales to form the field population. The locations of Wolf-Rayet stars are found to correlate with the positions of many of the youngest regions, providing additional support for our ability to accurately estimate ages. We address the effects of spatial resolution on the measured colors, magnitudes, and age estimates. While individual stars can occasionally show measurable differences in the colors and magnitudes, the age estimates for entire regions are only slightly affected.

We present a new, approximate method of calculating the column density of ammonia in mapping observations of the 23 GHz inversion lines. The temperature regime typically found in star forming regions allows for the assumption of a slowly varying partition function for ammonia. It is therefore possible to determine the column density using only the (J=1,K=1) inversion transition rather than the typical combination of the (1,1) and (2,2) transitions, with additional uncertainties comparable to or less than typical observational error. The proposed method allows column density and mass estimates to be extended into areas of lower signal to noise ratio. We show examples of column density maps around a number of cores in the W3 and Perseus star-forming regions made using this approximation, along with a comparison to the corresponding results obtained using the full two-transition approach. We suggest that this method is a useful tool in studying the distribution of mass around YSOs, particularly in the outskirts of the protostellar envelope where the (2,2) ammonia line is often undetectable on the short timescales necessary for large area mapping.

Several scenarios have been suggested to explain the phase-space distribution of the Milky Way (MW) satellite galaxies in a disc of satellites (DoS). To quantitatively compare these different possibilities, a new method analysing angular momentum directions in modelled data is presented. It determines how likely it is to find sets of angular momenta as concentrated and as close to a polar orientation as is observed for the MW satellite orbital poles. The method can be easily applied to orbital pole data from different models. The observed distribution of satellite orbital poles is compared to published angular momentum directions of subhalos derived from six cosmological state-of-the-art simulations in the Aquarius project. This tests the possibility that filamentary accretion might be able to naturally explain the satellite orbits within the DoS. For the most likely alignment of main halo and MW disc spin, the probability to reproduce the MW satellite orbital pole properties turns out to be less than 0.5 per cent in Aquarius models. Even an isotropic distribution of angular momenta has a higher likelihood to produce the observed distribution. The two Via Lactea cosmological simulations give results similar to the Aquarius simulations. Comparing instead with numerical models of galaxy-interactions gives a probability of up to 90 per cent for some models to draw the observed distribution of orbital poles from the angular momenta of tidal debris. This indicates that the formation as tidal dwarf galaxies in a single encounter is a viable, if not the only, process to explain the phase-space distribution of the MW satellite galaxies.

Observations of binaries in clusters tend to be of visual binaries with separations of 10s – 100s au. Such binaries are ‘intermediates’ and their destruction or survival depends on the exact details of their individual dynamical history. We investigate the stochasticity of the destruction of such binaries and the differences between the initial and processed populations using N-body simulations. We concentrate on Orion Nebula Cluster-like clusters, where the observed binary separation distribution ranges from 62 – 620 au. We find that, starting from the same initial binary population in statistically identical clusters, the number of intermediate binaries that are destroyed after 1 Myr can vary by a factor of >2, and that the resulting separation distributions can be statistically completely different in initially substructured clusters. We also find that the mass ratio distributions are altered (destroying more low mass ratio systems), but not as significantly as the binary fractions or separation distributions. We conclude that finding very different intermediate (visual) binary populations in different clusters does not provide conclusive evidence that the initial populations were different.

Using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope, we have obtained high-resolution ultraviolet spectra of the C I absorption toward two stars behind the Local Leo Cold Cloud (LLCC). At a distance (~20 pc) that places it well inside the Local Bubble, the LLCC is the nearest example of the coldest known (T~20 K) diffuse interstellar clouds. The STIS measurements of the C I fine-structure excitation toward HD 85259 and HD 83023 indicate that the thermal gas pressure of the LLCC is much greater than that of the warm clouds in the Local Bubble. The mean LLCC pressure measured toward these two stars (60,000 cm$^{-3}$ K) implies an H I density of ~3000 cm$^{-3}$ and a cloud thickness of ~200 AU at the 20 K cloud temperature. Such a thin, cold, dense structure could arise at the collision interface between converging flows of warm gas. However, the measured LLCC pressure is appreciably higher than that expected in the colliding cloud interpretation given the velocity and column density constraints on warm clouds in the HD 85259 and HD 83023 sightlines. Additional STIS measurements of the Zn II, Ni II, and Cr II column densities toward HD 85259 indicate that the LLCC has a modest “warm cloud” dust depletion pattern consistent with its low dust-to-gas ratio determined from H I 21 cm and 100 micron observations. In support of the inferred sheet-like geometry for the LLCC, a multi-epoch comparison of the Na I absorption toward a high-proper-motion background star reveals a 40% column density variation indicative of LLCC Na I structure on a scale of ~50 AU.

Massive AGN-driven outflows are invoked by AGN-galaxy co-evolutionary models to suppress both star formation and black hole accretion. Massive molecular outflows have recently been revealed in some AGN hosts. However, the physical properties and structure of these AGN-driven molecular outflows are still poorly constrained. Here we present new IRAM PdBI observations of Mrk231, the closest quasar known, targeting both the CO(1-0) and CO(2-1) transitions. We detect broad wings in both transitions, tracing a massive molecular outflow with velocities up to 800 km/s. The wings are spatially resolved at high significance level (5-11 sigma), indicating that the molecular outflow extends on the kpc scale. The CO(2-1)/CO(1-0) ratio of the red broad wings is consistent with the ratio observed in the narrow core, while the blue broad wing is less excited than the core. The latter result suggests that quasar driven outflow models invoking shocks (which would predict higher gas excitation) are not appropriate to describe the bulk of the outflow in Mrk231. However, we note that within the central 700 pc the CO(2-1)/CO(1-0) ratio of the red wing is slightly, but significantly, higher than in the line core, suggesting that shocks may play a role in the central region. We also find that the average size of the outflow anticorrelates with the critical density of the transition used as a wind tracer. This indicates that, although diffuse and dense clumps coexist in the outflowing gas, dense outflowing clouds have shorter lifetime and that they evaporate into the diffuse component along the outflow or, more simply, that diffuse clouds are more efficiently accelerated to larger distances by radiation pressure.

I’m trying to summarize the science communicated via oral presentations and by posters at the IAU Symposium 287 “Cosmic Masers – from OH to H_0″, which took place from January 29 to February 3, 2012 in Stellenbosch, South Africa.

I’m trying to summarize the science communicated via oral presentations and by posters at the IAU Symposium 285 “Cosmic Masers – from OH to H_0″, which took place from January 29 to February 3, 2012 in Stellenbosch, South Africa.

We have analysed optical, near-, and mid-IR images of the bright planetary nebula (PN) IC 418. These images probe unambiguously for the first time a number of low surface brightness structures and shells around the bright main nebula, including radial filaments or rays, a system of three concentric rings, and two detached haloes with sizes ~ 150″ and 220″\times250″, respectively. The main nebula is slightly off-centered with respect to the elliptical outer halo. The time-lapse between the two haloes is 10,000-50,000 yr, whereas the time-lapse between the three concentric rings is ~ 630 yr. We emphasize the advantages of near- and mid-IR imaging for the detection of faint structures and shells around optically bright nebulae.

A search for emission lines in foreground galaxies in quasar spectra (z(gal) < z(QSO)) of the Sloan Digital Sky Survey (SDSS) data release 5 (DR5) reveals 23 examples of quasars shining through low redshift, foreground galaxies at small impact parameters (< 10 kpc). About 74,000 quasar spectra were examined by searching for narrow H{\alpha} emission lines at z < 0.38, at a flux level greater than 5 \times 10^-17 ergs cm^-2 s^-1, then confirming that other expected emission lines of the H II regions in the galaxy are detected. The galaxies were deblended from the quasar images to get colors and morphologies. For cases that allow the galaxy and the quasar to be deblended, the galaxies are blue (0.95 <(u-r)< 1.95). Extinction and reddening through the galaxies is determined from the (g-i) color excesses of the quasars. These reddening values are compared with the flux ratio of H{\alpha} to H{\beta}, which reflect the extinction for an undetermined fraction of the sightline through each galaxy. No trends were found relating E(B-V)_(g-i), impact parameter (b), and (u-r) for the galaxies or between E(B-V) derived from (g-i) and that derived from H{\alpha}/H{\beta}. Comparison with previous studies of quasar absorption systems indicate our sample is more reddened, suggesting disk-dominated absorber galaxies. Measurement or limits on galactic, interstellar Ca II and Na I absorption lines are given from the quasar spectrum. No trends were found relating Ca II equivalent width (W (Ca II)) or Na I equivalent width (W (Na I)) to b, but a correlation of r_s = -0.77 ({\alpha} = 0.05) was found relating W (Ca II) and E(B-V)(g-i) .

Stellar streams provide unique probes of galactic potentials, with the longer streams normally providing the cleaner measurements. In this paper, we show an example of a short tidal stream that is particularly sensitive to the shape of the Milky Way’s dark matter halo: the globular cluster tidal stream NGC 5466. This stream has an interesting deviation from a smooth orbit at its western edge. We show that such a deviation favours an underlying oblate or triaxial halo (irrespective of plausible variations in the Milky Way disc properties and the specific halo parametrisation chosen); spherical or prolate halo shapes can be excluded at a high confidence level. Therefore, more extensive data sets along the NGC 5466 tidal stream promise strong constraints on the Milky Way halo shape.

The Alcock-Paczynski (AP) effect uses the fact that, when analyzed with the correct geometry, we should observe structure that is statistically isotropic in the Universe. For structure undergoing cosmological expansion with the background, this constrains the product of the Hubble parameter and the angular diameter distance. However, the expansion of the Universe is inhomogeneous and local curvature depends on density. We argue that this distorts the AP effect on small scales. After analyzing the dynamics of galaxy pairs in the Millennium simulation, we find an interplay between peculiar velocities, galaxy properties and local density that affects how pairs trace cosmological expansion. We find that only low mass, isolated galaxy pairs trace the average expansion with a minimum “correction” for peculiar velocities. Other pairs require larger, more cosmology and redshift dependent peculiar velocity corrections and, in the small-separation limit of being bound in a collapsed system, do not carry cosmological information.

We have observed a sample of 19 carbon stars in the Sculptor, Carina, Fornax, and Leo I dwarf spheroidal galaxies with the Infrared Spectrograph on the Spitzer Space Telescope. The spectra show significant quantities of dust around the carbon stars in Sculptor, Fornax, and Leo I, but little in Carina. Previous comparisons of carbon stars with similar pulsation properties in the Galaxy and the Magellanic Clouds revealed no evidence that metallicity affected the production of dust by carbon stars. However, the more metal-poor stars in the current sample appear to be generating less dust. These data extend two known trends to lower metallicities. In more metal-poor samples, the SiC dust emission weakens, while the acetylene absorption strengthens. The bolometric magnitudes and infrared spectral properties of the carbon stars in Fornax are consistent with metallicities more similar to carbon stars in the Magellanic Clouds than in the other dwarf spheroidals in our sample. A study of the carbon budget in these stars reinforces previous considerations that the dredge-up of sufficient quantities of carbon from the stellar cores may trigger the final superwind phase, ending a star’s lifetime on the asymptotic giant branch.

We have observed a sample of 19 carbon stars in the Sculptor, Carina, Fornax, and Leo I dwarf spheroidal galaxies with the Infrared Spectrograph on the Spitzer Space Telescope. The spectra show significant quantities of dust around the carbon stars in Sculptor, Fornax, and Leo I, but little in Carina. Previous comparisons of carbon stars with similar pulsation properties in the Galaxy and the Magellanic Clouds revealed no evidence that metallicity affected the production of dust by carbon stars. However, the more metal-poor stars in the current sample appear to be generating less dust. These data extend two known trends to lower metallicities. In more metal-poor samples, the SiC dust emission weakens, while the acetylene absorption strengthens. The bolometric magnitudes and infrared spectral properties of the carbon stars in Fornax are consistent with metallicities more similar to carbon stars in the Magellanic Clouds than in the other dwarf spheroidals in our sample. A study of the carbon budget in these stars reinforces previous considerations that the dredge-up of sufficient quantities of carbon from the stellar cores may trigger the final superwind phase, ending a star’s lifetime on the asymptotic giant branch.

We report the discovery of a new ultra-faint globular cluster in the constellation of Ursa Minor, based on stellar photometry from the MegaCam imager at the Canada-France-Hawaii Telescope (CFHT). We find that this cluster, Munoz 1, is located at a distance of 45 +/- 5 kpc and at a projected distance of only 45 arcmin from the center of the Ursa Minor dSph galaxy. Using a Maximum Likelihood technique we measure a half-light radius of 0.5 arcmin, or equivalently 7 pc and an ellipticity consistent with being zero. We estimate its absolute magnitude to be M_V=-0.4 +/- 0.9, which corresponds to L_V=120 (+160, -65) L_sun and we measure a heliocentric radial velocity of -137 +/- 4 km/s based on Keck/DEIMOS spectroscopy. This new satellite is separate from Ursa Minor by ~30 kpc and 110 km/s suggesting the cluster is not obviously associated with the dSph, despite the very close angular separation. Based on its photometric properties and structural parameters we conclude that Munoz 1 is a new ultra-faint stellar cluster. Along with Segue 3 this is one of the faintest stellar clusters known to date.

Theoretical studies of structure formation find an inverse proportionality between the concentration of dark matter haloes and virial mass. This trend has been recently confirmed for virial masses Mvir > ~6e12 Msun by the observation of the X-ray emission from the hot halo gas. We present an alternative approach to this problem, exploring the concentration of dark matter haloes over galaxy scales on a sample of 18 early-type systems. Our c-Mvir relation is consistent with the X-ray analysis, extending towards lower virial masses, covering the range from ~4e11 Msun up to 5e12 Msun. A combination of the lensing analysis along with photometric data allows us to constrain the baryon fraction within a few effective radii, which is compared with prescriptions for adiabatic contraction (AC) of the dark matter haloes. We find that the standard methods for AC are strongly disfavored, requiring additional mechanisms — such as mass loss during the contraction process — to play a role during the phases following the collapse of the haloes.

We present the discovery of a source with broadband infrared photometric characteristics similar to Sakurai’s Object. WISE J180956.27-330500.2 (hereafter, J1810-3305) shows very red WISE colors, but a very blue 2MASS [K] vs. WISE [W1 (3.4 micron)] color. It was not visible during the IRAS era, but now has a 12 micron flux well above the IRAS point source catalog detection limit. There are also indications of variability in historical optical photographic plates, as well as in multi-epoch AKARI mid-infrared measurements. The broadband infrared spectral energy distribution, post-IRAS brightening and multiwavelength variability are all characteristics also shared by Sakurai’s Object – a post asymptotic giant branch (post-AGB) star which underwent a late thermal pulse and recently ejected massive envelopes of dust that are currently expanding and cooling. Optical progenitor colors suggest that J1810-3305 may have been of late spectral class. Its dramatic infrared brightening, and the detection of a late-type optical counterpart are consistent with a scenario in which we have caught an extremely massive dust ejection event (in 1998 or shortly before) during the thermal pulse of an AGB star, thus providing a unique opportunity to observe stellar evolution in this phase. J1810-3305 is the only source in the entire WISE preliminary data release with similar infrared SED and variability, emphasizing the rarity of such sources. Confirmation of its nature is of great importance.

This work aims to search for the observational footprints of the interactions between the interstellar medium (ISM) and stellar winds in the Wolf-Rayet (WR) nebula NGC6888 in order to understand its ionization structure, chemical composition, and kinematics. We have collected a set of integral field spectroscopy observations across NGC6888, obtained with PPAK in the optical range performing both 2D and 1D analyses. Attending to the 2D analysis in the northeast part of NGC6888, we have generated maps of the extinction structure and electron density. We produced statistical frequency distributions of the radial velocity and diagnostic diagrams. We have found that the spectra of a localized region to the southwest of this pointing can be represented well by shock models. Furthermore, we performed a thorough study of integrated spectra in nine regions over the whole nebula. We derived electron densities ranging from <100 to 360 cm^(-3). The electron temperature varies from ~7700 K to ~10200 K. A strong variation of up to a factor 10 between different regions in the nitrogen abundance has been found: N/H appears lower than the solar abundance in those positions observed at the edges and very enhanced in the observed inner parts. Oxygen appears slightly underabundant with respect to solar value, whereas the helium abundance is found to be above it. Finally, we provide a scenario for the evolution of NGC6888 to explain the features observed. This scheme consists of a structure of multiple shells.

We present continuum high resolution Submillimeter Array (SMA) observations of the transition disk object RX J1633.9-2442, which is located in the Ophiuchus molecular cloud and has recently been identified as a likely site of ongoing giant planet formation. The observations were taken at 340 GHz (880 micron) with the SMA in its most extended configuration, resulting in an angular resolution of 0.3″ (35 AU at the distance of the target). We find that the disk is highly inclined (i ~50 deg) and has an inner cavity ~25 AU in radius, which is clearly resolved by our observations. We simultaneously model the entire optical to millimeter wavelength spectral energy distribution (SED) and SMA visibilities of RX J1633.9-2442 in order to constrain the structure of its disk. We find that an empty cavity ~25 AU in radius is inconsistent with the excess emission observed at 12, 22, and 24 micron. Instead, the mid-IR excess can be modeled by either a narrow, optically thick ring at ~10 AU or an optically thin region extending from ~7 AU to ~25 AU. The inner disk (r < 5 AU) is mostly depleted of small dust grains as attested by the lack of detectable near-IR excess. We also present deep Keck aperture masking observations in the near-IR, which rule out the presence of a companion up to 500 times fainter than the primary star (in K-band) for projected separations in the 5-20 AU range. We argue that the complex structure of the RX J1633.9-2442 disk is best explained by multiple planets embedded within the disk. We also suggest that the properties and incidence of objects such as RX J1633.9-2442, T Cha, and LkCa 15 (and those of the companions recently identified to these two latter objects) are most consistent with the runaway gas accretion phase of the core accretion model, when giant planets gain their envelopes and suddenly become massive enough to open wide gaps in the disk.

We have used a combination of broad-band near-infrared and optical Johnson-Cousins photometry to study the dust properties in the line of sight to the Galactic globular cluster M4. We have investigated the reddening effects in terms of absolute strength and variation across the cluster field, as well as the form of the reddening law defined by the type of dust. Here, we determine the ratio of absolute to selective extinction (R_V) in the line of sight towards M4, which is known to be a useful indicator for the type of dust and therefore characterizes the applicable reddening law. Our method is independent of age assumptions and appears to be significantly more precise and accurate than previous approaches. We obtain A_V/E(B-V)=3.76\pm0.07 (random error) for the dust in the line of sight to M4 for our set of filters. With this value, the distance to M4 is found to be 1.80\pm0.05 kpc (random error). A reddening map for M4 has been created, which reveals a spatial differential reddening of delta E(B-V)>0.2 mag across a field within 10′ around the cluster centre and a total mean reddening of E(B-V)=0.37\pm0.01. In order to provide accurate zero points for the extinction coefficients of our photometric filters, we investigated the impact of stellar parameters such as temperature, surface gravity and metallicity on the extinction properties in different bandpasses. Using both synthetic ATLAS9 spectra and observed spectral energy distributions, we found similar sized effects for the range of temperature, surface gravity, and metallicity typical of globular cluster stars: each causes a change of about 3% in the necessary correction factor for each filter combination. From our calculations, we provide extinction zero points for Johnson-Cousins and 2MASS filters, spanning a wide range of stellar parameters and dust types, suited for accurate, object-specific extinction corrections.

We report the results of a search for molecular oxygen (O2) toward the Orion Bar, a prominent photodissociation region at the southern edge of the HII region created by the luminous Trapezium stars. We observed the spectral region around the frequency of the O2 N_J = 3_3 – 1_2 transition at 487 GHz and the 5_4 – 3_4 transition at 774 GHz using the Heterodyne Instrument for the Far Infrared on the Herschel Space Observatory. Neither line was detected, but the 3sigma upper limits established here translate to a total line-of-sight O2 column density < 1.5 10^16 cm^-2 for an emitting region whose temperature is between 30K and 250 K, or < 1 10^16 cm^-2 if the O2 emitting region is primarily at a temperature of ~< 100 K. Because the Orion Bar is oriented nearly edge-on relative to our line of sight, the observed column density is enhanced by a factor estimated to be between 4 and 20 relative to the face-on value. Our upper limits imply that the face-on O2 column density is less than 4 10^15 cm^-2, a value that is below, and possibly well below, model predictions for gas with a density of 10^4 – 10^5 cm^-3 exposed to a far ultraviolet flux 10^4 times the local value, conditions inferred from previous observations of the Orion Bar. The discrepancy might be resolved if: (1) the adsorption energy of O atoms to ice is greater than 800 K; (2) the total face-on Av of the Bar is less than required for O2 to reach peak abundance; (3) the O2 emission arises within dense clumps with a small beam filling factor; or, (4) the face-on depth into the Bar where O2 reaches its peak abundance, which is density dependent, corresponds to a sky position different from that sampled by our Herschel beams.

The massive binary system Eta Carinae and the surrounding HII complex, the Carina Nebula, are potential particle acceleration sites from which very-high-energy (VHE; E > 100 GeV) \gamma-ray emission could be expected. This paper presents data collected during VHE \gamma-ray observations with the H.E.S.S. telescope array from 2004 to 2010, which cover a full orbit of Eta Carinae. In the 33.1-hour data set no hint of significant \gamma-ray emission from Eta Carinae has been found and an upper limit on the \gamma-ray flux of 7.7 x 10-13 ph cm-2 s-1 (99% confidence level) is derived above the energy threshold of 470 GeV. Together with the detection of high-energy (HE; 0.1 GeV > E > 100 GeV) \gamma-ray emission by the Fermi-LAT up to 100 GeV, and assuming a continuation of the average HE spectral index into the VHE domain, these results imply a cut-off in the \gamma-ray spectrum between the HE and VHE \gamma-ray range. This could be caused either by a cut-off in the accelerated particle distribution or by severe \gamma-\gamma\ absorption losses in the wind collision region. Furthermore, the search for extended \gamma-ray emission from the Carina Nebula resulted in an upper limit on the \gamma-ray flux of 4.2 x 10-12 ph cm-2 s-1 (99% confidence level). The derived upper limit of ~23 on the cosmic-ray enhancement factor is compared with results found for the old-age mixed-morphology supernova remnant W 28.

We have carried out survey observations toward the Galactic plane at l~38 deg in the 12CO and 13CO J=1-0 lines using the Nobeyama Radio Observatory 45-m telescope. A wide area (0.8 x 0.8 deg) was mapped with high spatial resolution (17″). The line of sight samples the gas in both the Sagittarius arm and the inter-arm regions. The present observations reveal how the structure and physical conditions vary across a spiral arm. We classify the molecular gas in the line of sight into two distinct components based on its appearance: the bright and compact B component and the fainter and diffuse (i.e., more extended) D component. The B component is predominantly seen at the spiral arm velocities, while the D component dominates at the inter-arm velocities and is also found at the spiral arm velocities. We introduce the brightness distribution function and the brightness distribution index (BDI, which indicates the dominance of the B component) in order to quantify the map’s appearance. The radial velocities of BDI peaks coincide with those of high 12CO J=3-2/12CO J=1-0 intensity ratio (i.e., warm gas) and H II regions, and tend to be offset from the line brightness peaks at lower velocities (i.e., presumably downstream side of the arm). Our observations reveal that the gas structure at small scales changes across a spiral arm: bright and spatially confined structures develop in a spiral arm, leading to star formation at downstream side, while extended emission dominates in the inter-arm region.

Recent observational results indicate that the functional shape of the spatially-resolved star formation-molecular gas density relation depends on the spatial scale considered. These results may indicate a fundamental role of sampling effects on scales that are typically only a few times larger than those of the largest molecular clouds. To investigate the impact of this effect, we construct simple models for the distribution of molecular clouds in a typical star-forming spiral galaxy, and, assuming a power-law relation between SFR and cloud mass, explore a range of input parameters. We confirm that the slope and the scatter of the simulated SFR-molecular gas surface density relation depend on the size of the sub-galactic region considered, due to stochastic sampling of the molecular cloud mass function, and the effect is larger for steeper relations between SFR and molecular gas. There is a general trend for all slope values to tend to ~unity for region sizes larger than 1-2 kpc, irrespective of the input SFR-cloud relation. The region size of 1-2 kpc corresponds to the area where the cloud mass function becomes fully sampled. We quantify the effects of selection biases in data tracing the SFR, either as thresholds (i.e., clouds smaller than a given mass value do not form stars) or backgrounds (e.g., diffuse emission unrelated to current star formation is counted towards the SFR). Apparently discordant observational results are brought into agreement via this simple model, and the comparison of our simulations with data for a few galaxies supports a steep (>1) power law index between SFR and molecular gas.

We have previously reported on chemical abundance trends with evolutionary state in the globular cluster NGC 6397 discovered in analyses of spectra taken with FLAMES at the VLT. Here, we reinvestigate the FLAMES-UVES sample of 18 stars, ranging from just above the turnoff point (TOP) to the red giant branch below the bump. Inspired by new calibrations of the infrared flux method, we adopt a set of hotter temperature scales. Chemical abundances are determined for six elements (Li, Mg, Ca, Ti, Cr, and Fe). Signatures of cluster-internal pollution are identified and corrected for in the analysis of Mg. On the modified temperature scales, evolutionary trends in the abundances of Mg and Fe are found to be significant at the 2{\sigma} and 3{\sigma} levels, respectively. The detailed evolution of abundances for all six elements agrees with theoretical isochrones, calculated with effects of atomic diffusion and a weak to moderately strong efficiency of turbulent mixing. The age of these models is compatible with the external determination from the white dwarf cooling sequence. We find that the abundance analysis cannot be reconciled with the strong turbulent-mixing efficiency inferred elsewhere for halo field stars. A weak mixing efficiency reproduces observations best, indicating a diffusion-corrected primordial lithium abundance of log {\epsilon}(Li) = 2.57 +- 0.10. At 1.2{\sigma}, this value agrees well with WMAP-calibrated Big-Bang nucleosynthesis predictions.

Deep, near-infrared JHK-maps were observed for 10 nearby, grand-design, spiral galaxies using HAWK-I/VLT to study the distribution of young stellar clusters in them and thereby determine whether strong spiral perturbations can influence star formation. Complete, magnitude-limited candidate lists of star-forming complexes were obtained by searching within the K-band maps. The properties of the complexes were derived from (H-K)-(J-H) diagrams including the identification of the youngest complexes (i.e. <7 Myr) and the estimation of their extinction. Young stellar clusters with ages <7 Myr have significant internal extinction in the range of Av=3-7m, while older ones typically have Av<1m. The cluster luminosity function (CLF) is well-fitted by a power law with an exponent of around -2 and displays no evidence of a high luminosity cut-off. The brightest cluster complexes in the disk reach luminosities of Mk = -15.5m or estimated masses of 10^6 Mo. At radii with a strong, two-armed spiral pattern, the star formation rate in the arms is higher by a factor of 2-5 than in the inter-arm regions. The CLF in the arms is also shifted towards brighter Mk by at least 0.4m. We also detect clusters with colors compatible with Large Magellanic Cloud intermediate age clusters and Milky Way globular clusters. The (J-K)-Mk diagram of several galaxies shows, for the brightest clusters, a clear separation between young clusters that are highly attenuated by dust and older ones with low extinction. The gap in the (J-K)-Mk diagrams implies that there has been a rapid expulsion of dust at an age around 7 Myr, possibly triggered by supernovae. Strong spiral perturbations concentrate the formation of clusters in the arm regions and shifts their CLF towards brighter magnitudes.

We present the first numerical simulations that self-consistently follow the formation of dense molecular clouds in colliding flows. Our calculations include a time-dependent model for the H2 and CO chemistry that runs alongside a detailed treatment of the dominant heating and cooling processes in the ISM. We adopt initial conditions characteristic of the warm neutral medium and study two different flow velocities – a slow flow with v = 6.8 km/s and a fast flow with v = 13.6 km/s. The clouds formed by the collision of these flows form stars, with star formation beginning after 16 Myr in the case of the slower flow, but after only 4.4 Myr in the case of the faster flow. In both flows, the formation of CO-dominated regions occurs only around 2 Myr before the onset of star formation. Prior to this, the clouds produce very little emission in the J = 1 -> 0 transition line of CO, and would probably not be identified as molecular clouds in observational surveys. In contrast, our models show that H2-dominated regions can form much earlier, with the timing depending on the details of the flow. In the case of the slow flow, small pockets of gas become fully molecular around 10 Myr before star formation begins, while in the fast flow, the first H2-dominated regions occur around 3 Myr before the first prestellar cores form. Our results are consistent with models of molecular cloud formation in which the clouds are dominated by “dark” molecular gas for a considerable proportion of their assembly history.

We present the statistical properties of molecular clumps in the Galactic center 50 km s$^{-1}$ molecular cloud (GCM-0.02-0.07) based on observations of the CS $J=1-0$ emission line with the Nobeyama Millimeter Array. In the cloud, 37 molecular clumps with local thermal equilibrium (LTE) masses of $2\times10^2-6\times10^3 M_\odot$ were identified by using the {\it clumpfind} algorithm. The velocity widths of the molecular clumps are about five-fold those of Galactic disk molecular clouds with the same radius. The virial-theorem masses are three-fold the LTE masses. The mass and size spectra can be described by power laws of $dN/dM\propto M^{-2.6\pm0.1}$ ($M\gtrsim 900M_\odot$) and $dN/dR\propto R^{-5.9\pm0.3}$ ($R\gtrsim 0.35$ pc), respectively. The statistical properties of the region interacting with the Sgr A East shell and those of the non-interacting part of the cloud are significantly different. The interaction probably makes the mass function steeper, from $dN/dM\propto M^{-2.0\pm0.1}$ in the non-interacting part to $dN/dM\propto M^{-4.0\pm0.2}$ in the interacting region. On the other hand, the interaction presumably truncates the size spectrum on the larger side of $R\sim 0.4$ pc.

In this paper, we investigate the structure of condensation fronts from warm diffuse gas to cold neutral medium (CNM) under the plane parallel geometry. The solutions have two parameters, the pressure of the CNM and the mass flux across the transition front, and their ranges are much wider than previously thought. First, we consider the pressure range where the three phases, the CNM, the unstable phase, and the warm neutral medium, can coexist in the pressure equilibrium. In a wide range of the mass flux, we find solutions connecting the CNM and the unstable phase. Moreover, we find solutions in larger pressure range where there is only one thermal equilibrium state or the CNM. These solutions can be realized in shock-compressed regions that are promising sites of molecular cloud formation. We also find remarkable properties in our solutions. Heat conduction becomes less important with increasing mass flux, and the thickness of the transition layer is characterized by the cooling length instead of the Field length.

In this paper, we investigate the structure of condensation fronts from warm diffuse gas to cold neutral medium (CNM) under the plane parallel geometry. The solutions have two parameters, the pressure of the cold neutral medium (CNM) and the mass flux across the transition front and their ranges are much wider than previous thought. First, we consider the pressure range where the three phases, the CNM, the unstable phase, and the warm neutral medium, can coexist in the pressure equilibrium. In a wide range of the mass flux, we find solutions connecting the CNM and unstable phase. Moreover, we find solutions in larger pressure range where there is only one thermal equilibrium state, or the CNM. These solutions can be realized in shock-compressed regions that is promising sites of the molecular cloud formation. We also find remarkable properties in our solutions. The heat conduction becomes less important with increasing mass flux, and the thickness of the transition layer is characterized by the cooling length instead of the Field length.

We discuss the relevance of directional detection experiments in the post-discovery era and propose a method to extract the local dark matter phase space distribution from directional data. The first feature of this method is a parameterization of the dark matter distribution function in terms of integrals of motion, which can be analytically extended to infer properties of the global distribution if certain equilibrium conditions hold. The second feature of our method is a decomposition of the distribution function in moments of a model independent basis, with minimal reliance on the ansatz for its functional form. We illustrate our method using the Via Lactea II N-body simulation as well as an analytical model for the dark matter halo. We conclude that O(1000) events are necessary to measure deviations from the Standard Halo Model and constrain or measure the presence of anisotropies.

We test two different methods of using near-infrared extinction to estimate distances to dark clouds in the first quadrant of the Galaxy using large near infrared (2MASS and UKIDSS) surveys. VLBI parallax measurements of masers around massive young stars provide the most direct and bias-free measurement of the distance to these dark clouds. We compare the extinction distance estimates to these maser parallax distances. We also compare these distances to kinematic distances, including recent re-calibrations of the Galactic rotation curve. The extinction distance methods agree with the maser parallax distances (within the errors) between 66% and 100% of the time (depending on method and input survey) and between 85% and 100% of the time outside of the crowded Galactic center. Although the sample size is small, extinction distance methods reproduce maser parallax distances better than kinematic distances; furthermore, extinction distance methods do not suffer from the kinematic distance ambiguity. This validation gives us confidence that these extinction methods may be extended to additional dark clouds where maser parallaxes are not available.

We present a new method to determine the age spread of resolved stellar populations in a starburst cluster. The method relies on a two-step process. In the first step, kinematic members of the cluster are identified based on multi-epoch astrometric monitoring. In the second step, a Bayesian analysis is carried out, comparing the observed photometric sequence of cluster members with sets of theoretical isochrones. When applying this methodology to optical and near-infrared high angular resolution Hubble Space Telescope (HST) and adaptive optics observations of the ~5 Myr old starburst cluster Westerlund 1 and ~2 Myr old starburst cluster NGC 3603 YC, we derive upper limits for the age spreads of 0.4 and 0.1 Myr, respectively. The results strongly suggest that star formation in these starburst clusters happened almost instantaneously.

We perform the largest currently available set of direct N-body calculations of young star cluster models to study the dynamical influence, especially through the ejections of the most massive star in the cluster, on the current relation between the maximum-stellar-mass and the star-cluster-mass. We vary several initial parameters such as the initial half-mass radius of the cluster, the initial binary fraction, and the degree of initial mass segregation. Two different pairing methods are used to construct massive binaries for more realistic initial conditions of massive binaries. We find that lower mass clusters (= 1000 Msun), no most-massive star escapes the cluster within 3 Myr regardless of the initial conditions if clusters have initial half-mass radii, r_0.5, >= 0.8 pc. However, a few of the initially smaller sized clusters (r_0.5 = 0.3 pc), which have a higher density, eject their most massive star within 3 Myr. If clusters form with a compact size and their massive stars are born in a binary system with a mass-ratio biased towards unity, the probability that the mass of the most massive star in the cluster changes due to the ejection of the initially most massive star can be as large as 20 per cent. Stellar collisions increase the maximum-stellar-mass in a large number of clusters when clusters are relatively dense (M_ecl >= 10^3 Msun and r_0.5 = 0.3 pc) and binary-rich. Overall, we conclude that dynamical effects hardly influence the observational maximum-stellar-mass — cluster mass relation.

We used the N-body code of Hernquist and Ostriker (1992) to build a dozen cuspy ({\gamma}\approx 1) triaxial models of stellar systems through dissipationless collapses of initially spherical distributions of 10^6 particles. We chose four sets of initial conditions that resulted in models morphologically resembling E2, E3, E4 and E5 galaxies, respectively. Within each set, three different seed numbers were selected for the random number generator used to create the initial conditions, so that the three models of each set are statistically equivalent. We checked the stability of our models using the values of their central densities and of their moments of inertia, which turned out to be very constant indeed. The changes of those values were all less than 3 per cent over one Hubble time and, moreover, we show that the most likely cause of those changes are relaxation effects in the numerical code. We computed the six Lyapunov exponents of nearly 5,000 orbits in each model in order to recognize regular, partially and fully chaotic orbits. All the models turned out to be highly chaotic, with less than 25 per cent of their orbits being regular. We conclude that it is quite possible to obtain cuspy triaxial stellar models that contain large fractions of chaotic orbits and are highly stable. The difficulty to build such models with the method of Schwarzschild (1979) should be attributed to the method itself and not to physical causes.

We present a comprehensive study of the periodic variations observed in OGLE I-band light curves of SMC Be/X-ray binaries, discovering new optical periodicities in 9 systems. We find that these periodicities derive from a number of mechanisms, notably disturbance of the decretion disk on the orbital period of the system, and aliased non-radial pulsations. We develop metrics that allow these mechanisms to be distinguished on the basis of the shape of the folded optical light curve, and use these metrics to categorise the periodicities present in \sim 50 SMC binary systems. We conclude that extreme care must be taken in the interpretation of the OGLE light curves since only around 30% of the periodicities present can be unambiguously attributed to orbital periods.

The disk around the Herbig Ae star HD\,169142 was imaged and resolved at 18.8 and 24.5\,$\mu$m using Subaru/COMICS. We interpret the observations using a 2D radiative transfer model and find evidence for the presence of a large gap. The MIR images trace dust that emits at the onset of the strong rise in the spectral energy distribution (SED) at 20\,$\mu$m, therefore are very sensitive to the location and characteristics of the inner wall of the outer disk and its dust. We determine the location of the wall to be 23$^{+3}_{-5}$\,AU from the star. An extra component of hot dust must exist close to the star. We find that a hydrostatic optically thick inner disk does not produce enough flux in the NIR and an optically thin geometrically thick component is our solution to fit the SED. Considering the recent findings of gaps and holes in a number of Herbig Ae/Be group I disks, we suggest that such disk structures may be common in group I sources. Classification as group I should be considered a support for classification as a transitional disk, though improved imaging surveys are needed to support this speculation.

Blue (metal-poor) globular clusters are observed to have half-light radii that are ~20% larger than their red (metal-rich) counterparts. The origin of this enhancement is not clear and differences in either the luminosity function or in the actual size of the clusters have been proposed. I analyze a set of dynamically self-consistent Monte Carlo globular cluster simulations to determine the origin of this enhancement. I find that my simulated blue clusters have larger half-light radii due to differences in the luminosity functions of metal-poor and metal-rich stars. I find that the blue clusters can also be physically larger, but only if they have a substantial number of black holes heating their central regions. In this case the difference between half-light radii is significantly larger than observed. I conclude that the observed difference in half-light radii between red and blue globular clusters is due to differences in their luminosity functions and that half-light radius is not a reliable proxy for cluster size.

We present high-resolution I-band imaging of the core of the globular cluster M15 obtained at the 2.5 m Nordic Optical Telescope with FastCam, a low readout noise L3CCD based instrument. Short exposure times (30 ms) were used to record 200000 images (512 x 512 pixels each) over a period of 2 hours 43 min. The lucky imaging technique was then applied to generate a final image of the cluster centre with FWHM ~ 0″.1 and 13″ x 13″ FoV. We obtained a catalogue of objects in this region with a limiting magnitude of I=19.5. I-band photometry and astrometry are reported for 1181 stars. This is the deepest I-band observation of the M15 core at this spatial resolution. Simulations show that crowding is limiting the completeness of the catalogue. At shorter wavelengths, a similar number of objects has been reported using HST/WFPC observations of the same field. The cross-match with the available HST catalogues allowed us to produce colour-magnitude diagrams where we identify new Blue Straggler star candidates and previously known stars of this class.

The purpose of this study is to analyze the dynamical role of a radiation field on the growth rate of the unstable Kelvin – Helmholtz (KH) perturbations. As a first step toward this purpose, the analyze is done in a general way, irrespective of applying the model to a specific astronomical system. The transition zone between the two layers of the fluid is ignored. Then, we perform a linear analysis and by imposing suitable boundary conditions and considering a radiation field, we obtain appropriate dispersion relation. Unstable modes are studied by solving the dispersion equation numerically, and then growth rates of them are obtained. By analyzing our dispersion relation, we show that for a wide range of the input parameters, the radiation field has a destabilizing effect on KH instability. In eruptions of the galaxies or supermassive stars, the radiation field is dynamically important and because of the enhanced KH growth rates in the presence of the radiation; these eruptions can inject more momentum and energy into their environment and excite more turbulent motions.

With the Infrared Camera on board AKARI, we carried out near-infrared (2.5-5.0 micron) spectroscopy of the central kiloparsec region of the barred spiral galaxy, NGC1097, categorized as Seyfert 1 with a circumnuclear starburst ring. Our observations mapped the area of ~50″*10″ with the resolution of ~5″, covering about a half of the ring and the galactic center. As a result, we spatially resolve the starburst ring in the polycyclic aromatic hydrocarbon 3.3 micron, the aliphatic hydrocarbon 3.4-3.6 micron features, and the hydrogen Br alpha 4.05 micron emission. They exhibit spatial distributions significantly different from each other, indicating that the environments vary considerably around the ring. In particular, the aliphatic features are enhanced near the bar connecting the ring with the nucleus, where the structure of hydrocarbon grains seems to be relatively disordered. Near the center, the continuum emission and the CO/SiO absorption features are strong, which indicates that the environments inside the ring are dominated by old stellar populations. The near-infrared spectra do not show any evidence for the presence of nuclear activity.

Since globular clusters (GCs) are old, low-N systems their dynamics is widely believed to be fully dominated by collisional two-body processes, and their surface brightness profiles are fit by King models. However, for many GCs, especially those with HST-resolved central regions, and `extra-tidal’ features, King models provide poor fits. We suggest that this is partly because collisionless dynamics is also important and contribute to shaping the cluster properties. We show using time-scale and length-scale arguments that except for the very centers of clusters, collisionless dynamics should be more important than collisional. We then fit 38 GCs analyzed by Noyola and Gebhardt (2006) with (collisional) King and (collisionless) DARKexp models over the full available radial range, and find that the latter provide a better fit to 29 GCs; for six of these the fit is at least ~5x better in term of rms. DARKexp models are theoretically derived maximum entropy equilibrium states of self-gravitating collisionless systems and have already been shown to fit the results of dark matter N-body simulations. (We do not attempt fits with ad hoc fitting functions.)

In the present article, we present a galactic gravitational model of three degrees of freedom, in order to investigate and reveal the behavior of orbits in a binary quasar system. The two quasars are hosted in a pair of interacting disk galaxies. We study in detail the regular or chaotic character of motion, in two different cases: the time independent model in both 2D and 3D dynamical systems and the evolving 3D model. Our numerical calculations, indicate that a large fraction of orbits in the 2D system are chaotic in the time independent case. A careful analysis suggest that several Lindblad resonances are also responsible for the chaotic motion of stars in both hosts galaxies. In the time dependent system, we follow the evolution of 3D orbits in our dynamical model, as the two interacting host galaxies develop dense and massive quasars in their cores, by mass transportation from the disks to their nuclei. In this interesting case, there are orbits that change their orbital character from regular to chaotic and vise versa and also orbits that maintain their character during the galactic evolution. These results strongly indicate that the ordered or chaotic nature of 3D orbits depends, not only in the galactic interaction but also in the presence of the quasars in the galactic cores of the host galaxies. The outcomes derived from our dynamical model are compared with observational data. Some theoretical arguments to support the numerically derived outcomes are presented, both in 2D and 3D systems, while a comparison with similar earlier work is also made.

We use distant Blue Horizontal Branch stars with Galactocentric distances 16 < r/kpc < 48 as kinematic tracers of the Milky Way dark halo. We model the tracer density as an oblate, power-law embedded within a spherical power-law potential. Using a distribution function method, we estimate the overall power-law potential and the velocity anisotropy of the halo tracers. We measure the slope of the potential to be gamma ~ 0.4 and the overall mass within 50 kpc is ~ 4 x 10^11 M_sol. The tracer velocity anisotropy is radially biased with beta ~ 0.5, which is in good agreement with local solar neighbourhood studies. Our results provide an accurate outer circular velocity profile for the Milky Way and suggest a relatively high concentration dark matter halo (c_vir ~ 20).

We examine the \alpha-element abundance ratio, [\alpha/Fe], of 5620 stars, observed by the Sloan Extension for Galactic Understanding and Exploration survey in the region 6 kpc < R < 16 kpc, 0.15 kpc < |Z| < 1.5 kpc, as a function of Galactocentric radius R and distance from the Galactic plane |Z|. Our results show that the high-\alpha\ thick disk population has a short scale length (L_thick ~ 1.8 kpc) compared to the low-\alpha population, which is typically associated with the thin disk. We find that the fraction of high-\alpha\ stars in the inner disk increases at large |Z|, and that high-\alpha\ stars lag in rotation compared to low-\alpha\ stars. In contrast, the fraction of high-\alpha\ stars in the outer disk is low at all |Z|, and high- and low-\alpha\ stars have similar rotational velocities up to 1.5 kpc from the plane. We interpret these results to indicate that different processes were responsible for the high-\alpha\ populations in the inner and outer disk. The high-\alpha\ population in the inner disk has a short scale length and large scale height, consistent with a scenario in which the thick disk forms during an early gas-rich accretion phase. Stars far from the plane in the outer disk may have reached their current locations through heating by minor mergers. The lack of high-\alpha\ stars at large R and |Z| also places strict constraints on the strength of radial migration via transient spiral structure.

It has been known for a long time that the satellite galaxies of the Milky Way (MW) show a significant amount of phase-space correlation, they are distributed in a highly inclined Disc of Satellites (DoS). We have extended the previous studies on the DoS by analysing for the first time the orientations of streams of stars and gas, and the distributions of globular clusters within the halo of the MW. It is shown that the spatial distribution of MW globular clusters classified as young halo clusters (YH GC) is very similar to the DoS, while 7 of the 14 analysed streams align with the DoS. The probability to find the observed clustering of streams is only 0.3 per cent when assuming isotropy. The MW thus is surrounded by a vast polar structure (VPOS) of subsystems (satellite galaxies, globular clusters and streams), spreading from Galactocentric distances as small as 10 kpc out to 250 kpc. These findings demonstrate that a near-isotropic infall of cosmological sub-structure components onto the MW is essentially ruled out because a large number of infalling objects would have had to be highly correlated, to a degree not natural for dark matter sub-structures. The majority of satellites, streams and YH GCs had to be formed as a correlated population. This is possible in tidal tails consisting of material expelled from interacting galaxies. We discuss the tidal scenario for the formation of the VPOS, including successes and possible challenges. The potential consequences of the MW satellites being tidal dwarf galaxies are severe. If all the satellite galaxies and YH GCs have been formed in an encounter between the young MW and another gas-rich galaxy about 10-11 Gyr ago, then the MW does not have any luminous dark-matter substructures and the missing satellites problem becomes a catastrophic failure of the standard cosmological model.

These lectures attempt to expose the most important ideas, which have been proposed to explain the formation of stars with particular emphasis on the formation of brown dwarfs and low-mass stars. We first describe the important physical processes which trigger the collapse of a self-gravitating piece of fluid and regulate the star formation rate in molecular clouds. Then we review the various theories which have been proposed along the years to explain the origin of the stellar initial mass function paying particular attention to four models, namely the competitive accretion and the theories based respectively on stopped accretion, MHD shocks and turbulent dispersion. As it is yet unsettled whether the brown dwarfs form as low-mass stars, we present the theory of brown dwarfs based on disk fragmentation stressing all the uncertainties due to the radiative feedback and magnetic field. Finally, we describe the results of large scale simulations performed to explain the collapse and fragmentation of molecular clouds.

We have selelcted 556 Red Horizontal Branch (RHB) stars along the streams of the Sagittarius dwarf galaxy (Sgr) from SDSS DR7 spectroscopic data using a theoretical model. The metallicity and \alpha-elements distributions are investigated for stars in the Sgr streams and for Galactic stars at the same locations. We find that the Sgr stars have two peaks in the metallicity distribution while the Galactic stars have a more prominent metal-poor peak. Meanwhile, [\alpha/Fe] ratios of the Sgr stars are lower than those of the Galactic stars. Among the Sgr stars, we find a difference in the metallicity distribution between the leading and trailing arms of the Sgr tidal tails. The metallicity and [\alpha/Fe] distribution of the leading arm is similar to that of the Galaxy. The trailing arm is composed mainly of a metal rich component and [\alpha/Fe] is obviously lower than that of the Galactic stars. The metallicity gradient is -(1.8 \pm 0.3)\times10^{-3} dex degree^{-1} in the first wrap of the trailing arm and -(1.5 \pm 0.4)\times10^{-3} dex degree^{-1} in the first wrap of the leading arm. No significant gradient exists along the second wraps of the leading or trailing arms. It seems that the Sgr dwarf galaxy initially lost the metal poor component in the second wrap (older) arms due to the tidal force of our Galaxy and then the metal rich component is disrupted in the first wrap (younger) arms. Finally, we found that the velocity dispersion of the trailing arm from 88^\circ<\Lambda_\odot<112^\circ is \sigma = 9.808 \pm 1.0 km s^{-1}, which is consistent with previous work in the literature.

In this paper we examine the effect of the formation and evolution of the disk galaxy on the distribution of dark halo matter. We have made simulations of isolated dark matter (DM) halo and two component (DM + baryons). N-body technique was used for stellar and DM particles and TVD MUSCL scheme for gas-dynamic simulations. The simulations include the processes of star formation, stellar feedback, heating and cooling of the interstellar medium. The results of numerical experiments with high spatial resolution let us to conclude in two main findings. First, accounting of star formation and supernova feedback resolves the so-called problem of cusp in distribution of dark matter predicted by cosmological simulations. Second, the interaction of dark matter with dynamic substructures of stellar and gaseous galactic disk (e.g., spiral waves, bar) has an impact on the shape of the dark halo. In particular, the in-plane distribution of dark matter is more symmetric in runs, where the baryonic component was taken into account.

We report the initial results from an ongoing multi-year spectroscopic survey of novae in M33. The survey resulted in the spectroscopic classification of six novae (M33N 2006-09a, 2007-09a, 2009-01a, 2010-10a, 2010-11a, and 2011-12a) and a determination of rates of decline (t_2 times) for four of them (2006-09a, 2007-09a, 2009-01a, and 2010-10a). When these data are combined with existing spectroscopic data for two additional M33 novae (2003-09a and 2008-02a) we find that 5 of the 8 novae with available spectroscopic class appear to be members of either the He/N or Fe IIb (hybrid) classes, with only two clear members of the Fe II spectroscopic class. This initial finding is very different from what would be expected based on the results for M31 and the Galaxy where Fe II novae dominate, and the He/N and Fe IIb classes together make up only ~20% of the total. It is plausible that the increased fraction of He/N and Fe IIb novae observed in M33 thus far may be the result of the younger stellar population that dominates this galaxy, which is expected to produce novae that harbor generally more massive white dwarfs than those typically associated with novae in M31 or the Milky Way.

We report the initial results from an ongoing multi-year spectroscopic survey of novae in M33. The survey resulted in the spectroscopic classification of six novae (M33N 2006-09a, 2007-09a, 2009-01a, 2010-10a, 2010-11a, and 2011-12a) and a determination of rates of decline (t_2 times) for four of them (2006-09a, 2007-09a, 2009-01a, and 2010-10a). When these data are combined with existing spectroscopic data for two additional M33 novae (2003-09a and 2008-02a) we find that 5 of the 8 novae with available spectroscopic class appear to be members of either the He/N or Fe IIb (hybrid) classes, with only two clear members of the Fe II spectroscopic class. This initial finding is very different from what would be expected based on the results for M31 and the Galaxy where Fe II novae dominate, and the He/N and Fe IIb classes together make up only ~20% of the total. It is plausible that the increased fraction of He/N and Fe IIb novae observed in M33 thus far may be the result of the younger stellar population that dominates this galaxy, which is expected to produce novae that harbor generally more massive white dwarfs than those typically associated with novae in M31 or the Milky Way.

We report measurement of trigonometric parallax of IRAS 05168+3634 with VERA. The parallax is 0.532 +/- 0.053 mas, corresponding to a distance of 1.88+0.21/-0.17 kpc. This result is significantly smaller than the previous distance estimate of 6 kpc based on kinematic distance. This drastic change in the source distance revises not only physical parameters of IRAS 05168+3634, but also its location of the source, placing it in the Perseus arm rather than the Outer arm. We also measure proper motions of the source. A combination of the distance and the proper motions with systemic velocity yields rotation velocity ({\Theta}) of 227+9/-11 km s-1 at the source, assuming {\Theta}0 = 240 km s-1. Our result combined with previous VLBI results for six sources in the Perseus arm indicates that the sources rotate systematically slower than the Galactic rotation velocity at the LSR. In fact, we show observed disk peculiar motions averaged over the seven sources in the Perseus arm as (Umean, Vmean) = (11 +/- 3, -17 +/- 3) km s-1, indicating that these seven sources are systematically moving toward the Galactic center, and lag behind the Galactic rotation.

The blue compact dwarf galaxy I Zw 18 is one of the most metal poor systems known in the local Universe (12 + log(O/H) $=$ 7.17). In this work we study I Zw 18 using data from {\it Spitzer}, {\it Herschel Space Telescope} and IRAM Plateau de Bure Interferometer. Our data set includes the most sensitive maps of I Zw 18, to date, in both, the far infrared and the CO $J=1\rightarrow0$ transition. We use dust emission models to derive a dust mass upper limit of only M$_{dust}\leq1.1\times10^4$ M$_{\odot}$ ($3\sigma$ limit). This upper limit is driven by the non-detection at 160 $\mu$m, and it is a factor of 4-10 times smaller than previous estimates (depending upon the model used). We also estimate an upper limit to the total dust-to-gas mass ratio of M$_{Dust}$/M$_{gas}\leq5.0\times10^{-5}$. If a linear correlation between the dust-to-gas mass ratio and metallicity (measure as O/H) were to hold, we would expect a ratio of 3.9$\times10^{-4}$. We also show that the infrared SED is similar to that of starbursting systems.

Observations towards IRC +10216 of CS, SiO, SiS, NaCl, KCl, AlCl, AlF, and NaCN have been carried out with the IRAM 30-m telescope in the 80-357.5 GHz frequency range. A large number of rotational transitions covering a wide range of energy levels, including highly excited vibrational states, are detected in emission and serve to trace different regions of the envelope. Radiative transfer calculations based on the LVG formalism have been performed to derive molecular abundances from the innermost out to the outer layers. The excitation calculations include infrared pumping to excited vibrational states and inelastic collisions, for which up-to-date rate coefficients for rotational and, in some cases, ro-vibrational transitions are used. We find that in the inner layers CS, SiO, and SiS have abundances relative to H$_2$ of 4e-6, 1.8e-7, and 3e-6, respectively, and that CS and SiS have significant lower abundances in the outer envelope, which implies that they actively contribute to the formation of dust. Moreover, in the inner layers, the amount of sulfur and silicon in gas phase molecules is only 27 % for S and 5.6 % for Si, implying that these elements have already condensed onto grains, most likely in the form of MgS and SiC. Metal-bearing molecules lock up a relatively small fraction of metals, although our results indicate that NaCl, KCl, AlCl, AlF, and NaCN, despite their refractory character, are not significantly depleted in the cold outer layers. In these regions a few percent of the metals Na, K, and Al survive in the gas phase, either in atomic or molecular form, and are therefore available to participate in the gas phase chemistry in the outer envelope.

We study the regular or chaotic nature of motion in a disk galaxy with a dense nucleus and an asymmetric dark halo. Two cases, the 2D model and the 3D model, are investigated. In the 2D model, a considerable fraction of the phase plane is covered by chaotic orbits. Two factors seem to be responsible for the chaotic motion: (i) the dense nucleus and (ii) the asymmetries in the dark halo. Our numerical experiments suggest, that there are several chaotic components on the Poincare phase plane. Different chaotic components are induced by the asymmetries in the halo. Each chaotic component seems to have a different value of the Lyapunov Characteristic Exponent, for small values of the asymmetry parameter and a unique LCE for larger values of the asymmetry parameter. A comparison of the present results with outcomes from previous work is also presented.

Context. In photo-dissociation regions (PDRs), Polycyclic Aromatic Hydrocarbons (PAHs) could be produced by evaporation of Very Small Grains (VSGs) by the impinging UV radiation field from a nearby star. Aims. We investigate quantitatively the transition zone between evaporating Very Small Grains (eVSGs) and PAHs in several PDRs. Methods. We study the relative contribution of PAHs and eVSGs to the mid-IR emission in a wide range of excitation conditions. We fit the observed mid-IR emission of PDRs by using a set of template band emission spectra of PAHs, eVSGs and gas lines. The fitting tool PAHTAT (PAH Toulouse Astronomical Templates) is made available to the community as an IDL routine. From the results of the fit, we derive the fraction of carbon f_eVSG locked in eVSGs and compare it to the intensity of the local UV radiation field. Results. We show a clear decrease of f_eVSG with increasing intensity of the local UV radiation field, which supports the scenario of photo-destruction of eVSGs. Conversely, this dependence can be used to quantify the intensity of the UV radiation field for different PDRs, including non resolved ones. Conclusions. PAHTAT can be used to trace the intensity of the local UV radiation field in regions where eVSGs evaporate, which correspond to relatively dense (nH = [100, 10^5 ] cm-3) and UV irradiated PDRs (G0 = [100, 5x10^4]) where H2 emits in rotational lines.

Context. In photo-dissociation regions (PDRs), Polycyclic Aromatic Hydrocarbons (PAHs) could be produced by evaporation of Very Small Grains (VSGs) by the impinging UV radiation field from a nearby star. Aims. We investigate quantitatively the transition zone between evaporating Very Small Grains (eVSGs) and PAHs in several PDRs. Methods. We study the relative contribution of PAHs and eVSGs to the mid-IR emission in a wide range of excitation conditions. We fit the observed mid-IR emission of PDRs by using a set of template band emission spectra of PAHs, eVSGs and gas lines. The fitting tool PAHTAT (PAH Toulouse Astronomical Templates) is made available to the community as an IDL routine. From the results of the fit, we derive the fraction of carbon feVSG locked in eVSGs and compare it to the intensity of the local UV radiation field. Results. We show a clear decrease of feVSG with increasing intensity of the local UV radiation field, which supports the scenario of photo-destruction of eVSGs. Conversely, this dependence can be used to quantify the intensity of the UV radiation field for different PDRs, including non resolved ones. Conclusions. PAHTAT can be used to trace the intensity of the local UV radiation field in regions where eVSGs evaporate, which correspond to relatively dense (nH = [100, 105 ] cm-3) and UV irradiated PDRs (G0 = [100, 5x10^4]) where H2 emits in rotational lines.

We investigate the possibility that the observed rotation of galaxies can be accounted for by invoking a massive baryonic disc with no need for non-baryonic dark matter or a massive halo. There are 5 primary reasons for suggesting this: 1. there are well known disc surface mass density distributions that naturally produce the observed rotation curves of galaxies. 2. there are a number of rotation curve `puzzles’ that cannot be explained by a massive dark matter halo i.e. the success of maximum disc fitting, HI gas scaling to the observed rotation, the disc/halo conspiracy and the interpretation of the Tully-Fisher relation. 3. recent 21cm observations show an almost constant HI surface density and a distinct `cut-off’ or edge to galactic discs. We explain this constant surface density in terms of either an optical depth effect or the onset of molecular gas formation and hence the possibility of considerably more gas existing in galaxies. We suggest that the HI cut-off does indeed mark the edge of the galactic disc. 4. there have been an increasing number of recent observations that imply that X_CO may be ten times higher in the outer Galaxy. This `dark’ gas may provide adequate mass to account for galaxy rotation. 5. we show that the additional baryonic mass required to account for the rotation of galaxies is just that required to reconcile observed baryons with those predicted by big bang nucleosynthesis. Mestel discs can be used to straight forwardly explain the scaling laws of galaxies, particularly the observed relation between rotation velocity and radius and the oft used Tully-Fisher relation. We discuss observations of the baryonic content of galactic discs and where sufficient `hidden’ baryons might be found to account for the rotation.