Table of contents

Volume 214

Number 2, October 2014

Previous issue Next issue

15

, , , and

Using a compilation of 25 studies from the literature, we investigate the evolution of the star-forming galaxy (SFG) main sequence (MS) in stellar mass and star formation rate (SFR) out to z ∼ 6. After converting all observations to a common set of calibrations, we find a remarkable consensus among MS observations (∼0.1 dex 1σ interpublication scatter). By fitting for time evolution of the MS in bins of constant mass, we deconvolve the observed scatter about the MS within each observed redshift bin. After accounting for observed scatter between different SFR indicators, we find the width of the MS distribution is ∼0.2 dex and remains constant over cosmic time. Our best fits indicate the slope of the MS is likely time-dependent, with our best-fit log SFR(M*, t) = (0.84 ± 0.02 − 0.026 ± 0.003 × t)log M* − (6.51 ± 0.24 − 0.11 ± 0.03 × t), where t is the age of the universe in Gyr. We use our fits to create empirical evolutionary tracks in order to constrain MS galaxy star formation histories (SFHs), finding that (1) the most accurate representations of MS SFHs are given by delayed-τ models, (2) the decline in fractional stellar mass growth for a "typical" MS galaxy today is approximately linear for most of its lifetime, and (3) scatter about the MS can be generated by galaxies evolving along identical evolutionary tracks assuming an initial 1σ spread in formation times of ∼1.4 Gyr.

16

, , and

We propose a novel numerical method for solving multi-dimensional, special relativistic Boltzmann equations for neutrinos coupled with hydrodynamics equations. This method is meant to be applied to simulations of core-collapse supernovae. We handle special relativity in a non-conventional way, taking account of all orders of v/c. Consistent treatment of the advection and collision terms in the Boltzmann equations has been a challenge, which we overcome by employing two different energy grids: Lagrangian remapped and laboratory fixed grids. We conduct a series of basic tests and perform a one-dimensional simulation of core-collapse, bounce, and shock-stall for a 15 M progenitor model with a minimum but essential set of microphysics. We demonstrate in the latter simulation that our new code is capable of handling all phases in core-collapse supernova. For comparison, a non-relativistic simulation is also conducted with the same code, and we show that they produce qualitatively wrong results in neutrino transfer. Finally, we discuss a possible incorporation of general relativistic effects into our method.

17

, , , , , , , , and

Orbital monitoring of M-type binaries is essential for constraining their fundamental properties. This is particularly useful in young systems, where the extended pre-main-sequence evolution can allow for precise isochronal dating. Here, we present the continued astrometric monitoring of the more than 200 binaries of the AstraLux Large Multiplicity Survey, building both on our previous work, archival data, and new astrometric data spanning the range of 2010–2012. The sample is very young overall—all included stars have known X-ray emission, and a significant fraction (18%) of them have recently also been identified as members of young moving groups in the solar neighborhood. We identify ∼30 targets that both have indications of being young and for which an orbit either has been closed or appears possible to close in a reasonable time frame (a few years to a few decades). One of these cases, GJ 4326, is, however, identified as probably being substantially older than has been implied from its apparent moving group membership, based on astrometric and isochronal arguments. With further astrometric monitoring, these targets will provide a set of empirical isochrones, against which theoretical isochrones can be calibrated, and which can be used to evaluate the precise ages of nearby young moving groups.

18

, , , , and

New experimental absolute atomic transition probabilities are reported for 203 lines of V ii. Branching fractions are measured from spectra recorded using a Fourier transform spectrometer and an echelle spectrometer. The branching fractions are normalized with radiative lifetime measurements to determine the new transition probabilities. Generally good agreement is found between this work and previously reported V ii transition probabilities. Two spectrometers, independent radiometric calibration methods, and independent data analysis routines enable a reduction in systematic uncertainties, in particular those due to optical depth errors. In addition, new hyperfine structure constants are measured for selected levels by least squares fitting line profiles in the FTS spectra. The new V ii data are applied to high resolution visible and UV spectra of the Sun and metal-poor star HD 84937 to determine new, more accurate V abundances. Lines covering a range of wavelength and excitation potential are used to search for non-LTE effects. Very good agreement is found between our new solar photospheric V abundance, log ε(V) = 3.95 from 15 V ii lines, and the solar-system meteoritic value. In HD 84937, we derive [V/H] = −2.08 from 68 lines, leading to a value of [V/Fe] = 0.24.

19

, , , , , and

We have surveyed a complete extent of Leo A—an apparently isolated gas-rich low-mass dwarf irregular galaxy in the Local Group. The B, V, and I passband CCD images (typical seeing ∼0farcs8) were obtained with the Subaru Telescope equipped with the Suprime-Cam mosaic camera. The wide-field (20' × 24') photometry catalog of 38,856 objects (V ∼ 16–26 mag) is presented. This survey is also intended to serve as "a finding chart" for future imaging and spectroscopic observation programs of Leo A.

20

, , , , and

We present the J- and H-band source catalog covering the AKARI North Ecliptic Pole field. Filling the gap between the optical data from other follow-up observations and mid-infrared (MIR) data from AKARI, our near-infrared (NIR) data provides contiguous wavelength coverage from optical to MIR. For the J- and H-band imaging, we used the FLoridA Multi-object Imaging Near-ir Grism Observational Spectrometer on the Kitt Peak National Observatory 2.1 m telescope covering a 5.1 deg2 area down to a 5σ depth of ∼21.6 mag and ∼21.3 mag (AB) for the J and H bands with an astrometric accuracy of 0farcs14 and 0farcs17 for 1σ in R.A. and decl. directions, respectively. We detected 208,020 sources for the J band and 203,832 sources for the H band. This NIR data is being used for studies including the analysis of the physical properties of infrared sources such as stellar mass and photometric redshifts, and will be a valuable data set for various future missions.

21

, , , , and

Six low-mass embedded sources (L1489, L1551-IRS5, TMR1, TMC1-A, L1527, and TMC1) in Taurus have been observed with Herschel-PACS to cover the full spectrum from 50 to 210 μm as part of the Herschel key program, "Dust, Ice, and Gas In Time." The relatively low intensity of the interstellar radiation field surrounding Taurus minimizes contamination of the [C ii] emission associated with the sources by diffuse emission from the cloud surface, allowing study of the [C ii] emission from the source. In several sources, the [C ii] emission is distributed along the outflow, as is the [O i] emission. The atomic line luminosities correlate well with each other, as do the molecular lines, but the atomic and molecular lines correlate poorly. The relative contribution of CO to the total gas cooling is constant at ∼30%, while the cooling fraction by H2O varies from source to source, suggesting different shock properties resulting in different photodissociation levels of H2O. The gas with a power-law temperature distribution with a moderately high density can reproduce the observed CO fluxes, indicative of CO close to LTE. However, H2O is mostly subthermally excited. L1551-IRS5 is the most luminous source (Łbol = 24.5 L) and the [O i] 63.1 μm line accounts for more than 70% of its FIR line luminosity, suggesting complete photodissociation of H2O by a J shock. In L1551-IRS5, the central velocity shifts of the [O i] line, which exceed the wavelength calibration uncertainty (∼70 km s−1) of PACS, are consistent with the known redshifted and blueshifted outflow direction.

22

, , and

We develop new hyperon equation of state (EoS) tables for core-collapse supernova simulations and neutron stars. These EoS tables are based on a density-dependent relativistic hadron field theory where baryon–baryon interaction is mediated by mesons, using the parameter set DD2 for nucleons. Furthermore, light and heavy nuclei along with interacting nucleons are treated in the nuclear statistical equilibrium model of Hempel and Schaffner-Bielich which includes excluded volume effects. Of all possible hyperons, we consider only the contribution of Λs. We have developed two variants of hyperonic EoS tables: in the npΛϕ case the repulsive hyperon–hyperon interaction mediated by the strange ϕ meson is taken into account, and in the npΛ case it is not. The EoS tables for the two cases encompass a wide range of densities (10−12 to ∼1 fm−3), temperatures (0.1 to 158.48 MeV), and proton fractions (0.01 to 0.60). The effects of Λ hyperons on thermodynamic quantities such as free energy per baryon, pressure, or entropy per baryon are investigated and found to be significant at higher densities. The cold, β-equilibrated EoS (with the crust included self-consistently) results in a 2.1 M maximum mass neutron star for the npΛϕ case, whereas that for the npΛ case is 1.95 M. The npΛϕ EoS represents the first supernova EoS table involving hyperons that is directly compatible with the recently measured 2 M neutron stars.

23

, , , , and

As a step toward a comprehensive overview of the infrared (IR) diagnostics of the central engines and host galaxies of quasars at low redshift, we present SpitzerSpace Telescope spectroscopic (5–40 μm) and photometric (24, 70, and 160 μm) measurements of all Palomar–Green (PG) quasars at z < 0.5 and Two Micron All Sky Survey (2MASS) quasars at z < 0.3. We supplement these data with Herschel measurements at 160 μm. The sample is composed of 87 optically selected PG quasars and 52 near-IR-selected 2MASS quasars. Here we present the data, measure the prominent spectral features, and separate emission due to star formation from that emitted by the dusty circumnuclear torus. We find that the mid-IR (5–30 μm) spectral shape for the torus is largely independent of quasar IR luminosity with scatter in the spectral energy distribution (SED) shape of ≲0.2 dex. Except for the silicate features, no large difference is observed between PG (unobscured—silicate emission) and 2MASS (obscured—silicate absorption) quasars. Only mild silicate features are observed in both cases. When in emission, the peak wavelength of the silicate feature tends to be longer than 9.7 μm, possibly indicating effects on grain properties near the active galactic nucleus. The IR color is shown to correlate with the equivalent width of the aromatic features, indicating that the slope of the quasar mid- to far-IR SED is to first order driven by the fraction of radiation from star formation in the IR bands.

24

, , , , , , , , , et al

The 3D-HST and CANDELS programs have provided WFC3 and ACS spectroscopy and photometry over ≈900 arcmin2 in five fields: AEGIS, COSMOS, GOODS-North, GOODS-South, and the UKIDSS UDS field. All these fields have a wealth of publicly available imaging data sets in addition to the Hubble Space Telescope(HST) data, which makes it possible to construct the spectral energy distributions (SEDs) of objects over a wide wavelength range. In this paper we describe a photometric analysis of the CANDELS and 3D-HST HST imaging and the ancillary imaging data at wavelengths 0.3–8 μm. Objects were selected in the WFC3 near-IR bands, and their SEDs were determined by carefully taking the effects of the point-spread function in each observation into account. A total of 147 distinct imaging data sets were used in the analysis. The photometry is made available in the form of six catalogs: one for each field, as well as a master catalog containing all objects in the entire survey. We also provide derived data products: photometric redshifts, determined with the EAZY code, and stellar population parameters determined with the FAST code. We make all the imaging data that were used in the analysis available, including our reductions of the WFC3 imaging in all five fields. 3D-HST is a spectroscopic survey with the WFC3 and ACS grisms, and the photometric catalogs presented here constitute a necessary first step in the analysis of these grism data. All the data presented in this paper are available through the 3D-HST Web site (http://3dhst.research.yale.edu).

25

, , , , , and

We present new calculations of Rosseland and Planck gaseous mean opacities relevant to the atmospheres of giant planets and ultracool dwarfs. Such calculations are used in modeling the atmospheres, interiors, formation, and evolution of these objects. Our calculations are an expansion of those presented in Freedman et al. to include lower pressures, finer temperature resolution, and also the higher metallicities most relevant for giant planet atmospheres. Calculations span 1 μbar to 300 bar, and 75–4000 K, in a nearly square grid. Opacities at metallicities from solar to 50 times solar abundances are calculated. We also provide an analytic fit to the Rosseland mean opacities over the grid in pressure, temperature, and metallicity. In addition to computing mean opacities at these local temperatures, we also calculate them with weighting functions up to 7000 K, to simulate the mean opacities for incident stellar intensities, rather than locally thermally emitted intensities. The chemical equilibrium calculations account for the settling of condensates in a gravitational field and are applicable to cloud-free giant planet and ultracool dwarf atmospheres, but not circumstellar disks. We provide our extensive opacity tables for public use.

26

, , , , and

New red and violet system line lists for the CN isotopologues 13C14N and 12C15N have been generated. These new transition data are combined with those previously derived for 12C14N, and applied to the determination of CNO abundances in the solar photosphere and in four red giant stars: Arcturus, the bright, very low-metallicity star HD 122563, and the carbon-enhanced metal-poor stars HD 196944 and HD 201626. When both red and violet system lines are detectable in a star, their derived N abundances are in good agreement. The mean N abundances determined in this work are also generally in accord with published values.

27

, , , , , , , , , et al

Recently the number of main-sequence and subgiant stars exhibiting solar-like oscillations that are resolved into individual mode frequencies has increased dramatically. While only a few such data sets were available for detailed modeling just a decade ago, the Kepler mission has produced suitable observations for hundreds of new targets. This rapid expansion in observational capacity has been accompanied by a shift in analysis and modeling strategies to yield uniform sets of derived stellar properties more quickly and easily. We use previously published asteroseismic and spectroscopic data sets to provide a uniform analysis of 42 solar-type Kepler targets from the Asteroseismic Modeling Portal. We find that fitting the individual frequencies typically doubles the precision of the asteroseismic radius, mass, and age compared to grid-based modeling of the global oscillation properties, and improves the precision of the radius and mass by about a factor of three over empirical scaling relations. We demonstrate the utility of the derived properties with several applications.

28

and

Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) are methods used to stochastically solve the radiative transport and diffusion equations, respectively. These methods combine into a hybrid transport-diffusion method we refer to as IMC–DDMC. We explore a multigroup IMC–DDMC scheme that in DDMC, combines frequency groups with sufficient optical thickness. We term this procedure "opacity regrouping." Opacity regrouping has previously been applied to IMC–DDMC calculations for problems in which the dependence of the opacity on frequency is monotonic. We generalize opacity regrouping to non-contiguous groups and implement this in SuperNu, a code designed to do radiation transport in high-velocity outflows with non-monotonic opacities. We find that regrouping of non-contiguous opacity groups generally improves the speed of IMC–DDMC radiation transport. We present an asymptotic analysis that informs the nature of the Doppler shift in DDMC groups and summarize the derivation of the Gentile–Fleck factor for modified IMC–DDMC. We test SuperNu using numerical experiments including a quasi-manufactured analytic solution, a simple 10 group problem, and the W7 problem for Type Ia supernovae. We find that opacity regrouping is necessary to make our IMC–DDMC implementation feasible for the W7 problem and possibly Type Ia supernova simulations in general. We compare the bolometric light curves and spectra produced by the SuperNu and PHOENIX radiation transport codes for the W7 problem. The overall shape of the bolometric light curves are in good agreement, as are the spectra and their evolution with time. However, for the numerical specifications we considered, we find that the peak luminosity of the light curve calculated using SuperNu is ∼10% less than that calculated using PHOENIX.