Table of contents

Transiting planet surveys like Kepler have provided a wealth of information on the distribution of planetary radii, particularly for the new populations of super-Earth- and sub-Neptune-sized planets. In order to aid in the physical interpretation of these radii, we compute model radii for low-mass rocky planets with hydrogen–helium envelopes. We provide model radii for planets 1–20 M_⊕, with envelope fractions 0.01%–20%, levels of irradiation 0.1–1000 times Earth's, and ages from 100 Myr to 10 Gyr. In addition we provide simple analytic fits that summarize how radius depends on each of these parameters. Most importantly, we show that at fixed H/He envelope fraction, radii show little dependence on mass for planets with more than ∼1% of their mass in their envelope. Consequently, planetary radius is to a first order a proxy for planetary composition, i.e., H/He envelope fraction, for Neptune- and sub-Neptune-sized planets. We recast the observed mass–radius relationship as a mass–composition relationship and discuss it in light of traditional core accretion theory. We discuss the transition from rocky super-Earths to sub-Neptune planets with large volatile envelopes. We suggest ∼1.75 R_⊕ as a physically motivated dividing line between these two populations of planets. Finally, we discuss these results in light of the observed radius occurrence distribution found by Kepler.

Water vapor has been detected in protoplanetary disks. In this work, we model the distribution of water vapor in protoplanetary disks with a thermo-chemical code. For a set of parameterized disk models, we calculate the distribution of dust temperature and radiation field of the disk with a Monte Carlo method, and then solve the gas temperature distribution and chemical composition. The radiative transfer includes detailed treatment of scattering by atomic hydrogen and absorption by water of Lyα photons, since the Lyα line dominates the UV spectrum of accreting young stars. In a fiducial model, we find that warm water vapor with temperature around 300 K is mainly distributed in a small and well-confined region in the inner disk. The inner boundary of the warm water region is where the shielding of UV field due to dust and water itself become significant. The outer boundary is where the dust temperature drops below the water condensation temperature. A more luminous central star leads to a more extended distribution of warm water vapor, while dust growth and settling tends to reduce the amount of warm water vapor. Based on typical assumptions regarding the elemental oxygen abundance and the water chemistry, the column density of warm water vapor can be as high as 10²² cm⁻². A small amount of hot water vapor with temperature higher than ∼300 K exists in a more extended region in the upper atmosphere of the disk. Cold water vapor with temperature lower than 100 K is distributed over the entire disk, produced by photodesorption of the water ice.

In the local universe, there is good evidence that, at a given stellar mass M, the gas-phase metallicity Z is anti-correlated with the star formation rate (SFR) of the galaxies. It has also been claimed that the resulting Z(M, SFR) relation is invariant with redshift—the so-called "fundamental metallicity relation" (FMR). Given a number of difficulties in determining metallicities, especially at higher redshifts, the form of the Z(M, SFR) relation and whether it is really independent of redshift is still very controversial. To explore this issue at z > 2, we used VLT-SINFONI and Subaru-MOIRCS near-infrared spectroscopy of 20 zCOSMOS-deep galaxies at 2.1 < z < 2.5 to measure the strengths of up to five emission lines: [O ii] λ3727, Hβ, [O iii] λ5007, Hα, and [N ii] λ6584. This near-infrared spectroscopy enables us to derive O/H metallicities, and also SFRs from extinction corrected Hα measurements. We find that the mass–metallicity relation (MZR) of these star-forming galaxies at z ≈ 2.3 is lower than the local Sloan Digital Sky Survey (SDSS) MZR by a factor of three to five, a larger change than found by Erb et al. using [N ii]/Hα-based metallicities from stacked spectra. We discuss how the different selections of the samples and metallicity calibrations used may be responsible for this discrepancy. The galaxies show direct evidence that the SFR is still a second parameter in the MZR at these redshifts. However, determining whether the Z(M, SFR) relation is invariant with epoch depends on the choice of extrapolation used from local samples, because z > 2 galaxies of a given mass have much higher SFRs than the local SDSS galaxies. We find that the zCOSMOS galaxies are consistent with a non-evolving FMR if we use the physically motivated formulation of the Z(M, SFR) relation from Lilly et al., but not if we use the empirical formulation of Mannucci et al.

Burst oscillations during the rising phases of thermonuclear X-ray bursts are usually believed to originate from flame spreading on the neutron star surface. However, the decrease of fractional oscillation amplitude with rise time, which provides a main observational support for the flame spreading model, have so far been reported from only a few bursts. Moreover, the non-detection and intermittent detections of rise oscillations from many bursts are not yet understood considering the flame spreading scenario. Here, we report the decreasing trend of fractional oscillation amplitude from an extensive analysis of a large sample of Rossi X-ray Timing Explorer Proportional Counter Array bursts from 10 neutron star low-mass X-ray binaries. This trend is 99.99% significant for the best case, which provides, to the best of our knowledge, by far the strongest evidence of such a trend. Moreover, it is important to note that an opposite trend is not found in any of the bursts. The concave shape of the fractional amplitude profiles for all the bursts suggests latitude-dependent flame speeds, possibly due to the effects of the Coriolis force. We also systematically study the roles of low fractional amplitude and low count rate for non-detection and intermittent detections of rise oscillations, and attempt to understand them within the flame spreading scenario. Our results support a weak turbulent viscosity for flame spreading, and imply that burst rise oscillations originate from an expanding hot spot, thus making these oscillations a more reliable tool to constrain the neutron star equations of state.

Electrons accelerated in solar flares are injected into the surrounding plasma, where they are subjected to the influence of collisional (Coulomb) energy losses. Their evolution is modeled by a partial differential equation describing continuity of electron number. In a recent paper, Dobranskis & Zharkova claim to have found an "updated exact analytical solution" to this continuity equation. Their solution contains an additional term that drives an exponential decrease in electron density with depth, leading them to assert that the well-known solution derived by Brown, Syrovatskii & Shmeleva, and many others is invalid. We show that the solution of Dobranskis & Zharkova results from a fundamental error in the application of the method of characteristics and is hence incorrect. Further, their comparison of the "new" analytical solution with numerical solutions of the Fokker–Planck equation fails to lend support to their result. We conclude that Dobranskis & Zharkova's solution of the universally accepted and well-established continuity equation is incorrect, and that their criticism of the correct solution is unfounded. We also demonstrate the formal equivalence of the approaches of Syrovatskii & Shmeleva and Brown, with particular reference to the evolution of the electron flux and number density (both differential in energy) in a collisional thick target. We strongly urge use of these long-established, correct solutions in future works.

Attempts to explain the source of r-process elements in our solar system (S.S.) by particular astrophysical sites still face entwined uncertainties, stemming from the extrapolation of nuclear properties far from stability, inconsistent sources of different properties (e.g., nuclear masses and β-decay properties), and the poor understanding of astrophysical conditions, which are hard to disentangle. In this paper we present results from the investigation of r-process in the high-entropy wind (HEW) of core-collapse supernovae (here chosen as one of the possible scenarios for this nucleosynthesis process), using new nuclear-data input calculated in a consistent approach, for masses and β-decay properties from the new finite-range droplet model FRDM(2012). The accuracy of the new mass model is 0.56 MeV with respect to AME2003, to which it was adjusted. We compare the new HEW r-process abundance pattern to the latest S.S. r-process residuals and to our earlier calculations with the nuclear-structure quantities based on FRDM(1992). Substantial overall and specific local improvements in the calculated pattern of the r-process between A ≃ 110 and ²⁰⁹Bi, as well as remaining deficiencies, are discussed in terms of the underlying spherical and deformed shell structure far from stability.

We present extensive optical and near-infrared photometric and spectroscopic observations of the stripped-envelope supernova SN 2010as. Spectroscopic peculiarities such as initially weak helium features and low expansion velocities with a nearly flat evolution place this object in the small family of events previously identified as transitional Type Ib/c supernovae (SNe). There is ubiquitous evidence of hydrogen, albeit weak, in this family of SNe, indicating that they are in fact a peculiar kind of Type IIb SNe that we name "flat-velocity Type IIb. The flat-velocity evolution—which occurs at different levels between 6000 and 8000 km s⁻¹ for different SNe—suggests the presence of a dense shell in the ejecta. Despite the spectroscopic similarities, these objects show surprisingly diverse luminosities. We discuss the possible physical or geometrical unification picture for such diversity. Using archival Hubble Space Telescope images, we associate SN 2010as with a massive cluster and derive a progenitor age of ≈6 Myr, assuming a single star-formation burst, which is compatible with a Wolf–Rayet progenitor. Our hydrodynamical modeling, on the contrary, indicates that the pre-explosion mass was relatively low, ≈4 M_☉. The seeming contradiction between a young age and low pre-SN mass may be solved by a massive interacting binary progenitor.

We analyze the physical conditions of the cool, photoionized (T ∼10⁴ K) circumgalactic medium (CGM) using the COS-Halos suite of gas column density measurements for 44 gaseous halos within 160 kpc of L ∼ L* galaxies at z ∼ 0.2. These data are well described by simple photoionization models, with the gas highly ionized (n_H ii/n_H ≳ 99%) by the extragalactic ultraviolet background. Scaling by estimates for the virial radius, R_vir, we show that the ionization state (tracked by the dimensionless ionization parameter, U) increases with distance from the host galaxy. The ionization parameters imply a decreasing volume density profile n_H = (10^{−4.2 ± 0.25})(R/R_vir)^{−0.8 ± 0.3}. Our derived gas volume densities are several orders of magnitude lower than predictions from standard two-phase models with a cool medium in pressure equilibrium with a hot, coronal medium expected in virialized halos at this mass scale. Applying the ionization corrections to the H i column densities, we estimate a lower limit to the cool gas mass $M_{\rm CGM}^{\rm cool} \gt 6.5 \times 10^{10}$ M_☉ for the volume within R < R_vir. Allowing for an additional warm-hot, O vi-traced phase, the CGM accounts for at least half of the baryons purported to be missing from dark matter halos at the 10¹² M_☉ scale.

We perform hybrid simulations of a super-Alfvénic quasi-parallel shock, driven by a coronal mass ejection (CME), propagating in the outer coronal/solar wind at distances of between 3 to 6 solar radii. The hybrid treatment of the problem enables the study of the shock propagation on the ion timescale, preserving ion kinetics and allowing for a self-consistent treatment of the shock propagation and particle acceleration. The CME plasma drags the embedded magnetic field lines stretching from the sun, and propagates out into interplanetary space at a greater velocity than the in situ solar wind, driving the shock, and producing very energetic particles. Our results show that electromagnetic Alfvén waves are generated at the shock front. The waves propagate upstream of the shock and are produced by the counter-streaming ions of the solar wind plasma being reflected at the shock. A significant fraction of the particles are accelerated in two distinct phases: first, particles drift from the shock and are accelerated in the upstream region, and second, particles arriving at the shock get trapped and are accelerated at the shock front. A fraction of the particles diffused back to the shock, which is consistent with the Fermi acceleration mechanism.

We revisit the evidence for the contribution of the long-lived radioactive nuclides ⁴⁴Ti, ⁵⁵Fe, ⁵⁶Co, ⁵⁷Co, and ⁶⁰Co to the UVOIR light curve of SN 1987A. We show that the V-band luminosity constitutes a roughly constant fraction of the bolometric luminosity between 900 and 1900 days, and we obtain an approximate bolometric light curve out to 4334 days by scaling the late time V-band data by a constant factor where no bolometric light curve data is available. Considering the five most relevant decay chains starting at ⁴⁴Ti, ⁵⁵Co, ⁵⁶Ni, ⁵⁷Ni, and ⁶⁰Co, we perform a least squares fit to the constructed composite bolometric light curve. For the nickel isotopes, we obtain best fit values of M(⁵⁶Ni) = (7.1 ± 0.3) × 10⁻² M_☉ and M(⁵⁷Ni) = (4.1 ± 1.8) × 10⁻³ M_☉. Our best fit ⁴⁴Ti mass is M(⁴⁴Ti) = (0.55 ± 0.17) × 10⁻⁴ M_☉, which is in disagreement with the much higher (3.1 ± 0.8) × 10⁻⁴ M_☉ recently derived from INTEGRAL observations. The associated uncertainties far exceed the best fit values for ⁵⁵Co and ⁶⁰Co and, as a result, we only give upper limits on the production masses of M(⁵⁵Co) < 7.2 × 10⁻³ M_☉ and M(⁶⁰Co) < 1.7 × 10⁻⁴ M_☉. Furthermore, we find that the leptonic channels in the decay of ⁵⁷Co (internal conversion and Auger electrons) are a significant contribution and constitute up to 15.5% of the total luminosity. Consideration of the kinetic energy of these electrons is essential in lowering our best fit nickel isotope production ratio to [⁵⁷Ni/⁵⁶Ni] = 2.5 ± 1.1, which is still somewhat high but is in agreement with gamma-ray observations and model predictions.

For the first time, we reveal large amounts of cold molecular gas in a ram-pressure-stripped tail, out to a large "intracluster" distance from the galaxy. With the Actama Pathfinder EXperiment (APEX) telescope, we have detected ¹²CO(2–1) emission corresponding to more than 10⁹ M_☉ of H₂ in three Hα bright regions along the tail of the Norma cluster galaxy ESO 137-001, out to a projected distance of 40 kpc from the disk. ESO 137-001 has an 80 kpc long and bright X-ray tail associated with a shorter (40 kpc) and broader tail of numerous star forming H ii regions. The amount of ∼1.5 × 10⁸ M_☉ of H₂ found in the most distant region is similar to molecular masses of tidal dwarf galaxies, though the standard Galactic CO-to-H₂ factor could overestimate the H₂ content. Along the tail, we find the amount of molecular gas to drop, while masses of the X-ray-emitting and diffuse ionized components stay roughly constant. Moreover, the amounts of hot and cold gas are large and similar, and together nearly account for the missing gas from the disk. We find a very low SFE (τ_dep > 10¹⁰ yr) in the stripped gas in ESO 137-001 and suggest that this is due to a low average gas density in the tail, or turbulent heating of the interstellar medium that is induced by a ram pressure shock. The unprecedented bulk of observed H₂ in the ESO 137-001 tail suggests that some stripped gas may survive ram pressure stripping in the molecular phase.

Sunspots are a canonical marker of the Sun's internal magnetic field which flips polarity every ∼22 yr. The principal variation of sunspots, an ∼11 yr variation, modulates the amount of the magnetic field that pierces the solar surface and drives significant variations in our star's radiative, particulate, and eruptive output over that period. This paper presents observations from the Solar and Heliospheric Observatory and Solar Dynamics Observatory indicating that the 11 yr sunspot variation is intrinsically tied to the spatio-temporal overlap of the activity bands belonging to the 22 yr magnetic activity cycle. Using a systematic analysis of ubiquitous coronal brightpoints and the magnetic scale on which they appear to form, we show that the landmarks of sunspot cycle 23 can be explained by considering the evolution and interaction of the overlapping activity bands of the longer-scale variability.

We study the Hα and Ca ii 8542 Å line spectra of four typical Ellerman bombs (EBs) in the active region NOAA 11765 on 2013 June 6, observed with the Fast Imaging Solar Spectrograph installed at the 1.6 m New Solar Telescope at Big Bear Solar Observatory. Considering that EBs may occur in a restricted region in the lower atmosphere, and that their spectral lines show particular features, we propose a two-cloud model to fit the observed line profiles. The lower cloud can account for the wing emission, and the upper cloud is mainly responsible for the absorption at line center. After choosing carefully the free parameters, we get satisfactory fitting results. As expected, the lower cloud shows an increase of the source function, corresponding to a temperature increase of 400–1000 K in EBs relative to the quiet Sun. This is consistent with previous results deduced from semi-empirical models and confirms that local heating occurs in the lower atmosphere during the appearance of EBs. We also find that the optical depths can increase to some extent in both the lower and upper clouds, which may result from either direct heating in the lower cloud, or illumination by an enhanced radiation on the upper cloud. The velocities derived from this method, however, are different from those obtained using the traditional bisector method, implying that one should be cautious when interpreting this parameter. The two-cloud model can thus be used as an efficient method to deduce the basic physical parameters of EBs.

We present results from the first Suzaku observation of the high-mass X-ray binary 4U 1538−522. The broadband spectral coverage of Suzaku allows for a detailed spectral analysis, characterizing the cyclotron resonance scattering feature at 23.0 ± 0.4  keV and the iron Kα line at 6.426 ± 0.008 keV, as well as placing limits on the strengths of the iron Kβ line and the iron K edge. We track the evolution of the spectral parameters both in time and in luminosity, notably finding a significant positive correlation between cyclotron line energy and luminosity. A dip and spike in the light curve is shown to be associated with an order-of-magnitude increase in column density along the line of sight, as well as significant variation in the underlying continuum, implying the accretion of a overdense region of a clumpy stellar wind. We also present a phase-resolved analysis, with most spectral parameters of interest showing significant variation with phase. Notably, both the cyclotron line energy and the iron Kα line intensity vary significantly with phase, with the iron line intensity significantly out of phase with the pulse profile. We discuss the implications of these findings in the context of recent work in the areas of accretion column physics and cyclotron resonance scattering feature formation.

We performed a systematic spectroscopic observation of a protocluster at z = 6.01 in the Subaru Deep Field. We took spectroscopy for all 53 i' dropout galaxies down to z' = 27.09 mag in/around the protocluster region. From these observations, we confirmed that 28 galaxies are at z ∼ 6, 10 of which are clustered in a narrow redshift range of Δz < 0.06. To trace the evolution of this primordial structure, we applied the same i' dropout selection and the same overdensity measurements used in the observations to a semi-analytic model built upon the Millennium Simulation. We obtain a relation between the significance of overdensities observed at z ∼ 6 and the predicted dark matter halo mass at z = 0. This protocluster with 6σ overdensity is expected to grow into a galaxy cluster with a mass of ∼5 × 10¹⁴ M_☉ at z = 0. Ten galaxies within 10 comoving Mpc of the overdense region can, with more than an 80% probability, merge into a single dark matter halo by z = 0. No significant differences appeared in UV and Lyα luminosities between the protocluster and field galaxies, suggesting that this protocluster is still in the early phase of cluster formation before the onset of any obvious environmental effects. However, further observations are required to study other properties, such as stellar mass, dust, and age. We do find that galaxies tend to be in close pairs in this protocluster. These pair-like subgroups will coalesce into a single halo and grow into a more massive structure. We may witness an onset of cluster formation at z ∼ 6 toward a cluster as seen in local universe.

The Cl 1604 supercluster at z ∼ 0.9 is one of the most extensively studied high-redshift large-scale structures, with more than 500 spectroscopically confirmed members. It consists of eight clusters and groups, with members numbering from a dozen to nearly a hundred, providing a broad range of environments for investigating the large-scale environmental effects on galaxy evolution. Here we examine the properties of 48 post-starburst galaxies in Cl 1604, comparing them to other galaxy populations in the same supercluster. Incorporating photometry from ground-based optical and near-infrared imaging, along with Spitzer mid-infrared observations, we derive stellar masses for all Cl 1604 members. The colors and stellar masses of the K+A galaxies support the idea that they are progenitors of red sequence galaxies. Their morphologies, residual star formation rates, and spatial distributions suggest that galaxy mergers may be the principal mechanism producing post-starburst galaxies. Interaction between galaxies and the dense intracluster medium (ICM) is also effective, but only in the cores of dynamically evolved clusters. The prevalence of post-starburst galaxies in clusters correlates with the dynamical state of the host cluster, as both galaxy mergers and the dense ICM produce post-starburst galaxies. We also investigate the incompleteness and contamination of K+A samples selected by means of Hδ and [O ii] equivalent widths. K+A samples may be up to ∼50% incomplete due to the presence of LINERs/Seyferts, and up to ∼30% of K+A galaxies could have substantial star formation activity.

Gas-giant planets emit a large fraction of their light in the mid-infrared (≳3 μm), where photometry and spectroscopy are critical to our understanding of the bulk properties of extrasolar planets. Of particular importance are the L- and M-band atmospheric windows (3–5 μm), which are the longest wavelengths currently accessible to ground-based, high-contrast imagers. We present binocular LBT adaptive optics (AO) images of the HR 8799 planetary system in six narrow-band filters from 3 to 4 μm, and a Magellan AO image of the 2M1207 planetary system in a broader 3.3 μm band. These systems encompass the five known exoplanets with luminosities consistent with L → T transition brown dwarfs. Our results show that the exoplanets are brighter and have shallower spectral slopes than equivalent temperature brown dwarfs in a wavelength range that contains the methane fundamental absorption feature (spanned by the narrow-band filters and encompassed by the broader 3.3 μm filter). For 2M1207 b, we find that thick clouds and non-equilibrium chemistry caused by vertical mixing can explain the object's appearance. For the HR 8799 planets, we present new models that suggest the atmospheres must have patchy clouds, along with non-equilibrium chemistry. Together, the presence of a heterogeneous surface and vertical mixing presents a picture of dynamic planetary atmospheres in which both horizontal and vertical motions influence the chemical and condensate profiles.

We calculate the mass loss driven by magnetohydrodynamic (MHD) waves from hot Jupiters by using MHD simulations in one-dimensional flux tubes. If a gaseous planet has a magnetic field, MHD waves are excited by turbulence at the surface, dissipate in the upper atmosphere, and drive gas outflows. Our calculation shows that mass-loss rates are comparable to the observed mass-loss rates of hot Jupiters; therefore, it is suggested that gas flow driven by MHD waves can play an important role in the mass loss from gaseous planets. The mass-loss rate varies dramatically with the radius and mass of a planet: a gaseous planet with a small mass but an inflated radius produces a very large mass-loss rate. We also derive an analytical expression for the dependence of mass-loss rate on planet radius and mass that is in good agreement with the numerical calculation. The mass-loss rate also depends on the amplitude of the velocity dispersion at the surface of a planet. Thus, we expect to infer the condition of the surface and the internal structure of a gaseous planet from future observations of mass-loss rate from various exoplanets.

We report the detection of a dispersed fast radio burst (FRB) in archival intermediate-latitude Parkes Radio Telescope data. The burst appears to be of the same physical origin as the four purported extragalactic FRBs reported by Thornton et al. This burst's arrival time precedes the Thornton et al. bursts by 10 years. We consider that this survey, and many other archival low-latitude (|gb| < 30°) pulsar surveys, have been searched for FRBs but produced fewer detections than the comparatively brief Thornton et al. search. Such a rate dependence on Galactic position could provide critical supporting evidence for an extragalactic origin for FRBs. To test this, we form an analytic expression to account for Galactic position and survey setup in FRB rate predictions. Employing a sky temperature, scattering, and dispersion model of the Milky Way, we compute the expected number of FRBs if they are isotropically distributed on the sky with respect to the Galactic position (i.e., local), and if they are of extragalactic origin. We demonstrate that the relative detection rates reject a local origin with a confidence of 99.96% (∼3.6σ). The extragalactic predictions provide a better agreement; however, there are still strong discrepancies with the low-latitude detection rate at a confidence of 99.69% (∼2.9σ). However, for the extragalactic population, the differences in predicted versus detected population may be accounted for by a number of factors, which we discuss.

We explore XMM-Newton observations of the eclipsing super-soft X-ray source CAL 87 in order to map the accretion structures of the system. Indirect imaging techniques were applied in X-ray light curves to provide eclipse maps. The surface brightness distribution exhibits an extended and symmetric emission, and a feature is revealed from the hardest X-rays that is likely due to a bright spot. A rate of $\dot{P} = (+6\pm 2) \times 10^{-10}$ for changes in the orbital period of the system was derived from the eclipses. There is no significant variation of the emission lines even during eclipses, arguing that the lines are formed in an extended region. The continuum emission dominates the decrease in flux that is observed during eclipses. The O viii Lyα line reveals a broadening velocity that is estimated to be 365$^{+65}_{-69}$ km s⁻¹ (at 1σ), marginal evidence for asymmetry in its profile, and sometimes shows evidence of double-peaked emission. Together, the results support that the wind-driven mass transfer scenario is running in CAL 87.

A long-standing question is whether radiative cooling can lead to local condensation of cold gas in the hot atmospheres of galaxies and galaxy clusters. We address this problem by studying the nature of local instabilities in rotating, stratified, weakly magnetized, optically thin plasmas in the presence of radiative cooling and anisotropic thermal conduction. For both axisymmetric and nonaxisymmetric linear perturbations, we provide general equations which can be applied locally to specific systems to establish whether they are unstable and, in case of instability, to determine the kind of evolution (monotonically growing or overstable) and the growth rates of the unstable modes. We present results for models of rotating plasmas representative of Milky-Way-like galaxy coronae and cool-cores of galaxy clusters. We show that the unstable modes arise from a combination of thermal, magnetothermal, magnetorotational, and heat-flux-driven buoyancy instabilities. Local condensation of cold clouds tends to be hampered in cluster cool cores, while it is possible under certain conditions in rotating galactic coronae. If the magnetic field is sufficiently weak, then the magnetorotational instability is dominant even in these pressure-supported systems.

The analysis of Ca-K line spectra as a function of solar latitude and integrated over longitudes taken during the period of 1989–2011 indicates that the K₁ width attains maximum amplitude at various latitude belts at different phases of the solar cycle. The FWHM of the K₁ distribution at different latitudes shows that its width varies by about 30% for the equatorial belt (<30°) and 11% for the polar region (>70°) latitudes. Interestingly, the K₁ width varies by ∼6% only in the 40–60° latitude belts during the solar cycle. The analysis of cross-correlation coefficients of the K₁ width between the 35° latitude and other latitude belts as a function of phase differences indicates that the activity representing toroidal field shifted at a uniform rate of about 5.1 m s⁻¹ in the northern hemisphere from mid-latitudes toward the equator. In the southern hemisphere, activity shifted at a faster rate, ∼14 m s⁻¹, in the beginning of the cycle and the speed decreased with time, yielding an average speed of 7.5 m s⁻¹ toward the equator. The shift of activity in the higher latitude belts showed complex behavior, indicating poleward and equatorward migration. These findings, especially the fewer variations in mid-latitude belts as compared to polar regions, asymmetry in the speed of the shift in the activity in both hemispheres, and complex variation in the direction of the shift in the activity representing poloidal fields in mid-latitude belts, will have an important implication on the modeling of the solar dynamo.

The Fe xiii 10747 and 10798 Å lines observed in the solar corona are sensitive to the coronal magnetic field in such a way that, in principle, the full vector field at a point on the line of sight can be inferred from their combined polarization signals. This paper presents analytical inversion formulae for the field parameters and analyzes the uncertainty of magnetic field measurements made from such observations, assuming emission dominated by a single region along the line of sight. We consider the case of the current Coronal Multi-channel Polarimeter (CoMP) instrument as well as the future Coronal Solar Magnetism Observatory (COSMO) and Advanced Technology Solar Telescope (ATST) instruments. Uncertainties are estimated with a direct analytic inverse and with a Markov Chain Monte Carlo algorithm. We find that (in effect) two components of the vector field can be recovered with CoMP, and well recovered with COSMO or ATST, but that the third component can only be recovered when the solar magnetic field is strong and optimally oriented.

We describe the methodology to include nonlinear evolution, including tidal effects, in the computation of subhalo distribution properties in both cold (CDM) and warm (WDM) dark matter universes. Using semi-analytic modeling, we include effects from dynamical friction, tidal stripping, and tidal heating, allowing us to dynamically evolve the subhalo distribution. We calibrate our nonlinear evolution scheme to the CDM subhalo mass function in the Aquarius N-body simulation, producing a subhalo mass function within the range of simulations. We find tidal effects to be the dominant mechanism of nonlinear evolution in the subhalo population. Finally, we compute the subhalo mass function for m_χ = 1.5 keV WDM including the effects of nonlinear evolution, and compare radial number densities and mass density profiles of subhalos in CDM and WDM models. We show that all three signatures differ between the two dark matter models, suggesting that probes of substructure may be able to differentiate between them.

Cosmological constraints from X-ray and microwave observations of galaxy clusters are subjected to systematic uncertainties. Non-thermal pressure support due to internal gas motions in galaxy clusters is one of the major sources of astrophysical uncertainties. Using a mass-limited sample of galaxy clusters from a high-resolution hydrodynamical cosmological simulation, we characterize the non-thermal pressure fraction profile and study its dependence on redshift, mass, and mass accretion rate. We find that the non-thermal pressure fraction profile is universal across redshift when galaxy cluster radii are defined with respect to the mean matter density of the universe instead of the commonly used critical density. We also find that the non-thermal pressure is predominantly radial, and the gas velocity anisotropy profile exhibits strong universality when galaxy cluster radii are defined with respect to the mean matter density of the universe. However, we find that the non-thermal pressure fraction is strongly dependent on the mass accretion rate of the galaxy cluster. We provide fitting formulae for the universal non-thermal pressure fraction and velocity anisotropy profiles of gas in galaxy clusters, which should be useful in modeling astrophysical uncertainties pertinent to using galaxy clusters as cosmological probes.

We motivate the Ĝ infrared search for extraterrestrial civilizations with large energy supplies. We discuss some philosophical difficulties of the search for extraterrestrial intelligence (SETI), and how communication SETI circumvents them. We review "Dysonian SETI," the search for artifacts of alien civilizations, and find that it is highly complementary to traditional communication SETI; the two together might succeed where either one alone has not. We discuss the argument of Hart that spacefaring life in the Milky Way should be either galaxy-spanning or non-existent, and examine a portion of his argument that we call the "monocultural fallacy." We discuss some rebuttals to Hart that invoke sustainability and predict long Galaxy colonization timescales. We find that the maximum Galaxy colonization timescale is actually much shorter than previous work has found (<10⁹ yr), and that many "sustainability" counter-arguments to Hart's thesis suffer from the monocultural fallacy. We extend Hart's argument to alien energy supplies and argue that detectably large energy supplies can plausibly be expected to exist because life has the potential for exponential growth until checked by resources or other limitations, and intelligence implies the ability to overcome such limitations. As such, if Hart's thesis is correct, then searches for large alien civilizations in other galaxies may be fruitful; if it is incorrect, then searches for civilizations within the Milky Way are more likely to succeed than Hart argued. We review some past Dysonian SETI efforts and discuss the promise of new mid-infrared surveys, such as that of WISE.

We describe the framework and strategy of the Ĝ infrared search for extraterrestrial civilizations with large energy supplies, which will use the wide-field infrared surveys of WISE and Spitzer to search for these civilizations' waste heat. We develop a formalism for translating mid-infrared photometry into quantitative upper limits on extraterrestrial energy supplies. We discuss the likely sources of false positives, how dust can and will contaminate our search, and prospects for distinguishing dust from alien waste heat. We argue that galaxy-spanning civilizations may be easier to distinguish from natural sources than circumstellar civilizations (i.e., Dyson spheres), although GAIA will significantly improve our capability to identify the latter. We present a zeroth order null result of our search based on the WISE all-sky catalog: we show, for the first time, that Kardashev Type iii civilizations (as Kardashev originally defined them) are very rare in the local universe. More sophisticated searches can extend our methodology to smaller waste heat luminosities, and potentially entirely rule out (or detect) both Kardashev Type iii civilizations and new physics that allows for unlimited "free" energy generation.

Massive stars that end their lives with helium cores in the range of 35–65 M_☉ are known to produce repeated thermonuclear outbursts due to a recurring pair-instability. In some of these events, solar masses of material are ejected in repeated outbursts of several × 10⁵⁰ erg each. Collisions between these shells can sometimes produce very luminous transients that are visible from the edge of the observable universe. Previous one-dimensional (1D) studies of these events produce thin, high-density shells as one ejection plows into another. Here, in the first multi-dimensional simulations of these collisions, we show that the development of a Rayleigh–Taylor instability truncates the growth of the high-density spike and drives mixing between the shells. The progenitor is a 110 M_☉ solar-metallicity star that was shown in earlier work to produce a superluminous supernova. The light curve of this more realistic model has a peak luminosity and duration that are similar to those of 1D models but a structure that is smoother.

Type Iax supernovae (SNe Iax) are thermonuclear explosions that are related to SNe Ia, but are physically distinct. The most important differences are that SNe Iax have significantly lower luminosity (1%–50% that of typical SNe Ia), lower ejecta mass (∼0.1–0.5 M_☉), and may leave a bound remnant. The most extreme SN Iax is SN 2008ha, which peaked at M_V = −14.2 mag, about 5 mag below that of typical SNe Ia. Here, we present Hubble Space Telescope (HST) images of UGC 12682, the host galaxy of SN 2008ha, taken 4.1 yr after the peak brightness of SN 2008ha. In these deep, high-resolution images, we detect a source coincident (0.86 HST pixels; 0043; 1.1σ) with the position of SN 2008ha with M_F814W = −5.4 mag. We determine that this source is unlikely to be a chance coincidence, but that scenario cannot be completely ruled out. If this source is directly related to SN 2008ha, it is either the luminous bound remnant of the progenitor white dwarf (WD) or its companion star. The source is consistent with being an evolved >3 M_☉ initial mass star, and is significantly redder than the SN Iax 2012Z progenitor system, the first detected progenitor system for a thermonuclear SN. If this source is the companion star for SN 2008ha, there is a diversity in SN Iax progenitor systems, perhaps related to the diversity in SN Iax explosions. If the source is the bound remnant of the WD, it must have expanded significantly. Regardless of the nature of this source, we constrain the progenitor system of SN 2008ha to have an age of <80 Myr.

NASA's Wide-field Infrared Survey Explorer (WISE) spacecraft has been brought out of hibernation and has resumed surveying the sky at 3.4 and 4.6 μm. The scientific objectives of the NEOWISE reactivation mission are to detect, track, and characterize near-Earth asteroids and comets. The search for minor planets resumed on 2013 December 23, and the first new near-Earth object (NEO) was discovered 6 days later. As an infrared survey, NEOWISE detects asteroids based on their thermal emission and is equally sensitive to high and low albedo objects; consequently, NEOWISE-discovered NEOs tend to be large and dark. Over the course of its three-year mission, NEOWISE will determine radiometrically derived diameters and albedos for ∼2000 NEOs and tens of thousands of Main Belt asteroids. The 32 months of hibernation have had no significant effect on the mission's performance. Image quality, sensitivity, photometric and astrometric accuracy, completeness, and the rate of minor planet detections are all essentially unchanged from the prime mission's post-cryogenic phase.

Re-brightening bumps are frequently observed in gamma-ray burst afterglows. Many scenarios have been proposed to interpret the origin of these bumps, of which a blast wave encountering a density-jump in the circumburst environment has been questioned by recent works. We develop a set of differential equations to calculate the relativistic outflow encountering the density-jump by extending the work of Huang et al. This approach is a semi-analytic method and is very convenient. Our results show that late high-amplitude bumps cannot be produced under common conditions, rather only a short plateau may emerge even when the encounter occurs at an early time (<10⁴ s). In general, our results disfavor the density-jump origin for those observed bumps, which is consistent with the conclusion drawn from full hydrodynamics studies. The bumps thus should be caused by other scenarios.

We perform cosmological hydrodynamics simulations with non-equilibrium primordial chemistry to obtain 59 minihalos that host first stars. The obtained minihalos are used as the initial conditions of local three-dimensional radiation hydrodynamics simulations to investigate the formation of the first stars. We find that two-thirds of the minihalos host multiple stars, while the other third has single stars. The mass of the stars found in our simulations are in the range of 1 M_☉ ≲ M ≲ 300 M_☉, peaking at several× 10 M_☉. Most of the very massive stars of ≳ 140 M_☉ are born as single stars, although not all of the single stars are very massive. We also find a few stars of ≲ 1 M_☉ that are kicked by the gravitational three body interactions to the position distant from the center of mass. The frequency that a star forming minihalo contains a binary system is ∼50%. We also investigate the abundance pattern of the stellar remnants by summing up the contributions from the first stars in the simulations. Consequently, the pattern is compatible with that of the low metallicity damped Lyα systems or the extremely metal-poor (EMP) stars, if the mass spectrum obtained in our experiment is shifted to the low mass side by 0.2 dex. If we consider the case that an EMP star is born in the remnant of the individual minihalo without mixing with others, the chemical signature of the pair instability supernova is more prominent, because most of them are born as single stars.

In a sample of 60 quasars selected from the Sloan Digital Sky Survey with at least six-epoch spectroscopy, we investigate the variability of emission lines and continuum luminosity at various aspects. A strong anti-correlation between the variability and continuum luminosity at 2500 Å is found for the sample, which is consistent with previous works. In individual sources, we find that half of the sample objects follow the trend of being bluer when brighter, while the remaining half follow the redder-when-brighter (RWB) trend. Although the mechanism for RWB is unclear, the effects of host galaxy contribution due to seeing variations cannot be completely ruled out. As expected from the photoionization model, the positive correlations between the broad emission line and continuum luminosity are found in most individual sources, as well as for the whole sample. We confirm the Baldwin effect in most individual objects and the whole sample, while a negative Baldwin effect is also found in several quasars, which can be at least partly (if not all) due to the host galaxy contamination. We find positive correlations between the broad emission line luminosity and line width in most individual quasars, as well as the whole sample, implying a line base that is more variable than the line core.

We present ALMA observations of the [C ii] line and far-infrared (FIR) continuum of a normally star-forming galaxy in the reionization epoch, the z = 6.96 Lyα emitter (LAE) IOK-1. Probing to sensitivities of σ_line = 240 μJy beam⁻¹ (40 km s⁻¹ channel) and σ_cont = 21 μJy beam⁻¹, we found the galaxy undetected in both [C ii] and continuum. Comparison of ultraviolet (UV)–FIR spectral energy distribution (SED) of IOK-1, including our ALMA limit, with those of several types of local galaxies (including the effects of the cosmic microwave background, CMB, on the FIR continuum) suggests that IOK-1 is similar to local dwarf/irregular galaxies in SED shape rather than highly dusty/obscured galaxies. Moreover, our 3σ FIR continuum limit, corrected for CMB effects, implies intrinsic dust mass M_dust < 6.4 × 10⁷ M_☉, FIR luminosity L_FIR < 3.7 × 10¹⁰ L_☉ (42.5–122.5 μm), total IR luminosity L_IR < 5.7 × 10¹⁰ L_☉ (8–1000 μm), and dust-obscured star formation rate (SFR) < 10 M_☉ yr⁻¹, if we assume that IOK-1 has a dust temperature and emissivity index typical of local dwarf galaxies. This SFR is 2.4 times lower than one estimated from the UV continuum, suggesting that <29% of the star formation is obscured by dust. Meanwhile, our 3σ [C ii] flux limit translates into [C ii] luminosity, L_[C ii] < 3.4 × 10⁷ L_☉. Locations of IOK-1 and previously observed LAEs on the L_[C ii] versus SFR and L_[C ii]/L_FIR versus L_FIR diagrams imply that LAEs in the reionization epoch have significantly lower gas and dust enrichment than AGN-powered systems and starbursts at similar/lower redshifts, as well as local star-forming galaxies.

SPTpol, POLARBEAR, and BICEP2 have recently measured the cosmic microwave background (CMB) B-mode polarization in various sky regions of several tens of square degrees and obtained BB power spectra in the multipole range 20–3000, detecting the components due to gravitational lensing and to inflationary gravitational waves. We analyze jointly the results of these three experiments and propose modifications to their analyses of the spectra to include in the model, in addition to the gravitational lensing and the inflationary gravitational wave components, and also the effects induced by the cosmic polarization rotation (CPR), if it exists within current upper limits. Although in principle our analysis would also lead to new constraints on CPR, in practice these can only be given on its fluctuations 〈δα²〉, since constraints on its mean angle are inhibited by the derotation which is applied by current CMB polarization experiments, in order to cope with the insufficient calibration of the polarization angle. The combined data fits from all three experiments (with 29% CPR–SPTpol correlation, depending on the theoretical model) gives the constraint 〈δα²〉^1/2 < 27.3 mrad (156), with r = 0.194 ± 0.033. These results show that the present data are consistent with no CPR detection and the constraint on CPR fluctuation is about 15. This method of constraining the CPR is new, is complementary to previous tests, which use the radio and optical/UV polarization of radio galaxies and the CMB E-mode polarization, and adds a new constraint for the sky areas observed by SPTpol, POLARBEAR, and BICEP2.

An analytical comparative study of a two-fluid and a gyrokinetic model of kinetic Alfvén waves (KAWs) is presented for various solar and space plasma environments. Based on the linear KAW dispersion relation for gyrokinetics (Howes et al. 2006), the wave group velocity and electromagnetic polarizations are obtained analytically. Then the gyrokinetic wave properties are compared with those of the two-fluid model. The results show that both models agree well with each other not only in the long wavelength regime (≫ the ion gyroradius ρ_i) for all cases considered, but also in wavelengths ∼ρ_i and ≪ρ_i (still much larger than the electron gyroscale) for a moderate or low (≲ 1) and a high (≫1) ion/electron temperature ratio T_0i/T_0e, respectively. However, the fluid model calculations deviate strongly from the gyrokinetic model at scales <ρ_i for a relatively low T_0i/T_0e due to the electron gyroradius effect. Meanwhile, the plasma β_i can make the gyrokinetic dispersion relation of KAWs become complex and sometimes have an oscillation-like structure. With the inherent simplicity of the fluid theory, these results may improve our understanding of the applicability of the two-fluid model, and may have important implications for computer simulation studies of KAWs in the solar and space plasma surroundings.

Based on high-resolution optical spectra obtained with ESPaDOnS at Canada–France–Hawaii Telescope, we determine fundamental parameters (T_eff, R, L_bol, log g,  and metallicity) for 59 candidate members of nearby young kinematic groups. The candidates were identified through the BANYAN Bayesian inference method of Malo et al., which takes into account the position, proper motion, magnitude, color, radial velocity, and parallax (when available) to establish a membership probability. The derived parameters are compared to Dartmouth magnetic evolutionary models and field stars with the goal of constraining the age of our candidates. We find that, in general, low-mass stars in our sample are more luminous and have inflated radii compared to older stars, a trend expected for pre-main-sequence stars. The Dartmouth magnetic evolutionary models show a good fit to observations of field K and M stars, assuming a magnetic field strength of a few kG, as typically observed for cool stars. Using the low-mass members of the β Pictoris moving group, we have re-examined the age inconsistency problem between lithium depletion age and isochronal age (Hertzspring–Russell diagram). We find that the inclusion of the magnetic field in evolutionary models increases the isochronal age estimates for the K5V–M5V stars. Using these models and field strengths, we derive an average isochronal age between 15 and 28 Myr and we confirm a clear lithium depletion boundary from which an age of 26 ± 3 Myr is derived, consistent with previous age estimates based on this method.

We report deep EVN and eMERLIN observations of the Type Ia SN 2014J in the nearby galaxy M82. Our observations represent, together with JVLA observations of SNe 2011fe and 2014J, the most sensitive radio studies of Type Ia SNe ever. By combining data and a proper modeling of the radio emission, we constrain the mass-loss rate from the progenitor system of SN 2014J to $\dot{M} \lesssim 7.0\times 10^{-10} {\,{M_{\odot } \,\rm yr^{-1}}}$ (for a wind speed of 100 km s⁻¹). If the medium around the supernova is uniform, then n_ISM ≲ 1.3  cm⁻³, which is the most stringent limit for the (uniform) density around a Type Ia SN. Our deep upper limits favor a double-degenerate (DD) scenario—involving two WD stars—for the progenitor system of SN 2014J, as such systems have less circumstellar gas than our upper limits. By contrast, most single-degenerate (SD) scenarios, i.e., the wide family of progenitor systems where a red giant, main-sequence, or sub-giant star donates mass to an exploding WD, are ruled out by our observations. (While completing our work, we noticed that a paper by Margutti et al. was submitted to The Astrophysical Journal. From a non-detection of X-ray emission from SN 2014J, the authors obtain limits of $\dot{M} \lesssim 1.2 \times 10^{-9}$ M_☉ yr⁻¹ (for a wind speed of 100 km s⁻¹) and n_ISM ≲ 3.5  cm⁻³, for the ρ∝r⁻² wind and constant density cases, respectively. As these limits are less constraining than ours, the findings by Margutti et al. do not alter our conclusions. The X-ray results are, however, important to rule out free–free and synchrotron self-absorption as a reason for the radio non-detections.) Our estimates on the limits on the gas density surrounding SN2011fe, using the flux density limits from Chomiuk et al., agree well with their results. Although we discuss the possibilities of an SD scenario passing observational tests, as well as uncertainties in the modeling of the radio emission, the evidence from SNe 2011fe and 2014J points in the direction of a DD scenario for both.

We carry out high-speed photometry on 20 of the shortest-period, detached white dwarf binaries known and discover systems with eclipses, ellipsoidal variations (due to tidal deformations of the visible white dwarf), and Doppler beaming. All of the binaries contain low-mass white dwarfs with orbital periods of less than four hr. Our observations identify the first eight tidally distorted white dwarfs, four of which are reported for the first time here. We use these observations to place empirical constraints on the mass–radius relationship for extremely low-mass (⩽0.30 M_☉) white dwarfs. We also detect Doppler beaming in several of these binaries, which confirms their high-amplitude radial-velocity variability. All of these systems are strong sources of gravitational radiation, and long-term monitoring of those that display ellipsoidal variations can be used to detect spin-up of the tidal bulge due to orbital decay.

We study the quasi-static evolution of coronal magnetic fields constructed from the non-linear force-free field (NLFFF) approximation aiming to understand the relation between the magnetic field topology and ribbon emission during an X1.5 flare in active region (AR) NOAA 11166. The flare with a quasi-elliptical and two remote ribbons occurred on 2011 March 9 at 23:13 UT over a positive flux region surrounded by negative flux at the center of the bipolar AR. Our analysis of the coronal magnetic structure with potential and NLFFF solutions unveiled the existence of a single magnetic null point associated with a fan–spine topology and is co-spatial with the hard X-ray source. The footpoints of the fan separatrix surface agree with the inner edge of the quasi-elliptical ribbon and the outer spine is linked to one of the remote ribbons. During the evolution, the slow footpoint motions stressed the field lines along the polarity inversion line and caused electric current layers in the corona around the fan separatrix surface. These current layers trigger magnetic reconnection as a consequence of dissipating currents, which are visible as cusp-shaped structures at lower heights. The reconnection process reorganized the magnetic field topology whose signatures are observed at the separatrices/quasi-separatrix layer structure in both the photosphere and the corona during the pre-to-post flare evolution. In agreement with previous numerical studies, our results suggest that the line-tied footpoint motions perturb the fan-spine system and cause null point reconnection, which eventually causes the flare emission at the footpoints of the field lines.

The Solar Optical Telescope on board Hinode observed a sunspot (AR 11836) with two light bridges (LBs) on 2013 August 31. We analyzed a two-hour Ca ii H emission intensity data set and detected strong five-minute oscillation power on both LBs and in the inner penumbra. The time–distance plot reveals that the five-minute oscillation phase does not vary significantly along the thin bridge, indicating that the oscillations are likely to originate from underneath it. The slit taken along the central axis of the wide LB exhibits a standing wave feature. However, at the center of the wide bridge, the five-minute oscillation power is found to be stronger than at its sides. Moreover, the time–distance plot across the wide bridge exhibits a herringbone pattern that indicates a counter-stream of two running waves, which originated at the bridge's sides. Thus, the five-minute oscillations on the wide bridge also resemble the properties of running penumbral waves. The five-minute oscillations are suppressed in the umbra, while the three-minute oscillations occupy all three cores of the sunspot's umbra, separated by the LBs. The three-minute oscillations were found to be in phase at both sides of the LBs. This may indicate that either LBs do not affect umbral oscillations, or that umbral oscillations at different umbral cores share the same source. It also indicates that LBs are rather shallow objects situated in the upper part of the umbra. We found that umbral flashes (UFs) follow the life cycles of umbral oscillations with much larger amplitudes. They cannot propagate across LBs. UFs dominate the three-minute oscillation power within each core; however, they do not disrupt the phase of umbral oscillation.

In many gamma-ray bursts a distinct blackbody spectral component is present, which is attributed to the emission from the photosphere of a relativistically expanding plasma. The properties of this component (temperature and flux) can be linked to the properties of the outflow and have been presented in the case where there is no sub-photospheric dissipation and the photosphere is in coasting phase. First, we present the derivation of the properties of the outflow for finite winds, including when the photosphere is in the accelerating phase. Second, we study the effect of localized sub-photospheric dissipation on the estimation of the parameters. Finally, we apply our results to GRB 090902B. We find that during the first epoch of this burst the photosphere is most likely to be in the accelerating phase, leading to smaller values of the Lorentz factor than the ones previously estimated. For the second epoch, we find that the photosphere is likely to be in the coasting phase.

We present a study of the discrete clouds and filaments in the Magellanic Stream using a new high-resolution survey of neutral hydrogen (H i) conducted with the H75 array of the Australia Telescope Compact Array, complemented by single-dish data from the Parkes Galactic All-Sky Survey. From the individual and combined data sets, we have compiled a catalog of 251 clouds and listed their basic parameters, including a morphological description useful for identifying cloud interactions. We find an unexpectedly large number of head–tail clouds in the region. The implication for the formation mechanism and evolution is discussed. The filaments appear to originate entirely from the Small Magellanic Cloud and extend into the northern end of the Magellanic Bridge.

Numerical studies of primordial star formation suggest that the first stars in the universe may have been very massive. Stellar models indicate that non-rotating Population III stars with initial masses of 140–260 M_☉ die as highly energetic pair-instability supernovae. We present new two-dimensional simulations of primordial pair-instability supernovae done with the CASTRO code. Our simulations begin at earlier times than previous multidimensional models, at the onset of core contraction, to capture any dynamical instabilities that may be seeded by core contraction and explosive burning. Such instabilities could enhance explosive yields by mixing hot ash with fuel, thereby accelerating nuclear burning, and affect the spectra of the supernova by dredging up heavy elements from greater depths in the star at early times. Our grid of models includes both blue supergiants and red supergiants over the range in progenitor mass expected for these events. We find that fluid instabilities driven by oxygen and helium burning arise at the upper and lower boundaries of the oxygen shell ∼20–100 s after core bounce. Instabilities driven by burning freeze out after the SN shock exits the helium core. As the shock later propagates through the hydrogen envelope, a strong reverse shock forms that drives the growth of Rayleigh–Taylor instabilities. In red supergiant progenitors, the amplitudes of these instabilities are sufficient to mix the supernova ejecta.

We present optical spectroscopy of galaxies in clusters detected through the Sunyaev–Zel'dovich (SZ) effect with the South Pole Telescope (SPT). We report our own measurements of 61 spectroscopic cluster redshifts, and 48 velocity dispersions each calculated with more than 15 member galaxies. This catalog also includes 19 dispersions of SPT-observed clusters previously reported in the literature. The majority of the clusters in this paper are SPT-discovered; of these, most have been previously reported in other SPT cluster catalogs, and five are reported here as SPT discoveries for the first time. By performing a resampling analysis of galaxy velocities, we find that unbiased velocity dispersions can be obtained from a relatively small number of member galaxies (≲ 30), but with increased systematic scatter. We use this analysis to determine statistical confidence intervals that include the effect of membership selection. We fit scaling relations between the observed cluster velocity dispersions and mass estimates from SZ and X-ray observables. In both cases, the results are consistent with the scaling relation between velocity dispersion and mass expected from dark-matter simulations. We measure a ∼30% log-normal scatter in dispersion at fixed mass, and a ∼10% offset in the normalization of the dispersion–mass relation when compared to the expectation from simulations, which is within the expected level of systematic uncertainty.

The asteroid (31) Euphrosyne is the largest body of its namesake family, and it contains more than 99% of the family mass. Among large asteroid families, the Euphrosyne group is peculiar because of its quite steep size–frequency distribution (SFD), significantly depleted in large- and medium-sized asteroids (8 < D < 12 km). The current steep SFD of the Euphrosyne family has been suggested to be the result of a grazing impact in which only the farthest, smallest members failed to accrete. The Euphrosyne family is, however, also very peculiar because of its dynamics: near its center it is crossed by the ν₆ = g − g₆ linear secular resonance, and it hosts the largest population (140 bodies) of asteroids in ν₆ antialigned librating states (or Tina-like asteroids) in the main belt. In this work we investigated the orbital evolution of newly obtained members of the dynamical family, with an emphasis on its interaction with the ν₆ resonance. Because of its unique resonant configuration, large- and medium-sized asteroids tend to migrate away from the family orbital region faster than small-sized objects, which were ejected farther away from the family center. As a consequence, the SFD of the Euphrosyne family becomes steeper in time with a growing depletion in the number of the largest family members. We estimate that the current SFD could be attained from a typical, initial SFD on timescales of 500 Myr, consistent with estimates of the family age obtained with other independent methods.

We use hydrodynamic simulations to study the effect of spiral arms on the star formation rate (SFR) in nuclear rings of barred-spiral galaxies. We find that spiral arms can be an efficient means of gas transport from the outskirts to the central parts, provided that the arms are rotating slower than the bar. While the ring star formation in models with no arms or corotating arms is active only during around the bar growth phase, arm-driven gas accretion both significantly enhances and prolongs the ring star formation in models with slow-rotating arms. The arm-enhanced SFR is larger by a factor of ∼3–20 than in the no-arm model, with larger values corresponding to stronger and slower arms. Arm-induced mass inflows also make dust lanes stronger. Nuclear rings in slow-arm models are ∼45% larger than in the no-arm counterparts. Star clusters that form in a nuclear ring exhibit an age gradient in the azimuthal direction only when the SFR is small, whereas no notable age gradient is found in the radial direction for models with arm-induced star formation.

The search for diffuse non-thermal inverse Compton (IC) emission from galaxy clusters at hard X-ray energies has been undertaken with many instruments, with most detections being either of low significance or controversial. Because all prior telescopes sensitive at E > 10 keV do not focus light and have degree-scale fields of view, their backgrounds are both high and difficult to characterize. The associated uncertainties result in lower sensitivity to IC emission and a greater chance of false detection. In this work, we present 266 ks NuSTAR observations of the Bullet cluster, which is detected in the energy range 3–30 keV. NuSTAR's unprecedented hard X-ray focusing capability largely eliminates confusion between diffuse IC and point sources; however, at the highest energies, the background still dominates and must be well understood. To this end, we have developed a complete background model constructed of physically inspired components constrained by extragalactic survey field observations, the specific parameters of which are derived locally from data in non-source regions of target observations. Applying the background model to the Bullet cluster data, we find that the spectrum is well—but not perfectly—described as an isothermal plasma with kT = 14.2 ± 0.2 keV. To slightly improve the fit, a second temperature component is added, which appears to account for lower temperature emission from the cool core, pushing the primary component to kT ∼ 15.3 keV. We see no convincing need to invoke an IC component to describe the spectrum of the Bullet cluster, and instead argue that it is dominated at all energies by emission from purely thermal gas. The conservatively derived 90% upper limit on the IC flux of 1.1 × 10⁻¹² erg s⁻¹ cm⁻² (50–100 keV), implying a lower limit on B ≳ 0.2 μG, is barely consistent with detected fluxes previously reported. In addition to discussing the possible origin of this discrepancy, we remark on the potential implications of this analysis for the prospects for detecting IC in galaxy clusters in the future.

We present a detailed study of the complete evolution of a coronal mass ejection (CME). We have tracked the evolution of both the ejecta and its shock, and further fit the evolution of the fronts to a simple but physics-based analytical model. This study focuses on the CME initiated on the Sun on 2012 July 12 and arriving at the Earth on 2012 July 14. Shock and ejecta fronts were observed by white light images, as well as in situ by the Advanced Composition Explorer satellite. We find that the propagation of the two fronts is not completely dependent upon one another, but can each be modeled in the heliosphere with a drag model that assumes the dominant force of affecting CME evolution to be the aerodynamic drag force of the ambient solar wind. Results indicate that the CME ejecta front undergoes a more rapid deceleration than the shock front within 50 R_☉ and therefore the propagation of the two fronts is not completely coupled in the heliosphere. Using the graduated cylindrical shell model, as well as data from time-elongation stack plots and in situ signatures, we show that the drag model can accurately describe the behavior of each front, but is more effective with the ejecta. We also show that without the in situ data, based on measurements out to 80 R_☉ combined with the general values for drag model parameters, the arrival of both the shock and ejecta can be predicted within four hours of arrival.

Magnetic reconnection in the corona results in contracting flare loops, releasing energy into plasma heating and shocks. The hydrodynamic shocks produced in this manner drive thermal conduction fronts (TCFs) which transport energy into the chromosphere and drive upflows (evaporation) and downflows (condensation) in the cooler, denser footpoint plasma. Observations have revealed that certain properties of the transition point between evaporation and condensation (the "flow reversal point" or FRP), such as temperature and velocity-temperature derivative at the FRP, vary between different flares. These properties may provide a diagnostic tool to determine parameters of the coronal energy release mechanism and the loop atmosphere. In this study, we develop a one-dimensional hydrodynamical flare loop model with a simplified three-region atmosphere (chromosphere/transition region/corona), with TCFs initiated by shocks introduced in the corona. We investigate the effect of two different flare loop parameters (post-shock temperature and transition region temperature ratio) on the FRP properties. We find that both of the evaporation characteristics have scaling-law relationships to the varied flare parameters, and we report the scaling exponents for our model. This provides a means of using spectroscopic observations of the chromosphere as quantitative diagnostics of flare energy release in the corona.

The Faraday dispersion function (FDF), which can be derived from an observed polarization spectrum by Faraday rotation measure synthesis, is a profile of polarized emissions as a function of Faraday depth. We study intrinsic FDFs along sight lines through face-on Milky Way like galaxies by means of a sophisticated galactic model incorporating three-dimensional MHD turbulence, and investigate how much information the FDF intrinsically contains. Since the FDF reflects distributions of thermal and cosmic-ray electrons as well as magnetic fields, it has been expected that the FDF could be a new probe to examine internal structures of galaxies. We, however, find that an intrinsic FDF along a sight line through a galaxy is very complicated, depending significantly on actual configurations of turbulence. We perform 800 realizations of turbulence and find no universal shape of the FDF even if we fix the global parameters of the model. We calculate the probability distribution functions of the standard deviation, skewness, and kurtosis of FDFs and compare them for models with different global parameters. Our models predict that the presence of vertical magnetic fields and the large-scale height of cosmic-ray electrons tend to make the standard deviation relatively large. In contrast, the differences in skewness and kurtosis are relatively less significant.

M74 (NGC 628) is a famous face-on spiral galaxy, hosting three core-collapse supernovae (SNe): SN Ic 2002ap, SN II-P 2003gd, and SN II-P 2013ej. However, its distance is not well known. We present a distance estimation for this galaxy based on the Tip of the Red Giant Branch (TRGB) method. We obtain photometry of the resolved stars in the arm-free region of M74 from F555W and F814W images in the Hubble Space Telescope archive. The color–magnitude diagram of the resolved stars shows a dominant red giant branch (RGB) as well as blue main sequence stars, red helium burning stars, and asymptotic giant branch stars. The I-band luminosity function of the RGB stars shows the TRGB to be at I_TRGB = 26.13 ± 0.03 mag, and T_RGB = 25.97 ± 0.03. From this, we derive the distance modulus to M74 to be 30.04 ± 0.04 (random) ± 0.12 (systematic) (corresponding to a linear distance of 10.19 ± 0.14 ± 0.56 Mpc). With this distance estimate, we calibrate the standardized candle method for SNe II-P. From the absolute magnitudes of SN 2003gd, we derive a value of the Hubble constant, H₀ = 72 ± 6 (random) ± 7 (systematic) km s⁻¹ Mpc⁻¹. It is similar to recent estimates based on the luminosity calibration of Type Ia supernovae.

We report an observational estimate of the rate of stellar tidal disruption flares (TDFs) in inactive galaxies based on a successful search for these events among transients in galaxies using archival Sloan Digital Sky Survey (SDSS) multi-epoch imaging data (Stripe 82). This search yielded 186 nuclear flares in galaxies, 2 of which are excellent TDF candidates. Because of the systematic nature of the search, the very large number of galaxies, the long time of observation, and the fact that non-TDFs were excluded without resorting to assumptions about TDF characteristics, this study provides an unparalleled opportunity to measure the TDF rate. To compute the rate of optical stellar tidal disruption events, we simulate our entire pipeline to obtain the efficiency of detection. The rate depends on the light curves of TDFs, which are presently still poorly constrained. Using only the observed part of the SDSS light curves gives a model-independent upper limit to the optical TDF rate, $\dot{N}<2\times 10^{-4}\,{\rm yr}^{-1}\,{\rm galaxy}^{-1}$ (90% CL), under the assumption that the SDSS TDFs are representative examples. We develop three empirical models of the light curves based on the two SDSS light curves and two more recent and better-sampled Pan-STARRS TDF light curves, leading to our best estimate of the rate: $\dot{N}_{\rm TDF} = (1.5\hbox{--}2.0)_{-1.3}^{+2.7} \times 10^{-5} \,{\rm yr}^{-1}\, {\rm galaxy}^{-1}$. We explore the modeling uncertainties by considering two theoretically motivated light curve models, as well as two different relationships between black hole mass and galaxy luminosity, and two different treatments of the cutoff in the visibility of TDFs at large M_BH. From this we conclude that these sources of uncertainty are not significantly larger than the statistical ones. Our results are applicable for galaxies hosting black holes with mass in the range of a few 10⁶–10⁸ M_☉, and translates to a volumetric TDF rate of (4–8) × 10^{−8 ± 0.4} yr⁻¹ Mpc⁻³, with the statistical uncertainty in the exponent.

Quasars are variable on timescales from days to years in UV/optical and generally appear bluer while they brighten. The physics behind the variations in fluxes and colors remains unclear. Using Sloan Digital Sky Survey g- and r-band photometric monitoring data for quasars in Stripe 82, we find that although the flux variation amplitude increases with timescale, the color variability exhibits the opposite behavior. The color variability of quasars is prominent at timescales as short as ∼10 days, but gradually reduces toward timescales up to years. In other words, the variable emission at shorter timescales is bluer than that at longer timescales. This timescale dependence is clearly and consistently detected at all redshifts from z = 0 to 3.5; thus, it cannot be due to contamination to broadband photometry from emission lines that do not respond to fast continuum variations. The discovery directly rules out the possibility that simply attributes the color variability to contamination from a non-variable redder component such as the host galaxy. It cannot be interpreted as changes in global accretion rate either. The thermal accretion disk fluctuation model is favored in the sense that fluctuations in the inner, hotter region of the disk are responsible for short-term variations, while longer-term and stronger variations are expected from the larger and cooler disk region. An interesting implication is that one can use quasar variations at different timescales to probe disk emission at different radii.

High Reynolds-number turbulence driven by stellar convection in main-sequence stars generates stochastic gravitational radiation. We calculate the wave-strain power spectral density as a function of the zero-age main-sequence mass for an individual star and for an isotropic, universal stellar population described by the Salpeter initial mass function and redshift-dependent Hopkins–Beacom star formation rate. The spectrum is a broken power law, which peaks near the turnover frequency of the largest turbulent eddies. The signal from the Sun dominates the universal background. For the Sun, the far-zone power spectral density peaks at S(f_peak) = 5.2 × 10⁻⁵² Hz⁻¹ at frequency f_peak = 2.3 × 10⁻⁷ Hz. However, at low observing frequencies f < 3 × 10⁻⁴ Hz, the Earth lies inside the Sun's near zone and the signal is amplified to S_near(f_peak) = 4.1 × 10⁻²⁷ Hz⁻¹ because the wave strain scales more steeply with distance (∝d⁻⁵) in the near zone than in the far zone (∝d⁻¹). Hence the Solar signal may prove relevant for pulsar timing arrays. Other individual sources and the universal background fall well below the projected sensitivities of the Laser Interferometer Space Antenna and next-generation pulsar timing arrays. Stellar convection sets a fundamental noise floor for more sensitive stochastic gravitational-wave experiments in the more distant future.

We present new K-band spectroscopy of the UY Aur binary star system. Our data are the first to show H₂ emission in the spectrum of UY Aur A and the first to spectrally resolve the Brγ line in the spectrum of UY Aur B. We see an increase in the strength of the Brγ line in UY Aur A and a decrease in Brγ and H₂ line luminosity for UY Aur B compared to previous studies. Converting Brγ line luminosity to accretion rate, we infer that the accretion rate onto UY Aur A has increased by 2 × 10⁻⁹M_☉ yr⁻¹ per year since a rate of zero was observed in 1994. The Brγ line strength for UY Aur B has decreased by a factor of 0.54 since 1994, but the K-band flux has increased by 0.9 mag since 1998. The veiling of UY Aur B has also increased significantly. These data evince a much more luminous disk around UY Aur B. If the lower Brγ luminosity observed in the spectrum of UY Aur B indicates an intrinsically smaller accretion rate onto the star, then UY Aur A now accretes at a higher rate than UY Aur B. However, extinction at small radii or mass pile-up in the circumstellar disk could explain decreased Brγ emission around UY Aur B even when the disk luminosity implies an increased accretion rate. In addition to our scientific results for the UY Aur system, we discuss a dedicated pipeline we have developed for the reduction of echelle-mode data from the ARIES spectrograph.

Observations of novae at radio frequencies provide us with a measure of the total ejected mass, density profile, and kinetic energy of a nova eruption. The radio emission is typically well characterized by the free–free emission process. Most models to date have assumed spherical symmetry for the eruption, although for as long as there have been radio observations of these systems, it has been known that spherical eruptions are too simplistic a geometry. In this paper, we build bipolar models of the nova eruption, assuming the free–free process, and show the effects of varying different parameters on the radio light curves. The parameters considered include the ratio of the minor- to major-axis, the inclination angle, and shell thickness. We also show the uncertainty introduced when fitting spherical-model synthetic light curves to bipolar-model synthetic light curves. We find that the optically thick phase rises with the same power law (S_ν∝t²) for both the spherical and bipolar models. In the bipolar case, there is a "plateau" phase—depending on the thickness of the shell as well as the ratio of the minor- to major-axis—before the final decline, which follows the same power law (S_ν∝t⁻³) as in the spherical case. Finally, fitting spherical models to the bipolar-model synthetic light curves requires, in the worst-case scenario, doubling the ejected mass, more than halving the electron temperature, and reducing the shell thickness by nearly a factor of 10. This implies that in some systems we have been over-predicting the ejected masses and under-predicting the electron temperature of the ejecta.

N70 and N185 are two large (⩾100 pc in diameter) shell nebulae in the Large Magellanic Cloud (LMC). Their high [S ii]/Hα ratios rival those of supernova remnants (SNRs), but they are not confirmed as SNRs. To study their physical nature, we have obtained XMM-Newton X-ray observations and high-dispersion long-slit echelle spectroscopic observations of these two nebulae. The X-ray spectra of both nebulae can be well interpreted with an optically thin thermal (∼0.2 keV) plasma with the average LMC abundance in a collisional ionization equilibrium. N70 encompasses the OB association LH114. Although N70 has a modest expansion velocity and essentially thermal radio emission, its diffuse X-ray luminosity (∼6.1 × 10³⁵ erg s⁻¹) is higher than that from a quiescent superbubble with N70's density, size, and expansion velocity; thus, N70 is most likely a superbubble that is recently energized by an interior SNR. N185 does not contain any known OB association, and its X-ray luminosity is an order of magnitude lower than expected if it is a quiescent superbubble. N185 has nonthermal radio emission and has high-velocity material expanding at nearly 200 km s⁻¹, similar to many known SNRs in the LMC. Its X-ray luminosity (∼1.9 × 10³⁵ erg s⁻¹) is also consistent with that of an evolved SNR. We therefore suggest that N185 is energized by a recent supernova.

We study the dynamics of the giant elliptical galaxy M87 from the central to the outermost regions with the made-to-measure (M2M) method. We use a new catalog of 922 globular cluster line-of-sight velocities extending to a projected radius of 180 kpc (equivalent to 25 M87 effective radii), and SAURON integral field unit data within the central 2.4 kpc. There are 263 globular clusters, mainly located beyond 40 kpc, newly observed by the Next Generation Virgo Survey. For the M2M modeling, the gravitational potential is taken as a combination of a luminous matter potential with a constant stellar mass-to-light ratio and a dark matter potential modeled as a logarithmic potential. Our best-fit dynamical model returns a stellar mass-to-light ratio in the I band of M/L_I = 6.0 ± 0.3 $M_{\odot } \,L_{\odot }^{-1}$ with a dark matter potential scale velocity of 591 ± 50 km s⁻¹ and scale radius of 42 ± 10 kpc. We determine the total mass of M87 within 180 kpc to be (1.5  ±  0.2) × 10¹³ M_☉. The mass within 40 kpc is smaller than previous estimates determined using globular cluster kinematics that did not extend beyond ∼45 kpc. With our new globular cluster velocities at much larger radii, we see that globular clusters around 40 kpc show an anomalously large velocity dispersion which affected previous results. The mass we derive is in good agreement with that inferred from ROSAT X-ray observation out to 180 kpc. Within 30 kpc our mass is also consistent with that inferred from Chandra and XMM-Newton X-ray observations, while within 120 kpc it is about 20% smaller. The model velocity dispersion anisotropy β parameter for the globular clusters in M87 is small, varying from −0.2 at the center to 0.2 at ∼40 kpc, and gradually decreasing to zero at ∼120 kpc.

We report a detection of a faint near-Earth asteroid (NEA) using our synthetic tracking technique and the CHIMERA instrument on the Palomar 200 inch telescope. With an apparent magnitude of 23 (H = 29, assuming detection at 20 lunar distances), the asteroid was moving at 632 day⁻¹ and was detected at a signal-to-noise ratio (S/N) of 15 using 30 s of data taken at a 16.7 Hz frame rate. The detection was confirmed by a second observation 77 minutes later at the same S/N. Because of its high proper motion, the NEA moved 7 arcsec over the 30 s of observation. Synthetic tracking avoided image degradation due to trailing loss that affects conventional techniques relying on 30 s exposures; the trailing loss would have degraded the surface brightness of the NEA image on the CCD down to an approximate magnitude of 25 making the object undetectable. This detection was a result of our 12 hr blind search conducted on the Palomar 200 inch telescope over two nights, scanning twice over six (53 × 0046) fields. Detecting only one asteroid is consistent with Harris's estimates for the distribution of the asteroid population, which was used to predict a detection of 1.2 NEAs in the H-magnitude range 28–31 for the two nights. The experimental design, data analysis methods, and algorithms are presented. We also demonstrate milliarcsecond-level astrometry using observations of two known bright asteroids on the same system with synthetic tracking. We conclude by discussing strategies for scheduling observations to detect and characterize small and fast-moving NEAs using the new technique.

We quantify the effects of refraction in transit transmission spectroscopy on spectral absorption features and on temporal variations that could be used to obtain altitude-dependent spectra for planets orbiting stars of different stellar types. We validate our model against altitude-dependent transmission spectra of the Earth from ATMOS and against lunar eclipse spectra from Pallé et al. We perform detectability studies to show the potential effects of refraction on hypothetical observations of Earth analogs with the James Webb Space Telescope NIRSPEC. Due to refraction, there will be a maximum tangent pressure level that can be probed during transit for each given planet–star system. We show that because of refraction, for an Earth-analog planet orbiting in the habitable zone of a Sun-like star only the top 0.3 bars of the atmosphere can be probed, leading to a decrease in the signal-to-noise ratio (S/N) of absorption features by 60%, while for an Earth-analog planet orbiting in the habitable zone of an M5V star it is possible to probe almost the entire atmosphere with minimal decreases in S/N. We also show that refraction can result in temporal variations in the transit transmission spectrum which may provide a way to obtain altitude-dependent spectra of exoplanet atmospheres. Additionally, the variations prior to ingress and subsequent to egress provide a way to probe pressures greater than the maximum tangent pressure that can be probed during transit. Therefore, probing the maximum range of atmospheric altitudes, and in particular the near-surface environment of an Earth-analog exoplanet, will require looking at out-of-transit refracted light in addition to the in-transit spectrum.

We report on the design and performance of the Bicep2 instrument and on its three-year data set. Bicep2 was designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 1°–5°(ℓ = 40–200), near the expected peak of the B-mode polarization signature of primordial gravitational waves from cosmic inflation. Measuring B-modes requires dramatic improvements in sensitivity combined with exquisite control of systematics. The Bicep2 telescope observed from the South Pole with a 26 cm aperture and cold, on-axis, refractive optics. Bicep2 also adopted a new detector design in which beam-defining slot antenna arrays couple to transition-edge sensor (TES) bolometers, all fabricated on a common substrate. The antenna-coupled TES detectors supported scalable fabrication and multiplexed readout that allowed Bicep2 to achieve a high detector count of 500 bolometers at 150 GHz, giving unprecedented sensitivity to B-modes at degree angular scales. After optimization of detector and readout parameters, Bicep2 achieved an instrument noise-equivalent temperature of $15.8\ \mu \mathrm{K}\sqrt{\mathrm{s}}$. The full data set reached Stokes Q and U map depths of 87.2 nK in square-degree pixels (52 μK) over an effective area of 384 deg² within a 1000 deg² field. These are the deepest CMB polarization maps at degree angular scales to date. The power spectrum analysis presented in a companion paper has resulted in a significant detection of B-mode polarization at degree scales.

We performed sub-parsec (∼0.06 pc) scale simulations of two idealized molecular clouds with different masses undergoing a collision. Gas clumps with densities greater than 10⁻²⁰ g cm⁻³ (0.3 × 10⁴ cm⁻³) were identified as pre-stellar cores and tracked throughout the simulation. The colliding system showed a partial gas arc morphology with core formation in the oblique shock front at the collision interface. These characteristics support NANTEN observations of objects suspected to be colliding giant molecular clouds (GMCs). We investigated the effect of turbulence and collision speed on the resulting core population and compared the cumulative mass distribution to cores in observed GMCs. Our results suggest that a faster relative velocity increases the number of cores formed but that cores grow via accretion predominately while in the shock front, leading to a slower shock being more important for core growth. The core masses obey a power-law relation with index γ = −1.6, in good agreement with observations. This suggests that core production through collisions should follow a similar mass distribution as quiescent formation, albeit at a higher mass range. If cores can be supported against collapse during their growth, then the estimated ram pressure from gas infall is of the right order to counter the radiation pressure and form a star of 100 M_☉.

We report the detection of an anomalous accretion flare in the tight eccentric pre-main-sequence binary system DQ Tau. In a multi-epoch survey consisting of randomly acquired low- to moderate-resolution near-infrared spectra obtained over a period of almost 10 yr, we detect a significant and simultaneous brightening of four standard accretion indicators (Ca ii infrared triplet, the Paschen and Brackett series H i lines, and He i 1.083 μm), on back-to-back nights (ϕ = 0.372 and 0.433) with the flare increasing in strength as the system approached apastron (ϕ = 0.5). The mass accretion rate measured for the anomalous flare is nearly an order of magnitude stronger than the average quiescent rate. While previous observations established that frequent, periodic accretion flares phased with periastron passages occur in this system, these data provide evidence that orbitally modulated accretion flares occur near apastron, when the stars make their closest approach to the circumbinary disk. The timing of the flare suggests that this outburst is due to interactions of the stellar cores (or the highly truncated circumstellar disks) with material in non-axisymmetric structures located at the inner edge of the circumbinary disk. We also explore the optical/infrared spectral type mismatch previously observed for T Tauri stars (TTSs) and successfully model the shape of the spectra from 0.8 to 1.0 μm and the strengths of the TiO and FeH bands as manifestations of large cool spots on the surfaces of the stellar companions in DQ Tau. These findings illustrate that a complete model of near-infrared spectra of many TTSs must include parameters for spot filling factors and temperatures.

The radii of debris disks and the sizes of their dust grains are important tracers of the planetesimal formation mechanisms and physical processes operating in these systems. Here we use a representative sample of 34 debris disks resolved in various Herschel Space Observatory (Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA) programs to constrain the disk radii and the size distribution of their dust. While we modeled disks with both warm and cold components, and identified warm inner disks around about two-thirds of the stars, we focus our analysis only on the cold outer disks, i.e., Kuiper-belt analogs. We derive the disk radii from the resolved images and find a large dispersion for host stars of any spectral class, but no significant trend with the stellar luminosity. This argues against ice lines as a dominant player in setting the debris disk sizes, since the ice line location varies with the luminosity of the central star. Fixing the disk radii to those inferred from the resolved images, we model the spectral energy distribution to determine the dust temperature and the grain size distribution for each target. While the dust temperature systematically increases toward earlier spectral types, the ratio of the dust temperature to the blackbody temperature at the disk radius decreases with the stellar luminosity. This is explained by a clear trend of typical sizes increasing toward more luminous stars. The typical grain sizes are compared to the radiation pressure blowout limit s_blow that is proportional to the stellar luminosity-to-mass ratio and thus also increases toward earlier spectral classes. The grain sizes in the disks of G- to A-stars are inferred to be several times s_blow at all stellar luminosities, in agreement with collisional models of debris disks. The sizes, measured in the units of s_blow, appear to decrease with the luminosity, which may be suggestive of the disk's stirring level increasing toward earlier-type stars. The dust opacity index β ranges between zero and two, and the size distribution index q varies between three and five for all the disks in the sample.

The majority of massive stars are formed in binary systems. It is hence reasonable to expect that most core-collapse supernovae (CCSNe) take place in binaries and the existence of a companion star may leave some imprints in observed features. Having this in mind, we have conducted two-dimensional hydrodynamical simulations of the collisions of CCSNe ejecta with the companion star in an almost-equal-mass (∼10 M_☉) binary to find out possible consequences of such events. In particular we pay attention to the amount of mass removed and its dependence on the binary separation. In contrast to the previous surmise, we find that the companion mass is stripped not by momentum transfer but by shock heating. Up to 25% of the original mass can be removed for the closest separations and the removed mass decreases as M_ub∝a^−4.3 with the binary separation a. By performing some experimental computations with artificially modified densities of incident ejecta, we show that if the velocity of ejecta is fixed, the density of incident ejecta is the single important parameter that actually determines the removed mass as $M_{{\rm ub}} \propto \rho _{{\rm ej}} ^{1.4}$. On the other hand, another set of simulations with modified velocities of incident ejecta demonstrate that the strength of the forward shock, which heats up the stellar material and causes the mass loss of the companion star, is actually the key parameter for the removed mass.

Kepler data from G-, K-, and M-type stars are used to study conditions that lead to superflares with energies above 10³⁴ erg. From the 117,661 stars included, 380 show superflares with a total of 1690 such events. We study whether parameters, like effective temperature or rotation rate, have any effect on the superflare occurrence rate or energy. With increasing effective temperature we observe a decrease in the superflare rate, which is analogous to the previous findings of a decrease in dynamo activity with increasing effective temperature. For slowly rotating stars, we find a quadratic increase of the mean occurrence rate with the rotation rate up to a critical point, after which the rate decreases linearly. Motivated by standard dynamo theory, we study the behavior of the relative starspot coverage, approximated as the relative brightness variation. For faster rotating stars, an increased fraction of stars shows higher spot coverage, which leads to higher superflare rates. A turbulent dynamo is used to study the dependence of the Ohmic dissipation as a proxy of the flare energy on the differential rotation or shear rate. The resulting statistics of the dissipation energy as a function of dynamo number is similar to the observed flare statistics as a function of the inverse Rossby number and shows similarly strong fluctuations. This supports the idea that superflares might well be possible for solar-type G stars.

In this study, we present high-resolution millimeter observations of the dust and gas disk of the T Tauri star AS 205 N and its companion, AS 205 S, obtained with the Atacama Large Millimeter Array. The gas disk around AS 205 N, for which infrared emission spectroscopy demonstrates significant deviations from Keplerian motion that has been interpreted as evidence for a disk wind, also displays significant deviations from Keplerian disk emission in the observations presented here. Detections near both AS 205 N and S are obtained in 1.3 mm continuum, ¹²CO 2–1, ¹³CO 2–1, and C¹⁸O 2–1. The ¹²CO emission is extended up to ∼2'' from AS 205 N, and both ¹²CO and ¹³CO display deviations from Keplerian rotation at all angular scales. Two possible explanations for these observations hold up best to close scrutiny—tidal interaction with AS 205 S or disk winds (or a combination of the two)—and we discuss these possibilities in some detail.

Motivated by its important role in the collisional growth of dust particles in protoplanetary disks, we investigate the probability distribution function (PDF) of the relative velocity of inertial particles suspended in turbulent flows. Using the simulation from our previous work, we compute the relative velocity PDF as a function of the friction timescales, τ_p1 and τ_p2, of two particles of arbitrary sizes. The friction time of the particles included in the simulation ranges from 0.1τ_η to 54T_L, where τ_η and T_L are the Kolmogorov time and the Lagrangian correlation time of the flow, respectively. The relative velocity PDF is generically non-Gaussian, exhibiting fat tails. For a fixed value of τ_p1, the PDF shape is the fattest for equal-size particles (τ_p2 = τ_p1), and becomes thinner at both τ_p2 < τ_p1 and τ_p2 > τ_p1. Defining f as the friction time ratio of the smaller particle to the larger one, we find that, at a given f in (1/2) ≲ f ≲ 1, the PDF fatness first increases with the friction time τ_{p, h} of the larger particle, peaks at τ_{p, h} ≃ τ_η, and then decreases as τ_{p, h} increases further. For 0 ⩽ f ≲ (1/4), the PDF becomes continuously thinner with increasing τ_{p, h}. The PDF is nearly Gaussian only if τ_{p, h} is sufficiently large (≫T_L). These features are successfully explained by the Pan & Padoan model. Using our simulation data and some simplifying assumptions, we estimated the fractions of collisions resulting in sticking, bouncing, and fragmentation as a function of the dust size in protoplanetary disks, and argued that accounting for non-Gaussianity of the collision velocity may help further alleviate the bouncing barrier problem.

We derive the conductivity tensor for axisymmetric perturbations of a hot, collisionless, and charge-neutral plasma in the shearing sheet approximation. Our results generalize the well-known linear Vlasov theory for uniform plasmas to differentially rotating plasmas and can be used for wide range of kinetic stability calculations. We apply these results to the linear theory of the magneto-rotational instability (MRI) in collisionless plasmas. We show analytically and numerically how the general kinetic theory results derived here reduce in appropriate limits to previous results in the literature, including the low-frequency guiding center (or "kinetic MHD") approximation, Hall magnetohydrodynamics (MHD), and the gyro-viscous approximation. We revisit the cold plasma model of the MRI and show that, contrary to previous results, an initially unmagnetized collisionless plasma is linearly stable to axisymmetric perturbations in the cold plasma approximation. In addition to their application to astrophysical plasmas, our results provide a useful framework for assessing the linear stability of differentially rotating plasmas in laboratory experiments.

The nature of the progenitor systems of type Ia supernovae (SNe Ia) is still unclear. One way to distinguish between the single-degenerate scenario and double-degenerate scenario for their progenitors is to search for the surviving companions (SCs). Using a technique that couples the results from multi-dimensional hydrodynamics simulations with calculations of the structure and evolution of main-sequence- (MS-) and helium-rich SCs, the color and magnitude of MS- and helium-rich SCs are predicted as functions of time. The SC candidates in Galactic type Ia supernova remnants (Ia SNR) and nearby extragalactic Ia SNRs are discussed. We find that the maximum detectable distance of MS SCs (helium-rich SCs) is 0.6–4 Mpc (0.4–16 Mpc), if the apparent magnitude limit is 27 in the absence of extinction, suggesting that the Large and Small Magellanic Clouds and the Andromeda Galaxy are excellent environments in which to search for SCs. However, only five Ia SNRs have been searched for SCs, showing little support for the standard channels in the singe-degenerate scenario. To better understand the progenitors of SNe Ia, we encourage the search for SCs in other nearby Ia SNRs.

The Sloan Digital Sky Survey (SDSS) Galaxy Zoo project has revealed a number of spectacular galaxies possessing extended emission-line regions (EELRs), the most famous being Hanny's Voorwerp galaxy. We present another EELR object discovered in the SDSS endeavor: the Teacup active galactic nucleus (AGN). Nicknamed for its EELR, which has a "handle"-like structure protruding 15 kpc into the northeast quadrant of the galaxy. We analyze the physical conditions of this galaxy with long-slit, ground-based spectroscopy from the Lowell, Lick, and KPNO observatories. With the Lowell 1.8 m Perkin's telescope we took multiple observations at different offset positions, allowing us to recover spatially resolved spectra across the galaxy. Line diagnostics indicate the ionized gas is photoionized primarily by the AGN. Additionally we are able to derive the hydrogen density from the [S ii] λ6716/λ6731 ratio. We generated two-component photoionization models for each spatially resolved Lowell spectrum. These models allow us to calculate the AGN bolometric luminosity seen by the gas at different radii from the nuclear center of the Teacup. Our results show a drop in bolometric luminosity by more than two orders of magnitude from the EELR to the nucleus, suggesting that the AGN has decreased in luminosity by this amount in a continuous fashion over 46,000 yr, supporting the case for a dying AGN in this galaxy independent of any IR based evidence. We demonstrate that spatially resolved photoionization modeling could be applied to EELRs to investigate long timescale variability.

Photophoretic motion can transport illuminated particles in protoplanetary disks. In a previous paper, we focused on the modeling of steady state photophoretic forces based on the compositions derived from tomography and heat transfer. Here, we present microgravity experiments which deviate significantly from the steady state calculations of the first paper. The experiments on average show a significantly smaller force than predicted with a large variation in absolute photophoretic force and in the direction of motion with respect to the illumination. Time-dependent modeling of photophoretic forces for heat-up and rotation shows that the variations in strength and direction observed can be well explained by the particle reorientation in the limited experiment time of a drop tower experiment. In protoplanetary disks, random rotation subsides due to gas friction on short timescales and the results of our earlier paper hold. Rotation has a significant influence in short duration laboratory studies. Observing particle motion and rotation under the influence of photophoresis can be considered as a basic laboratory analog experiment to Yarkovsky and YORP effects.

We investigated the magnitude–phase relation of (162173) 1999 JU3, a target asteroid for the JAXA Hayabusa 2 sample return mission. We initially employed the International Astronomical Union's H–G formalism but found that it fits less well using a single set of parameters. To improve the inadequate fit, we employed two photometric functions: the Shevchenko and Hapke functions. With the Shevchenko function, we found that the magnitude–phase relation exhibits linear behavior in a wide phase angle range (α = 5°–75°) and shows weak nonlinear opposition brightening at α < 5°, providing a more reliable absolute magnitude of H_V = 19.25 ± 0.03. The phase slope (0.039 ± 0.001 mag deg⁻¹) and opposition effect amplitude (parameterized by the ratio of intensity at α = 03 to that at α = 5°, I(03)/I(5°) = 1.31 ± 0.05) are consistent with those of typical C-type asteroids. We also attempted to determine the parameters for the Hapke model, which are applicable for constructing the surface reflectance map with the Hayabusa 2 onboard cameras. Although we could not constrain the full set of Hapke parameters, we obtained possible values, w = 0.041, g = −0.38, B₀ = 1.43, and h = 0.050, assuming a surface roughness parameter $\bar{\theta }$ = 20°. By combining our photometric study with a thermal model of the asteroid, we obtained a geometric albedo of p_v = 0.047 ± 0.003, phase integral q = 0.32 ± 0.03, and Bond albedo A_B = 0.014 ± 0.002, which are commensurate with the values for common C-type asteroids.

We investigate the relationships between stellar mass, gas-phase oxygen abundance (metallicity), star formation rate (SFR), and dust content of star-forming galaxies at z ∼ 1.6 using Subaru/FMOS spectroscopy in the COSMOS field. The mass–metallicity (MZ) relation at z ∼ 1.6 is steeper than the relation observed in the local universe. The steeper MZ relation at z ∼ 1.6 is mainly due to evolution in the stellar mass where the MZ relation begins to turnover and flatten. This turnover mass is 1.2 dex larger at z ∼ 1.6. The most massive galaxies at z ∼ 1.6 (∼10¹¹ M_☉) are enriched to the level observed in massive galaxies in the local universe. The MZ relation we measure at z ∼ 1.6 supports the suggestion of an empirical upper metallicity limit that does not significantly evolve with redshift. We find an anti-correlation between metallicity and SFR for galaxies at a fixed stellar mass at z ∼ 1.6, which is similar to trends observed in the local universe. We do not find a relation between stellar mass, metallicity, and SFR that is independent of redshift; rather, our data suggest that there is redshift evolution in this relation. We examine the relation between stellar mass, metallicity, and dust extinction, and find that at a fixed stellar mass, dustier galaxies tend to be more metal rich. From examination of the stellar masses, metallicities, SFRs, and dust extinctions, we conclude that stellar mass is most closely related to dust extinction.

We utilize 16 band Hubble Space Telescope (HST) observations of 18 lensing clusters obtained as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program to search for z ∼ 6–8 galaxies. We report the discovery of 204, 45, and 13 Lyman-break galaxy candidates at z ∼ 6, z ∼ 7, and z ∼ 8, respectively, identified from purely photometric redshift selections. This large sample, representing nearly an order of magnitude increase in the number of magnified star-forming galaxies at z ∼ 6–8 presented to date, is unique in that we have observations in four WFC3/UVIS UV, seven ACS/WFC optical, and all five WFC3/IR broadband filters, which enable very accurate photometric redshift selections. We construct detailed lensing models for 17 of the 18 clusters to estimate object magnifications and to identify two new multiply lensed z ≳ 6 candidates. The median magnifications over the 17 clusters are 4, 4, and 5 for the z ∼ 6, z ∼ 7, and z ∼ 8 samples, respectively, over an average area of 4.5 arcmin² per cluster. We compare our observed number counts with expectations based on convolving "blank" field UV luminosity functions through our cluster lens models and find rough agreement down to ∼27 mag, where we begin to suffer significant incompleteness. In all three redshift bins, we find a higher number density at brighter observed magnitudes than the field predictions, empirically demonstrating for the first time the enhanced efficiency of lensing clusters over field surveys. Our number counts also are in general agreement with the lensed expectations from the cluster models, especially at z ∼ 6, where we have the best statistics.

We have monitored 12 intrinsic narrow absorption lines (NALs) in five quasars and seven mini-broad absorption lines (mini-BALs) in six quasars for a period of 4–12 yr (1–3.5 yr in the quasar rest-frame). We present the observational data and the conclusions that follow immediately from them, as a prelude to a more detailed analysis. We found clear variability in the equivalent widths (EWs) of the mini-BAL systems but no easily discernible changes in their profiles. We did not detect any variability in the NAL systems or in narrow components that are often located at the center of mini-BAL profiles. Variations in mini-BAL EWs are larger at longer time intervals, reminiscent of the trend seen in variable BALs. If we assume that the observed variations result from changes in the ionization state of the mini-BAL gas, we infer lower limits to the gas density ∼10³–10⁵ cm⁻³ and upper limits on the distance of the absorbers from the central engine of the order of a few kiloparsecs. Motivated by the observed variability properties, we suggest that mini-BALs can vary because of fluctuations of the ionizing continuum or changes in partial coverage while NALs can vary primarily because of changes in partial coverage.

Radiative feedback produced by stellar populations played a vital role in early structure formation. In particular, photons below the Lyman limit can escape the star-forming regions and produce a background ultraviolet (UV) flux, which consequently may influence the pristine halos far away from the radiation sources. These photons can quench the formation of molecular hydrogen by photodetachment of H⁻. In this study, we explore the impact of such UV radiation on fragmentation in massive primordial halos of a few times 10⁷M_☉. To accomplish this goal, we perform high resolution cosmological simulations for two distinct halos and vary the strength of the impinging background UV field in units of J₂₁ assuming a blackbody radiation spectrum with a characteristic temperature of T_rad = 10⁴ K. We further make use of sink particles to follow the evolution for 10,000 yr after reaching the maximum refinement level. No vigorous fragmentation is observed in UV-illuminated halos while the accretion rate changes according to the thermal properties. Our findings show that a few 10²–10⁴ solar mass protostars are formed when halos are irradiated by J₂₁ = 10–500 at z > 10 and suggest a strong relation between the strength of the UV flux and mass of a protostar. This mode of star formation is quite different from minihalos, as higher accretion rates of about 0.01–0.1 M_☉ yr⁻¹ are observed by the end of our simulations. The resulting massive stars are potential cradles for the formation of intermediate-mass black holes at earlier cosmic times and contribute to the formation of a global X-ray background.

Transit and radial velocity searches are two techniques for identifying nearby extrasolar planets to Earth that transit bright stars. Identifying a robust sample of these exoplanets around bright stars for detailed atmospheric characterization is a major observational undertaking. In this study we describe a framework that answers the question of whether a transit or radial velocity survey is more efficient at finding transiting exoplanets given the same amount of observing time. Within the framework we show that a transit survey's window function can be approximated using the hypergeometric probability distribution. We estimate the observing time required for a transit survey to find a transiting Earth-sized exoplanet in the habitable zone (HZ) with an emphasis on late-type stars. We also estimate the radial velocity precision necessary to detect the equivalent HZ Earth-mass exoplanet that also transits when using an equal amount of observing time as the transit survey. We find that a radial velocity survey with σ_rv ∼ 0.6 m s⁻¹ precision has comparable efficiency in terms of observing time to a transit survey with the requisite photometric precision σ_phot ∼ 300 ppm to find a transiting Earth-sized exoplanet in the HZ of late M dwarfs. For super-Earths, a σ_rv ∼ 2.0 m s⁻¹ precision radial velocity survey has comparable efficiency to a transit survey with σ_phot ∼ 2300 ppm.