Table of contents

Gamma-ray burst GRB 140430A was detected by the Swift satellite and observed promptly with the imaging polarimeter RINGO3 mounted on the Liverpool Telescope, with observations beginning while the prompt γ-ray emission was still ongoing. In this paper, we present densely sampled (10-s temporal resolution) early optical light curves (LCs) in 3 optical bands and limits to the degree of optical polarization. We compare optical, X-ray, and gamma-ray properties and present an analysis of the optical emission during a period of high-energy flaring. The complex optical LC cannot be explained merely with a combination of forward and reverse shock emission from a standard external shock, implying additional contribution of emission from internal shock dissipation. We estimate an upper limit for time averaged optical polarization during the prompt phase to be as low as P < 12% (1σ). This suggests that the optical flares and early afterglow emission in this GRB are not highly polarized. Alternatively, time averaging could mask the presence of otherwise polarized components of distinct origin at different polarization position angles.

We present a comprehensive nucleosynthesis study of the neutrino-driven wind in the aftermath of a binary neutron star merger. Our focus is the initial remnant phase when a massive central neutron star is present. Using tracers from a recent hydrodynamical simulation, we determine total masses and integrated abundances to characterize the composition of unbound matter. We find that the nucleosynthetic yields depend sensitively on both the life time of the massive neutron star and the polar angle. Matter in excess of up to 9 × 10⁻³M_⊙ becomes unbound until ∼200 ms. Due to electron fractions of Y_e ≈ 0.2–0.4, mainly nuclei with mass numbers A < 130 are synthesized, complementing the yields from the earlier dynamic ejecta. Mixing scenarios with these two types of ejecta can explain the abundance pattern in r-process enriched metal-poor stars. Additionally, we calculate heating rates for the decay of the freshly produced radioactive isotopes. The resulting light curve peaks in the blue band after about 4 hr. Furthermore, high opacities due to heavy r-process nuclei in the dynamic ejecta lead to a second peak in the infrared after 3–4 days.

G73.9+0.9 has been classified as a probable shell-type supernova remnant, though it has also been suggested that it could have a pulsar wind nebula (PWN). Here, a broadband model of the non-thermal emission of G73.9+0.9 from radio to gamma rays is presented. The model includes a new gamma-ray observation obtained by the analysis of seven years of data from the Fermi/LAT telescope. Above 200 MeV, the source is detected with a significance of 13σ and the spectrum of the radiation is best described by a power law with an index of ∼2.5. The leptonic mechanisms are hard to reconcile with the measured radio and gamma-ray spectral energy distribution. A PWN origin for the high-energy emission is also not very likely, due to the lack of detection of pulsars and of X-ray emission in the region, as well as from the shape of the gamma-ray spectrum. Given the possibility that the object is interacting with molecular clouds, a hadronic origin of the high-energy emission is more likely, and the spectral properties of the cosmic rays responsible for this radiation are derived.

Many sulfur-bearing species have been detected in different astronomical environments and have allowed us to derive important information about the chemical and physical composition of interstellar regions. In particular, these species have also been shown to trace and probe hot-core environment time evolution. Among the most prominent sulfur-bearing molecules is SO, the sulfur monoxide radical, one of the more ubiquitous and abundant, which is also observed in its isotopic substituted species such as ³⁴SO and S¹⁸O. Due to the importance of this simple diatomic system, and in order to face the challenge of modern radioastronomical facilities, an extension to the THz range of the rare isotopologues of sulfur monoxide has been performed. High-resolution rotational molecular spectroscopy has been employed to extend the available data set of four isotopic species, SO, ³⁴SO, S¹⁷O, and S¹⁸O, up to the 1.5 THz region. The frequency coverage and spectral resolution of our measurements allowed a better constraint of the molecular constants of the four species considered, specifically focusing on the two oxygen-substituted isotopologues. Our measurements were also employed in an isotopically invariant fit including all of the available pure rotational and ro-vibrational transitions for all of the SO isotopologues, thus enabling accurate predictions of the rotational transitions at higher frequencies. We also provide comparisons with recent works performed on the same system, demonstrating the quality of our experiment and the improvement of the data sets for all of the species considered. Transition frequencies for this system can now be used with confidence by the astronomical community well into the THz spectral region.

We investigate heating and acceleration of protons from a thermal gas with a generic diffusion and acceleration model, and subject to Coulomb scattering and energy loss, as was done by Petrosian & East for electrons. As protons gain energy their loss to electrons becomes important. Thus, we need to solve the coupled proton–electron kinetic equation. We numerically solve the coupled Fokker–Planck equations and compute the time evolution of the spectra of both particles. We show that this can lead to a quasi-thermal component plus a high-energy nonthermal tail. We determine the evolution of the nonthermal tail and the quasi-thermal component. The results may be used to explore the possibility of inverse bremsstrahlung radiation as a source of hard X-ray emissions from hot sources such as solar flares, accretion disk coronas, and the intracluster medium of galaxy clusters. We find that the emergence of nonthermal protons is accompanied by excessive heating of the entire plasma, unless the turbulence needed for scattering and acceleration is steeper than Kolmogorov and the acceleration parameters, the duration of the acceleration, and/or the initial distributions are significantly fine-tuned. These results severely constrain the feasibility of the nonthermal inverse bremsstrahlung process producing hard X-ray emissions. However, the nonthermal tail may be the seed particles for further re-acceleration to relativistic energies, say by a shock. In the Appendix we present some tests of the integrity of the algorithm used and present a new formula for the energy loss rate due to inelastic proton–proton interactions.

The distribution of galaxies displays anisotropy on different scales and it is often referred to as galaxy alignment. To understand the origin of galaxy alignments on small scales, one must investigate how galaxies were accreted in the early universe and quantify their primordial anisotropy at the time of accretion. In this paper we use N-body simulations to investigate the accretion of subhalos, focusing on their alignment with halo shape and the orientation of mass distribution on the large scale, defined using the Hessian matrix of the density field. The large/small (e1/e3) eigenvalues of the Hessian matrix define the fast/slow collapse direction of matter on the large scale. We find that: (1) the halo major axis is well aligned with the e3 (slow collapse) direction, and it is stronger for massive halos; (2) subhalos are predominantly accreted along the major axis of the host halo, and the alignment increases with the host halo mass. Most importantly, this alignment is universal; (3) accretion of subhalos with respect to the e3 direction is not universal. In massive halos, subhalos are accreted along the e3 (even more strongly than the alignment with the halo major axis), but in low-mass halos subhalos are accreted perpendicular to e3. The transitional mass is lower at high redshift. The last result well explains the puzzling correlation (both in recent observations and simulations) that massive galaxies/halos have their spin perpendicular to the filament, and the spin of low-mass galaxies/halos is slightly aligned with the filament, under the assumption that the orbital angular momentum of subhalos is converted to halo spin.

Using spectroscopically selected galaxies from the Sloan Digital Sky Survey we present a detection of reddening effects from the circumgalactic medium of galaxies which we attribute to an extended distribution of dust. We detect the mean change in the colors of "standard crayons" correlated with the presence of foreground galaxies at $z\sim 0.05$ as a function of angular separation. Following Peek & Graves, we create standard crayons using passively evolving galaxies corrected for Milky Way reddening and color-redshift trends, leading to a sample with as little as 2% scatter in color. We devise methods to ameliorate possible systematic effects related to the estimation of colors, and we find an excess reddening induced by foreground galaxies at a level ranging from 10 to 0.5 mmag on scales ranging from 30 kpc to 1 Mpc. We attribute this effect to a large-scale distribution of dust around galaxies similar to the findings of Ménard et al. We find that circumgalactic reddening is a weak function of stellar mass over the range $6\times {10}^{9}\;{M}_{\odot }$–$6\times {10}^{10}\;{M}_{\odot }$ and note that this behavior appears to be consistent with recent results on the distribution of metals in the gas phase. We also find that circumgalactic reddening has no detectable dependence on the specific star formation rate of the host galaxy.

Models of the dynamical evolution of the early solar system that follow the dispersal of the gaseous protoplanetary disk have been widely successful in reconstructing the current orbital configuration of the giant planets. Statistically, some of the most successful dynamical evolution simulations have initially included a hypothetical fifth giant planet, of ice giant (IG) mass, which gets ejected by a gas giant during the early solar system's proposed instability phase. We investigate the likelihood of an IG ejection (IGE) event by either Jupiter or Saturn through constraints imposed by the current orbits of their wide-separation regular satellites Callisto and Iapetus, respectively. We show that planetary encounters that are sufficient to eject an IG often provide excessive perturbations to the orbits of Callisto and Iapetus, making it difficult to reconcile a planet ejection event with the current orbit of either satellite. Quantitatively, we compute the likelihood of reconciling a regular Jovian satellite orbit with the current orbit of Callisto following an IGE by Jupiter of ∼42%, and conclude that such a large likelihood supports the hypothesis of a fifth giant planet's existence. A similar calculation for Iapetus reveals that it is much more difficult for Saturn to have ejected an IG and reconciled a Kronian satellite orbit with that of Iapetus (likelihood ∼1%), although uncertainties regarding the formation of Iapetus, with its unusual orbit, complicates the interpretation of this result.

What can we tell about exoplanet habitability if currently only the stellar properties, planet radius, and the incoming stellar flux are known? A planet is in the habitable zone (HZ) if it harbors liquid water on its surface. The HZ is traditionally conceived as a sharp region around stars because it is calculated for one planet with specific properties. Such an approach is limiting because the planet's atmospheric and geophysical properties, which influence the presence of liquid water on the surface, are currently unknown but expected to be diverse. A statistical HZ description is outlined that does not favor one planet type. Instead, the stellar and planet properties are treated as random variables, and a continuous range of planet scenarios is considered. Various probability density functions are assigned to each random variable, and a combination of Monte Carlo sampling and climate modeling is used to generate synthetic exoplanet populations with known surface climates. Then, the properties of the subpopulation bearing liquid water is analyzed. Given our current observational knowledge, the HZ takes the form of a weakly constrained but smooth probability function. The HZ has an inner edge, but a clear outer edge is not seen. Currently only optimistic upper limits can be derived for the potentially observable HZ occurrence rate. Finally, we illustrate through an example how future data on exoplanet atmospheres will help to narrow down the probability that an exoplanet harbors liquid water, and we identify the greatest observational challenge in the way of finding a habitable exoplanet.

Motivated by the discovery of prolate rotation of stars in Andromeda II (And II), a dwarf spheroidal companion of M31, we study its origin via mergers of disky dwarf galaxies. We simulate merger events between two identical dwarfs changing the initial inclination of their disks with respect to the orbit and the amount of orbital angular momentum. On radial orbits, the amount of prolate rotation in the merger remnants correlates strongly with the inclination of the disks and is well understood as due to the conservation of the angular momentum component of the disks along the merger axis. For non-radial orbits, prolate rotation may still be produced if the orbital angular momentum is initially not much larger than the intrinsic angular momentum of the disks. The orbital structure of the remnants with significant rotation is dominated by box orbits in the center and long-axis tubes in the outer parts. The frequency analysis of stellar orbits in the plane perpendicular to the major axis reveals the presence of two families roughly corresponding to inner and outer long-axis tubes. The fraction of inner tubes is largest in the remnant forming from disks that are initially oriented most vertically, and is responsible for the boxy shape of the galaxy. We conclude that prolate rotation results from mergers with a variety of initial conditions and no fine tuning is necessary to reproduce this feature. We compare the properties of our merger remnants to those of dwarfs resulting from the tidal stirring scenario and the data for And II.

Following our previous work, which related generic features in the sky-averaged (global) 21-cm signal to properties of the intergalactic medium, we now investigate the prospects for constraining a simple galaxy formation model with current and near-future experiments. Markov-Chain Monte Carlo fits to our synthetic data set, which includes a realistic galactic foreground, a plausible model for the signal, and noise consistent with 100 hr of integration by an ideal instrument, suggest that a simple four-parameter model that links the production rate of Lyα, Lyman-continuum, and X-ray photons to the growth rate of dark matter halos can be well-constrained (to ∼0.1 dex in each dimension) so long as all three spectral features expected to occur between 40 ≲ ν/MHz ≲ 120 are detected. Several important conclusions follow naturally from this basic numerical result, namely that measurements of the global 21-cm signal can in principle (i) identify the characteristic halo mass threshold for star formation at all redshifts z ≳ 15, (ii) extend z ≲ 4 upper limits on the normalization of the X-ray luminosity star formation rate (L_X–SFR) relation out to z ∼ 20, and (iii) provide joint constraints on stellar spectra and the escape fraction of ionizing radiation at z ∼ 12. Though our approach is general, the importance of a broadband measurement renders our findings most relevant to the proposed Dark Ages Radio Explorer, which will have a clean view of the global 21-cm signal from ∼40 to 120 MHz from its vantage point above the radio-quiet, ionosphere-free lunar far-side.

Hot methane is found in many "cool" sub-stellar astronomical sources including brown dwarfs and exoplanets, as well as in combustion environments on Earth. We report on the first high-resolution laboratory absorption spectra of hot methane at temperatures up to 1200 K. Our observations are compared to the latest theoretical spectral predictions and recent brown dwarf spectra. The expectation that millions of weak absorption lines combine to form a continuum, not seen at room temperature, is confirmed. Our high-resolution transmittance spectra account for both the emission and absorption of methane at elevated temperatures. From these spectra, we obtain an empirical line list and continuum that is able to account for the absorption of methane in high temperature environments at both high and low resolution. Great advances have recently been made in the theoretical prediction of hot methane, and our experimental measurements highlight the progress made and the problems that still remain.

${\mathcal{T}}$-REx (Tau Retrieval of Exoplanets) is a novel, fully Bayesian atmospheric retrieval code custom built for extrasolar atmospheres. In Waldmann et al., the transmission spectroscopic case was introduced, and here we present the emission spectroscopy spectral retrieval for the ${\mathcal{T}}$-REx framework. Compared to transmission spectroscopy, the emission case is often significantly more degenerate due to the need to retrieve the full atmospheric temperature–pressure (TP) profile. This is particularly true in the case of current measurements of exoplanetary atmospheres, which are either of low signal-to-noise, low spectral resolution, or both. We present a new way of combining two existing approaches to the modeling of the said TP profile: (1) the parametric profile, where the atmospheric TP structure is analytically approximated by a few model parameters, (2) the layer-by-layer approach, where individual atmospheric layers are modeled. Both of these approaches have distinct advantages and disadvantages in terms of convergence properties and potential model biases. The ${\mathcal{T}}$-REx hybrid model presented here is a new two-stage TP profile retrieval, which combines the robustness of the analytic solution with the accuracy of the layer-by-layer approach. The retrieval process is demonstrated using simulations of the hot-Jupiter WASP-76b and the hot-super-Earth 55 Cnc e as well as the secondary eclipse measurements of HD 189733b.

We present a survey of 41 Kepler Objects of Interest (KOIs) for exomoons using Bayesian photodynamics, more than tripling the number of KOIs surveyed with this technique. We find no compelling evidence for exomoons although 13 KOIs yield spurious detections driven by instrumental artifacts, stellar activity, and/or perturbations from unseen bodies. Regarding the latter, we find seven KOIs exhibiting >5 σ evidence of transit timing variations, including the "mega-Earth" Kepler-10c, likely indicating an additional planet in that system. We exploit the moderately large sample of 57 unique KOIs surveyed to date to infer several useful statistics. For example, although there is a diverse range in sensitivities, we find that we are sensitive to Pluto–Charon mass-ratio systems for ≃40% of KOIs studied and Earth–Moon mass-ratios for 1 in 8 cases. In terms of absolute mass, our limits probe down to 1.7 Ganymede masses, with a sensitivity to Earth-mass moons for 1 in 3 cases studied and to the smallest moons capable of sustaining an Earth-like atmosphere (0.3 M_⨁) for 1 in 4. Despite the lack of positive detections to date, we caution against drawing conclusions yet, since our most interesting objects remain under analysis. Finally, we point out that had we searched for the photometric transit signals of exomoons alone, rather than using photodynamics, we estimate that 1 in 4 KOIs would have erroneously been concluded to harbor exomoons due to residual time correlated noise in the Kepler data, posing a serious problem for alternative methods.

GJ1214b is a warm sub-Neptune transiting in front of a nearby M dwarf star. Recent observations indicate the presence of high and thick clouds or haze whose presence requires strong atmospheric mixing. In order to understand the transport and distribution of such clouds/haze, we study the atmospheric circulation and the vertical mixing of GJ1214b with a 3D General Circulation Model for cloud-free hydrogen-dominated atmospheres (metallicity of 1, 10, and 100 times the solar value) and for a water-dominated atmosphere. We analyze the effect of the atmospheric metallicity on the thermal structure and zonal winds. We also analyze the zonal mean meridional circulation and show that it corresponds to an anti-Hadley circulation in most of the atmosphere with upwelling at mid-latitude and downwelling at the equator on average. This circulation must be present on a large range of synchronously rotating exoplanets with a strong impact on cloud formation and distribution. Using simple tracers, we show that vertical winds on GJ1214b can be strong enough to loft micrometric particles and that the anti-Hadley circulation leads to a minimum of tracers at the equator. We find that the strength of the vertical mixing increases with metallicity. We derive 1D equivalent eddy diffusion coefficients and find simple parametrizations from ${K}_{{zz}}=7\times {10}^{2}\times {P}_{\mathrm{bar}}^{-0.4}\;{{\rm{m}}}^{2}\;{{\rm{s}}}^{-1}$ for solar metallicity to ${K}_{{zz}}=3\times {10}^{3}\times {P}_{\mathrm{bar}}^{-0.4}\;{{\rm{m}}}^{2}\;{{\rm{s}}}^{-1}$ for the 100× solar metallicity. These values should favor an efficient formation of photochemical haze in the upper atmosphere of GJ1214b.

We have measured electron impact ionization for Fe⁷⁺ from the ionization threshold up to 1200 eV. The measurements were performed using the TSR heavy ion storage ring. The ions were stored long enough prior to measurements to remove most metastables, resulting in a beam of 94% ground-level ions. Comparing with the previously recommended atomic data, we find that the Arnaud & Raymond cross section is up to about 40% larger than our measurement, with the largest discrepancies below about 400 eV. The cross section of Dere agrees to within 10%, which is about the magnitude of the experimental uncertainties. The remaining discrepancies between our measurement and the Dere calculations are likely due to shortcomings in the theoretical treatment of the excitation-autoionization contribution.

Lieu et al. have recently claimed that it is possible to substantially improve the sensitivity of radio-astronomical observations. In essence, their proposal is to make use of the intensity of the photon shot noise as a measure of the photon arrival rate. Lieu et al. provide a detailed quantum-mechanical calculation of a proposed measurement scheme that uses two detectors and conclude that this scheme avoids the sensitivity degradation that is associated with photon bunching. If correct, this result could have a profound impact on radio astronomy. Here I present a detailed analysis of the sensitivity attainable using shot-noise measurement schemes that use either one or two detectors, and demonstrate that neither scheme can avoid the photon bunching penalty. I perform both semiclassical and fully quantum calculations of the sensitivity, obtaining consistent results, and provide a formal proof of the equivalence of these two approaches. These direct calculations are furthermore shown to be consistent with an indirect argument based on a correlation method that establishes an independent limit to the sensitivity of shot-noise measurement schemes. Furthermore, these calculations are directly applicable to the regime of interest identified by Lieu et al. Collectively, these results conclusively demonstrate that the photon-bunching sensitivity penalty applies to shot-noise measurement schemes just as it does to ordinary photon counting, in contradiction to the fundamental claim made by Lieu et al. The source of this contradiction is traced to a logical fallacy in their argument.

The redshifted 21 cm line of neutral hydrogen (H i), potentially observable at low radio frequencies (∼50–200 MHz), is a promising probe of the physical conditions of the intergalactic medium during Cosmic Dawn and the Epoch of Reionization (EoR). The sky-averaged H i signal is expected to be extremely weak (∼100 mK) in comparison to the Galactic foreground emission (∼10⁴ K). Moreover, the sky-averaged spectra measured by ground-based instruments are affected by chromatic propagation effects (∼tens of kelvin) originating in the ionosphere. We analyze data collected with the upgraded Broadband Instrument for Global Hydrogen Reionization Signal system deployed at the Murchison Radio-astronomy Observatory to assess the significance of ionospheric effects on the detection of the global EoR signal. The ionospheric effects identified in these data are, particularly during nighttime, dominated by absorption and emission. We measure some properties of the ionosphere, such as the electron temperature (T_e ≈ 470 K at nighttime), magnitude, and variability of optical depth (τ_{100 MHz} ≈ 0.01 and δτ ≈ 0.005 at nighttime). According to the results of a statistical test applied on a large data sample, very long integrations (∼100 hr collected over approximately 2 months) lead to increased signal-to-noise ratio even in the presence of ionospheric variability. This is further supported by the structure of the power spectrum of the sky temperature fluctuations, which has flicker noise characteristics at frequencies ≳10⁻⁵ Hz, but becomes flat below ≈10⁻⁵ Hz. Hence, we conclude that the stochastic error introduced by the chromatic ionospheric effects tends to zero in an average. Therefore, the ionospheric effects and fluctuations are not fundamental impediments preventing ground-based instruments from integrating down to the precision required by global EoR experiments, provided that the ionospheric contribution is properly accounted for in the data analysis.

High-resolution galactic neutral hydrogen (HI) data obtained with the Green Bank Telescope (GBT) over 56 square degrees of sky around l = 132°, b = 25° are compared with small-scale structure in the Cosmic Microwave Background observed by PLANCK, specifically at 143 and 857 GHz, as well as with 100 μm observations from the IRIS survey. The analysis uses data in 13 2° × 2° sub-areas found in the IRSA database at IPAC. The results confirm what has been reported previously; nearby galactic HI features and high-frequency continuum sources believed to be cosmological are in fact clearly associated. While several attempts strongly suggest that the associations are statistically significant, the key to understanding the phenomenon lies in the fact that in any given area HI is associated with cirrus dust at certain HI velocities and with 143 GHz features at different velocities. At the same time, for the 13 sub-areas studied, there is very little overlap between the dust and 143 GHz features. The data do not imply that the HI itself gives rise to the high-frequency continuum emission. Rather, they appear to indicate undiagnosed brightness enhancements indirectly associated with the HI. If low density interstellar electrons concentrated into clumps, or observed in directions where their integrated line-of-sight column densities are greater than the background in a manner similar to the phenomena that give rise to structure in diffuse HI structure, they will profoundly affect attempts to create a foreground electron mask used for processing PLANCK as well as WMAP data.

We analyzed the orientation of the sample of ACO galaxy clusters. We examined the alignment in a subsample of 1056 galaxy structures taken from the Panko–Flin (2006) Catalog with known BM morphological types. We were looking for a correlation between the orientation of the cluster and the positions of neighboring clusters. The Binggeli effect (the excess of small values of the Δθ angles between the direction toward neighboring clusters and the cluster position angle) is observed, having a range up to about 45 h⁻¹ Mpc. The strongest effect was found for elongated BM type I clusters. This is probably connected with the origins of the supergiant galaxy and with cluster formation along a long filament or plane in a supercluster.

We present a model for the evolution of the galaxy ultraviolet (UV) luminosity function (LF) across cosmic time where star formation is linked to the assembly of dark matter halos under the assumption of a mass-dependent, but redshift-independent, efficiency. We introduce a new self-consistent treatment of the halo star formation history, which allows us to make predictions at z > 10 (lookback time ≲500 Myr), when growth is rapid. With a calibration at a single redshift to set the stellar-to-halo mass ratio, and no further degrees of freedom, our model captures the evolution of the UV LF over all available observations (0 ≲ z ≲ 10). The significant drop in luminosity density of currently detectable galaxies beyond z ∼ 8 is explained by a shift of star formation toward less massive, fainter galaxies. Assuming that star formation proceeds down to atomic cooling halos, we derive a reionization optical depth $\tau ={0.056}_{-0.010}^{+0.007},$ fully consistent with the latest Planck measurement, implying that the universe is fully reionized at $z={7.84}_{-0.98}^{+0.65}.$ In addition, our model naturally produces smoothly rising star formation histories for galaxies with L ≲ L_* in agreement with observations and hydrodynamical simulations. Before the epoch of reionization at z > 10 we predict the LF to remain well-described by a Schechter function, but with an increasingly steep faint-end slope (α ∼ −3.5 at z ∼ 16). Finally, we construct forecasts for surveys with James Webb Space Telescope (JWST) and Wide-field Infrared Survey Telescope (WFIRST) and predict that galaxies out to z ∼ 14 will be observed. Galaxies at z > 15 will likely be accessible to JWST and WFIRST only through the assistance of strong lensing magnification.

Understanding whether Helium can sediment to the core of galaxy clusters is important for a number of problems in cosmology and astrophysics. All current models addressing this question are one-dimensional and do not account for the fact that magnetic fields can effectively channel ions and electrons, leading to anisotropic transport of momentum, heat, and particle diffusion in the weakly collisional intracluster medium (ICM). This anisotropy can lead to a wide variety of instabilities, which could be relevant for understanding the dynamics of heterogeneous media. In this paper, we consider the radial temperature and composition profiles as obtained from a state-of-the-art Helium sedimentation model and analyze its stability properties. We find that the associated radial profiles are unstable to different kinds of instabilities depending on the magnetic field orientation at all radii. The fastest growing modes are usually related to generalizations of the magnetothermal instability (MTI) and the heat-flux-driven buoyancy instability which operate in heterogeneous media. We find that the effect of sedimentation is to increase (decrease) the predicted growth rates in the inner (outer) cluster region. The unstable modes grow quickly compared to the sedimentation timescale. This suggests that the composition gradients as inferred from sedimentation models, which do not fully account for the anisotropic character of the weakly collisional environment, might not be very robust. Our results emphasize the subtleties involved in understanding the gas dynamics of the ICM and argue for the need of a comprehensive approach to address the issue of Helium sedimentation beyond current models.

In this paper we study a key phase in the formation of massive galaxies: the transition of star-forming galaxies into massive (M_stars ∼ 10¹¹M_⊙), compact (r_e ∼ 1 kpc) quiescent galaxies, which takes place from z ∼ 3 to z ∼ 1.5. We use HST grism redshifts and extensive photometry in all five 3D-HST/CANDELS fields, more than doubling the area used previously for such studies, and combine these data with Keck MOSFIRE and NIRSPEC spectroscopy. We first confirm that a population of massive, compact, star-forming galaxies exists at z ≳ 2, using K-band spectroscopy of 25 of these objects at 2.0 < z < 2.5. They have a median [N ii]/Hα ratio of 0.6, are highly obscured with SFR(tot)/SFR(Hα) ∼10, and have a large range of observed line widths. We infer from the kinematics and spatial distribution of Hα that the galaxies have rotating disks of ionized gas that are a factor of ∼2 more extended than the stellar distribution. By combining measurements of individual galaxies, we find that the kinematics are consistent with a nearly Keplerian fall-off from V_rot ∼ 500 km s⁻¹ at 1 kpc to V_rot ∼ 250 km s⁻¹ at 7 kpc, and that the total mass out to this radius is dominated by the dense stellar component. Next, we study the size and mass evolution of the progenitors of compact massive galaxies. Even though individual galaxies may have had complex histories with periods of compaction and mergers, we show that the population of progenitors likely followed a simple inside-out growth track in the size–mass plane of ${\rm{\Delta }}\mathrm{log}{r}_{{\rm{e}}}\sim 0.3{\rm{\Delta }}\mathrm{log}{M}_{{\rm{stars}}}.$ This mode of growth gradually increases the stellar mass within a fixed physical radius, and galaxies quench when they reach a stellar density or velocity dispersion threshold. As shown in other studies, the mode of growth changes after quenching, as dry mergers take the galaxies on a relatively steep track in the size–mass plane.

The detection of an excess of emission at microwave frequencies with respect to the predicted free–free emission has been reported for several Galactic H ii regions. Here, we investigate the case of RCW 49, for which the Cosmic Background Imager tentatively (∼3σ) detected Anomalous Microwave Emission (AME) at 31 GHz on angular scales of 7'. Using the Australia Telescope Compact Array, we carried out a multi-frequency (5, 19, and 34 GHz) continuum study of the region, complemented by observations of the H109α radio recombination line. The analysis shows that: (1) the spatial correlation between the microwave and IR emission persists on angular scales from 34 to 04, although the degree of the correlation slightly decreases at higher frequencies and on smaller angular scales; (2) the spectral indices between 1.4 and 5 GHz are globally in agreement with optically thin free–free emission, however, ∼30% of these are positive and much greater than −0.1, consistent with a stellar wind scenario; and (3) no major evidence for inverted free–free radiation is found, indicating that this is likely not the cause of the Anomalous Emission in RCW 49. Although our results cannot rule out the spinning dust hypothesis to explain the tentative detection of AME in RCW 49, they emphasize the complexity of astronomical sources that are very well known and studied, such as H ii regions, and suggest that, at least in these objects, the reported excess of emission might be ascribed to alternative mechanisms such as stellar winds and shocks.

We present the initial results of our investigation of the star-forming complex W49, one of the youngest and most luminous massive star-forming regions in our Galaxy. We used Spitzer/Infrared Array Camera (IRAC) data to investigate massive star formation with the primary objective of locating a representative set of protostars and the clusters of young stars that are forming around them. We present our source catalog with the mosaics from the IRAC data. In this study we used a combination of IRAC, MIPS, Two Micron All Sky Survey, and UKIRT Deep Infrared Sky Survey (UKIDSS) data to identify and classify the young stellar objects (YSOs). We identified 232 Class 0/I YSOs, 907 Class II YSOs, and 74 transition disk candidate objects using color–color and color–magnitude diagrams. In addition, to understand the evolution of star formation in W49, we analyzed the distribution of YSOs in the region to identify clusters using a minimal spanning tree method. The fraction of YSOs that belong to clusters with ≥7 members is found to be 52% for a cutoff distance of 96'', and the ratio of Class II/I objects is 2.1. We compared the W49 region to the G305 and G333 star-forming regions and concluded that W49 has the richest population, with seven subclusters of YSOs.

We present a sample of 20 runaway M dwarf candidates (RdMs) within 1 kpc of the Sun whose Galactocentric (GC) velocities exceed 400 km s⁻¹. The candidates were selected from the Sloan Digital Sky Survey (SDSS) DR7 M Dwarf Catalog of West et al. Our RdMs have SDSS+USNO-B proper motions that are consistent with those recorded in the PPMXL, LSPM, and combined Wide-field Infrared Survey Explorer+SDSS+Two-micron All-sky Survey catalogs. Sixteen RdMs are classified as dwarfs, while the remaining four RdMs are subdwarfs. We model the Galactic potential using a bulge-disk-halo profile. Our fastest RdM, with a GC velocity of 658.5 ± 236.9 km s⁻¹, is a possible hypervelocity candidate, as it is unbound in 77% of our simulations. About half of our RdMs have kinematics that are consistent with ejection from the Galactic center. Seven of our RdMs have kinematics consistent with an ejection scenario from M31 or M32 to within 2σ, although our distance-limited survey makes such a realization unlikely. No more than four of our RdMs may have originated from the Leo stream. We propose that to within measurement errors, most of our bound RdMs are likely disk runaways or halo objects, and may have been accelerated through a series of multi-body interactions within the Galactic disk or possibly supernovae explosions.

At a projected distance of ∼26 pc from Sgr A*, the Arches cluster provides insight into star formation in the extreme Galactic center (GC) environment. Despite its importance, many key properties, such as the cluster's internal structure and orbital history, are not well known. We present an astrometric and photometric study of the outer region of the Arches cluster (R > 625) using Hubble Space Telescope WFC3IR. Using proper motions, we calculate membership probabilities for stars down to F153M = 20 mag (∼2.5 M_⊙) over a 120'' × 120'' field of view, an area 144 times larger than previous astrometric studies of the cluster. We construct the radial profile of the Arches to a radius of 75'' (∼3 pc at 8 kpc), which can be well described by a single power law. From this profile we place a 3σ lower limit of 2.8 pc on the observed tidal radius, which is larger than the predicted tidal radius (1–2.5 pc). Evidence of mass segregation is observed throughout the cluster, and no tidal tail structures are apparent along the orbital path. The absence of breaks in the profile suggests that the Arches has not likely experienced its closest approach to the GC between ∼0.2 and 1 Myr ago. If accurate, this constraint indicates that the cluster is on a prograde orbit and is located in front of the sky plane that intersects Sgr A*. However, further simulations of clusters in the GC potential are required to interpret the observed profile with more confidence.

In the current era of large surveys and massive data sets, autoclassification of astrophysical sources using intelligent algorithms is becoming increasingly important. In this paper we present the catalog of variable sources in the Third XMM-Newton Serendipitous Source catalog (3XMM) autoclassified using the Random Forest machine learning algorithm. We used a sample of manually classified variable sources from the second data release of the XMM-Newton catalogs (2XMMi-DR2) to train the classifier, obtaining an accuracy of ∼92%. We also evaluated the effectiveness of identifying spurious detections using a sample of spurious sources, achieving an accuracy of ∼95%. Manual investigation of a random sample of classified sources confirmed these accuracy levels and showed that the Random Forest machine learning algorithm is highly effective at automatically classifying 3XMM sources. Here we present the catalog of classified 3XMM variable sources. We also present three previously unidentified unusual sources that were flagged as outlier sources by the algorithm: a new candidate supergiant fast X-ray transient, a 400 s X-ray pulsar, and an eclipsing 5 hr binary system coincident with a known Cepheid.

Spectral energy distributions for 132 classical and type II Cepheids were searched for evidence of excess flux above the photospheric level in the mid-infrared. Eight of them were found to have unambiguously strong excess emission while a further 13 showed evidence of weak emission. The presence of emission appears to be unrelated to either the pulsational amplitude or the effective temperature while strong emission is limited to stars with periods longer than 11 days, with a single exception. For the stars with strong emission we attempted to fit the energy distribution with a stellar wind model. No acceptable fit could be found for silicate grains. With graphite or iron grains we could only obtain an acceptable fit if the maximum dust temperature was significantly lower than the condensation temperature. We conclude that the excess emission is not evidence of mass loss.

We present ultraviolet (UV) and optical photometry and spectra of the 1999aa-like supernova (SN) iPTF14bdn. The UV data were observed using the Swift Ultraviolet/Optical Telescope and constitute the first UV spectral series of a 1999aa-like SN. From the photometry, we measure Δm₁₅(B) = 0.84 ± 0.05 mag and blue UV colors at epochs earlier than −5 days. The spectra show that the early-time blue colors are the result of less absorption between 2800−3200 Å than is present in normal SNe Ia. Using model spectra fits of the data at −10 and +10 days, we identify the origin of this spectral feature to be a temperature effect in which doubly ionized iron group elements create an opacity "window." We determine that the detection of high temperatures and large quantities of iron group elements at early epochs imply the mixing of a high Ni mass into the outer layers of the SN ejecta. We also identify the source of the I-band secondary maximum in iPTF14bdn to be the decay of Fe iii to Fe ii, as is seen in normal SNe Ia.

We present a sample of 11 M31 Cepheids in stellar clusters, derived from the overlap of the Panchromatic Hubble Andromeda Treasury cluster catalog and the Pan-STARRS1 (PS1) disk Cepheid catalog. After identifying the PS1 Cepheids in the Hubble Space Telescope (HST) catalog, we calibrate the PS1 mean magnitudes using the higher resolution HST photometry, revealing up to 1 mag offsets due to crowding effects in the ground-based catalog. We measure ages of the clusters by performing single-age stellar population fits to their color–magnitude diagrams excluding their Cepheids. From these cluster age measurements, we derive an empirical period–age relation which agrees well with the existing literature values. By confirming this relation for M31 Cepheids, we justify its application in high-precision pointwise age estimation across M31.

The Extreme-ultraviolet Imaging Spectrometer (EIS) on the Hinode spacecraft observed flare footpoint regions coincident with a surge for an M3.7 flare observed on 2011 September 25 at N12 E33 in active region 11302. The flare was observed in spectral lines of O vi, Fe x, Fe xii, Fe xiv, Fe xv, Fe xvi, Fe xvii, Fe xxiii, and Fe xxiv. The EIS observations were made coincident with hard X-ray bursts observed by RHESSI. Overlays of the RHESSI images on the EIS raster images at different wavelengths show a spatial coincidence of features in the RHESSI images with the EIS upflow and downflow regions, as well as loop-top or near-loop-top regions. A complex array of phenomena were observed, including multiple evaporation regions and the surge, which was also observed by the Solar Dynamics Observatory/Atmospheric Imaging Assembly telescopes. The slit of the EIS spectrometer covered several flare footpoint regions from which evaporative upflows in Fe xxiii and Fe xxiv lines were observed with Doppler speeds greater than 500 km s⁻¹. For ions such as Fe xv both evaporative outflows (∼200 km s⁻¹) and downflows (∼30–50 km s⁻¹) were observed. Nonthermal motions from 120 to 300 km s⁻¹ were measured in flare lines. In the surge, Doppler speeds are found from about 0 to over 250 km s⁻¹ in lines from ions such as Fe xiv. The nonthermal motions could be due to multiple sources slightly Doppler-shifted from each other or turbulence in the evaporating plasma. We estimate the energetics of the hard X-ray burst and obtain a total flare energy in accelerated electrons of ≥7 × 10²⁸ erg. This is a lower limit because only an upper limit can be determined for the low-energy cutoff to the electron spectrum. We find that detailed modeling of this event would require a multithreaded model owing to its complexity.

Quasi-periodic propagating intensity disturbances have been observed in large coronal loops in extreme ultraviolet images over a decade, and are widely accepted to be slow magnetosonic waves. However, spectroscopic observations from Hinode/EIS revealed their association with persistent coronal upflows, making this interpretation debatable. We perform a 2.5D magnetohydrodynamic simulation to imitate the chromospheric evaporation and the following reflected patterns in a flare loop. Our model encompasses the corona, transition region, and chromosphere. We demonstrate that the quasi periodic propagating intensity variations captured by the synthesized Solar Dynamics Observatory/Atmospheric Imaging Assembly 131, 94 Å emission images match the previous observations well. With particle tracers in the simulation, we confirm that these quasi periodic propagating intensity variations consist of reflected slow mode waves and mass flows with an average speed of 310 km s⁻¹ in an 80 Mm length loop with an average temperature of 9 MK. With the synthesized Doppler shift velocity and intensity maps of the Solar and Heliospheric Observatory/Solar Ultraviolet Measurement of Emitted Radiation Fe xix line emission, we confirm that these reflected slow mode waves are propagating waves.

The O i 135.56 nm line is covered by NASA's Interface Region Imaging Spectrograph (IRIS) small explorer mission which studies how the solar atmosphere is energized. We study here the formation and diagnostic potential of this line by means of non-local thermodynamic equilibrium modeling employing both 1D semi-empirical and 3D radiation magnetohydrodynamic models. We study the basic formation mechanisms and derive a quintessential model atom that incorporates essential atomic physics for the formation of the O i 135.56 nm line. This atomic model has 16 levels and describes recombination cascades through highly excited levels by effective recombination rates. The ionization balance O i/O ii is set by the hydrogen ionization balance through charge exchange reactions. The emission in the O i 135.56 nm line is dominated by a recombination cascade and the line is optically thin. The Doppler shift of the maximum emission correlates strongly with the vertical velocity in its line forming region, which is typically located at 1.0–1.5 Mm height. The total intensity of the line emission is correlated with the square of the electron density. Since the O i 135.56 nm line is optically thin, the width of the emission line is a very good diagnostic of non-thermal velocities. We conclude that the O i 135.56 nm line is an excellent probe of the middle chromosphere, and compliments other powerful chromospheric diagnostics of IRIS such as the Mg ii h & k lines and the C ii lines around 133.5 nm.

The 3-part appearance of many coronal mass ejections (CMEs) arising from erupting filaments emerges from a large magnetic flux tube structure, consistent with the form of the erupting filament system. Other CMEs arising from erupting filaments lack a clear 3-part structure and reasons for this have not been researched in detail. This paper aims to further establish the link between CME structure and the structure of the erupting filament system and to investigate whether CMEs which lack a 3-part structure have different eruption characteristics. A survey is made of 221 near-limb filament eruptions observed from 2013 May 03 to 2014 June 30 by Extreme UltraViolet (EUV) imagers and coronagraphs. Ninety-two filament eruptions are associated with 3-part structured CMEs, 41 eruptions are associated with unstructured CMEs. The remaining 88 are categorized as failed eruptions. For 34% of the 3-part CMEs, processing applied to EUV images reveals the erupting front edge is a pre-existing loop structure surrounding the filament, which subsequently erupts with the filament to form the leading bright front edge of the CME. This connection is confirmed by a flux-rope density model. Furthermore, the unstructured CMEs have a narrower distribution of mass compared to structured CMEs, with total mass comparable to the mass of 3-part CME cores. This study supports the interpretation of 3-part CME leading fronts as the outer boundaries of a large pre-existing flux tube. Unstructured (non 3-part) CMEs are a different family to structured CMEs, arising from the eruption of filaments which are compact flux tubes in the absence of a large system of enclosing closed field.

After two Atacama Large Millimeter/submillimeter Array (ALMA) observing cycles, only a handful of [C ii] 158 μm emission line searches in z > 6 galaxies have reported a positive detection, questioning the applicability of the local [C ii]–star formation rate (SFR) relation to high-z systems. To investigate this issue we use the Vallini et al. (V13) model,based on high-resolution, radiative transfer cosmological simulations to predict the [C ii] emission from the interstellar medium of a z ≈ 7 (halo mass M_h = 1.17 × 10¹¹M_⊙) galaxy. We improve the V13 model by including (a) a physically motivated metallicity (Z) distribution of the gas, (b) the contribution of photodissociation regions (PDRs), and (c) the effects of cosmic microwave background (CMB) on the [C ii] line luminosity. We study the relative contribution of diffuse neutral gas to the total [C ii] emission (F_diff/F_tot) for different SFR and Z values. We find that the [C ii] emission arises predominantly from PDRs: regardless of the galaxy properties, F_diff/F_tot ≤ 10%, since at these early epochs the CMB temperature approaches the spin temperature of the [C ii] transition in the cold neutral medium (T_CMB ∼ ${T}_{s}^{{\rm{CNM}}}$ ∼ 20 K). Our model predicts a high-z [C ii]–SFR relation, consistent with observations of local dwarf galaxies (0.02 < Z/Z_⊙ < 0.5). The [C ii] deficit suggested by actual data (L_Cii < 2.0 × 10⁷L_⊙ in BDF3299 at z ≈ 7.1) if confirmed by deeper ALMA observations, can be ascribed to negative stellar feedback disrupting molecular clouds around star formation sites. The deviation from the local [C ii]–SFR would then imply a modified Kennicutt–Schmidt relation in z > 6 galaxies. Alternatively/in addition, the deficit might be explained by low gas metallicities (Z < 0.1 Z_⊙).

We present a detailed gravitational lens model of the galaxy cluster RCS2 J232727.6-020437. Due to cosmological dimming of cluster members and intra-cluster light, its high redshift (z = 0.6986) makes it ideal for studying high-redshift background galaxies. Using new Advanced Camera for Surveys and WFC3/IR Hubble Space Telescope data, we identify 16 new multiple images belonging to 6 distinct source galaxies. From Multi-object Spectrometer for Infrared Exploration (MOSFIRE) follow-up, we identify a strong emission line in the spectrum of one multiple image, measuring the redshift of that system to z = 2.083. With a highly magnified (μ ≳ 3) source plane area of ∼0.9 arcmin² at z = 7, RCS2 J232727.6-020437 has a lensing efficiency comparable to the Hubble Frontier Fields clusters. We discover four highly magnified z ∼ 7 candidate Lyman-break galaxies behind the cluster. Correcting for magnification, we find that all four galaxies are fainter than $0.5{L}_{\star }.$ One candidate is detected at >10σ in both Spitzer/IRAC [3.6] and [4.5] channels. A spectroscopic follow-up with MOSFIRE does not result in the detection of the Lyα emission line from any of the four candidates. From the MOSFIRE spectra, we place median upper limits on the Lyα flux of $(3-11)\times {10}^{-19}\;\mathrm{erg}\;{{\rm{s}}}^{-1}\;{\mathrm{cm}}^{-2}$ (5σ).

We present Sedonu, a new open source, steady-state, special relativistic Monte Carlo (MC) neutrino transport code, available at bitbucket.org/srichers/sedonu. The code calculates the energy- and angle-dependent neutrino distribution function on fluid backgrounds of any number of spatial dimensions, calculates the rates of change of fluid internal energy and electron fraction, and solves for the equilibrium fluid temperature and electron fraction. We apply this method to snapshots from two-dimensional simulations of accretion disks left behind by binary neutron star mergers, varying the input physics and comparing to the results obtained with a leakage scheme for the cases of a central black hole and a central hypermassive neutron star. Neutrinos are guided away from the densest regions of the disk and escape preferentially around 45° from the equatorial plane. Neutrino heating is strengthened by MC transport a few scale heights above the disk midplane near the innermost stable circular orbit, potentially leading to a stronger neutrino-driven wind. Neutrino cooling in the dense midplane of the disk is stronger when using MC transport, leading to a globally higher cooling rate by a factor of a few and a larger leptonization rate by an order of magnitude. We calculate neutrino pair annihilation rates and estimate that an energy of 2.8 × 10⁴⁶ erg is deposited within 45° of the symmetry axis over 300 ms when a central BH is present. Similarly, 1.9 × 10⁴⁸ erg is deposited over 3 s when an HMNS sits at the center, but neither estimate is likely to be sufficient to drive a gamma-ray burst jet.

We describe directed searches for continuous gravitational waves (GWs) in data from the sixth Laser Interferometer Gravitational-wave Observatory (LIGO) science data run. The targets were nine young supernova remnants not associated with pulsars; eight of the remnants are associated with non-pulsing suspected neutron stars. One target's parameters are uncertain enough to warrant two searches, for a total of 10. Each search covered a broad band of frequencies and first and second frequency derivatives for a fixed sky direction. The searches coherently integrated data from the two LIGO interferometers over time spans from 5.3–25.3 days using the matched-filtering ${\mathcal{F}}$-statistic. We found no evidence of GW signals. We set 95% confidence upper limits as strong (low) as 4 × 10⁻²⁵ on intrinsic strain, 2 × 10⁻⁷ on fiducial ellipticity, and 4 × 10⁻⁵ on r-mode amplitude. These beat the indirect limits from energy conservation and are within the range of theoretical predictions for neutron-star ellipticities and r-mode amplitudes.

Rotation evolution of late-type stars is dominated by magnetic braking and the underlying factors that control this angular momentum loss are important for the study of stellar spin-down. In this work, we study angular momentum loss as a function of two different aspects of magnetic activity using a calibrated Alfvén wave-driven magnetohydrodynamic wind model: the strengths of magnetic spots and their distribution in latitude. By driving the model using solar and modified solar surface magnetograms, we show that the topology of the field arising from the net interaction of both small-scale and large-scale field is important for spin-down rates and that angular momentum loss is not a simple function of large scale magnetic field strength. We find that changing the latitude of magnetic spots can modify mass and angular momentum loss rates by a factor of two. The general effect that causes these differences is the closing down of large-scale open field at mid- and high-latitudes by the addition of the small-scale field. These effects might give rise to modulation of mass and angular momentum loss through stellar cycles, and present a problem for ab initio attempts to predict stellar spin-down based on wind models. For all the magnetogram cases considered here, from dipoles to various spotted distributions, we find that angular momentum loss is dominated by the mass loss at mid-latitudes. The spin-down torque applied by magnetized winds therefore acts at specific latitudes and is not evenly distributed over the stellar surface, though this aspect is unlikely to be important for understanding spin-down and surface flows on stars.

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at λ = 0.9, 2.8, 8.0, 9.8 mm for DoAr 25 and at λ = 1.3, 2.8, 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, β, is consistent with the observed dust emission in both disks, with low-β in the inner disk and high-β in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (a_max), we constrain radial variations of a_max in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our a_max(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions.

We present the first results from our survey of the star-forming complex W3, combining VRI photometry with multiobject spectroscopy to identify and characterize the high-mass stellar population across the region. With 79 new spectral classifications, we bring the total number of spectroscopically confirmed O- and B-type stars in W3 to 105. We find that the high-mass slope of the mass function in W3 is consistent with a Salpeter IMF, and that the extinction toward the region is best characterized by an R_V of approximately 3.6. B-type stars are found to be more widely dispersed across the W3 giant molecular cloud (GMC) than previously realized: they are not confined to the high-density layer (HDL) created by the expansion of the neighboring W4 H ii region into the GMC. This broader B-type population suggests that star formation in W3 began spontaneously up to 8–10 Myr ago, although at a lower level than the more recent star formation episodes in the HDL. In addition, we describe a method of optimizing sky subtraction for fiber spectra in regions of strong and spatially variable nebular emission.

A joint analysis is done of the radio and X-ray observations of SN 1993J. It is argued that neither synchrotron cooling behind the forward shock nor thermal cooling behind the reverse shock is supported by observations. In order for adiabatic models to be consistent, a reinterpretation of the radius of the spatially resolved very long baseline interferometry-source (VLBI) is needed during the first few hundred days. Instead of reflecting the position of the forward shock, it is then associated with the expansion of the Rayleigh–Taylor unstable region emanating from the contact discontinuity. Although observations imply a constant ratio between the energy densities in magnetic fields and relativistic electrons, they do not appear to scale individually with the thermal energy density behind the forward shock; rather, in adiabatic models, the evolution of the magnetic field strength is best understood as scaling inversely with the supernova radius.

We present the discovery of three new Milky Way satellites from our search for compact stellar overdensities in the photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS 1, or PS1) 3π survey. The first satellite, Laevens 3, is located at a heliocentric distance of d = 67 ± 3 kpc. With a total magnitude of M_V = −4.4 ± 0.3 and a half-light radius of r_h = 7 ± 2 pc, its properties resemble those of outer halo globular clusters. The second system, Draco II/Laevens 4, is a closer and fainter satellite (d ∼ 20 kpc, M_V = −2.9 ± 0.8), whose uncertain size (${r}_{h}={19}_{-6}^{+8}\;\mathrm{pc}$) renders its classification difficult without kinematic information; it could either be a faint and extended globular cluster or a faint and compact dwarf galaxy. The third satellite, Sagittarius II/Laevens 5 (Sgr II), has an ambiguous nature, as it is either the most compact dwarf galaxy or the most extended globular cluster in its luminosity range (${r}_{h}={37}_{-8}^{+9}\;\mathrm{pc}$ and M_V = −5.2 ± 0.4). At a heliocentric distance of 67 ± 5 kpc, this satellite lies intriguingly close to the expected location of the trailing arm of the Sagittarius stellar stream behind the Sagittarius dwarf spheroidal galaxy (Sgr dSph). If confirmed through spectroscopic follow up, this connection would locate this part of the trailing arm of the Sagittarius stellar stream that has so far gone undetected. It would further suggest that Sgr II was brought into the Milky Way halo as a satellite of the Sgr dSph.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 870 μm (345 GHz) data for 49 high-redshift (0.47 < z < 2.85), luminous ($11.7\lt \mathrm{log}({L}_{{\rm{bol}}}/{L}_{\odot })\lt 14.2$) radio-powerful active galactic nuclei (AGNs), obtained to constrain cool dust emission from starbursts concurrent with highly obscured radiative-mode black hole (BH) accretion in massive galaxies that possess a small radio jet. The sample was selected from the Wide-field Infrared Survey Explorer with extremely steep (red) mid-infrared colors and with compact radio emission from NVSS/FIRST. Twenty-six sources are detected at 870 μm, and we find that the sample has large mid- to far-infrared luminosity ratios, consistent with a dominant and highly obscured quasar. The rest-frame 3 GHz radio powers are $24.7\lt \mathrm{log}({P}_{\text{3.0 GHz}}/{{\rm{W}}\;\mathrm{Hz}}^{-1})\lt 27.3,$ and all sources are radio-intermediate or radio-loud. BH mass estimates are 7.7 < log(M_BH/M_⊙) < 10.2. The rest-frame 1–5 μm spectral energy distributions are very similar to the "Hot DOGs" (hot dust-obscured galaxies), and steeper (redder) than almost any other known extragalactic sources. ISM masses estimated for the ALMA-detected sources are 9.9 < log (M_ISM/M_⊙) < 11.75 assuming a dust temperature of 30 K. The cool dust emission is consistent with star formation rates reaching several thousand M_⊙ yr⁻¹, depending on the assumed dust temperature, but we cannot rule out the alternative that the AGN powers all the emission in some cases. Our best constrained source has radiative transfer solutions with approximately equal contributions from an obscured AGN and a young (10–15 Myr) compact starburst.

We present a study exploring the nature and properties of the circumgalactic medium (CGM) and its connection to the atomic gas content in the interstellar medium (ISM) of galaxies as traced by the H i 21 cm line. Our sample includes 45 low-z (0.026–0.049) galaxies from the GALEX Arecibo SDSS Survey (Galaxy Evolution Explorer/Arecibo/Sloan Digital Sky Survey). Their CGM was probed via absorption in the spectra of background quasi-stellar objects at impact parameters of 63–231 kpc. The spectra were obtained with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. We detected neutral hydrogen (Lyα absorption lines) in the CGM of 92% of the galaxies. We find that the radial profile of the CGM as traced by the Lyα equivalent width can be fit as an exponential with a scale length of roughly the virial radius of the dark matter halo. We found no correlation between the orientation of the sightline relative to the galaxy's major axis and the Lyα equivalent width. The velocity spread of the circumgalactic gas is consistent with that seen in the atomic gas in the ISM. We find a strong correlation (99.8% confidence) between the gas fraction (M(H i)/M_⋆) and the impact-parameter-corrected Lyα equivalent width. This is stronger than the analogous correlation between corrected Lyα equivalent width and specific star formation rate (SFR)/M_⋆ (97.5% confidence). These results imply a physical connection between the H i disk and the CGM, which is on scales an order of magnitude larger. This is consistent with the picture in which the H i disk is nourished by accretion of gas from the CGM.

Many parameters constraining the spectral appearance of exoplanets are still poorly understood. We therefore study the properties of irradiated exoplanet atmospheres over a wide parameter range including metallicity, C/O ratio, and host spectral type. We calculate a grid of 1D radiative-convective atmospheres and emission spectra. We perform the calculations with our new Pressure–Temperature Iterator and Spectral Emission Calculator for Planetary Atmospheres (PETIT) code, assuming chemical equilibrium. The atmospheric structures and spectra are made available online. We find that atmospheres of planets with C/O ratios ∼1 and ${T}_{{\rm{eff}}}$ ≳ 1500 K can exhibit inversions due to heating by the alkalis because the main coolants CH₄, H₂O, and HCN are depleted. Therefore, temperature inversions possibly occur without the presence of additional absorbers like TiO and VO. At low temperatures we find that the pressure level of the photosphere strongly influences whether the atmospheric opacity is dominated by either water (for low C/O) or methane (for high C/O), or both (regardless of the C/O). For hot, carbon-rich objects this pressure level governs whether the atmosphere is dominated by methane or HCN. Further we find that host stars of late spectral type lead to planetary atmospheres which have shallower, more isothermal temperature profiles. In agreement with prior work we find that for planets with ${T}_{{\rm{eff}}}\lt 1750$ K the transition between water or methane dominated spectra occurs at C/O ∼ 0.7, instead of ∼1, because condensation preferentially removes oxygen.

We present 50 nights of polarimetric observations of HD 189733 in the B band using the POLISH2 aperture-integrated polarimeter at the Lick Observatory Shane 3-m telescope. This instrument, commissioned in 2011, is designed to search for Rayleigh scattering from short-period exoplanets due to the polarized nature of scattered light. Since these planets are spatially unresolvable from their host stars, the relative contribution of the planet-to-total system polarization is expected to vary with an amplitude of the order of 10 parts per million (ppm) over the course of the orbit. Non-zero and also variable at the 10 ppm level, the inherent polarization of the Lick 3-m telescope limits the accuracy of our measurements and currently inhibits conclusive detection of scattered light from this exoplanet. However, the amplitude of observed variability conservatively sets a 99.7% confidence upper limit to the planet-induced polarization of the system of 60 ppm in the B band, which is consistent with a previous upper limit from the POLISH instrument at the Palomar Observatory 5-m telescope. A physically motivated Rayleigh scattering model, which includes the depolarizing effects of multiple scattering, is used to conservatively set a 99.7% confidence upper limit to the geometric albedo of HD 189733b of A_g < 0.40. This value is consistent with the value ${A}_{g}=0.226\pm 0.091$ derived from occultation observations with Hubble Space Telescope STIS, but it is inconsistent with the large ${A}_{g}=0.61\pm 0.12$ albedo reported by Berdyugina et al.

The magnetic field measured by Voyager 1 prior to its heliocliff encounter on 2012.65 showed an unexpectedly complex transition from the primarily azimuthal inner-heliosheath field to the draped interstellar field tilted by some 20° to the nominal azimuthal direction. Most prominent were two regions of enhanced magnetic field strength depleted in energetic charged particles of heliospheric origin. These regions were interpreted as magnetic flux tubes connected to the outer heliosheath that provided a path for the particles to escape. Despite large increases in strength, the field's direction did not change appreciably at the boundaries of these flux tubes. Rather, the field's direction changed gradually over several months prior to the heliocliff crossing. It is shown theoretically that the heliopause, as a pressure equilibrium layer, can become unstable to interchange of magnetic fields between the inner and the outer heliosheaths. The curvature of magnetic field lines and the anti-sunward gradient in plasma kinetic pressure provide conditions favorable for an interchange. Magnetic shear between the heliosheath and the interstellar fields reduces the growth rates, but does not fully stabilize the heliopause against perturbations propagating in the latitudinal direction. The instability could create a transition layer permeated by magnetic flux tubes, oriented parallel to each other and alternately connected to the heliosheath or the interstellar regions.

This is the second paper in a series where we build a self-consistent model to simulate the mass-loss process of a close-orbit magnetized giant exoplanet, so-called hot Jupiter (HJ). In this paper we generalize the hydrodynamic (HD) model of an HJ's expanding hydrogen atmosphere, proposed in the first paper, to include the effects of intrinsic planetary magnetic field. The proposed self-consistent axisymmetric 2D magnetohydrodynamics model incorporates radiative heating and ionization of the atmospheric gas, basic hydrogen chemistry for the appropriate account of major species composing HJ's upper atmosphere and related radiative energy deposition, and ${{\rm{H}}}_{3}^{+}$ and Lyα cooling processes. The model also takes into account a realistic solar-type X-ray/EUV spectrum for calculation of intensity and column density distribution of the radiative energy input, as well as gravitational and rotational forces acting in a tidally locked planet–star system. An interaction between the expanding atmospheric plasma and an intrinsic planetary magnetic dipole field leads to the formation of a current-carrying magnetodisk that plays an important role for topology and scaling of the planetary magnetosphere. A cyclic character of the magnetodisk behavior, composed of consequent phases of the disk formation followed by the magnetic reconnection with the ejection of a ring-type plasmoid, has been discovered and investigated. We found that the mass-loss rate of an HD 209458b analog planet is weakly affected by the equatorial surface field <0.3 G, but is suppressed by an order of magnitude at the field of 1 G.

We perform a multi-wavelength polarimetric study of the quasar CTA 102 during an extraordinarily bright γ-ray outburst detected by the Fermi Large Area Telescope in 2012 September–October when the source reached a flux of F_{>100 MeV} = 5.2 ± 0.4 × 10⁻⁶ photons cm⁻² s⁻¹. At the same time, the source displayed an unprecedented optical and near-infrared (near-IR) outburst. We study the evolution of the parsec-scale jet with ultra-high angular resolution through a sequence of 80 total and polarized intensity Very Long Baseline Array images at 43 GHz, covering the observing period from 2007 June to 2014 June. We find that the γ-ray outburst is coincident with flares at all the other frequencies and is related to the passage of a new superluminal knot through the radio core. The powerful γ-ray emission is associated with a change in direction of the jet, which became oriented more closely to our line of sight (θ ∼ 12) during the ejection of the knot and the γ-ray outburst. During the flare, the optical polarized emission displays intra-day variability and a clear clockwise rotation of electric vector position angles (EVPAs), which we associate with the path followed by the knot as it moves along helical magnetic field lines, although a random walk of the EVPA caused by a turbulent magnetic field cannot be ruled out. We locate the γ-ray outburst a short distance downstream of the radio core, parsecs from the black hole. This suggests that synchrotron self-Compton scattering of NIR to ultraviolet photons is the probable mechanism for the γ-ray production.

Nonhelical shear dynamos are studied with a particular focus on the possibility of coherent dynamo action. The primary results—serving as a follow up to the results of Squire & Bhattacharjee—pertain to the "magnetic shear-current effect" as a viable mechanism to drive large-scale magnetic field generation. This effect raises the interesting possibility that the saturated state of the small-scale dynamo could drive large-scale dynamo action, and is likely to be important in the unstratified regions of accretion disk turbulence. In this paper, the effect is studied at low Reynolds numbers, removing the complications of small-scale dynamo excitation and aiding analysis by enabling the use of quasi-linear statistical simulation methods. In addition to the magnetically driven dynamo, new results on the kinematic nonhelical shear dynamo are presented. These illustrate the relationship between coherent and incoherent driving in such dynamos, demonstrating the importance of rotation in determining the relative dominance of each mechanism.

Calibrating the photometric redshifts of ≳10⁹ galaxies for upcoming weak lensing cosmology experiments is a major challenge for the astrophysics community. The path to obtaining the required spectroscopic redshifts for training and calibration is daunting, given the anticipated depths of the surveys and the difficulty in obtaining secure redshifts for some faint galaxy populations. Here we present an analysis of the problem based on the self-organizing map, a method of mapping the distribution of data in a high-dimensional space and projecting it onto a lower-dimensional representation. We apply this method to existing photometric data from the COSMOS survey selected to approximate the anticipated Euclid weak lensing sample, enabling us to robustly map the empirical distribution of galaxies in the multidimensional color space defined by the expected Euclid filters. Mapping this multicolor distribution lets us determine where—in galaxy color space—redshifts from current spectroscopic surveys exist and where they are systematically missing. Crucially, the method lets us determine whether a spectroscopic training sample is representative of the full photometric space occupied by the galaxies in a survey. We explore optimal sampling techniques and estimate the additional spectroscopy needed to map out the color–redshift relation, finding that sampling the galaxy distribution in color space in a systematic way can efficiently meet the calibration requirements. While the analysis presented here focuses on the Euclid survey, similar analysis can be applied to other surveys facing the same calibration challenge, such as DES, LSST, and WFIRST.

SCORCH (Simulations and Constructions of the Reionization of Cosmic Hydrogen) is a new project to study the Epoch of Reionization (EoR). In this first paper, we probe the connection between observed high-redshift galaxies and simulated dark matter halos to better understand the primary source of ionizing radiation. High-resolution N-body simulations are run to quantify the abundance of dark matter halos as a function of mass M, accretion rate $\dot{M},$ and redshift z. A new fit for the halo mass function dn/dM is ≈20% more accurate at the high-mass end. A novel approach is used to fit the halo accretion rate function ${dn}/d\dot{M}$ in terms of the halo mass function. Abundance matching against the observed galaxy luminosity function is used to estimate the luminosity–mass relation and the luminosity–accretion-rate relation. The inferred star formation efficiency is not monotonic with M nor $\dot{M},$ but reaches a maximum value at a characteristic mass $\sim 2\times {10}^{11}\ {M}_{\odot }$ and a characteristic accretion rate $\sim 6\times {10}^{2}\ {M}_{\odot }\;{{\rm{yr}}}^{-1}$ at z ≈ 6. We find a universal EoR luminosity–accretion-rate relation and construct a fiducial model for the galaxy luminosity function. The Schechter parameters evolve such that ${\phi }_{\star }$ decreases, ${M}_{\star }$ is fainter, and α is steeper at higher redshifts. We forecast for the upcoming James Webb Space Telescope and show that with apparent magnitude limit ${m}_{{\rm{AB}}}\approx 31$ (32), it can observe $\gtrsim 11$ (24) unlensed galaxies per square degree per unit redshift at least down to ${M}_{\star }$ at z ≲ 13 (14).

Radioactive nuclei play an important role in planetary evolution by providing an internal heat source, which affects planetary structure and helps facilitate plate tectonics. A minimum level of nuclear activity is thought to be necessary—but not sufficient—for planets to be habitable. Extending previous work that focused on short-lived nuclei, this paper considers the delivery of long-lived radioactive nuclei to circumstellar disks in star forming regions. Although the long-lived nuclear species are always present, their abundances can be enhanced through multiple mechanisms. Most stars form in embedded cluster environments, so that disks can be enriched directly by intercepting ejecta from supernovae within the birth clusters. In addition, molecular clouds often provide multiple episodes of star formation, so that nuclear abundances can accumulate within the cloud; subsequent generations of stars can thus receive elevated levels of radioactive nuclei through this distributed enrichment scenario. This paper calculates the distribution of additional enrichment for ⁴⁰K, the most abundant of the long-lived radioactive nuclei. We find that distributed enrichment is more effective than direct enrichment. For the latter mechanism, ideal conditions lead to about 1 in 200 solar systems being directly enriched in ⁴⁰K at the level inferred for the early solar nebula (thereby doubling the abundance). For distributed enrichment from adjacent clusters, about 1 in 80 solar systems are enriched at the same level. Distributed enrichment over the entire molecular cloud is more uncertain, but can be even more effective.

Many works have attempted to investigate the astrophysical origin of the neutron-capture elements in the metal-poor star HD 140283. However, no definite conclusions have been drawn. In this work, using the abundance-decomposed approach, we find that the metal-poor star HD 140283 is a weak r-process star. Although this star is a weak r-process star, its Ba abundance mainly originates from the main r-process. This is the reason that the ratio [Ba/Eu]$=\;-0.58\pm 0.15$ for HD 140283 is close to the ratio of the main r-process. Based on the comparison of the abundances in the six-weak r-process stars, we find that their element abundances possess a robust nature. On the other hand, we find that the robust nature of the abundance of the extreme main r-process stars ([r/Fe] $\geqslant$ 1.5) can be extended to the lighter neutron-capture elements. Furthermore, the abundance characteristics of the weak r-process and main r-process are investigated. The abundance robustness of the two category r-process stars could be used as the constraint of the r-process theory and could be used to investigate the astrophysical origins of the elements in the metal-poor stars and population I stars.

In this work, we updated the catalog of Galactic Cepheids with 24 μm photometry by cross-matching the positions of known Galactic Cepheids to the recently released MIPSGAL point source catalog. We have added 36 new sources featuring MIPSGAL photometry in our analysis, thus increasing the existing sample to 65. Six different sources of compiled Cepheid distances were used to establish a 24 μm period–luminosity (P–L) relation. Our recommended 24 μm P–L relation is ${M}_{24\mu {\rm{m}}}=-3.18(\pm 0.10)\mathrm{log}P-2.46(\pm 0.10),$ with an estimated intrinsic dispersion of 0.20 mag, and is derived from 58 Cepheids exhibiting distances based on a calibrated Wesenheit function. The slopes of the P–L relations were steepest when tied solely to the 10 Cepheids exhibiting trigonometric parallaxes from the Hubble Space Telescope and Hipparcos. Statistical tests suggest that these P–L relations are significantly different from those associated with other methods of distance determination, and simulations indicate that difference may arise from the small sample size.

We have observed and spatially resolved a set of seven A-type stars in the nearby Ursa Major moving group with the Classic, CLIMB, and PAVO beam combiners on the Center for High Angular Resolution Astronomy Array. At least four of these stars have large rotational velocities ($v\mathrm{sin}i\;\gtrsim$ 170 $\mathrm{km}\;{{\rm{s}}}^{-1}$) and are expected to be oblate. These interferometric measurements, the stars' observed photometric energy distributions, and $v\mathrm{sin}i$ values are used to computationally construct model oblate stars from which stellar properties (inclination, rotational velocity, and the radius and effective temperature as a function of latitude, etc.) are determined. The results are compared with MESA stellar evolution models to determine masses and ages. The value of this new technique is that it enables the estimation of the fundamental properties of rapidly rotating stars without the need to fully image the star. It can thus be applied to stars with sizes comparable to the interferometric resolution limit as opposed to those that are several times larger than the limit. Under the assumption of coevality, the spread in ages can be used as a test of both the prescription presented here and the MESA evolutionary code for rapidly rotating stars. With our validated technique, we combine these age estimates and determine the age of the moving group to be 414 ± 23 Myr, which is consistent with, but much more precise than previous estimates.

We have explored the relationship between hard X-ray (HXR) emissions and Doppler velocities caused by the chromospheric evaporation in two X1.6 class solar flares on 2014 September 10 and October 22, respectively. Both events display double ribbons and the Interface Region Imaging Spectrograph slit is fixed on one of their ribbons from the flare onset. The explosive evaporations are detected in these two flares. The coronal line of Fe xxi 1354.09 Å shows blueshifts, but the chromospheric line of C i 1354.29 Å shows redshifts during the impulsive phase. The chromospheric evaporation tends to appear at the front of the flare ribbon. Both Fe xxi and C i display their Doppler velocities with an "increase-peak-decrease" pattern that is well related to the "rising-maximum-decay" phase of HXR emissions. Such anti-correlation between HXR emissions and Fe xxi Doppler shifts and correlation with C i Doppler shifts indicate the electron-driven evaporation in these two flares.

We report on an eruption involving a relatively rare filament–filament interaction on 2013 June 21, observed by SDO and STEREO-B. The two filaments were separated in height with a "double-decker" configuration. The eruption of the lower filament began simultaneously with a descent of the upper filament, resulting in a convergence and direct interaction of the two filaments. The interaction was accompanied by the heating of surrounding plasma and an apparent crossing of a loop-like structure through the upper filament. The subsequent coalescence of the filaments drove a bright front ahead of the erupting structures. The whole process was associated with a C3.0 flare followed immediately by an M2.9 flare. Shrinking loops and descending dark voids were observed during the M2.9 flare at different locations above a C-shaped flare arcade as part of the energy release, giving us unique insight into the flare dynamics.

We present a comprehensive method to analyze small-scale heating events in detail in a 3D magnetohydrodynamics simulation for quiet-Sun conditions. The method determines the number, volume, and some general geometric properties of the small-scale heating events at different instants in a simulation with a volume of 16 × 8 × 16 Mm³, spanning from the top of the convection zone to the corona. We found that there are about 10⁴ small-scale heating events at any instant above the simulated area of 128 Mm². They occur mainly at heights between 1.5 and 3.0 Mm. We determine the average value of their projected vertical extent, which ranges from 375 to 519 km over time, and we show that height, volume, and energy distribution of the events at any instant resemble power laws. Finally, we demonstrate that larger heating structures are a combination of much smaller heating events and that small-scale heating events dissipate enough energy to maintain the coronal energetic balance at any instant.

We present new integral field spectroscopy of the gravitationally lensed broad absorption line (BAL) quasar H1413+117, covering the ultraviolet restframe spectral range. We observe strong microlensing signatures in lensed image D, and we use this microlensing to simultaneously constrain both the broad emission and broad absorption line gas. The wavelength independence of image D magnifications across the broad emission lines (BELs) indicates a lower limit on the broad emission line region (BELR) size equal to the Einstein radius (ER) of the system: ≳11 ${(\langle M\rangle /{M}_{\odot })}^{0.5}$ lt-day for a lens redshift of 1.4 and ≳15 ${(\langle M\rangle /{M}_{\odot })}^{0.5}$ lt-day for z_L = 0.94. Lensing simulations verify that the observed wavelength independence is very unlikely for BELRs with significant velocity stratification at size scales below an ER. We perform spectral decomposition to derive the intrinsic BEL and continuum spectrum, subject to BAL absorption. We reconstruct the intrinsic BAL absorption profile, whose features allow us to constrain outflow kinematics in the context of a disk-wind model. We find a very sharp, blueshifted onset of absorption of 1500 km s⁻¹ in both C iv and N v, which may correspond to an inner edge of a disk-wind's radial outflow. The lower ionization Si iv and Al iii have higher-velocity absorption onsets, consistent with a decreasing ionization parameter with radius in an accelerating outflow. There is evidence of strong absorption in the BEL component, which indicates a high covering factor for absorption over two orders of magnitude in outflow radius.

The Fermi satellite detected GeV flashes from cosmological gamma-ray bursts (GRBs). In two GRBs, an optical counterpart of the GeV flash was detected. Such flashes are predicted by the model of a blast wave running into a medium loaded with copious ${e}^{\pm }$ pairs. Here we examine a sample of seven bursts with the best GeV+optical data and further test the model. We find that the observed light curves are in agreement with the theoretical predictions, which allows us to measure three parameters for each burst: the Lorentz factor of the explosion, its isotropic kinetic energy, and the external density. With the possible exception of GRB 090510 (the only short burst in the sample), the ambient medium is consistent with a wind from a Wolf–Rayet progenitor. The wind density parameter $A=\rho {r}^{2}$ varies in the sample around 10¹¹ g cm⁻¹. The initial Lorentz factor of the blast wave varies from 200 to 540, and correlates with the burst luminosity. Radiative efficiency of the prompt emission varies between 0.1 and 0.8. For the two bursts with a detected optical flash, GRB 120711A and GRB 130427A, we also estimate the magnetization of the external blast wave. Remarkably, despite its small number of free parameters, the model reproduces the entire optical light curve of GRB 120711A (with its sharp peak, fast decay, plateau, and break) as well as the GeV data. The spectra of GeV flashes are predicted to extend above 0.1 TeV, where they can be detected by ground-based Cherenkov telescopes.

We use two-dimensional relativistic hydrodynamic numerical calculations to show that highly collimated relativistic jets can be produced in neutron star merger models of short-duration gamma-ray bursts (GRBs) without the need for a highly directed engine or a large net magnetic flux. Even a hydrodynamic engine generating a very wide sustained outflow on small scales can, in principle, produce a highly collimated relativistic jet, facilitated by a dense surrounding medium that provides a cocoon surrounding the jet core. An oblate geometry to the surrounding gas significantly enhances the collimation process. Previous numerical simulations have shown that the merger of two neutron stars produces an oblate, expanding cloud of dynamical ejecta. We show that this gas can efficiently collimate the central engine power much like the surrounding star does in long-duration GRB models. For typical short-duration GRB central engine parameters, we find jets with opening angles of an order of 10° in which a large fraction of the total outflow power of the central engine resides in highly relativistic material. These results predict large differences in the opening angles of outflows from binary neutron star mergers versus neutron star–black hole mergers.

We present high-precision timing observations spanning up to nine years for 37 millisecond pulsars monitored with the Green Bank and Arecibo radio telescopes as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We describe the observational and instrumental setups used to collect the data, and methodology applied for calculating pulse times of arrival; these include novel methods for measuring instrumental offsets and characterizing low signal-to-noise ratio timing results. The time of arrival data are fit to a physical timing model for each source, including terms that characterize time-variable dispersion measure and frequency-dependent pulse shape evolution. In conjunction with the timing model fit, we have performed a Bayesian analysis of a parameterized timing noise model for each source, and detect evidence for excess low-frequency, or "red," timing noise in 10 of the pulsars. For 5 of these cases this is likely due to interstellar medium propagation effects rather than intrisic spin variations. Subsequent papers in this series will present further analysis of this data set aimed at detecting or limiting the presence of nanohertz-frequency gravitational wave signals.

Unexpected structure in images of astronomical sources often presents itself upon visual inspection of the image, but such apparent structure may either correspond to true features in the source or be due to noise in the data. This paper presents a method for testing whether inferred structure in an image with Poisson noise represents a significant departure from a baseline (null) model of the image. To infer image structure, we conduct a Bayesian analysis of a full model that uses a multiscale component to allow flexible departures from the posited null model. As a test statistic, we use a tail probability of the posterior distribution under the full model. This choice of test statistic allows us to estimate a computationally efficient upper bound on a p-value that enables us to draw strong conclusions even when there are limited computational resources that can be devoted to simulations under the null model. We demonstrate the statistical performance of our method on simulated images. Applying our method to an X-ray image of the quasar 0730+257, we find significant evidence against the null model of a single point source and uniform background, lending support to the claim of an X-ray jet.

We report first results from an ongoing monitoring campaign to measure time delays between the six images of the quasar SDSS J2222+2745, gravitationally lensed by a galaxy cluster. The time delay between A and B, the two most highly magnified images, is measured to be ${\tau }_{{\rm{AB}}}=47.7\pm 6.0$ days (95% confidence interval), consistent with previous model predictions for this lens system. The strong intrinsic variability of the quasar also allows us to derive a time delay value of ${\tau }_{{\rm{CA}}}=722\pm 24$ days between image C and A, in spite of modest overlap between their light curves in the current data set. Image C, which is predicted to lead all the other lensed quasar images, has undergone a sharp, monotonic flux increase of 60%–75% during 2014. A corresponding brightening is firmly predicted to occur in images A and B during 2016. The amplitude of this rise indicates that time delays involving all six known images in this system, including those of the demagnified central images D–F, will be obtainable from further ground-based monitoring of this system during the next few years.

We have discovered 11 ultra-faint (r ≲ 22.1) low surface brightness (LSB, central surface brightness 23 ≲ μ_r ≲ 26) dwarf galaxy candidates in one deep Virgo field of just 576 arcmin² obtained by the Large Binocular Camera at the Large Binocular Telescope. Their association with the Virgo cluster is supported by their distinct position in the central surface brightness—total magnitude plane with respect to the background galaxies of similar total magnitude. They have typical absolute magnitudes and scale sizes, if at the distance of Virgo, in the range −13 ≲ M_r ≲ −9 and 250 ≲ r_s ≲ 850 pc, respectively. Their colors are consistent with a gradually declining star formation history with a specific star formation rate of the order of 10⁻¹¹ yr⁻¹, i.e., 10 times lower than that of main sequence star-forming galaxies. They are older than the cluster formation age and appear to be regular in morphology. They represent the faintest extremes of the population of low luminosity LSB dwarfs that has recently been detected in wider surveys of the Virgo cluster. Thanks to the depth of our observations, we are able to extend the Virgo luminosity function down to M_r ∼ −9.3 (corresponding to total masses M ∼ 10⁷M_⊙), finding an average faint-end slope α ≃ −1.4. This relatively steep slope puts interesting constraints on the nature of the dark matter and, in particular, on warm dark matter (WDM) often invoked to solve the overprediction of the dwarf number density by the standard cold dark matter scenario. We derive a lower limit on the WDM particle mass >1.5 keV.

In this paper we constrain four alternative models to the late cosmic acceleration in the universe: Chevallier–Polarski–Linder (CPL), interacting dark energy (IDE), Ricci holographic dark energy (HDE), and modified polytropic Cardassian (MPC). Strong lensing (SL) images of background galaxies produced by the galaxy cluster Abell 1689 are used to test these models. To perform this analysis we modify the LENSTOOL lens modeling code. The value added by this probe is compared with other complementary probes: Type Ia supernovae (SN Ia), baryon acoustic oscillations (BAO), and cosmic microwave background (CMB). We found that the CPL constraints obtained for the SL data are consistent with those estimated using the other probes. The IDE constraints are consistent with the complementary bounds only if large errors in the SL measurements are considered. The Ricci HDE and MPC constraints are weak, but they are similar to the BAO, SN Ia, and CMB estimations. We also compute the figure of merit as a tool to quantify the goodness of fit of the data. Our results suggest that the SL method provides statistically significant constraints on the CPL parameters but is weak for those of the other models. Finally, we show that the use of the SL measurements in galaxy clusters is a promising and powerful technique to constrain cosmological models. The advantage of this method is that cosmological parameters are estimated by modeling the SL features for each underlying cosmology. These estimations could be further improved by SL constraints coming from other galaxy clusters.

Without any doubt, solar flaring loops possess a multithread internal structure that is poorly resolved, and there are no means to observe heating episodes and thermodynamic evolution of the individual threads. These limitations cause fundamental problems in numerical modeling of flaring loops, such as selection of a structure and a number of threads, and an implementation of a proper model of the energy deposition process. A set of one-dimensional (1D) hydrodynamic and two-dimensional (2D) magnetohydrodynamic models of a flaring loop are developed to compare energy redistribution and plasma dynamics in the course of a prototypical solar flare. Basic parameters of the modeled loop are set according to the progenitor M1.8 flare recorded in AR 10126 on 2002 September 20 between 09:21 UT and 09:50 UT. The nonideal 1D models include thermal conduction and radiative losses of the optically thin plasma as energy-loss mechanisms, while the nonideal 2D models take into account viscosity and thermal conduction as energy-loss mechanisms only. The 2D models have a continuous distribution of the parameters of the plasma across the loop and are powered by varying in time and space along and across the loop heating flux. We show that such 2D models are an extreme borderline case of a multithread internal structure of the flaring loop, with a filling factor equal to 1. Nevertheless, these simple models ensure the general correctness of the obtained results and can be adopted as a correct approximation of the real flaring structures.

Most coronal physicists now seem to agree that loops are composed of tangled magnetic strands and have both isothermal and multithermal cross-field temperature distributions. As yet, however, there is no information on the relative importance of each of these categories, and we do not know how common one is with respect to the other. In this paper, we investigate these temperature properties for all loop segments visible in the 171-Å image of AR 11294, which was observed by the Atmospheric Imaging Assembly (AIA) on 2011 September 15. Our analysis revealed 19 loop segments, but only 2 of these were clearly isothermal. Six additional segments were effectively isothermal, that is, the plasma emission to which AIA is sensitive could not be distinguished from isothermal emission, within measurement uncertainties. One loop had both isothermal transition region and multithermal coronal solutions. Another five loop segments require multithermal plasma to reproduce the AIA observations. The five remaining loop segments could not be separated reliably from the background in the crucial non-171-Å AIA images required for temperature analysis. We hope that the direction of coronal heating models and the efforts modelers spend on various heating scenarios will be influenced by these results.

The Earth appears non-chondritic in its abundances of refractory lithophile elements, posing a significant problem for our understanding of its formation and evolution. It has been suggested that this non-chondritic composition may be explained by collisional erosion of differentiated planetesimals of originally chondritic composition. In this work, we present N-body simulations of terrestrial planet formation that track the growth of planetary embryos from planetesimals. We simulate evolution through the runaway and oligarchic growth phases under the Grand Tack model and in the absence of giant planets. These simulations include a state-of-the-art collision model that allows multiple collision outcomes, such as accretion, erosion, and bouncing events, and enables tracking of the evolving core mass fraction of accreting planetesimals. We show that the embryos grown during this intermediate stage of planet formation exhibit a range of core mass fractions, and that with significant dynamical excitation, enough mantle can be stripped from growing embryos to account for the Earth's non-chondritic Fe/Mg ratio. We also find that there is a large diversity in the composition of remnant planetesimals, with both iron-rich and silicate-rich fragments produced via collisions.

We present studies of C/2015 D1 (SOHO), the first sunskirting comet ever seen from ground stations over the past half century. The Solar and Heliospheric Observatory (SOHO) witnessed its peculiar light curve with a huge dip followed by a flare-up around perihelion: the dip was likely caused by sublimation of olivines, directly evidenced by a coincident temporary disappearance of the tail. The flare-up likely reflects a disintegration event, which we suggest was triggered by intense thermal stress established within the nucleus interior. Photometric data reveal an increasingly dusty coma, indicative of volatile depletion. A catastrophic mass-loss rate of ∼10⁵ kg s⁻¹ around perihelion was seen. Ground-based Xingming Observatory spotted the post-perihelion debris cloud. Our morphological simulations of post-perihelion images find newly released dust grains of size a ≳ 10 μm in radius; however, a temporal increase in a_min was also witnessed, possibly owing to swift dispersions of smaller grains swept away by radiation forces without replenishment. Together with the fading profile of the light curve, a power-law dust size distribution with index γ = 3.2 ± 0.1 is derived. We detected no active remaining cometary nuclei over ∼0.1 km in radius in post-perihelion images acquired at Lowell Observatory. Applying a radial nongravitational parameter, ${{\mathcal{A}}}_{1}=\left(1.209\pm 0.118\right)\times {10}^{-6}$ AU day⁻², from an isothermal water–ice sublimation model to the SOHO astrometry significantly reduces residuals and sinusoidal trends in the orbit determination. The nucleus mass ∼10⁸–10⁹ kg and the radius ∼50–150 m (bulk density ρ_d = 0.4 g cm⁻³ assumed) before the disintegration are deduced from the photometric data; consistent results were determined from the nongravitational effects.

We perform three-dimensional radiation hydrodynamic simulations of the structure and dynamics of the radiation-dominated envelopes of massive stars at the location of the iron opacity peak. One-dimensional hydrostatic calculations predict an unstable density inversion at this location, whereas our simulations reveal a complex interplay of convective and radiative transport whose behavior depends on the ratio of the photon diffusion time to the dynamical time. The latter is set by the ratio of the optical depth per pressure scale height, ${\tau }_{0},$ to ${\tau }_{{\rm{c}}}=c/{c}_{{\rm{g}}},$ where ${c}_{{\rm{g}}}\approx 50\;\mathrm{km}\;{{\rm{s}}}^{-1}$ is the isothermal sound speed in the gas alone. When ${\tau }_{0}\gg {\tau }_{{\rm{c}}},$ convection reduces the radiation acceleration and removes the density inversion. The turbulent energy transport in the simulations agrees with mixing length theory and provides its first numerical calibration in the radiation-dominated regime. When ${\tau }_{0}\ll {\tau }_{{\rm{c}}},$ convection becomes inefficient and the turbulent energy transport is negligible. The turbulent velocities exceed c_g, driving shocks and large density fluctuations that allow photons to preferentially diffuse out through low-density regions. However, the effective radiation acceleration is still larger than the gravitational acceleration so that the time average density profile contains a modest density inversion. In addition, the simulated envelope undergoes large-scale oscillations with periods of a few hours. The turbulent velocity field may affect the broadening of spectral lines and therefore stellar rotation measurements in massive stars, while the time variable outer atmosphere could lead to variations in their mass loss and stellar radius.

We measure the binarity of detached M dwarfs in the Kepler field with orbital periods in the range of 1–90 days. Kepler's photometric precision and nearly continuous monitoring of stellar targets over time baselines ranging from 3 months to 4 years make its detection efficiency for eclipsing binaries nearly complete over this period range and for all radius ratios. Our investigation employs a statistical framework akin to that used for inferring planetary occurrence rates from planetary transits. The obvious simplification is that eclipsing binaries have a vastly improved detection efficiency that is limited chiefly by their geometric probabilities to eclipse. For the M-dwarf sample observed by the Kepler Mission, the fractional incidence of eclipsing binaries implies that there are ${0.11}_{-0.04}^{+0.02}$ close stellar companions per apparently single M dwarf. Our measured binarity is higher than previous inferences of the occurrence rate of close binaries via radial velocity techniques, at roughly the 2σ level. This study represents the first use of eclipsing binary detections from a high quality transiting planet mission to infer binary statistics. Application of this statistical framework to the eclipsing binaries discovered by future transit surveys will establish better constraints on short-period M+M binary rate, as well as binarity measurements for stars of other spectral types.

The formation of planetesimals requires that primordial dust grains grow from micron- to kilometer-sized bodies. Dust traps caused by gas pressure maxima have been proposed as regions where grains can concentrate and grow fast enough to form planetesimals, before radially migrating onto the star. We report new VLA Ka and Ku observations of the protoplanetary disk around the Herbig Ae/Be star MWC 758. The Ka image shows a compact emission region in the outer disk, indicating a strong concentration of big dust grains. Tracing smaller grains, archival ALMA data in band 7 continuum shows extended disk emission with an intensity maximum to the northwest of the central star, which matches the VLA clump position. The compactness of the Ka emission is expected in the context of dust trapping, as big grains are trapped more easily than smaller grains in gas pressure maxima. We develop a nonaxisymmetric parametric model inspired by a steady-state vortex solution with parameters adequately selected to reproduce the observations, including the spectral energy distribution. Finally, we compare the radio continuum with SPHERE scattered light data. The ALMA continuum spatially coincides with a spiral-like feature seen in scattered light, while the VLA clump is offset from the scattered light maximum. Moreover, the ALMA map shows a decrement that matches a region devoid of scattered polarized emission. Continuum observations at a different wavelength are necessary to conclude whether the VLA-ALMA difference is an opacity or a real dust segregation.

Diffuse radio emission in galaxy clusters is known to be related to cluster mass and cluster dynamical state. We collect the observed fluxes of radio halos, relics, and mini-halos for a sample of galaxy clusters from the literature, and calculate their radio powers. We then obtain the values of cluster mass or mass proxies from previous observations, and also obtain the various dynamical parameters of these galaxy clusters from optical and X-ray data. The radio powers of relics, halos, and mini-halos are correlated with the cluster masses or mass proxies, as found by previous authors, while the correlations concerning giant radio halos are in general the strongest. We found that the inclusion of dynamical parameters as the third dimension can significantly reduce the data scatter for the scaling relations, especially for radio halos. We therefore conclude that the substructures in X-ray images of galaxy clusters and the irregular distributions of optical brightness of member galaxies can be used to quantitatively characterize the shock waves and turbulence in the intracluster medium responsible for re-accelerating particles to generate the observed diffuse radio emission. The power of radio halos and relics is correlated with cluster mass proxies and dynamical parameters in the form of a fundamental plane.

We present the results from a stellar population modeling analysis of a sample of 162 z = 4.5 and 14 z = 5.7 Lyα emitting galaxies (LAEs) in the Boötes field, using deep Spitzer/IRAC data at 3.6 and 4.5 μm from the Spitzer Lyα Survey, along with Hubble Space Telescope NICMOS and WFC3 imaging at 1.1 and 1.6 μm for a subset of the LAEs. This represents one of the largest samples of high-redshift LAEs imaged with Spitzer IRAC. We find that 30/162 (19%) of the z = 4.5 LAEs and 9/14 (64%) of the z = 5.7 LAEs are detected at ≥3σ in at least one IRAC band. Individual z = 4.5 IRAC-detected LAEs have a large range of stellar mass, from 5 × 10⁸–10¹¹${M}_{\odot }$. One-third of the IRAC-detected LAEs have older stellar population ages of 100 Myr–1 Gyr, while the remainder have ages <100 Myr. A stacking analysis of IRAC-undetected LAEs shows this population to be primarily low mass (8–20 × 10⁸${M}_{\odot }$) and young (64–570 Myr). We find a correlation between stellar mass and the dust-corrected ultraviolet-based star formation rate (SFR) similar to that at lower redshifts, in that higher mass galaxies exhibit higher SFRs. However, the z = 4.5 LAE correlation is elevated 4–5 times in SFR compared to continuum-selected galaxies at similar redshifts. The exception is the most massive LAEs which have SFRs similar to galaxies at lower redshifts suggesting that they may represent a different population of galaxies than the traditional lower-mass LAEs, perhaps with a different mechanism promoting Lyα photon escape.