We have developed our original made-to-measure (M2M) algorithm, PRIMAL, with the aim of modelling the Galactic disc from upcoming Gaia data. From a Milky Way like N-body disc galaxy simulation, we have created mock Gaia data using M0III stars as tracers, taking into account extinction and the expected Gaia errors. In PRIMAL, observables calculated from the N-body model are compared with the target stars, at the position of the target stars. Using PRIMAL, the masses of the N-body model particles are changed to reproduce the target mock data, and the gravitational potential is automatically adjusted by the changing mass of the model particles. We have also adopted a new resampling scheme for the model particles to keep the mass resolution of the N-body model relatively constant. We have applied PRIMAL to this mock Gaia data and we show that PRIMAL can recover the structure and kinematics of a Milky Way like barred spiral disc, along with the apparent bar structure and pattern speed of the bar despite the galactic extinction and the observational errors.
The measurement of Zeeman splitting in spectral lines---both in emission and absorption---can provide direct estimates of the magnetic field strength and direction in atomic and molecular clouds, both in our own Milky Way and in external galaxies. This method will probe the magnetic field in the warm and cold neutral components of the interstellar medium, providing a complement to the extensive SKA Faraday studies planning to probe the field in the ionized components.
We estimated the scale-length of the thin disc with the J and W1 magnitudes of the most probable Red Clump (RC) stars in the Galactic plane, $-0\overset{^\circ}.5 \leq b \leq +0\overset{^\circ}.5$, in 19 equal sized fields with consecutive Galactic longitudes which cover the interval $90^\circ \leq l \leq 270^\circ$. Our results are constrained with respect to the solar space density ($D^*=5.95$), which indicates that the radial variation of the density is lower for higher Galactocentric distances. The scale-length of the thin disc is 2 kpc for the fields in the Galactic anticentre direction or close to this direction, while it decreases continuously in the second and third quadrants reaching to a lower limit of $h$ = 1.6 kpc at the Galactic longitudes $l$ = 90$^\circ$ and $l$ = 270$^\circ$. The distribution of the scale-length in 19 fields is consistent with the predictions from the Galaxia model and its variation with longitude is probably due to the inhomogeneity structure of the disc caused by the accreted material or other features such as warp and flare.
Gas outflows are believed to play a pivotal role in shaping galaxies, as they regulate both star formation and black hole growth. Despite their ubiquitous presence, the origin and the acceleration mechanism of such powerful and extended winds is not yet understood. Direct observations of the cold gas component in objects with detected outflows at other wavelengths are needed to assess the impact of the outflow on the host galaxy interstellar medium (ISM). We observed with the Plateau de Bure Interferometer an obscured quasar at z~1.5, XID2028, for which the presence of an ionised outflow has been unambiguously signalled by NIR spectroscopy. The detection of CO(3-2) emission in this source allows us to infer the molecular gas content and compare it to the ISM mass derived from the dust emission. We then analyze the results in the context of recent insights on scaling relations, which describe the gas content of the overall population of star-forming galaxies at a similar redshifts. The Star formation efficiency (~100) and gas mass (M_gas=2.1-9.5x10^{10} M_sun) inferred from the CO(3-2) line depend on the underlying assumptions on the excitation of the transition and the CO-to-H2 conversion factor. However, the combination of this information and the ISM mass estimated from the dust mass suggests that the ISM/gas content of XID2028 is significantly lower than expected for its observed M$_\star$, sSFR and redshift, based on the most up-to-date calibrations (with gas fraction <20% and depletion time scale <340 Myr). Overall, the constraints we obtain from the far infrared and millimeter data suggest that we are observing QSO feedback able to remove the gas from the host
High-precision astrometry throughout the Local Group is a unique capability of the Hubble Space Telescope (HST), with potential for transformative science, including constraining the nature of dark matter, probing the epoch of reionization, and understanding key physics of galaxy evolution. While Gaia will provide unparalleled astrometric precision for bright stars in the inner halo of the Milky Way, HST is the only current mission capable of measuring accurate proper motions for systems at greater distances (> 80 kpc), which represents the vast majority of galaxies in the Local Group. The next generation of proper-motion measurements will require long time baselines, spanning many years to decades and possibly multiple telescopes, combining HST with the James Webb Space Telescope (JWST) or the Wide-Field Infrared Survey Telescope (WFIRST). However, the current HST allocation process is not conducive to such multi-cycle/multi-mission science, which will bear fruit primarily over many years. We propose an HST astrometry initiative to enable long-time-baseline, multi-mission science, which we suggest could be used to provide comprehensive kinematic measurements of all dwarf galaxies and high surface-density stellar streams in the Local Group with HST's Advanced Camera for Surveys (ACS) or Wide Field Camera 3 (WFC3). Such an initiative not only would produce forefront scientific results within the next 5 years of HST's life, but also would serve as a critical anchor point for future missions to obtain unprecedented astrometric accuracy, ensuring that HST leaves a unique and lasting legacy for decades to come.
Using cosmological hydrodynamics simulations we study the angular momentum content of the simulated galaxies in relation with their morphological type. We found that not only the angular momentum of the disk component follow the expected theoretical relation, Mo, Mao & Whiye (1998), but also the spheroidal one, with a gap due to its lost of angular momentum, in agreement with Fall & Romanowsky (2013),. We also found that the galaxy size can plot in one general relation, despite the morphological type, as found by Kravtsov (2013).
We present spectra of the nuclear regions of 50 nearby (D = 1 - 92 Mpc, median = 20 Mpc) galaxies of morphological types E to Sm. The spectra, obtained with the Gemini Near-IR Spectrograph on the Gemini North telescope, cover a wavelength range of approximately 0.85-2.5 microns at R~1300--1800. There is evidence that most of the galaxies host an active galactic nucleus (AGN), but the range of AGN luminosities (log (L2-10 keV [erg/s]) = 37.0-43.2) in the sample means that the spectra display a wide variety of features. Some nuclei, especially the Seyferts, exhibit a rich emission-line spectrum. Other objects, in particular the type 2 Low Ionisation Nuclear Emission Region galaxies, show just a few, weak emission lines, allowing a detailed view of the underlying stellar population. These spectra display numerous absorption features sensitive to the stellar initial mass function, as well as molecular bands arising in cool stars, and many other atomic absorption lines. We compare the spectra of subsets of galaxies known to be characterised by intermediate-age and old stellar populations, and find clear differences in their absorption lines and continuum shapes. We also examine the effect of atmospheric water vapor on the signal-to-noise ratio achieved in regions between the conventional NIR atmospheric windows, of potential interest to those planning observations of redshifted emission lines or other features affected by telluric H2O. Further exploitation of this data set is in progress, and the reduced spectra and data reduction tools are made available to the community.
We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $\lambda^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile.
We present the results of a targeted search for intervening HI absorption in six nearby, gas-rich galaxies using the Australia Telescope Compact Array (ATCA). The sightlines searched have impact parameters of 10-20 kpc. By targeting nearby galaxies we are also able to map their HI emission, allowing us to directly relate the absorption-line detection rate to the extended HI distribution. The continuum sightlines intersect the HI disk in four of the six galaxies, but no intervening absorption was detected. Of these four galaxies, we find that three of the non-detections are the result of the background source being too faint. In the fourth case we find that the ratio of the spin temperature to the covering factor ($T_{\mathrm{S}}/f$) must be much higher than expected ($\gtrsim$5700 K) in order to explain the non-detection. We discuss how the structure of the background continuum sources may have affected the detection rate of HI absorption in our sample, and the possible implications for future surveys. Future work including an expanded sample, and VLBI observations, would allow us to better investigate the expected detection rate, and influence of background source structure, on the results of future surveys.
We utilize magnetohydrodynamic (MHD) simulations to develop a numerical model for GMC-GMC collisions between nearly magnetically critical clouds. The goal is to determine if, and under what circumstances, cloud collisions can cause pre-existing magnetically subcritical clumps to become supercritical and undergo gravitational collapse. We first develop and implement new photodissociation region (PDR) based heating and cooling functions that span the atomic to molecular transition, creating a multiphase ISM and allowing modeling of non-equilibrium temperature structures. Then in 2D and with ideal MHD, we explore a wide parameter space of magnetic field strength, magnetic field geometry, collision velocity, and impact parameter, and compare isolated versus colliding clouds. We find factors of ~2-3 increase in mean clump density from typical collisions, with strong dependence on collision velocity and magnetic field strength, but ultimately limited by flux-freezing in 2D geometries. For geometries enabling flow along magnetic field lines, greater degrees of collapse are seen. We discuss observational diagnostics of cloud collisions, focussing on 13CO(J=2-1), 13CO(J=3-2), and 12CO(J=8-7) integrated intensity maps and spectra, which we synthesize from our simulation outputs. We find the ratio of J=8-7 to lower-J emission is a powerful diagnostic probe of GMC collisions.
We conducted a continuous survey of infrared and visual images of 18020 2MFGC galaxies which were selected on an automatic basis from 1.64 mln extended objects of the 2MASS XSC catalog based on the ratio of the infrared axes a/b>=3. This work aims to exclude "false" objects from the list of flat galaxies. Having observed more than 80 thousand images in different filters, we were able to detect 1512 such objects (8.4% of the total number). We found 23 galaxies duplicated in 2MASS, which have two 2MFGC numbers correspondingly, and three flat galaxies which are not included in other catalogs and are located close to three "false" galaxies. Galaxies with magnitudes fainter than K_s =13 mag compose the main part of the excluded objects. They show small angular sizes, low surface brightnesses and concentration ratios. The results of the work in the form of the 2MFGC table with notes are given in the astronomical databases VizieR, NED, HyperLeda.
The motivation for this paper is to review the limits set on the MACHO content of the Galactic halo by microlensing experiments in the direction of the Large Magellanic Cloud. This has been prompted by recent measurements of the Galactic rotation curve, which suggest that the limits have been biassed by the assumption of an over-massive halo. The paper first discusses the security of the detection efficiency calculations which are central to deriving the MACHO content of the Galactic halo. It then sets out to compare the rotation curves from various halo models with recent observations, with a view to establishing what limits can be put on an all-MACHO halo. The main thrust of the paper is to investigate whether lighter halo models which are consistent with microlensing by an all-MACHO halo are also consistent with recent measures of the Galactic rotation curve. In this case the population of bodies discovered by the MACHO collaboration would make up the entire dark matter content of the Galactic halo. The main result of this paper is that it is easy to find low mass halo models consistent with the observed Galactic rotation curve, which also imply an optical depth to microlensing similar to that found by the MACHO collaboration. This means that all-MACHO halos cannot be ruled out on the basis of their observations. In conclusion, limits placed on the MACHO content of the Galactic halo from microlensing surveys in the Magellanic Clouds are inconsistent and model dependent, and do not provide a secure basis for rejecting an all-MACHO halo.
It is common to attribute a flat rotation curve to our Galaxy. However Galazutdinov et al. (2015) in a recent paper have obtained a Keplerian rotation curve for outer parts of the Galaxy. They have calculated the distances from equivalent widths of interstellar CaII lines. The radial velocity was also measured on the interstellar CaII absorption line. We explain this unexpected result assuming elliptical orbits of some objects in our Galaxy. The application of formulas derived with the assumption of circular orbits to elliptical ones mimics the flat rotation curve.
Tidal interactions of galaxies in galaxy clusters have been proposed as one potential explanation of the morphology-density relation at low masses. The efficiency of tidal transformation is expected to depend strongly on the orbit of a galaxy within the cluster halo. The orbit determines both the strength of the cluster's global tidal field and the probability of encounters with other cluster members. Additionally, we investigate the effect of an inclination between the disk of the infalling galaxy and its orbital plane. We also compare our results to observational data. We find that galaxies that entered a cluster from the outskirts are unlikely to be significantly transformed. Closer to the cluster centre, tidal interactions are a more efficient mechanism for producing harassed galaxies. But the inclination of the disk can reduce the mass loss significantly, yet it amplifies the thickening.
Cheung (2007) identified a sample of 100 candidate X-shaped radio galaxies
using the NRAO FIRST survey; these are small-axial-ratio extended radio sources
with off-axis emission. Here we present radio images of 52 of these sources
that have been made from archival Very Large Array data with resolution of
about 1 arcsec. Fifty-one of the 52 were observed at 1.4 GHz, seven were
observed at 1.4 GHz and 5 GHz, and one was observed only at 5 GHz. We also
present overlays of the SDSS red images for 41 of the sources, and DSS II
overlays for the remainder. Optical counterparts have been identified for most
sources, but there remain a few empty fields.
Our higher resolution VLA images along with FIRST survey images of the
sources in the sample reveal that extended extragalactic radio sources with
small axial ratios are largely (60%) cases of double radio sources with twin
lobes that have off-axis extensions, usually with inversion-symmetric
structure. The available radio images indicate that at most 20% sources might
be genuine X-shaped radio sources that could have formed by a restarting of
beams in a new direction following an interruption and axis flip. The remaining
20% are in neither of these categories. The implications of this result for the
gravitational wave background are discussed in Roberts, Saripalli, and
Subrahmanyan (2015).
Coalescence of super massive black holes (SMBH's) in galactic mergers is potentially the dominant contributor to the low frequency gravitational wave background (GWB). It was proposed by Merritt and Ekers (2002) that X-shaped radio galaxies are signposts of such coalescences, and that their abundance might be used to predict the magnitude of the gravitational wave background. In Roberts et al. (2015) we present radio images of all 52 X-shaped radio source candidates out of the sample of 100 selected by Cheung (2007) for which archival VLA data were available. These images indicate that at most 21% of the candidates might be genuine X-shaped radio sources that were formed by a restarting of beams in a new direction following a major merger. This suggests that fewer than 1.3% of extended radio sources appear to be candidates for genuine axis reorientations, much smaller than the 7% suggested by Leahy and Parma (1992). Thus the associated gravitational wave background may be substantially smaller than previous estimates. These results can be used to normalize detailed calculations of the SMBH coalescence rate and the GWB.
We present a sub-50 pc-scale analysis of the gravitational lens system SDP.81 at redshift 3.042 using Atacama Large submillimetre/Millimetre Array (ALMA) science verification data. These were taken at 236 and 290 GHz using baselines up to 15 km, giving unprecedented insight into the structure of a high-redshift sub-mm galaxy. At mm-wavelengths, the observed system comprises four images in a cusp configuration with an extended, low surface brightness Einstein ring. We model both the mass distribution of the gravitational lensing galaxy and the pixelated surface brightness distribution of the (unlensed) background source using a novel Bayesian technique that fits the data directly in visibility space. We find the mm-wavelength dust emission to be magnified by a factor of u = 17.6 +/- 0.4. The total star-formation rate of the galaxy is 315 +/- 60 M_sol / yr after correcting for the lensing magnification. Our pixelated reconstruction shows the dust emission from SDP.81 to be non-uniform, composed of multiple regions that are heated both by diffuse and by strongly clumped star-formation. We find a possible variation in the spectral slope between the different star-forming regions, which is presumably due a range of dust temperatures within the source. The highest surface brightness region is a ~1.9 x 0.7 kpc disk-like structure, which is surrounded by extended star formation at 20-30 M_sol / yr / kpc^2. The disk contains three compact regions exceeding 120 M_sol / yr / kpc^2, with a maximum of 190 +/- 20 M_sol / yr / kpc^2. This upper limit is below the expectation for Eddington-limited star formation within a radiation-pressure supported starburst.
Reverberation-mapping-based scaling relations are often used to estimate the masses of black holes from single-epoch spectra of AGN. While the radius-luminosity relation that is the basis of these scaling relations is determined using reverberation mapping of the H$\beta$ line in nearby AGN, the scaling relations are often extended to use other broad emission lines, such as MgII, in order to get black hole masses at higher redshifts when H$\beta$ is redshifted out of the optical waveband. However, there is no radius-luminosity relation determined directly from MgII. Here, we present an attempt to perform reverberation mapping using MgII in the well-studied nearby Seyfert 1, NGC 5548. We used Swift to obtain UV grism spectra of NGC 5548 once every two days from April to September 2013. Concurrent photometric UV monitoring with Swift provides a well determined continuum lightcurve that shows strong variability. The MgII emission line, however, is not strongly correlated with the continuum variability, and there is no significant lag between the two. We discuss these results in the context of using MgII scaling relations to estimate high-redshift black hole masses.
We analyze the time evolution of the luminosity of a cluster of Population III protostars formed in the early universe. We argue from the Jeans criterion that primordial gas can collapse to form a cluster of first stars that evolve relatively independently of one another (i.e., with negligible gravitational interaction). We model the collapse of individual protostellar clumps using 2+1D nonaxisymmetric numerical hydrodynamics simulations. Each collapse produces a protostar surrounded by a massive disk (i.e., $M_{\rm disk} / M_{*} \gtrsim 0.1$), whose evolution we follow for a further 30--40 kyr. Gravitational instabilities result in the fragmentation and the formation of gravitationally bound clumps within the disk. The accretion of these fragments by the host protostar produces accretion and luminosity bursts on the order of $10^6\,\LSun$. Within the cluster, we show that a simultaneity of such events across several protostellar cluster members can elevate the cluster luminosity to 5--10${\times}$ greater than expected, and that the cluster spends $\sim15\%$ of it's star-forming history at these levels. This enhanced luminosity effect is particularly enabled in clusters of modest size with $\simeq$ 10--20 members. In one such instance, we identify a confluence of burst events that raise the luminosity to nearly $1000{\times}$ greater than the cluster mean luminosity, resulting in $L > 10^8\,\LSun$. This phenomenon arises solely through the gravitational-instability--driven episodic fragmentation and accretion that characterizes this early stage of protostellar evolution.
To better characterize the global star formation (SF) activity in a galaxy, one needs to know not only the star formation rate (SFR) but also the rest-frame, far-infrared (FIR) color (e.g., the 60-to-100 $\mu$m color, $C(60/100)$] of the dust emission. The latter probes the average intensity of the dust heating radiation field and scales statistically with the effective SFR surface density in star-forming galaxies including (ultra-)luminous infrared galaxies [(U)LIRGs]. To this end, we exploit here a new spectroscopic approach involving only two emission lines: CO\,(7$-$6) at 372 $\mu$m and [NII] at 205 $\mu$m. For local (U)LIRGs, the ratios of the CO (7$-$6) luminosity ($L_{\rm CO\,(7-6)}$) to the total infrared luminosity ($L_{\rm IR}$; 8$-$1000 $\mu$m) are fairly tightly distributed (to within $\sim$0.12 dex) and show little dependence on $C(60/100)$. This makes $L_{\rm CO\,(7-6)}$ a good SFR tracer, which is less contaminated by active galactic nuclei (AGN) than $L_{\rm IR}$ and may also be much less sensitive to metallicity than $L_{\rm CO\,(1-0)}$. Furthermore, the logarithmic [NII] 205 $\mu$m to CO (7$-$6) luminosity ratio is fairly steeply (at a slope of $\sim$$-1.4$) correlated with $C(60/100)$, with a modest scatter ($\sim$0.23 dex). This makes it a useful estimator on $C(60/100)$ with an implied uncertainty of $\sim$0.15 [or $\lesssim$4 K in the dust temperature ($T_{\rm dust}$) in the case of a graybody emission with $T_{\rm dust} \gtrsim 30$ K and a dust emissivity index $\beta \ge 1$]. Our locally calibrated SFR and $C(60/100)$ estimators are shown to be consistent with the published data of (U)LIRGs of $z$ up to $\sim$6.5.
The first N-body simulation of interacting galaxies, even producing spiral arms, was performed by Erik Holmberg in Lund (1941), not with a numerical computer, but by his arrangement of movable light-bulbs and photocells to measure the luminosity at each bulb and thereby estimate the gravitational force. A decade later, and with the arrival of the first programable computers, computations of galactic dynamics were performed, which were later transferred into a N-body simulation movie. I present here the background details for this work with a description of the important elements to note in the movie which may be retrieved at this http URL .
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
We study the kinematics of a local sample of stars, located within a cylinder of 500 pc radius centered on the Sun, in the RAVE dataset. We find clear asymmetries in the $v_R$-$v_\phi$ velocity distributions of thin and thick disk stars: here are more stars moving radially outwards for low azimuthal velocities and more radially inwards for high azimuthal velocities. Such asymmetries have been previously reported for the thin disk as being due to the Galactic bar, but this is the first time that the same type of structures are seen in the thick disk. Our findings imply that the velocities of thick disk stars should no longer be described by Schwarzschild's, multivariate Gaussian or purely axisymmetric distributions. Furthermore, the nature of previously reported substructures in the thick disk needs to be revisited as these could be associated with dynamical resonances rather than to accretion events. It is clear that dynamical models of the Galaxy must fit the 3D velocity distributions of the disks, rather than the projected 1D, if we are to understand the Galaxy fully.
The ionising continuum from active galactic nuclei (AGN) is fundamental for interpreting their broad emission lines and understanding their impact on the surrounding gas. Furthermore, it provides hints on how matter accretes onto supermassive black holes. Using HST's Wide Field Camera 3 we have constructed the first stacked ultraviolet (rest-frame wavelengths 600-2500\AA) spectrum of 53 luminous quasars at z=2.4, with a state-of-the-art correction for the intervening Lyman forest and Lyman continuum absorption. The continuum slope ($f_\nu \propto \nu^{\alpha_\nu}$) of the full sample shows a break at ~912\AA\ with spectral index $\alpha_\nu=-0.61\pm0.01$ at $\lambda>912$\AA\ and a softening at shorter wavelengths ($\alpha_\nu=-1.70 \pm 0.61$ at $\lambda\leq 912$\AA). Our analysis proves that a proper intergalactic medium absorption correction is required to establish the intrinsic continuum emission of quasars. We interpret our average ultraviolet spectrum in the context of photoionisation, accretion disk models, and quasar contribution to the ultraviolet background. We find that observed broad line ratios are consistent with those predicted assuming an ionising slope of $\alpha_\mathrm{ion}=$-2.0, similar to the observed ionising spectrum in the same wavelength range. The continuum break and softening are consistent with accretion disk plus X-ray corona models when black hole spin is taken into account. Our spectral energy distribution yields a 30% increase to previous estimates of the specific quasar emissivity, such that quasars may contribute significantly to the total specific Lyman limit emissivity estimated from the Ly$\alpha$ forest at z<3.2.
Cluster mergers may play a fundamental role in the formation and evolution of cluster galaxies. Stroe et al. (2014) revealed unexpected over-densities of candidate H$\alpha$ emitters near the ~1 Mpc-wide shock fronts of the massive (~2x10$^{15}$M$_{\odot}$) "Sausage" merging cluster, CIZA J2242.8+5301. We used Keck/DEIMOS and WHT/AF2 to confirm 83 H$\alpha$ emitters in and around the merging cluster. We find that cluster star-forming galaxies in the hottest X-ray gas and/or in the cluster sub-cores (away from the shock fronts) show high [SII]6716/[SII]6761 and high [SII]6716/H$\alpha$, implying very low electron densities (<30x lower than all other star-forming galaxies outside the cluster) and significant contribution from supernovae, respectively. All cluster star-forming galaxies near the cluster centre show evidence of significant outflows (blueshifted Na D~200-300km/s), likely driven by supernovae. Strong outflows are also found for the cluster H$\alpha$ AGN. H$\alpha$ star-forming galaxies in the merging cluster follow the z~0 mass-metallicity relation, showing systematically higher metallicity (~0.15-0.2 dex) than H$\alpha$ emitters outside the cluster (projected R>2.5 Mpc). This suggests that the shock front may have triggered remaining metal-rich gas which galaxies were able to retain into forming stars. Our observations show that the merger of impressively massive (~10$^{15}$M$_\odot$) clusters can provide the conditions for significant star-formation and AGN activity, but, as we witness strong feedback by star-forming galaxies and AGN (and given how massive the merging cluster is), such sources will likely quench in a few 100 Myrs.
We have used the publicly released Dark Energy Survey data to hunt for new satellites of the Milky Way in the Southern hemisphere. Our search yielded a large number of promising candidates. In this paper, we announce the discovery of 9 new unambiguous ultra-faint objects, whose authenticity can be established with the DES data alone. Based on the morphological properties, three of the new satellites are dwarf galaxies, one of which is located at the very outskirts of the Milky Way, at a distance of 380 kpc. The remaining 6 objects have sizes and luminosities comparable to the Segue 1 satellite and can not be classified straightforwardly without follow-up spectroscopic observations. The satellites we have discovered cluster around the LMC and the SMC. We show that such spatial distribution is unlikely under the assumption of isotropy, and, therefore, conclude that at least some of the new satellites must have been associated with the Magellanic Clouds in the past.
We analyze the low-redshift (z~0.2) circumgalactic medium by comparing absorption-line data from the COS-Halos Survey to absorption around a matched galaxy sample from two cosmological hydrodynamic simulations. The models include different prescriptions for galactic outflows, namely hybrid energy/momentum driven wind (ezw), and constant winds (cw). We extract for comparison direct observables including equivalent widths, covering factors, ion ratios, and kinematics. Both wind models are generally in good agreement with these observations for HI and certain low ionization metal lines, but show poorer agreement with higher ionization metal lines including SiIII and OVI that are well-observed by COS-Halos. These discrepancies suggest that both wind models predict too much cool, metal-enriched gas and not enough hot gas, and/or that the metals are not sufficiently well-mixed. This may reflect our model assumption of ejecting outflows as cool and unmixing gas. Our ezw simulation includes a heuristic prescription to quench massive galaxies by super-heating its ISM gas, which we show yields sufficient low ionisation absorption to be broadly consistent with observations, but also substantial OVI absorption that is inconsistent with data, suggesting that gas around quenched galaxies in the real Universe does not cool. At impact parameters of <50 kpc, recycling winds dominate the absorption of low ions and even HI, while OVI almost always arises from metals ejected longer than 1 Gyr ago. The similarity between the wind models is surprising, since we show that they differ substantially in their predicted amount and phase distribution of halo gas. We show that this similarity owes mainly to our comparison here at fixed stellar mass rather than at fixed halo mass in our previous works, which suggests that CGM properties are more closely tied to the stellar mass of galaxies rather than halo mass.
We present an X-ray and multiwavelength study of 33 weak emission-line quasars (WLQs) and 18 quasars that are analogs of the extreme WLQ, PHL 1811, at z ~ 0.5-2.9. New Chandra 1.5-9.5 ks exploratory observations were obtained for 32 objects while the others have archival X-ray observations. Significant fractions of these luminous type 1 quasars are distinctly X-ray weak compared to typical quasars, including 16 (48%) of the WLQs and 17 (94%) of the PHL 1811 analogs with average X-ray weakness factors of 17 and 39, respectively. We measure a relatively hard ($\Gamma=1.16_{-0.32}^{+0.37}$) effective power-law photon index for a stack of the X-ray weak subsample, suggesting X-ray absorption, and spectral analysis of one PHL 1811 analog, J1521+5202, also indicates significant intrinsic X-ray absorption. We compare composite SDSS spectra for the X-ray weak and X-ray normal populations and find several optical-UV tracers of X-ray weakness; e.g., Fe II rest-frame equivalent width and relative color. We describe how orientation effects under our previously proposed "shielding-gas" scenario can likely unify the X-ray weak and X-ray normal populations. We suggest that the shielding gas may naturally be understood as a geometrically thick inner accretion disk that shields the broad line region from the ionizing continuum. If WLQs and PHL 1811 analogs have very high Eddington ratios, the inner disk could be significantly puffed up (e.g., a slim disk). Shielding of the broad emission-line region by a geometrically thick disk may have a significant role in setting the broad distributions of C IV rest-frame equivalent width and blueshift for quasars more generally.
The discovery of multiple stellar populations in globular clusters has implications for all the aspects of the study of these stellar systems. In this paper, by means of N-body simulations, we study the evolution of binary stars in multiple-population clusters and explore the implications of the initial differences in the spatial distribution of different stellar populations for the evolution and survival of their binary stars. Our simulations show that initial differences between the spatial distribution of first-generation (FG) and second-generation (SG) stars can leave a fingerprint in the current properties of the binary population. SG binaries are disrupted more efficiently than those of the FG population resulting in a global SG binary fraction smaller than that of the FG. As for surviving binaries, dynamical evolution produces a difference between the SG and the FG binary binding energy distribution with the SG population characterized by a larger fraction of high binding energy (more bound) binaries. We have also studied the dependence of the binary properties on the distance from the cluster centre. Although the global binary fraction decreases more rapidly for the SG population, the local binary fraction measured in the cluster inner regions may still be dominated by SG binaries. The extent of the differences between the surviving FG and SG binary binding energy distribution also varies radially within the cluster and is larger in the cluster inner regions.
Most of the massive elliptical galaxies in the universe stopped forming stars billions of years ago, even though plenty of hot gas remains available for star formation. Here we present compelling evidence indicating that quenching of star formation depends on both black-hole feedback and Type Ia supernova heating. We analyze Chandra X-ray observations of ten massive ellipticals, five with extended, potentially star-forming multiphase gas and five single-phase ellipticals with no star formation. The ratio of cooling time to freefall time at 1--10 kpc in the multiphase galaxies is tc/tff ~10, indicating that precipitation-driven feedback limits cooling but does not eliminate condensation. In the same region of the single-phase galaxies, the radial profiles of gas entropy are consistent with a thermally stable (tc/tff > 20) supernova-driven outflow that sweeps stellar ejecta out of the galaxy. However, in one of those single-phase ellipticals (NGC 4261) we find tc/tff < 10 at < 300 pc. Notably, its jets are ~50 times more powerful than in the other nine ellipticals, in agreement with models indicating that precipitation near the black hole should switch its fueling mode from Bondi-like accretion to cold chaotic accretion. We conclude by hypothesizing that black-hole feedback outbursts shut off star formation in massive elliptical galaxies by switching on supernova sweeping capable of raising the entropy of the outflowing stellar ejecta to > 100 keV cm^2.
Using a sample of 69,919 red giants from the SDSS-III/APOGEE Data Release 12, we measure the distribution of stars in the [$\alpha$/Fe] vs. [Fe/H] plane and the metallicity distribution functions (MDF) across an unprecedented volume of the Milky Way disk, with radius $3<R<15$ kpc and height $|z|<2$ kpc. Stars in the inner disk ($R<5$ kpc) lie along a single track in [$\alpha$/Fe] vs. [Fe/H], starting with $\alpha$-enhanced, metal-poor stars and ending at [$\alpha$/Fe]$\sim0$ and [Fe/H]$\sim+0.4$. At larger radii we find two distinct sequences in [$\alpha$/Fe] vs. [Fe/H] space, with a roughly solar-$\alpha$ sequence that spans a decade in metallicity and a high-$\alpha$ sequence that merges with the low-$\alpha$ sequence at super-solar [Fe/H]. The location of the high-$\alpha$ sequence is nearly constant across the disk, however there are very few high-$\alpha$ stars at $R>11$ kpc. The peak of the midplane MDF shifts to lower metallicity at larger $R$, reflecting the Galactic metallicity gradient. Most strikingly, the shape of the midplane MDF changes systematically with radius, with a negatively skewed distribution at $3<R<7$ kpc, to a roughly Gaussian distribution at the solar annulus, to a positively skewed shape in the outer Galaxy. For stars with $|z|>1$ kpc or [$\alpha$/Fe]$>0.18$, the MDF shows little dependence on $R$. The positive skewness of the outer disk MDF may be a signature of radial migration; we show that blurring of stellar populations by orbital eccentricities is not enough to explain the reversal of MDF shape but a simple model of radial migration can do so.
We report the redshift of an unlensed, highly obscured submillimetre galaxy (SMG), HS1700.850.1, the brightest SMG (S850um =19.1 mJy) detected in the JCMT/SCUBA-2 Baryonic Structure Survey, based on the detection of its CO line emission. Using the IRAM PdBI-WIDEX with 3.6GHz band width, we serendipitously detect an emission line at 150.6 GHz. Confirmation of the identification of this line as CO(5-4) at z = 2.816 comes from a search over 14.5 GHz in the 3-mm and 2-mm atmospheric windows, meaning that it does not reside in the z~2.30 proto-cluster in this field. Measurement of the 870um source size (<0.74") from the Sub-Millimeter Array (SMA) confirms a compact emission in a S870um =14.5mJy, LIR~10^13 Lsun component, suggesting an Eddington-limited starburst. We use the double-peaked CO line profile measurements along with the SMA size constraints to study the gas dynamics of a HyLIRG, estimating the gas and dynamical masses of HS1700.850.1. While HS1700.850.1 is one of the most extreme galaxies known in the Universe, we find that it occupies a relative void in the Lyman-Break Galaxy distribution in this field. Comparison with other extreme objects at similar epochs (HyLIRG Quasars), and cosmological simulations, suggests such an anti-bias of bright SMGs could be relatively common, with the brightest SMGs rarely occupying the most overdense regions at z=2-4.
The Taurus-Auriga molecular complex (TMC) is the main laboratory for the study of low mass star formation. The density and properties of interstellar dust are expected to vary across the TMC. These variations trace important processes such as dust nucleation or the magnetic field coupling with the cloud. In this article, we show how the combination of near ultraviolet (NUV) and infrared (IR) photometry can be used to derive the strength of the 2175 \AA\ bump and thus any enhancement in the abundance of small dust grains and PAHs in the dust grains size distribution. This technique is applied to the envelope of the TMC, mapped by the GALEX All Sky Survey (AIS). UV and IR photometric data have been retrieved from the GALEX-AIS and the 2MASS catalogues. NUV and K-band star counts have been used to identify the areas in the cloud envelope where the 2175 \AA\ bump is weaker than in the diffuse ISM namely, the low column density extensions of L1495, L1498 and L1524 in Taurus, L1545, L1548, L1519, L1513 in Auriga and L1482-83 in the California region. This finding agrees with previous results on dust evolution derived from Spitzer data and suggests that dust grains begin to decouple from the environmental galactic magnetic field already in the envelope.
The spectral index of synchrotron emission is an important parameter in understanding the properties of cosmic ray electrons (CREs) and the interstellar medium (ISM). We determine the synchrotron spectral index ($\alpha_{\rm nt}$) of four nearby star-forming galaxies, namely NGC 4736, NGC 5055, NGC 5236 and NGC 6946 at sub-kpc linear scales. The $\alpha_{\rm nt}$ was determined between 0.33 and 1.4 GHz for all the galaxies. We find the spectral index to be flatter ($\gtrsim -0.7$) in regions with total neutral (atomic + molecular) gas surface density, $\Sigma_{\rm gas} \gtrsim \rm 50~M_\odot pc^{-2}$, typically in the arms and inner parts of the galaxies. In regions with $\Sigma_{\rm gas} \lesssim \rm 50~M_\odot pc^{-2}$, especially in the interarm and outer regions of the galaxies, the spectral index steepens sharply to $<-1.0$. The flattening of $\alpha_{\rm nt}$ is unlikely to be caused due to thermal free--free absorption at 0.33 GHz. Our result is consistent with the scenario where the CREs emitting at frequencies below $\sim0.3$ GHz are dominated by bremsstrahlung and/or ionization losses. For denser medium ($\Sigma_{\rm gas} \gtrsim \rm 200~M_\odot pc^{-2}$), having strong magnetic fields ($\sim 30~\mu$G), $\alpha_{\rm nt}$ is seen to be flatter than $-0.5$, perhaps caused due to ionization losses. We find that, due to the clumpy nature of the ISM, such dense regions cover only a small fraction of the galaxy ($\lesssim5$ percent). Thus, the galaxy-integrated spectrum may not show indication of such loss mechanisms and remain a power-law over a wide range of radio frequencies (between $\sim 0.1$ to 10 GHz).
A growing body of evidence has been supporting the existence of so-called "dark molecular gas" (DMG), which is invisible in the most common tracer of molecular gas, i.e., CO rotational emission. DMG is believed to be the main gas component of the intermediate extinction region between A$\rm_v$$\sim$0.05-2, roughly corresponding to the self-shielding threshold of H$_2$ and $^{13}$CO. To quantify DMG relative to HI and CO, we are pursuing three observational techniques, namely, HI self-absorption, OH absorption, and TeraHz C$^+$ emission. In this paper, we focus on preliminary results from a CO and OH absorption survey of DMG candidates. Our analysis show that the OH excitation temperature is close to that of the Galactic continuum background and that OH is a good DMG tracer co-existing with molecular hydrogen in regions without CO. Through systematic "absorption mapping" by Square Kilometer Array (SKA) and ALMA, we will have unprecedented, comprehensive knowledge of the ISM components including DMG in terms of their temperature and density, which will impact our understanding of galaxy evolution and star formation profoundly.
As part of a long-term project to revisit the kinematics and dynamics of the large disc galaxies of the Local Group, we present the first deep, wide-field (42' x 56') 3D-spectroscopic survey of the ionized gas disc of Messier 33. Fabry-Perot interferometry has been used to map its Ha distribution and kinematics at unprecedented angular resolution (<3'') and resolving power (12600), with the 1.6m telescope at the Observatoire du Mont Megantic. The ionized gas distribution follows a complex, large-scale spiral structure, unsurprisingly coincident with the already-known spiral structures of the neutral and molecular gas discs. The kinematical analysis of the velocity field shows that the rotation center of the Ha disc is distant from the photometric center by 170 pc (sky projected distance) and that the kinematical major-axis position angle and disc inclination are in excellent agreement with photometric values. The Ha rotation curve agrees very well with the HI rotation curves for 0 < R < 6.5 kpc, but the Ha velocities are 10-20 km/s higher for R > 6.5 kpc. The reason for this discrepancy is not well understood. The velocity dispersion profile is relatively flat around 16 km/s, which is at the low end of velocity dispersions of nearby star-forming galactic discs. A strong relation is also found between the Ha velocity dispersion and the Ha intensity. Mass models were obtained using the Ha rotation curve but, as expected, the dark matter halo's parameters are not very well constrained since the optical rotation curve only extends out to 8 kpc.
We present JCMT SCUBA-2 850microns submillimetre (submm) observations of 30 mid-infrared (mid-IR) luminous AGN, detected jointly by the WISE all-sky IR survey and the NVSS/FIRST radio survey. These rare sources are selected by their extremely red mid-infrared spectral energy distributions (SEDs) and compact radio counterparts. Further investigations show that they are highly obscured, have abundant warm AGN-heated dust and are thought to be experiencing intense AGN feedback. These galaxies appear to be consistent with an AGN-dominated galaxy, and could be a transient phase of merging galaxies. When comparing the number of submm galaxies (SMGs) detected serendipitously in the surrounding 1.5-arcmin to those in blank-field submm surveys, there is a very significant overdensity, of order 5, but no sign of radial clustering centred at our primary objects. The WISE/radio-selected AGN thus reside in 10-Mpc-scale overdense environments, that could be forming in pre-viralised clusters of galaxies. WISE/radio-selected AGNs appear to be the strongest signposts of high-density regions of active, luminous and dusty galaxies. SCUBA-2 850microns observations indicate that their submm fluxes are low compared to many popular AGN SED templates, hence the WISE/radio-selected AGNs have either less cold and/or more warm dust emission than normally assumed for typical AGN. Most of the targets are not detected, only four targets are detected at SCUBA-2 850microns, and have total IR luminosities >= 10^13 L_solar, if their redshifts are consistent with the subset of the 10 SCUBA-2 undetected targets with known redshifts, z ~ 0.44 - 2.86.
We report the discovery of eight new Milky Way companions in ~1,800 deg^2 of optical imaging data collected during the first year of the Dark Energy Survey (DES). Each system is identified as a statistically significant over-density of individual stars consistent with the expected isochrone and luminosity function of an old and metal-poor stellar population. The objects span a wide range of absolute magnitudes (M_V from -2.2 mag to -7.4 mag), physical sizes (10 pc to 170 pc), and heliocentric distances (30 kpc to 330 kpc). Based on the low surface brightnesses, large physical sizes, and/or large Galactocentric distances of these objects, several are likely to be new ultra-faint satellite galaxies of the Milky Way and/or Magellanic Clouds. We introduce a likelihood-based algorithm to search for and characterize stellar over-densities, as well as identify stars with high satellite membership probabilities. We also present completeness estimates for detecting ultra-faint galaxies of varying luminosities, sizes, and heliocentric distances in the first-year DES data.
We introduce a Bayesian method for fitting faint, resolved stellar spectra in order to obtain simultaneous estimates of redshift and stellar-atmospheric parameters. We apply the method to thousands of spectra---covering 5160-5280 Angs. at resolution R~20,000---that we have acquired with the MMT/Hectochelle fibre spectrograph for red-giant and horizontal branch candidates along the line of sight to the Milky Way's dwarf spheroidal satellite in Draco. The observed stars subtend an area of ~4 deg^2, extending ~3 times beyond Draco's nominal `tidal' radius. For each spectrum we tabulate the first four moments---central value, variance, skewness and kurtosis---of posterior probability distribution functions representing estimates of the following physical parameters: line-of-sight velocity v_los, effective temperature (T_eff), surface gravity (logg) and metallicity ([Fe/H]). After rejecting low-quality measurements, we retain a new sample consisting of 2813 independent observations of 1565 unique stars, including 1879 observations for 631 stars with (as many as 13) repeat observations. Parameter estimates have median random errors of sigma_{v_los}=0.88 km/s, sigma_{T_eff}=162 K, sigma_logg=0.37 dex and sigma_[Fe/H]=0.20 dex. Our estimates of physical parameters distinguish ~470 likely Draco members from interlopers in the Galactic foreground.
About one third of early-type barred galaxies host small-scale secondary bars. The formation and evolution of such double-barred galaxies remain far from being well understood. In order to understand the formation of such systems, we explore a large parameter space of isolated pure-disk simulations. We show that a dynamically cool inner disk embedded in a hotter outer disk can naturally generate a steady secondary bar while the outer disk forms a large-scale primary bar. The independent bar instabilities of inner and outer disks result in long-lived double-barred structures whose dynamical properties are comparable with observations. This formation scenario indicates that the secondary bar might form from the general bar instability, the same as the primary bar. Under some circumstances, the interaction of the bars and the disk leads to the two bars aligning or single, nuclear, bars only. Simulations that are cool enough of the center to experience clump instabilities may also generate steady double-barred galaxies. In this case, the secondary bars are "fast", i.e., the bar length is close to the co-rotation radius. This is the first time that double-barred galaxies containing a fast secondary bar are reported. Previous orbit-based studies had suggested that fast secondary bars are not dynamically possible.
Dust lanes, nuclear rings, and nuclear spirals are typical gas structures in the inner region of barred galaxies Their shapes and properties are linked to the physical parameters of the host galaxy. We use high-resolution hydrodynamical simulations to study 2D gas flows in simple barred galaxy models. The nuclear rings formed in our simulations can be divided into two groups: one group is nearly round and the other is highly elongated. We find that roundish rings may not form when the bar pattern speed is too high or the bulge central density is too low. We also study the periodic orbits in our galaxy models, and find that the concept of inner Lindblad resonance (ILR) may be generalized by the extent of $x_2$ orbits. All roundish nuclear rings in our simulations settle in the range of $x_2$ orbits (or ILRs). However, knowing the resonances is insufficient to pin down the exact location of these nuclear rings. We suggest that the backbone of round nuclear rings is the $x_2$ orbital family, i.e. round nuclear rings are allowed only in the radial range of $x_2$ orbits. A round nuclear ring forms exactly at the radius where the residual angular momentum of infalling gas balances the centrifugal force, which can be described by a parameter $f_{\rm ring}$ measured from the rotation curve. The gravitational torque on gas in high pattern speed models is larger, leading to a smaller ring size than in the low pattern speed models. Our result may have important implications for using nuclear rings to measure the parameters of real barred galaxies with 2D gas kinematics.
Under cosmic irradiation, the interstellar water ice mantles evolve towards a compact amorphous state. Crystalline ice amorphisation was previously monitored mainly in the keV to hundreds of keV ion energies. We experimentally investigate heavy ion irradiation amorphisation of crystalline ice, at high energies closer to true cosmic rays, and explore the water-ice sputtering yield. We irradiated thin crystalline ice films with MeV to GeV swift ion beams, produced at the GANIL accelerator. The ice infrared spectral evolution as a function of fluence is monitored with in-situ infrared spectroscopy (induced amorphisation of the initial crystalline state into a compact amorphous phase). The crystalline ice amorphisation cross-section is measured in the high electronic stopping-power range for different temperatures. At large fluence, the ice sputtering is measured on the infrared spectra, and the fitted sputtering-yield dependence, combined with previous measurements, is quadratic over three decades of electronic stopping power. The final state of cosmic ray irradiation for porous amorphous and crystalline ice, as monitored by infrared spectroscopy, is the same, but with a large difference in cross-section, hence in time scale in an astrophysical context. The cosmic ray water-ice sputtering rates compete with the UV photodesorption yields reported in the literature. The prevalence of direct cosmic ray sputtering over cosmic-ray induced photons photodesorption may be particularly true for ices strongly bonded to the ice mantles surfaces, such as hydrogen-bonded ice structures or more generally the so-called polar ices.
The Wide-field Nearby Galaxy-cluster Survey (WINGS) is a wide-field
multi-wavelength survey of X-ray selected clusters at z =0.04-0.07. The
original 34'x34' WINGS field-of- view has now been extended to cover a 1 sq.deg
field with both photometry and spectroscopy. In this paper we present the
Johnson B and V-band OmegaCAM/VST observations of 46 WINGS clusters, together
with the data reduction, data quality and Sextractor photometric catalogs.
With a median seeing of 1arcs in both bands, our 25-minutes exposures in each
band typically reach the 50% completeness level at V=23.1 mag. The quality of
the astrometric and photometric accuracy has been verified by comparison with
the 2MASS as well as with SDSS astrometry, and SDSS and previous WINGS imaging.
Star/galaxy separation and sky-subtraction procedure have been tested comparing
with previous WINGS data.
The Sextractor photometric catalogues are publicly available at the CDS, and
will be included in the next release of the WINGS database on the VO together
with the OmegaCAM reduced images. These data form the basis for a large ongoing
spectroscopic campaign with AAOmega/AAT and is being employed for a variety of
studies. [abridged]
We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.04 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using an extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at angular resolutions as fine as 23 milliarcseconds (mas; corresponding to an un-magnified spatial scale of 180 pc at z=3.042). The ALMA images clearly show two main gravitational arc components with emission tracing a radius of 1.5 arcseconds. We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO data has an angular resolution of 170 mas and the emission is found to broadly trace the gravitational arc structures. We detect H2O line emission but only using the shortest baselines. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the increase in angular resolution. Finally, we detect weak unresolved continuum emission at all three observed frequencies from a position that is spatially coincident with the centre of the foreground lensing galaxy.
We present the spectroscopic redshift catalog from a wide-field survey of the fields of 28 galaxy-mass strong gravitational lenses. We discuss the acquisition and reduction of the survey data, collected over 40 nights of 6.5m MMT and Magellan time, employing four different multi-object spectrographs. We determine that no biases are introduced by combining datasets obtained with different instrument/spectrograph combinations. Special care is taken to determine redshift uncertainties using repeat observations. The redshift catalog consists of 9768 new and unique galaxy redshifts. 82.4% of the catalog redshifts are between z=0.1 and z=0.7, and the catalog median redshift is z=0.36. The data from this survey will be used to study the lens environments and line-of-sight structures to gain a better understanding of the effects of large scale structure on lens statistics and lens-derived parameters.
Context. Stars are born together from giant molecular clouds and, if we
assume that the priors were chemically homogeneous and well-mixed, we expect
them to share the same chemical composition. Most of the stellar aggregates are
disrupted while orbiting the Galaxy and most of the dynamic information is
lost, thus the only possibility of reconstructing the stellar formation history
is to analyze the chemical abundances that we observe today.
Aims. The chemical tagging technique aims to recover disrupted stellar
clusters based merely on their chemical composition. We evaluate the viability
of this technique to recover co-natal stars that are no longer gravitationally
bound.
Methods. Open clusters are co-natal aggregates that have managed to survive
together. We compiled stellar spectra from 31 old and intermediate-age open
clusters, homogeneously derived atmospheric parameters, and 17 abundance
species, and applied machine learning algorithms to group the stars based on
their chemical composition. This approach allows us to evaluate the viability
and efficiency of the chemical tagging technique.
Results. We found that stars at different evolutionary stages have distinct
chemical patterns that may be due to NLTE effects, atomic diffusion, mixing,
and biases. When separating stars into dwarfs and giants, we observed that a
few open clusters show distinct chemical signatures while the majority show a
high degree of overlap. This limits the recovery of co-natal aggregates by
applying the chemical tagging technique. Nevertheless, there is room for
improvement if more elements are included and models are improved.
We investigate the luminosity and redshift dependence of the quasar continuum by means of composite spectrum using a large non-BAL radio-quiet quasar sample drawn from the Sloan Digital Sky Survey. Quasar continuum slopes in the UV-Opt band are measured at two different wavelength ranges, i.e., $\alpha_{\nu12}$ ($1000\sim 2000 \rm\AA$) and $\alpha_{\nu24}$ ($2000 \sim 4000 \rm\AA$) derived from power law fitting. Generally, the UV spectra slope becomes harder (higher $\alpha_{\nu}$) towards higher bolometric luminosity. On the other hand, when quasars are further grouped into luminosity bins, we find both $\alpha_{\nu12}$ and $\alpha_{\nu24}$ show significant anti-correlation with redshift (i.e., quasar continuum becomes redder towards higher redshift). We suggest that the cosmic dust extinction is very likely the cause of this observed $\alpha_\nu-z$ relation. We build a simple cosmic dust extinction model to quantify the observed reddening tendency and find an effective dust density $n\sigma_v \sim 10^{-5}h~\rm Mpc^{-1}$ at $z<1.5$. The other possibilities that could produce such a reddening effect have also been discussed.
Recent discoveries of super-massive black holes at high redshifts indicate a possible tension with the standard Lambda CDM paradigm of early universe cosmology which has difficulties in explaining the origin of the required nonlinear compact seeds which trigger the formation of these super-massive black holes. Here we show that cosmic string loops which result from a scaling solution of strings formed during a phase transition in the very early universe lead to an additional source of compact seeds. The number density of string-induced seeds dominates at high redshifts and can help trigger the formation of the observed super-massive black holes.
Measuring radio source counts is critical for characterizing new extragalactic populations, brings a wealth of science within reach and will inform forecasts for SKA and its pathfinders. Yet there is currently great debate (and few measurements) about the behaviour of the 1.4-GHz counts in the microJy regime. One way to push the counts to these levels is via 'stacking', the covariance of a map with a catalogue at higher resolution and (often) a different wavelength. For the first time, we cast stacking in a fully bayesian framework, applying it to (i) the SKADS simulation and (ii) VLA data stacked at the positions of sources from the VIDEO survey. In the former case, the algorithm recovers the counts correctly when applied to the catalogue, but is biased high when confusion comes into play. This needs to be accounted for in the analysis of data from any relatively-low-resolution SKA pathfinders. For the latter case, the observed radio source counts remain flat below the 5-sigma level of 85 microJy as far as 40 microJy, then fall off earlier than the flux hinted at by the SKADS simulations and a recent P(D) analysis (which is the only other measurement from the literature at these flux-density levels, itself extrapolated in frequency). Division into galaxy type via spectral-energy distribution reveals that normal spiral galaxies dominate the counts at these fluxes.
Due to their proximity, high dark matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels.
Hierarchical structure formation implies that the number of subhalos within a dark matter halo depends not only on halo mass, but also on the formation history of the halo. This dependence on the formation history, which is highly correlated with halo concentration, can account for the super-Poissonian scatter in subhalo occupation at a fixed halo mass that has been previously measured in simulations. Here we propose a model to predict the subhalo abundance function for individual host halos, that incorporates both halo mass and concentration. We combine results of cosmological simulations with a new suite of zoom-in simulations of Milky Way-mass halos to calibrate our model. We show the model can successfully reproduce the mean and the scatter of subhalo occupation in these simulations. The implications of this correlation between subhalo abundance and halo concentration are further investigated. We also discuss cases in which inferences about halo properties can be affected if this correlation between subhalo abundance and halo concentration is ignored; in these cases our model would give a more accurate inference. We propose that with future deep surveys, satellite occupation in the low-mass regime can be used to verify the existence of halo assembly bias.
Heat input roughly balances radiative cooling in the gaseous cores of galaxy clusters even when the central cooling time is short, implying that cooling triggers a feedback loop that maintains thermal balance. Furthermore, cores with short cooling times tend to have multiphase structure, suggesting that the intracluster medium (ICM) becomes locally thermally unstable for cooling times < 1 Gyr. In this work, we use 2D and 3D hydrodynamic simulations to study the onset of condensation in idealized galaxy-cluster cores. In particular, we look at how the condensation process depends on the ratio of cooling time to freefall time and on the geometry of the gravitational potential. We conclude that the ICM can always evolve to a state in which condensation occurs if given enough time, but that an initial timescale ratio tcool /tff < 10 is needed for thermal instability to grow quickly enough to affect realistic cluster cores within a timescale that is relevant for cosmological structure formation. We find that instability leads to convection and that perturbations continue to grow while the gas convects. Condensation occurs when the timescale ratio in the low-entropy tail of the perturbation distribution drops below tcool /tff < 3, even if the volume-averaged timescale ratio is substantially greater. In our simulations, the geometry of the gravitational potential does not have a strong effect on thermal stability. Finally, we find that if condensation is powering feedback, a conversion efficiency of around 10^-3 for converting the condensed mass into thermal energy is sufficient to maintain thermal balance in the ICM.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
In this paper, we present the properties of 10 halo globular clusters with luminosities $L\simeq 5-7\times 10^5{L_\odot}$ in the Local Group galaxy M33 using the images of {\it Hubble Space Telescope} Wide Field Planetary Camera 2 in the F555W and F814W bands. We obtained ellipticities, position angles and surface brightness profiles for them. In general, the ellipticities of M33 sample clusters are similar to those of M31 clusters. The structural and dynamical parameters are derived by fitting the profiles to three different models combined with mass-to-light ratios ($M/L$ values) from population-synthesis models. The structural parameters include core radii, concentration, half-light radii {\bf and} central surface brightness. The dynamical parameters include the integrated cluster mass, integrated binding energy, central surface mass density {\bf and} predicted line-of-sight velocity dispersion at the cluster center. The velocity dispersions of four clusters predicted here agree well with the observed dispersions by Larsen et al. The results here showed that the majority of the sample halo globular clusters are well fitted by King model as well as by Wilson model, and better than by S\'ersic model. In general, the properties of clusters in M33, M31 and the Milky Way fall in the same regions of parameter spaces. The tight correlations of cluster properties indicate a "fundamental plane" for clusters, which reflects some universal physical conditions and processes operating at the epoch of cluster formation.
Numerical simulations of active galactic nuclei (AGN) feedback in cool-core galaxy clusters have successfully avoided classical cooling flows, but often produce too much cold gas. We perform adaptive mesh simulations that include momentum-driven AGN feedback, self-gravity, star formation and stellar feedback, focusing on the interplay between cooling, AGN heating and star formation in an isolated cool-core cluster. Cold clumps triggered by AGN jets and turbulence form filamentary structures tens of kpc long. This cold gas feeds both star formation and the supermassive black hole (SMBH), triggering an AGN outburst that increases the entropy of the ICM and reduces its cooling rate. Within 1-2 Gyr, star formation completely consumes the cold gas, leading to a brief shutoff of the AGN. The ICM quickly cools and redevelops multiphase gas, followed by another cycle of star formation/AGN outburst. Within 6.5 Gyr, we observe three such cycles. There is good agreement between our simulated cluster and the observations of cool-core clusters. ICM cooling is dynamically balanced by AGN heating, and a cool-core appearance is preserved. The minimum cooling time to free-fall time ratio typically varies between a few and $\gtrsim 20$. The star formation rate (SFR) covers a wide range, from 0 to a few hundred $\rm M_{\odot}\, yr^{-1}$, with an average of $\sim 40 \,\rm M_{\odot}\, yr^{-1}$. The instantaneous SMBH accretion rate shows large variations on short timescales, but the average value correlates well with the SFR. Simulations without stellar feedback or self-gravity produce qualitatively similar results, but a lower SMBH feedback efficiency (0.1% compared to 1%) results in too many stars.
We present results from a blind survey to identify strong gravitational lenses among the population of low-redshift early-type galaxies. The SINFONI Nearby Elliptical Lens Locator Survey (SNELLS) uses integral-field infrared spectroscopy to search for lensed emission line sources behind massive lens candidates at $z$<0.055. From 27 galaxies observed, we have recovered one previously-known lens (ESO325-G004) at $z$=0.034, and discovered two new systems, at $z$=0.031 and $z$=0.052. All three lens galaxies have high velocity dispersions (\sigma>300 km/s) and \alpha-element abundances ([Mg/Fe]>0.3). From the lensing configurations we derive total J-band mass-to-light ratios of 1.8$\pm$0.1, 2.1$\pm$0.1 and 1.9$\pm$0.2 within the $\sim$2 kpc Einstein radius. Correcting for estimated dark-matter contributions, and comparing to stellar population models with a Milky Way (Kroupa) initial mass function (IMF), we determine the "mass excess factor", \alpha. Assuming the lens galaxies have "old" stellar populations (10$\pm$1 Gyr), the average IMF mass factor is $\langle\alpha\rangle$=1.10$\pm$0.08$\pm$0.10, where the first error is random and the second is systematic. If we instead fit the stellar populations from 6dF optical survey spectra, all three galaxies are consistent with being old, but the age errors are 3-4 Gyr, due to limited signal-to-noise. The IMF constraints are therefore looser in this case, with $\langle\alpha\rangle$ = $1.23^{+0.16}_{-0.13}\pm{0.10}$. Our results are thus consistent with a Kroupa IMF (\alpha=1.00) on average, and strongly reject very heavy IMFs with \alpha>2. A Salpeter IMF (\alpha=1.55) is inconsistent at the 3.5$\sigma$ level if the galaxies are old, but cannot be excluded using age constraints derived from the currently-available optical spectra.
We study the Kennicutt-Schmidt relation between average star formation rate and average cold gas surface density in the Hi dominated ISM of nearby spiral and dwarf irregular galaxies. We divide the galaxies into grid cells varying from sub-kpc to tens of kpc in size. Grid-cell measurements of low SFRs using H-alpha emission can be biased and scatter may be introduced because of non-uniform sampling of the IMF or because of stochastically varying star formation. In order to alleviate these issues, we use far-ultraviolet emission to trace SFR, and we sum up the fluxes from different bins with the same gas surface density to calculate the average $\Sigma_{SFR}$ at a given value of $\Sigma_{gas}$. We study the resulting Kennicutt-Schmidt relation in 400 pc, 1 kpc and 10 kpc scale grids in nearby massive spirals and in 400 pc scale grids in nearby faint dwarf irregulars. We find a relation with a power law slope of 1.5 in the HI-dominated regions for both kinds of galaxies. The relation is offset towards longer gas consumption timescales compared to the molecular hydrogen dominated centres of spirals, but the offset is an order-of-magnitude less than that quoted by earlier studies. Our results lead to the surprising conclusion that conversion of gas to stars is independent of metallicity in the HI dominated regions of star-forming galaxies. Our observed relations are better fit by a model of star formation based on thermal and hydrostatic equilibrium in the ISM, in which feedback driven turbulence sets the thermal pressure.
We present a detailed gravitational lens model of the galaxy cluster RCS2 J232727.6-020437. Due to cosmological dimming of cluster members and ICL, its high redshift ($z=0.6986$) makes it ideal for studying background galaxies. Using new ACS and WFC3/IR HST data, we identify 16 multiple images. From MOSFIRE follow up, we identify a strong emission line in the spectrum of one multiple image, likely confirming the redshift of that system to $z=2.083$. With a highly magnified ($\mu\gtrsim2$) source plane area of $\sim0.7$ arcmin$^2$ at $z=7$, RCS2 J232727.6-020437 has a lensing efficiency comparable to the Hubble Frontier Fields clusters. We discover four highly magnified $z\sim7$ candidate Lyman-break galaxies behind the cluster, one of which may be multiply-imaged. Correcting for magnification, we find that all four galaxies are fainter than $0.5 L_{\star}$. One candidate is detected at ${>10\sigma}$ in both Spitzer/IRAC [3.6] and [4.5] channels. A spectroscopic follow-up with MOSFIRE does not result in the detection of the Lyman-alpha emission line from any of the four candidates. From the MOSFIRE spectra we place median upper limits on the Lyman-alpha flux of $5-14 \times 10^{-19}\, \mathrm{erg \,\, s^{-1} cm^{-2}}$ ($5\sigma$).
Sulphur is an important, volatile alpha element but its role in the Galactic chemical evolution is still uncertain. We derive the S abundances in RGB stars in three Galactic globular clusters (GC) that cover a wide metallicity range (-2.3<[Fe/H]<-1.2), namely M4, M22, and M30. The halo field stars show a large scatter in the [S/Fe] ratio in this metallicity span, which is inconsistent with canonical chemical evolution models. To date, very few measurements of [S/Fe] exist for stars in GCs, which are good tracers of the chemical enrichment of their environment. However, some light and alpha elements show star-to-star variations within individual GCs and it is yet unclear whether sulphur also varies between GC stars. We used the the infrared spectrograph CRIRES to obtain high-resolution (R~50000), high signal-to-noise (SNR~200 per px) spectra in the region of the S I multiplet 3 at 1045 nm for 15 GC stars selected from the literature (6 stars in M4, 6 stars in M22 and 3 stars in M30). Multiplet 3 is better suited for S abundance derivation than the more commonly used lines of multiplet 1 at 920 nm, since its lines are not blended by telluric absorption or other stellar features at low metallicity. We used spectral synthesis to derive the [S/Fe] ratio of the stars assuming local thermodynamic equilibrium (LTE). We find mean [S/Fe] = 0.58 +/- 0.01 +/- 0.20 dex (statistical and systematic error) for M4, [S/Fe] = 0.57+/-0.01+/-0.19 dex for M22, and [S/Fe] = 0.55+/-0.02+/-0.16 dex for M30. The negative NLTE corrections are estimated to be in the order of the systematic uncertainties. With the tentative exception of two stars with measured high S abundances, we conclude that sulphur behaves like a typical alpha element in the studied Galactic GCs, showing enhanced abundances with respect to the solar value at metallicities below [Fe/H] = -1.0 dex without a considerable spread.
We examine the performance of four different methods which are used to measure mass segregation in star-forming regions: the radial variation of the mass function $\mathcal{M}_{\rm MF}$; the minimum spanning tree-based $\Lambda_{\rm MSR}$ method; the local surface density $\Sigma_{\rm LDR}$ method; and the $\Omega_{\rm GSR}$ technique, which isolates groups of stars and determines whether the most massive star in each group is more centrally concentrated than the average star. All four methods have been proposed in the literature as techniques for quantifying mass segregation, yet they routinely produce contradictory results as they do not all measure the same thing. We apply each method to synthetic star-forming regions to determine when and why they have shortcomings. When a star-forming region is smooth and centrally concentrated, all four methods correctly identify mass segregation when it is present. However, if the region is spatially substructured, the $\Omega_{\rm GSR}$ method fails because it arbitrarily defines groups in the hierarchical distribution, and usually discards positional information for many of the most massive stars in the region. We also show that the $\Lambda_{\rm MSR}$ and $\Sigma_{\rm LDR}$ methods can sometimes produce apparently contradictory results, because they use different definitions of mass segregation. We conclude that only $\Lambda_{\rm MSR}$ measures mass segregation in the classical sense (without the need for defining the centre of the region), although $\Sigma_{\rm LDR}$ does place limits on the amount of previous dynamical evolution in a star-forming region.
We present the J and H-band source catalog covering the AKARI North Ecliptic Pole field. Filling the gap between the optical data from other follow-up observations and mid-infrared (MIR) data from AKARI, our near-infrared (NIR) data provides contiguous wavelength coverage from optical to MIR. For the J and H-band imaging, we used the FLoridA Multi-object Imaging Near-ir Grism Observational Spectrometer (FLAMINGOS) on the Kitt Peak National Observatory 2.1m telescope covering a 5.1 deg2 area down to a 5 sigma depth of ~21.6 mag and ~21.3 mag (AB) for J and H-band with an astrometric accuracy of 0.14" and 0.17" for 1 sigma in R.A. and Decl. directions, respectively. We detected 208,020 sources for J-band and 203,832 sources for H-band. This NIR data is being used for studies including analysis of the physical properties of infrared sources such as stellar mass and photometric redshifts, and will be a valuable dataset for various future missions.
Submillimetre (submm) observations of WISE-selected, dusty, luminous, high-redshift galaxies have revealed intriguing overdensities around them on arcmin scales. They could be the best signposts of overdense environments on the sky.
We estimate the size and distribution of the parent populations for the 6 largest (at least 20 stars in the Solar neighborhood) chemical groups identified in the Chemical Tagging experiment by Mitschang et al.~2014. Stars in the abundance groups tend to lie near a boundary in angular momentum versus eccentricity space where the probability is highest for a star to be found in the Solar neighborhood and where orbits have apocenter approximately equal to the Sun's galactocentric radius. Assuming that the parent populations are uniformly distributed at all azimuthal angles in the Galaxy, we estimate that the parent populations of these abundance groups contain at least 200,000 members. The spread in angular momentum of the groups implies that the assumption of a uniform azimuthal distribution only fails for the two youngest groups and only for the highest angular momentum stars in them. The parent populations of three thin disk groups have narrow angular momentum distributions, but tails in the eccentricity and angular momentum distributions suggest that only a small fraction of stars have migrated and increased in eccentricity. In contrast, the parent populations of the thick disk groups exhibit both wide angular momentum and eccentricity distributions implying that both heating and radial migration has taken place.
Aims. We explore the problem of the site of production of Eu. We use also the information present in the observed spread in the Eu abundances in the early Galaxy, not only its average trend. Moreover, we extend to other heavy elements (Ba, Sr, Rb, Zr) our investigations to provide additional constraints to our results. Methods. We adopt a stochastic chemical evolution model taking into account inhomogeneous mixing. The adopted yields of Eu from neutron star mergers (NSM) and from core-collapse supernovae (SNII) are those that are able to explain the average [Eu/Fe]-[Fe/H] trend observed for solar neighborhood stars, in the framework of a well-tested homogeneous model for the chemical evolution of the MilkyWay. Rb, Sr, Zr, and Ba are produced by both the s- and r-process. The s-process contribution by spinstars is the same as in our previous papers. Results. NSM that merge in less than 10 Myr or NSM combined with a source of r-process generated by massive stars can explain the spread of [Eu/Fe] in the Galactic halo. The combination of r-process production by NSM and s-process production by spinstars is able to reproduce the available observational data for Sr, Zr and Ba. We also show the first predictions for Rb in the Galactic halo. Conclusions. We confirm previous results that either NSM with very short time scale or both NSM and at least a fraction of SNII should have contributed to the synthesis of Eu in the Galaxy. The r-process production by NSM - complemented by an s-process production by spinstars - provide results compatible with our previous findings based on other r-process sites. We critically discuss the weak and strong points of both NSM and SNII scenarios for producing Eu and eventually suggest that the best solution is probably a mixed one in which both sources produce Eu. In fact, this scenario better reproduces the scatter observed in all the studied elements. [abridged]
The whole Sloan Digital Sky Survey (SDSS, 14,555 deg^2 has been searched for intergalactic globular clusters (IGCs) in the Local Group (LG). Using optical, infrared, and ultraviolet photometric selection criteria and photometric redshifts, the 2.1x10^8 of objects in the SDSS Galaxy Catalogue were reduced to only 183,791 brighter than r_o = 19 that might be GCs. Visual examination of their SDSS images recovered 84 percent of the confirmed GCs in M31 and M33 and yielded 17 new GC candidates, 5 of them of high confidence, which we could confirm as GCs in MegaPrime images from the Canada, France, Hawaii Telescope. These 5 GCs are within M31's halo, but the other 12 candidates are not close to LG galaxies or galaxies within 3 Mpc of the LG. Even though this search covers only one-third of the sky and some GCs could have been missed, it suggests that the LG does not contain a large population of IGCs more luminous than Mv ~ -6. In the direction of the M81 Group, the search yielded five candidate GCs, probable members of that group.
The 9.7-micron silicate absorption profile in the interstellar medium provides important information on the physical and chemical composition of interstellar dust grains. Measurements in the Milky Way have shown that the profile in the diffuse interstellar medium is very similar to the amorphous silicate profiles found in circumstellar dust shells around late M stars, and narrower than the silicate profile in denser star-forming regions. Here, we investigate the silicate absorption profile towards the very heavily obscured nucleus of NGC 4418, the galaxy with the deepest known silicate absorption feature, and compare it to the profiles seen in the Milky Way. Comparison between the 8-13 micron spectrum obtained with TReCS on Gemini and the larger aperture spectrum obtained from the Spitzer archive indicates that the former isolates the nuclear emission, while Spitzer detects low surface brightness circumnuclear diffuse emission in addition. The silicate absorption profile towards the nucleus is very similar to that in the diffuse ISM in the Milky Way with no evidence of spectral structure from crystalline silicates or silicon carbide grains.
The high quality light curves of Kepler space telescope make it possible to analyze the optical variability of AGNs with an unprecedented time resolution. Studying the asymmetry in variations could give independent constraints on the physical models for AGN variability. In this paper, we use Kepler observations of 19 sources to perform analyses on the variability asymmetry of AGNs. We apply smoothing-correction to light curves to deduct the bias to high frequency variability asymmetry, caused by long term variations which are poorly sampled due to the limited length of light curves. A parameter $\beta$ based on structure functions is introduced to quantitively describe the asymmetry and its uncertainty is measured using extensive Monte-Carlo simulations. Individual sources show no evidence of asymmetry at timescales of $1\sim20$ days and there is not a general trend toward positive or negative asymmetry over the whole sample. Stacking data of all 19 AGNs, we derive averaged $\overline{\beta}$ of 0.00$\pm$0.03 and -0.02$\pm$0.04 over timescales of 1$\sim$5 days and 5$\sim$20 days, respectively, statistically consistent with zero. Quasars and Seyfert galaxies show similar asymmetry parameters. Our results indicate that short term optical variations in AGNs are highly symmetric.
We study the properties of gas in and around 10^12 solar mass halos at z=2 using a suite of high-resolution cosmological hydrodynamic 'zoom' simulations. We quantify the thermal and dynamical structure of these gaseous reservoirs in terms of their mean radial distributions and angular variability along different sightlines. With each halo simulated at three levels of increasing resolution, the highest reaching a baryon mass resolution of ~10,000 solar masses, we study the interaction of filamentary inflow and the quasi-static hot halo atmosphere. We highlight the discrepancy between the spatial resolution available in the halo gas as opposed to within the galaxy itself, and find that stream morphologies become increasingly complex at higher resolution, with large coherent flows revealing density and temperature structure at progressively smaller scales. Moreover, multiple gas components co-exist at the same radius within the halo, making radially averaged analyses misleading. This is particularly true where the hot, quasi-static, high entropy halo atmosphere interacts with cold, rapidly inflowing, low entropy accretion. We investigate the process of gas virialization and identify different regimes for the heating of gas as it accretes from the intergalactic medium. Haloes at this mass have a well-defined virial shock, associated with a sharp jump in temperature and entropy at ~1.25 r_vir. The presence, radius, and radial width of this boundary feature, however, vary not only from halo to halo, but also as a function of angular direction, covering roughly ~85% of the 4pi sphere. Our findings are relevant for the proper interpretation of observations pertaining to the circumgalactic medium, including evidence for large amounts of cold gas surrounding massive haloes at intermediate redshifts.
We present a new symplectic numerical integrator designed for collisional gravitational $N$-body problems which makes use of Kepler solvers. The integrator is also reversible and conserves 9 integrals of motion of the $N$-body problem to machine precision. The integrator is second order, but the order can easily be increased by the method of \citeauthor{yos90}. We use fixed time step in all tests studied in this paper to ensure preservation of symplecticity. We study small $N$ collisional problems and perform comparisons with typically used integrators. In particular, we find comparable or better performance when compared to the 4th order Hermite method and much better performance than adaptive time step symplectic integrators introduced previously. The integrator is a promising tool in collisional gravitational dynamics. We plan larger $N$ tests of the method in future work.
A sample of 46 nearby clusters observed with Chandra is analyzed to produce radial density, temperature, entropy and metallicity profiles, as well as other morphological measurements. The entropy profiles are computed to larger radial extents than in previous Chandra cluster sample analyses. We find that the iron mass fraction measured in the inner 0.15 R500 shows a larger dispersion across the sample of low-mass clusters, than it does for the sample of high-mass clusters. We interpret this finding as the result of the mixing of more haloes in large clusters than in small clusters, which leads to an averaging of the metal content in the large clusters, and thus less dispersion of metallicity for high-mass clusters. This interpretation lends support to the idea that the low-entropy, metal-rich gas of merging haloes reaches clusters' centers, which explains observations of Core-Collapse Supernova products metallicity peaks, and which is seen in hydrodynamical simulations. The gas in these merging haloes would have to reach the centers of clusters without mixing in the outer regions, in order to support our interpretation. On the other hand, metallicity dispersion does not change with mass in the outer regions of clusters, suggesting that most of the outer metals come from a source with a more uniform metallicity level, such as during pre-enrichment. We also measure a correlation between the metal content in low-mass clusters and the degree to which their Intra-Cluster Medium (ICM) is morphologically disturbed, as measured by centroid shift. This suggests an alternative interpretation of the large width of the metallicity distribution in low-mass clusters, whereby a metallicity boost in the center of low-mass clusters is induced as a transitional state, during mergers.
Selection of extragalactic point sources, e.g. QSOs, is often hampered by significant selection effects causing existing samples to have rather complex selection functions. We explore whether a purely astrometric selection of extragalactic point sources, e.g. QSOs, is feasible with the ongoing Gaia mission. Such a selection would be interesting as it would be unbiased in terms of colours of the targets and hence would allow selection also with colours in the stellar sequence. We have analyzed a total of 18 representative regions of the sky by using \textit{GUMS}, the simulator prepared for ESAs Gaia mission, both in the range of $12\le G \le 20$ mag and $12\le G \le 18$ mag. For each region we determine the density of apparently stationary stellar sources, i.e. sources for which Gaia cannot measure a significant proper motion. The density is contrasted with the density of extragalactic point sources, e.g. QSOs, in order to establish in which celestial directions a pure astrometric selection is feasible. When targeting regions at galactic latitude $|b| \ge 30^\mathrm{o}$ the ratio of QSOs to apparently stationary stars is above 50\% and when observing towards the poles the fraction of QSOs goes up to about $\sim80$\%. We show that the proper motions from the proposed Gaia successor mission in about 20 years would dramatically improve these results at all latitudes. Detection of QSOs solely from zero proper motion, unbiased by any assumptions on spectra, might lead to the discovery of new types of QSOs or new classes of extragalactic point sources.
We present a multi-wavelength observational study of a low-mass star-forming region, L1251-C, with observational results at wavelengths from the near-infrared to the millimeter. Spitzer Space Telescope observations confirmed that IRAS 22343+7501 is a small group of protostellar objects. The extended emission to east-west direction with its intensity peak at the center of L1251A has been detected at 350 and 850 micron with the CSO and JCMT telescopes, tracing dense envelope materials around L1251A. The single-dish data from the KVN and TRAO telescopes show inconsistencies between the intensity peaks of several molecular line emission and that of the continuum emission, suggesting complex distributions of molecular abundances around L1251A. The SMA interferometer data, however, show intensity peaks of CO 2-1 and 13CO 2-1 located at the position of IRS 1, which is both the brightest source in IRAC image and the weakest source in the 1.3 mm dust continuum map. IRS 1 is the strongest candidate for the driving source of the newly detected compact CO 2-1 outflow. Over the whole region (14' by 14') of L125l-C, 3 Class I and 16 Class II sources have been detected, including three YSOs in L1251A. A comparison with the average projected distance among 19 YSOs in L1251-C and that among 3 YSOs in L1251A suggests L1251-C is an example of low-mass cluster formation, where protostellar objects are forming in a small group.
Sufficiently massive clumps of molecular gas collapse under self-gravity and fragment to spawn a cluster of stars that have a range of masses. We investigate observationally the early stages of formation of a stellar cluster in a massive filamentary infrared dark cloud, G28.34+0.06 P1, in the 1.3mm continuum and spectral line emission using the ALMA. Sensitive continuum data reveal further fragmentation in five dusty cores at a resolution of several 10^3 AU. Spectral line emission from C18O, CH3OH, 13CS, H2CO and N2D+ are detected for the first time toward these dense cores. We found that three cores are chemically more evolved as compared with the other two; interestingly though, all of them are associated with collimated outflows as suggested by evidence from the CO, SiO, CH3OH, H2CO and SO emissions. The parsec-scale kinematics in NH3 exhibit velocity gradients along the filament, consistent with accretion flows toward the clumps and cores. The moderate luminosity and the chemical signatures indicate that the five cores harbor low- to intermediate-mass protostars that likely become massive ones at the end of the accretion. Despite the fact that the mass limit reached by the 1\sigma dust continuum sensitivity is 30 times lower than the thermal Jeans mass, there is a lack of a distributed low-mass protostellar population in the clump. Our observations indicate that in a protocluster, low-mass stars form at a later stage after the birth of more massive protostars.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
We examine the effects that the modelling of a Boxy/Peanut (B/P) bulge will have on the estimates of the stellar gravitational potential, forces, orbital structure and bar strength of barred galaxies. We present a method for obtaining the potential of disc galaxies from surface density images, assuming a vertical density distribution (height function), which is let to vary with position, thus enabling it to represent the geometry of a B/P. We construct a B/P height function after the results from a high-resolution, N-body+SPH simulation of an isolated galaxy and compare the resulting dynamical model to those obtained with the commonly used, position-independent "flat" height functions. We show that methods that do not allow for a B/P can induce errors in the forces in the bar region of up to 40% and demonstrate that this has a significant impact on the orbital structure of the model, which in turn determines its kinematics and morphology. Furthermore, we show that the bar strength is reduced in the presence of a B/P. We conclude that neglecting the vertical extent of a B/P can introduce considerable errors in the dynamical modelling. We also examine the errors introduced in the model due to uncertainties in the parameters of the B/P and show that even for generous but realistic values of the uncertainties, the error will be noticeably less than that of not modelling a B/P bulge at all.
Spiral structure is the most distinctive feature of disk galaxies and yet debate persists about which theory of spiral structure is the correct one. Many versions of the density wave theory demand that the pitch angle is uniquely determined by the distribution of mass in the bulge and disk of the galaxy. We present evidence that the tangent of the pitch angle of logarithmic spiral arms in disk galaxies correlates strongly with the density of neutral atomic hydrogen in the disk and with the central stellar bulge mass of the galaxy. These three quantities, when plotted against each other, form a planar relationship which, we argue should be fundamental to our understanding of spiral structure in disk galaxies. We further argue that any successful theory of spiral structure must be able to explain this relationship.
A number of stellar sources have been advocated as the origin of the enriched material required to explain the abundance anomalies seen in ancient globular clusters (GCs). Most studies to date have compared the yields from potential sources (asymptotic giant branch stars (AGBs), fast rotating massive stars (FRMS), high mass interacting binaries (IBs), and very massive stars (VMS)) with observations of specific elements that are observed to vary from star-to-star in GCs, focussing on extreme GCs such as NGC 2808, which display large He variations. However, a consistency check between the results of fitting extreme cases with the requirements of more typical clusters, has rarely been done. Such a check is particularly timely given the constraints on He abundances in GCs now available. Here we show that all of the popular enrichment sources fail to reproduce the observed trends in GCs, focussing primarily on Na, O and He. In particular, we show that any model that can fit clusters like NGC 2808, will necessarily fail (by construction) to fit more typical clusters like 47 Tuc or NGC 288. All sources severely over-produce He for most clusters. Additionally, given the large differences in He spreads between clusters, but similar spreads observed in Na--O, only sources with large degrees of stochasticity in the resulting yields will be able to fit the observations. We conclude that no enrichment source put forward so far (AGBs, FRMS, IBs, VMS - or combinations thereof) is consistent with the observations of GCs. Finally, the observed trends of increasing [N/Fe] and He spread with increasing cluster mass cannot be resolved within a self-enrichment framework, without further exacerbating the mass budget problem.
We present the abundance analysis of 82 red giant branch stars in the dense, metal-poor globular cluster NGC 6093 (M 80), the largest sample of stars analyzed in this way for this cluster. From high resolution UVES spectra of 14 stars and intermediate resolution GIRAFFE spectra for the other stars we derived abundances of O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Ba, La, Ce, Pr, Nd, Sm, Eu. On our UVES metallicity scale the mean metal abundance of M 80 is [Fe/H]=-1.791+/-0.006+/-0.076 (+/-statistical +/-systematic error) with rms=0.023 (14 stars). M 80 shows star to star variations in proton-capture elements, and the extension of the Na-O anticorrelation perfectly fit the relations with (i) total cluster mass, (ii) horizontal branch morphology, and (iii) cluster concentration previously found by our group. The chemistry of multiple stellar populations in M 80 does not look extreme. The cluster is also a typical representative of halo globular clusters for what concerns the pattern of alpha-capture and Fe-group elements. However we found that a significant contribution from the s-process is required to account for the distribution of neutron-capture elements. A minority of stars in M 80 seem to exhibit slightly enhanced abundances of s-process species, compatible with those observed in M 22 and NGC 1851, although further confirmation from larger samples is required.
We report the discovery of rapid variations of a high-velocity CIV broad absorption line trough in the quasar SDSS J141007.74+541203.3. This object was intensively observed in 2014 as a part of the Sloan Digital Sky Survey Reverberation Mapping Project, during which 32 epochs of spectroscopy were obtained with the Baryon Oscillation Spectroscopic Survey spectrograph. We observe significant (>4sigma) variability in the equivalent width of the broad (~4000 km/s wide) CIV trough on rest-frame timescales as short as 1.20 days (~29 hours), the shortest broad absorption line variability timescale yet reported. The equivalent width varied by ~10% on these short timescales, and by about a factor of two over the duration of the campaign. We evaluate several potential causes of the variability, concluding that the most likely cause is a rapid response to changes in the incident ionizing continuum. If the outflow is at a radius where the recombination rate is higher than the ionization rate, the timescale of variability places a lower limit on the density of the absorbing gas of n_e > 3.9 x 10^5 cm^-3. The broad absorption line variability characteristics of this quasar are consistent with those observed in previous studies of quasars, indicating that such short-term variability may in fact be common and thus can be used to learn about outflow characteristics and contributions to quasar/host-galaxy feedback scenarios.
To provide a quantitative cosmological context to ongoing observational work on the formation histories and location of compact massive galaxies, we locate and study a sample of exceptionally compact systems in the Bolshoi simulation, using the dark matter structural parameters from a real, compact massive galaxy (NGC1277) as a basis for our working criteria. We find that over 80% of objects in this nominal compact category are substructures of more massive groups or clusters, and that the probability of a given massive substructure being this compact increases significantly with the mass of the host structure; rising to ~30% for the most massive clusters in the simulation. Tracking the main progenitors of this subsample back to z=2, we find them all to be distinct structures with scale radii and densities representative of the population as a whole at this epoch. What does characterise their histories, in addition to mostly becoming substructures, is that they have almost all experienced below-average mass accretion since z=2; a third of them barely retaining, or even losing mass during the intervening 10 Gyr.
We explore the connection between the black hole mass and its relativistic jet for a sample of radio-loud AGN (z < 1), in which the relativistic jet parameters are well estimated by means of long term monitoring with the 14m Mets\"ahovi millimeter wave telescope and the Very Long Base-line Array (VLBA). NIR host galaxy images taken with the NOTCam on the Nordic Optical Telescope (NOT) and retrieved from the 2MASS all-sky survey allowed us to perform a detailed surface brightness decomposition of the host galaxies in our sample and to estimate reliable black hole masses via their bulge luminosities. We present early results on the correlations between black hole mass and the relativistic jet parameters. Our preliminary results suggest that the more massive the black hole is, the faster and the more luminous jet it produces.
Our work is based on the "Bluedisk" project, a program to map the neutral gas in a sample of 25 HI-rich spirals and a similar number of control galaxies with the Westerbork Synthesis Radio Telescope (WSRT). In this paper we focus on the HI properties of the galaxies in the environment of our targeted galaxies. In total, we extract 65 galaxies from the WSRT cubes with stellar masses between $10^8M_{\odot}$ and $10^{11}M_{\odot}$. Most of these galaxies are located on the same HI mass-size relation and "HI-plane" as normal spiral galaxies. We find that companions around HI-rich galaxies tend to be HI-rich as well and to have larger R90,HI/R50,HI. This suggests a scenario of "HI conformity", similar to the colour conformity found by Weinmann et al. (2006): galaxies tend to adopt the HI properties of their neighbours. We visually inspect the outliers from the HI mass-size relation and galaxies which are offset from the HI plane and find that they show morphological and kinematical signatures of recent interactions with their environment. We speculate that these outliers have been disturbed by tidal or ram-pressure stripping processes, or in a few cases, by accretion events.
We present a survey of star clusters in the halo of IC 10, a starburst galaxy in the Local Group based on Subaru R band images and NOAO Local Group Survey UBVRI images. We find five new star clusters. All these star clusters are located far from the center of IC 10, while previously known star clusters are mostly in the main body. Interestingly the distribution of these star clusters shows an asymmetrical structure elongated along the east and south-west direction. We derive UBVRI photometry of 66 star clusters including these new star clusters as well as previously known star clusters. Ages of the star clusters are estimated from the comparison of their UBVRI spectral energy distribution with the simple stellar population models. We find that the star clusters in the halo are all older than 1 Gyr, while those in the main body have various ages from very young (several Myr) to old (>1 Gyr). The young clusters (<10 Myr) are mostly located in the H{\alpha} emission regions and are concentrated on a small region at 2' in the south-east direction from the galaxy center, while the old clusters are distributed in a wider area than the disk. Intermediate-age clusters (~100 Myr) are found in two groups. One is close to the location of the young clusters and the other is at ~4' from the location of the young clusters. The latter may be related with past merger or tidal interaction.
The presence of a helical magnetic field threading the jet of an Active Galactic Nucleus (AGN) should give rise to a gradient in the observed Faraday rotation measure (RM) across the jet, due to the associated systematic change in the line-of-sight magnetic field. Reports of observations of transverse RM gradients across AGN jets have appeared in the literature starting from 2002, but concerns were raised about the resolution required for these gradients to be reliable, and there was a lack of a full understanding of the best approach to accurate estimation of the uncertainties of local RM values. These questions have now been resolved by recent Monte Carlo simulations carried out by various groups, enabling both a verification of previously published results and reliable analyses of new data. We consider here RM gradients across the jet structures of 15 AGN, some previously published in the refereed literature but without a correct and complete error analysis, and some published for the first time here, all of which have monotonic transverse RM gradients with significances of at least 3 sigma.
ALMA is revolutionizing the way we study and understand the astrophysics of galaxies, both as a whole and individually. By exploiting its unique sensitivity and resolution to make spatially and spectrally resolved images of the gas and dust in the interstellar medium (ISM), ALMA can reveal new information about the relationship between stars and gas, during and between galaxies' cycles of star formation and AGN fueling. However, this can only be done for a modest number of targets, and thus works in the context of large samples drawn from other surveys, while providing parallel deep imaging in small fields around. Recent ALMA highlights are reviewed, and some areas where ALMA will potentially make great contributions in future are discussed.
We have observed 15 red giant stars in the relatively massive, metal-poor globular cluster NGC 4833 using the MIKE spectrograph at Magellan. We calculate stellar parameters for each star and perform a standard abundance analysis to derive abundances of 43 species of 39 elements, including 20 elements heavier than the iron group. We derive <[Fe/H]> = -2.25 +/- 0.02 from Fe I lines and <[Fe/H> = -2.19 +/- 0.013 from Fe II lines. We confirm earlier results that found no internal metallicity spread in NGC 4833, and there are no significant star-to-star abundance dispersions among any elements in the iron group (19 <= Z <= 30). We recover the usual abundance variations among the light elements C, N, O, Na, Mg, Al, and possibly Si. The heavy-element distribution reflects enrichment by r-process nucleosynthesis ([Eu/Fe] = +0.36 +/- 0.03), as found in many other metal-poor globular clusters. We investigate small star-to-star variations found among the neutron-capture elements, and we conclude that these are probably not real variations. Upper limits on the Th abundance, log epsilon (Th/Eu) < -0.47 +/- 0.09, indicate that NGC 4833, like other globular clusters where Th has been studied, did not experience a so-called "actinide boost."
A galaxy group catalog is built from the sample of the 2MASS Redshift Survey almost complete to Ks=11.75 over 91% of the sky. Constraints in the construction of the groups were provided by scaling relations determined by close examination of well defined groups with masses between 10^11 and 10^15 Msun. Group masses inferred from Ks luminosities are statistically in agreement with masses calculated from application of the virial theorem. While groups have been identified over the full redshift range of the sample, the properties of the nearest and farthest groups are uncertain and subsequent analysis has only considered groups with velocities between 3,000 and 10,000 km/s. The 24,044 galaxies in this range are identified with 13,607 entities, 3,461 of them with two or more members. A group mass function is constructed. The Sheth-Tormen formalism provides a good fit to the shape of the mass function for group masses above 6/h x 10^12 Msun but the count normalization is poor. Summing all the mass associated with the galaxy groups between 3,000 and 10,000 km/s gives a density of collapsed matter as a fraction of the critical density of Omega_collapsed = 0.16.
The sensitivity of pulsar timing arrays to gravitational waves is, at some level, limited by timing noise. Red timing noise - the stochastic wandering of pulse arrival times with a red spectrum - is prevalent in slow-spinning pulsars and has been identified in many millisecond pulsars. Phenomenological models of timing noise, such as from superfluid turbulence, suggest that the timing noise spectrum plateaus below some critical frequency, $f_c$, potentially aiding the hunt for gravitational waves. We examine this effect for individual pulsars by calculating minimum observation times, $T_{\rm min}(f_c)$, over which the gravitational wave signal becomes larger than the timing noise plateau. We do this in two ways: 1) in a model-independent manner, and 2) by using the superfluid turbulence model for timing noise as an example to illustrate how neutron star parameters can be constrained. We show that the superfluid turbulence model can reproduce the data qualitatively from a number of pulsars observed as part of the Parkes Pulsar Timing Array. We further show how a value of $f_c$, derived either through observations or theory, can be related to $T_{\rm min}$. This provides a diagnostic whereby the usefulness of timing array pulsars for gravitational-wave detection can be quantified.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
We demonstrate that high abundances of water vapor could have existed in extremely low metallicity ($10^{-3}$ solar) partially shielded gas, during the epoch of first metal enrichment of the interstellar medium of galaxies at high redshifts.
We present a new potential-density pair designed to model nearly isothermal star clusters (and similar self-gravitating systems) with a central core and an outer turnover radius, beyond which density falls off as $r^{-4}$. In the intermediate zone, the profile is similar to that of an isothermal sphere (density $\rho \propto r^{-2}$), somewhat less steep than the King 1962 profile, and with the advantage that many dynamical quantities can be written in a simple closed form. We derive analytic expressions for the cluster binding energy, central velocity dispersion, and escape velocity, and apply these to create toy models for cluster core collapse and evaporation. We rederive classical results for evaporating, collapsing, and quasi-equilibrium (heated) clusters, and fit our projected surface brightness profiles to observed globular and open clusters. We find that the quality of the fit is generally at least as good as that for the surface brightness profiles of King 1962. This model can be used for convenient computation of the dynamics and evolution of globular and nuclear star clusters.
The evolution and distribution of the angular momentum of dark matter halos have been discussed in several studies over the last decades. In particular, the idea arose that angular momentum conservation should allow to infer the total angular momentum of the entire dark matter halo from measuring the angular momentum of the baryonic component, which is populating the center of the halo, especially for disk galaxies. To test this idea and to understand the connection between the angular momentum of the dark matter halo and its galaxy, we use the Magneticum Simulations. We successfully produce populations of spheroidal and disk galaxies self-consistently. Thus, we are able to study the dependance of galactic properties on their morphology. We find that: (I) The specific angular momentum of stars in disk and spheroidal galaxies as function of their stellar mass compares well with observational results; (II) The specific angular momentum of the stars in disk galaxies is slightly smaller compared to the specific angular momentum of the cold gas, in good agreement with observations; (III) Simulations including the baryonic component show a dichotomy in the specific stellar angular momentum distribution when splitting the galaxies according to their morphological type. This dichotomy can also be seen in the spin parameter, where disk galaxies populate halos with slightly larger spin compared to spheroidal galaxies; (IV) Disk galaxies preferentially populate halos in which the angular momentum vector of the dark matter component in the central part shows a better alignment to the angular momentum vector of the entire halo; (V) The specific angular momentum of the cold gas in disk galaxies is approximately 40 percent smaller than the specific angular momentum of the total dark matter halo and shows a significant scatter.
We have carried out 13CO (J=2-1) observations of the active star-forming region N159 West in the LMC with ALMA. We have found that the CO distribution at a sub-pc scale is highly elongated with a small width. These elongated clouds called "filaments" show straight or curved distributions with a typical width of 0.5-1.0 pc and a length of 5-10 pc. All the known infrared YSOs are located toward the filaments. We have found broad CO wings of two molecular outflows toward young high-mass stars in N159W-N and N159W-S, whose dynamical timescale is ~10^4 yrs. This is the first discovery of protostellar outflow in external galaxies. For N159W-S which is located toward an intersection of two filaments we set up a hypothesis that the two filaments collided with each other ~10^5 yrs ago and triggered formation of the high-mass star having ~37Mo. The colliding clouds show significant enhancement in linewidth in the intersection, suggesting excitation of turbulence in the shocked interface layer between them as is consistent with the magneto-hydro-dynamical numerical simulations (Inoue & Fukui 2013). This turbulence increases the mass accretion rate to ~10^-4 Mo yr^-1, which is required to overcome the stellar feedback to form the high-mass star.
A deep, wide-field, near-infrared imaging survey was used to construct an extinction map of the southeastern part of the California Molecular Cloud (CMC) with $\sim$ 0.5 arc min resolution. The same region was also surveyed in the $^{12}$CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1) emission lines at the same angular resolution. Strong spatial variations in the abundances of $^{13}$CO and C$^{18}$O were found to be correlated with variations in gas temperature, consistent with temperature dependent CO depletion/desorption on dust grains. The $^{13}$CO to C$^{18}$O abundance ratio was found to increase with decreasing extinction, suggesting selective photodissociation of C$^{18}$O by the ambient UV radiation field. The cloud averaged X-factor is found to be $<$X$_{\rm CO}$$>$ $=$ 2.53 $\times$ 10$^{20}$ ${\rm cm}^{-2}~({\rm K~km~s}^{-1})^{-1}$, somewhat higher than the Milky Way average. On sub-parsec scales we find no single empirical value of the X-factor that can characterize the molecular gas in cold (T$_{\rm k}$ $\lesssim$ 15 K) regions, with X$_{\rm CO}$ $\propto$ A$_{\rm V}$$^{0.74}$ for A$_{\rm V}$ $\gtrsim$ 3 magnitudes. However in regions containing relatively hot (T$_{\rm ex}$ $\gtrsim$ 25 K) gas we find a clear correlation between W($^{12}$CO) and A$_{\rm V}$ over a large (3 $\lesssim$ A$_{\rm V}$ $\lesssim$ 25 mag) extinction range. This suggests a constant X$_{\rm CO}$ $=$ 1.5 $\times$ 10$^{20}$ ${\rm cm}^{-2}~({\rm K~km~s}^{-1})^{-1}$ for the hot gas, a lower value than either the average for the CMC or Milky Way. We find a correlation between X$_{\rm CO}$ and T$_{\rm ex}$ with X$_{\rm CO}$ $\propto$ T$_{\rm ex}$$^{-0.7}$ suggesting that the global X-factor of a cloud may depend on the relative amounts of hot gas within it.
This is the third in a series of papers reporting on a large reverberation-mapping campaign aimed to study the properties of active galactic nuclei (AGNs) with high accretion rates. We present new results on the variability of the optical Fe II emission lines in 10 AGNs observed by the Yunnan Observatory 2.4m telescope during 2012--2013. We detect statistically significant time lags, relative to the AGN continuum, in nine of the sources. This accurate measurement is achieved by using a sophisticated spectral fitting scheme that allows for apparent flux variations of the host galaxy, and several narrow lines, due to the changing observing conditions. Six of the newly detected lags are indistinguishable from the Hbeta lags measured in the same sources. Two are significantly longer and one is slightly shorter. Combining with Fe II lags reported in previous studies, we find a Fe II radius--luminosity relationship similar to the one for Hbeta, although our sample by itself shows no clear correlation. The results support the idea that Fe II emission lines originate in photoionized gas which, for the majority of the newly reported objects, is indistinguishable from the Hbeta-emitting gas. We also present a tentative correlation between the lag and intensity of Fe II and Hbeta and comment on its possible origin.
A family of spherical halo models with flat circular velocity curves is presented. This includes models in which the rotation curve has a finite central value but declines outwards (like the Jaffe model). It includes models in which the rotation curve is rising in the inner parts, but flattens asymptotically (like the Binney model). The family encompasses models with both finite and singular (cuspy) density profiles. The self-consistent distribution function depending on binding energy $E$ and angular momentum $L$ is derived and the kinematical properties of the models discussed. These really describe the properties of the total matter (both luminous and dark). For comparison with observations, it is better to consider tracer populations of stars. These can be used to represent elliptical galaxies or the spheroidal components of spiral galaxies. Accordingly, we study the properties of tracers with power-law or Einasto profiles moving in the doubloon potential. Under the assumption of spherical alignment, we provide a simple way to solve the Jeans equations for the velocity dispersions. This choice of alignment is supported by observations on the stellar halo of the Milky Way. Power-law tracers have prolate spheroidal velocity ellipsoids everywhere. However, this is not the case for Einasto tracers, for which the velocity ellipsoids change from prolate to oblate spheroidal near the pole. Asymptotic forms of the velocity distributions close to the escape speed are also derived, with an eye to application to the high velocity stars in the Milky Way. Power-law tracers have power-law or Maxwellian velocity distributions tails, whereas Einasto tracers have super-exponential cut-offs.
[Abridged] Fullerenes have been recently detected in various circumstellar and interstellar environments, raising the question of their formation pathway. It has been proposed that they can form by the photo-chemical processing of large polycyclic aromatic hydrocarbons (PAHs). Following our previous work on the evolution of PAHs in the NGC 7023 reflection nebula, we evaluate, using photochemical modeling, the possibility that the PAH C$_{66}$H$_{20}$ (i.e. circumovalene) can lead to the formation of C$_{60}$ upon irradiation by ultraviolet photons. The chemical pathway involves full dehydrogenation, folding into a floppy closed cage and shrinking of the cage by loss of C$_2$ units until it reaches the symmetric C$_{60}$ molecule. At 10" from the illuminating star and with realistic molecular parameters, the model predicts that 100\% of C$_{66}$H$_{20}$ is converted into C$_{60}$ in $\sim$ 10$^5$ years, a timescale comparable to the age of the nebula. Shrinking appears to be the kinetically limiting step of the whole process. Hence, PAHs larger than C$_{66}$H$_{20}$ are unlikely to contribute significantly to the formation of C$_{60}$, while PAHs containing between 60 and 66 C atoms should contribute to the formation of C$_{60}$ with shorter timescales, and PAHs containing less than 60 C atoms will be destroyed. Assuming a classical size distribution for the PAH precursors, our model predicts absolute abundances of C$_{60}$ are up to several $10^{-4}$ of the elemental carbon, i.e. less than a percent of the typical interstellar PAH abundance, which is consistent with observational studies. According to our model, once formed, C$_{60}$ can survive much longer than other fullerenes because of the remarkable stability of the \cs molecule at high internal energies.Hence, a natural consequence is that \cs is more abundant than other fullerenes in highly irradiated environments.
The geometry and physics of the spiral structure of the giant Hubble type Sc
galaxy NGC 309 is studied. A schematic of two patterns with three arms in each
is suggested for the blue spiral. The red and blue patterns form together a
grand design with two-fold symmetry. A possible gas-dynamics explanation of the
phenomenon is suggested which shows how the two-arm red spiral may induce the
formation of the six-arm coherent blue spiral.
Key words: galaxies: individual (NGC 309) -- galaxies: spiral
In this paper, we apply a method identified by Puerari & Dottori (1997) to find the corotation radii (CR) in spiral galaxies. We apply our method to 57 galaxies, 17 of which have already have their CR locations determined using other methods. The method we adopted entails taking Fourier transforms along radial cuts in the u, g, r, i, and z wavebands and comparing the phase angles as a function of radius between them. The radius at which the phase angles cross indicates the location of the corotation radius. We then calculated the relative bar pattern speed, $\mathcal{R}$, and classified the bar as "fast", where $\mathcal{R} < 1.4$, slow, where $\mathcal{R} \geq 1.4$, or intermediate, where the errors on $\mathcal{R}$ are consistent with the bar being "slow" or "fast". For the 17 galaxies that had their CR locations previously measured, we found that our results were consistent with the values of $\mathcal{R}$ obtained by the computer simulations of Rautiainen, Salo & Laurikainen (2008). For the larger sample, our results indicate that 34 out of 57 galaxies (~60%) have fast bars. We discuss these results in the context of its implications for dark matter concentrations in disk galaxies. We also discuss these results in the context of different models for spiral structure in disk galaxies.
The Swift AGN and Cluster Survey (SACS) uses 125 deg^2 of Swift XRT serendipitous fields with variable depths surrounding gamma-ray bursts to provide a medium depth (4e-15 erg/s/cm^2) and area survey filling the gap between deep, narrow Chandra/XMM-Newton surveys and wide, shallow ROSAT surveys. Here we present a catalog of 22,563 point sources and 442 extended sources and examine the number counts of the AGN and galaxy cluster populations. SACS provides excellent constraints on the AGN number counts at the bright end with negligible uncertainties due to cosmic variance, and these constraints are consistent with previous measurements. We use Wise mid-infrared (MIR) colors to classify the sources. For AGN we can roughly separate the point sources into MIR-red and MIR-blue AGN, finding roughly equal numbers of each type in the soft X-ray band (0.5-2 keV), but fewer MIR-blue sources in the hard X-ray band (2-8 keV). The cluster number counts, with 5% uncertainties from cosmic variance, are also consistent with previous surveys but span a much larger continuous flux range. Deep optical or IR follow-up observations of this cluster sample will significantly increase the number of higher redshift (z > 0.5) X-ray-selected clusters.
Links to: arXiv, form interface, find, astro-ph, recent, 1503, contact, help (Access key information)
