Table of contents

It is shown here how new experimental data, for the electrical properties of solid CO, can be used to fill important gaps in our understanding of the evolution of prestellar cores. Dust grains with a mantle of CO lead to a reduction in the degree of ionization in these cores by a factor of between 5 and 6. The lifetimes for expulsion of magnetic fields from cores, a process generally necessary for gravitational collapse, are reduced from current estimates of several megayears, by a similar factor. This removes a major inconsistency, since lifetimes now tally with typical ages of prestellar cores of a few hundred thousand to 10⁶ yr, derived from observations. With the reduced timescales, cores also escape disruption by Galactic supernova remnants. Our results provide a natural mechanism for the generation of so-called magnetically supercritical cores, in which the magnetic field alone cannot prevent gravitational collapse. In addition, we find a minimum value for the density of prestellar cores of ≥(1.1 ± 0.1) × 10⁴ H₂ cm⁻³, in agreement with observations.

Dense plasma fragments were observed to fall back on the solar surface by the Solar Dynamics Observatory after an eruption on 2011 June 7, producing strong extreme-ultraviolet brightenings. Previous studies investigated impacts in regions of weak magnetic field. Here we model the $\sim 300$ km s⁻¹ impact of fragments channelled by the magnetic field close to active regions. In the observations, the magnetic channel brightens before the fragment impact. We use a 3D-MHD model of spherical blobs downfalling in a magnetized atmosphere. The blob parameters are constrained from the observation. We run numerical simulations with different ambient densitie and magnetic field intensities. We compare the model emission in the 171 Å channel of the Atmospheric Imaging Assembly with the observed one. We find that a model of downfall channelled in an ∼1 MK coronal loop confined by a magnetic field of ∼10–20 G, best explains qualitatively and quantitatively the observed evolution. The blobs are highly deformed and further fragmented when the ram pressure becomes comparable to the local magnetic pressure, and they are deviated to be channelled by the field because of the differential stress applied by the perturbed magnetic field. Ahead of them, in the relatively dense coronal medium, shock fronts propagate, heat, and brighten the channel between the cold falling plasma and the solar surface. This study shows a new mechanism that brightens downflows channelled by the magnetic field, such as in accreting young stars, and also works as a probe of the ambient atmosphere, providing information about the local plasma density and magnetic field.

We present a multi-wavelength polarimetric and spectral study of the M87 jet obtained at sub-arcsecond resolution between 2002 and 2008. The observations include multi-band archival VLA polarimetry data sets along with Hubble Space Telescope (HST) imaging polarimetry. These observations have better angular resolution than previous work by factors of 2–3 and in addition, allow us to explore the time domain. These observations envelop the huge flare in HST-1 located 086 from the nucleus. The increased resolution enables us to view more structure in each knot, showing several resolved sub-components. We also see apparent helical structure in the polarization vectors in several knots, with polarization vectors turning either clockwise or counterclockwise near the flux maxima in various places as well as showing filamentary undulations. Some of these characteristics are correlated with flux and polarization maxima while others are not. We also examine the total flux and fractional polarization and look for changes in both radio and optical since the observations of Perlman et al. (1999) and test them against various models based on shocks and instabilities in the jet. Our results are broadly consistent with previous spine-sheath models and recollimation shock models; however, they require additional combinations of features to explain the observed complexity, e.g., shearing of magnetic field lines near the jet surface and compression of the toroidal component near shocks. In particular, in many regions we find apparently helical features both in total flux and polarization. We discuss the physical interpretation of these features.

We present new SOFIA-FORCAST observations obtained in 2016 February of the archetypal outbursting low-mass young stellar object FU Orionis, and we compare the continuum, solid-state, and gas properties with mid-infrared data obtained at the same wavelengths in 2004 with Spitzer-IRS. In this study, we conduct the first mid-infrared spectroscopic comparison of an FUor over a long time period. Over a 12-year period, UBVR monitoring indicates that FU Orionis has continued its steady decrease in overall brightness by ∼14%. We find that this decrease in luminosity occurs only at wavelengths ≲20 μm. In particular, the continuum shortward of the silicate emission complex at 10 μm exhibits a ∼12% (∼3σ) drop in flux density but no apparent change in slope; both the Spitzer and SOFIA spectra are consistent with a 7200 K blackbody. Additionally, the detection of water absorption is consistent with the Spitzer spectrum. The silicate emission feature at 10 μm continues to be consistent with unprocessed grains, unchanged over 12 years. We conclude that either the accretion rate in FU Orionis has decreased by ∼12–14% over this time baseline or the inner disk has cooled, but the accretion disk remains in a superheated state outside the innermost region.

CO₂ ice is an important reservoir of carbon and oxygen in star- and planet-forming regions. Together with water and CO, CO₂ sets the physical and chemical characteristics of interstellar icy grain mantles, including desorption and diffusion energies for other ice constituents. A detailed understanding of CO₂ ice spectroscopy is a prerequisite to characterize CO₂ interactions with other volatiles both in interstellar ices and in laboratory experiments of interstellar ice analogs. We report laboratory spectra of the CO₂ longitudinal optical (LO) phonon mode in pure CO₂ ice and in CO₂ ice mixtures with H₂O, CO, and O₂ components. We show that the LO phonon mode position is sensitive to the mixing ratio of various ice components of astronomical interest. In the era of the James Webb Space Telescope, this characteristic could be used to constrain interstellar ice compositions and morphologies. More immediately, LO phonon mode spectroscopy provides a sensitive probe of ice mixing in the laboratory and should thus enable diffusion measurements with higher precision than has been previously possible.

The effects of a non-gradient flux term originating from the motion of convective elements with entropy perturbations of either sign are investigated and incorporated into a modified version of stellar mixing length theory (MLT). Such a term, first studied by Deardorff in the meteorological context, might represent the effects of cold intense downdrafts caused by the rapid cooling in the granulation layer at the top of the convection zone of late-type stars. These intense downdrafts were first seen in the strongly stratified simulations of Stein & Nordlund in the late 1980s. These downdrafts transport heat nonlocally, a phenomenon referred to as entropy rain. Moreover, the Deardorff term can cause upward enthalpy transport even in a weakly Schwarzschild-stably stratified layer. In that case, no giant cell convection would be excited. This is interesting in view of recent observations, which could be explained if the dominant flow structures were of small scale even at larger depths. To study this possibility, three distinct flow structures are examined: one in which convective structures have similar size and mutual separation at all depths, one in which the separation increases with depth, but their size is still unchanged, and one in which both size and separation increase with depth, which is the standard flow structure. It is concluded that the third possibility with fewer and thicker downdrafts in deeper layers remains the most plausible, but it may be unable to explain the suspected absence of large-scale flows with speeds and scales expected from MLT.

Knowledge of galaxy evolution rests on cross-sectional observations of different objects at different times. Understanding of galaxy evolution rests on longitudinal interpretations of how these data relate to individual objects moving through time. The connection between the two is often assumed to be clear, but we use a simple "physics-free" model to show that it is not and that exploring its nuances can yield new insights. Comprising nothing more than 2094 loosely constrained lognormal star formation histories (SFHs), the model faithfully reproduces the following data it was not designed to match: stellar mass functions at $z\leqslant 8;$ the slope of the star formation rate/stellar mass relation (the SFR "Main Sequence") at $z\leqslant 6;$ the mean $\mathrm{sSFR}(\equiv \mathrm{SFR}/{M}_{* })$ of low-mass galaxies at $z\leqslant 7;$ "fast-" and "slow-track" quenching; downsizing; and a correlation between formation timescale and $\mathrm{sSFR}({M}_{* },t)$ similar to results from simulations that provides a natural connection to bulge growth. We take these findings—which suggest that quenching is the natural downturn of all SFHs affecting galaxies at rates/times correlated with their densities—to mean that: (1) models in which galaxies are diversified on Hubble timescales by something like initial conditions rival the dominant grow-and-quench framework as good descriptions of the data; or (2) absent spatial information, many metrics of galaxy evolution are too undiscriminating—if not inherently misleading—to confirm a unique explanation. We outline future tests of our model but stress that, even if ultimately incorrect, it illustrates how exploring different paradigms can aid learning and, we hope, more detailed modeling efforts.

The typical optical–UV continuum slopes observed in many type-1 active galactic nuclei (AGNs) are redder than expected from thin accretion disk (AD) models. A possible resolution to this conundrum is that many AGNs are reddened by dust along the line of sight. To explore this possibility, we stack 5000 SDSS AGNs with luminosity $L\approx {10}^{45}\,\mathrm{erg}\,{{\rm{s}}}^{-1}$ and redshift $z\sim 0.4$ in bins of optical continuum slope ${\alpha }_{\mathrm{opt}}$ and width of the broad Hβ emission line. We measure the equivalent width (EW) of the NaID absorption feature in each stacked spectrum. We find a linear relation between ${\alpha }_{\mathrm{opt}}$ and EW(NaID), such that EW(NaID) increases as ${\alpha }_{\mathrm{opt}}$ becomes redder. In the bin with the smallest Hβ width, objects with the bluest slopes, which are similar to AD predictions, are found to have $\mathrm{EW}(\mathrm{NaID})=0$, supporting the line of sight dust hypothesis. This conclusion is also supported by the dependence of the Hα/Hβ line ratio on ${\alpha }_{\mathrm{opt}}$. The implied relationship between continuum slope and dust reddening is given by ${E}_{{\rm{B}}-{\rm{V}}}\approx 0.2\cdot (-0.1-{\alpha }_{\mathrm{opt}})$, and the implied reddening of a typical type-1 AGN with ${\alpha }_{\mathrm{opt}}=-0.5$ is ${E}_{{\rm{B}}-{\rm{V}}}\approx 0.08\,\mathrm{mag}$. Photoionization calculations show that the line of sight dusty gas responsible for reddening is too ionized to produce the observed sodium features. Therefore, we argue that the sodium absorption arises in regions of the host ISM that are shielded from the AGN radiation along lines of sight to the stars, and the correlation with ${\alpha }_{\mathrm{opt}}$ arises since ISM columns along shielded and non-shielded sightlines are correlated. This scenario is supported by the similarity of the relation between ${E}_{{\rm{B}}-{\rm{V}}}$ and the Na i column implied by our results with the relation in the Milky Way found by previous studies.

At present, the Babcock–Leighton flux transport solar dynamo models appear to be the most promising models for explaining diverse observational aspects of the sunspot cycle. The success of these flux transport dynamo models is largely dependent upon a single-cell meridional circulation with a deep equatorward component at the base of the Sun's convection zone. However, recent observations suggest that the meridional flow may in fact be very shallow (confined to the top 10% of the Sun) and more complex than previously thought. Taken together, these observations raise serious concerns on the validity of the flux transport paradigm. By accounting for the turbulent pumping of magnetic flux, as evidenced in magnetohydrodynamic simulations of solar convection, we demonstrate that flux transport dynamo models can generate solar-like magnetic cycles even if the meridional flow is shallow. Solar-like periodic reversals are recovered even when meridional circulation is altogether absent. However, in this case, the solar surface magnetic field dynamics does not extend all the way to the polar regions. Very importantly, our results demonstrate that the Parker–Yoshimura sign rule for dynamo wave propagation can be circumvented in Babcock–Leighton dynamo models by the latitudinal component of turbulent pumping, which can generate equatorward propagating sunspot belts in the absence of a deep, equatorward meridional flow. We also show that variations in turbulent pumping coefficients can modulate the solar cycle amplitude and periodicity. Our results suggest the viability of an alternate magnetic flux transport paradigm—mediated via turbulent pumping—for sustaining solar-stellar dynamo action.

Stars ejected from the Galactic Center can be used to place important constraints on the Milky Way potential. Since existing hypervelocity stars are too distant to accurately determine orbits, we have conducted a search for nearby candidates using full three-dimensional velocities. Since the efficacy of such studies is often hampered by deficiencies in proper motion catalogs, we have chosen to utilize the reliable, high-precision Sloan Digital Sky Survey (SDSS) Stripe 82 proper motion catalog. Although we do not find any candidates which have velocities in excess of the escape speed, we identify 226 stars on orbits that are consistent with Galactic Center ejection. This number is significantly larger than what we would expect for halo stars on radial orbits and cannot be explained by disk or bulge contamination. If we restrict ourselves to metal-rich stars, we find 29 candidates with [Fe/H] > −0.8 dex and 10 with [Fe/H] > −0.6 dex. Their metallicities are more consistent with what we expect for bulge ejecta, and so we believe these candidates are especially deserving of further study. We have supplemented this sample using our own radial velocities, developing an algorithm to use proper motions for optimizing candidate selection. This technique provides considerable improvement on the blind spectroscopic sample of SDSS, being able to identify candidates with an efficiency around 20 times better than a blind search.

The baryonic Tully–Fisher relation (BTFR) is both a valuable observational tool and a critical test of galaxy formation theory. We explore the systematic uncertainty in the slope and the scatter of the observed line-width BTFR utilizing homogeneously measured, unresolved H i observations for 930 isolated galaxies. We measure a fiducial relation of $\,{\mathrm{log}}_{10}\,{M}_{\mathrm{baryon}}=3.24\,{\mathrm{log}}_{10}\,{V}_{\mathrm{rot}}\,+\,3.21$ with observed scatter of 0.25 dex over a baryonic mass range of 10^7.4 to 10^11.3${M}_{\odot }$ where ${V}_{\mathrm{rot}}$ is measured from 20% H i line widths. We then conservatively vary the definitions of ${M}_{\mathrm{baryon}}$ and ${V}_{\mathrm{rot}}$, the sample definition and the linear fitting algorithm. We obtain slopes ranging from 2.64 to 3.53 and scatter measurements ranging from 0.14 to 0.41 dex, indicating a significant systematic uncertainty of 0.25 in the BTFR slope derived from unresolved H i line widths. We next compare our fiducial slope to literature measurements, where reported slopes range from 3.0 to 4.3 and scatter is either unmeasured, immeasurable, or as large as 0.4 dex. Measurements derived from unresolved H i line widths tend to produce slopes of 3.3, while measurements derived strictly from resolved asymptotic rotation velocities tend to produce slopes of 3.9. The single largest factor affecting the BTFR slope is the definition of rotation velocity. The sample definition, the mass range and the linear fitting algorithm also significantly affect the measured BTFR. We find that galaxies in our sample with ${V}_{\mathrm{rot}}\lt 100$ km s⁻¹ are consistent with the line-width BTFR of more massive galaxies, but these galaxies drive most of the observed scatter. It is critical when comparing predictions to an observed BTFR that the rotation velocity definition, the sample selection and the fitting algorithm are similarly defined. We recommend direct statistical comparisons between data sets with commensurable properties as opposed to simply comparing BTFR power-law fits.

We present results from a comprehensive submillimeter spectral survey toward the source Orion South, based on data obtained with the Heterodyne Instrument for the Far-Infrared instrument on board the Herschel Space Observatory, covering the frequency range of 480 to 1900 GHz. We detect 685 spectral lines with signal-to-noise ratios (S/Ns) > 3σ, originating from 52 different molecular and atomic species. We model each of the detected species assuming conditions of Local Thermodynamic Equilibrium. This analysis provides an estimate of the physical conditions of Orion South (column density, temperature, source size, and V_LSR). We find evidence for three different cloud components: a cool (T_ex ∼ 20–40 K), spatially extended (>60''), and quiescent (ΔV_FWHM ∼ 4 km s⁻¹) component; a warmer (T_ex ∼ 80–100 K), less spatially extended (∼30''), and dynamic (ΔV_FWHM ∼ 8 km s⁻¹) component, which is likely affected by embedded outflows; and a kinematically distinct region (T_ex > 100 K; V_LSR ∼ 8 km s⁻¹), dominated by emission from species that trace ultraviolet irradiation, likely at the surface of the cloud. We find little evidence for the existence of a chemically distinct "hot-core" component, likely due to the small filling factor of the hot core or hot cores within the Herschel beam. We find that the chemical composition of the gas in the cooler, quiescent component of Orion South more closely resembles that of the quiescent ridge in Orion-KL. The gas in the warmer, dynamic component, however, more closely resembles that of the Compact Ridge and Plateau regions of Orion-KL, suggesting that higher temperatures and shocks also have an influence on the overall chemistry of Orion South.

Motivated by recent results in stellar evolution that predict the existence of hybrid white dwarf (WD) stars with a C–O core inside an O–Ne shell, we simulate thermonuclear (Type Ia) supernovae from these hybrid progenitors. We use the FLASH code to perform multidimensional simulations in the deflagration-to-detonation transition (DDT) explosion paradigm. Our hybrid progenitor models were produced with the MESA stellar evolution code and include the effects of the Urca process, and we map the progenitor model to the FLASH grid. We performed a suite of DDT simulations over a range of ignition conditions consistent with the progenitor's thermal and convective structure assuming multiple ignition points. To compare the results from these hybrid WD stars to previous results from C–O WDs, we construct a set of C–O WD models with similar properties and similarly simulate a suite of explosions. We find that despite significant variability within each suite, trends distinguishing the explosions are apparent in their ${}^{56}\mathrm{Ni}$ yields and the kinetic properties of the ejecta. We compare our results with other recent work that studies explosions from these hybrid progenitors.

Motivated by the significant interaction of convection, rotation, and magnetic field in many astrophysical objects, we investigate the interplay between large-scale flows driven by rotating convection and an imposed magnetic field. We utilize a simple model in two dimensions comprised of a plane layer that is rotating about an axis inclined to gravity. It is known that this setup can result in strong mean flows; we numerically examine the effect of an imposed horizontal magnetic field on such flows. We show that increasing the field strength in general suppresses the time-dependent mean flows, but in some cases it organizes them, leading to stronger time-averaged flows. Furthermore, we discuss the effect of the field on the correlations responsible for driving the flows and the competition between Reynolds and Maxwell stresses. A change in behavior is observed when the (fluid and magnetic) Prandtl numbers are decreased. In the smaller Prandtl number regime, it is shown that significant mean flows can persist even when the quenching of the overall flow velocity by the field is relatively strong.

The origin of the broad line region (BLR) in active galaxies remains unknown. It seems to be related to the underlying accretion disk, but an efficient mechanism is required to raise the material from the disk surface without giving signatures of the outflow that are too strong in the case of the low ionization lines. We discuss in detail two proposed mechanisms: (1) radiation pressure acting on dust in the disk atmosphere creating a failed wind and (2) the gravitational instability of the underlying disk. We compare the predicted location of the inner radius of the BLR in those two scenarios with the observed position obtained from the reverberation studies of several active galaxies. The failed dusty outflow model well represents the observational data while the predictions of the self-gravitational instability are not consistent with observations. The issue that remains is why do we not see any imprints of the underlying disk instability in the BLR properties.

We use the TESIS EUV telescope to study the current sheet signatures observed during flux rope eruption. The special feature of the TESIS telescope was its ability to image the solar corona up to a distance of 2 ${R}_{\odot }$ from the Sun's center in the Fe 171 Å line. The Fe 171 Å line emission illuminates the magnetic field lines, and the TESIS images reveal the coronal magnetic structure at high altitudes. The analyzed coronal mass ejection (CME) had a core with a spiral—flux rope—structure. The spiral shape indicates that the flux rope radius varied along its length. The flux rope had a complex temperature structure: cold legs (70,000 K, observed in He 304 Å line) and a hotter core (0.7 MK, observed in Fe 171 Å line). Such a structure contradicts the common assumption that the CME core is a cold prominence. When the CME impulsively accelerated, a dark double Y-structure appeared below the flux rope. The Y-structure timing, location, and morphology agree with the previously performed MHD simulations of the current sheet. We interpreted the Y-structure as a hot envelope of the current sheet and hot reconnection outflows. The Y-structure had a thickness of 6.0 Mm. Its length increased over time from 79 Mm to more than 411 Mm.

Blazars radiate from radio through gamma-ray frequencies and thereby make ideal targets for multifrequency studies. Such studies allow the properties of the emitting jet to be constrained. 3C 279 is among the most notable blazars and therefore subject to extensive multifrequency campaigns. We report the results of a campaign ranging from near-IR to gamma-ray energies that targeted an outburst of 3C 279 in 2015 June. The campaign pivots around the detection in only 50 ks by INTEGRAL, whose IBIS/ISGRI data pin down the high-energy component of the spectral energy distribution (SED) between Swift-XRT data and Fermi-LAT data. The overall SED from near-IR to gamma rays can be well represented by either a leptonic or a lepto-hadronic radiation transfer model. Even though the data are equally well represented by the two models, their inferred parameters challenge the physical conditions in the jet. In fact, the leptonic model requires parameters with a magnetic field far below equipartition with the relativistic particle energy density. In contrast, equipartition may be achieved with the lepto-hadronic model, although this implies an extreme total jet power close to the Eddington luminosity.

Magnetic fields (B-fields) play a key role in the formation and evolution of protoplanetary disks, but their properties are poorly understood due to the lack of observational constraints. Using CanariCam at the 10.4 m Gran Telescopio Canarias, we have mapped out the mid-infrared polarization of the protoplanetary disk around the Herbig Ae star AB Aur. We detect ∼0.44% polarization at 10.3 μm from AB Aur's inner disk (r < 80 au), rising to ∼1.4% at larger radii. Our simulations imply that the mid-infrared polarization of the inner disk arises from dichroic emission of elongated particles aligned in a disk B-field. The field is well ordered on a spatial scale, commensurate with our resolution (∼50 au), and we infer a poloidal shape tilted from the rotational axis of the disk. The disk of AB Aur is optically thick at 10.3 μm, so polarimetry at this wavelength is probing the B-field near the disk surface. Our observations therefore confirm that this layer, favored by some theoretical studies for developing magneto-rotational instability and its resultant viscosity, is indeed very likely to be magnetized. At radii beyond ∼80 au, the mid-infrared polarization results primarily from scattering by dust grains with sizes up to ∼1 μm, a size indicating both grain growth and, probably, turbulent lofting of the particles from the disk mid-plane.

Early quiescent galaxies at $z\sim 2$ are known to be remarkably compact compared to their nearby counterparts. Possible progenitors of these systems include galaxies that are structurally similar, but are still rapidly forming stars. Here, we present Karl G. Jansky Very Large Array (VLA) observations of the CO(1–0) line toward three such compact, star-forming galaxies (SFGs) at $z\sim 2.3$, significantly detecting one. The VLA observations indicate baryonic gas fractions $\gtrsim 5$ times lower and gas depletion timescales $\gtrsim 10$ times shorter than normal, extended massive SFGs at these redshifts. At their current star formation rates, all three objects will deplete their gas reservoirs within 100 Myr. These objects are among the most gas-poor objects observed at $z\gt 2$, and are outliers from standard gas scaling relations, a result that remains true regardless of assumptions about the CO–H₂ conversion factor. Our observations are consistent with the idea that compact, SFGs are in a rapid state of transition to quiescence in tandem with the buildup of the $z\sim 2$ quenched population. In the detected compact galaxy, we see no evidence of rotation or that the CO-emitting gas is spatially extended relative to the stellar light. This casts doubt on recent suggestions that the gas in these compact galaxies is rotating and significantly extended compared to the stars. Instead, we suggest that, at least for this object, the gas is centrally concentrated, and only traces a small fraction of the total galaxy dynamical mass.

The merger of two white dwarfs is expected to result in a central fast-rotating core surrounded by a debris disk, in which magnetorotational instabilities give rise to a hot magnetized corona and a magnetized outflow. The dissipation of magnetic energy via reconnection could lead to the acceleration of cosmic-rays (CRs) in the expanding material, which would result in high energy neutrinos. We discuss the possibility of using these neutrino signals as probes of the outflow dynamics, magnetic energy dissipation rate, and CR acceleration efficiency. Importantly, the accompanying high-energy gamma-rays are absorbed within these sources because of the large optical depth, so these neutrino sources can be regarded as hidden cosmic-ray accelerators that are consistent with the non-detection of gamma-rays with Fermi-LAT. While the CR generation rate is highly uncertain, if it reaches $\sim {10}^{45}\,\mathrm{erg}\,{\mathrm{Mpc}}^{-3}\,{\mathrm{yr}}^{-1}$, the diffuse neutrino flux could contribute a substantial fraction of the IceCube observations. We also evaluate the prospect of observing individual merger events, which provides a means for testing such sources in the future.

We report the discovery of a new ultra-faint dwarf satellite companion of the Milky Way (MW) based on the early survey data from the Hyper Suprime-Cam Subaru Strategic Program. This new satellite, Virgo I, which is located in the constellation of Virgo, has been identified as a statistically significant (5.5σ) spatial overdensity of star-like objects with a well-defined main sequence and red giant branch in the color–magnitude diagram. The significance of this overdensity increases to 10.8σ when the relevant isochrone filter is adopted for the search. Based on the distribution of the stars around the likely main-sequence turnoff at r ∼ 24 mag, the distance to Virgo I is estimated as 87 kpc, and its most likely absolute magnitude calculated from a Monte Carlo analysis is M_V = −0.8 ± 0.9 mag. This stellar system has an extended spatial distribution with a half-light radius of ${38}_{-11}^{+12}$ pc, which clearly distinguishes it from a globular cluster with comparable luminosity. Thus, Virgo I is one of the faintest dwarf satellites known and is located beyond the reach of the Sloan Digital Sky Survey. This demonstrates the power of this survey program to identify very faint dwarf satellites. This discovery of Virgo I is based only on about 100 square degrees of data, thus a large number of faint dwarf satellites are likely to exist in the outer halo of the MW.

When a circumbinary disk surrounds a binary whose secondary's mass is at least $\sim {10}^{-2}\times$ the primary's mass, a nearly empty cavity with radius a few times the binary separation is carved out of the disk. Narrow streams of material pass from the inner edge of the circumbinary disk into the domain of the binary itself, where they eventually join onto the small disks orbiting the members of the binary. Using data from three-dimensional magnetohydrodynamics simulations of this process, we determine the luminosity of these streams; it is mostly due to weak laminar shocks, and is in general only a few percent of the luminosity of adjacent regions of either the circumbinary disk or the "mini-disks." This luminosity therefore hardly affects the deficit in the thermal continuum predicted on the basis of a perfectly dark gap region.

We present a detailed study of the formation of an inverse S-shaped filament prior to its eruption in active region NOAA 11884 from 2013 October 31 to November 2. In the initial stage, clockwise rotation of a small positive sunspot around the main negative trailing sunspot formed a curved filament. Then the small sunspot cancelled with the negative magnetic flux to create a longer active-region filament with an inverse S-shape. At the cancellation site a brightening was observed in UV and EUV images and bright material was transferred to the filament. Later the filament erupted after cancellation of two opposite polarities below the upper part of the filament. Nonlinear force-free field extrapolation of vector photospheric fields suggests that the filament may have a twisted structure, but this cannot be confirmed from the current observations.

We have investigated the photo-stability of pristine and super-hydrogenated pyrene cations (${{\rm{C}}}_{16}{{\rm{H}}}_{10+m}^{+},m=0,6,{\rm{}}\,{\rm{or}}\,{\rm{}}16$) by means of gas-phase action spectroscopy. Optical absorption spectra and photo-induced dissociation mass spectra are presented. By measuring the yield of mass-selected photo-fragment ions as a function of laser pulse intensity, the number of photons (and hence the energy) needed for fragmentation of the carbon backbone was determined. Backbone fragmentation of pristine pyrene ions (${{\rm{C}}}_{16}{{\rm{H}}}_{10}^{+}$) requires absorption of three photons of energy just below 3 eV, whereas super-hydrogenated hexahydropyrene (C₁₆H${}_{16}^{+}$) must absorb two such photons and fully hydrogenated hexadecahydropyrene (C₁₆H${}_{26}^{+}$) only a single photon. These results are consistent with previously reported dissociation energies for these ions. Our experiments clearly demonstrate that the increased heat capacity from the additional hydrogen atoms does not compensate for the weakening of the carbon backbone when pyrene is hydrogenated. In photodissociation regions, super-hydrogenated Polycyclic Aromatic Hydrocarbons (PAHs) have been proposed to serve as catalysts for H₂ formation. Our results indicate that carbon backbone fragmentation may be a serious competitor to H₂ formation at least for small hydrogenated PAHs like pyrene.

Irradiation of solid nitrogen at 4 K with far-ultraviolet light from a synchrotron caused excitation to the upper state of the Vegard–Kaplan (VK) system; the emission in that system was simultaneously recorded in wavelength region 200–440 nm. The lifetimes of emission lines for VK (0, 1) to (0, 12) transitions were measured in the range of 2.12 ∼ 2.65 s. The threshold wavelength to observe the VK emission was 175.0 ± 3.5 nm, corresponding to energy 7.08 ± 0.14 eV. This investigation of the generation of icy VK nitrogen enhances our understanding of its photochemistry in space.

We measure the Sunyaev–Zel'dovich (SZ) signal toward a set of 47 clusters with a median mass of 9.5 × 10¹⁴M_☉ and a median redshift of 0.40 using data from Planck and the ground-based Bolocam receiver. When Planck XMM-like masses are used to set the scale radius ${\theta }_{{\rm{s}}}$, we find consistency between the integrated SZ signal, ${Y}_{5{\rm{R}}500}$, derived from Bolocam and Planck based on generalized Navarro, Frenk, and White model fits using A10 shape parameters, with an average ratio of 1.069 ± 0.030 (allowing for the ≃5% Bolocam flux calibration uncertainty). We also perform a joint fit to the Bolocam and Planck data using a modified A10 model with the outer logarithmic slope β allowed to vary, finding β = 6.13 ± 0.16 ± 0.76 (measurement error followed by intrinsic scatter). In addition, we find that the value of β scales with mass and redshift according to $\beta \propto {M}^{0.077\pm 0.026}\times {(1+z)}^{-0.06\pm 0.09}$. This mass scaling is in good agreement with recent simulations. We do not observe the strong trend of β with redshift seen in simulations, though we conclude that this is most likely due to our sample selection. Finally, we use Bolocam measurements of Y₅₀₀ to test the accuracy of the Planck completeness estimate. We find consistency, with the actual number of Planck detections falling approximately 1σ below the expectation from Bolocam. We translate this small difference into a constraint on the effective mass bias for the Planck cluster cosmology results, with $(1-b)=0.93\pm 0.06$.

This study entails the third part of a global flare energetics project, in which Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) data of 191 M and X-class flare events from the first 3.5 years of the Solar Dynamics Observatory mission are analyzed. We fit a thermal and a nonthermal component to RHESSI spectra, yielding the temperature of the differential emission measure (DEM) tail, the nonthermal power-law slope and flux, and the thermal/nonthermal cross-over energy e_co. From these parameters, we calculate the total nonthermal energy E_nt in electrons with two different methods: (1) using the observed cross-over energy e_co as low-energy cutoff, and (2) using the low-energy cutoff e_wt predicted by the warm thick-target bremsstrahlung model of Kontar et al. Based on a mean temperature of T_e = 8.6 MK in active regions, we find low-energy cutoff energies of ${e}_{\mathrm{wt}}=6.2\pm 1.6\,\mathrm{keV}$ for the warm-target model, which is significantly lower than the cross-over energies ${e}_{\mathrm{co}}=21\pm 6\,\mathrm{keV}$. Comparing with the statistics of magnetically dissipated energies E_mag and thermal energies E_th from the two previous studies, we find the following mean (logarithmic) energy ratios with the warm-target model: ${E}_{\mathrm{nt}}=0.41\ {E}_{\mathrm{mag}}$, ${E}_{\mathrm{th}}=0.08\ {E}_{\mathrm{mag}}$, and ${E}_{\mathrm{th}}=0.15\ {E}_{\mathrm{nt}}$. The total dissipated magnetic energy exceeds the thermal energy in 95% and the nonthermal energy in 71% of the flare events, which confirms that magnetic reconnection processes are sufficient to explain flare energies. The nonthermal energy exceeds the thermal energy in 85% of the events, which largely confirms the warm thick-target model.

We numerically investigate global properties of rotating neutron stars (NSs) using the allowed band of QCD equations of state derived by Kurkela et al. This band is constrained by chiral effective theory at low densities and perturbative QCD at high densities, and is thus, in essence, a controlled constraint from first-principles physics. Previously, this band of equations of state was used to investigate non-rotating NSs only; in this work, we extend these results to any rotation frequency below the mass-shedding limit. We investigate mass–radius curves, allowed mass–frequency regions, radius–frequency curves for a typical $1.4{M}_{\odot }$ star, and the values of the moment of inertia of the double pulsar PSR J0737-3039A, a pulsar for which the moment of inertia may be constrained observationally in a few years. We present limits on observational data coming from these constraints, and identify values of observationally relevant parameters that would further constrain the allowed region for the QCD equation of state. We also discuss how much this region would be constrained by a measurement of the moment of inertia of the double pulsar PSR J0737-3039A.

The average star formation rate (SFR) in galaxies has been declining since the redshift of 2. A fraction of galaxies quench and become quiescent. We constrain two key properties of the quenching process: the quenching timescale and the quenching rate among galaxies. We achieve this by analyzing the galaxy number density profile in NUV−u color space and the distribution in NUV−u versus u − i color–color diagram with a simple toy-model framework. We focus on galaxies in three mass bins between 10¹⁰ and 10^10.6M_⊙. In the NUV−u versus u − i color–color diagram, the red u − i galaxies exhibit a different slope from the slope traced by the star-forming galaxies. This angled distribution and the number density profile of galaxies in NUV−u space strongly suggest that the decline of the SFR in galaxies has to accelerate before they turn quiescent. We model this color–color distribution with a two-phase exponential decline star formation history. The models with an e-folding time in the second phase (the quenching phase) of 0.5 Gyr best fit the data. We further use the NUV−u number density profile to constrain the quenching rate among star-forming galaxies as a function of mass. Adopting an e-folding time of 0.5 Gyr in the second phase (or the quenching phase), we found the quenching rate to be 19%/Gyr, 25%/Gyr and 33%/Gyr for the three mass bins. These are upper limits of the quenching rate as the transition zone could also be populated by rejuvenated red-sequence galaxies.

An important source of opacity in exoplanet atmospheres at short visible and near-UV wavelengths is Rayleigh scattering of light on molecules. It is accompanied by a related, albeit weaker process—Raman scattering. We analyze the signatures of Raman scattering imprinted in the reflected light and the geometric albedo of exoplanets, which could provide information about atmospheric properties. Raman scattering affects the geometric albedo spectra of planets in the following ways. First, it causes filling-in of strong absorption lines in the incident radiation, thus producing sharp peaks in the albedo. Second, it shifts the wavelengths of spectral features in the reflected light causing the so-called Raman ghost lines. Raman scattering can also cause a broadband reduction of the albedo due to wavelength shifting of a stellar spectrum with red spectral index. Observing the Raman peaks in the albedo could be used to measure the column density of gas, thus providing constraints on the presence of clouds in the atmosphere. Observing the Raman ghost lines could be used to spectroscopically identify the main scatterer in the atmosphere, even molecules like H₂ or N₂, which do not have prominent spectral signatures in the optical wavelength range. If detected, ghost lines could also provide information about the temperature of the atmosphere. In this paper, we investigate the effects of Raman scattering in hydrogen- and nitrogen-dominated atmospheres. We analyze the feasibility of detecting the signatures of Raman scattering with the existing and future observational facilities, and of using these signatures as probes of exoplanetary atmospheres.

We incorporate our experimentally derived thermal rate coefficients for C + ${{\rm{H}}}_{3}^{+}$ forming CH⁺ and CH₂⁺ into a commonly used astrochemical model. We find that the Arrhenius–Kooij equation typically used in chemical models does not accurately fit our data and instead we use a more versatile fitting formula. At a temperature of 10 K and a density of 10⁴ cm⁻³, we find no significant differences in the predicted chemical abundances, but at higher temperatures of 50, 100, and 300 K we find up to factor of 2 changes. In addition, we find that the relatively small error on our thermal rate coefficients, ∼15%, significantly reduces the uncertainties on the predicted abundances compared to those obtained using the currently implemented Langevin rate coefficient with its estimated factor of 2 uncertainty.

Using a one-dimensional hybrid expanding box model, we investigate properties of the solar wind in the outer heliosphere. We assume a proton–electron plasma with a strictly transverse ambient magnetic field and, aside from the expansion, we take into account the influence of a continuous injection of cold pick-up protons through the charge-exchange process between the solar wind protons and hydrogen of interstellar origin. The injected cold pick-up protons form a ring distribution function, which rapidly becomes unstable, and generate Alfvén cyclotron waves. The Alfvén cyclotron waves scatter pick-up protons to a spherical shell distribution function that thickens over that time owing to the expansion-driven cooling. The Alfvén cyclotron waves heat solar wind protons in the perpendicular direction (with respect to the ambient magnetic field) through cyclotron resonance. At later times, the Alfvén cyclotron waves become parametrically unstable and the generated ion-acoustic waves heat protons in the parallel direction through Landau resonance. The resulting heating of the solar wind protons is efficient on the expansion timescale.

We report the direct imaging detection of a low-mass companion to a young, moderately active star V450 And, that was previously identified with the radial velocity (RV) method. The companion was found in high-contrast images obtained with the Subaru Telescope equipped with the HiCIAO camera and AO188 adaptive optics system. From the public ELODIE and SOPHIE archives we extracted available high-resolution spectra and RV measurements, along with RVs from the Lick planet search program. We combined our multi-epoch astrometry with these archival, partially unpublished RVs, and found that the companion is a low-mass star, not a brown dwarf, as previously suggested. We found the best-fitting dynamical masses to be ${m}_{1}={1.141}_{-0.091}^{+0.037}$ and ${m}_{2}={0.279}_{-0.020}^{+0.023}$${M}_{\odot }$. We also performed spectral analysis of the SOPHIE spectra with the iSpec code. Hipparcos time-series photometry shows a periodicity of P = 5.743 day, which is also seen in the SOPHIE spectra as an RV modulation of the star A. We interpret it as being caused by spots on the stellar surface, and the star to be rotating with the given period. From the rotation and level of activity, we found that the system is ${380}_{-100}^{+220}$ Myr old, consistent with an isochrone analysis (${220}_{-90}^{+2120}$ Myr). This work may serve as a test case for future studies of low-mass stars, brown dwarfs, and exoplanets by combination of RV and direct imaging data.

NASA's Kepler Mission uncovered a wealth of planetary systems, many with planets on short-period orbits. These short-period systems reside around 50% of Sun-like stars and are similarly prevalent around M dwarfs. Their formation and subsequent evolution is the subject of active debate. In this paper, we simulate late-stage, in situ planet formation across a grid of planetesimal disks with varying surface density profiles and total mass. We compare simulation results with observable characteristics of the Kepler sample. We identify mixture models with different primordial planetesimal disk properties that self-consistently recover the multiplicity, radius, period and period ratio, and duration ratio distributions of the Kepler planets. We draw three main conclusions. (1) We favor a "frozen-in" narrative for systems of short-period planets, in which they are stable over long timescales, as opposed to metastable. (2) The "Kepler dichotomy," an observed phenomenon of the Kepler sample wherein the architectures of planetary systems appear to either vary significantly or have multiple modes, can naturally be explained by formation within planetesimal disks with varying surface density profiles. Finally, (3) we quantify the nature of the "Kepler dichotomy" for both GK stars and M dwarfs, and find that it varies with stellar type. While the mode of planet formation that accounts for high multiplicity systems occurs in 24% ± 7% of planetary systems orbiting GK stars, it occurs in 63% ± 16% of planetary systems orbiting M dwarfs.

AE Aquarii (AE Aqr) is a cataclysmic binary hosting one of the fastest rotating (${P}_{\mathrm{spin}}$ = 33.08 s) white dwarfs (WDs) known. Based on seven years of Fermi Large Area Telescope (LAT) Pass 8 data, we report on a deep search for gamma-ray emission from AE Aqr. Using X-ray observations from ASCA, XMM-Newton, Chandra, Swift, Suzaku, and NuSTAR, spanning 20 years, we substantially extend and improve the spin ephemeris of AE Aqr. Using this ephemeris, we searched for gamma-ray pulsations at the spin period of the WD. No gamma-ray pulsations were detected above 3σ significance. Neither phase-averaged gamma-ray emission nor gamma-ray variability of AE Aqr is detected by Fermi LAT. We impose the most restrictive upper limit to the gamma-ray flux from AE Aqr to date: $1.3\times {10}^{-12}$ erg cm⁻² s⁻¹ in the 100 MeV–300 GeV energy range, providing constraints on models.

A sufficiently large temperature anisotropy can sometimes drive various types of electromagnetic plasma micro-instabilities, which can play an important role in the dynamics of relativistic pair plasmas in space, astrophysics, and laboratory environments. Here, we provide a detailed description of the cyclotron instability of parallel propagating electromagnetic waves in relativistic pair plasmas on the basis of a relativistic anisotropic distribution function. Using plasma kinetic theory and particle-in-cell simulations, we study the influence of the relativistic temperature and the temperature anisotropy on the collective and noncollective modes of these plasmas. Growth rates and dispersion curves from the linear theory show a good agreement with simulations results.

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850 μm continuum observations of the original Lyα Blob (LAB) in the SSA22 field at z = 3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with a total flux density of S₈₅₀ = 1.68 ± 0.06 mJy, corresponding to a star-formation rate of ∼150 M_⊙ yr⁻¹. The submillimeter sources are associated with several faint (m ≈ 27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (λ ≈ 5850 Å). One of these companions is spectroscopically confirmed with the Keck Multi-Object Spectrometer For Infra-Red Exploration to lie within 20 projected kpc and 250 km s⁻¹ of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high-resolution cosmological zoom simulation of a 10¹³M_⊙ halo at z = 3, including stellar, dust, and Lyα radiative transfer, we can model the ALMA+STIS observations and demonstrate that Lyα photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Lyα emission with similar surface brightness and morphology to observed giant LABs.

A planet having protective ozone within the collimated beam of a gamma-ray burst (GRB) may suffer ozone depletion, potentially causing a mass extinction event to existing life on a planet's surface and oceans. We model the dangers of long GRBs to planets in the Milky Way and utilize a static statistical model of the Galaxy, which matches major observable properties, such as the inside-out star formation history (SFH), metallicity evolution, and three-dimensional stellar number density distribution. The GRB formation rate is a function of both the SFH and metallicity. However, the extent to which chemical evolution reduces the GRB rate over time in the Milky Way is still an open question. Therefore, we compare the damaging effects of GRBs to biospheres in the Milky Way using two models. One model generates GRBs as a function of the inside-out SFH. The other model follows the SFH, but generates GRB progenitors as a function of metallicity, thereby favoring metal-poor host regions of the Galaxy over time. If the GRB rate only follows the SFH, the majority of the GRBs occur in the inner Galaxy. However, if GRB progenitors are constrained to low-metallicity environments, then GRBs only form in the metal-poor outskirts at recent epochs. Interestingly, over the past 1 Gyr, the surface density of stars (and their corresponding planets), which survive a GRB is still greatest in the inner galaxy in both models. The present-day danger of long GRBs to life at the solar radius (R_⊙ = 8 kpc) is low. We find that at least ∼65% of stars survive a GRB over the past 1 Gyr. Furthermore, when the GRB rate was expected to have been enhanced at higher redshifts, such as z ≳ 0.5, our results suggest that a large fraction of planets would have survived these lethal GRB events.

A galaxy group catalog is constructed from the 2MASS Redshift Survey (2MRS) with the use of a halo-based group finder. The halo mass associated with a group is estimated using a "GAP" method based on the luminosity of the central galaxy and its gap with other member galaxies. Tests using mock samples show that this method is reliable, particularly for poor systems containing only a few members. On average, 80% of all the groups have completeness $\gt 0.8$, and about 65% of the groups have zero contamination. Halo masses are estimated with a typical uncertainty of $\sim 0.35\,\mathrm{dex}$. The application of the group finder to the 2MRS gives 29,904 groups from a total of 43,246 galaxies at $z\leqslant 0.08$, with 5286 groups having two or more members. Some basic properties of this group catalog is presented, and comparisons are made with other group catalogs in overlap regions. With a depth to $z\sim 0.08$ and uniformly covering about 91% of the whole sky, this group catalog provides a useful database to study galaxies in the local cosmic web, and to reconstruct the mass distribution in the local universe.

We perform two-dimensional axially symmetric radiation hydrodynamic simulations to assess the impact of outflows and radiative force feedback from massive protostars by varying when the protostellar outflow starts, and to determine the ratio of ejection to accretion rates and the strength of the wide-angle disk wind component. The star-formation efficiency, i.e., the ratio of final stellar mass to initial core mass, is dominated by radiative forces and the ratio of outflow to accretion rates. Increasing this ratio has three effects. First, the protostar grows slower with a lower luminosity at any given time, lowering radiative feedback. Second, bipolar cavities cleared by the outflow become larger, further diminishing radiative feedback on disk and core scales. Third, the higher momentum outflow sweeps up more material from the collapsing envelope, decreasing the protostar's potential mass reservoir via entrainment. The star-formation efficiency varies with the ratio of ejection to accretion rates from 50% in the case of very weak outflows to as low as 20% for very strong outflows. At latitudes between the low-density bipolar cavity and the high-density accretion disk, wide-angle disk winds remove some of the gas, which otherwise would be part of the accretion flow onto the disk; varying the strength of these wide-angle disk winds, however, alters the final star-formation efficiency by only ±6%. For all cases, the opening angle of the bipolar outflow cavity remains below 20° during early protostellar accretion phases, increasing rapidly up to 65° at the onset of radiation pressure feedback.

The composition of a planet's atmosphere is determined by its formation, evolution, and present-day insolation. A planet's spectrum therefore may hold clues on its origins. We present a "chain" of models, linking the formation of a planet to its observable present-day spectrum. The chain links include (1) the planet's formation and migration, (2) its long-term thermodynamic evolution, (3) a variety of disk chemistry models, (4) a non-gray atmospheric model, and (5) a radiometric model to obtain simulated spectroscopic observations with James Webb Space Telescope and ARIEL. In our standard chemistry model the inner disk is depleted in refractory carbon as in the Solar System and in white dwarfs polluted by extrasolar planetesimals. Our main findings are: (1) envelope enrichment by planetesimal impacts during formation dominates the final planetary atmospheric composition of hot Jupiters. We investigate two, under this finding, prototypical formation pathways: a formation inside or outside the water iceline, called "dry" and "wet" planets, respectively. (2) Both the "dry" and "wet" planets are oxygen-rich (C/O < 1) due to the oxygen-rich nature of the solid building blocks. The "dry" planet's C/O ratio is <0.2 for standard carbon depletion, while the "wet" planet has typical C/O values between 0.1 and 0.5 depending mainly on the clathrate formation efficiency. Only non-standard disk chemistries without carbon depletion lead to carbon-rich C/O ratios >1 for the "dry" planet. (3) While we consistently find C/O ratios <1, they still vary significantly. To link a formation history to a specific C/O, a better understanding of the disk chemistry is thus needed.

Previous detections of methyl and ethyl formate make other small substituted formates potential candidates for observation in the interstellar medium. Among them, vinyl formate is one of the simplest unsaturated carboxylic ester. The aim of this work is to provide direct experimental frequencies of the ground vibrational state of vinyl formate in a large spectral range for astrophysical use. The room-temperature rotational spectrum of vinyl formate has been measured from 80 to 360 GHz and analyzed in terms of Watson's semirigid rotor Hamiltonian. Two thousand six hundred transitions within J = 3–88 and K_a = 0–28 were assigned to the most stable conformer of vinyl formate and a new set of spectroscopic constants was accurately determined. Spectral features of vinyl formate were then searched for in Orion KL, Sgr B2(N), B1-b, and TMC-1 molecular clouds. Upper limits to the column density of vinyl formate are provided.

Massive stars are key players in the evolution of galaxies, yet their formation pathway remains unclear. In this work, we use data from several galaxy-wide surveys to build an unbiased data set of ∼600 massive young stellar objects, ∼200 giant molecular clouds (GMCs), and ∼100 young (<10 Myr) optical stellar clusters (SCs) in the Large Magellanic Cloud. We employ this data to quantitatively study the location and clustering of massive star formation and its relation to the internal structure of GMCs. We reveal that massive stars do not typically form at the highest column densities nor centers of their parent GMCs at the ∼6 pc resolution of our observations. Massive star formation clusters over multiple generations and on size scales much smaller than the size of the parent GMC. We find that massive star formation is significantly boosted in clouds near SCs. However, whether a cloud is associated with an SC does not depend on either the cloud's mass or global surface density. These results reveal a connection between different generations of massive stars on timescales up to 10 Myr. We compare our work with Galactic studies and discuss our findings in terms of GMC collapse, triggered star formation, and a potential dichotomy between low- and high-mass star formation.

The neutron star ocean is a plasma of ions and electrons that extends from the base of the neutron star's envelope to a depth where the plasma crystallizes into a solid crust. During an accretion outburst in an X-ray transient, material accumulates in the envelope of the neutron star primary. This accumulation compresses the neutron star's outer layers and induces nuclear reactions in the ocean and crust. Accretion-driven heating raises the ocean's temperature and increases the frequencies of $g$-modes in the ocean; when accretion halts, the ocean cools and ocean $g$-mode frequencies decrease. If the observed low-frequency quasi-periodic oscillations on accreting neutron stars are $g$-modes in the ocean, the observed quasi-periodic oscillation frequencies will increase during the outburst—reaching a maximum when the ocean temperature reaches steady state—and subsequently decrease during quiescence. For time-averaged accretion rates during outbursts between $\langle \dot{M}\rangle =0.1\mbox{--}1.0\,{\dot{M}}_{\mathrm{Edd}}$ the predicted $g$-mode fundamental n = 1 l = 2 frequency is between ≈3–7 Hz for slowly rotating neutron stars. Accreting neutron stars that require extra shallow heating, such as the Z-sources MAXI J0556-332, MXB 1659-29, and XTE J1701-462, have predicted $g$-mode fundamental frequencies between ≈3–16 HZ. Therefore, observations of low-frequency quasi-periodic oscillations between $\approx 8\mbox{--}16\,\mathrm{Hz}$ in these sources, or in other transients that require shallow heating, will support a $g$-mode origin for the observed quasi-periodic oscillations.

We report on the results from a large observational campaign on the bare Seyfert galaxy Ark 120, jointly carried out in 2014 with XMM-Newton, Chandra, and NuSTAR. The fortunate line of sight to this source, devoid of any significant absorbing material, provides an incomparably clean view to the nuclear regions of an active galaxy. Here we focus on the analysis of the iron fluorescence features, which form a composite emission pattern in the 6–7 keV band. The prominent Kα line from neutral iron at 6.4 keV is resolved in the Chandra High-Energy Transmission Grating spectrum to a full-width at half maximum of ${4700}_{-1500}^{+2700}$ km s⁻¹, consistent with an origin from the optical broad-line region. Excess components are detected on both sides of the narrow Kα line: the red one (6.0–6.3 keV) clearly varies in strength in about one year, and hints at the presence of a broad, mildly asymmetric line from the accretion disk; the blue one (6.5–7.0 keV), instead, is likely a blend of different contributions, and appears to be constant when integrated over long enough exposures. However, the Fe K excess emission map computed over the 7.5 days of the XMM-Newton monitoring shows that both the red and blue features are actually highly variable on timescales of ∼10–15 hr, suggesting that they might arise from short-lived hotspots on the disk surface, located within a few tens of gravitational radii from the central supermassive black hole and possibly illuminated by magnetic reconnection events. Any alternative explanation would still require a highly dynamic, inhomogeneous disk/coronal system, involving clumpiness and/or instability.

We propose to use the flux variability of lensed quasar images induced by gravitational microlensing to measure the transverse peculiar velocity of lens galaxies over a wide range of redshift. Microlensing variability is caused by the motions of the observer, the lens galaxy (including the motion of the stars within the galaxy), and the source. Hence, its frequency is directly related to the galaxy's transverse peculiar velocity. The idea is to count time-event rates (e.g., peak or caustic crossing rates) in the observed microlensing light curves of lensed quasars that can be compared with model predictions for different values of the transverse peculiar velocity. To compensate for the large timescale of microlensing variability, we propose to count and model the number of events in an ensemble of gravitational lenses. We develop the methodology to achieve this goal and apply it to an ensemble of 17 lensed quasar systems. In spite of the shortcomings of the available data, we have obtained tentative estimates of the peculiar velocity dispersion of lens galaxies at z ∼ 0.5, ${\sigma }_{\mathrm{pec}}(0.53\pm 0.18)\simeq (638\pm 213)\sqrt{\langle m\rangle /0.3\,{M}_{\odot }}\,\mathrm{km}\,{{\rm{s}}}^{-1}$. Scaling at zero redshift, we derive ${\sigma }_{\mathrm{pec}}(0)\simeq (491\pm 164)\sqrt{\langle m\rangle /0.3\,{M}_{\odot }}\,\mathrm{km}\,{{\rm{s}}}^{-1}$, consistent with peculiar motions of nearby galaxies and with recent N-body nonlinear reconstructions of the Local Universe based on ΛCDM. We analyze the different sources of uncertainty of the method and find that for the present ensemble of 17 lensed systems the error is dominated by Poisson noise, but that for larger ensembles the impact of the uncertainty on the average stellar mass may be significant.

The detection of periodicity in the broadband non-thermal emission of blazars has so far been proven to be elusive. However, there are a number of scenarios that could lead to quasi-periodic variations in blazar light curves. For example, an orbital or thermal/viscous period of accreting matter around central supermassive black holes could, in principle, be imprinted in the multi-wavelength emission of small-scale blazar jets, carrying such crucial information about plasma conditions within the jet launching regions. In this paper, we present the results of our time series analysis of the ∼9.2 yr long, and exceptionally well-sampled, optical light curve of the BL Lac object OJ 287. The study primarily used the data from our own observations performed at the Mt. Suhora and Kraków Observatories in Poland, and at the Athens Observatory in Greece. Additionally, SMARTS observations were used to fill some of the gaps in the data. The Lomb–Scargle periodogram and the weighted wavelet Z-transform methods were employed to search for possible quasi-periodic oscillations in the resulting optical light curve of the source. Both methods consistently yielded a possible quasi-periodic signal around the periods of ∼400 and ∼800 days, the former with a significance (over the underlying colored noise) of $\geqslant 99 \%$. A number of likely explanations for this are discussed, with preference given to a modulation of the jet production efficiency by highly magnetized accretion disks. This supports previous findings and the interpretation reported recently in the literature for OJ 287 and other blazar sources.

We used ultra-deep J and K_s images secured with the near-infrared (NIR) GSAOI camera assisted by the multi-conjugate adaptive optics system GeMS at the GEMINI South Telescope in Chile, to obtain a (K_s, J − K_s) color–magnitude diagram (CMD) for the bulge globular cluster NGC 6624. We obtained the deepest and most accurate NIR CMD from the ground for this cluster, by reaching K_s ∼ 21.5, approximately 8 mag below the horizontal branch level. The entire extension of the Main Sequence (MS) is nicely sampled and at K_s ∼ 20 we detected the so-called MS "knee" in a purely NIR CMD. By taking advantage of the exquisite quality of the data, we estimated the absolute age of NGC 6624 (t_age = 12.0 ± 0.5 Gyr), which turns out to be in good agreement with previous studies in the literature. We also analyzed the luminosity and mass functions of MS stars down to M ∼ 0.45 M_⊙, finding evidence of a significant increase of low-mass stars at increasing distances from the cluster center. This is a clear signature of mass segregation, confirming that NGC 6624 is in an advanced stage of dynamical evolution.

We present the results of a new spectroscopic survey for dusty intervening absorption systems, particularly damped Lyα absorbers (DLAs), toward reddened quasars. The candidate quasars are selected from mid-infrared photometry from the Wide-field Infrared Survey Explorer combined with optical and near-infrared photometry. Out of 1073 candidates, we secure low-resolution spectra for 108 using the Nordic Optical Telescope on La Palma, Spain. Based on the spectra, we are able to classify 100 of the 108 targets as quasars. A large fraction (50%) is observed to have broad absorption lines (BALs). Moreover, we find six quasars with strange breaks in their spectra, which are not consistent with regular dust reddening. Using template fitting, we infer the amount of reddening along each line of sight ranging from A(V) ≈ 0.1 to 1.2 mag (assuming a Small Magellanic Cloud extinction curve). In four cases, the reddening is consistent with dust exhibiting the 2175 Å feature caused by an intervening absorber, and for two of these, an Mg ii absorption system is observed at the best-fit absorption redshift. In the rest of the cases, the reddening is most likely intrinsic to the quasar. We observe no evidence for dusty DLAs in this survey. However, the large fraction of BAL quasars hampers the detection of absorption systems. Out of the 50 non-BAL quasars, only 28 have sufficiently high redshift to detect Lyα in absorption.

We present the results of a search for new circumstellar disks around low-mass stars and brown dwarfs with spectral types >K5 that are confirmed or candidate members of nearby young moving groups. Our search input sample was drawn from the BANYAN surveys of Malo et al. and Gagné et al. Two Micron All-Sky Survey and Wide-field Infrared Survey Explorer data were used to detect near- to mid-infrared excesses that would reveal the presence of circumstellar disks. A total of 13 targets with convincing excesses were identified: 4 are new and 9 were already known in the literature. The new candidates are 2MASS J05010082–4337102 (M4.5), J08561384–1342242 (${\rm{M}}8\gamma$), J12474428–3816464 (${\rm{M}}9\gamma$), and J02265658–5327032 (${\rm{L}}0\delta$); they are candidate members of the TW Hya ($\sim 10\pm 3$ Myr), Columba (∼${42}_{-4}^{+6}\,$Myr), and Tucana-Horologium ($\sim 45\pm 4$ Myr) associations, with masses of 120 and 13–18 ${M}_{\mathrm{Jup}}$. The M8–L0 objects in Columba and Tucana-Horologium are potentially among the first substellar disk systems aged ∼40 Myr. Estimates of the new candidates' mean disk temperatures and fractional luminosities are in the ranges ∼$135\mbox{--}520\,{\rm{K}}$ and $0.021\mbox{--}0.15$, respectively. New optical spectroscopy of J0501–4337 reveals strong Hα emission, possibly indicating ongoing accretion, provides a detection of lithium absorption, and shows a radial velocity measurement that is consistent with a membership to Columba. We also present a near-infrared spectrum of J0226–5327 that reveals Paschen β emission and shows signs of low surface gravity, consistent with accretion from a disk and a young age.

Low-resolution Spitzer spectral map data (>1700 spectra) of ten reflection nebulae (RNe) fields are analyzed using the data and tools available through the NASA Ames PAH IR Spectroscopic Database. The PAH emission is broken down into PAH charge state using a database fitting approach. Here, the physics of the PAH emission process is taken into account and uses target appropriate parameters, e.g., a stellar radiation model for the exciting star. The breakdown results are combined with results derived using the traditional PAH band strength approach, which interprets particular PAH band strength ratios as proxies for the PAH charge state, e.g., the 6.2/11.2 μm PAH band strength ratio. These are successfully calibrated against their database equivalent; the PAH ionized fraction (f_i). The PAH ionized fraction is converted into the PAH ionization parameter, which relates the PAH ionized fraction to the strength of the radiation field, gas temperature and electron density. The behavior of the 12.7 μm PAH band is evaluated as a tracer for PAH ionization and erosion. The plot of the 8.6 versus 11.2 μm PAH band strength for the northwest photo-dominated region (PDR) in NGC 7023 is shown to be a robust diagnostic template for the PAH ionized fraction. Remarkably, most of the other RNe fall within the limits set by NGC 7023. Finally, PAH spectroscopic templates are constructed and verified as principal components. Template spectra derived from NGC 7023 and NGC 2023 compare extremely well with each other, with those derived for NGC 7023 successfully reproducing the PAH emission observed from NGC 2023.

We present the first results from an ongoing survey to characterize the circumgalactic medium (CGM) of massive high-redshift galaxies detected as submillimeter galaxies (SMGs). We constructed a parent sample of 163 SMG–QSO pairs with separations less than ∼36'' by cross-matching far-infrared-selected galaxies from Herschel with spectroscopically confirmed QSOs. The Herschel sources were selected to match the properties of the SMGs. We determined the sub-arcsecond positions of six Herschel sources with the Very Large Array and obtained secure redshift identification for three of those with near-infrared spectroscopy. The QSO sightlines probe transverse proper distances of 112, 157, and 198 kpc at foreground redshifts of 2.043, 2.515, and 2.184, respectively, which are comparable to the virial radius of the ∼10¹³M_⊙ halos expected to host SMGs. High-quality absorption-line spectroscopy of the QSOs reveals systematically strong H i Lyα absorption around all three SMGs, with rest-frame equivalent widths of ∼2–3 Å. However, none of the three absorbers exhibit compelling evidence for optically thick H i gas or metal absorption, in contrast to the dominance of strong neutral absorbers in the CGM of luminous z ∼ 2 QSOs. The low covering factor of optically thick H i gas around SMGs tentatively indicates that SMGs may not have as prominent cool gas reservoirs in their halos as the coeval QSOs and that they may inhabit less massive halos than previously thought.

Recent observations of the photosphere using high spatial and temporal resolution show small dynamic features at or below the current resolving limits. A new pixel dynamics method has been developed to analyze spectral profiles and quantify changes in line displacement, width, asymmetry, and peakedness of photospheric absorption lines. The algorithm evaluates variations of line profile properties in each pixel and determines the statistics of such fluctuations averaged over all pixels in a given region. The method has been used to derive statistical characteristics of pixel fluctuations in observed quiet-Sun regions, an active region with no eruption, and an active region with an ongoing eruption. Using Stokes I images from the Vector Spectromagnetograph (VSM) of the Synoptic Optical Long-term Investigations of the Sun (SOLIS) telescope on 2012 March 13, variations in line width and peakedness of Fe i 6301.5 Å are shown to have a distinct spatial and temporal relationship with an M7.9 X-ray flare in NOAA 11429. This relationship is observed as stationary and contiguous patches of pixels adjacent to a sunspot exhibiting intense flattening in the line profile and line-center displacement as the X-ray flare approaches peak intensity, which is not present in area scans of the non-eruptive active region. The analysis of pixel dynamics allows one to extract quantitative information on differences in plasma dynamics on sub-pixel scales in these photospheric regions. The analysis can be extended to include the Stokes parameters and study signatures of vector components of magnetic fields and coupled plasma properties.

Terrestrial exoplanets in the canonical habitable zone may have a variety of initial water fractions due to random volatile delivery by planetesimals. If the total planetary water complement is high, the entire surface may be covered in water, forming a "waterworld." On a planet with active tectonics, competing mechanisms act to regulate the abundance of water on the surface by determining the partitioning of water between interior and surface. Here we explore how the incorporation of different mechanisms for the degassing and regassing of water changes the volatile evolution of a planet. For all of the models considered, volatile cycling reaches an approximate steady state after $\sim 2\ \mathrm{Gyr}$. Using these steady states, we find that if volatile cycling is either solely dependent on temperature or seafloor pressure, exoplanets require a high abundance ($\gtrsim 0.3 \%$ of total mass) of water to have fully inundated surfaces. However, if degassing is more dependent on seafloor pressure and regassing mainly dependent on mantle temperature, the degassing rate is relatively large at late times and a steady state between degassing and regassing is reached with a substantial surface water fraction. If this hybrid model is physical, super-Earths with a total water fraction similar to that of the Earth can become waterworlds. As a result, further understanding of the processes that drive volatile cycling on terrestrial planets is needed to determine the water fraction at which they are likely to become waterworlds.

The origin of hydrogen in chondritic components is poorly understood. Their isotopic composition is heavier than the solar nebula gas. In addition, in most meteorites, hydrous silicates are found to be lighter than the coexisting organic matter. Ionizing irradiation recently emerged as an efficient hydrogen fractionating process in organics, but its effect on H-bearing silicates remains essentially unknown. We report the evolution of the D/H of hydrous silicates experimentally irradiated by electrons. Thin films of amorphous silica, amorphous "serpentine," and pellets of crystalline muscovite were irradiated at 4 and 30 keV. For all samples, irradiation leads to a large hydrogen loss correlated with a moderate deuterium enrichment of the solid residue. The entire data set can be described by a Rayleigh distillation. The calculated fractionation factor is consistent with a kinetically controlled fractionation during the loss of hydrogen. Furthermore, for a given ionizing condition, the deuteration of the silicate residues is much lower than the deuteration measured on irradiated organic macromolecules. These results provide firm evidence of the limitations of ionizing irradiation as a driving mechanism for D-enrichment of silicate materials. The isotopic composition of the silicate dust cannot rise from a protosolar to a chondritic signature during solar irradiations. More importantly, these results imply that irradiation of the disk naturally induces a strong decoupling of the isotopic signatures of coexisting organics and silicates. This decoupling is consistent with the systematic difference observed between the heavy organic matter and the lighter water typically associated with minerals in the matrix of most carbonaceous chondrites.

We present a robust measurement of the rest-frame UV luminosity function (LF) and its evolution during the peak epoch of cosmic star formation at $1\lt z\lt 3$. We use our deep near-ultraviolet imaging from WFC3/UVIS on the Hubble Space Telescope and existing Advanced Camera for Surveys (ACS)/WFC and WFC3/IR imaging of three lensing galaxy clusters, Abell 2744 and MACS J0717 from the Hubble Frontier Field survey and Abell 1689. Combining deep UV imaging and high magnification from strong gravitational lensing, we use photometric redshifts to identify 780 ultra-faint galaxies with ${M}_{\mathrm{UV}}\lt -12.5$ AB mag at $1\lt z\lt 3$. From these samples, we identified five new, faint, multiply imaged systems in A1689. We run a Monte Carlo simulation to estimate the completeness correction and effective volume for each cluster using the latest published lensing models. We compute the rest-frame UV LF and find the best-fit faint-end slopes of $\alpha =-1.56\pm 0.04$, $\alpha =-1.72\pm 0.04$, and $\alpha =-1.94\pm 0.06$ at $1.0\lt z\lt 1.6$, $1.6\lt z\lt 2.2$, and $2.2\lt z\lt 3.0$, respectively. Our results demonstrate that the UV LF becomes steeper from $z\sim 1.3$ to $z\sim 2.6$ with no sign of a turnover down to ${M}_{\mathrm{UV}}=-14$ AB mag. We further derive the UV LFs using the Lyman break "dropout" selection and confirm the robustness of our conclusions against different selection methodologies. Because the sample sizes are so large and extend to such faint luminosities, the statistical uncertainties are quite small, and systematic uncertainties (due to the assumed size distribution, for example) likely dominate. If we restrict our analysis to galaxies and volumes above $\gt 50 \%$ completeness in order to minimize these systematics, we still find that the faint-end slope is steep and getting steeper with redshift, though with slightly shallower (less negative) values ($\alpha =-1.55\pm 0.06$, −1.69 ± 0.07, and −1.79 ± 0.08 for $z\sim 1.3$, 1.9, and 2.6, respectively). Finally, we conclude that the faint star-forming galaxies with UV magnitudes of $-18.5\lt {M}_{\mathrm{UV}}\lt -12.5$ covered in this study produce the majority (55%–60%) of the unobscured UV luminosity density at $1\lt z\lt 3$.

Intermittency of heating in weakly collisional plasma turbulence is an active subject of research, with significant potential impact on understanding of the solar wind, solar corona, and astrophysical plasmas. Recent studies suggest a role of vorticity in plasma heating. In magnetohydrodynamics small-scale vorticity is generated near current sheets and this effect persists in kinetic plasma, as demonstrated here with hybrid and fully kinetic particle-in-cell simulations. Furthermore, vorticity enhances local kinetic effects, with a generalized resonance condition selecting sign-dependent enhancements or reductions of proton heating and thermal anisotropy. In such plasmas heating is correlated with vorticity and current density, but more strongly with vorticity. These results help explain several prior results that find kinetic effects and energization near to, but not centered on, current sheets. Evidently intermittency in kinetic plasma involves multiple physical quantities, and the associated coherent structures and nonthermal effects are closely related.

Because WISE J085510.83−071442.5 (hereafter WISE 0855-0714) is the coldest known brown dwarf (∼250 K) and one of the Sun's closest neighbors (2.2 pc), it offers a unique opportunity to study a planet-like atmosphere in an unexplored regime of temperature. To detect and characterize inhomogeneities in its atmosphere (e.g., patchy clouds, hot spots), we have performed time-series photometric monitoring of WISE 0855-0714 at 3.6 and 4.5 μm with the Spitzer Space Telescope during two 23 hr periods that were separated by several months. For both bands, we have detected variability with peak-to-peak amplitudes of 4%–5% and 3%–4% in the first and second epochs, respectively. The light curves are semiperiodic in the first epoch for both bands, but they are more irregular in the second epoch. Models of patchy clouds have predicted a large increase in mid-infrared (mid-IR) variability amplitudes (for a given cloud covering fraction) with the appearance of water ice clouds at T_eff < 375 K, so if such clouds are responsible for the variability of WISE 0855-0714, then its small amplitudes of variability indicate a very small deviation in cloud coverage between hemispheres. Alternatively, the similarity in mid-IR variability amplitudes between WISE 0855-0714 and somewhat warmer T and Y dwarfs may suggest that they share a common origin for their variability (i.e., not water clouds). In addition to our variability data, we have examined other constraints on the presence of water ice clouds in the atmosphere of WISE 0855-0714, including the recent mid-IR spectrum from Skemer et al. (2016). We find that robust evidence of such clouds is not yet available.

We report on low-frequency radio (85–35 MHz) spectral observations of four different type II radio bursts, which exhibited fundamental-harmonic emission and split-band structure. Each of the bursts was found to be closely associated with a whitelight coronal mass ejection (CME) close to the Sun. We estimated the coronal magnetic field strength from the split-band characteristics of the bursts, by assuming a model for the coronal electron density distribution. The choice of the model was constrained, based on the following criteria: (1) when the radio burst is observed simultaneously in the upper and lower bands of the fundamental component, the location of the plasma level corresponding to the frequency of the burst in the lower band should be consistent with the deprojected location of the leading edge (LE) of the associated CME; (2) the drift speed of the type II bursts derived from such a model should agree closely with the deprojected speed of the LE of the corresponding CMEs. With the above conditions, we find that: (1) the estimated field strengths are unique to each type II burst, and (2) the radial variation of the field strength in the different events indicate a pattern. It is steepest for the case where the heliocentric distance range over which the associated burst is observed is closest to the Sun, and vice versa.

We report on a search for radio transients at 340 MHz with the Jansky Very Large Array Low-band Ionosphere and Transient Experiment (VLITE). Between 2015 July 29 and September 27, operating in commensal mode, VLITE imaged approximately 2800 pointings covering 12,000 deg² on the sky, sampling timescales ranging from tens of seconds to several hours on a daily basis. In addition, between 2015 February 25 and May 9, VLITE observed 55 epochs of roughly 2–4 hr each toward the COSMOS field. Using existing radio source catalogs, we have searched all of the daily VLITE images for transients, while for the COSMOS field we compared individual images and the summed image to search for new sources in repeated observations of the same field. The wide range of timescales makes VLITE sensitive to both coherent and incoherent transient source classes. No new transients are found, allowing us to set stringent upper limits on transients at milli-jansky levels and at low frequencies where comparatively few such surveys have been carried out to date. An all-sky isotropic surface density of bursting radio transients with similar rates, durations, and intensities as the unusual transient GCRT J1745−3009, discovered in wide-field monitoring toward the Galactic center, is ruled out with high confidence. The resulting non-detections allows us to argue that this is a coherent source, whose properties most resemble the growing class of nulling pulsars. We end with a discussion of the future prospects for the detection of transients by VLITE and other experiments.

Interstellar pickup ions (PUIs) play a significant part in mediating the solar wind (SW) interaction with the interstellar medium. In this paper, we examine the details of spatial variation of the PUI velocity distribution function (VDF) in the SW by solving the PUI transport equation. We assume the PUI distribution is isotropic resulting from strong pitch-angle scattering by wave–particle interaction. A three-dimensional model combining the MHD treatment of the background SW and neutrals with a kinetic treatment of PUIs throughout the heliosphere and the surrounding local interstellar medium has been developed. The model generates PUI power-law tails via second-order Fermi process. We analyze how PUIs transform across the heliospheric termination shock and obtain the PUI phase space distribution in the inner heliosheath including continuing velocity diffusion. Our simulated PUI spectra are compared with observations made by New Horizons, Ulysses, Voyager 1, 2, and Cassini, and a satisfactory agreement is demonstrated. Some specific features in the observations, for example, a cutoff of PUI VDF at v = V_SW and a f ∝ v⁻⁵ tail in the reference frame of the SW, are well represented by the model.

We present optical and near-infrared spectroscopy of WISEA J061543.91−124726.8, which we rediscovered as a high motion object in the AllWISE survey. The spectra of this object are unusual; while the red optical (λ > 7000 Å) and near-infrared spectra exhibit characteristic TiO, VO, and H₂O bands of a late-M dwarf, the blue portion of its optical spectrum shows a significant excess of emission relative to late-M-type templates. The excess emission is relatively featureless, with the exception of a prominent and very broad Na i D doublet. We find that no single, ordinary star can reproduce these spectral characteristics. The most likely explanation is an unresolved binary system of an M7 dwarf and a cool white dwarf. The flux of a cool white dwarf drops in the optical red and near-infrared, due to collision-induced absorption, thus allowing the flux of a late-M dwarf to show through. This scenario, however, does not explain the Na D feature, which is unlike that of any known white dwarf, but which could perhaps be explained via unusual abundance or pressure conditions.

The transport of energetic electrons in a solar flare is modeled using a time-dependent one-dimensional Fokker–Planck code that incorporates asymmetric magnetic convergence. We derive the temporal and spectral evolution of the resulting hard X-ray (HXR) emission in the conjugate chromospheric footpoints, assuming thick target photon production, and characterize the time evolution of the numerically simulated footpoint asymmetry and its relationship to the photospheric magnetic configuration. The thick target HXR asymmetry in the conjugate footpoints is found to increase with magnetic field ratio as expected. However, we find that the footpoint HXR asymmetry saturates for conjugate footpoint magnetic field ratios ≥4. This result is borne out in a direct comparison with observations of 44 double-footpoint flares. The presence of such a limit has not been reported before, and may serve as both a theoretical and observational benchmark for testing a range of particle transport and flare morphology constraints, particularly as a means to differentiate between isotropic and anisotropic particle injection.

The double adiabatic expansion of the nearly collisionless solar wind plasma creates conditions for the firehose instability to develop and efficiently prevent the further increase of the plasma temperature in the direction parallel to the interplanetary magnetic field. The conditions imposed by the firehose instability have been extensively studied using idealized approaches that ignore the mutual effects of electrons and protons. Recently, more realistic approaches have been proposed that take into account the interplay between electrons and protons,  unveiling new regimes of the parallel oscillatory modes. However, for oblique wave propagation the instability develops distinct branches that grow much faster and may therefore be more efficient than the parallel firehose instability in constraining the temperature anisotropy of the plasma particles. This paper reports for the first time on the effects of electron plasma properties on the oblique proton firehose (PFH) instability and provides a comprehensive vision of the entire unstable wave-vector spectrum, unifying the proton and the smaller electron scales. The plasma β and temperature anisotropy regimes considered here are specific for the solar wind and magnetospheric conditions, and enable the electrons and protons to interact via the excited electromagnetic fluctuations. For the selected parameters, simultaneous electron and PFH instabilities can be observed with a dispersion spectrum of the electron firehose (EFH) extending toward the proton scales. Growth rates of the PFH instability are markedly boosted by the anisotropic electrons, especially in the oblique direction where the EFH growth rates are orders of magnitude higher.

In this paper, we report our multiwavelength observations of the C3.1 circular-ribbon flare SOL2015-10-16T10:20 in active region (AR) 12434. The flare consisted of a circular flare ribbon (CFR), an inner flare ribbon (IFR) inside it, and a pair of short parallel flare ribbons (PFRs). The PFRs located to the north of the IFR were most striking in the Interface Region Imaging Spectrograph (IRIS) 1400 and 2796 Å images. For the first time, we observed the circular-ribbon flare in the Ca ii H line of the Solar Optical Telescope on board Hinode, which has a similar shape as observed in the Atmospheric Imaging Assembly 1600 Å on board the Solar Dynamic Observatory (SDO). Photospheric line-of-sight magnetograms from the Helioseismic and Magnetic Imager on board SDO show that the flare is associated with positive polarities with a negative polarity inside. The IFR and CFR were cospatial with the negative polarity and positive polarities, implying the existence of a magnetic null point (${\boldsymbol{B}}=0$) and a dome-like spine–fan topology. During the impulsive phase of the flare, "two-step" raster observations of IRIS with a cadence of 6 s and an exposure time of 2 s showed plasma downflow at the CFR in the Si ivλ1402.77 line ($\mathrm{log}T\approx 4.8$), suggesting chromospheric condensation. The downflow speeds first increased rapidly from a few km s⁻¹ to the peak values of 45–52 km s⁻¹, before decreasing gradually to the initial levels. The decay timescales of condensation were 3–4 minutes, indicating ongoing magnetic reconnection. Interestingly, the downflow speeds are positively correlated with the logarithm of the Si iv line intensity and time derivative of the GOES soft X-ray (SXR) flux in 1–8 Å. The radio dynamic spectra are characterized by a type III radio burst associated with the flare, which implies that the chromospheric condensation was most probably driven by nonthermal electrons. Using an analytical expression and the peak Doppler velocity, we derive the lower limit of energy flux of the precipitating electrons, i.e., 0.65 × 10¹⁰ erg cm⁻² s⁻¹. The Si iv line intensity and SXR derivative show quasi-periodic pulsations with periods of 32–42 s, which are likely caused by intermittent null-point magnetic reconnections modulated by the fast wave propagating along the fan surface loops at a phase speed of 950–1250 km s⁻¹. Periodic accelerations and precipitations of the electrons result in periodic heating observed in the Si iv line and SXR.

We examine the long-term time evolution (1965–2015) of the relationships between solar wind proton temperature (T_p) and speed (V_p) and between the proton density (n_p) and speed using OMNI solar wind observations taken near Earth. We find a long-term decrease in the proton temperature–speed (T_p–V_p) slope that lasted from 1972 to 2010, but has been trending upward since 2010. Since the solar wind proton density–speed (n_p–V_p) relationship is not linear like the T_p–V_p relationship, we perform power-law fits for n_p–V_p. The exponent (steepness in the n_p–V_p relationship) is correlated with the solar cycle. This exponent has a stronger correlation with current sheet tilt angle than with sunspot number because the sunspot number maxima vary considerably from cycle to cycle and the tilt angle maxima do not. To understand this finding, we examined the average n_p for different speed ranges, and found that for the slow wind n_p is highly correlated with the sunspot number, with a lag of approximately four years. The fast wind n_p variation was less, but in phase with the cycle. This phase difference may contribute to the n_p–V_p exponent correlation with the solar cycle. These long-term trends are important since empirical formulas based on fits to T_p and V_p data are commonly used to identify interplanetary coronal mass ejections, but these formulas do not include any time dependence. Changes in the solar wind density over a solar cycle will create corresponding changes in the near-Earth space environment and the overall extent of the heliosphere.

We present optical long-slit observations of the complete sample of 71 Type 2 active galactic nuclei (AGNs) with double-peaked narrow emission lines at z < 0.1 in the Sloan Digital Sky Survey. Double-peaked emission lines are produced by a variety of mechanisms including disk rotation, kiloparsec-scale dual AGNs, and narrow-line region (NLR) kinematics (outflows or inflows). We develop a novel kinematic classification technique to determine the nature of these objects using long-slit spectroscopy alone. We determine that 86% of the double-peaked profiles are produced by moderate-luminosity AGN outflows, 6% are produced by rotation, and 8% are ambiguous. While we are unable to directly identify dual AGNs with long-slit data alone, we explore their potential kinematic classifications with this method. We also find a positive correlation between the NLR size and luminosity of the AGN NLRs (R${}_{\mathrm{NLR}}\propto \,{L}_{[{\rm{O}}\,{\rm{III}}]}^{0.21\pm 0.05}$), indicating a clumpy two-zone ionization model for the NLR.

The G7 IIIa single-lined spectroscopic binary, β Her, is studied with high-resolution, high-signal-to-noise spectra taken over 10 seasons from 23MR2000 to 10MY2009. Absolute radial velocities, corrected for convective blueshifts, are determined and new orbital parameters are derived. Line-depth ratios are used to measure temperature variation ∼2 K. A Fourier analysis is done for the line broadening, yielding a projected rotation velocity of 3.27 ± 0.20 km s⁻¹ and a radial–tangential macroturbulence dispersion of 6.43 ± 0.08 km s⁻¹. The "C" shaped bisector of Fe iλ6253 has its blue-most point at a relative flux level of 0.52, consistent with what is expected from β Her's absolute magnitude. The third-signature plot indicates granulation velocities 20% larger than the Sun's. Mapping the λ6253 line bisector onto the third-signature curve results in a flux deficit of 12.6 ± 1.0% that can be interpreted as arising from a temperature difference between granules and inter-granular lanes of 132 K. The flux deficit peaks near 5.5 km s⁻¹ from the line center, suggesting the velocity difference between granules and lanes is ∼20% larger than that found for recently analyzed K giants.

Galaxy mergers are important events that can determine the fate of a galaxy by changing its morphology, star formation activity and mass growth. Merger systems have commonly been identified from their disturbed morphologies, and we now can employ integral field spectroscopy to detect and analyze the impact of mergers on stellar kinematics as well. We visually classified galaxy morphology using deep images (${\mu }_{{\rm{r}}}=28\,\mathrm{mag}\,{\mathrm{arcsec}}^{-2}$) taken by the Blanco 4 m telescope at the Cerro Tololo Inter-American Observatory. In this paper we investigate 63 bright (${M}_{{\rm{r}}}\lt -19.3$) spectroscopically selected galaxies in Abell 119, of which 53 are early type and 20 show a disturbed morphology by visual inspection. A misalignment between the major axes in the photometric image and the kinematic map is conspicuous in morphologically disturbed galaxies. Our sample is dominated by early-type galaxies, yet it shows a surprisingly tight Tully–Fisher relation except for the morphologically disturbed galaxies which show large deviations. Three out of the eight slow rotators in our sample are morphologically disturbed. The morphologically disturbed galaxies are generally more asymmetric, visually as well as kinematically. Our findings suggest that galaxy interactions, including mergers and perhaps fly-bys, play an important role in determining the orientation and magnitude of a galaxy's angular momentum.

We present the measurement of the projected and redshift-space two-point correlation function (2pcf) of the new catalog of Chandra COSMOS-Legacy active galactic nucleus (AGN) at 2.9 ≤ z ≤ 5.5 ($\langle {L}_{\mathrm{bol}}\rangle \,\sim$ 10⁴⁶ erg s⁻¹) using the generalized clustering estimator based on phot-z probability distribution functions in addition to any available spec-z. We model the projected 2pcf, estimated using π_max = 200 h⁻¹ Mpc with the two-halo term and we derive a bias at z ∼ 3.4 equal to b = ${6.6}_{-0.55}^{+0.60}$, which corresponds to a typical mass of the hosting halos of log M_h = ${12.83}_{-0.11}^{+0.12}$ h⁻¹M_⊙. A similar bias is derived using the redshift-space 2pcf, modeled including the typical phot-z error σ_z = 0.052 of our sample at z ≥ 2.9. Once we integrate the projected 2pcf up to π_max = 200 h⁻¹ Mpc, the bias of XMM and Chandra COSMOS at z = 2.8 used in Allevato et al. is consistent with our results at higher redshifts. The results suggest only a slight increase of the bias factor of COSMOS AGNs at z ≳ 3 with the typical hosting halo mass of moderate-luminosity AGNs almost constant with redshift and equal to log M_h = ${12.92}_{-0.18}^{+0.13}$ at z = 2.8 and log M_h = ${12.83}_{-0.11}^{+0.12}$ at z ∼ 3.4, respectively. The observed redshift evolution of the bias of COSMOS AGNs implies that moderate-luminosity AGNs still inhabit group-sized halos at z ≳ 3, but slightly less massive than observed in different independent studies using X-ray AGNs at z ≤ 2.

It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudo-spectral code Dedalus. Because the flame propagation timescale is much longer than the convection timescale, we approximate the flame as fixed in space, and only consider its effects on the buoyancy of the fluid. By evolving a passive scalar field, we derive a turbulent chemical diffusivity produced by the convection as a function of height, ${D}_{{\rm{t}}}(z)$. Convection can stall a flame if the chemical mixing timescale, set by the turbulent chemical diffusivity, ${D}_{{\rm{t}}}$, is shorter than the flame propagation timescale, set by the thermal diffusivity, κ, i.e., when ${D}_{{\rm{t}}}\gt \kappa$. However, we find ${D}_{{\rm{t}}}\lt \kappa$ for most of the flame because convective plumes are not dense enough to penetrate into the flame. Extrapolating to realistic stellar conditions, this implies that convective mixing cannot stall a carbon flame and that "hybrid carbon–oxygen–neon" white dwarfs are not a typical product of stellar evolution.

The outer heliosphere is a dynamic region shaped largely by the interaction between the solar wind and the interstellar medium. While interplanetary magnetic field and plasma observations by the Voyager spacecraft have significantly improved our understanding of this vast region, modeling the outer heliosphere still remains a challenge. We simulate the three-dimensional, time-dependent solar wind flow from 1 to 80 astronomical units (au), where the solar wind is assumed to be supersonic, using a two-fluid model in which protons and interstellar neutral hydrogen atoms are treated as separate fluids. We use 1 day averages of the solar wind parameters from the OMNI data set as inner boundary conditions to reproduce time-dependent effects in a simplified manner which involves interpolation in both space and time. Our model generally agrees with Ulysses data in the inner heliosphere and Voyager data in the outer heliosphere. Ultimately, we present the model solar wind parameters extracted along the trajectory of the New Horizons spacecraft. We compare our results with in situ plasma data taken between 11 and 33 au and at the closest approach to Pluto on 2015 July 14.

Previous studies have shown that the radiation emitted by a rapidly rotating magnetar embedded in a young supernova can greatly amplify its luminosity. These one-dimensional studies have also revealed the existence of an instability arising from the piling up of radiatively accelerated matter in a thin dense shell deep inside the supernova. Here, we examine the problem in two dimensions and find that, while instabilities cause mixing and fracture this shell into filamentary structures that reduce the density contrast, the concentration of matter in a hollow shell persists. The extent of the mixing depends upon the relative energy input by the magnetar and the kinetic energy of the inner ejecta. The light curve and spectrum of the resulting supernova will be appreciably altered, as will the appearance of the supernova remnant, which will be shellular and filamentary. A similar pile up and mixing might characterize other events where energy is input over an extended period by a centrally concentrated source, e.g., a pulsar, radioactive decay, a neutrino-powered wind, or colliding shells. The relevance of our models to the recent luminous transient ASASSN-15lh is briefly discussed.

Coronal supra-arcade downflows (SADs) are observed as dark trails descending toward hot turbulent-fan-shaped regions. Due to the large temperature values and gradients in these fan regions, the thermal conduction (TC) should be very efficient. While several models have been proposed to explain the triggering and the evolution of SADs, none of these scenarios address a systematic consideration of TC. Thus, we accomplish this task numerically simulating the evolution of SADs within this framework. That is, SADs are conceived as voided (subdense) cavities formed by nonlinear waves triggered by downflowing bursty localized reconnection events in a perturbed hot fan. We generate a properly turbulent fan, obtained by a stirring force that permits control of the energy and vorticity input in the medium where SADs develop. We include anisotropic TC and consider plasma properties consistent with observations. Our aim is to study whether it is possible to prevent SADs from vanishing by thermal diffusion. We find that this will be the case, depending on the turbulence parameters, in particular if the magnetic field lines are able to envelope the voided cavities, thermally isolating them from the hot environment. Velocity shear perturbations that are able to generate instabilities of the Kelvin–Helmholtz type help to produce magnetic islands, extending the lifetime of SADs.

In a recent work by Sun et al., the color variation of quasars, namely the bluer-when-brighter trend, was found to be timescale dependent using the SDSS $g/r$ band light curves in Stripe 82. Such timescale dependence, i.e., bluer variation at shorter timescales, supports the thermal fluctuation origin of the UV/optical variation in quasars, and can be modeled well with the inhomogeneous accretion disk model. In this paper, we extend the study to much shorter wavelengths in the rest frame (down to extreme UV) using GALaxy Evolution eXplorer (GALEX) photometric data of quasars collected in two ultraviolet bands (near-UV and far-UV). We develop Monte Carlo simulations to correct for possible biases due to the considerably larger photometric uncertainties in the GALEX light curves (particularly in the far-UV, compared with the SDSS $g/r$ bands), which otherwise could produce artificial results. We securely confirm the previously discovered timescale dependence of the color variability with independent data sets and at shorter wavelengths. We further find that the slope of the correlation between the amplitude of the color variation and timescale appears even steeper than predicted by the inhomogeneous disk model, which assumes that disk fluctuations follow a damped random walk (DRW) process. The much flatter structure function observed in the far-UV compared with that at longer wavelengths implies deviation from the DRW process in the inner disk, where rest-frame extreme UV radiation is produced.

The power mechanism and accretion geometry for low-power FR 1 radio galaxies are poorly understood in comparison to those for Seyfert galaxies and QSOs. In this paper, we use the diagnostic power of the Lyα recombination line observed using the Cosmic Origins Spectrograph (COS) aboard the Hubble Space Telescope (HST) to investigate the accretion flows in three well-known, nearby FR 1s: M87, NGC 4696, and Hydra A. The Lyα emission line's luminosity, velocity structure, and the limited knowledge of its spatial extent provided by COS are used to assess conditions within a few parsecs of the supermassive black hole in these radio-mode active galactic nuclei. We observe strong Lyα emission in all three objects with total luminosity similar to that seen in BL Lacertae objects. M87 shows a complicated emission-line profile in Lyα, which varies spatially across the COS aperture and possibly temporally over several epochs of observation. In both NGC 4696 and M87, the Lyα luminosities ∼10⁴⁰ erg s⁻¹ are closely consistent with the observed strength of the ionizing continuum in Case B recombination theory and with the assumption of a near-unity covering factor. It is possible that the Lyα-emitting clouds are ionized largely by beamed radiation associated with the jets. Long-slit UV spectroscopy can be used to test this hypothesis. Hydra A and the several BL Lac objects studied in this and previous papers have Lyα luminosities larger than M87 but their extrapolated, nonthermal continua are so luminous that they overpredict the observed strength of Lyα, a clear indicator of relativistic beaming in our direction. Given their substantial space density (∼4 × 10⁻³ Mpc⁻³), the unbeamed Lyman continuum radiation of FR 1s may make a substantial minority contribution (∼10%) to the local UV background if all FR 1s are similar to M87 in ionizing flux level.

We discuss the intensity ratio of the O iv line at 1401.16 Å to the Si iv line at 1402.77 Å in Interface Region Imaging Spectrograph (IRIS) spectra. This intensity ratio is important if it can be used to measure high electron densities that cannot be measured using line intensity ratios of two different O iv lines from the multiplet within the IRIS wavelength range. Our discussion is in terms of considerably earlier observations made from the Skylab manned space station and other spectrometers on orbiting spacecraft. The earlier data on the O iv and Si iv ratio and other intersystem line ratios not available to IRIS are complementary to IRIS data. In this paper, we adopt a simple interpretation based on electron density. We adopt a set of assumptions and calculate the electron density as a function of velocity in the Si iv line profiles of two explosive events. At zero velocity the densities are about 2–3 × 10¹¹ cm⁻³, and near 200 km s⁻¹ outflow speed the densities are about 10¹² cm⁻³. The densities increase with outflow speed up to about 150 km s⁻¹ after which they level off. Because of the difference in the temperature of formation of the two lines and other possible effects such as non-ionization equilibrium, these density measurements do not have the precision that would be available if there were some additional lines near the formation temperature of O iv.

Until recently, only a handful of dusty, star-forming galaxies (DSFGs) were known at z > 4, most of them significantly amplified by gravitational lensing. Here, we have increased the number of such DSFGs substantially, selecting galaxies from the uniquely wide 250, 350, and 500 μm Herschel-ATLAS imaging survey on the basis of their extremely red far-infrared colors and faint 350 and 500 μm flux densities, based on which, they are expected to be largely unlensed, luminous, rare, and very distant. The addition of ground-based continuum photometry at longer wavelengths from the James Clerk Maxwell Telescope and the Atacama Pathfinder Experiment allows us to identify the dust peak in their spectral energy distributions (SEDs), with which we can better constrain their redshifts. We select the SED templates that are best able to determine photometric redshifts using a sample of 69 high-redshift, lensed DSFGs, then perform checks to assess the impact of the CMB on our technique, and to quantify the systematic uncertainty associated with our photometric redshifts, σ = 0.14 (1 + z), using a sample of 25 galaxies with spectroscopic redshifts, each consistent with our color selection. For Herschel-selected ultrared galaxies with typical colors of S₅₀₀/S₂₅₀ ∼ 2.2 and S₅₀₀/S₃₅₀ ∼ 1.3 and flux densities, S₅₀₀ ∼ 50 mJy, we determine a median redshift, ${\hat{z}}_{\mathrm{phot}}=3.66$, an interquartile redshift range, 3.30–4.27, with a median rest-frame 8–1000 μm luminosity, ${\hat{L}}_{\mathrm{IR}}$, of 1.3 ×  10¹³L_⊙. A third of the galaxies lie at z > 4, suggesting a space density, ρ_z > 4, of ≈6 × 10⁻⁷ Mpc⁻³. Our sample contains the most luminous known star-forming galaxies, and the most overdense cluster of starbursting proto-ellipticals found to date.

Although there has been much progress in understanding how galaxies evolve, we still do not understand how and when they stop forming stars and become quiescent. We address this by applying our galaxy spectral energy distribution models, which incorporate physically motivated star formation histories (SFHs) from cosmological simulations, to a sample of quiescent galaxies at $0.2\lt z\lt 2.1$. A total of 845 quiescent galaxies with multi-band photometry spanning rest-frame ultraviolet through near-infrared wavelengths are selected from the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) data set. We compute median SFHs of these galaxies in bins of stellar mass and redshift. At all redshifts and stellar masses, the median SFHs rise, reach a peak, and then decline to reach quiescence. At high redshift, we find that the rise and decline are fast, as expected, because the universe is young. At low redshift, the duration of these phases depends strongly on stellar mass. Low-mass galaxies ($\mathrm{log}({M}_{* }/{M}_{\odot })\sim 9.5$) grow on average slowly, take a long time to reach their peak of star formation ($\gtrsim 4$ Gyr), and then the declining phase is fast ($\lesssim 2$ Gyr). Conversely, high-mass galaxies ($\mathrm{log}({M}_{* }/{M}_{\odot })\sim 11$) grow on average fast ($\lesssim 2$ Gyr), and, after reaching their peak, decrease the star formation slowly ($\gtrsim 3$). These findings are consistent with galaxy stellar mass being a driving factor in determining how evolved galaxies are, with high-mass galaxies being the most evolved at any time (i.e., downsizing). The different durations we observe in the declining phases also suggest that low- and high-mass galaxies experience different quenching mechanisms, which operate on different timescales.

Recent observations show that white dwarfs (WDs) in cataclysmic variables (CVs) have an average mass significantly higher than isolated WDs and WDs in post-common envelope binaries (PCEBs), which are thought to be the progenitors of CVs. This suggests that either the WDs have grown in mass during the PCEB/CV evolution or the binaries with low-mass WDs are unable to evolve to be CVs. In this paper, we calculate the evolution of accreting WD binaries with updated hydrogen accumulation efficiency and angular momentum loss (AML) prescriptions. We show that thermal-timescale mass transfer is not effective in changing the average WD mass distribution. The WD mass discrepancy is most likely related to unstable mass transfer in WD binaries, in which an efficient mechanism of AML is required.

Debris disks often take the form of eccentric rings with azimuthal asymmetries in surface brightness. Such disks are often described as showing "pericenter glow," an enhancement of the disk brightness in regions nearest the central star. At long wavelengths, however, the disk apocenters should appear brighter than their pericenters: in the long-wavelength limit, we find that the apocenter/pericenter flux ratio scales as $1+e$ for disk eccentricity e. We produce new models of this "apocenter glow" to explore its causes and wavelength dependence and study its potential as a probe of dust grain properties. Based on our models, we argue that several far-infrared and (sub)millimeter images of the Fomalhaut and  Eridani debris rings obtained with Herschel, JCMT, SHARC II, ALMA, and ATCA should be reinterpreted as suggestions or examples of apocenter glow. This reinterpretation yields new constraints on the disks' dust grain properties and size distributions.

Until now, systematic errors in strong gravitational lens modeling have been acknowledged but have never been fully quantified. Here, we launch an investigation into the systematics induced by constraint selection. We model the simulated cluster Ares 362 times using random selections of image systems with and without spectroscopic redshifts and quantify the systematics using several diagnostics: image predictability, accuracy of model-predicted redshifts, enclosed mass, and magnification. We find that for models with >15 image systems, the image plane rms does not decrease significantly when more systems are added; however, the rms values quoted in the literature may be misleading as to the ability of a model to predict new multiple images. The mass is well constrained near the Einstein radius in all cases, and systematic error drops to <2% for models using >10 image systems. Magnification errors are smallest along the straight portions of the critical curve, and the value of the magnification is systematically lower near curved portions. For >15 systems, the systematic error on magnification is ∼2%. We report no trend in magnification error with the fraction of spectroscopic image systems when selecting constraints at random; however, when using the same selection of constraints, increasing this fraction up to ∼0.5 will increase model accuracy. The results suggest that the selection of constraints, rather than quantity alone, determines the accuracy of the magnification. We note that spectroscopic follow-up of at least a few image systems is crucial because models without any spectroscopic redshifts are inaccurate across all of our diagnostics.

I explore whether close-in super-Earths were formed as rocky bodies that failed to grow fast enough to become the cores of gas giants before the natal protostellar disk dispersed. I model the failed cores' inward orbital migration in the low-mass or type I regime to stopping points at distances where the tidal interaction with the protostellar disk applies zero net torque. The three kinds of migration traps considered are those due to the dead zone's outer edge, the ice line, and the transition from accretion to starlight as the disk's main heat source. As the disk disperses, the traps move toward final positions near or just outside 1 au. Planets at this location exceeding about 3 M_⊕ open a gap, decouple from their host traps, and migrate inward in the high-mass or type II regime to reach the vicinity of the star. I synthesize the population of planets that formed in this scenario, finding that a fraction of the observed super-Earths could have been failed cores. Most super-Earths that formed this way have more than 4 M_⊕, so their orbits when the disks dispersed were governed by type II migration. These planets have solid cores surrounded by gaseous envelopes. Their subsequent photoevaporative mass loss is most effective for masses originally below about 6 M_⊕. The failed core scenario suggests a division of the observed super-Earth mass–radius diagram into five zones according to the inferred formation history.

Direct imaging of exoplanets by reflected starlight is extremely challenging due to the large luminosity ratio to the primary star. Wave-front control is a critical technique to attenuate the speckle noise in order to achieve an extremely high contrast. We present a phase quantization study of a spatial light modulator (SLM) for wave-front control to meet the contrast requirement of detection of a terrestrial planet in the habitable zone of a solar-type star. We perform the numerical simulation by employing the SLM with different phase accuracy and actuator numbers, which are related to the achievable contrast. We use an optimization algorithm to solve the quantization problems that is matched to the controllable phase step of the SLM. Two optical configurations are discussed with the SLM located before and after the coronagraph focal plane mask. The simulation result has constrained the specification for SLM phase accuracy in the above two optical configurations, which gives us a phase accuracy of 0.4/1000 and 1/1000 waves to achieve a contrast of 10⁻¹⁰. Finally, we have demonstrated that an SLM with more actuators can deliver a competitive contrast performance on the order of 10⁻¹⁰ in comparison to that by using a deformable mirror.

We analyze the relationships between atomic, neutral hydrogen (H i) and star formation (SF) in the 12 low-mass SHIELD galaxies. We compare high spectral (∼0.82 km s⁻¹ ch⁻¹) and spatial resolution (physical resolutions of 160–640 pc) H i imaging from the VLA with Hα and far-ultraviolet imaging. We quantify the degree of co-spatiality between star-forming regions and regions of high H i column densities. We calculate the global star formation efficiencies (SFE; ${{\rm{\Sigma }}}_{\mathrm{SFR}}$ / ${{\rm{\Sigma }}}_{{\rm{H}}{\rm{I}}}$) and examine the relationships among the SFE and H i mass, H i column density, and star formation rate (SFR). The systems are consuming their cold neutral gas on timescales of order a few gigayears. While we derive an index for the Kennicutt–Schmidt relation of N ≈ 0.68 ± 0.04 for the SHIELD sample as a whole, the values of N vary considerably from system to system. By supplementing SHIELD results with those from other surveys, we find that H i mass and UV-based SFR are strongly correlated over five orders of magnitude. Identification of patterns within the SHIELD sample allows us to bin the galaxies into three general categories: (1) mainly co-spatial H i and SF regions, found in systems with the highest peak H i column densities and highest total H i masses; (2) moderately correlated H i and SF regions, found in systems with moderate H i column densities; and (3) obvious offsets between H i and SF peaks, found in systems with the lowest total H i masses. SF in these galaxies is dominated by stochasticity and random fluctuations in their ISM.

In this article, we compare optical light curves of two SN2002es-like Type Ia supernovae (SNe), iPTF14atg and iPTF14dpk, from the intermediate Palomar Transient Factory. Although the two light curves resemble each other around and after maximum, they show distinct early-phase rise behavior in the r-band. On the one hand, iPTF14atg revealed a slow and steady rise that lasted for 22 days with a mean rise rate of 0.2–0.3 mag day^–1, before it reached the R-band peak (−18.05 mag). On the other hand, iPTF14dpk rose rapidly to −17 mag within a day of discovery with a rise rate $\gt 1.8\,{\rm{mag}}\,{{\rm{day}}}^{-1}$, and then rose slowly to its peak (−18.19 mag) with a rise rate similar to iPTF14atg. The apparent total rise time of iPTF14dpk is therefore only 16 days. We show that emission from iPTF14atg before −17 days with respect to its maximum can be entirely attributed to radiation produced by collision between the SN and its companion star. Such emission is absent from iPTF14dpk probably because of an unfavored viewing angle, provided that SN2002es-like events arise from the same progenitor channel. We further show that an SN2002es-like SN may experience a dark phase after the explosion but before its radioactively powered light curve becomes visible. This dark phase may be lit by radiation from supernova–companion interaction.

We present new deep UBVRI images and high-resolution multi-object optical spectroscopy of the young (∼6–10 Myr old), relatively nearby (800 pc) open cluster IC 2395. We identify nearly 300 cluster members and use the photometry to estimate their spectral types, which extend from early B to middle M. We also present an infrared imaging survey of the central region using the IRAC and MIPS instruments on board the Spitzer Space Telescope, covering the wavelength range from 3.6 to 24 μm. Our infrared observations allow us to detect dust in circumstellar disks originating over a typical range of radii from ∼0.1 to ∼10 au from the central star. We identify 18 Class II, 8 transitional disk, and 23 debris disk candidates, respectively, 6.5%, 2.9%, and 8.3% of the cluster members with appropriate data. We apply the same criteria for transitional disk identification to 19 other stellar clusters and associations spanning ages from ∼1 to ∼18 Myr. We find that the number of disks in the transitional phase as a fraction of the total with strong 24 μm excesses ([8] – [24] ≥ 1.5) increases from (8.4 ± 1.3)% at ∼3 Myr to (46 ± 5)% at ∼10 Myr. Alternative definitions of transitional disks will yield different percentages but should show the same trend.

We present theoretical constraints for the formation of the newly discovered dark star clusters (DSCs) with high mass-to-light (${ \mathcal M }/{ \mathcal L }$) ratios, from Taylor et al. These compact stellar systems photometrically resemble globular clusters (GCs) but have dynamical ${ \mathcal M }/{ \mathcal L }$ ratios of ∼10–100, closer to the expectations for dwarf galaxies. The baryonic properties of the DSCs suggest that their host dark matter halos likely virialized at high redshift with ${ \mathcal M }\gt {10}^{8}{M}_{\odot }$. We use a new set of high-resolution N-body simulations of Centaurus A to determine whether there is a set of z = 0 subhalos whose properties are in line with these observations. While we find such a set of subhalos, when we extrapolate the dark matter density profiles into the inner 20 pc, no dark matter halo associated with Centaurus A in our simulations, at any redshift, can replicate the extremely high central mass densities of the DSCs. Among the most likely options for explaining 10⁵–${10}^{7}{M}_{\odot }$ within subhalos of 10 pc diameter is the presence of a central massive black hole (BH). We therefore propose that the DSCs are remnant cusps of stellar systems surrounding the central BHs of dwarf galaxies that have been almost completely destroyed by interactions with Centaurus A.

We present kinematic analyses of the 12 galaxies in the "Survey of H i in Extremely Low-mass Dwarfs" (SHIELD). We use multi-configuration interferometric observations of the H i 21 cm emission line from the Karl G. Jansky Very Large Array (VLA)²²to produce image cubes at a variety of spatial and spectral resolutions. Both two- and three-dimensional fitting techniques are employed in an attempt to derive inclination-corrected rotation curves for each galaxy. In most cases, the comparable magnitudes of velocity dispersion and projected rotation result in degeneracies that prohibit unambiguous circular velocity solutions. We thus make spatially resolved position–velocity cuts, corrected for inclination using the stellar components, to estimate the circular rotation velocities. We find ${v}_{\mathrm{circ}}$$\leqslant$ 30 km s⁻¹ for the entire survey population. Baryonic masses are calculated using single-dish H i fluxes from Arecibo and stellar masses derived from HST and Spitzer imaging. Comparison is made with total dynamical masses estimated from the position–velocity analysis. The SHIELD galaxies are then placed on the baryonic Tully–Fisher relation. There exists an empirical threshold rotational velocity, V${}_{\mathrm{rot}}$ < 15 km s⁻¹, below which current observations cannot differentiate coherent rotation from pressure support. The SHIELD galaxies are representative of an important population of galaxies whose properties cannot be described by current models of rotationally dominated galaxy dynamics.

Tidal tails are created in major mergers involving disk galaxies. It remains to be explored how the tidal tails trace the assembly history of massive galaxies. We identify a sample of 461 merging galaxies with long tidal tails, from 35,076 galaxies mass-complete at ${M}_{\star }\geqslant {10}^{9.5}\,{M}_{\odot }$ and $0.2\leqslant z\leqslant 1$, based on Hubble Space Telescope/ACS F814W imaging data and public catalogs of the COSMOS field. The long tails refer to those with length equal to or greater than the diameter of their host galaxies. The mergers with tidal tails are selected using our novel ${A}_{{\rm{O}}}-{D}_{{\rm{O}}}$ technique for strong asymmetric features, along with visual examination. Our results show that the fraction of tidal-tailed mergers evolves mildly with redshift, as $\sim {(1+z)}^{2.0\pm 0.4}$, and becomes relatively higher in less-massive galaxies, out to z = 1. With a timescale of 0.5 Gyr for the tidal-tailed mergers, we obtain that the occurrence rate of such mergers follows $0.01\pm 0.007{(1+z)}^{2.3\pm 1.4}$ Gyr⁻¹, and corresponds to ∼0.3 events since z = 1,  as well as roughly one-third of the total budget of major mergers from the literature. For disk-involved major mergers, nearly half of them have undergone a phase with long tidal tails.

Solving the mystery of the origin of chondrules is one of the most elusive goals in the field of meteoritics. Recently, the idea of planet(esimal) collisions releasing splashes of lava droplets, long considered out of favor, has been reconsidered as a possible origin of chondrules by several papers. One of the main problems with this idea is the lack of quantitative and simple models that can be used to test this scenario by directly comparing to the many known observables of chondrules. In Paper I of this series, we presented a simple thermal evolution model of a spherically symmetric expanding cloud of molten lava droplets that is assumed to emerge from a collision between two planetesimals. The production of lava could be either because the two planetesimals were already in a largely molten (or almost molten) state due to heating by ²⁶Al, or due to impact jetting at higher impact velocities. In the present paper, number II of this series, we use this model to calculate whether or not volatile elements such as Na and K will remain abundant in these droplets or whether they will get depleted due to evaporation. The high density of the droplet cloud (e.g., small distance between adjacent droplets) causes the vapor to quickly reach saturation pressure and thus shuts down further evaporation. We show to what extent, and under which conditions, this keeps the abundances of these elements high, as is seen in chondrules. We find that for most parameters of our model (cloud mass, expansion velocity, initial temperature) the volatile elements Mg, Si, and Fe remain entirely in the chondrules. The Na and K abundances inside the droplets will initially stay mostly at their initial values due to the saturation of the vapor pressure, but at some point start to drop due to the cloud expansion. However, as soon as the temperature starts to decrease, most or all of the vapor recondenses again. At the end, the Na and K elements retain most of their initial abundances, albeit occasionally somewhat reduced, depending on the parameters of the expanding cloud model. These findings appear to be qualitatively consistent with the analysis of Semarkona Type II chondrules by Hewins et al. who found evidence for sodium evaporation followed by recondensation.

The Neutron-star Interior Composition Explorer is an X-ray astrophysics payload that will be placed on the International Space Station. Its primary science goal is to measure with high accuracy the pulse profiles that arise from the non-uniform thermal surface emission of rotation-powered pulsars. Modeling general relativistic effects on the profiles will lead to measuring the radii of these neutron stars and to constraining their equation of state. Achieving this goal will depend, among other things, on accurate knowledge of the source, sky, and instrument backgrounds. We use here simple analytic estimates to quantify the level at which these backgrounds need to be known in order for the upcoming measurements to provide significant constraints on the properties of neutron stars. We show that, even in the minimal-information scenario, knowledge of the background at a few percent level for a background-to-source countrate ratio of 0.2 allows for a measurement of the neutron star compactness to better than 10% uncertainty for most of the parameter space. These constraints improve further when more realistic assumptions are made about the neutron star emission and spin, and when additional information about the source itself, such as its mass or distance, are incorporated.

I model the effect of rapid stellar rotation on a planet's insolation. Fast-rotating stars have induced pole-to-equator temperature gradients (known as gravity darkening) of up to several thousand Kelvin that affect the star's luminosity and peak emission wavelength as a function of latitude. When orbiting such a star, a planet's annual insolation can strongly vary depending on its orbital inclination. Specifically, inclined orbits result in temporary exposure to the star's hotter poles. I find that gravity darkening can drive changes in a planet's equilibrium temperature of up to ∼15% due to increased irradiance near the stellar poles. This effect can also vary a planet's exposure to UV radiation by up to ∼80% throughout its orbit as it is exposed to an irradiance spectrum corresponding to different stellar effective temperatures over time.

The key elements of the Babcock–Leighton dynamos are the generation of poloidal field through decay and the dispersal of tilted bipolar active regions and the generation of toroidal field through the observed differential rotation. These models are traditionally known as flux transport dynamo models as the equatorward propagations of the butterfly wings in these models are produced due to an equatorward flow at the bottom of the convection zone. Here we investigate the role of downward magnetic pumping near the surface using a kinematic Babcock–Leighton model. We find that the pumping causes the poloidal field to become predominately radial in the near-surface shear layer, which allows the negative radial shear to effectively act on the radial field to produce a toroidal field. We observe a clear equatorward migration of the toroidal field at low latitudes as a consequence of the dynamo wave even when there is no meridional flow in the deep convection zone. Both the dynamo wave and the flux transport type solutions are thus able to reproduce some of the observed features of the solar cycle including the 11-year periodicity. The main difference between the two types of solutions is the strength of the Babcock–Leighton source required to produce the dynamo action. A second consequence of the magnetic pumping is that it suppresses the diffusion of fields through the surface, which helps to allow an 11-year cycle at (moderately) larger values of magnetic diffusivity than have previously been used.

We present cosmological parameter constraints obtained from galaxy clusters identified by their Sunyaev–Zel'dovich effect signature in the 2500 square-degree South Pole Telescope Sunyaev Zel'dovich (SPT-SZ) survey. We consider the 377 cluster candidates identified at $z\gt 0.25$ with a detection significance greater than five, corresponding to the 95% purity threshold for the survey. We compute constraints on cosmological models using the measured cluster abundance as a function of mass and redshift. We include additional constraints from multi-wavelength observations, including Chandra X-ray data for 82 clusters and a weak lensing-based prior on the normalization of the mass-observable scaling relations. Assuming a spatially flat ΛCDM cosmology, we combine the cluster data with a prior on H₀ and find ${\sigma }_{8}=0.784\pm 0.039$ and ${{\rm{\Omega }}}_{m}=0.289\pm 0.042$, with the parameter combination ${\sigma }_{8}{({{\rm{\Omega }}}_{m}/0.27)}^{0.3}=0.797\pm 0.031$. These results are in good agreement with constraints from the cosmic microwave background (CMB) from SPT, WMAP, and Planck, as well as with constraints from other cluster data sets. We also consider several extensions to ΛCDM, including models in which the equation of state of dark energy w, the species-summed neutrino mass, and/or the effective number of relativistic species (${N}_{\mathrm{eff}}$) are free parameters. When combined with constraints from the Planck CMB, H₀, baryon acoustic oscillation, and SNe, adding the SPT cluster data improves the w constraint by 14%, to $w=-1.023\pm 0.042$.