MID-IR luminosities and UV/optical star formation rates at z < 1.4

Ultraviolet (UV) nonionizing continuum and mid-infrared (IR) emission constitute the basis of two widely used star formation (SF) indicators at intermediate and high redshifts. We study 2430 galaxies with z < 1.4 in the Extended Groth Strip with deep MIPS 24 μm observations from FIDEL, spectroscopy from DEEP2, and UV, optical, and near-IR photometry from the AEGIS. The data are coupled with dust-reddened stellar population models and Bayesian spectral energy distribution (SED) fitting to estimate dust-corrected star formation rates (SFRs). In order to probe the dust heating from stellar populations of various ages, the derived SFRs were averaged over various timescales—from 100 Myr for "current" SFR (corresponding to young stars) to 1-3 Gyr for long-timescale SFRs (corresponding to the light-weighted age of the dominant stellar populations). These SED-based UV/optical SFRs are compared to total IR luminosities extrapolated from 24 μm observations, corresponding to 10-18 μm rest frame. The total IR luminosities are in the range of normal star-forming galaxies and luminous IR galaxies (10^10-10^12 L_☉). We show that the IR luminosity can be estimated from the UV and optical photometry to within a factor of 2, implying that most z < 1.4 galaxies are not optically thick. We find that for the blue, actively star-forming galaxies the correlation between the IR luminosity and the UV/optical SFR shows a decrease in scatter when going from shorter to longer SFR-averaging timescales. We interpret this as the greater role of intermediate age stellar populations in heating the dust than what is typically assumed. Equivalently, we observe that the IR luminosity is better correlated with dust-corrected optical luminosity than with dust-corrected UV light. We find that this holds over the entire redshift range. Many so-called green valley galaxies are simply dust-obscured actively star-forming galaxies. However, there exist 24 μm detected galaxies, some with L_IR>10^11 L_☉, yet with little current SF. For them a reasonable amount of dust absorption of stellar light (but presumably higher than in nearby early-type galaxies) is sufficient to produce the observed levels of IR, which includes a large contribution from intermediate and old stellar populations. In our sample, which contains very few ultraluminous IR galaxies, optical and X-ray active galactic nuclei do not contribute on average more than ~50% to the mid-IR luminosity, and we see no evidence for a large population of "IR excess" galaxies.

Testing quasar unification: radiative transfer in clumpy winds

Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to 1043 erg s-1. Here, we introduce a simple treatment of clumping, and find that a filling factor of ˜0.01 moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Lyα and C IV 1550 Å at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question.

The microphysics of collisionless shock waves

Collisionless shocks, that is shocks mediated by electromagnetic processes, are customary in space physics and in astrophysics. They are to be found in a great variety of objects and environments: magnetospheric and heliospheric shocks, supernova remnants, pulsar winds and their nebulæ, active galactic nuclei, gamma-ray bursts and clusters of galaxies shock waves. Collisionless shock microphysics enters at different stages of shock formation, shock dynamics and particle energization and/or acceleration. It turns out that the shock phenomenon is a multi-scale non-linear problem in time and space. It is complexified by the impact due to high-energy cosmic rays in astrophysical environments. This review adresses the physics of shock formation, shock dynamics and particle acceleration based on a close examination of available multi-wavelength or in situ observations, analytical and numerical developments. A particular emphasis is made on the different instabilities triggered during the shock formation and in association with particle acceleration processes with regards to the properties of the background upstream medium. It appears that among the most important parameters the background magnetic field through the magnetization and its obliquity is the dominant one. The shock velocity that can reach relativistic speeds has also a strong impact over the development of the micro-instabilities and the fate of particle acceleration. Recent developments of laboratory shock experiments has started to bring some new insights in the physics of space plasma and astrophysical shock waves. A special section is dedicated to new laser plasma experiments probing shock physics.

Spectroscopic Studies of Fermi LAT Blazars

Blazar research offers a view to one of the most energetic physical processes known to man. The high-energy end of blazar emission has been probed by the Fermi satellite mission since 2008, and it has catalogued more than a thousand gamma- ray bright blazars. However, a large fraction of these sources have no spectroscopic classification at lower energies. In this thesis, optical spectra for sixteen Fermi blazar candidates are published. The optical spectroscopic data have been observed with the Nordic Optical Telescope on the island of La Palma, Spain, during the summer of 2015. The ALFOSC instrument was used, with exposure times from 800 to 3000 seconds per target, yielding signal- to-noise ratios from 10 to 38. All of the sixteen targets show a flat, featureless optical spectrum, characteristic to BL Lacertae objects. The spectra of two targets contain faint emission features, and faint absorption features are seen in three targets. However, none of the features could be reliably identified. Therefore all of the targets are classified as BL Lacertae objects. This classification is supported by the statistical distribution of Fermi -selected active galactic nuclei; more than half of the identified Fermi AGN are BL Lacs. However, the classification of this sample could be improved further with a new observing campaign. This is especially true for the objects with uncertain spectral features.

Examining the growth and evolution of magnetic fields in clusters of galaxies: a numerical perspective

Magnetic fields are ubiquitous in galaxy cluster atmospheres and have a variety of astrophysical and cosmological consequences. Magnetic fields can contribute to the pressure support of clusters, affect thermal conduction, and modify the evolution of bubbles driven by active galactic nuclei. However, we currently do not fully understand the origin and evolution of these fields throughout cosmic time. Furthermore, we do not have a general understanding of the relationship between magnetic field strength and topology and other cluster properties, such as mass and X-ray luminosity. We can now begin to answer some of these questions using large-scale cosmological magnetohydrodynamic (MHD) simulations of the formation of galaxy clusters including the seeding and growth of magnetic fields. Using large-scale cosmological simulations with the FLASH code combined with a simplified model of the acceleration of cosmic rays responsible for the generation of radio halos, we find that the galaxy cluster frequency distribution and expected number counts of radio halos from upcoming low-frequency sur- veys are strongly dependent on the strength of magnetic fields. Thus, a more complete understanding of the origin and evolution of magnetic fields is necessary to understand and constrain models of diffuse synchrotron emission from clusters. One favored model for generating magnetic fields is through the amplification of weak seed fields in active galactic nuclei (AGN) accretion disks and their subsequent injection into cluster atmospheres via AGN-driven jets and bubbles. However, current large-scale cosmological simulations cannot directly include the physical processes associated with the accretion and feedback processes of AGN or the seeding and merging of the associated SMBHs. Thus, we must include these effects as subgrid models. In order to carefully study the growth of magnetic fields in clusters via AGN-driven outflows, we present a systematic study of SMBH and AGN subgrid models. Using dark-matter only cosmological simulations, we find that many important quantities, such as the relationship between SMBH mass and galactic bulge velocity dispersion and the merger rate of black holes, are highly sensitive to the subgrid model assumptions of SMBHs. In addition, using MHD calculations of an isolated cluster, we find that magnetic field strengths, extent, topology, and relationship to other gas quantities such as temperature and density are also highly dependent on the chosen model of accretion and feedback. We use these systematic studies of SMBHs and AGN inform and constrain our choice of subgrid models, and we use those results to outline a fully cosmological MHD simulation to study the injection and growth of magnetic fields in clusters of galaxies. This simulation will be the first to study the birth and evolution of magnetic fields using a fully closed accretion-feedback cycle, with as few assumptions as possible and a clearer understanding of the effects of the various parameter choices.

Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, 6 x 10(19) eV. The anisotropy was measured by the fraction of arrival directions that are less than 3.1 degrees from the position of an active galactic nucleus within 75 Mpc (using the Veron-Cetty and Veron 12th catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is (38(-6)(+7))%, compared with 21% expected for isotropic cosmic rays. This is down from the early estimate of (69-(+11)(13))%. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation. (C) 2010 Elsevier B.V. All rights reserved.

The Ultraviolet-Bright, Slowly Declining Transient PS1-11af as a Partial Tidal Disruption Event

We present the Pan-STARRS1 discovery of the long-lived and blue transient PS1-11af, which was also detected by Galaxy Evolution Explorer with coordinated observations in the near-ultraviolet (NUV) band. PS1-11af is associated with the nucleus of an early type galaxy at redshift z = 0.4046 that exhibits no evidence for star formation or active galactic nucleus activity. Four epochs of spectroscopy reveal a pair of transient broad absorption features in the UV on otherwise featureless spectra. Despite the superficial similarity of these features to P-Cygni absorptions of supernovae (SNe), we conclude that PS1-11af is not consistent with the properties of known types of SNe. Blackbody fits to the spectral energy distribution are inconsistent with the cooling, expanding ejecta of a SN, and the velocities of the absorption features are too high to represent material in homologous expansion near a SN photosphere. However, the constant blue colors and slow evolution of the luminosity are similar to previous optically selected tidal disruption events (TDEs). The shape of the optical light curve is consistent with models for TDEs, but the minimum accreted mass necessary to power the observed luminosity is only 0.002 M, which points to a partial disruption model. A full disruption model predicts higher bolometric luminosities, which would require most of the radiation to be emitted in a separate component at high energies where we lack observations. In addition, the observed temperature is lower than that predicted by pure accretion disk models for TDEs and requires reprocessing to a constant, lower temperature. Three deep non-detections in the radio with the Very Large Array over the first two years after the event set strict limits on the production of any relativistic outflow comparable to Swift J1644+57, even if off-axis.

PCA Tomography: how to extract information from data cubes

Astronomy has evolved almost exclusively by the use of spectroscopic and imaging techniques, operated separately. With the development of modern technologies, it is possible to obtain data cubes in which one combines both techniques simultaneously, producing images with spectral resolution. To extract information from them can be quite complex, and hence the development of new methods of data analysis is desirable. We present a method of analysis of data cube (data from single field observations, containing two spatial and one spectral dimension) that uses Principal Component Analysis (PCA) to express the data in the form of reduced dimensionality, facilitating efficient information extraction from very large data sets. PCA transforms the system of correlated coordinates into a system of uncorrelated coordinates ordered by principal components of decreasing variance. The new coordinates are referred to as eigenvectors, and the projections of the data on to these coordinates produce images we will call tomograms. The association of the tomograms (images) to eigenvectors (spectra) is important for the interpretation of both. The eigenvectors are mutually orthogonal, and this information is fundamental for their handling and interpretation. When the data cube shows objects that present uncorrelated physical phenomena, the eigenvector`s orthogonality may be instrumental in separating and identifying them. By handling eigenvectors and tomograms, one can enhance features, extract noise, compress data, extract spectra, etc. We applied the method, for illustration purpose only, to the central region of the low ionization nuclear emission region (LINER) galaxy NGC 4736, and demonstrate that it has a type 1 active nucleus, not known before. Furthermore, we show that it is displaced from the centre of its stellar bulge.

Cometary evolution : clues on physical properties from chondritic interplanetary dust particles

The degree of diversity or similarity detected in comets depends primarily on the lifetimes of the individual cometary nuclei at the time of analysis. It is inherent in our understanding of cometary orbital dynamics and the seminal model of comet origins that cometary evolution is the natural order of events in our Solar System. Thus, predictions of cometary behaviour in terms of bulk physical, mineralogical or chemical parameters should contain an appreciation of temporal variation(s). Previously, Rietmeijer and Mackinnon [1987] developed mineralogical bases for the chemical evolution of cometary nuclei primarily with regard to the predominantly silicate fraction of comet nuclei. We suggested that alteration of solids in cometary nuclei should be expected and that indications of likely reactants and products can be derived from judicious comparison with terrestrial diagenetic environments which include hydrocryogenic and low-temperature aqueous alterations. In a further development of this concept, Rietmeijer [1988] provides indirect evidence for the formation of sulfides and oxides in comet nuclei. Furthermore, Rietmeijer [1988] noted that timescales for hydrocryogenic and low-temperature reactions involving liquid water are probably adequate for relatively mature comets, e.g. P/comet Halley. In this paper, we will address the evolution of comet nuclei physical parameters such as solid particle grain size, porosity and density. In natural environments, chemical evolution (e.g. mineral reactions) is often accompanied by changes in physical properties. These concurrent changes are well-documented in the terrestrial geological literature, especially in studies of sediment diagenesis and we suggest that similar basic principles apply within the upper few meters of active comet nuclei. The database for prediction of comet nuclei physical parameters is, in principle, the same as used for the proposition of chemical evolution. We use detailed mineralogical studies of chondritic interplanetary dust particles (IDPS) as a guide to the likely constitution of mature comets traversing the inner Solar System. While there is, as yet, no direct proof that a specific sub-group or type of chondritic IDP is derived from a specific comet, it is clear that these particles are extraterrestrial in origin and that a certain portion of the interplanetary flux received by the Earth is cometary in origin. Two chondritic porous (CP) IDPS, sample numbers W701OA2 and W7029CI, from the Johnson Space Center Cosmic Dust Collection have been selected for this study of putative cometary physical parameters. This particular type of particle is considered a likely candidate for a cometary origin on the basis of mineralogy, bulk composition and morphology. While many IDPs have been subjected to intensive study over the past decade, we can develop a physical parameter model on only these two CP IDPs because few others have been studied in sufficient detail.

Dark Spots, Bubbles, and Shells in the Lobes of Extragalactic Radio Sources

Based on maps of the extragalactic radio sources Cyg A, Her A, Cen A, 3C 277.3 and others, arguments are given that the twin-jets from the respective active galactic nucleus ram their channels repeatedly through thin, massive shells. The jets are thereby temporarily choked and blow radio bubbles. Warm shell matter in the cocoon shows up radio-dark through electron-scattering.

Starburst in the ultraluminous galaxy Arp 220: Constraints from observations of radio recombination lines and continuum

We present observations of radio recombination lines (RRL) from the starburst galaxy Arp 220 at 8.1 GHz (H92 alpha) and 1.4 GHz (H167 alpha and H165 alpha) and at 84 GHz (H42 alpha), 96 GHz (H40 alpha) and 207 GHz (H31 alpha) using the Very Large Array and the IRAM 30 m telescope, respectively. RRLs were detected at all the frequencies except 1.4 GHz, where a sensitive upper limit was obtained. We also present continuum flux measurements at these frequencies as well as at 327 MHz made with the VLA. The continuum spectrum, which has a spectral index alpha similar to -0.6 (S-nu proportional to nu(alpha)) between 5 and 10 GHz, shows a break near 1.5 GHz, a prominent turnover below 500 MHz, and a flatter spectral index above 50 GHz. We show that a model with three components of ionized gas with different densities and area covering factors can consistently explain both RRL and continuum data. The total mass of ionized gas in the three components is 3.2 x 10(7) M., requiring 3 x 10(5) O5 stars with a total Lyman continuum production rate N-Lyc similar to 1.3 x 10(55) photons s(-1). The ratio of the expected to observed Br alpha and Br gamma fluxes implies a dust extinction A(V) similar to 45 mag. The derived Lyman continuum photon production rate implies a continuous star formation rate (SFR) averaged over the lifetime of OB stars of similar to 240 M yr(-1). The Lyman continuum photon Production rate of similar to 3% associated with the high-density H II regions implies a similar SFR at recent epochs (t < 10(5) yr). An alternative model of high-density gas, which cannot be excluded on the basis of the available data, predicts 10 times higher SFR at recent epochs. If confirmed, this model implies that star formation in Arp 220 consists of multiple starbursts of very high SFR (few times 10(3) M. yr(-1)) and short duration (similar to 10(5) yr). The similarity of IR excess, L-IR/L-Ly alpha similar to 24, in Arp 220 to values observed in starburst galaxies shows that most of the high luminosity of Arp 220 is due to the ongoing starburst rather than to a hidden active galactic nucleus (AGN). A comparison of the IR excesses in Arp 220, the Galaxy, and M33 indicates that the starburst in Arp 220 has an initial mass function that is similar to that in normal galaxies and has a duration longer than 107 yr. If there was no infall of gas during this period, then the star formation efficiency (SFE) in Arp 220 is similar to 50%. The high SFR and SFE in Arp 220 is consistent with their known dependences on mass and density of gas in star-forming regions of normal galaxies.

Hot gaseous coronae of early-type galaxies and their radio luminosity function

Recent X-ray observations have revealed that early-type galaxies (which usually produce extended double radio sources) generally have hot gaseous haloes extending up to approx102kpc1,2. Moreover, much of the cosmic X-ray background radiation is probably due to a hotter, but extremely tenuous, intergalactic medium (IGM)3. We have presented4–7 an analytical model for the propagation of relativistic beams from galactic nuclei, in which the beams' crossing of the pressure-matched interface between the IGM and the gaseous halo, plays an important role. The hotspots at the ends of the beams fade quickly when their advance becomes subsonic with respect to the IGM. This model has successfully predicted (for typical double radio sources) the observed8 current mean linear-size (approx2Dsime350 kpc)4,5, the observed8–11 decrease in linear-size with cosmological redshift4–6 and the slope of the linear-size versus radio luminosity10,12–14 relation6. We have also been able to predict the redshift-dependence of observed numbers and radio luminosities of giant radio galaxies7,15. Here, we extend this model to include the propagation of somewhat weaker beams. We show that the observed flattening of the local radio luminosity function (LRLF)16–20 for radio luminosity Papproximately 1024 W Hz-1 at 1 GHz can be explained without invoking ad hoc a corresponding break in the beam power function Phi(Lb), because the heads of the beams with Lb < 1025 W Hz-1 are decelerated to sonic velocity within the halo itself, which leads to a rapid decay of radio luminosity and a reduced contribution of these intrinsically weaker sources to the observed LRLF.

Seyfert Galaxies: Nuclear Radio Structure And Unification

A radio study of a carefully selected sample of 20 Seyfert galaxies that are matched in orientation-independent parameters, which are measures of intrinsic active galactic nucleus power and host galaxy properties, is presented to test the predictions of the unified scheme hypothesis. Our sample sources have core flux densities greater than 8 mJy at 5 GHz on arcsec scales due to the feasibility requirements. These simultaneous parsec-scale and kiloparsec-scale radio observations reveal (1) that Seyfert 1 and Seyfert 2 galaxies have an equal tendency to show compact radio structures on milliarcsecond scales, (2) the distributions of parsec-scale and kiloparsec-scale radio luminosities are similar for both Seyfert 1 and Seyfert 2 galaxies, (3) there is no evidence for relativistic beaming in Seyfert galaxies, (4) similar distributions of source spectral indices in spite of the fact that Seyferts show nuclear radio flux density variations, and (5) the distributions of the projected linear size for Seyfert 1 and Seyfert 2 galaxies are not significantly different as would be expected in the unified scheme. The latter could be mainly due to a relatively large spread in the intrinsic sizes. We also find that a starburst alone cannot power these radio sources. Finally, an analysis of the kiloparsec-scale radio properties of the CfA Seyfert galaxy sample shows results consistent with the predictions of the unified scheme.

Lopsidedness and sloshing in centres of advanced mergers of galaxies

We measure the non-axisymmetry in the luminosity distribution in the central few kpc of a sample of advanced mergers of galaxies, by analyzing their 2MASS images. All mergers show a high central asymmetry: the centres of isophotes show a striking sloshing pattern with a spatial variation of upto 30% within the central 1 kpc; and the Fourier amplitude for lopsidedness (m = 1) shows high values upto 0.2 within the central 5 kpc. The central asymmetry is estimated to be long-lived, lasting for ~ a few Gyr or ~ 100 local dynamical timescales. This will significantly affect the dynamical evolution of this region, by helping fuel the central active galactic nucleus, and also by causing the secular growth of the bulge driven by lopsidedness.

Unstable m=1 modes of counter-rotating Keplerian discs

We study the linear m= 1 counter-rotating instability in a two-component, nearly Keplerian disc. Our goal is to understand these slow modes in discs orbiting massive black holes in galactic nuclei. They are of interest not only because they are of large spatial scale and can hence dominate observations but also because they can be growing modes that are readily excited by accretion events. Self-gravity being non-local, the eigenvalue problem results in a pair of coupled integral equations, which we derive for a two-component softened gravity disc. We solve this integral eigenvalue problem numerically for various values of mass fraction in the counter-rotating component. The eigenvalues are in general complex, being real only in the absence of the counter-rotating component, or imaginary when both components have identical surface density profiles. Our main results are as follows: (i) the pattern speed appears to be non-negative, with the growth (or damping) rate being larger for larger values of the pattern speed; (ii) for a given value of the pattern speed, the growth (or damping) rate increases as the mass in the counter-rotating component increases; (iii) the number of nodes of the eigenfunctions decreases with increasing pattern speed and growth rate. Observations of lopsided brightness distributions would then be dominated by modes with the least number of nodes, which also possess the largest pattern speeds and growth rates.