Table of contents

The anti-glitch of magnetar 1E 2259+586 is analyzed theoretically. An enhanced particle wind during the observational interval takes away additional rotational energy of the neutron star, which will result in a net spin-down of the magnetar, i.e., an anti-glitch. In the wind braking scenario of the anti-glitch, there are several predictions: (1) a radiative event will always accompany the anti-glitch, (2) there will be a decrease/variation of the braking index after the anti-glitch, and (3) the anti-glitch is just a period of enhanced spin-down. If there are enough timing observations, a period of enhanced spin-down is expected instead of an anti-glitch. Applications to previous timing events of SGR 1900+14 and PSR J1846−0258 are also included. It is shown that current timing events of 1E 2259+586, SGR 1900+14, and PSR J1846−0258 can be understood safely in the wind braking model. The enhanced spin-down and absence of an anti-glitch before the giant flare of SGR 1806−20 is consistent with the wind braking scenario.

During a stellar tidal disruption event (TDE), an accretion disk forms as stellar debris returns to the disruption site and circularizes. Rather than being confined within the circularizing radius, the disk can spread to larger radii to conserve angular momentum. A spreading disk is a source of matter for re-accretion at rates that may exceed the later stellar fallback rate, although a disk wind can suppress its contribution to the central black hole accretion rate. A spreading disk is detectible through a break in the central accretion rate history or, at longer wavelengths, by its own emission. We model the evolution of TDE disk size and accretion rate by accounting for the time-dependent fallback rate, for the influence of wind losses in the early advective stage, and for the possibility of thermal instability for accretion rates intermediate between the advection-dominated and gas-pressure-dominated states. The model provides a dynamic basis for modeling TDE light curves. All or part of a young TDE disk will precess as a solid body because of the Lense–Thirring effect, and precession may manifest itself as a quasi-periodic modulation of the light curve. The precession period increases with time. Applying our results to the jetted TDE candidate Swift J1644+57, whose X-ray light curve shows numerous quasi-periodic dips, we argue that the data best fit a scenario in which a main-sequence star was fully disrupted by an intermediate mass black hole on an orbit significantly inclined from the black hole equator, with the apparent jet shutoff at t = 500 days corresponding to a disk transition from the advective state to the gas-pressure-dominated state.

The universe is expected to be permeated by a stochastic background of gravitational radiation of astrophysical and cosmological origin. This background is capable of exciting oscillations in solar-like stars. Here we show that solar-like oscillators can be employed as giant hydrodynamical detectors for such a background in the μHz to mHz frequency range, which has remained essentially unexplored until today. We demonstrate this approach by using high-precision radial velocity data for the Sun to constrain the normalized energy density of the stochastic gravitational-wave background around 0.11 mHz. These results open up the possibility for asteroseismic missions like CoRoT and Kepler to probe fundamental physics.

In the heliosheath beyond the termination shock, low energy (<0.5 keV) neutral atoms are created by charge exchange with interstellar neutrals. Detecting these neutrals from Earth's orbit is difficult because their flux is reduced substantially by ionization losses as they propagate from about 100 to 1 AU and because there are a variety of other signals and backgrounds that compete with this weak signal. Observations from IBEX-Lo and -Hi from two opposing vantage points in Earth's orbit established a lower energy limit of about 0.1 keV on measurements of energetic neutral atoms (ENAs) from the heliosphere and the form of the energy spectrum from about 0.1 to 6 keV in two directions in the sky. Below 0.1 keV, the detailed ENA spectrum is not known, and IBEX provides only upper limits on the fluxes. However, using some assumptions and taking constraints on the spectrum into account, we find indications that the spectrum turns over at an energy between 0.1 and 0.2 keV.

We present a 4 deg² weak gravitational lensing survey of subhalos in the very nearby Coma cluster using the Subaru/Suprime-Cam. The large apparent size of cluster subhalos allows us to measure the mass of 32 subhalos detected in a model-independent manner, down to the order of 10⁻³ of the virial mass of the cluster. Weak-lensing mass measurements of these shear-selected subhalos enable us to investigate subhalo properties and the correlation between subhalo masses and galaxy luminosities for the first time. The mean distortion profiles stacked over subhalos show a sharply truncated feature which is well-fitted by a Navarro–Frenk–White (NFW) mass model with the truncation radius, as expected due to tidal destruction by the main cluster. We also found that subhalo masses, truncation radii, and mass-to-light ratios decrease toward the cluster center. The subhalo mass function, dn/dln M_sub, in the range of 2 orders of magnitude in mass, is well described by a single power law or a Schechter function. Best-fit power indices of $1.09^{+0.42}_{-0.32}$ for the former model and $0.99_{-0.23}^{+0.34}$ for the latter, are in remarkable agreement with slopes of ∼0.9–1.0 predicted by the cold dark matter paradigm. The tangential distortion signals in the radial range of 0.02–2 h⁻¹ Mpc from the cluster center show a complex structure which is well described by a composition of three mass components of subhalos, the NFW mass distribution as a smooth component of the main cluster, and a lensing model from a large scale structure behind the cluster. Although the lensing signals are 1 order of magnitude lower than those for clusters at z ∼ 0.2, the total signal-to-noise ratio, S/N = 13.3, is comparable, or higher, because the enormous number of background source galaxies compensates for the low lensing efficiency of the nearby cluster.

A bright ray that developed in the wake of a fast coronal mass ejection (CME) on 2005 September 7 presents a unique opportunity to study the early development and physical characteristics of a reconnecting current sheet (CS). Polarization brightness images from the Mk4 K-Coronameter at the Mauna Loa Solar Observatory are used to determine the structure of the ray along its axis low in the corona as it progressed outward. Coverage of the early development of the ray out to ∼1.3 R_☉ for a period of ∼27 hr after the start of the event enables for the first time in white light a measurement of a CME CS from the top of the arcade to the base of the flux rope. Measured widths of the ray are combined to obtain the kinematics of the upper and lower "Y"-points described in reconnection flux-rope models such as that of Lin & Forbes. The time dependence of these points are used to derive values for the speed and acceleration of the growth of the CS. We note the appearance of a large structure which increases in size as it expands outward in the early development of the ray and an apparent oscillation with a period of ∼0.5 hr in the position angle of the ray.

We measure quasar variability using the Panoramic Survey Telescope and Rapid Response System 1 Survey (Pan-STARRS1 or PS1) and the Sloan Digital Sky Survey (SDSS) and establish a method of selecting quasars via their variability in 10⁴ deg² surveys. We use 10⁵ spectroscopically confirmed quasars that have been well measured in both PS1 and SDSS and take advantage of the decadal timescales that separate SDSS measurements and PS1 measurements. A power law model fits the data well over the entire time range tested, 0.01–10 yr. Variability in the current PS1–SDSS data set can efficiently distinguish between quasars and nonvarying objects. It improves the purity of a griz quasar color cut from 4.1% to 48% while maintaining 67% completeness. Variability will be very effective at finding quasars in data sets with no u band and in redshift ranges where exclusively photometric selection is not efficient. We show that quasars' rest-frame ensemble variability, measured as a root mean squared in Δ magnitudes, is consistent with V(z, L, t) = A₀(1 + z)^0.37(L/L₀)^−0.16(t/1 yr)^0.246, where L₀ = 10⁴⁶ erg s⁻¹ and A₀ = 0.190, 0.162, 0.147, or 0.141 in the g_P1, r_P1, i_P1, or z_P1filter, respectively. We also fit across all four filters and obtain median variability as a function of z, L, and λ as V(z, L, λ, t) = 0.079(1 + z)^0.15(L/L₀)^−0.2(λ/1000 nm)^−0.44(t/1 yr)^0.246.

Galactic superwinds may be driven by very hot outflows generated by overlapping supernovae within the host galaxy. We use the Chevalier & Clegg (CC85) wind model and the observed correlation between X-ray luminosities of galaxies and their star formation rates (SFRs) to constrain the mass-loss rates ($\dot{M}_{\rm hot}$) across a wide range of SFRs, from dwarf starbursts to ultraluminous infrared galaxies. We show that for fixed thermalization and mass-loading efficiencies, the X-ray luminosity of the hot wind scales as L_X∝SFR², significantly steeper than is observed for star-forming galaxies: L_X∝SFR. Using this difference, we constrain the mass-loading and thermalization efficiency of hot galactic winds. For reasonable values of the thermalization efficiency (≲ 1) and for SFR ≳ 10 M_☉ yr⁻¹ we find that $\dot{M}_{\rm hot}/{\rm SFR}\lesssim \,1$, which is significantly lower than required by integrated constraints on the efficiency of stellar feedback in galaxies and potentially too low to explain observations of winds from rapidly star-forming galaxies. In addition, we highlight the fact that heavily mass-loaded winds cannot be described by the adiabatic CC85 model because they become strongly radiative.

Turbulent motions are essential to the mixing of entrained fluids and are also capable of amplifying weak initial magnetic fields by small-scale dynamo action. Here we perform a systematic study of turbulent mixing in magnetized media, using three-dimensional magnetohydrodynamic simulations that include a scalar concentration field. We focus on how mixing depends on the magnetic Prandtl number, Pm, from 1 to 4 and the Mach number, ${\mathcal {M}},$ from 0.3 to 2.4. For all subsonic flows, we find that the velocity power spectrum has a k^−5/3 slope in the early kinematic phase, but steepens due to magnetic back reactions as the field saturates. The scalar power spectrum, on the other hand, flattens compared to k^−5/3 at late times, consistent with the Obukohov–Corrsin picture of mixing as a cascade process. At higher Mach numbers, the velocity power spectrum also steepens due to the presence of shocks, and the scalar power spectrum again flattens accordingly. Scalar structures are more intermittent than velocity structures in subsonic turbulence, whereas for supersonic turbulence, velocity structures appear more intermittent than the scalars only in the kinematic phase. Independent of the Mach number of the flow, scalar structures are arranged in sheets in both the kinematic and saturated phases of the magnetic field evolution. For subsonic turbulence, scalar dissipation is hindered in the strong magnetic field regions, probably due to Lorentz forces suppressing the buildup of scalar gradients, whereas for supersonic turbulence, scalar dissipation increases monotonically with increasing magnetic field strength. At all Mach numbers, mixing is significantly slowed by the presence of dynamically important small-scale magnetic fields, implying that mixing in the interstellar medium and in galaxy clusters is less efficient than modeled in hydrodynamic simulations.

Direct N-body simulations of globular clusters in a realistic Milky-Way-like potential are carried out using the code NBODY6 to determine the impact of the host galaxy disk mass and geometry on the survival of star clusters. A relation between disk mass and star-cluster dissolution timescale is derived. These N-body models show that doubling the mass of the disk from 5 × 10¹⁰ M_☉ to 10 × 10¹⁰ M_☉ halves the dissolution time of a satellite star cluster orbiting the host galaxy at 6 kpc from the galactic center. Different geometries in a disk of identical mass can determine either the survival or dissolution of a star cluster orbiting within the inner 6 kpc of the galactic center. Furthermore, disk geometry has measurable effects on the mass loss of star clusters up to 15 kpc from the galactic center. N-body simulations performed with a fine output time step show that at each disk crossing the outer layers of star clusters experiences an increase in velocity dispersion of ∼5% of the average velocity dispersion in the outer section of star clusters. This leads to an enhancement of mass loss—a clearly discernable effect of disk shocking. By running models with different inclinations, we determine that star clusters with an orbit that is perpendicular to the Galactic plane have larger mass loss rates than do clusters that evolve in the Galactic plane or in an inclined orbit.

For most planets in the range of radii from 1 to 4 R_⊕, water is a major component of the interior composition. At high pressure H₂O can be solid, but for larger planets, like Neptune, the temperature can be too high for this. Mass and age play a role in determining the transition between solid and fluid (and mixed) water-rich super-Earth. We use the latest high-pressure and ultra-high-pressure phase diagrams of H₂O, and by comparing them with the interior adiabats of various planet models, the temperature evolution of the planet interior is shown, especially for the state of H₂O. It turns out that the bulk of H₂O in a planet's interior may exist in various states such as plasma, superionic, ionic, Ice VII, Ice X, etc., depending on the size, age, and cooling rate of the planet. Different regions of the mass–radius phase space are also identified to correspond to different planet structures. In general, super-Earth-size planets (isolated or without significant parent star irradiation effects) older than about 3 Gyr would be mostly solid.

We continue our exploration of the collective properties of neutron-star X-ray binaries in the stellar fields (i.e., outside globular clusters) of normal galaxies. In Paper I of this series, we considered high-mass X-ray binaries (HMXBs). In this paper (Paper II), we consider low-mass X-ray binaries (LMXBs), whose evolutionary scenario is very different from that of HMXBs. We consider the evolution of primordial binaries up to the stage where the neutron star just formed in the supernova explosion of the primary is in a binary with its low-mass, unevolved companion, and this binary has circularized tidally, producing what we call a pre-low-mass X-ray binary (pre-LMXB). We study the constraints on the formation of such pre-LMXBs in detail (since these are low-probability events), and calculate their collective properties and formation rates. To this end, we first consider the changes in the binary parameters in the various steps involved, viz., the common-envelope phase, the supernova, and the tidal evolution. This naturally leads to a clarification of the constraints. We then describe our calculation of the evolution of the distributions of primordial binary parameters into those of pre-LMXB parameters, following the standard evolutionary scenario for individual binaries. We display the latter as both bivariate and monovariate distributions, discuss their essential properties, and indicate the influences of some essential factors on these. Finally, we calculate the formation rate of these pre-LMXBs. The results of this paper will be used in a subsequent one to compute the expected X-ray luminosity function of LMXBs.

There are a number of key structures that can be used to reveal the formation and modification history of organic matter in the cosmos. For instance, the susceptibility of organic matter to heat is well documented and the relative thermal stabilities of different isomers can be used as cosmothermometers. Yet despite being an important variable, no previously recognized organic marker of pressure exists. The absence of a pressure marker is unfortunate considering our ability to effectively recognize extraterrestrial organic structures both remotely and in the laboratory. There are a wide variety of pressures in cosmic settings that could potentially be reflected by organic structures. Therefore, to develop an organic cosmic pressure marker, we have used state-of-the-art diamond anvil cell (DAC) and synchrotron-source Fourier transform infrared (FTIR) spectroscopy to reveal the effects of pressure on the substitution patterns for representatives of the commonly encountered methyl substituted naphthalenes, specifically the dimethylnaphthalenes. Interestingly, although temperature and pressure effects are concordant for many isomers, pressure appears to have the opposite effect to heat on the final molecular architecture of the 1,5-dimethylnaphthalene isomer. Our data suggest the possibility of the first pressure parameter or "cosmo-barometer" (1,5-dimethylnaphthalene/total dimethylnaphthalenes) that can distinguish pressure from thermal effects. Information can be obtained from the new pressure marker either remotely by instrumentation on landers or rovers or directly by laboratory measurement, and its use has relevance for all cases where organic matter, temperature, and pressure interplay in the cosmos.

We present the search for the [C ii] emission line in three z > 6.5 Lyα emitters (LAEs) and one J-dropout galaxy using the Combined Array for Research in Millimeter-wave Astronomy and the Plateau de Bure Interferometer. We observed three bright z ∼ 6.5–7 LAEs discovered in the Subaru Deep Field (SDF) and the multiple imaged lensed z ∼ 11 galaxy candidate found behind the galaxy cluster MACSJ0647.7+7015. For the LAEs IOK-1 (z = 6.965), SDF J132415.7+273058 (z = 6.541), and SDF J132408.3+271543 (z = 6.554) we find upper limits for the [C ii] line luminosity of <2.05, <4.52, and <10.56 × 10⁸ L_☉, respectively. We find upper limits to the far-IR (FIR) luminosity of the galaxies using a spectral energy distribution template of the local galaxy NGC 6946 and taking into account the effects of the cosmic microwave background on the millimeter observations. For IOK-1, SDF J132415.7+273058, and SDF J132408.3+271543 we find upper limits for the FIR luminosity of <2.33, 3.79, and 7.72 × 10¹¹ L_☉, respectively. For the lensed galaxy MACS0647-JD, one of the highest-redshift galaxy candidates to date with $z_{\rm ph}=10.7^{+0.6}_{-0.4}$, we put an upper limit in the [C ii] emission of <1.36 × 10⁸ × (μ/15)⁻¹ L_☉ and an upper limit in the FIR luminosity of <6.1 × 10¹⁰ × (μ/15)⁻¹ L_☉ (where μ is the magnification factor). We explore the different conditions relevant for the search for [C ii] emission in high-redshift galaxies as well as the difficulties for future observations with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Cerro Chajnantor Atacama Telescope (CCAT).

Sub-parsec binary massive black holes (BBHs) have long been thought to exist in many QSOs but remain observationally elusive. In this paper, we propose a novel method to probe sub-parsec BBHs through microlensing of lensed QSOs. If a QSO hosts a sub-parsec BBH in its center, it is expected that the BBH is surrounded by a circumbinary disk, each component of the BBH is surrounded by a small accretion disk, and a gap is opened by the secondary component in between the circumbinary disk and the two small disks. Assuming such a BBH structure, we generate mock microlensing light curves for some QSO systems that host BBHs with typical physical parameters. We show that microlensing light curves of a BBH QSO system at the infrared–optical–UV bands can be significantly different from those of corresponding QSO system with a single massive black hole (MBH), mainly because of the existence of the gap and the rotation of the BBH (and its associated small disks) around the center of mass. We estimate the half-light radii of the emission region at different wavelengths from mock light curves and find that the obtained half-light radius versus wavelength relations of BBH QSO systems can be much flatter than those of single MBH QSO systems at a wavelength range determined by the BBH parameters, such as the total mass, mass ratio, separation, accretion rates, etc. The difference is primarily due to the existence of the gap. Such unique features on the light curves and half-light radius–wavelength relations of BBH QSO systems can be used to select and probe sub-parsec BBHs in a large number of lensed QSOs to be discovered by current and future surveys, including the Panoramic Survey Telescope and Rapid Response System, the Large Synoptic Survey telescope, and Euclid.

So far, the straight flux tube model proposed by Edwin & Roberts is the most commonly used tool in practical coronal seismology, in particular, to infer values of the (coronal) magnetic field from observed, standing kink mode oscillations. In this paper, we compare the period predicted by this basic model with three-dimensional (3D) numerical simulations of standing kink mode oscillations, as the period is a crucial parameter in the seismological inversion to determine the magnetic field. We perform numerical simulations of standing kink modes in both straight and curved 3D coronal loops and consider excitation by internal and external drivers. The period of oscillation for the displacement of dense coronal loops is determined by the loop length and the kink speed, in agreement with the estimate based on analytical theory for straight flux tubes. For curved coronal loops embedded in a magnetic arcade and excited by an external driver, a secondary mode with a period determined by the loop length and external Alfvén speed is also present. When a low number of oscillations is considered, these two periods can result in a single, non-resolved (broad) peak in the power spectrum, particularly for low values of the density contrast for which the two periods will be relatively similar. In that case (and for this particular geometry), the presence of this additional mode would lead to ambiguous seismological estimates of the magnetic field strength.

In this work, UV and white-light (WL) coronagraphic data are combined to derive the full set of plasma physical parameters along the front of a shock driven by a coronal mass ejection. Pre-shock plasma density, shock compression ratio, speed, and inclination angle are estimated from WL data, while pre-shock plasma temperature and outflow velocity are derived from UV data. The Rankine-Hugoniot (RH) equations for the general case of an oblique shock are then applied at three points along the front located between 2.2 and 2.6 R_☉ at the shock nose and at the two flanks. Stronger field deflection (by ∼46°), plasma compression (factor ∼2.7), and heating (factor ∼12) occur at the nose, while heating at the flanks is more moderate (factor 1.5–3.0). Starting from a pre-shock corona where protons and electrons have about the same temperature (T_p ∼ T_e ∼ 1.5 × 10⁶ K), temperature increases derived with RH equations could better represent the proton heating (by dissipation across the shock), while the temperature increase implied by adiabatic compression (factor ∼2 at the nose, ∼1.2–1.5 at the flanks) could be more representative of electron heating: the transit of the shock causes a decoupling between electron and proton temperatures. Derived magnetic field vector rotations imply a draping of field lines around the expanding flux rope. The shock turns out to be super-critical (sub-critical) at the nose (at the flanks), where derived post-shock plasma parameters can be very well approximated with those derived by assuming a parallel (perpendicular) shock.

We present an alternative model of coronal-loop oscillations, which considers that the waves are trapped in a two-dimensional waveguide formed by the entire arcade of field lines. This differs from the standard one-dimensional model which treats the waves as the resonant oscillations of just the visible bundle of field lines. Within the framework of our two-dimensional model, the two types of oscillations that have been observationally identified, flare-induced waves and "decayless" oscillations, can both be attributed to MHD fast waves. The two components of the signal differ only because of the duration and spatial extent of the source that creates them. The flare-induced waves are generated by strong localized sources of short duration, while the decayless background can be excited by a continuous, stochastic source. Further, the oscillatory signal arising from a localized, short-duration source can be interpreted as a pattern of interference fringes produced by waves that have traveled diverse routes of various pathlengths through the waveguide. The resulting amplitude of the fringes slowly decays in time with an inverse square root dependence. The details of the interference pattern depend on the shape of the arcade and the spatial variation of the Alfvén speed. The rapid decay of this wave component, which has previously been attributed to physical damping mechanisms that remove energy from resonant oscillations, occurs as a natural consequence of the interference process without the need for local dissipation.

Since exoplanets were detected using the radial velocity method, they have revealed a diverse distribution of orbital configurations. Among these are planets in highly eccentric orbits (e > 0.5). Most of these systems consist of a single planet but several have been found to also contain a longer period planet in a near-circular orbit. Here we use the latest Keplerian orbital solutions to investigate four known systems which exhibit this extreme eccentricity diversity; HD 37605, HD 74156, HD 163607, and HD 168443. We place limits on the presence of additional planets in these systems based on the radial velocity residuals. We show that the two known planets in each system exchange angular momentum through secular oscillations of their eccentricities. We calculate the amplitude and timescale for these eccentricity oscillations and associated periastron precession. We further demonstrate the effect of mutual orbital inclinations on the amplitude of high-frequency eccentricity oscillations. Finally, we discuss the implications of these oscillations in the context of possible origin scenarios for unequal eccentricities.

We present 13 Type Ia supernovae (SNe Ia) observed in the rest-frame near-infrared (NIR) from 0.02 < z < 0.09 with the WIYN High-resolution Infrared Camera on the WIYN 3.5 m telescope. With only one to three points per light curve and a prior on the time of maximum from the spectrum used to type the object, we measure an H-band dispersion of spectroscopically normal SNe Ia of 0.164 mag. These observations continue to demonstrate the improved standard brightness of SNe Ia in an H band, even with limited data. Our sample includes two SNe Ia at z ∼ 0.09, which represent the most distant rest-frame NIR H-band observations published to date. This modest sample of 13 NIR SNe Ia represent the pilot sample for "SweetSpot"—a 3 yr NOAO Survey program that will observe 144 SNe Ia in the smooth Hubble flow. By the end of the survey we will have measured the relative distance to a redshift of z ∼ 0.05%–1%. Nearby Type Ia supernova (SN Ia) observations such as these will test the standard nature of SNe Ia in the rest-frame NIR, allow insight into the nature of dust, and provide a critical anchor for future cosmological SN Ia surveys at higher redshift.

Pulsars, if existing and detectable in the immediate vicinity of the massive black hole (MBH) in the Galactic center (GC), may be used as a superb tool to probe both the environment and the metric of the central MBH. The recent discovery of a magnetized pulsar in the GC suggests that many more pulsars should exist near the MBH. In this paper, we estimate the number and the orbital distribution of pulsars in the vicinity of the MBH in the GC by assuming that the pulsar progenitors, similar to the GC S-stars, were captured to orbits tightly bound to the MBH through the tidal breakup of stellar binaries. We use the current observations on both the GC S-stars and the hypervelocity stars to calibrate the injection rate(s) of and the dynamical model(s) for the stellar binaries. By including the relaxation processes, supernova kicks, and gravitational wave radiation in our simulations, we estimate that ∼97–190 (9–14) pulsars may presently orbit the central MBH with semimajor axes ⩽4000 AU (⩽1000 AU), which is compatible with the current observational constraints on the number of the GC pulsars. The semimajor axis and the pericenter distance of the pulsar closest to the central MBH are probably in the range of ∼120–460 AU and ∼2–230 AU, respectively. Future telescopes, such as the Square Kilometer Array, may be able to detect a significant number of pulsars with semimajor axis smaller than a few thousand AU in the GC. Long-term monitoring of these pulsars would be helpful in constraining both the environment and the metric of the central MBH. Our preferred model also results in about ten hyperfast pulsars with velocity ≳ 1500 km s⁻¹ moving away from the Milky Way.

The formation of hypercompact (HC) H ii regions is an important stage in massive star formation. Spectral line and continuum observations can explore its dynamic conditions. We present high angular resolution observations carried out with the Submillimeter Array (SMA) and the Very Large Array (VLA) toward the HC H ii region G35.58-0.03. With the 1.3 mm SMA and 1.3 cm VLA, we detected a total of about 25 transitions of 8 different species and their isotopologues (CO, CH₃CN, SO₂, CH₃CCH, OCS, CS, H, and NH₃). G35.58-0.03 consists of an HC H ii core with electron temperature $T_e^*\geqq 5500$ K, emission measure EM ≈1.9 × 10⁹ pc cm⁻⁶, local volume electron density n_e = 3.3 × 10⁵ cm⁻³, and a same width of radio recombination line FWHM ≈ 43.2 km s⁻¹ for both H30α and H38β at its intrinsic core size ∼3714 AU. The H30α line shows evidence of an ionized outflow driving a molecular outflow. Based on the derived Lyman continuum flux, there should be an early-type star equivalent to O6.5 located inside the H ii region. From the continuum spectral energy distribution from 3.6 cm, 2.0 cm, 1.3 cm, 1.3 mm and 0.85 mm to 0.45 mm, we distinguished the free–free emission (25% ∼ 55%) from the warm dust component (75% ∼ 45%) at 1.3 mm. The molecular envelope shows evidence of infall and outflow with an infall rate of 0.05 M_☉ yr⁻¹ and a mass loss rate of 5.2 × 10⁻³ M_☉ yr⁻¹. The derived momentum (∼0.05 M_☉ km s⁻¹) is consistent between the infalling and outflowing gas per year. It is suggested that the infall is predominant and the envelope mass of dense core is increasing rapidly, but the accretion in the inner part might already be halted.

We present joint constraints on the distribution of Mg ii absorption around high redshift galaxies obtained by combining two orthogonal probes, the integrated Mg ii absorption seen in stacked background galaxy spectra and the distribution of parent galaxies of individual strong Mg ii systems as seen in the spectra of background quasars. We present a suite of models that can be used to predict, for different two- and three-dimensional distributions, how the projected Mg ii absorption will depend on a galaxy's apparent inclination, the impact parameter b and the azimuthal angle between the projected vector to the line of sight and the projected minor axis. In general, we find that variations in the absorption strength with azimuthal angles provide much stronger constraints on the intrinsic geometry of the Mg ii absorption than the dependence on the inclination of the galaxies. In addition to the clear azimuthal dependence in the integrated Mg ii absorption that we reported earlier in Bordoloi et al., we show that strong equivalent width Mg ii absorbers (W_r(2796) ⩾ 0.3 Å) are also asymmetrically distributed in azimuth around their host galaxies: 72% of the absorbers in Kacprzak et al., and 100% of the close-in absorbers within 35 kpc of the center of their host galaxies, are located within 50° of the host galaxy's projected semi minor axis. It is shown that either composite models consisting of a simple bipolar component plus a spherical or disk component, or a single highly softened bipolar distribution, can well represent the azimuthal dependencies observed in both the stacked spectrum and quasar absorption-line data sets within 40 kpc. Simultaneously fitting both data sets, we find that in the composite model the bipolar cone has an opening angle of ∼100° (i.e., confined to within 50° of the disk axis) and contains about two-thirds of the total Mg ii absorption in the system. The single softened cone model has an exponential fall off with azimuthal angle with an exponential scale length in opening angle of about 45°. We conclude that the distribution of Mg ii gas at low impact parameters is not the same as that found at high impact parameters. At larger impact parameters beyond 40 kpc, there is evidence for a much more symmetric distribution, significantly different from that closer in to the galaxies.

We investigate the accretion of solid materials onto circumplanetary disks from heliocentric orbits rotating in protoplanetary disks, which is a key process for the formation of regular satellite systems. In the late stage of the gas-capturing phase of giant planet formation, the accreting gas from protoplanetary disks forms circumplanetary disks. Since the accretion flow toward the circumplanetary disks affects the particle motion through gas drag force, we use hydrodynamic simulation data for the gas drag term to calculate the motion of solid materials. We consider a wide range of size for the solid particles (10⁻²–10⁶ m), and find that the accretion efficiency of the solid particles peaks around 10 m sized particles because energy dissipation of drag with circum-planetary disk gas in this size regime is most effective. The efficiency for particles larger than 10 m becomes lower because gas drag becomes less effective. For particles smaller than 10 m, the efficiency is lower because the particles are strongly coupled with the background gas flow, which prevents particles from accretion. We also find that the distance from the planet where the particles are captured by the circumplanetary disks is in a narrow range and well described as a function of the particle size.

Only a very small fraction of the asteroid population at size scales comparable to the object that exploded over Chelyabinsk, Russia has been discovered to date, and physical properties are poorly characterized. We present previously unreported detections of 105 close approaching near-Earth objects (NEOs) by the Wide-field Infrared Survey Explorer (WISE) mission's NEOWISE project. These infrared observations constrain physical properties such as diameter and albedo for these objects, many of which are found to be smaller than 100 m. Because these objects are intrinsically faint, they were detected by WISE during very close approaches to the Earth, often at large apparent on-sky velocities. We observe a trend of increasing albedo with decreasing size, but as this sample of NEOs was discovered by visible light surveys, it is likely that selection biases against finding small, dark NEOs influence this finding.

We present an all-sky sample of 984 candidate intermediate-mass Galactic star-forming regions that are color selected from the Infrared Astronomical Satellite (IRAS) Point Source Catalog and morphologically classify each object using mid-infrared Wide-field Infrared Survey Explorer (WISE) images. Of the 984 candidates, 616 are probable star-forming regions (62.6%), 128 are filamentary structures (13.0%), 39 are point-like objects of unknown nature (4.0%), and 201 are galaxies (20.4%). We conduct a study of four of these regions, IRAS 00259+5625, IRAS 00420+5530, IRAS 01080+5717, and IRAS 05380+2020, at Galactic latitudes |b| > 5° using optical spectroscopy from the Wyoming Infrared Observatory, along with near-infrared photometry from the Two-Micron All Sky Survey, to investigate their stellar content. New optical spectra, color–magnitude diagrams, and color–color diagrams reveal their extinctions, spectrophotometric distances, and the presence of small stellar clusters containing 20–78 M_☉ of stars. These low-mass diffuse star clusters contain ∼65–250 stars for a typical initial mass function, including one or more mid-B stars as their most massive constituents. Using infrared spectral energy distributions we identify young stellar objects near each region and assign probable masses and evolutionary stages to the protostars. The total infrared luminosity lies in the range 190–960 L_☉, consistent with the sum of the luminosities of the individually identified young stellar objects.

The Schmidt–Kennicutt relation links the surface densities of gas to the star formation rate in galaxies. The physical origin of this relation, and in particular its break, i.e., the transition between an inefficient regime at low gas surface densities and a main regime at higher densities, remains debated. Here, we study the physical origin of the star formation relations and breaks in several low-redshift galaxies, from dwarf irregulars to massive spirals. We use numerical simulations representative of the Milky Way and the Large and Small Magellanic Clouds with parsec up to subparsec resolution, and which reproduce the observed star formation relations and the relative variations of the star formation thresholds. We analyze the role of interstellar turbulence, gas cooling, and geometry in drawing these relations at 100 pc scale. We suggest in particular that the existence of a break in the Schmidt–Kennicutt relation could be linked to the transition from subsonic to supersonic turbulence and is independent of self-shielding effects. With this transition being connected to the gas thermal properties and thus to the metallicity, the break is shifted toward high surface densities in metal-poor galaxies, as observed in dwarf galaxies. Our results suggest that together with the collapse of clouds under self-gravity, turbulence (injected at galactic scale) can induce the compression of gas and regulate star formation.

In contrast to massive, bulge hosting galaxies, very few supermassive black holes (SMBHs) are known in either low-mass or bulgeless galaxies. Such a population could provide clues to the origins of SMBHs and to secular pathways for their growth. Using the all-sky Wide-field Infrared Survey Explorer (WISE ) survey, and bulge-to-disk decompositions from the Sloan Digital Sky Survey (SDSS) Data Release 7, we report the discovery of a population of local (z < 0.3) bulgeless disk galaxies with extremely red mid-infrared colors which are highly suggestive of a dominant active galactic nucleus (AGN), despite having no optical AGN signatures in their SDSS spectra. Using various mid-infrared selection criteria from the literature, there are between 30 and over 300 bulgeless galaxies with possible AGNs. Other known scenarios that can heat the dust to high temperatures do not appear to explain the observed colors of this sample. If these galaxies are confirmed to host AGNs, this study will provide a breakthrough in characterizing the properties of SMBHs in the low bulge mass regime and in understanding their relation with their host galaxies. Mid-infrared selection identifies AGNs that dominate their host galaxy's emission and therefore reveal a different AGN population than that uncovered by optical studies. We find that the fraction of all galaxies identified as candidate AGNs by WISE is highest at lower stellar masses and drops dramatically in higher mass galaxies, in striking contrast to the findings from optical studies.

We measured the trigonometric parallax of an H₂O maser source associated with the massive star-forming region NGC 6334I(North), hereafter as NGC 6334I(N), with the VLBI Exploration of Radio Astrometry. The derived annual parallax is 0.789 ± 0.161 mas, corresponding to a distance of 1.26$^{+0.33}_{-0.21}$ kpc. Using the new distance, we recalculated the physical parameters (masses and luminosities) of the dust continuum cores in the region, and the revised parameters are only ∼50% of their originally reported values. We also traced 23 relative proper motions of the H₂O masers associated with SMA1 (central millimeter source in the region) between epochs, which exhibit an average amplitude of maser proper motion of ∼2.03 mas yr⁻¹ (∼12.22 km s⁻¹), tracing a bipolar outflow. The bipolar outflow structure extends through ∼600 mas (∼720 AU), with a dynamical timescale of ∼295 yr. Using an expanding flow model, we derived the most plausible dynamical center of the outflow, pointing to SMA1b (1.3 cm and 7 mm continuum source) as the outflow driving source. Based on our results and other existing parallax results, we derive the pitch angles of the Sagittarius arm as 62 ± 54 along the Galactic longitude range of ∼07 − ∼401 assuming a perfect logarithmic spiral. We derived the peculiar motion of NGC 6334I(N) to be −4 ± 1 km s⁻¹ toward the Galactic center, 8 ± 2 km s⁻¹ in the direction of the Galactic rotation, and 25 ± 2 km s⁻¹ toward the Galactic north pole.

In this paper, the nonlinear evolution of a bistable interstellar medium is investigated using two-dimensional simulations with a realistic cooling rate, thermal conduction, and physical viscosity. The calculations are performed using periodic boundary conditions without any external dynamical forcing. As the initial condition, a spatially uniform unstable gas under thermal equilibrium is considered. At the initial stage, the unstable gas quickly segregates into two phases: cold neutral medium (CNM) and warm neutral medium (WNM). Then, self-sustained turbulence with velocity dispersion of 0.1–0.2 km s⁻¹ is observed in which the CNM moves around in the WNM. We find that the interfacial medium (IFM) between the CNM and WNM plays an important role in sustaining the turbulence. The self-sustaining mechanism can be divided into two steps. First, thermal conduction drives fast flows streaming into concave CNM surfaces toward the WNM. The kinetic energy of the fast flows in the IFM is incorporated into that of the CNM through the phase transition. Second, turbulence inside the CNM deforms interfaces and forms other concave CNM surfaces, leading to fast flows in the IFM. This drives the first step again and a cycle is established by which turbulent motions are self-sustained.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was created by adding polarimetric capability to the BLAST experiment that was flown in 2003, 2005, and 2006. BLASTPol inherited BLAST's 1.8 m primary and its Herschel/SPIRE heritage focal plane that allows simultaneous observation at 250, 350, and 500 μm. We flew BLASTPol in 2010 and again in 2012. Both were long duration Antarctic flights. Here we present polarimetry of the nearby filamentary dark cloud Lupus I obtained during the 2010 flight. Despite limitations imposed by the effects of a damaged optical component, we were able to clearly detect submillimeter polarization on degree scales. We compare the resulting BLASTPol magnetic field map with a similar map made via optical polarimetry. (The optical data were published in 1998 by J. Rizzo and collaborators.) The two maps partially overlap and are reasonably consistent with one another. We compare these magnetic field maps to the orientations of filaments in Lupus I, and we find that the dominant filament in the cloud is approximately perpendicular to the large-scale field, while secondary filaments appear to run parallel to the magnetic fields in their vicinities. This is similar to what is observed in Serpens South via near-IR polarimetry, and consistent with what is seen in MHD simulations by F. Nakamura and Z. Li.

The modeling of turbulence, whether it be numerical or analytical, is a difficult challenge. Turbulence is amenable to analysis with linear theory if it is subject to rapid distortions, i.e., motions occurring on a timescale that is short compared to the timescale for nonlinear interactions. Such an approach (referred to as rapid distortion theory) could prove useful for understanding aspects of astrophysical turbulence, which is often subject to rapid distortions, such as supernova explosions or the free-fall associated with gravitational instability. As a proof of principle, a particularly simple problem is considered here: the evolution of vorticity due to a planar rarefaction in an ideal gas. Analytical solutions are obtained for incompressive modes having a wave vector perpendicular to the distortion; as in the case of gradient-driven instabilities, these are the modes that couple most strongly to the mean flow. Vorticity can either grow or decay in the wake of a rarefaction front, and there are two competing effects that determine which outcome occurs: entropy fluctuations couple to the mean pressure gradient to produce vorticity via baroclinic effects, whereas vorticity is damped due to the conservation of angular momentum as the fluid expands. Whether vorticity grows or decays depends upon the ratio of entropic to vortical fluctuations at the location of the front; growth occurs if this ratio is of order unity or larger. In the limit of purely entropic fluctuations in the ambient fluid, a strong rarefaction generates vorticity with a turbulent Mach number on the order of the rms of the ambient entropy fluctuations. The analytical results are shown to compare well with results from two- and three-dimensional numerical simulations. Analytical solutions are also derived in the linear regime of Reynolds-averaged turbulence models. This highlights an inconsistency in standard turbulence models that prevents them from accurately capturing the physics of rarefaction–turbulence interaction. In addition to providing physical insight, the solutions derived here can be used to verify algorithms of both the Reynolds-averaged and direct numerical simulation variety. Finally, dimensional analysis of the equations indicates that rapid distortion of turbulence can give rise to two distinct regimes in the turbulent spectrum: a distortion range at large scales where linear distortion effects dominate, and an inertial range at small scales where nonlinear effects dominate.

Observations of occultations of bright γ-ray sources by the Sun may reveal predicted pair halos around blazars and/or new physics, such as, e.g., hypothetical light dark matter particles—axions. We use Fermi Gamma-Ray Space Telescope (Fermi) data to analyze four occultations of blazar 3C 279 by the Sun on October 8 each year from 2008 to 2011. A combined analysis of the observations of these occultations allows a point-like source at the position of 3C 279 to be detected with significance of ≈3σ, but does not reveal any significant excess over the flux expected from the quiescent Sun. The likelihood ratio test rules out complete transparency of the Sun to the blazar γ-ray emission at a 3σ confidence level.

Within the next five years, it is expected that the Advanced LIGO/Virgo network will have reached a sensitivity sufficient to enable the routine detection of gravitational waves. Beyond the initial detection, the scientific promise of these instruments relies on the effectiveness of our physical parameter estimation capabilities. A major part of this effort has been toward the detection and characterization of gravitational waves from compact binary coalescence, e.g., the coalescence of binary neutron stars. While several previous studies have investigated the accuracy of parameter estimation with advanced detectors, the majority have relied on approximation techniques such as the Fisher Matrix which are insensitive to the non-Gaussian nature of the gravitational wave posterior distribution function. Here we report average statistical uncertainties that will be achievable for strong detection candidates (S/N = 20) over a comprehensive sample of source parameters. We use the Markov Chain Monte Carlo based parameter estimation software developed by the LIGO/Virgo Collaboration with the goal of updating the previously quoted Fisher Matrix bounds. We find the recovery of the individual masses to be fractionally within 9% (15%) at the 68% (95%) credible intervals for equal-mass systems, and within 1.9% (3.7%) for unequal-mass systems. We also find that the Advanced LIGO/Virgo network will constrain the locations of binary neutron star mergers to a median uncertainty of 5.1 deg² (13.5 deg²) on the sky. This region is improved to 2.3 deg² (6 deg²) with the addition of the proposed LIGO India detector to the network. We also report the average uncertainties on the luminosity distances and orbital inclinations of strong detections that can be achieved by different network configurations.

Although the mechanisms responsible for heating the Sun's corona and accelerating the solar wind are still being actively investigated, it is largely accepted that photospheric motions provide the energy source and that the magnetic field must play a key role in the process. Verdini et al. presented a model for heating and accelerating the solar wind based on the turbulent dissipation of Alfvén waves. We first use a time-dependent model of the solar wind to reproduce one of Verdini et al.'s solutions; then, we extend its application to the case where the energy equation includes thermal conduction and radiation losses, and the upper chromosphere is part of the computational domain. Using this model, we explore the parameter space and describe the characteristics of a fast solar wind solution. We discuss how this formulation may be applied to a three-dimensional MHD model of the corona and solar wind.

We report results of three-dimensional magnetohydrodynamical (MHD) simulations of global accretion disks threaded with weak vertical magnetic fields. We perform the simulations in the spherical coordinates with different temperature profiles and accordingly different rotation profiles. In the cases with a spatially constant temperature, because the rotation frequency is vertically constant in the equilibrium condition, general properties of the turbulence excited by magnetorotational instability are quantitatively similar to those obtained in local shearing box simulations. On the other hand, in the cases with a radially variable temperature profile, the vertical differential rotation, which is inevitable in the equilibrium condition, winds up the magnetic field lines in addition to the usual radial differential rotation. As a result, the coherent wound magnetic fields contribute to the Maxwell stress in the surface regions. We obtain nondimensional density and velocity fluctuations ∼0.1–0.2 at the midplane. The azimuthal power spectra of the magnetic fields show shallower slopes, ∼m⁰ − m⁻¹, than those of velocity and density. The Poynting flux associated with the MHD turbulence drives intermittent and structured disk winds as well as sound-like waves toward the midplane. The mass accretion mainly occurs near the surfaces, and the gas near the midplane slowly moves outward in the time domain of the present simulations. The vertical magnetic fields are also dragged inward in the surface regions, while they stochastically move outward and inward around the midplane. We also discuss an observational implication of induced spiral structure in the simulated turbulent disks.

Understanding how black holes accrete and supply feedback to their environment is one of the outstanding challenges of modern astrophysics. Swift J1910.2–0546 is a candidate black hole low-mass X-ray binary that was discovered in 2012 when it entered an accretion outburst. To investigate the binary configuration and the accretion morphology, we monitored the evolution of the outburst for ≃3 months at X-ray, UV, optical (B, V, R, I), and near-infrared (J, H, K) wavelengths using Swift and SMARTS. The source evolved from a hard to a soft X-ray spectral state with a relatively cold accretion disk that peaked at ≃0.5 keV. A Chandra/HETG spectrum obtained during this soft state did not reveal signatures of an ionized disk wind. Both the low disk temperature and the absence of a detectable wind could indicate that the system is viewed at relatively low inclination. The multi-wavelength light curves revealed two notable features that appear to be related to X-ray state changes. First, a prominent flux decrease was observed in all wavebands ≃ 1–2 weeks before the source entered the soft state. This dip occurred in (0.6–10 keV) X-rays ≃ 6 days later than at longer wavelengths, which could possibly reflect the viscous timescale of the disk. Second, about two weeks after the source transitioned back into the hard state, the UV emission significantly increased while the X-rays steadily decayed. We discuss how these observations may reflect changes in the accretion morphology, perhaps related to the quenching/launch of a jet or the collapse/recovery of a hot flow.

We perform a systematic analysis of neutron star radius constraints from five quiescent low-mass X-ray binaries and examine how they depend on measurements of their distances and amounts of intervening absorbing material, as well as their assumed atmospheric compositions. We construct and calibrate to published results a semi-analytic model of the neutron star atmosphere which approximates these effects for the predicted masses and radii. Starting from mass and radius probability distributions established from hydrogen-atmosphere spectral fits of quiescent sources, we apply this model to compute alternate sets of probability distributions. We perform Bayesian analyses to estimate neutron star mass–radius curves and equation of state (EOS) parameters that best-fit each set of distributions, assuming the existence of a known low-density neutron star crustal EOS, a simple model for the high-density EOS, causality, and the observation that the neutron star maximum mass exceeds 2 M_☉. We compute the posterior probabilities for each set of distance measurements and assumptions about absorption and composition. We find that, within the context of our assumptions and our parameterized EOS models, some absorption models are disfavored. We find that neutron stars composed of hadrons are favored relative to those with exotic matter with strong phase transitions. In addition, models in which all five stars have hydrogen atmospheres are found to be weakly disfavored. Our most likely models predict neutron star radii that are consistent with current experimental results concerning the nature of the nucleon–nucleon interaction near the nuclear saturation density.

We present observations of the eccentric γ-ray binary B1259–63/LS 2883 with the ChandraX-ray Observatory. The images reveal a variable, extended (about 4'', or ∼1000 times the binary orbit size) structure, which appears to be moving away from the binary with the velocity of 0.05 of the speed of light. The observed emission is interpreted as synchrotron radiation from relativistic particles supplied by the pulsar. However, the fast motion through the circumbinary medium would require the emitting cloud to be loaded with a large amount of baryonic matter. Alternatively, the extended emission can be interpreted as a variable extrabinary shock in the pulsar wind outflow launched near binary apastron. The resolved variable X-ray nebula provides an opportunity to probe pulsar winds and their interaction with stellar winds in a previously inaccessible way.

X-ray emission from the supermassive binary system η Car declines sharply around periastron. This X-ray minimum has two distinct phases—the lowest flux phase in the first ∼3 weeks and a brighter phase thereafter. In 2009, the ChandraX-ray Observatory monitored the first phase five times and found the lowest observed flux at ∼1.9 × 10⁻¹² erg cm⁻² s⁻¹ (3–8 keV). The spectral shape changed such that the hard band above ∼4 keV dropped quickly at the beginning and the soft band flux gradually decreased to its lowest observed value in ∼2 weeks. The hard band spectrum had begun to recover by that time. This spectral variation suggests that the shocked gas producing the hottest X-ray gas near the apex of the wind–wind collision (WWC) is blocked behind the dense inner wind of the primary star, which later occults slightly cooler gas downstream. Shocked gas previously produced by the system at earlier orbital phases is suggested to produce the faint residual X-ray emission seen when the emission near the apex is completely blocked by the primary wind. The brighter phase is probably caused by the re-appearance of the WWC plasma, whose emissivity significantly declined during the occultation. We interpret this to mean that the X-ray minimum is produced by a hybrid mechanism of an occultation and a decline in the emissivity of the WWC shock. We constrain timings of superior conjunction and periastron based on these results.

We have compiled photometry at 3.4, 4.6, 12, and 22 μm from the all-sky survey performed by the Wide-field Infrared Survey Explorer (WISE) for all known members of the Taurus complex of dark clouds. Using these data and photometry from the Spitzer Space Telescope, we have identified members with infrared excess emission from circumstellar disks and have estimated the evolutionary stages of the detected disks, which include 31 new full disks and 16 new candidate transitional, evolved, evolved transitional, and debris disks. We have also used the WISE All-Sky Source Catalog to search for new disk-bearing members of Taurus based on their red infrared colors. Through optical and near-infrared spectroscopy, we have confirmed 26 new members with spectral types of M1–M7. The census of disk-bearing stars in Taurus should now be largely complete for spectral types earlier than ∼M8 (M ≳ 0.03 M_☉).

The far-ultraviolet (FUV; 912–1700 Å) radiation field from accreting central stars in classical T Tauri systems influences the disk chemistry during the period of giant planet formation. The FUV field may also play a critical role in determining the evolution of the inner disk (r < 10 AU), from a gas- and dust-rich primordial disk to a transitional system where the optically thick warm dust distribution has been depleted. Previous efforts to measure the true stellar+accretion-generated FUV luminosity (both hot gas emission lines and continua) have been complicated by a combination of low-sensitivity and/or low-spectral resolution and did not include the contribution from the bright Lyα emission line. In this work, we present a high-resolution spectroscopic study of the FUV radiation fields of 16 T Tauri stars whose dust disks display a range of evolutionary states. We include reconstructed Lyα line profiles and remove atomic and molecular disk emission (from H₂ and CO fluorescence) to provide robust measurements of both the FUV continuum and hot gas lines (e.g., Lyα, N v, C iv, He ii) for an appreciable sample of T Tauri stars for the first time. We find that the flux of the typical classical T Tauri star FUV radiation field at 1 AU from the central star is ∼10⁷ times the average interstellar radiation field. The Lyα emission line contributes an average of 88% of the total FUV flux, with the FUV continuum accounting for an average of 8%. Both the FUV continuum and Lyα flux are strongly correlated with C iv flux, suggesting that accretion processes dominate the production of both of these components. On average, only ∼0.5% of the total FUV flux is emitted between the Lyman limit (912 Å) and the H₂ (0–0) absorption band at 1110 Å. The total and component-level high-resolution radiation fields are made publicly available in machine-readable format.

We present measurements of the luminosity and color-dependence of galaxy clustering at 0.2 < z < 1.0 in the Prism Multi-object Survey. We quantify the clustering with the redshift-space and projected two-point correlation functions, ξ(r_p, π) and w_p(r_p), using volume-limited samples constructed from a parent sample of over ∼130, 000 galaxies with robust redshifts in seven independent fields covering 9 deg² of sky. We quantify how the scale-dependent clustering amplitude increases with increasing luminosity and redder color, with relatively small errors over large volumes. We find that red galaxies have stronger small-scale (0.1 Mpc h⁻¹ < r_p < 1 Mpc h⁻¹) clustering and steeper correlation functions compared to blue galaxies, as well as a strong color dependent clustering within the red sequence alone. We interpret our measured clustering trends in terms of galaxy bias and obtain values of b_gal ≈ 0.9–2.5, quantifying how galaxies are biased tracers of dark matter depending on their luminosity and color. We also interpret the color dependence with mock catalogs, and find that the clustering of blue galaxies is nearly constant with color, while redder galaxies have stronger clustering in the one-halo term due to a higher satellite galaxy fraction. In addition, we measure the evolution of the clustering strength and bias, and we do not detect statistically significant departures from passive evolution. We argue that the luminosity– and color–environment (or halo mass) relations of galaxies have not significantly evolved since z ∼ 1. Finally, using jackknife subsampling methods, we find that sampling fluctuations are important and that the COSMOS field is generally an outlier, due to having more overdense structures than other fields; we find that "cosmic variance" can be a significant source of uncertainty for high-redshift clustering measurements.

We present new Very Large Array 6 cm H₂CO observations toward four extragalactic radio continuum sources (B0212+735, 3C 111, NRAO 150, and BL Lac) to explore the structure of foreground Galactic clouds as revealed by absorption variability. This project adds a new epoch in the monitoring observations of the sources reported by Marscher and collaborators in the mid-1990s. Our new observations confirm the monotonic increase in H₂CO absorption strength toward NRAO 150. We do not detect significant variability of our 2009 spectra with respect to the 1994 spectra of 3C111, B0212+735, and BL Lac; however, we find significant variability of the 3C111 2009 spectrum with respect to archive observations conducted in 1991 and 1992. Our analysis supports that changes in absorption lines could be caused by chemical and/or geometrical gradients in the foreground clouds and not necessarily by small-scale (∼10 AU) high-density molecular clumps within the clouds.

Infrared emission is an invaluable tool for quantifying star formation in galaxies. Because the 8 μm polycyclic aromatic hydrocarbon (PAH) emission has been found to correlate with other well-known star formation tracers, it has widely been used as a star formation rate (SFR) tracer. There are, however, studies that challenge the accuracy and reliability of the 8 μm emission as a SFR tracer. Our study, part of the Herschel (Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA) M33 Extended Survey (HERM33ES) open time key program, aims at addressing this issue by analyzing the infrared emission from the nearby spiral galaxy M33 at the high spatial scale of ∼75 pc. Combining data from the Herschel Space Observatory and the Spitzer Space Telescope, we find that the 8 μm emission is better correlated with the 250 μm emission, which traces cold interstellar gas, than with the 24 μm emission. Furthermore, the L(8)/L(250) ratio is more tightly correlated with the 3.6 μm emission, a tracer of evolved stellar populations and stellar mass, than with a combination of Hα and 24 μm emission, a tracer of SFR. The L(8)/L(24) ratio is highly depressed in 24 μm luminous regions, which correlate with known H ii regions. We also compare our results with the dust emission models by Draine & Li. We confirm that the depression of 8 μm PAH emission near star-forming regions is higher than what is predicted by models; this is possibly an effect of increased stellar radiation from young stars destroying the dust grains responsible for the 8 μm emission as already suggested by other authors. We find that the majority of the 8 μm emission is fully consistent with heating by the diffuse interstellar medium, similar to what recently determined for the dust emission in M31 by Draine et al. We also find that the fraction of 8 μm emission associated with the diffuse interstellar radiation field ranges between ∼60% and 80% and is 40% larger than the diffuse fraction at 24 μm.

How and where cosmic rays are produced, and how they diffuse through various turbulent media, represent fundamental problems in astrophysics with far-reaching implications, both in terms of our theoretical understanding of high-energy processes in the Milky Way and beyond, and the successful interpretation of space-based and ground based GeV and TeV observations. For example, recent and ongoing detections, e.g., by Fermi (in space) and HESS (in Namibia), of γ-rays produced in regions of dense molecular gas hold important clues for both processes. In this paper, we carry out a comprehensive numerical investigation of relativistic particle acceleration and transport through turbulent magnetized environments in order to derive broadly useful scaling laws for the energy diffusion coefficients.

The molecular structure of methanediol has been investigated by means of quantum chemical calculations. Two conformers, corresponding to C₂ and C_s symmetries, respectively, were considered. The C₂ conformer is found to lie about 1.7 (at 298 K) or 2.3 (at 0 K) kcal mol⁻¹ below the C_s conformer. Predictions for their rotational constants, vibrational frequencies, IR intensities, and dipole moments have been provided. The lowest-lying isomer has a very low dipole moment, around 0.03 D, whereas the C_s conformer has a relatively high dipole moment, namely, 2.7 D. The barrier for the C_s→C₂ process is predicted to be around 0.7–1 kcal mol⁻¹. Based on the energetic results the proportion of the C_s conformer is likely to be negligible under low temperature conditions, such as in the interstellar medium. Therefore, it is predicted that detection by radioastronomy of methanediol would be rather unlikely.

The Hubble Space Telescope/NICMOS transmission spectrum measurements of HD 189733b that suggest the detection of methane (CH₄) in an exoplanet atmosphere have been a source of significant controversy. With what is probably the best analyzed exoplanet spectroscopy data set to date, different teams, using different methods, have claimed evidence both contradicting and supporting the original findings. Here, we report results from a uniform spectral retrieval analysis of the three independent published spectra together with null hypothesis testing. Based on Bayesian model comparison, we find that two of the three spectra show strong evidence (⩾3.6σ) for the detection of molecular features mainly due to water and methane while the third is consistent with a weak molecular detection at the 2.2σ level. We interpret the agreement in the spectral modulation established by previous authors and the atmospheric retrieval results presented here, as a confirmation of the original detection of molecular absorbers in the atmosphere of HD 189733b.

Observing the Sun at high time and spatial scales is a step toward understanding the finest and fundamental scales of heating events in the solar corona. The high-resolution coronal (Hi-C) instrument has provided the highest spatial and temporal resolution images of the solar corona in the EUV wavelength range to date. Hi-C observed an active region on 2012 July 11 that exhibits several interesting features in the EUV line at 193 Å. One of them is the existence of short, small brightenings "sparkling" at the edge of the active region; we call these EUV bright dots (EBDs). Individual EBDs have a characteristic duration of 25 s with a characteristic length of 680 km. These brightenings are not fully resolved by the SDO/AIA instrument at the same wavelength; however, they can be identified with respect to the Hi-C location of the EBDs. In addition, EBDs are seen in other chromospheric/coronal channels of SDO/AIA, which suggests a temperature between 0.5 and 1.5 MK. Based on their frequency in the Hi-C time series, we define four different categories of EBDs: single peak, double peak, long duration, and bursty. Based on a potential field extrapolation from an SDO/HMI magnetogram, the EBDs appear at the footpoints of large-scale, trans-equatorial coronal loops. The Hi-C observations provide the first evidence of small-scale EUV heating events at the base of these coronal loops, which have a free magnetic energy of the order of 10²⁶ erg.

A study of the kinematics and arrival times of coronal mass ejections (CMEs) at Earth, derived from time-elongation maps (J-maps) constructed from STEREO/heliospheric imager (HI) observations, provides an opportunity to understand the heliospheric evolution of CMEs in general. We implement various reconstruction techniques, based on the use of time-elongation profiles of propagating CMEs viewed from single or multiple vantage points, to estimate the dynamics of three geo-effective CMEs. We use the kinematic properties, derived from analysis of the elongation profiles, as inputs to the Drag Based Model for the distance beyond which the CMEs cannot be tracked unambiguously in the J-maps. The ambient solar wind into which these CMEs, which travel with different speeds, are launched, is different. Therefore, these CMEs will evolve differently throughout their journey from the Sun to 1 AU. We associate the CMEs, identified and tracked in the J-maps, with signatures observed in situ near 1 AU by the WIND spacecraft. By deriving the kinematic properties of each CME, using a variety of existing methods, we assess the relative performance of each method for the purpose of space weather forecasting. We discuss the limitations of each method, and identify the major constraints in predicting the arrival time of CMEs near 1 AU using HI observations.

Data from the Princess Sirindhorn Neutron Monitor at Doi Inthanon, Thailand, with a vertical cutoff rigidity of 16.8 GV, were utilized to determine the diurnal anisotropy (DA) of Galactic cosmic rays (GCRs) near Earth during solar minimum conditions between 2007 November and 2010 November. We identified trains of enhanced DA over several days, which often recur after a solar rotation period (∼27 days). By investigating solar coronal holes as identified from synoptic maps and solar wind parameters, we found that the intensity and anisotropy of cosmic rays are associated with the high-speed streams (HSSs) in the solar wind, which are in turn related to the structure and evolution of coronal holes. An enhanced DA was observed after the onset of some, but not all, HSSs. During time periods of recurrent trains, the DA was often enhanced or suppressed according to the sign of the interplanetary magnetic field B, which suggests a contribution from a mechanism involving a southward gradient in the GCR density, n, and a gradient anisotropy along B × ∇n. In one non-recurrent and one recurrent sequence, an HSS from an equatorial coronal hole was merged with that from a trailing mid-latitude extension of a polar coronal hole, and the slanted HSS structure in space with suppressed GCR density can account for the southward GCR gradient. We conclude that the gradient anisotropy is a source of temporary changes in the GCR DA under solar minimum conditions, and that the latitudinal GCR gradient can sometimes be explained by the coronal hole morphology.

Several lines of argument support the existence of a link between activity at the nuclei of galaxies, in the form of an accreting supermassive black hole, and star formation activity in these galaxies. Radio jets have long been argued to be an ideal mechanism that allows active galactic nuclei (AGNs) to interact with their host galaxies and affect star formation. We use a sample of radio sources in the North Ecliptic Pole (NEP) field to study the nature of this putative link, by means of spectral energy distribution (SED) fitting. We employ the excellent spectral coverage of the AKARI infrared space telescope and the rich ancillary data available in the NEP to build SEDs extending from UV to far-IR wavelengths. We find a significant AGN component in our sample of relatively faint radio sources (<mJy). A positive correlation is found between the luminosity of the AGN component and that of star formation in the host galaxy, independent of the radio luminosity. In contrast, for narrow redshift and AGN luminosity ranges, we find that increasing radio luminosity leads to a decrease in the specific star formation rate. The most radio-loud AGNs are found to lie on the main sequence of star formation for their respective redshifts. For the first time, we potentially see such a two-sided feedback process in the same sample. We discuss the possible suppression of star formation, but not total quenching, in systems with strong radio jets, that supports the maintenance nature of feedback from radio AGN jets.

We extend our previous model-independent determination of the intergalactic background light, based purely on galaxy survey data, out to a wavelength of 5 μm. Our approach enables us to constrain the range of photon densities, based on the uncertainties from observationally determined luminosity densities and colors. We further determine a 68% confidence upper and lower limit on the opacity of the universe to γ-rays up to energies of 1.6/(1 + z) TeV. A comparison of our lower limit redshift-dependent opacity curves to the opacity limits derived from the results of both ground-based air Cerenkov telescope and Fermi-LAT observations of PKS 1424+240 allows us to place a new upper limit on the redshift of this source, independent of IBL modeling.

We used the continuous-time random-walk Monte Carlo technique to study anew the formation of H₂ on the surfaces of interstellar dust grains in diffuse interstellar clouds. For our study, we considered three different grain materials, olivine (a polycrystalline silicate), amorphous silicate, and amorphous carbon, as well as a grain temperature that depends on granular size. For some runs, we included temperature fluctuations. Four different granular surfaces were used, one "flat" with one type of binding site due to physisorption, one "rough" with five different types of physisorption binding sites due to lateral forces, and two with sites for chemisorption, one in which chemisorption sites are entered through precursor physisorption sites, and one in which chemisorption is direct but occurs with a barrier for the adsorption of the first hydrogen atom. We found that on flat and rough olivine surfaces, molecular hydrogen is formed at low efficiencies, with smaller grains contributing very little despite their large numbers due to high temperatures. For flat amorphous carbon and amorphous silicate surfaces, the efficiency increases, reaching unity for the largest grains. For models with barrierless chemisorption, the efficiency of formation of H₂ is near unity at all grain sizes considered, while for direct chemisorption via a barrier, we found efficiencies of 0.13–0.6 depending upon the barrier, but independent of grain size. Treating the flat olivine and amorphous silicate surfaces with temperature fluctuations increases the efficiency of H₂ formation.

The evolutionary connection between nuclear starbursts and active galactic nuclei (AGNs) in luminous infrared galaxies (LIRGs) and ultraluminous infrared galaxies (ULIRGs), which result from galaxy interactions and mergers and produce the bulk of their radiation as infrared (IR) emission, is not well understood. To this effort, we present and examine new spectrophotometric data for five U/LIRGs (10¹¹ < L_IR < 10¹³L_☉) within the IRAS 2 Jy Redshift Survey with 0.05 ≲ z ≲ 0.07. We show that our sample consists almost entirely of composite objects—thus hosting both a nuclear starburst and an AGN—using the BPT diagrams. We then show that for our sample of U/LIRGs the properties that describe their nuclear starbursts and AGNs (e.g., star formation rate, L[O iii], optical D parameter, D4000, and EW(Hδ)) are independent of one another, ensuring that no biases affect correlations between these parameters and the object locations on the BPT diagrams. Finally, we derive evolutionary paths on the BPT diagram involving ${[{\rm N\,{\scriptsize{II}}}]/{\rm H}\alpha }$ that are based on how these parameters vary between two U/LIRGs positioned at the end-points of these paths. The U/LIRGs at the end-points of a given path represent the beginning and end states of a U/LIRG evolving along that path. These paths may be able to specifically explain how all local U/LIRGs evolve along the BPT diagram, and serve as a starting point for future quantitative analysis on the evolution of U/LIRGs.

The gamma-ray-detected blazar 3C 454.3 exhibits dramatic flux and polarization variations in the optical and near-infrared bands. In 2010 December, the object emitted a very bright outburst. We monitored it for approximately four years (including the 2010 outburst) by optical and near-infrared photopolarimetry. During the 2010 outburst, the object emitted two rapid, redder brightenings, at which the polarization degrees (PDs) in both bands increased significantly and the bands exhibited a frequency-dependent polarization. The observed frequency-dependent polarization leads us to propose that the polarization vector is composed of two vectors. Therefore, we separate the observed polarization vectors into short- and long-term components that we attribute to the emissions of the rapid brightenings and the outburst that varied the timescale of days and months, respectively. The estimated PD of the short-term component is greater than the maximum observed PD and is close to the theoretical maximum PD. We constrain the bulk Lorentz factors and inclination angles between the jet axis and the line of sight from the estimated PDs. In this case, the inclination angle of the emitting region of short-term component from the first rapid brightening should be equal to 90°, because the estimated PD of the short-term component was approximately equal to the theoretical maximum PD. Thus, the Doppler factor at the emitting region of the first rapid brightening should be equal to the bulk Lorentz factor.

Our current view of galaxies considers them as systems of stars and gas embedded in extended halos of dark matter, much of it formed by the infall of smaller systems at earlier times. The true extent of a galaxy remains poorly determined, with the "virial radius" (R_vir) providing a characteristic separation between collapsed structures in dynamical equilibrium and external infalling matter. Other physical estimates of the extent of gravitational influence include the gravitational radius, gas accretion radius, and "galactopause" arising from outflows that stall at 100–200 kpc over a range of outflow parameters and confining gas pressures. Physical criteria are proposed to define bound structures, including a more realistic definition of R_vir(M_*, M_h, z_a) for stellar mass M_* and halo mass M_h, half of which formed at "assembly redshifts" ranging from z_a ≈ 0.7–1.3. We estimate the extent of bound gas and dark matter around L* galaxies to be ∼200 kpc. The new virial radii, with mean 〈R_vir〉 ≈ 200 kpc, are 40%–50% smaller than values estimated in recent Hubble Space Telescope/Cosmic Origins Spectrograph detections of H i and O vi absorbers around galaxies. In the new formalism, the Milky Way stellar mass, log M_* = 10.7 ± 0.1, would correspond to $R_{\rm vir} = 153^{+25}_{-16}$ kpc for half-mass halo assembly at z_a = 1.06 ± 0.03. The frequency per unit redshift of low-redshift O vi absorption lines in QSO spectra suggests absorber sizes ∼150 kpc when related to intervening 0.1L* galaxies. This formalism is intended to clarify semantic differences arising from observations of extended gas in galactic halos, circumgalactic medium (CGM), and filaments of the intergalactic medium (IGM). Astronomers should refer to bound gas in the galactic halo or CGM, and unbound gas at the CGM–IGM interface, on its way into the IGM.

We present a measurement of the angular bispectrum of the millimeter-wave sky in observing bands centered at roughly 95, 150, and 220 GHz, on angular scales of 1' ≲ θ ≲ 10' (multipole number 1000 ≲ l ≲ 10,000). At these frequencies and angular scales, the main contributions to the bispectrum are expected to be the thermal Sunyaev–Zel'dovich (tSZ) effect and emission from extragalactic sources, predominantly dusty, star-forming galaxies (DSFGs) and active galactic nuclei. We measure the bispectrum in 800 deg² of three-band South Pole Telescope data, and we use a multi-frequency fitting procedure to separate the bispectrum of the tSZ effect from the extragalactic source contribution. We simultaneously detect the bispectrum of the tSZ effect at >10σ, the unclustered component of the extragalactic source bispectrum at >5σ in each frequency band, and the bispectrum due to the clustering of DSFGs—i.e., the clustered cosmic infrared background (CIB) bispectrum—at >5σ. This is the first reported detection of the clustered CIB bispectrum. We use the measured tSZ bispectrum amplitude, compared to model predictions, to constrain the normalization of the matter power spectrum to be σ₈ = 0.787 ± 0.031 and to predict the amplitude of the tSZ power spectrum at l = 3000. This prediction improves our ability to separate the thermal and kinematic contributions to the total SZ power spectrum. The addition of bispectrum data improves our constraint on the tSZ power spectrum amplitude by a factor of two compared to power spectrum measurements alone and demonstrates a preference for a nonzero kinematic SZ (kSZ) power spectrum, with a derived constraint on the kSZ amplitude at l = 3000 of A_kSZ = 2.9 ± 1.6 μK², or A_kSZ = 2.6 ± 1.8 μK² if the default A_kSZ > 0 prior is removed.

We present SDO/AIA observations of an eruptive X-class flare of 2012 July 12, and compare its evolution with the predictions of a three-dimensional (3D) numerical simulation. We focus on the dynamics of flare loops that are seen to undergo slipping reconnection during the flare. In the Atmospheric Imaging Assembly (AIA) 131 Å observations, lower parts of 10 MK flare loops exhibit an apparent motion with velocities of several tens of km s⁻¹ along the developing flare ribbons. In the early stages of the flare, flare ribbons consist of compact, localized bright transition-region emission from the footpoints of the flare loops. A differential emission measure analysis shows that the flare loops have temperatures up to the formation of Fe xxiv. A series of very long, S-shaped loops erupt, leading to a coronal mass ejection observed by STEREO. The observed dynamics are compared with the evolution of magnetic structures in the "standard solar flare model in 3D." This model matches the observations well, reproducing the apparently slipping flare loops, S-shaped erupting loops, and the evolution of flare ribbons. All of these processes are explained via 3D reconnection mechanisms resulting from the expansion of a torus-unstable flux rope. The AIA observations and the numerical model are complemented by radio observations showing a noise storm in the metric range. Dm-drifting pulsation structures occurring during the eruption indicate plasmoid ejection and enhancement of the reconnection rate. The bursty nature of radio emission shows that the slipping reconnection is still intermittent, although it is observed to persist for more than an hour.

Large-scale plasma flows in the Sun's convection zone likely play a major role in solar dynamics on decadal timescales. In particular, quantifying meridional motions is a critical ingredient for understanding the solar cycle and the transport of magnetic flux. Because the signal of such features can be quite small in deep solar layers and be buried in systematics or noise, the true meridional velocity profile has remained elusive. We perform time–distance helioseismology measurements on several years worth of Global Oscillation Network Group Doppler data. A spherical harmonic decomposition technique is applied to a subset of acoustic modes to measure travel-time differences to try to obtain signatures of meridional flows throughout the solar convection zone. Center-to-limb systematics are taken into account in an intuitive yet ad hoc manner. Travel-time differences near the surface that are consistent with a poleward flow in each hemisphere and are similar to previous work are measured. Additionally, measurements in deep layers near the base of the convection zone suggest a possible equatorward flow, as well as partial evidence of a sign change in the travel-time differences at mid-convection zone depths. This analysis on an independent data set using different measurement techniques strengthens recent conclusions that the convection zone may have multiple "cells" of meridional flow. The results may challenge the common understanding of one large conveyor belt operating in the solar convection zone. Further work with helioseismic inversions and a careful study of systematic effects are needed before firm conclusions of these large-scale flow structures can be made.

Voyager 1 (V1) has been observing interstellar magnetic fields for more than one year beginning ≈2012/209, when V1 crossed a current sheet, a "CS0" having the structure of a tangential discontinuity. The inclination of this current sheet is consistent with an interstellar magnetic field B draped on a blunt heliopause. Two other current sheets (sector boundaries) were observed at ≈2012/167 and ≈2011/276 with high inclinations (99° ± 10° and 89° ± 10°, respectively). From 2013.0 to ≈2013.6, the difference between the azimuthal angle λ of B from the Parker spiral angle at the latitude 346 of V1 was λ − λ_P = 22° ± 3° and the corresponding difference of the elevation angle δ was δ − δ_P = 23° ± 8°. During 2012, the deviation from the Parker spiral angle was somewhat smaller. The interstellar magnetic field has a "west to east polarity," opposite to the direction of planetary motions. The magnitude of B varied smoothly in the range 0.38–0.59 nT with an average B = 0.486 ± 0.045 after 2012/237.7. The transition from heliosheath to interstellar magnetic fields is related to a "two-step" increase in the cosmic ray intensity observed by V1 from ≈2012.30 to ≈2012.65. The first step increase began near the end of an unusual "away-polarity" sector, and it reached a plateau when V1 moved into a "toward-polarity" sector that ended at CS0. The second step increase began slowly after V1 crossed CS0, and it ended abruptly at 2012/237.728.

Massive stars can alter physical conditions and properties of their ambient interstellar dust grains via radiative heating and shocks. The H ii regions in the Large Magellanic Cloud (LMC) offer ideal sites to study the stellar energy feedback effects on dust because stars can be resolved, and the galaxy's nearly face-on orientation allows us to unambiguously associate H ii regions with their ionizing massive stars. The Spitzer Space Telescope survey of the LMC provides multi-wavelength (3.6–160 μm) photometric data of all H ii regions. To investigate the evolution of dust properties around massive stars, we have analyzed spatially resolved IR dust emission from two classical H ii regions (N63 and N180) and two simple superbubbles (N70 and N144) in the LMC. We produce photometric spectral energy distributions (SEDs) of numerous small subregions for each region based on its stellar distributions and nebular morphologies. We use DustEM dust emission model fits to characterize the dust properties. Color–color diagrams and model fits are compared with the radiation field (estimated from photometric and spectroscopic surveys). Strong radial variations of SEDs can be seen throughout the regions, reflecting the available radiative heating. Emission from very small grains drastically increases at locations where the radiation field is the highest, while polycyclic aromatic hydrocarbons (PAHs) appear to be destroyed. PAH emission is the strongest in the presence of molecular clouds, provided that the radiation field is low.

We present far-infrared and submillimeter images of the η Crv debris disk system obtained with Herschel and SCUBA-2, as well as Hubble Space Telescope visible and near-infrared coronagraphic images. In the 70 μm Herschel image, we clearly separate the thermal emission from the warm and cold belts in the system, find no evidence for a putative dust population located between them, and precisely determine the geometry of the outer belt. We also find marginal evidence for azimuthal asymmetries and a global offset of the outer debris ring relative to the central star. Finally, we place stringent upper limits on the scattered light surface brightness of the outer ring. Using radiative transfer modeling, we find that it is impossible to account for all observed properties of the system under the assumption that both rings contain dust populations with the same properties. While the outer belt is in reasonable agreement with the expectations of steady-state collisional cascade models, albeit with a minimum grain size that is four times larger than the blow-out size, the inner belt appears to contain copious amounts of small dust grains, possibly below the blow-out size. This suggests that the inner belt cannot result from a simple transport of grains from the outer belt and rather supports a more violent phenomenon as its origin. We also find that the emission from the inner belt has not declined over three decades, a much longer timescale than its dynamical timescale, which indicates that the belt is efficiently replenished.

We conducted multi-element isotopic analyses of 11 presolar silicate grains from the Acfer 094 meteorite having unusual O isotopic compositions. Eight grains are ¹⁸O-rich, one is ¹⁶O-rich, and two are extremely ¹⁷O-rich. We constrained the grains' stellar sources by measuring their Si and Mg isotopic ratios, and also the ⁵⁴Fe/⁵⁶Fe and ⁵⁷Fe/⁵⁶Fe ratios for five grains. The Mg and Fe isotopic measurements were conducted after surrounding matrix grains were removed for more accurate ratios. Most of the ¹⁸O-rich silicates had anomalous Mg isotopic ratios, and their combined isotopic constraints are consistent with origins in low-mass Type II supernovae (SNe II) rather than high-metallicity stars. The isotopic ratios of the ¹⁶O-rich silicate are also consistent with an SN origin. Mixing small amounts of interior stellar material with the stellar envelope replicated all measured isotopic ratios except for ²⁹Si/²⁸Si and ⁵⁴Fe/⁵⁶Fe in some grains. The ²⁹Si/²⁸Si ratios of all SN-derived grains are matched by doubling the ²⁹Si yield in the Ne- and Si-burning zones. The ⁵⁴Fe/⁵⁶Fe ratios of the grains imply elemental fractionation in the Si/S zone, or introduction of isotopically solar Fe by secondary processing. The two highly ¹⁷O-rich silicates exhibited significant ²⁵Mg and/or ²⁶Mg enrichments and their isotopic ratios are best explained by strong dilution of 1.15 M_☉ CO nova matter. We estimate that ∼12% and 1% of presolar silicates have SN and nova origins, respectively, similar to presolar SiC and oxides. This implies that asymptotic giant branch stars are the dominant dust producers in the galaxy.

We have examined Ulysses magnetic field data using dynamic spectrogram techniques that compute wave amplitude, polarization, and direction of propagation over a broad range of frequencies and time. Events were identified that showed a strong polarization signature and an enhancement of power above the local proton gyrofrequency. We perform a statistical study of 502 wave events in an effort to determine when, where, and why they are observed. Most notably, we find that waves arising from newborn interstellar pickup ions are relatively rare and difficult to find. The quantities normally employed in theories of wave growth are neutral atom density and quantities related to their ionization and the subsequent dynamics such as wind speed, solar wind flux, and magnetic field orientation. We find the observations of waves to be largely uncorrelated to these quantities except for mean field direction where quasi-radial magnetic fields are favored and solar wind proton flux where wave observations appear to be favored by low flux conditions which runs contrary to theoretical expectations of wave generation. It would appear that an explanation based on source physics and instability growth rates alone is not adequate to account for the times when these waves are seen.

We have obtained WIYN and SOAR i' images of BL Lacertae objects and used these to detect or constrain the flux of the host galaxy. Under common standard candle assumptions, these data provide estimates of, or lower bounds on, the redshift. Our targets are a set of flat-spectrum radio counterparts of high flux Fermi Large Area Telescope sources, with sensitive spectral observations showing them to be continuum-dominated BL Lac objects. In this sample, 5 of 11 BL Lac objects yielded significant host detections, with standard candle redshifts z = 0.13–0.58. Our estimates and lower bounds are generally in agreement with other redshifts estimates, although our z = 0.374 estimate for J0543−5532 implies a significantly sub-luminous host.

Using deep narrowband H₂S1 and K_s-band imaging data obtained with CFHT/WIRCam, we identify a sample of 56 Hα emission-line galaxies (ELGs) at z = 2.24 with the 5σ depths of H₂S1 = 22.8 and K_s = 24.8 (AB) over a 383 arcmin² area in the Extended Chandra Deep Field South. A detailed analysis is carried out with existing multi-wavelength data in this field. Three of the 56 Hα ELGs are detected in Chandra 4 Ms X-ray observations and two of them are classified as active galactic nuclei. The rest-frame UV and optical morphologies revealed by HST/ACS and WFC3 deep images show that nearly half of the Hα ELGs are either merging systems or have a close companion, indicating that the merging/interacting processes play a key role in regulating star formation at cosmic epoch z = 2–3. About 14% are too faint to be resolved in the rest-frame UV morphology due to high dust extinction. We estimate dust extinction from spectral energy distributions. We find that dust extinction is generally correlated with Hα luminosity and stellar mass. Our results suggest that Hα ELGs are representative of star-forming galaxies. Applying extinction corrections to individual objects, we examine the intrinsic Hα luminosity function (LF) at z = 2.24, obtaining a best-fit Schechter function characterized by a faint-end slope of α = − 1.3. This is shallower than the typical slope of α ≃ −1.6 in previous works based on constant extinction correction. We demonstrate that this difference is mainly due to the different extinction corrections. The proper extinction correction is thus the key to recovering the intrinsic LF as the extinction globally increases with Hα luminosity. Moreover, we find that our Hα LF mirrors the stellar mass function of star-forming galaxies at the same cosmic epoch. This finding indeed reflects the tight correlation between star formation rate and stellar mass for star-forming galaxies, i.e., the so-called main sequence.

The entire evolution of the Milky Way, including its mass-assembly and star-formation history, is imprinted onto the chemo-dynamical distribution function of its member stars, f(x, v, [X/H]), in the multi-dimensional phase space spanned by position, velocity, and elemental abundance ratios. In particular, the chemo-dynamical distribution functions for low-mass stars (e.g., G- or K-type dwarfs) are precious tracers of the earliest stages of the Milky Way's formation, since their main-sequence lifetimes approach or exceed the age of the universe. A basic tenet of essentially all previous analyses is that the stellar metallicity, usually parameterized as [Fe/H], is conserved over time for main-sequence stars (at least those that have not been polluted due to mass transfer from binary companions). If this holds true, any correlations between metallicity and kinematics for long-lived main-sequence stars of different masses, effective temperatures, or spectral types must strictly be the same, since they reflect the same mass-assembly and star-formation histories. By analyzing a sample of nearby metal-poor halo and thick-disk stars on the main sequence, taken from Data Release 8 of the Sloan Digital Sky Survey, we find that the median metallicity of G-type dwarfs is systematically higher (by about 0.2 dex) than that of K-type dwarfs having the same median rotational velocity about the Galactic center. If it can be confirmed, this finding may invalidate the long-accepted assumption that the atmospheric metallicities of long-lived stars are conserved over time.

We report on a survey for Hα bow shock emission around nearby γ-detected energetic pulsars. This survey adds three Balmer-dominated neutron star bow shocks to the six previously confirmed examples. In addition to the shock around Fermi pulsar PSR J1741−2054, we now report Hα structures around two additional γ-ray pulsars, PSR J2030+4415 and PSR J1509−5850. These are the first known examples of Hα nebulae with pre-ionization halos. With new measurements, we show that a simple analytic model can account for the angular size and flux of the bow shocks' apices. The latter, in particular, provides a new pulsar probe and indicates large moments of inertia and smaller distances than previously assumed in several cases. In particular, we show that the re-measured PSR J0437−4715 shock flux implies I = (1.7 ± 0.2) × 10⁴⁵/(f_H isin i) g cm². We also derive a distance d ≈ 0.72 kpc for the γ-ray only pulsar PSR J2030+4415 and revised distances for PSRs J1959+2048 (1.4 kpc) and J2555+6535 (∼1 kpc), smaller than the conventional DM-estimated values. Finally, we report upper limits for 94 additional LAT pulsars. An estimate of the survey sensitivity indicates that for a warm neutral medium filling factor ϕ_WNM ∼ 0.3 there should be a total of approximately nine Hα bow shocks in our LAT-targeted survey; given that seven such objects are now known, a much larger ϕ_WNM seems problematic.

In the spirit of minimal modeling of complex systems, we develop an idealized two-column model to investigate the climate of tidally locked terrestrial planets with Earth-like atmospheres in the habitable zone of M-dwarf stars. The model is able to approximate the fundamental features of the climate obtained from three-dimensional (3D) atmospheric general circulation model (GCM) simulations. One important reason for the two-column model's success is that it reproduces the high cloud albedo of the GCM simulations, which reduces the planet's temperature and delays the onset of a runaway greenhouse state. The two-column model also clearly illustrates a secondary mechanism for determining the climate: the nightside acts as a "radiator fin" through which infrared energy can be lost to space easily. This radiator fin is maintained by a temperature inversion and dry air on the nightside, and plays a similar role to the subtropics on modern Earth. Since one-dimensional radiative–convective models cannot capture the effects of the cloud albedo and radiator fin, they are systematically biased toward a narrower habitable zone. We also show that cloud parameters are the most important in the two-column model for determining the day–night thermal emission contrast, which decreases and eventually reverses as the stellar flux increases. This reversal is important because it could be detected by future extrasolar planet characterization missions, which would suggest that the planet has Earth-like water clouds and is potentially habitable.

The MEarth survey is a search for small rocky planets around the smallest, nearest stars to the Sun as identified by high proper motion with red colors. We augmented our planetary search time series with lower cadence astrometric imaging and obtained two million images of approximately 1800 stars suspected to be mid-to-late M dwarfs. We fit an astrometric model to MEarth's images for 1507 stars and obtained trigonometric distance measurements to each star with an average precision of 5 mas. Our measurements, combined with the Two Micron All Sky Survey photometry, allowed us to obtain an absolute K_s magnitude for each star. In turn, this allows us to better estimate the stellar parameters than those obtained with photometric estimates alone and to better prioritize the targets chosen to monitor at high cadence for planetary transits. The MEarth sample is mostly complete out to a distance of 25 pc for stars of type M5.5V and earlier, and mostly complete for later type stars out to 20 pc. We find eight stars that are within 10 pc of the Sun for which there did not exist a published trigonometric parallax distance estimate. We release with this work a catalog of the trigonometric parallax measurements for 1507 mid-to-late M dwarfs, as well as new estimates of their masses and radii.

Using a combination of high-resolution Hubble Space Telescope/Wide-Field and Planetary Camera-2 observations, we explore the physical properties of the stellar populations in two intermediate-age star clusters, NGC 1831 and NGC 1868, in the Large Magellanic Cloud based on their color–magnitude diagrams. We show that both clusters exhibit extended main-sequence turn offs. To explain the observations, we consider variations in helium abundance, binarity, age dispersions, and the fast rotation of the clusters' member stars. The observed narrow main sequence excludes significant variations in helium abundance in both clusters. We first establish the clusters' main-sequence binary fractions using the bulk of the clusters' main-sequence stellar populations ≳ 1 mag below their turn-offs. The extent of the turn-off regions in color–magnitude space, corrected for the effects of binarity, implies that age spreads of order 300 Myr may be inferred for both clusters if the stellar distributions in color–magnitude space were entirely due to the presence of multiple populations characterized by an age range. Invoking rapid rotation of the population of cluster members characterized by a single age also allows us to match the observed data in detail. However, when taking into account the extent of the red clump in color–magnitude space, we encounter an apparent conflict for NGC 1831 between the age dispersion derived from that based on the extent of the main-sequence turn off and that implied by the compact red clump. We therefore conclude that, for this cluster, variations in stellar rotation rate are preferred over an age dispersion. For NGC 1868, both models perform equally well.

Recent studies suggest that metal-poor stars enhanced in carbon but containing low levels of neutron-capture elements may have been among the first to incorporate the nucleosynthesis products of the first generation of stars. We have observed 16 stars with enhanced carbon or nitrogen using the MIKE Spectrograph on the Magellan Telescopes at Las Campanas Observatory and the Tull Spectrograph on the Smith Telescope at McDonald Observatory. We present radial velocities, stellar parameters, and detailed abundance patterns for these stars. Strontium, yttrium, zirconium, barium, europium, ytterbium, and other heavy elements are detected. In four stars, these heavy elements appear to have originated in some form of r-process nucleosynthesis. In one star, a partial s-process origin is possible. The origin of the heavy elements in the rest of the sample cannot be determined unambiguously. The presence of elements heavier than the iron group offers further evidence that zero-metallicity rapidly rotating massive stars and pair instability supernovae did not contribute substantial amounts of neutron-capture elements to the regions where the stars in our sample formed. If the carbon- or nitrogen-enhanced metal-poor stars with low levels of neutron-capture elements were enriched by products of zero-metallicity supernovae only, then the presence of these heavy elements indicates that at least one form of neutron-capture reaction operated in some of the first stars.

We investigate the hypothesis that γ-ray-quiet active galactic nuclei (AGNs) have a greater tendency for jet bending than γ-ray-loud AGNs, revisiting the analysis of Tingay et al. We perform a statistical analysis using a large sample of 351 radio-loud AGNs along with γ-ray identifications from the Fermi Large Area Telescope (LAT). Our results show no statistically significant differences in jet-bending properties between γ-ray-loud and γ-ray-quiet populations, indicating that jet bending is not a significant factor for γ-ray detection in AGNs.

As a planet ages, it cools and its radius shrinks at a rate set by the efficiency with which heat is transported from the interior out to space. The bottleneck for this transport is at the boundary between the convective interior and the radiative atmosphere; the opacity there sets the global cooling rate. Models of planetary evolution are often one dimensional (1D), such that the radiative–convective boundary (RCB) is defined by a single temperature, pressure, and opacity. In reality the spatially inhomogeneous stellar heating pattern and circulation in the atmosphere could deform the RCB, allowing heat from the interior to escape more efficiently through regions with lower opacity. We present an analysis of the degree to which the RCB could be deformed and the resultant change in the evolutionary cooling rate. In this initial work we calculate the upper limit for this effect by comparing an atmospheric structure in local radiative equilibrium to its 1D equivalent. We find that the cooling through an uneven RCB could be enhanced over cooling through a uniform RCB by as much as 10%–50%. We also show that the deformation of the RCB (and the enhancement of the cooling rate) increases with a greater incident stellar flux or a lower inner entropy. Our results indicate that this mechanism could significantly change a planet's thermal evolution, causing it to cool and shrink more quickly than would otherwise be expected. This may exacerbate the well-known difficulty in explaining the very large radii observed for some hot Jupiters.

We present an analysis of the effects of dissipational baryonic physics on the local dark matter (DM) distribution at the location of the Sun, with an emphasis on the consequences for direct detection experiments. Our work is based on a comparative analysis of two cosmological simulations with identical initial conditions of a Milky Way halo, one of which (Eris) is a full hydrodynamic simulation and the other (ErisDark) is a DM-only one. We find that in Eris two distinct processes lead to a 30% enhancement of DM in the disk plane at the location of the Sun: the accretion and disruption of satellites resulting in a DM component with net angular momentum, and the contraction of baryons pulling the DM into the disk plane without forcing it to co-rotate. Owing to its particularly quiescent merger history for dark halos of Milky Way mass, the co-rotating dark disk in Eris is less massive than what has been suggested by previous work, contributing only 9% of the local DM density. Yet, since the simulation results in a realistic Milky Way analog galaxy, its DM halo provides a plausible alternative to the Maxwellian standard halo model (SHM) commonly used in direct detection analyses. The speed distribution in Eris is broadened and shifted to higher speeds, compared to its DM-only twin simulation ErisDark. At high speeds f(v) falls more steeply in Eris than in ErisDark or the SHM, easing the tension between recent results from the CDMS-II and XENON100 experiments. The non-Maxwellian aspects of f(v) are still present, but much less pronounced in Eris than in the DM-only runs. The weak dark disk increases the time-averaged scattering rate by only a few percent at low recoil energies. On the high velocity tail, however, the increase in typical speeds due to baryonic contraction results in strongly enhanced mean scattering rates compared to ErisDark, although they are still suppressed compared to the SHM. Similar trends are seen regarding the amplitude of the annual modulation, while the modulated fraction is increased compared to the SHM and decreased compared to ErisDark.

A number of recent studies have proposed that the stellar initial mass function (IMF) of early type galaxies varies systematically as a function of galaxy mass, with higher mass galaxies having bottom-heavy IMFs. These bottom-heavy IMFs have more low-mass stars relative to the number of high mass stars, and therefore naturally result in proportionally fewer neutron stars (NSs) and black holes (BHs). In this paper, we specifically predict the variation in the number of BHs and NSs based on the power-law IMF variation required to reproduce the observed mass-to-light ratio trends with galaxy mass. We then test whether such variations are observed by studying the field low-mass X-ray binary (LMXB) populations of nearby early-type galaxies. In these binaries, an NS or BH accretes matter from a low-mass donor star. Their number is therefore expected to scale with the number of BHs and NSs present in a galaxy. We find that the number of LMXBs per K-band light is similar among the galaxies in our sample. These data therefore demonstrate the uniformity of the slope of the IMF from massive stars down to those now dominating the K-band light and are consistent with an invariant IMF. Our results are inconsistent with an IMF which varies from a Kroupa/Chabrier like IMF for low-mass galaxies to a steep power-law IMF (with slope x = 2.8) for high mass galaxies. We discuss how these observations constrain the possible forms of the IMF variations and how future Chandra observations can enable sharper tests of the IMF.

Two techniques, the Partial Variance of Increments (PVI) and the Local Intermittency Measure (LIM), have been applied and compared using MESSENGER magnetic field data in the solar wind at a heliocentric distance of about 0.3 AU. The spatial properties of the turbulent field at different scales, spanning the whole inertial range of magnetic turbulence down toward the proton scales have been studied. LIM and PVI methodologies allow us to identify portions of an entire time series where magnetic energy is mostly accumulated, and regions of intermittent bursts in the magnetic field vector increments, respectively. A statistical analysis has revealed that at small time scales and for high level of the threshold, the bursts present in the PVI and the LIM series correspond to regions of high shear stress and high magnetic field compressibility.

In the recent paper of Antiochos, a new concept for the injection of magnetic helicity into the solar corona by small-scale convective motions and its condensation onto polarity inversion lines (PILs) was developed. We investigate this concept through global simulations of the Sun's photospheric and coronal magnetic fields, and compare the results with the hemispheric pattern of solar filaments. Assuming that the vorticity of the cells is predominantly counterclockwise/clockwise in the northern/southern hemisphere, the convective motions inject negative/positive helicity into each hemisphere. The simulations show that: (1) on a north–south oriented PIL, both differential rotation and convective motions inject the same sign of helicity, which matches that required to reproduce the hemispheric pattern of filaments. (2) On a high-latitude east–west oriented polar crown or subpolar crown PIL, the vorticity of the cells has to be approximately 2–3 times greater than the local differential-rotation gradient in order to overcome the incorrect sign of helicity injection from differential rotation. (3) In the declining phase of the cycle, as a bipole interacts with the polar field, in some cases, helicity condensation can reverse the effect of differential rotation along the east–west lead arm but not in all cases. The results show that this newly developed concept of magnetic helicity injection and condensation, in conjunction with the mechanisms used in Yeates et al., is a viable explanation for the hemispheric pattern of filaments. Future observational studies should focus on examining the vorticity component within convective motions to determine both its magnitude and latitudinal variation relative to the differential-rotation gradient on the Sun.

We present the observation of a major solar eruption that is associated with fast sunspot rotation. The event includes a sigmoidal filament eruption, a coronal mass ejection, and a GOES X2.1 flare from NOAA active region 11283. The filament and some overlying arcades were partially rooted in a sunspot. The sunspot rotated at ∼10° hr⁻¹ during a period of 6 hr prior to the eruption. In this period, the filament was found to rise gradually along with the sunspot rotation. Based on the Helioseismic and Magnetic Imager observation, for an area along the polarity inversion line underneath the filament, we found gradual pre-eruption decreases of both the mean strength of the photospheric horizontal field (B_h) and the mean inclination angle between the vector magnetic field and the local radial (or vertical) direction. These observations are consistent with the pre-eruption gradual rising of the filament-associated magnetic structure. In addition, according to the nonlinear force-free field reconstruction of the coronal magnetic field, a pre-eruption magnetic flux rope structure is found to be in alignment with the filament, and a considerable amount of magnetic energy was transported to the corona during the period of sunspot rotation. Our study provides evidence that in this event sunspot rotation plays an important role in twisting, energizing, and destabilizing the coronal filament–flux rope system, and led to the eruption. We also propose that the pre-event evolution of B_h may be used to discern the driving mechanism of eruptions.

Images recorded by the X-ray Telescope on board the Hinode spacecraft are used to provide high-cadence observations of solar jetting activity. A selection of the brightest of these polar jets shows a positive correlation with high-speed responses traced into the interplanetary medium. LASCO C2 and STEREO COR2 coronagraph images measure the coronal response to some of the largest jets, and also the nearby background solar wind velocity, thereby giving a determination of their speeds that we compare with Hinode observations. When using the full Solar Mass Ejection Imager (SMEI) data set, we track these same high-speed solar jet responses into the inner heliosphere and from these analyses determine their mass, flow energies, and the extent to which they retain their identity at large solar distances.

We investigated the influence of velocity shear and a radial density profile on the spatial development of the current-driven (CD) kink instability along helically magnetized relativistic jets via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we use a nonperiodic computational box, the jet flow is initially established across the computational grid, and a precessional perturbation at the inlet triggers growth of the kink instability. If the velocity shear radius is located inside the characteristic radius of the helical magnetic field, a static nonpropagating CD kink is excited as the perturbation propagates down the jet. Temporal growth disrupts the initial flow across the computational grid not too far from the inlet. On the other hand, if the velocity shear radius is outside the characteristic radius of the helical magnetic field, the kink is advected with the flow and grows spatially down the jet. In this case, flow is maintained to much larger distances from the inlet. The effect of different radial density profiles is more subtle. When the density increases with radius, the kink appears to saturate by the end of the simulation without apparent disruption of the helical twist. This behavior suggests that relativistic jets consisting of a tenuous spine surrounded by a denser medium with a velocity shear radius outside the radius of the maximum toroidal magnetic field have a relatively stable configuration.

We propose a catastrophic eruption model for the enormous energy release of magnetars during giant flares, in which a toroidal and helically twisted flux rope is embedded within a force-free magnetosphere. The flux rope stays in stable equilibrium states initially and evolves quasi-statically. Upon the loss of equilibrium, the flux rope cannot sustain the stable equilibrium states and erupts catastrophically. During the process, the magnetic energy stored in the magnetosphere is rapidly released as the result of destabilization of global magnetic topology. The magnetospheric energy that could be accumulated is of vital importance for the outbursts of magnetars. We carefully establish the fully open fields and partially open fields for various boundary conditions at the magnetar surface and study the relevant energy thresholds. By investigating the magnetic energy accumulated at the critical catastrophic point, we find that it is possible to drive fully open eruptions for dipole-dominated background fields. Nevertheless, it is hard to generate fully open magnetic eruptions for multipolar background fields. Given the observational importance of the multipolar magnetic fields in the vicinity of the magnetar surface, it would be worthwhile to explore the possibility of the alternative eruption approach in multipolar background fields. Fortunately, we find that flux ropes may give rise to partially open eruptions in the multipolar fields, which involve only partial opening of background fields. The energy release fractions are greater for cases with central-arcaded multipoles than those with central-caved multipoles that emerged in background fields. Eruptions would fail only when the centrally caved multipoles become extremely strong.

A new mechanism to form a magnetic pressure supported, high temperature corona above the photosphere of an accretion disk is explored using three dimensional radiation magnetohydrodynamic (MHD) simulations. The thermal properties of the disk are calculated self-consistently by balancing radiative cooling through the surfaces of the disk with heating due to dissipation of turbulence driven by magneto-rotational instability (MRI). As has been noted in previous work, we find the dissipation rate per unit mass increases dramatically with height above the mid-plane, in stark contrast to the α-disk model which assumes this quantity is a constant. Thus, we find that in simulations with a low surface density (and therefore a shallow photosphere), the fraction of energy dissipated above the photosphere is significant (about 3.4% in our lowest surface density model), and this fraction increases as surface density decreases. When a significant fraction of the accretion energy is dissipated in the optically thin photosphere, the gas temperature increases substantially and a high temperature, magnetic pressure supported corona is formed. The volume-averaged temperature in the disk corona is more than 10 times larger than at the disk mid-plane. Moreover, gas temperature in the corona is strongly anti-correlated with gas density, which implies the corona formed by MRI turbulence is patchy. This mechanism to form an accretion disk corona may help explain the observed relation between the spectral index and luminosity from active galactic nucleus (AGNs), and the soft X-ray excess from some AGNs. It may also be relevant to spectral state changes in X-ray binaries.

We present an online catalog of distance determinations for 6036 K giants, most of which are members of the Milky Way's stellar halo. Their medium-resolution spectra from the Sloan Digital Sky Survey/Sloan Extension for Galactic Understanding and Exploration are used to derive metallicities and rough gravity estimates, along with radial velocities. Distance moduli are derived from a comparison of each star's apparent magnitude with the absolute magnitude of empirically calibrated color–luminosity fiducials, at the observed (g − r)₀ color and spectroscopic [Fe/H]. We employ a probabilistic approach that makes it straightforward to properly propagate the errors in metallicities, magnitudes, and colors into distance uncertainties. We also fold in prior information about the giant-branch luminosity function and the different metallicity distributions of the SEGUE K-giant targeting sub-categories. We show that the metallicity prior plays a small role in the distance estimates, but that neglecting the luminosity prior could lead to a systematic distance modulus bias of up to 0.25 mag, compared to the case of using the luminosity prior. We find a median distance precision of 16%, with distance estimates most precise for the least metal-poor stars near the tip of the red giant branch. The precision and accuracy of our distance estimates are validated with observations of globular and open clusters. The stars in our catalog are up to 125 kpc from the Galactic center, with 283 stars beyond 50 kpc, forming the largest available spectroscopic sample of distant tracers in the Galactic halo.

We present IRAM Plateau de Bure Interferometer observations of the CO(3–2) and CO(5–4) line transitions from an Lyα blob at z ∼ 2.7 in order to investigate the gas kinematics, determine the location of the dominant energy source, and study the physical conditions of the molecular gas. CO line and dust continuum emissions are detected at the location of a strong MIPS source that is offset by ∼15 from the Lyα peak. Neither of these emission components is resolved with the 17 beam, showing that the gas and dust are confined to within ∼7 kpc from this galaxy. No millimeter source is found at the location of the Lyα peak, ruling out a central compact source of star formation as the power source for the Lyα emission. Combined with a spatially resolved spectrum of Lyα and He ii, we constrain the kinematics of the extended gas using the CO emission as a tracer of the systemic redshift. Near the MIPS source, the Lyα profile is symmetric, and its line center agrees with that of the CO line, implying that there are no significant bulk flows and that the photo-ionization from the MIPS source might be the dominant source of the Lyα emission. In the region near the Lyα peak, the gas is slowly receding (∼100 km s⁻¹) with respect to the MIPS source, thus making the hyper-/superwind hypothesis unlikely. We find a sub-thermal line ratio between two CO transitions, I_CO(5-4)/I_CO(3-2) = 0.97 ± 0.21. This line ratio is lower than the average values found in high-z submillimeter galaxies and QSOs but is consistent with the value found in the Galactic center, suggesting that there is a large reservoir of low-density molecular gas that is spread over the MIPS source and its vicinity.

Broad infrared emission features (e.g., at 3.3, 6.2, 7.7, 8.6, and 11.3 μm) from the gas phase interstellar medium have long been attributed to polycyclic aromatic hydrocarbons (PAHs). A significant portion (10%–20%) of the Milky Way's carbon reservoir is locked in PAH molecules, which makes their characterization integral to our understanding of astrochemistry. In molecular clouds and the dense envelopes and disks of young stellar objects (YSOs), PAHs are expected to be frozen in the icy mantles of dust grains where they should reveal themselves through infrared absorption. To facilitate the search for frozen interstellar PAHs, laboratory experiments were conducted to determine the positions and strengths of the bands of pyrene mixed with H₂O and D₂O ices. The D₂O mixtures are used to measure pyrene bands that are masked by the strong bands of H₂O, leading to the first laboratory determination of the band strength for the CH stretching mode of pyrene in water ice near 3.25 μm. Our infrared band strengths were normalized to experimentally determined ultraviolet band strengths, and we find that they are generally ∼50% larger than those reported by Bouwman et al. based on theoretical strengths. These improved band strengths were used to reexamine YSO spectra published by Boogert et al. to estimate the contribution of frozen PAHs to absorption in the 5–8 μm spectral region, taking into account the strength of the 3.25 μm CH stretching mode. It is found that frozen neutral PAHs contain 5%–9% of the cosmic carbon budget and account for 2%–9% of the unidentified absorption in the 5–8 μm region.

We study the dust present in 56 Galactic planetary nebulae (PNe) through their iron depletion factors, their C/O abundance ratios (in 51 objects), and the dust features that appear in their infrared spectra (for 33 objects). Our sample objects have deep optical spectra of good quality, and most of them also have ultraviolet observations. We use these observations to derive the iron abundances and the C/O abundance ratios in a homogeneous way for all the objects. We compile detections of infrared dust features from the literature and we analyze the available Spitzer/IRS spectra. Most of the PNe have C/O ratios below one and show crystalline silicates in their infrared spectra. The PNe with silicates have C/O <1, with the exception of Cn 1-5. Most of the PNe with dust features related to C-rich environments (SiC or the 30 μm feature usually associated to MgS) have C/O ≳ 0.8. Polycyclic aromatic hydrocarbons are detected over the full range of C/O values, including 6 objects that also show silicates. Iron abundances are low in all the objects, implying that more than 90% of their iron atoms are deposited into dust grains. The range of iron depletions in the sample covers about two orders of magnitude, and we find that the highest depletion factors are found in C-rich objects with SiC or the 30 μm feature in their infrared spectra, whereas some of the O-rich objects with silicates show the lowest depletion factors.

We extend a model for turbulence–flame interactions (TFI) to consider astrophysical flames with a particular focus on combustion in Type Ia supernovae. The inertial range of the turbulent cascade is nearly always under-resolved in simulations of astrophysical flows, requiring the use of a model in order to quantify the effects of subgrid-scale wrinkling of the flame surface. We provide implementation details to extend a well-tested TFI model to low-Prandtl number flames for use in the compressible hydrodynamics code flash. A local, instantaneous measure of the turbulent velocity is calibrated for flash and verification tests are performed. Particular care is taken to consider the relation between the subgrid rms turbulent velocity and the turbulent flame speed, especially for high-intensity turbulence where the turbulent flame speed is not expected to scale with the turbulent velocity. Finally, we explore the impact of different TFI models in full-star, three-dimensional simulations of Type Ia supernovae.