We compare the dynamics of satellite galaxies in the Sloan Digital Sky Survey to simple models in order to test the hypothesis that a large fraction of satellites co-rotate in coherent planes. We confirm the previously-reported excess of co-rotating satellite pairs located near diametric opposition with respect to the host, but show that this signal is unlikely to be due to rotating discs (or planes) of satellites. In particular, no overabundance of co-rotating satellites pairs is observed within $\sim 20^{\circ}-50^{\circ}$ of direct opposition, as would be expected for planar distributions inclined relative to the line-of-sight. Instead, the excess co-rotation for satellite pairs within $\sim 10^{\circ}$ of opposition is consistent with random noise associated with undersampling of an underlying isotropic velocity distribution. We conclude that at most $10\%$ of the hosts in our sample harbor co-rotating satellite planes (as traced by the luminous satellite population).
We investigate the relation between star formation (SF) and black hole accretion luminosities, using a sample of 492 type-2 active galactic nuclei (AGNs) at z < 0.22, which are detected in the far-infrared (FIR) surveys with AKARI and Herschel. We adopt FIR luminosities at 90 and 100 um as SF luminosities, assuming the proposed linear proportionality of star formation rate with FIR luminosities. By estimating AGN luminosities from [OIII]5007 and [OI]6300 emission lines, we find a positive linear trend between FIR and AGN luminosities over a wide dynamical range. This result appears to be inconsistent with the recent reports that low-luminosity AGNs show essentially no correlation between FIR and X-ray luminosities, while the discrepancy is likely due to the Malmquist and sample selection biases. By analyzing the spectral energy distribution, we find that pure-AGN candidates, of which FIR radiation is thought to be AGN-dominated, show significantly low-SF activities. These AGNs hosted by low-SF galaxies are rare in our sample (~ 1%). However, the low fraction of low-SF AGN is possibly due to observational limitations since the recent FIR surveys are insufficient to examine the population of high-luminosity AGNs hosted by low-SF galaxies.
To investigate the disk formation and jet launch in protostars is crucial to comprehend the earliest stages of star and planet formation. We aim to constrain the properties of the molecular jet and the disk of the HH 212 protostellar system at unprecedented angular scales through ALMA observations of sulfur-bearing molecules, SO 9(8)-8(7), SO 10(11)-10(10), SO2 8(2,6)-7(1,7). SO 9(8)-8(7) and SO2 8(2,6)-7(1,7) show broad velocity profiles. At systemic velocity they probe the circumstellar gas and the cavity walls. Going from low to high blue-/red-shifted velocities the emission traces the wide-angle outflow and the fast (~100-200 km/s) and collimated (~90 AU) molecular jet revealing the inner knots with timescales <50 years. The jet transports a mass loss rate >0.2-2e-6 Msun/yr, implying high ejection efficiency (>0.03-0.3). The SO and SO2 abundances in the jet are ~1e-7-1e-6. SO 10(11)-10(10) emission is compact and shows small-scale velocity gradients indicating that it originates partly from the rotating disk previously seen in HCO+ and C17O, and partly from the base of the jet. The disk mass is >0.002-0.013 Msun, and the SO abundance in the disk is ~1e-8-1e-7. SO and SO2 are effective tracers of the molecular jet in the inner few hundreds AU from the protostar. Their abundances indicate that 1% - 40% of sulfur is in SO and SO2 due to shocks in the jet/outflow and/or to ambipolar diffusion at the wind base. The SO abundance in the disk is 3-4 orders of magnitude larger than in evolved protoplanetary disks. This may be due to an SO enhancement in the accretion shock at the envelope-disk interface or in spiral shocks if the disk is partly gravitationally unstable.
We have analyzed the archival data of FUV observations for the region of GSH 006-15+7, a large shell-like structure discovered by Moss et al. (2012) from the H I velocity maps. FUV emission is seen to be enhanced in the lower supershell region. The FUV emission is considered to come mainly from the scattering of interstellar photons by dust grains. A corresponding Monte Carlo simulation indicates that the distance to the supershell is 1300 +- 800 pc, which is similar to the previous estimation of 1500 +- 500 pc based on kinematic considerations. The spectrum at lower Galactic latitudes of the supershell exhibits molecular hydrogen fluorescence lines; a simulation model for this candidate photodissociation region (PDR) yields an H_2 column density of N(H_2) = 10^{18.0-20.0} cm^{-2} with a rather high total hydrogen density of n_H ~ 30 cm^{-3}.
By means of our own cosmological-hydrodynamical simulation and
semi-analytical model we studied galaxy population properties in clusters and
groups, spanning over 10 different bands from UV to NIR, and their evolution
since redshift z=2. We compare our results in terms of galaxy red/blue
fractions and luminous-to-faint ratio (LFR) on the Red Sequence (RS) with
recent observational data reaching beyond z=1.5. Different selection criteria
were tested in order to retrieve galaxies belonging to the RS: either by their
quiescence degree measured from their specific SFR ("Dead Sequence"), or by
their position in a colour-colour plane which is also a function of sSFR. In
both cases, the colour cut and the limiting magnitude threshold were let
evolving with redshift, in order to follow the natural shift of the
characteristic luminosity in the LF.
We find that the Butcher-Oemler effect is wavelength-dependent, with the
fraction of blue galaxies increasing steeper in optical colours than in NIR.
Besides, only when applying a lower limit in terms of fixed absolute magnitude,
a steep BO effect can be reproduced, while the blue fraction results less
evolving when selecting samples by stellar mass or an evolving magnitude limit.
We then find that also the RS-LFR behaviour, highly debated in the literature,
is strongly dependent on the galaxy selection function: in particular its very
mild evolution recovered when measured in terms of stellar mass, is in
agreement with values reported for some of the highest redshift confirmed
(proto)clusters. As to differences through environments, we find that normal
groups and (to a lesser extent) cluster outskirts present the highest values of
both star forming fraction and LFR at low z, while fossil groups and cluster
cores the lowest: this separation among groups begins after z~0.5, while
earlier all group star forming properties are undistinguishable.
The dust cloud around $\lambda$ Orionis is observed to be circularly symmetric with a large angular extent ($\approx$ 8 degrees). However, whether the three-dimensional (3D) structure of the cloud is shell- or ring-like has not yet been fully resolved. We study the 3D structure using a new approach that combines a 3D Monte Carlo radiative transfer model for ultraviolet (UV) scattered light and an inverse Abel transform, which gives a detailed 3D radial density profile from a two-dimensional column density map of a spherically symmetric cloud. By comparing the radiative transfer models for a spherical shell cloud and that for a ring cloud, we find that only the shell model can reproduce the radial profile of the scattered UV light, observed using the S2/68 UV observation, suggesting a dust shell structure. However, the inverse Abel transform applied to the column density data from the Pan-STARRS1 dust reddening map results in negative values at a certain radius range of the density profile, indicating the existence of additional, non-spherical clouds near the nebular boundary. The additional cloud component is assumed to be of toroidal ring shape; we subtracted from the column density to obtain a positive, radial density profile using the inverse Abel transform. The resulting density structure, composed of a toroidal ring and a spherical shell, is also found to give a good fit to the UV scattered light profile. We therefore conclude that the cloud around $\lambda$ Ori is composed of both ring and shell structures.
In a new release, Iocco, Pato, and Bertone in arXiv:1505.05181 analyze the consistency of Modified Newtonian Dynamics (MOND) with their compiled Milky Way data and baryonic mass distribution models. We contribute to the discussion by feeding the mean of the seven baryonic mass distribution models that they considered in their original paper into the MOND formula assuming the so-called simple interpolation function, and directly plotting these results on top of the compiled observational rotation curve data from their original paper. Although there is no reason to assume that the mean of the seven baryonic mass distribution models is more correct than any of the individual models, it is a reasonable choice to feed into the equations and one that is less subject to bias inherent in choosing an arbitrary individual model for the MOND analysis to compare to the data. We find that the mean baryonic model using MOND with the simple interpolation function provides a striking fit to the rotation curve observational data with no parameter adjustments required, and we believe that this demonstration is visually transparent and contributes to the valuable discussion on the subject provided by Iocco, Pato, and Bertone. Our results are consistent with the findings of McGaugh (2008), but for an average of many baryonic models instead of just one.
We examine the effect of Galactic diffractive interstellar scintillation as a means of explaining the reported deficit of Fast Radio Burst (FRB) detections at low Galactic latitude. We model the unknown underlying FRB flux density distribution as a power law with a rate scaling as $S_\nu^{-5/2+\delta}$ and account for the fact that the FRBs are detected at unknown positions within the telescope beam. We find that the event rate of FRBs located off the Galactic plane may be enhanced by a factor ~30-300% relative to objects near the Galactic plane without necessarily affecting the slope of the distribution. For FRBs whose flux densities are subject to relatively weak diffractive scintillation, as is typical for events detected at high Galactic latitudes, we demonstrate that an effect associated with Eddington bias is responsible for the enhancement. The magnitude of the enhancement increases with the steepness of the underlying flux density distribution, so that existing limits on the disparity in event rates between high and low Galactic latitudes suggest that the FRB population has a steep differential flux density distribution, scaling as $S_\nu^{-3.5}$ or steeper. Existing estimates of the event rate in the flux density range probed by the High Time Resolution Universe (HTRU) survey overestimate the true rate by a factor of ~3.
GLEAM, the GaLactic and Extragalactic All-sky MWA survey, is a survey of the entire radio sky south of declination +25 deg at frequencies between 72 and 231 MHz, made with the Murchison Widefield Array (MWA) using a drift scan method that makes efficient use of the MWA's very large field-of-view. We present the observation details, imaging strategies and theoretical sensitivity for GLEAM. The survey ran for two years, the first year using 40 kHz frequency resolution and 0.5 s time resolution; the second year using 10 kHz frequency resolution and 2 s time resolution. The resulting image resolution and sensitivity depends on observing frequency, sky pointing and image weighting scheme. At 154 MHz the image resolution is approximately 2.5 x 2.2/cos(DEC+26.7) arcmin with sensitivity to structures up to ~10 deg in angular size. We provide tables to calculate the expected thermal noise for GLEAM mosaics depending on pointing and frequency and discuss limitations to achieving theoretical noise in Stokes I images. We discuss challenges, and their solutions, that arise for GLEAM including ionospheric effects on source positions and linearly polarised emission, and the instrumental polarisation effects inherent to the MWA's primary beam.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
We use the Sternberg et al. (2014) theory for interstellar atomic to molecular (HI-to-H$_2$) conversion to analyze HI-to-H$_2$ transitions in five (low-mass) star-forming and dark regions in the Perseus molecular cloud, B1, B1E, B5, IC348, and NGC1333. The observed HI mass surface densities of 6.3 to 9.2 M$_{\odot}$ pc$^{-2}$ are consistent with HI-to-H$_2$ transitions dominated by HI-dust shielding in predominantly atomic envelopes. For each source, we constrain the dimensionless parameter $\alpha G$, and the effective ratio, $I_{\rm UV}/n$, of the FUV intensity to hydrogen gas density. We find $\alpha G$ values from 5.0 to 47.0, implying characteristic atomic hydrogen densities 11.8 to 1.0 cm$^{-3}$, for $I_{\rm UV} \approx 1$ appropriate for Perseus. Our analysis implies that the dusty HI shielding layers are probably multiphased, with thermally unstable UNM gas in addition to cold CNM within the 21 cm kinematic radius.
We study how X-rays from stellar binary systems and the hot intracluster medium (ICM) affect the radiative cooling rates of gas in galaxies. Our study uses a novel implementation of gas cooling in the moving-mesh hydrodynamics code \textsc{arepo}. X-rays from stellar binaries do not affect cooling at all as their emission spectrum is too hard to effectively couple with galactic gas. In contrast, X-rays from the ICM couple well with gas in the temperature range $10^4 - 10^6$ K. Idealised simulations show that the hot halo radiation field has minimal impact on the dynamics of cooling flows in clusters because of the high virial temperature ($> 10^7$K), making the interaction between the gas and incident photons very ineffective. Satellite galaxies in cluster environments, on the other hand, experience a high radiation flux due to the emission from the host halo. Low mass satellites ($< 10^{12}\rm{M_\odot}$) in particular have virial temperatures that are exactly in the regime where the effect of the radiation field is maximal. Idealised simulations of satellite galaxies including only the effect of host halo radiation (no ram pressure stripping or tidal effects) fields show a drastic reduction in the amount of cool gas formed ($\sim 40\%$) on a short timescale of about $0.5$ Gyrs. A galaxy merger simulation including all the other environmental quenching mechanisms, shows about $20\%$ reduction in the stellar mass of the satellite and about $\sim 30\%$ reduction in star formation rate after $1$ Gyr due to the host hot halo radiation field. These results indicate that the hot halo radiation fields potentially play an important role in quenching galaxies in cluster environments.
During a galaxy merger, the supermassive black hole (SMBH) in each galaxy is thought to sink to the center of the potential and form a supermassive black hole binary; this binary can eject stars via 3-body scattering, bringing the SMBHs ever closer. In a static spherical galaxy model, the binary stalls at a separation of about a parsec after ejecting all the stars in its loss cone -- this is the well-known final parsec problem. However it has been shown that SMBH binaries in non-spherical galactic nuclei harden at a nearly constant rate until reaching the gravitational wave regime. Here we use a suite of direct N-body simulations to follow SMBH binary evolution in both corotating and counterrotating flattened galaxy models. For N larger than 500K, we find that the evolution of the SMBH binary is convergent, and is independent of the particle number. Rotation in general increases the hardening rate of SMBH binaries even more effectively than galaxy geometry alone. SMBH binary hardening rates are similar for co- and counterrotating galaxies. In the corotating case, the center of mass of SMBH binary settles into an orbit that is in a corotation resonance with the background rotating model, and the coalescence time is roughly few hundred Myr faster than a non-rotating flattened model. We find that counterrotation drives SMBHs to coalesce on a nearly radial orbit promptly after forming a hard binary. We discuss the implications for gravitational wave astronomy, hypervelocity star production, and the effect on the structure of the host galaxy.
New observations indicate that ultrafaint dwarf galaxies (UFD) -- the least luminous systems bound by dark matter halos (<10^5 Lsun) -- may have formed before reionization. The extrapolated virial masses today are uncertain with estimates ranging from 10^8 Msun to 10^9 Msun. We show that the progenitor halo masses of UFDs can be as low as Mvir = 10^7 Msun. Under the right conditions, such a halo can survive the energy input of a supernova and its radiative progenitor. A clumpy medium is much less susceptible to both internal and external injections of energy. It is less prone to SN sweeping because the coupling efficiency of the explosive energy is much lower than for a diffuse ISM. With the aid of the 3D hydro/ionization code Fyris, we show that sufficient baryons are retained to form stars following a single supernova event in dark matter halos down to Mvir ~ 10^7 Msun with radiative cooling. The gas survives the SN explosion, is enriched with the abundance yields of the discrete events, and reaches surface densities where low mass stars can form. Our highest resolution simulations reveal why cooling is so effective in retaining gas compared to any other factor. In the early stages, the super-hot metal-enriched SN ejecta exhibit strong cooling, leading to much of the explosive energy being lost. Consistent with earlier work, the baryons do *not* survive in smooth or adiabatic models in the event of a supernova. The smallest galaxies carry signatures of the earliest epochs of star formation, which may distinguish a small primordial galaxy from one that was stripped down to its present size through tidal interaction. We discuss these results in the context of local UFDs and damped Ly-alpha systems (z~2) at very low metallicity ([Fe/H] ~ -3). We show that both classes of objects are consistent with primordial low-mass systems that have experienced only a few enrichment events.
Cold accretion is a primary growth mechanism of simulated galaxies, yet observational evidence of "cold flows" at redshifts where they should be most efficient ($z=2$-4) is scarce. In simulations, cold streams manifest as Lyman-limit absorption systems (LLSs) with low heavy-element abundances similar to those of the diffuse IGM. Here we report on an abundance survey of 17 H I-selected LLSs at $z=3.2$-4.4 which exhibit no metal absorption in SDSS spectra. Using medium-resolution spectra obtained at Magellan, we derive ionization-corrected metallicities (or limits) with a Markov-Chain Monte Carlo sampling that accounts for the large uncertainty in $N_{\rm HI}$ measurements typical of LLSs. The metal-poor LLS sample overlaps with the IGM in metallicity and is best described by a model where $71^{+13}_{-11}\%$ are drawn from the IGM chemical abundance distribution. These represent roughly half of all LLSs at these redshifts, suggesting that 28-40$\%$ of the general LLS population at $z\sim3.7$ could trace unprocessed gas. An ancillary sample of ten LLSs without any a priori metal-line selection is best fit with $48^{+14}_{-12}\%$ of metallicities drawn from the IGM. We compare these results with regions of a moving-mesh simulation; the simulation finds only half as many baryons in IGM-metallicity LLSs, and most of these lie beyond the virial radius of the nearest galaxy halo. A statistically significant fraction of all LLSs have low metallicity and therefore represent candidates for accreting gas; large-volume simulations can establish what fraction of these candidates actually lie near galaxies and the observational prospects for detecting the presumed hosts in emission.
MUSE observations of NGC5813 reveal a complex structure in the velocity dispersion map, previously hinted by SAURON observations. The structure is reminiscent of velocity dispersion maps of galaxies comprising two counter-rotating discs, and may explain the existence of the kinematically distinct core (KDC). Further evidence for two counter-rotating components comes from the analysis of the higher moments of the stellar line-of-sight velocity distributions and fitting MUSE spectra with two separate Gaussian line-of-sight velocity distributions. The emission-line kinematics show evidence of being linked to the present cooling flows and the buoyant cavities seen in X-rays. We detect ionised gas in a nuclear disc-like structure, oriented like the KDC, which is, however, not directly related to the KDC. We build an axisymmetric Schwarzschild dynamical model, which shows that the MUSE kinematics can be reproduced well with two counter-rotating orbit families, characterised by relatively low angular momentum components, but clearly separated in integral phase space and with radially varying contributions. The model indicates that the counter-rotating components in NGC5813 are not thin discs, but dynamically hot structures. Our findings give further evidence that KDCs in massive galaxies should not necessarily be considered as structurally or dynamically decoupled regions, but as the outcomes of the mixing of different orbital families, where the balance in the distribution of mass of the orbital families is crucial. We discuss the formation of the KDC in NGC5813 within the framework of gas accretion, binary mergers and formation of turbulent thick discs from cold streams at high redshift.
We report contributions to cosmic infrared background (CIB) intensities originating from known galaxies, and their companions, at submillimeter wavelengths. Using the publicly-available UltraVISTA catalog, and maps at 250, 350, and 500 {\mu}m from Herschel/SPIRE, we perform a novel measurement that exploits the fact that correlated sources will bias stacked flux densities if the resolution of the image is poor; i.e., we intentionally smooth the image - in effect degrading the angular resolution - before stacking and summing intensities. By smoothing the maps we are capturing the contribution of faint (undetected in K_S ~ 23.4) sources that are physically associated with the detected sources. We find that the cumulative CIB increases with increased smoothing, reaching 9.82 +- 0.78, 5.77 +- 0.43, and 2.32 +- 0.19 nWm^-2/sr at 250, 350, and 500 {\mu}m at 300 arcsec full width half maximum. This corresponds to a fraction of the fiducial CIB of 0.94 +- 0.23, 1.07 +- 0.31, and 0.97 +- 0.26 at 250, 350, and 500 {\mu}m, where the uncertainties are dominated by those of the absolute CIB. We then propose, with a simple model combining parametric descriptions for stacked flux densities and stellar mass functions, that emission from galaxies with log(M/Msun) > 8.5 can account for the entire measured total intensities, and argue against contributions from extended, diffuse emission. Finally, we discuss prospects for future survey instruments to improve the estimates of the absolute CIB levels, and observe any potentially remaining emission at z > 4.
We present the first cosmological galaxy evolved using the modern smoothed particle hydrodynamics (SPH) code GASOLINE2 with superbubble feedback. We show that superbubble-driven galactic outflows powered by Type II supernovae alone can produce $\rm{L^*}$ galaxies with flat rotation curves with circular velocities $\sim 200\; \rm{km/s}$, low bulge-to-disc ratios, and stellar mass fractions that match observed values from high redshift to the present. These features are made possible by the high mass loadings generated by the evaporative growth of superbubbles. Outflows are driven extremely effectively at high redshift, expelling gas at early times and preventing overproduction of stars before $z=2$. Centrally concentrated gas in previous simulations has often lead to unrealistically high bulge to total ratios and strongly peaked rotation curves. We show that supernova-powered superbubbles alone can produce galaxies that agree well with observed properties without the need for additional feedback mechanisms or increased feedback energy. We present additional results arising from properly modelled hot feedback.
Observations of comets and asteroids show that the Solar Nebula that spawned our planetary system was rich in water and organic molecules. Bombardment brought these organics to the young Earth's surface, seeding its early chemistry. Unlike asteroids, comets preserve a nearly pristine record of the Solar Nebula composition. The presence of cyanides in comets, including 0.01% of methyl cyanide (CH3CN) with respect to water, is of special interest because of the importance of C-N bonds for abiotic amino acid synthesis. Comet-like compositions of simple and complex volatiles are found in protostars, and can be readily explained by a combination of gas-phase chemistry to form e.g. HCN and an active ice-phase chemistry on grain surfaces that advances complexity[3]. Simple volatiles, including water and HCN, have been detected previously in Solar Nebula analogues - protoplanetary disks around young stars - indicating that they survive disk formation or are reformed in situ. It has been hitherto unclear whether the same holds for more complex organic molecules outside of the Solar Nebula, since recent observations show a dramatic change in the chemistry at the boundary between nascent envelopes and young disks due to accretion shocks[8]. Here we report the detection of CH3CN (and HCN and HC3N) in the protoplanetary disk around the young star MWC 480. We find abundance ratios of these N-bearing organics in the gas-phase similar to comets, which suggests an even higher relative abundance of complex cyanides in the disk ice. This implies that complex organics accompany simpler volatiles in protoplanetary disks, and that the rich organic chemistry of the Solar Nebula was not unique.
Symbiotic stars are interacting binary systems with the longest orbital periods. They are typically formed by a white dwarf, a red giant and a nebula. These objects are natural astrophysical laboratories for studying the evolution of binaries. Current estimates of the population of Milky Way symbiotic stars vary from 3000 up to 400000. However, the current census is less than 300. The Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST) survey can obtain hundreds of thousands of stellar spectra per year, providing a good opportunity to search for new symbiotic stars. In this work we detect 4 of such binaries among 4,147,802 spectra released by the LAMOST, of which two are new identifications. The first is LAMOST J12280490-014825.7, considered to be an S-type halo symbiotic star. The second is LAMOST J202629.80+423652.0, a D-type symbiotic star.
We have performed a new abundance analysis of Carina Red Giant (RG) stars from spectroscopic data collected with UVES (high resolution) and FLAMES/GIRAFFE (high and medium resolution) at ESO/VLT. The former sample includes 44 RGs, while the latter consists of 65 (high) and ~800 (medium resolution) RGs, covering a significant fraction of the galaxy's RG branch (RGB), and red clump stars. To improve the abundance analysis at the faint magnitude limit, the FLAMES/GIRAFFE data were divided into ten surface gravity and effective temperature bins. The spectra of the stars belonging to the same gravity/temperature bin were stacked. This approach allowed us to increase by at least a factor of five the signal-to-noise ratio in the faint limit (V>20.5mag). We took advantage of the new photometry index cU,B,I introduced by Monelli et al. (2014), as an age and probably a metallicity indicator, to split stars along the RGB. These two stellar populations display distinct [Fe/H] and [Mg/H] distributions: their mean Fe abundances are -2.15$\pm$0.06dex (sig=0.28), and -1.75$\pm$0.03dex (sig=0.21), respectively. The two iron distributions differ at the 75% level. This supports preliminary results by Lemasle et al. (2012) and by Monelli et al. (2014). Moreover, we found that the old and intermediate-age stellar populations have mean [Mg/H] abundances of -1.91$\pm$0.05dex (sig=0.22) and -1.35$\pm$0.03dex (sig=0.22); these differ at the 83% level. Carina's {\alpha}-element abundances agree, within 1sigma, with similar abundances for field Halo stars and for cluster (Galactic, Magellanic) stars. The same outcome applies to nearby dwarf spheroidals and ultra-faint dwarf galaxies, in the iron range covered by Carina stars. Finally, we found evidence of a clear correlation between Na and O abundances, thus suggesting that Carina's chemical enrichment history is quite different than in the globular clusters.
The Mid-IR colors ($F_{8}/F_{24}$) of galaxies together with their IR-UV luminosity correlations can be used to get some insight into the relative abundance of the different dust grain populations present in them. The ELAIS-N1 field contains thousands of galaxies which do not have optical spectra but have been observed in the Mid-IR by {\it Spitzer} and UV by {\it GALEX} making it ideal for these studies. As part of this work we have selected a sample of galaxies from the ELAIS-N1 field which have photometric observations in the MIR and UV as well as photometric redshifts from the SDSS database. We put the constraint that the redshifts are $\le$ 0.1, thereby giving us a total of 309 galaxies. We find that the majority of the galaxies in the sample are PAH dominated due to their high MIR flux ratio. We also find a reasonable correlation between the Mid-IR and the UV luminosities out of which the Mid-IR emission from PAHs at 8 $\mu$m is marginally better correlated than the 24 $\mu$m VSG emission with the UV luminosities. However, if we divide the sample based on their $F_{8}/F_{24}$ ratios which is also an indicator of metallicity, the MIR-UV correlation seems to increase with the $F_{8}/F_{24}$ ratio. But the MIR-UV correlations are not very different for the PAHs and the VSG population within the individual metallicity groups.
(abridged) Methods: We derive maps of submillimeter dust optical depth and effective dust temperature from Herschel data that were calibrated against Planck. After calibration, we then fit a modified blackbody to the long-wavelength Herschel data, using the Planck-derived dust opacity spectral index beta, derived on scales of 30' (or ~1 pc). We use this model to make predictions of the submillimeter flux density at 850 micron, and we compare these in turn with APEX-Laboca observations. Results: A comparison of the submillimeter dust optical depth and near-infrared extinction data reveals evidence for an increased submillimeter dust opacity at high column densities, interpreted as an indication of grain growth in the inner parts of the core. Additionally, a comparison of the Herschel dust model and the Laboca data reveals that the frequency dependence of the submillimeter opacity, described by the spectral index beta, does not change. A single beta that is only slightly different from the Planck-derived value is sufficient to describe the data, beta=1.53+/-0.07. We apply a similar analysis to Barnard 68, a core with significantly lower column densities than FeSt 1-457, and we do not find evidence for grain growth but also a single beta. Conclusions: While we find evidence for grain growth from the dust opacity in FeSt 1-457, we find no evidence for significant variations in the dust opacity spectral index beta on scales 0.02<x<1 pc (or 36"<x<30'). The correction to the Planck-derived dust beta that we find in both cases is on the order of the measurement error, not including any systematic errors, and it would thus be reasonable to directly apply the dust beta from the Planck all-sky dust model. As a corollary, reliable effective temperature maps can be derived which would be otherwise affected by beta variations.
We present the open-source Python code pyMCZ that determines oxygen abundance and its distribution from strong emission lines in the standard metallicity scales, based on the original IDL code of Kewley & Dopita (2002) with updates from Kewley & Ellison (2008), and expanded to include more recently developed scales. The standard strong-line diagnostics have been used to estimate the oxygen abundance in the interstellar medium through various emission line ratios in many areas of astrophysics, including galaxy evolution and supernova host galaxy studies. We introduce a Python implementation of these methods that, through Monte Carlo (MC) sampling, better characterizes the statistical reddening-corrected oxygen abundance confidence region. Given line flux measurements and their uncertainties, our code produces synthetic distributions for the oxygen abundance in up to 13 metallicity scales simultaneously, as well as for E(B-V), and estimates their median values and their 66% confidence regions. In addition, we provide the option of outputting the full MC distributions, and their kernel density estimates. We test our code on emission line measurements from a sample of nearby supernova host galaxies ($z<0.15$) and compare our metallicity results with those from previous methods. We show that our metallicity estimates are consistent with previous methods but yield smaller uncertainties. We also offer visualization tools to assess the spread of the oxygen abundance in the different scales, as well as the shape of the estimated oxygen abundance distribution in each scale, and develop robust metrics for determining the appropriate MC sample size. The code is open access and open source and can be found at https://github.com/nyusngroup/pyMCZ
We introduce Copernicus Complexio (COCO), a high-resolution cosmological N-body simulation of structure formation in the $\Lambda$cdm model. COCO follows an approximately spherical region of radius $\sim 17.4h^{-1}$Mpc, in which the particle mass is $1.1 \times 10^5h^{-1}M_{\odot}$, embedded in a much larger periodic cube followed at lower resolution. Thus, the resolution in the inner volume is 60 times better than in the Millennium-II simulation. COCO gives the dark matter halo mass function over eight orders of magnitude in halo mass; it forms $\sim 60$ halos of galactic size, each resolved with about 10 million particles. The concentration-mass relation of COCO halos deviates from a single power law for masses $M_{200}<$a few$\times 10^{8}h^{-1}M_{\odot}$, where it flattens in agreement with results by Sanchez-Conde et al. We confirm the power-law character of the subhalo mass function, $\overline N(>\mu)\propto\mu^{-s}$, down to a reduced subhalo mass $M_{sub}/M_{200}\equiv\mu=10^{-6}$, with a best-fit power-law index, $s=0.94$, for hosts of mass $\langle M_{200}\rangle=10^12h^{-1}M_{\odot}$, increasing very slowly with host mass. The host-mass invariance of the reduced maximum circular velocity function of subhaloes, $\nu\equiv V_{max}/V_{200}$, hinted at in previous simulations, is clearly demonstrated over five orders of magnitude in host mass. Similarly, we find that the average, normalized radial distribution of subhaloes is approximately universal (i.e. independent of subhalo mass), as previously suggested by the Aquarius simulations of individual halos. Finally, we find that at fixed physical subhalo size, subhaloes in lower mass hosts typically have lower central densities than those in higher mass hosts.
The very high energy Galactic gamma-ray sky is partially opaque in the (0.1-10) PeV energy range. In the light of the recently detected high energy neutrino flux by IceCube, a comparable very high energy gamma-ray flux is expected in any scenario with a sizable Galactic contribution to the neutrino flux. Here we elaborate on the peculiar energy and anisotropy features imposed upon these very high energy gamma-rays by the absorption on the cosmic microwave background photons and Galactic interstellar light. As a notable application of our considerations, we study the prospects of probing the PeV-scale decaying DM scenario, proposed as a possible source of IceCube neutrinos, by extensive air shower (EAS) cosmic ray experiments.
Thanks to their thermal emission, Tidal Disruption Events (TDEs) were detected regularly in the soft X-rays and sometimes in the optical. Only few of them have been detected at hard X-rays: two are high redshift beamed events, one occurred at the core of a nearby galaxy and the last one is of a different nature, involving a compact object in the Milky Way. The aims of this work are to obtain a first sample of hard X-ray selected un-beamed TDEs, to determine their frequency and to probe if TDEs are usually or exceptionally emitting at hard X-rays. We performed extensive search for hard X-ray flares at the positions of over 53000 galaxies up to a distance of 100 Mpc in the Swift BAT archive. Light curves were extracted and parametrized. The quiescent hard X-ray emission was used to exclude persistently active galactic nuclei. Significant flares from non-active galaxies were derived and checked for possible contamination. We found a sample of nine TDE candidates, which translates in a rate of $2 \times 10^{-5}$ galaxy$^{-1}$ yr$^{-1}$ above the BAT detection limit. This rate is consistent with these observed by XMM-Newton at soft X-rays and in the optical from SDSS observations, and expected from simulations. We conclude that hard X-ray emission should be ubiquitous in un-beamed TDEs and that electrons should be accelerated in their accretion flow.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
We present an observational constraint for the typical active galactic nucleus (AGN) phase lifetime. The argument is based on the time lag between an AGN central engine switching on and becoming visible in X-rays, and the time the AGN then requires to photoionize a large fraction of the host galaxy. Based on the typical light travel time across massive galaxies, and the observed fraction of X-ray selected AGN without AGN-photoionized narrow lines, we estimate that the AGN phase typically lasts $\sim10^{5}$ years. This lifetime is short compared to the total growth time of $10^{7}-10^{9}$ years estimated from e.g. the Soltan argument and implies that black holes grow via many such short bursts and that AGN therefore "flicker" on and off. We discuss some consequences of this flickering behavior for AGN feedback and the analogy of X-ray binaries and AGN lifecycles.
We present data and initial results from VLT/X-shooter emission-line spectroscopy of 96 GRB-selected galaxies at 0.1<z<3.6, the largest sample of GRB host spectroscopy available to date. The majority of our GRBs was detected by Swift and 76% are at 0.5<z<2.5 with a median z~1.6. Based on Balmer and/or forbidden lines of oxygen, nitrogen and neon, we measure systemic redshifts, star-formation rates, visual attenuations (A_V), oxygen abundances and emission-line widths (sigma). We find a strong change of the typical physical properties of GRB hosts with redshift. The median SFR, for example, increases from ~0.6 M_sun/yr at z~0.6 up to ~15 M_sun/yr at z~2. A higher ratio of [OIII]/[OII] at higher redshifts leads to an increasing distance of GRB-selected galaxies to the locus of local galaxies in the BPT diagram. Oxygen abundances of the galaxies are distributed between 12+log(O/H)=7.9 and 12+log(O/H)=9.0 with a median of 12+log(O/H)~8.5. The fraction of GRB-selected galaxies with super-solar metallicities is around 20% at z<1 in the adopted metallicity scale. This is significantly less than the fraction of star-formation in similar galaxies, illustrating that GRBs are scarce in high-metallicity environments. At z~3, sensitivity limits us to probing only the most luminous GRB hosts for which we derive metallicities of Z ~< 0.5 Z_sun. Together with a high incidence of galaxies with similar metallicity in our sample at z~1.5, this indicates that the metallicity dependence observed at low redshift will not be dominant at z~3.
(Abridged) We present a spectral analysis of a deep (220 ks) XMM-Newton observation of the Phoenix cluster (SPT-CL J2344-4243), which we also combine with Chandra archival ACIS-I data. We extract CCD and RGS X-ray spectra from the core region to search for the signature of cold gas, and constrain the mass deposition rate in the cooling flow which is thought to be responsible of the massive star formation episode observed in the BCG. We find an average mass deposition rate of $\dot M = 620 (-190 +200)_{stat} (-50 +150)_{syst} M_\odot$/yr in the temperature range 0.3-3.0 keV from MOS data. A temperature-resolved analysis shows that a significant amount of gas is deposited only above 1.8 keV, while upper limits of the order of hundreds of $M_\odot$/yr can be put in the 0.3-1.8 keV temperature range. From pn data we obtain $\dot M = 210 (-80 +85)_{stat} ( -35 +60)_{syst} M_\odot$/yr, and the upper limits from the temperature-resolved analysis are typically a factor of 3 lower than MOS data. In the RGS spectrum, no line emission from ionization states below Fe XXIII is seen above $12 \AA$, and the amount of gas cooling below $\sim 3$ keV has a best-fit value $\dot M = 122_{-122}^{+343}$ $M_{\odot}$/yr. In addition, our analysis of the FIR SED of the BCG based on Herschel data provides $SFR = (530 \pm 50) M_\odot$/yr, significantly lower than previous estimates by a factor 1.5. Current data are able to firmly identify substantial amount of cooling gas only above 1.8 keV in the core of the Phoenix cluster. While MOS data analysis is consistent with values as high as $\dot M \sim 1000$ within $1 \sigma$, pn data provide $\dot M < 500 M_\odot$ yr$^{-1}$ at $3\sigma$ c.l. at temperature below 1.8 keV. At present, this discrepancy cannot be explained on the basis of known calibration uncertainties or other sources of statistical noise.
Context. The LIRG Zw 049.057 contains a compact obscured nucleus where a considerable amount of the galaxy's luminosity is generated. This nucleus contains a dusty environment that is rich in molecular gas. One approach to probing this kind of environment and to revealing what is hidden behind the dust is to study the rotational lines of molecules that couple well with the IR radiation emitted by the dust. Methods. We observed Zw 049.057 with PACS and SPIRE onboard the Herschel Space Observatory in rotational lines of H2O, H218O, OH, 18OH, and [O I]. We modeled the unresolved core of the galaxy using a spherically symmetric radiative transfer code. Results. We present the full SPIRE FTS spectrum of Zw 049.057, along with relevant spectral scans in the PACS range. We find that a minimum of two different components (nuclear and extended) are required in order to account for the rich molecular line spectrum. The nuclear component has a radius of 10-30 pc, a very high infrared surface brightness (1e14 Lsun/kpc2), warm dust (Td > 100 K), and a very large H2 column density (NH2 = 1e24-1e25 cm-2). The modeling also indicates high nuclear H2O (5e-6) and OH (4e-6) abundances relative to H2 as well as a low 16O/18O-ratio of 50-100. We also find a prominent infall signature in the [O I] line. We tentatively detect a 500 km/s outflow in the H2O 313->202 line. Conclusions. The high surface brightness of the core indicates the presence of either a buried active galactic nucleus or a very dense nuclear starburst.nThe H2O abundance is comparable to that of other compact (ultra-)luminous infrared galaxies such as NGC 4418 and Arp 220 - and also to hot cores in the Milky Way. The enhancement of 18O is a possible indicator that the nucleus of Zw 049.057 is in a similar evolutionary stage as the nuclei of Arp 220 - and more advanced than NGC 4418.
The Herschel Space Observatory has had a tremendous impact on the study of extragalactic dust. Specifically, early-type galaxies (ETG) have been the focus of several studies. In this paper we combine results from two Herschel studies - a Virgo cluster study HeViCS and a broader, low-redshift H-ATLAS/GAMA study - and contrast the dust and associated properties for similar mass galaxies. This comparison is motivated by differences in results exhibited between multiple Herschel studies of early-type galaxies. A comparison between consistent modified blackbody derived dust mass is carried out, revealing strong differences between the two samples in both dust mass and dust-to-stellar mass ratio. In particular, the HeViCS sample lacks massive ETG with as high a specific dust content as found in H-ATLAS. This is most likely connected with the difference in environment for the two samples. We calculate nearest neighbour environment densities in a consistent way, showing that H-ATLAS ETG occupy sparser regions of the local Universe, whereas HeViCS ETG occupy dense regions. This is also true for ETG that are not Herschel-detected but are in the Virgo and GAMA parent samples. Spectral energy distributions are fit to the panchromatic data. From these we find that in H-ATLAS the specific star formation rate anticorrelates with stellar mass and reaches values as high as in our Galaxy. On the other hand HeViCS ETG appear to have little star formation. Based on the trends found here, H-ATLAS ETG are thought to have more extended star formation histories and a younger stellar population than HeViCS ETG.
The existence of a classical bulge in disk galaxies holds important clue to
the assembly history of galaxies. Finding observational evidence of very low
mass classical bulges particularly in barred galaxies including our Milky Way,
is a challenging task as the bar driven secular evolution might bring
significant dynamical change to these bulges alongside the stellar disk.
Using high-resolution N-body simulation, we show that if a cool stellar disk
is assembled around a non-rotating low-mass classical bulge, the disk rapidly
grows a strong bar within a few rotation time scales. Later, the bar driven
secular process transform the initial classical bulge into a flattened rotating
stellar system whose central part also have grown a bar-like component rotating
in sync with the disk bar. During this time, a boxy/peanut (hereafter, B/P)
bulge is formed via the buckling instability of the disk bar and the vertical
extent of this B/P bulge being slightly higher than that of the classical
bulge, it encompasses the whole classical bulge. The resulting composite bulge
appears to be both photometrically and kinematically identical to a B/P bulge
without any obvious signature of the classical component. Our analysis suggest
that many barred galaxies in the local universe might be hiding such low-mass
classical bulges. We suggest that stellar population and chemodynamical
analysis might be required in establishing the evidence for such low-mass
classical bulges.
Extensive N-body simulations are among the key means for the study of numerous astrophysical and cosmological phenomena, so various schemes are developed for possibly higher accuracy computations. We demonstrate the principal possibility for revealing the evolution of a perturbed Hamiltonian system with an accuracy independent on time. The method is based on the Laplace transform and the derivation and analytical solution of an evolution equation in the phase space for the resolvent and using computer algebra.
We study the angular clustering of $\sim 6\times 10^5$ NVSS sources on scales $\gtrsim 50 h^{-1}$ Mpc in the context of the $\Lambda$CDM scenario. The analysis partially relies on the redshift distribution of 131 radio galaxies, inferred from the Hercules and CENSORS survey, and an empirical fit to the stellar to halo mass (SHM) relation. For redshifts $z\lesssim 0.7$, the fraction of radio activity versus stellar mass evolves as $f_{_{\rm RL}}\sim M_*^{\alpha_0+\alpha_1 z}$ where $\alpha_0=2.529{\pm0.184}$ and $\alpha_1=1.854^{+0.708}_{-0.761}$. The estimate on $\alpha_0$ is largely driven by the results of Best et al. (2005), while the constraint on $\alpha_1$ is new. We derive a biasing factor $b(z=0.5)=2.093^{+0.164}_{-0.109}$ between radio galaxies and the underlying mass.The function $b(z)=0.33z^2+0.85 z +1.6$ fits well the redshift dependence. We also provide convenient parametric forms for the redshift dependent radio luminosity function, which are consistent with the redshift distribution and the NVSS source count versus flux.
Signs of damping wing absorption attenuating the Lyman-$\alpha$ emission line of the first known $z \sim 7$ quasar, ULAS J1120+0641, recently provided exciting evidence of a significantly neutral IGM. This long-awaited signature of reionization was inferred, in part, from a deficit of flux in the quasar's Lyman-$\alpha$ emission line based on predictions from a composite of lower-redshift quasars. The composite sample was chosen based on its C IV emission line properties; however, as the original study by Mortlock et al. noted, the composite contained a slight velocity offset in C IV compared to ULAS J1120+0641. Here we test whether this offset may be related to the predicted strength of the Lyman-$\alpha$ emission line. We confirm the significant ($\sim 10$ per cent at r.m.s.) scatter in Lyman-$\alpha$ flux for quasars of a given C IV velocity and equivalent width found by Mortlock et al. We further find that among lower-redshift objects chosen to more closely match the C IV properties of ULAS J1120+0641, its Lyman-$\alpha$ emission falls within the observed distribution of fluxes. Among lower-redshift quasars chosen to more closely match in C IV velocity and equivalent width, we find that ULAS J1120+0641 falls within the observed distribution of Lyman-$\alpha$ emission line strengths. This suggests that damping wing absorption may not be present, potentially weakening the case for neutral gas around this object. Larger samples of z$>$7 quasars may therefore be needed to establish a clearer picture of the IGM neutral fraction at these redshifts.
Upcoming HI surveys will deliver large datasets, and automated processing using the full 3-D information (two positional dimensions and one spectral dimension) to find and characterize HI objects is imperative. In this context, visualization is an essential tool for enabling qualitative and quantitative human control on an automated source finding and analysis pipeline. We discuss how Visual Analytics, the combination of automated data processing and human reasoning, creativity and intuition, supported by interactive visualization, enables flexible and fast interaction with the 3-D data, helping the astronomer to deal with the analysis of complex sources. 3-D visualization, coupled to modeling, provides additional capabilities helping the discovery and analysis of subtle structures in the 3-D domain. The requirements for a fully interactive visualization tool are: coupled 1-D/2-D/3-D visualization, quantitative and comparative capabilities, combined with supervised semi-automated analysis. Moreover, the source code must have the following characteristics for enabling collaborative work: open, modular, well documented, and well maintained. We review four state of-the-art, 3-D visualization packages assessing their capabilities and feasibility for use in the case of 3-D astronomical data.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
Estimating the number of microlensing events observed in different parts of the Galactic bulge is a crucial point in planning microlensing experiments. Reliable estimates are especially important if observing resources are scarce, as is the case for space missions: K2, WFIRST, and Euclid. Here we show that the number of detected events can be reliably estimated based on statistics of stars observed in targeted fields. The statistics can be estimated relatively easily, which makes presented method suitable for planning future microlensing experiments.
We develop a simple dynamical model for the evolution of gas in the centres of barred spiral galaxies, using the Milky Way's Central Molecular Zone (CMZ, i.e., the central few hundred pc) as a case study. We show that, in the presence of a galactic bar, gas in a disc in the central regions of a galaxy will be driven inwards by angular momentum transport induced by acoustic instabilities within the bar's inner Lindblad resonance. This transport process drives turbulence within the gas that temporarily keeps it strongly gravitationally stable and prevents the onset of rapid star formation. However, at some point the rotation curve must transition from approximately flat to approximately solid body, and the resulting reduction in shear reduces the transport rates and causes gas to build up, eventually producing a gravitationally-unstable region that is subject to rapid and violent star formation. For the observed rotation curve of the Milky Way, the accumulation happens $\sim 100$ pc from the centre of the Galaxy, in good agreement with the observed location of gas clouds and young star clusters in the CMZ. The characteristic timescale for gas accumulation and star formation is of order $10-20$ Myr. We argue that similar phenomena should be ubiquitous in other barred spiral galaxies.
We present the main steps that will be taken to extract H$\alpha$ emission flux from Javalambre Photometric Local Universe Survey (J-PLUS) photometric data. For galaxies with $z\lesssim0.015$, the H$\alpha$+[NII] emission is covered by the J-PLUS narrow-band filter $F660$. We explore three different methods to extract the H$\alpha$ + [NII] flux from J-PLUS photometric data: a combination of a broad-band and a narrow-band filter ($r'$ and $F660$), two broad-band and a narrow-band one ($r'$, $i'$ and $F660$), and a SED-fitting based method using 8 photometric points. To test these methodologies, we simulated J-PLUS data from a sample of 7511 SDSS spectra with measured H$\alpha$ flux. Based on the same sample, we derive two empirical relations to correct the derived H$\alpha$+[NII] flux from dust extinction and [NII] contamination. We find that the only unbiased method is the SED fitting based one. The combination of two filters underestimates the measurements of the H$\alpha$ + [NII] flux by a 28%, while the three filters method by a 9%. We study the error budget of the SED-fitting based method and find that, in addition to the photometric error, our measurements have a systematic uncertainty of a 4.3%. Several sources contribute to this uncertainty: differences between our measurement procedure and the one used to derive the spectroscopic values, the use of simple stellar populations as templates, and the intrinsic errors of the spectra, which were not taken into account. Apart from that, the empirical corrections for dust extinction and [NII] contamination add an extra uncertainty of 14%. Given the J-PLUS photometric system, the best methodology to extract H$\alpha$ + [NII] flux is the SED-fitting based one. Using this method, we are able to recover reliable H$\alpha$ fluxes for thousands of nearby galaxies in a robust and homogeneous way.
We present a simulation of the long-term evolution of a Population III supernova remnant in a cosmological minihalo. Employing passive Lagrangian tracer particles, we investigate how chemical stratification and anisotropy in the explosion can affect the abundances of the first low-mass, metal-enriched stars. We find that reverse shock heating can leave the inner mass shells at entropies too high to cool, leading to carbon-enhancement in the re-collapsing gas. This hydrodynamic selection effect could explain the observed incidence of carbon-enhanced metal-poor (CEMP) stars at low metallicity. We further explore how anisotropic ejecta distributions, recently seen in direct numerical simulations of core-collapse explosions, may translate to abundances in metal-poor stars. We find that some of the observed scatter in the Population II abundance ratios can be explained by an incomplete mixing of supernova ejecta, even in the case of only one contributing enrichment event. We demonstrate that the customary hypothesis of fully-mixed ejecta clearly fails if post-explosion hydrodynamics prefers the recycling of some nucleosynthetic products over others. Furthermore, to fully exploit the stellar-archaeological program of constraining the Pop III initial mass function from the observed Pop II abundances, considering these hydrodynamical transport effects is crucial. We discuss applications to the rich chemical structure of ultra-faint dwarf satellite galaxies, to be probed in unprecedented detail with upcoming spectroscopic surveys.
Context. The elliptical galaxy NGC 3923 is known to be surrounded by a number
of stellar shells, a probable remnant of an accreted galaxy. Despite its
uniqueness, the deepest images of its outskirts come from the 80s. B\'{i}lek et
al. (2014) predicted a new shell to lie in this region on the basis of the MOND
theory of modified dynamics.
Aims. To obtain the deepest image ever of the galaxy and to map the tidal
features in it.
Methods. The image of the galaxy was taken by the MegaCam camera at the
Canada-France-Hawaii Telescope in the $g'$ band. It reached the
surface-brightness limit of 29 mag arcsec$^{-2}$. Moreover, we reanalyze an
archival HST image of the galaxy.
Results. We detected up to 42 shells in NGC 3923. This is by far most of all
galaxies. We present the description of the shells and other tidal features in
the galaxy. A probable progenitor of some of these features was discovered. The
shell system likely originates from two or more progenitors. The predicted
shell was not detected, but we found that the prediction was based on incorrect
assumptions and poor data.
We recover spiral and barred spiral patterns in disk galaxy simulations with a Wave Dark Matter (WDM) background (also known as Scalar Field Dark Matter (SFDM), Ultra-Light Axion (ULA) dark matter, and Bose-Einstein Condensate (BEC) dark matter). Here we show how the interaction between a baryonic disk and its Dark Matter Halo triggers the formation of spiral structures when the halo is allowed to have a triaxial shape and angular momentum. This is a more realistic picture within the WDM model since a non-spherical rotating halo seems to be more natural. By performing hydrodynamic simulations, along with earlier test particles simulations, we demonstrate another important way in which wave dark matter is consistent with observations. The common existence of bars in these simulations is particularly noteworthy. This may have consequences when trying to obtain information about the dark matter distribution in a galaxy, the mere presence of spiral arms or a bar usually indicates that baryonic matter dominates the central region and therefore observations, like rotation curves, may not tell us what the DM distribution is at the halo center. But here we show that spiral arms and bars can develop in DM dominated galaxies with a central density core without supposing its origin on mechanisms intrinsic to the baryonic matter.
We investigate the dependence of gas kinematics and column densities in the MgII-absorbing circumgalactic medium on galaxy color, azimuthal angle, and inclination to trace baryon cycle processes. Our sample of 30 foreground isolated galaxies at $0.3<z_{\rm gal}<1.0$, imaged with the Hubble Space Telescope, are probed by background quasars within a projected distance of $20<D<110$ kpc. From the high-resolution ($\Delta v\simeq 6.6$ km s$^{-1}$) quasar spectra, we quantify the extent of the absorber velocity structure with pixel-velocity two-point correlation functions. Absorbers with the largest velocity dispersions are associated with blue, face-on ($i<57^{\circ}$) galaxies probed along the projected minor axis ($\Phi \geq 45^{\circ}$), while those with the smallest velocity dispersions belong to red, face-on galaxies along the minor axis. The velocity structure is similar for edge-on ($i \geq 57^{\circ}$) galaxies regardless of galaxy color or azimuthal angle, for red galaxies with azimuthal angle, and for blue and red galaxies probed along the projected major axis ($\Phi<45^{\circ}$). The cloud column densities for face-on galaxies and red galaxies are smaller than for edge-on galaxies and blue galaxies, respectively. These results are consistent with biconical outflows along the minor axis for star-forming galaxies and accreting and/or rotating gas, which is most easily observed in edge-on galaxies probed along the major axis. Gas entrained in outflows may be fragmented with large velocity dispersions, while gas accreting onto or rotating around galaxies may be more coherent due to large path lengths and smaller velocity dispersions. Quiescent galaxies may exhibit little-to-no outflows along the minor axis, while accretion/rotation may exist along the major axis.
We use high-resolution (0.5 arcsec) CO(2-1) observations performed with ALMA to trace the kinematics of the molecular gas in the Seyfert 2 galaxy IC5063. A fast outflow of molecular gas extends along the entire radio jet, with the highest outflow velocities about 0.5kpc from the nucleus, at the location of the brighter hot-spot in the W lobe. The data show that a massive, fast outflow with velocities up to 650 km/s of cold molecular gas is present, in addition to one detected earlier in warm H2, HI and ionised gas. Both the central AGN and the radio jet could energetically drive the outflow. However, the characteristics of the outflowing gas point to the radio jet being the main driver. This is important, because IC5063, although one of the most powerful Seyfert galaxies, is a relatively weak radio source (P = 3x10^23 W/Hz). All the observed characteristics can be described by a scenario of a radio plasma jet expanding into a clumpy medium, interacting directly with the clouds and inflating a cocoon that drives a lateral outflow into the interstellar medium. This model is consistent with results obtained by recent simulations such as those of Wagner et al.. A stronger, direct interaction between the jet and a gas cloud is present at the location of the brighter W lobe. Even assuming the most conservative values for the conversion factor CO-to-H2, the mass of the outflowing gas is between 1.9 and 4.8x10^7 Msun. These amounts are much larger than those of the outflow of warm gas (molecular and ionized) and somewhat larger than of the HI outflow. This suggests that most of the observed cold molecular outflow is due to fast cooling after being shocked. This gas is the end product of the cooling process. Our CO observations demonstrate that fast outflows of molecular gas can be driven by relativistic jets.
In a Universe where AGN feedback regulates star formation in massive galaxies, a strong correlation between these two quantities is expected. If the gas causing star formation is also responsible for feeding the central black hole, then a positive correlation is expected. If powerful AGNs are responsible for the star formation quenching, then a negative correlation is expected. Observations so far have mainly found a mild correlation or no correlation at all (i.e. a flat relation between star formation rate (SFR) and AGN luminosity), raising questions about the whole paradigm of "AGN feedback". In this paper, we report the predictions of the GALFORM semi-analytical model, which has a very strong coupling between AGN activity and quenching of star formation. The predicted SFR-AGN luminosity correlation appears negative in the low AGN luminosity regime, where AGN feedback acts, but becomes strongly positive in the regime of the brightest AGN. Our predictions reproduce reasonably well recent observations by Rosario et al., yet there is some discrepancy in the normalisation of the correlation at low luminosities and high redshifts. Though this regime could be strongly influenced by observational biases, we argue that the disagreement could be ascribed to the fact that GALFORM neglects AGN variability effects. Interestingly, the galaxies that dominate the regime where the observations imply a weak correlation are massive early-type galaxies that are subject to AGN feedback. Nevertheless, these galaxies retain high enough molecular hydrogen contents to maintain relatively high star formation rates and strong infrared emission.
We combine Herschel observations of a total of 12 sources to construct the most uniform survey of HF and H2O in our Galactic disk. Both molecules are detected in absorption along all sight lines. The high spectral resolution of the Heterodyne Instrument for the Far-Infrared (HIFI) allows us to compare the HF and H2O distributions in 47 diffuse cloud components sampling the disk. We find that the HF and H2O velocity distributions follow each other almost perfectly and establish that HF and H2O probe the same gas-phase volume. Our observations corroborate theoretical predictions that HF is a sensitive tracer of H2 in diffuse clouds, down to molecular fractions of only a few percent. Using HF to trace H2 in our sample, we find that the N(H2O)-to-N(HF) ratio shows a narrow distribution with a median value of 1.51. Our results further suggest that H2O might be used as a tracer of H2 -within a factor 2.5- in the diffuse interstellar medium. We show that the measured factor of ~2.5 variation around the median is driven by true local variations in the H2O abundance relative to H2 throughout the disk. The latter variability allows us to test our theoretical understanding of the chemistry of oxygen-bearing molecules in the diffuse gas. We show that both gas-phase and grain-surface chemistry are required to reproduce our H2O observations. This survey thus confirms that grain surface reactions can play a significant role in the chemistry occurring in the diffuse interstellar medium n_H < 1000 cm^-3.
We have performed multivariate statistical analyses of photometric and chemical abundance parameters of three large samples of stars in the globular cluster $\omega$ Centauri. The statistical analysis of a sample of 735 stars based on seven chemical abundances with the method of Maximum Parsimony (cladistics) yields the most promising results: seven groups are found, distributed along three branches with distinct chemical, spatial and kinematical properties. A progressive chemical evolution can be traced from one group to the next, but also within groups, suggestive of an inhomogeneous chemical enrichment of the initial interstellar matter. The adjustment of stellar evolution models shows that the groups with metallicities [Fe/H]\textgreater{}-1.5 are Helium-enriched, thus presumably of second generation. The spatial concentration of the groups increases with chemical evolution, except for two groups, which stand out in their other properties as well. The amplitude of rotation decreases with chemical evolution, except for two of the three metal-rich groups, which rotate fastest, as predicted by recent hydrodynamical simulations. The properties of the groups are interpreted in terms of star formation in gas clouds of different origins. In conclusion, our multivariate analysis has shown that metallicity alone cannot segregate the different populations of $\omega$ Centauri.
We measure the merger fraction of Type 2 radio-loud and radio-quiet active galactic nuclei at z>1 using new samples. The objects have HST images taken with WFC3 in the IR channel. These samples are compared to the 3CR sample of radio galaxies at z>1 and to a sample of non-active galaxies. We also consider lower redshift radio galaxies with HST observations and previous generation instruments (NICMOS and WFPC2). The full sample spans an unprecedented range in both redshift and AGN luminosity. We perform statistical tests to determine whether the different samples are differently associated with mergers. We find that all (92%) radio-loud galaxies at z>1 are associated with recent or ongoing merger events. Among the radio-loud population there is no evidence for any dependence of the merger fraction on either redshift or AGN power. For the matched radio-quiet samples, only 38% are merging systems. The merger fraction for the sample of non-active galaxies at z>1 is indistinguishable from radio-quiet objects. This is strong evidence that mergers are the triggering mechanism for the radio-loud AGN phenomenon and the launching of relativistic jets from supermassive black holes. We speculate that major BH-BH mergers play a major role in spinning up the central supermassive black holes in these objects.
Neutral hydrogen represents the major observable baryonic constituent of galaxies that fuels the formation of stars through the transformation in molecular hydrogen. The emission of the hydrogen recombination line Halpha is the most direct tracer of the process that transforms gas (fuel) into stars. We continue to present Halpha3 (acronym for Halpha-alpha-alpha), an extensive Halpha+[NII] narrow-band imaging campaign of galaxies selected from the HI Arecibo Legacy Fast ALFA Survey (ALFALFA), using the instrumentation available at the San Pedro Martir observatory (Mexico). In only four years since 2011 we were able to complete in 48 nights the Halpha imaging observations of 724 galaxies in the region of the Coma supercluster 10^h < R.A. <16^h; 24^o < Dec. <28^o and 3900<cz<9000 kms^{-1}. Of these, 603 are selected from the HI Arecibo Legacy Fast ALFA Survey (ALFALFA) and constitute a 97% complete sample. They provide for the first time a complete census of the massive star formation properties of local gas-rich galaxies belonging to different environments (cluster vs filaments), morphological type (spirals vs dwarf Irr), over a wide range of stellar mass (10^{8}-10^{11.5} Modot) in the Coma Supercluster. The present Paper V provides the Halpha data and the derived star formation rates for the observed galaxies.
We conduct a detailed lensing, dynamics and stellar population analysis of nine massive lens early-type galaxies (ETGs) from the X-Shooter Lens Survey (XLENS). Combining gravitational lensing constraints from HST imaging with spatially-resolved kinematics and line-indices constraints from VLT X-Shooter (XSH) spectra, we infer the low-mass slope and the low cut-off mass of the stellar Initial Mass Function (IMF): $x_{250}=2.37^{+0.12}_{-0.12}$ and $M_{{\rm low}, 250}= 0.131^{+0.023}_{-0.026}\, M_{\odot}$, respectively, for a reference point with $\sigma \equiv 250\, {{\rm kms}}^{-1}$ and R$_{{\rm eff}} \equiv 10$ kpc. All the XLENS systems are consistent with an IMF slope steeper than Milky Way-like. We find no significant correlations between IMF slope and any other quantity, except for an anti-correlation between total dynamical mass density and low-mass IMF slope at the 87% CL [$dx/d\log(\rho)$ = $ -0.19^{+0.15}_{-0.15}$]. This anti-correlation is consistent with the low redshift lenses found by Smith et al. (2015) that have high velocity dispersions and high stellar mass densities but surprisingly shallow IMF slopes.
We have derived the Galactic bulge initial mass function of the SWEEPS field down to 0.15 $M_{\odot}$, using deep photometry collected with the Advanced Camera for Surveys on the Hubble Space Telescope. Observations at several epochs, spread over 9 years, allowed us to separate the disk and bulge stars down to very faint magnitudes, $F814W \approx$ 26 mag, with a proper-motion accuracy better than 0.5 mas/yr (20 km/s). This allowed us to determine the initial mass function of the pure bulge component uncontaminated by disk stars for this low-reddening field in the Sagittarius window. In deriving the mass function, we took into account the presence of unresolved binaries, errors in photometry, distance modulus and reddening, as well as the metallicity dispersion and the uncertainties caused by adopting different theoretical color-temperature relations. We found that the Galactic bulge initial mass function can be fitted with two power laws with a break at $M \sim$ 0.56 $M_{\odot}$, the slope being steeper ($\alpha = -2.41\pm$0.50) for the higher masses, and shallower ($\alpha = -1.25\pm$0.20) for the lower masses. In the high-mass range, our derived mass function agrees well with the mass function derived for other regions of the bulge. In the low-mass range however, our mass function is slightly shallower, which suggests that separating the disk and bulge components is particularly important in the low-mass range. The slope of the bulge mass function is also similar to the slope of the mass function derived for the disk in the high-mass regime, but the bulge mass function is slightly steeper in the low-mass regime. We used our new mass function to derive stellar mass--to--light values for the Galactic bulge and we found $M/L_{F814W} =$ 2.2$\pm$0.3 and $M/L_{F606W} =$ 3.2$\pm$0.5.
We present an analysis of the diffuse soft X-ray emission from the nuclear region of M51 combining both XMM-Newton RGS and Chandra data. Most of the RGS spectrum of M51 can be fitted with a thermal model with a temperature of $\sim0.5$ keV except for the OVII triplet, which is forbidden-line dominated. The Fe L-shell lines peak around the southern cloud, where the OVIII and NVII Lya lines also peak. In contrast, the peak of the OVII forbidden line is about 10$"$ offset from that of the other lines, indicating that it is from a spatially distinct component. The spatial distribution of the OVII triplet mapped by the Chandra data shows that most of the OVII triplet flux is located at faint regions near edges, instead of the southern cloud where other lines peak. This distribution of the OVII triplet is inconsistent with the photoionization model. Other mechanisms that could produce the anomalous OVII triplet, including a recombining plasma and charge exchange X-ray emission, are discussed.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)
Evidence for non-zero mean stellar velocities in the direction perpendicular to the Galactic plane has been accumulating from various recent large spectroscopic surveys. Previous analytical and numerical work has shown that a "breathing mode" of the Galactic disc, similar to what is observed in the Solar vicinity, can be the natural consequence of a non-axisymmetric internal perturbation of the disc. Here we provide a general analytical framework, in the context of perturbation theory, allowing us to compute the vertical bulk motions generated by a single internal perturber (bar or spiral pattern). In the case of the Galactic bar, we show that these analytically predicted bulk motions are well in line with the outcome of a numerical simulation. The mean vertical motions induced by the Milky Way bar are small (mean velocity of less than 1 km/sec) and cannot be responsible alone for the observed breathing mode, but they are existing. Our analytical treatment is valid close to the plane for all the non-axisymmetric perturbations of the disc that can be described by small-amplitude Fourier modes. Further work should study how the coupling of multiple internal perturbers and external perturbers is affecting the present analytical results.
The evidence that stellar systems surrounding the Milky Way (MW) are distributed in a Vast Polar Structure (VPOS) may be observationally biased by satellites detected in surveys of the northern sky. The recent discoveries of more than a dozen new systems in the southern hemisphere thus constitute a critical test of the VPOS phenomenon. We report that the new objects are located close to the original VPOS, with half of the sample having offsets less than 20 kpc. The positions of the new satellite galaxy candidates are so well aligned that the orientation of the revised best-fitting VPOS structure is preserved to within 9 degrees and the VPOS flattening is almost unchanged (31 kpc height). Interestingly, the shortest distance of the VPOS plane from the MW center is now only 2.5 kpc, indicating that the new discoveries balance out the VPOS at the Galactic center. The vast majority of the MW satellites are thus consistent with sharing a similar orbital plane as the Magellanic Clouds, confirming a hypothesis proposed by Kunkel & Demers and Lynden-Bell almost 40 years ago. We predict the absolute proper motions of the new objects assuming they orbit within the VPOS. Independent of the VPOS results we also predict the velocity dispersions of the new systems under three distinct assumptions: that they (i) are dark-matter-free star clusters obeying Newtonian dynamics, (ii) are dwarf satellites lying on empirical scaling relations of galaxies in dark matter halos, and (iii) obey MOND.
We present results from a sample of XMM-Newton and Suzaku observations of interstellar clouds that cast shadows in the soft X-ray background (SXRB) - the first uniform analysis of such a sample from these missions. By fitting to the on- and off-shadow spectra, we separated the foreground and Galactic halo components of the SXRB. We tested different foreground models - two solar wind charge exchange (SWCX) models and a Local Bubble (LB) model. We also examined different abundance tables. We found that Anders & Grevesse (1989) abundances, commonly used in previous SXRB studies, may result in overestimated foreground brightnesses and halo temperatures. We also found that assuming a single solar wind ionization temperature for a SWCX model can lead to unreliable results. We compared our measurements of the foreground emission with predictions of the SWCX emission from a smooth solar wind, finding only partial agreement. Using available observation-specific SWCX predictions and various plausible assumptions, we placed an upper limit on the LB's OVII intensity of ~0.8 photons/cm^2/s/sr (90% confidence). Comparing the halo results obtained with SWCX and LB foreground models implies that, if the foreground is dominated by SWCX and is brighter than ~1.5e-12 erg/cm^2/s/deg^2 (0.4-1.0 keV), then using an LB foreground model may bias the halo temperature upward and the 0.5-2.0 keV surface brightness downward by ~(0.2-0.3)e6 K and ~(1-2)e-12 erg/cm^2/s/deg^2, respectively. Similarly, comparing results from different observatories implies that there may be uncertainties in the halo temperature and surface brightness of up to ~0.2e6 K and ~25%, respectively, in addition to the statistical uncertainties. These uncertainties or biases may limit the ability of X-ray measurements to discriminate between Galactic halo models.
We present a study of the M83 cluster population, covering the disc of the galaxy between radii of 0.45 and 4.5 kpc. We aim to probe the properties of the cluster population as a function of distance from the galactic centre. We observe a net decline in cluster formation efficiency ($\Gamma$, i.e. the amount of star formation happening in bound clusters) from about 19 % in the inner region to 7 % in the outer part of the galaxy. The recovered $\Gamma$ values within different regions of M83 follow the same $\Gamma$ versus star formation rate density relation observed for entire galaxies. We also probe the initial cluster mass function (ICMF) as a function of galactocentric distance. We observe a significant steepening of the ICMF in the outer regions (from $-1.90\pm0.11$ to $-2.70\pm0.14$) and for the whole galactic cluster population (slope of $-2.18\pm0.07$) of M83. We show that this change of slope reflects a more fundamental change of the 'truncation mass' at the high-mass end of the distribution. This can be modelled as a Schechter function of slope $-2$ with an exponential cut-off mass ($M_{\rm c}$) that decreases significantly from the inner to the outer regions (from 4.00 to $0.25\times 10^5$ M$_\odot$) while the galactic $M_{\rm c}$ is $\approx1.60\times10^5$ M$_\odot$. The trends in \Gamma and ICMF are consistent with the observed radial decrease of the $\Sigma({\rm H}_2)$, hence in gas pressure. As gas pressure declines cluster formation becomes less efficient. We conclude that the host galaxy environment appears to regulate 1) the fraction of stars locked in clusters; 2) the upper mass limit of the ICMF, consistently described by a near-universal slope $-2$ truncated at the high-mass end.
The Lya emission line of HI is intrinsically the brightest feature in the spectrum of astrophysical nebulae, making it a very attractive observational tool with which to survey galaxies. Moreover as a UV resonance line, Lya possesses several unique characteristics that make it useful to study the ISM and ionizing stellar population at all cosmic epochs. In this review I present a summary of Lya observations of galaxies in the nearby universe. At UV magnitudes reachable with current facilities, only ~5% of the local galaxy population shows a Lya equivalent width (EW_Lya) that exceeds 20\AA. This fraction increases dramatically at higher z, but only in the local universe can we study galaxies in detail and assemble unprecedented multi-wavelength datasets. I discuss many local Lya observations, showing that when galaxies show net Lya emission, they ubiquitously produce large halos of scattered Lya, that dominate the integrated luminosity. We discuss how global EW_Lya and the Lya escape fraction (fescLya) are higher (EW_Lya >~ 20\AA\ and fescLya> 10%) in galaxies that represent the less massive and younger end of the distributions for local objects. This is connected with various properties, such that Lya-emitters have lower metallicities (median value of 12+log(O/H) ~ 8.1) and dust reddening. However, the presence of galactic outflows is also vital to Doppler shift the Lya line out of resonance with the HI, as high EW_Lya is found only among galaxies with winds faster than ~50km/s. The evidence is then assembled into a coherent picture, and the requirement for star formation driven feedback is discussed in the context of an evolutionary sequence where the ISM is accelerated and/or subject to fluid instabilities, which reduce the scattering of Lya. Concluding remarks take the form of perspectives upon the most pressing questions that can be answered by observation.
An algorithm is described for evolving the phase-space density of stars or compact objects around a massive black hole at the center of a galaxy. The technique is based on numerical integration of the Fokker-Planck equation in energy-angular momentum space, f(E,L,t), and includes, for the first time, diffusion coefficients that describe the effects of both random and correlated encounters (resonant relaxation), as well as energy loss due to emission of gravitational waves. Destruction or loss of stars into the black hole are treated by means of a detailed boundary-layer analysis. Performance of the algorithm is illustrated by calculating two-dimensional, time-dependent and steady-state distribution functions and their corresponding loss rates.
We present Australia Telescope Compact Array radio data of the dwarf irregular galaxy ESO 324-G024 which is seen in projection against the giant, northern lobe of the radio galaxy Centaurus A (Cen A, NGC 5128). The distorted morphology and kinematics of ESO 324-G024, as observed in the 21 cm spectral line emission of neutral hydrogen, indicate disruptions by external forces. We investigate whether tidal interactions and/or ram pressure stripping are responsible for the formation of the HI tail stretching to the northeast of ESO 324-G024 with the latter being most probable. Furthermore, we closely analyze the sub-structure of Cen A's polarized radio lobes to ascertain whether ESO 324-G024 is located in front, within or behind the northern lobe. Our multi-wavelength, multi-component approach allows us to determine that ESO 324-G024 is most likely behind the northern radio lobe of Cen A. This result helps to constrain the orientation of the lobe, which is likely inclined to our line of sight by approximately 60 degrees if NGC 5128 and ESO 324-G024 are at the same distance.
We present MUFFIT, a new generic code optimized to retrieve the main stellar population parameters of galaxies in photometric multi-filter surveys, and we check its reliability and feasibility with real galaxy data from the ALHAMBRA survey. Making use of an error-weighted $\chi^2$-test, we compare the multi-filter fluxes of galaxies with the synthetic photometry of mixtures of two single stellar populations at different redshifts and extinctions, to provide through a Monte Carlo method the most likely range of stellar population parameters (mainly ages and metallicities), extinctions, redshifts, and stellar masses. To improve the diagnostic reliability, MUFFIT identifies and removes from the analysis those bands that are significantly affected by emission lines. We highlight that the retrieved age-metallicity locus for a sample of $z \le 0.22$ early-type galaxies in ALHAMBRA at different stellar mass bins are in very good agreement with the ones from SDSS spectroscopic diagnostics. Moreover, a one-to-one comparison between the redshifts, ages, metallicities, and stellar masses derived spectroscopically for SDSS and by MUFFIT for ALHAMBRA reveals good qualitative agreements in all the parameters. In addition, and using as input the results from photometric-redshift codes, MUFFIT improves the photometric-redshift accuracy by $\sim 10$-$20\%$, and it also detects nebular emissions in galaxies, providing physical information about their strengths. Our results show the potential of multi-filter galaxy data to conduct reliable stellar population studies with the appropiate analysis techniques, as MUFFIT.
We carry out three-dimensional hydrodynamic simulations of the supernova remnants (SNRs) produced inside molecular clouds (MCs) near their surface using the HLL code (Harten et al. 1983). We explore the dynamical evolution and the X-ray morphology of SNRs after breaking through the MC surface for ranges of the explosion depths below the surface and the density ratios of the clouds to the intercloud media (ICM). We find that if an SNR breaks out through an MC surface in its Sedov stage, the outermost dense shell of the remnant is divided into several layers. The divided layers are subject to the Rayleigh-Taylor instability and fragmented. On the other hand, if an SNR breaks through an MC after the remnant enters the snowplow phase, the radiative shell is not divided to layers. We also compare the predictions of previous analytic solutions for the expansion of SNRs in stratified media with our onedimensional simulations. Moreover, we produce synthetic X-ray surface brightness in order to research the center-bright X-ray morphology shown in thermal composite SNRs. In the late stages, a breakout SNR shows the center-bright X-ray morphology inside an MC in our results. We apply our model to the observational results of the X-ray morphology of the thermal composite SNR 3C 391.
We investigate effects of a global magnetic field on the dynamics of an ensemble of clumps within a magnetized advection-dominated accretion flow by neglecting interactions between the clumps and then solving the collisionless Boltzman equation. In the strong-coupling limit, in which the averaged radial and the rotational velocities of the clumps follow the ADAF dynamics, the averaged radial velocity square of the clumps is calculated analytically for different magnetic field configurations. The value of the averaged radial velocity square of the clumps increases with increasing the strength of the radial or vertical components of the magnetic field. But a purely toroidal magnetic field geometry leads to a reduction of the value of the averaged radial velocity square of the clumps at the inner parts with increasing the strength of this component. Moreover, dynamics of the clumps strongly depends on the amount of the advected energy so that the value of the averaged radial velocity square of the clumps increases in the presence of a global magnetic field as the flow becomes more advective.
We observed with the JVLA at 3.6 and 1.3 cm a sample of 11 proto-brown dwarf candidates in Taurus in a search for thermal radio jets driven by the most embedded brown dwarfs. We detected for the first time four thermal radio jets in proto-brown dwarf candidates. We compiled data from UKIDSS, 2MASS, Spitzer, WISE and Herschel to build the Spectral Energy Distribution (SED) of the objects in our sample, which are similar to typical Class~I SEDs of Young Stellar Objects (YSOs). The four proto-brown dwarf candidates driving thermal radio jets also roughly follow the well-known trend of centimeter luminosity against bolometric luminosity determined for YSOs, assuming they belong to Taurus, although they present some excess of radio emission compared to the known relation for YSOs. Nonetheless, we are able to reproduce the flux densities of the radio jets modeling the centimeter emission of the thermal radio jets using the same type of models applied to YSOs, but with corresponding smaller stellar wind velocities and mass-loss rates, and exploring different possible geometries of the wind or outflow from the star. Moreover, we also find that the modeled mass outflow rates for the bolometric luminosities of our objects agree reasonably well with the trends found between the mass outflow rates and bolometric luminosities of YSOs, which indicates that, despite the "excess" centimeter emission, the intrinsic properties of proto-brown dwarfs are consistent with a continuation of those of very low mass stars to a lower mass range. Overall, our study favors the formation of brown dwarfs as a scaled-down version of low-mass stars.
Selection effects can bedevil the inference of the properties of a population of astronomical catalogues, unavoidably biasing the observed catalogue. This is particularly true when mapping interstellar extinction in three dimensions: more extinguished stars are fainter and so generally less likely to appear in any magnitude limited catalogue of observations. This paper demonstrates how to account for this selection effect when mapping extinction, so that accurate and unbiased estimates of the true extinction are obtained. We advocate couching the description of the problem explicitly as a Poisson point process, which allows the likelihoods employed to be easily and correctly normalised in such a way that accounts for the selection functions applied to construct the catalogue of observations.
We investigate the variations of the radial O/H abundance gradients from planetary nebulae (PN) located at different distances from the galactic plane. In particular, we determine the abundance gradients at different heights from the plane in order to investigate a possible gradient inversion for the objects at larger distances from the plane. We consider a large sample of PN with known distances, so that the height relative to the galactic plane can be derived, and accurate abundances, so that the gradients can be determined.
We present Karl G. Jansky Very Large Array (VLA) observations of the CO ($J = 2 \rightarrow 1$) line emission towards the $z = 6.419$ quasar SDSS J$114816.64+525150.3$ (J$1148+5251$). The molecular gas is found to be marginally resolved with a major axis of $0.9"$ (consistent with previous size measurements of the CO ($J = 7 \rightarrow 6$) emission). We observe tentative evidence for extended line emission towards the south west on a scale of ~$1.4"$, but this is only detected at $3.3\sigma$ significance and should be confirmed. The position of the molecular emission region is in excellent agreement with previous detections of low frequency radio continuum emission as well as [C ii] line and thermal dust continuum emission. These CO ($J = 2 \rightarrow 1$) observations provide an anchor for the low excitation part of the molecular line SED. We find no evidence for extended low excitation component, neither in the spectral line energy distribution nor the image. We fit a single kinetic gas temperature model of 50 K. We revisit the gas and dynamical masses in light of this new detection of a low order transition of CO, and confirm previous findings that there is no extended reservoir of cold molecular gas in J$1148+5251$, and that the source departs substantially from the low $z$ relationship between black hole mass and bulge mass. Hence, the characteristics of J$1148+5251$ at $z = 6.419$ are very similar to $z$~$2$ quasars, in the lack of a diffuse cold gas reservoir and kpc-size compactness of the star forming region.
The Monte Carlo method is the most popular technique to perform radiative transfer simulations in a general 3D geometry. The algorithms behind and acceleration techniques for Monte Carlo radiative transfer are discussed extensively in the literature, and many different Monte Carlo codes are publicly available. On the contrary, the design of a suite of components that can be used for the distribution of sources and sinks in radiative transfer codes has received very little attention. The availability of such models, with different degrees of complexity, has many benefits. For example, they can serve as toy models to test new physical ingredients, or as parameterised models for inverse radiative transfer fitting. For 3D Monte Carlo codes, this requires algorithms to efficiently generate random positions from 3D density distributions. We describe the design of a flexible suite of components for the Monte Carlo radiative transfer code SKIRT. The design is based on a combination of basic building blocks (which can be either analytical toy models or numerical models defined on grids or a set of particles) and the extensive use of decorators that combine and alter these building blocks to more complex structures. For a number of decorators, e.g. those that add spiral structure or clumpiness, we provide a detailed description of the algorithms that can be used to generate random positions. Advantages of this decorator-based design include code transparency, the avoidance of code duplication, and an increase in code maintainability. Moreover, since decorators can be chained without problems, very complex models can easily be constructed out of simple building blocks. Finally, based on a number of test simulations, we demonstrate that our design using customised random position generators is superior to a simpler design based on a generic black-box random position generator.
Links to: arXiv, form interface, find, astro-ph, recent, 1505, contact, help (Access key information)