Table of contents

We carried out 2.5-dimensional resistive MHD simulations to study the formation mechanism of molecular loops observed by Fukui et al. in the Galactic central region. Since it is hard to form molecular loops by lifting up dense molecular gas, we study the formation mechanism of molecular gas in rising magnetic arcades. This model is based on the in situ formation model of solar prominences, in which prominences are formed by cooling instability in helical magnetic flux ropes formed by imposing converging and shearing motion at footpoints of the magnetic arch anchored to the solar surface. We extended this model to Galactic center scale (a few hundreds of parsecs). Numerical results indicate that magnetic reconnection taking place in the current sheet that formed inside the rising magnetic arcade creates dense blobs confined by the rising helical magnetic flux ropes. Thermal instability taking place in the flux ropes forms dense molecular filaments floating at high Galactic latitude. The mass of the filament increases with time and can exceed ${10}^{5}$${M}_{\odot }$.

An X1.6 flare occurred in active region AR 12192 on 2014 October 22 at 14:02 UT and was observed by Hinode, IRIS, SDO, and RHESSI. We analyze a bright kernel that produces a white light (WL) flare with continuum enhancement and a hard X-ray (HXR) peak. Taking advantage of the spectroscopic observations of IRIS and Hinode/EIS, we measure the temporal variation of the plasma properties in the bright kernel in the chromosphere and corona. We find that explosive evaporation was observed when the WL emission occurred, even though the intensity enhancement in hotter lines is quite weak. The temporal correlation of the WL emission, HXR peak, and evaporation flows indicates that the WL emission was produced by accelerated electrons. To understand the WL emission process, we calculated the energy flux deposited by non-thermal electrons (observed by RHESSI) and compared it to the dissipated energy estimated from a chromospheric line (Mg ii triplet) observed by IRIS. The deposited energy flux from the non-thermal electrons is about (3–7.7) × 10¹⁰ erg cm⁻² s⁻¹ for a given low-energy cutoff of 30–40 keV, assuming the thick-target model. The energy flux estimated from the changes in temperature in the chromosphere measured using the Mg ii subordinate line is about (4.6–6.7) × 10⁹ erg cm⁻² s⁻¹: ∼6%–22% of the deposited energy. This comparison of estimated energy fluxes implies that the continuum enhancement was directly produced by the non-thermal electrons.

I demonstrate four tight correlations of total baryonic mass, velocity, and radius for a set of nearby disk galaxies: the mass–velocity relation ${M}_{{\rm{t}}}\propto {V}^{4};$ the mass–radius relation ${M}_{{\rm{t}}}\propto {R}^{2}$; the radius–velocity relation $R\propto {V}^{2};$ and the mass–radius–velocity relation ${M}_{{\rm{t}}}\propto {{RV}}^{2}$. The mass–velocity relation is the familiar Baryonic Tully–Fisher relation, and versions of the other three relations, using magnitude rather than baryonic mass, are also well known. These four observed correlations follow from a pair of more fundamental relations. First, the centripetal acceleration at the edge of the stellar disk is proportional to the acceleration predicted by Newtonian physics, and second, this acceleration is a constant that is related to Milgrom's constant. The two primary relations can be manipulated algebraically to generate the four observed correlations and allow little room for dark matter inside the radius of the stellar disk. The primary relations do not explain the velocity of the outer gaseous disks of spiral galaxies, which do not trace the Newtonian gravitational field of the observed matter.

We study the link between baryons and dark matter (DM) in 240 galaxies with spatially resolved kinematic data. Our sample spans 9 dex in stellar mass and includes all morphological types. We consider (1) 153 late-type galaxies (LTGs; spirals and irregulars) with gas rotation curves from the SPARC database, (2) 25 early-type galaxies (ETGs; ellipticals and lenticulars) with stellar and H i data from ATLAS${}^{3{\rm{D}}}$ or X-ray data from Chandra, and (3) 62 dwarf spheroidals (dSphs) with individual-star spectroscopy. We find that LTGs, ETGs, and "classical" dSphs follow the same radial acceleration relation: the observed acceleration (${g}_{\mathrm{obs}}$) correlates with that expected from the distribution of baryons (${g}_{\mathrm{bar}}$) over 4 dex. The relation coincides with the 1:1 line (no DM) at high accelerations but systematically deviates from unity below a critical scale of ∼10⁻¹⁰ m s⁻². The observed scatter is remarkably small ($\lesssim 0.13$ dex) and largely driven by observational uncertainties. The residuals do not correlate with any global or local galaxy property (e.g., baryonic mass, gas fraction, and radius). The radial acceleration relation is tantamount to a natural law: when the baryonic contribution is measured, the rotation curve follows, and vice versa. Including ultrafaint dSphs, the relation may extend by another 2 dex and possibly flatten at ${g}_{\mathrm{bar}}\lesssim {10}^{-12}$ m s⁻², but these data are significantly more uncertain. The radial acceleration relation subsumes and generalizes several well-known dynamical properties of galaxies, like the Tully–Fisher and Faber–Jackson relations, the "baryon-halo" conspiracies, and Renzo's rule.

We report on the analysis of the 10–1000 TeV large-scale sidereal anisotropy of Galactic cosmic rays (GCRs) with the data collected by the Tibet Air Shower Array from 1995 October to 2010 February. In this analysis, we improve the energy estimate and extend the decl. range down to −30°. We find that the anisotropy maps above 100 TeV are distinct from that at a multi-TeV band. The so-called tail-in and loss-cone features identified at low energies get less significant, and a new component appears at ∼100 TeV. The spatial distribution of the GCR intensity with an excess (7.2σ pre-trial, 5.2σ post-trial) and a deficit (−5.8σ pre-trial) are observed in the 300 TeV anisotropy map, in close agreement with IceCube's results at 400 TeV. Combining the Tibet results in the northern sky with IceCube's results in the southern sky, we establish a full-sky picture of the anisotropy in hundreds of TeV band. We further find that the amplitude of the first order anisotropy increases sharply above ∼100 TeV, indicating a new component of the anisotropy. All these results may shed new light on understanding the origin and propagation of GCRs.

Polarized foreground emission is a potential contaminant of attempts to measure the fluctuation power spectrum of highly redshifted 21 cm H i emission from the epoch of reionization. Using the Donald C. Backer Precision Array for Probing the Epoch of Reionization, we present limits on the observed power spectra of all four Stokes parameters in two frequency bands, centered at 126 MHz (z = 10.3) and 164 MHz (z = 7.66), for a three-month observing campaign of a deployment involving 32 antennas, for which results on the unpolarized power spectrum have been reported at z = 7.7 (by Parsons et al.) and at $7.5\lt z\lt 10.5$ (by Jacobs et al.). The power spectra in this paper are processed in the same way as by those authors, and show no definitive detection of polarized power. This nondetection is consistent with what is known about polarized sources, combined with the suppression of polarized power by fluctuations in the ionospheric rotation measure, which can strongly affect Stokes Q and U. We are able to show that the net effect of polarized leakage is a negligible contribution at the levels of the limits reported by Parsons et al. and Jacobs et al.

Long-term radio monitoring of the broad absorption line (BAL) quasar Mrk 231 at 17.6 GHz detected a strong flare in 2015. This triggered four epochs of Very Long Baseline Array (VLBA) observations from 8.4 to 43 GHz as well as three epochs of X-ray observations with NuSTAR and two with XMM over a 15 week period. Two ejected components were detected by the VLBA observations. A conservative lower bound on the apparent speed of the first ejection is attained by assuming that it was ejected when the flare began, ${v}_{\mathrm{app}}\gt 3.15c$. Serendipitous far-UV Hubble Space Telescope observations combined with our long-term radio monitoring seem to indicate that episodes of relativistic ejections suppress flux that is emitted at wavelengths shortward of the peak of the far-UV spectral energy distribution, similar to what has been observed in radio-loud quasars. Episodes of strong jet production also seem to suppress the high-ionization BAL wind seen in weak jet states. We found a statistically significant increase ($\sim 25 \%$) of the 3–12 keV flux during the radio flare relative to a quiescent radio state. This is explained by an ultra-fast (∼0.06c) X-ray-absorbing photoionized wind that is significantly detected only in the low-radio state (similar to Galactic black holes). Mrk 231 is becoming more radio loud. We found that the putative parsec-scale radio lobe doubled in brightness in nine years. Furthermore, large flares are more frequent, with three major flares occurring at ∼2 year intervals.

We use the lensing potential map from Planck CMB lensing reconstruction analysis and the "Public Cosmic Void Catalog" to measure the stacked void lensing potential. We have made an attempt to fit the HSW void profile parameters from the stacked lensing potential. In this profile, four parameters are needed to describe the shape of voids with different characteristic radii R_V. However, we have found that after reducing the background noise by subtracting the average background, there is a residue lensing power left in the data. The inclusion of the environment shifting parameter, ${\gamma }_{V}$, is necessary to get a better fit to the data with the residue lensing power. We divide the voids into two redshift bins: cmass1 ($0.45\lt z\lt 0.5$) and cmass2 ($0.5\lt z\lt 0.6$). Our best-fit parameters are $\alpha =1.989\pm 0.149$, $\beta =12.61\pm 0.56$, ${\delta }_{c}=-0.697\pm 0.025$, ${R}_{S}/{R}_{V}=1.039\pm 0.030$, ${\gamma }_{v}=(-7.034\pm 0.150)\times {10}^{-2}$ for the cmass1 sample with 123 voids and $\alpha =1.956\pm 0.165$, $\beta \,=12.91\pm 0.60$, ${\delta }_{c}=-0.673\pm 0.027$, ${R}_{S}/{R}_{V}=1.115\pm 0.032$, ${\gamma }_{v}=(-4.512\pm 0.114)\times {10}^{-2}$ for the cmass2 sample with 393 voids at 68% C.L. The addition of the environment shifting parameter is consistent with the conjecture that the Sloan Digital Sky Survey voids reside in an underdense region.

Correcting Type Ia Supernova brightnesses for extinction by dust has proven to be a vexing problem. Here we study the dust foreground to the highly reddened SN 2012cu, which is projected onto a dust lane in the galaxy NGC 4772. The analysis is based on multi-epoch, spectrophotometric observations spanning from 3300–9200 Å, obtained by the Nearby Supernova Factory. Phase-matched comparison of the spectroscopically twinned SN 2012cu and SN 2011fe across 10 epochs results in the best-fit color excess of ($E(B-V)$, RMS) = (1.00, 0.03) and total-to-selective extinction ratio of (R_V, RMS) = (2.95, 0.08) toward SN 2012cu within its host galaxy. We further identify several diffuse interstellar bands and compare the 5780 Å band with the dust-to-band ratio for the Milky Way (MW). Overall, we find the foreground dust-extinction properties for SN 2012cu to be consistent with those of the MW. Furthermore, we find no evidence for significant time variation in any of these extinction tracers. We also compare the dust extinction curve models of Cardelli et al., O'Donnell, and Fitzpatrick, and find the predictions of Fitzpatrick fit SN 2012cu the best. Finally, the distance to NGC4772, the host of SN 2012cu, at a redshift of z = 0.0035, often assigned to the Virgo Southern Extension, is determined to be 16.6 ± 1.1 Mpc. We compare this result with distance measurements in the literature.

Type Ibn supernovae (SNe) are a small yet intriguing class of explosions whose spectra are characterized by low-velocity helium emission lines with little to no evidence for hydrogen. The prevailing theory has been that these are the core-collapse explosions of very massive stars embedded in helium-rich circumstellar material (CSM). We report optical observations of six new SNe Ibn: PTF11rfh, PTF12ldy, iPTF14aki, iPTF15ul, SN 2015G, and iPTF15akq. This brings the sample size of such objects in the literature to 22. We also report new data, including a near-infrared spectrum, on the Type Ibn SN 2015U. In order to characterize the class as a whole, we analyze the photometric and spectroscopic properties of the full Type Ibn sample. We find that, despite the expectation that CSM interaction would generate a heterogeneous set of light curves, as seen in SNe IIn, most Type Ibn light curves are quite similar in shape, declining at rates around 0.1 mag day⁻¹ during the first month after maximum light, with a few significant exceptions. Early spectra of SNe Ibn come in at least two varieties, one that shows narrow P Cygni lines and another dominated by broader emission lines, both around maximum light, which may be an indication of differences in the state of the progenitor system at the time of explosion. Alternatively, the spectral diversity could arise from viewing-angle effects or merely from a lack of early spectroscopic coverage. Together, the relative light curve homogeneity and narrow spectral features suggest that the CSM consists of a spatially confined shell of helium surrounded by a less dense extended wind.

Several pulsar wind nebulae (PWN) have been detected in the TeV band in the last decade. TeV emission is typically interpreted in a purely leptonic scenario, but this often requires that the magnetic field in the nebula be much lower than the equipartition value, as well as the assumption of an enhanced density of target radiation at IR frequencies. In this work, we consider the possibility that, in addition to the relativistic electrons and positrons, relativistic hadrons are also present in these nebulae. Assuming that some of the emitted TeV photons are of hadronic origin, we compute the associated flux of $\sim 1\mbox{--}100$ TeV neutrinos. We use IceCube non-detection to put constraints on the fraction of TeV photons that might be contributed by hadrons and estimate the number of neutrino events that can be expected from these sources in ANTARES and KM3Net.

Using the high-quality observations of the Solar Dynamics Observatory, we present the interaction of two filaments (F1 and F2) in a long filament channel associated with twin coronal mass ejections (CMEs) on 2016 January 26. Before the eruption, a sequence of rapid cancellation and emergence of the magnetic flux has been observed, which likely triggered the ascending of the west filament (F1). The east footpoints of rising F1 moved toward the east far end of the filament channel, accompanied by post-eruption loops and flare ribbons. This likely indicated a large-scale eruption involving the long filament channel, which resulted from the interaction between F1 and the east filament (F2). Some bright plasma flew over F2, and F2 stayed at rest during the eruption, likely due to the confinement of its overlying lower magnetic field. Interestingly, the impulsive F1 pushed its overlying magnetic arcades to form the first CME, and F1 finally evolved into the second CME after the collision with the nearby coronal hole. We suggest that the interaction of F1 and the overlying magnetic field of F2 led to the merging reconnection that forms a longer eruptive filament loop. Our results also provide a possible picture of the origin of twin CMEs and show that the large-scale magnetic topology of the coronal hole is important for the eventual propagation direction of CMEs.

Galaxy mass assembly is an end product of structure formation in the ΛCDM cosmology. As an extension of Lee & Yi, we investigate the assembly history of stellar components in galaxies as a function of halo environments and stellar mass using semi-analytic approaches. In our fiducial model, halo mass intrinsically determines the formation and assembly of the stellar mass. Overall, the ex situ fraction slowly increases in central galaxies with increasing halo mass but sharply increases for $\mathrm{log}{M}_{* }/{M}_{\odot }\gtrsim 11$. A similar trend is also found in satellite galaxies, which implies that mergers are essential to build stellar masses above $\mathrm{log}{M}_{* }/{M}_{\odot }\sim 11$. We also examine the time evolution of the contribution of mass growth channels. Mergers become the primary channel in the mass growth of central galaxies when their host halo mass begins to exceed $\mathrm{log}{M}_{200}/{M}_{\odot }\sim 13$. However, satellite galaxies seldom reach the merger-dominant phase despite their reduced star-formation activities due to environmental effects.

We present a principal component (PC) analysis of 23 line-of-sight parameters (including the strengths of 16 diffuse interstellar bands, DIBs) for a well-chosen sample of single-cloud sightlines representing a broad range of environmental conditions. Our analysis indicates that the majority (∼93%) of the variations in the measurements can be captured by only four parameters The main driver (i.e., the first PC) is the amount of DIB-producing material in the line of sight, a quantity that is extremely well traced by the equivalent width of the λ5797 DIB. The second PC is the amount of UV radiation, which correlates well with the λ5797/λ5780 DIB strength ratio. The remaining two PCs are more difficult to interpret, but are likely related to the properties of dust in the line of sight (e.g., the gas-to-dust ratio). With our PCA results, the DIBs can then be used to estimate these line-of-sight parameters.

The Telescope Array (TA) shows a 20° hotspot as well as an excess of ultra-high-energy cosmic-rays (UHECRs) above 50 EeV when compared with the Auger spectrum. We consider the possibility that both the TA excess and hotspot are due to a dominant source in the northern sky. We carry out detailed simulations of UHECR propagation in both the intergalactic medium and the Galaxy, using different values for the intergalactic magnetic field. We consider two general classes of sources: transients and steady, adopting a mixed UHECR composition that is consistent with the one found by Auger. The spatial location of the sources is drawn randomly. We generate Auger-like and TA-like data sets from which we determine the spectrum, the sky maps, and the level of anisotropy. We find that, while steady sources are favored over transients, it is unlikely to account for all the currently available observational data. While we reproduce fairly well the Auger spectrum for the vast majority of the simulated data sets, most of the simulated data sets with a spectrum compatible with that of TA (at most a few percent depending on density model tested) show a much stronger anisotropy than the one observed. We find that the rare cases in which both the spectrum and the anisotropy are consistent require a steady source within ∼10 Mpc, to account for the flux excess, and a strong extragalactic magnetic field ∼10 nG, to reduce the excessive anisotropy.

We present a detailed study of the rest-optical (3600–7000 Å) nebular spectra of ∼380 star-forming galaxies at $z\simeq 2\mbox{--}3$, obtained with Keck/Multi-object Spectrometer for Infrared Exploration (MOSFIRE) as part of the Keck Baryonic Structure Survey (KBSS). The KBSS-MOSFIRE sample is representative of star-forming galaxies at these redshifts, with stellar masses ${M}_{* }={10}^{9}\mbox{--}{10}^{11.5}$${M}_{\odot }$ and star formation rates SFR = 3–1000 ${M}_{\odot }$ yr⁻¹. We focus on robust measurements of many strong diagnostic emission lines for individual galaxies: [O ii]λλ3727, 3729, [Ne iii]λ3869, Hβ, [O iii]$\lambda \lambda$ 4960, 5008, [N ii]λλ 6549, 6585, Hα, and [S ii]λλ6718, 6732. Comparisons with observations of typical local galaxies from the Sloan Digital Sky Survey and between subsamples of KBSS-MOSFIRE show that high-redshift galaxies exhibit a number of significant differences in addition to the well-known offset in log([O iii]$\lambda 5008$/Hβ) and log([N ii]$\lambda 6585$/Hα). We argue that the primary difference between H ii regions in $z\sim 2.3$ galaxies and those at $z\sim 0$ is an enhancement in the degree of nebular excitation, as measured by [O iii]/Hβ and ${\rm{R}}23\equiv \mathrm{log}$[([O iii]$\lambda \lambda 4960,5008$+[O ii]$\lambda \lambda 3727,3729$)/Hβ]. At the same time, KBSS-MOSFIRE galaxies are ∼10 times more massive than $z\sim 0$ galaxies with similar ionizing spectra and have higher N/O (likely accompanied by higher O/H) at fixed excitation. These results indicate the presence of harder ionizing radiation fields at fixed N/O and O/H relative to typical $z\sim 0$ galaxies, consistent with Fe-poor stellar population models that include massive binaries, and highlight a population of massive, high-specific star formation rate galaxies at high redshift with systematically different star formation histories than galaxies of similar stellar mass today.

NGC 1448 is one of the nearest luminous galaxies (L_8–1000μm > 10⁹L_⊙) to ours (z = 0.00390), and yet the active galactic nucleus (AGN) it hosts was only recently discovered, in 2009. In this paper, we present an analysis of the nuclear source across three wavebands: mid-infrared (MIR) continuum, optical, and X-rays. We observed the source with the Nuclear Spectroscopic Telescope Array (NuSTAR), and combined these data with archival Chandra data to perform broadband X-ray spectral fitting (≈0.5–40 keV) of the AGN for the first time. Our X-ray spectral analysis reveals that the AGN is buried under a Compton-thick (CT) column of obscuring gas along our line of sight, with a column density of N_H(los) ≳ 2.5 × 10²⁴ cm⁻². The best-fitting torus models measured an intrinsic 2–10 keV luminosity of L${}_{2-10,\mathrm{int}}\,=$ (3.5–7.6) × 10⁴⁰ erg s⁻¹, making NGC 1448 one of the lowest luminosity CTAGNs known. In addition to the NuSTAR observation, we also performed optical spectroscopy for the nucleus in this edge-on galaxy using the European Southern Observatory New Technology Telescope. We re-classify the optical nuclear spectrum as a Seyfert on the basis of the Baldwin–Philips–Terlevich diagnostic diagrams, thus identifying the AGN at optical wavelengths for the first time. We also present high spatial resolution MIR observations of NGC 1448 with Gemini/T-ReCS, in which a compact nucleus is clearly detected. The absorption-corrected 2–10 keV luminosity measured from our X-ray spectral analysis agrees with that predicted from the optical [O iii]λ5007 Å emission line and the MIR 12 μm  continuum, further supporting the CT nature of the AGN.

Single, double, and triple photoionization of Ne⁺ ions by single photons have been investigated at the synchrotron radiation source PETRA III in Hamburg, Germany. Absolute cross-sections were measured by employing the photon–ion merged-beams technique. Photon energies were between about 840 and 930 eV, covering the range from the lowest-energy resonances associated with the excitation of one single K-shell electron up to double excitations involving one K- and one L-shell electron, well beyond the K-shell ionization threshold. Also, photoionization of neutral Ne was investigated just below the K edge. The chosen photon energy bandwidths were between 32 and 500 meV, facilitating the determination of natural line widths. The uncertainty of the energy scale is estimated to be 0.2 eV. For comparison with existing theoretical calculations, astrophysically relevant photoabsorption cross-sections were inferred by summing the measured partial ionization channels. Discussion of the observed resonances in the different final ionization channels reveals the presence of complex Auger-decay mechanisms. The ejection of three electrons from the lowest K-shell-excited Ne⁺($1s2{s}^{2}2{p}^{6}{}^{2}{{\rm{S}}}_{1/2}$) level, for example, requires cooperative interaction of at least four electrons.

We present near-infrared spectra for 144 candidate planetary systems identified during Campaigns 1–7 of the NASA K2 Mission. The goal of the survey was to characterize planets orbiting low-mass stars, but our Infrared Telescope Facility/SpeX and Palomar/TripleSpec spectroscopic observations revealed that 49% of our targets were actually giant stars or hotter dwarfs reddened by interstellar extinction. For the 72 stars with spectra consistent with classification as cool dwarfs (spectral types K3–M4), we refined their stellar properties by applying empirical relations based on stars with interferometric radius measurements. Although our revised temperatures are generally consistent with those reported in the Ecliptic Plane Input Catalog (EPIC), our revised stellar radii are typically 0.13 ${R}_{\odot }$ (39%) larger than the EPIC values, which were based on model isochrones that have been shown to underestimate the radii of cool dwarfs. Our improved stellar characterizations will enable more efficient prioritization of K2 targets for follow-up studies.

Recent observations have shown that a growing number of the most massive Galactic globular clusters contain multiple populations of stars with different [Fe/H] and neutron-capture element abundances. NGC 6273 has only recently been recognized as a member of this "iron-complex" cluster class, and we provide here a chemical and kinematic analysis of >300 red giant branch and asymptotic giant branch member stars using high-resolution spectra obtained with the Magellan–M2FS and VLT–FLAMES instruments. Multiple lines of evidence indicate that NGC 6273 possesses an intrinsic metallicity spread that ranges from about [Fe/H] = −2 to −1 dex, and may include at least three populations with different [Fe/H] values. The three populations identified here contain separate first (Na/Al-poor) and second (Na/Al-rich) generation stars, but a Mg–Al anti-correlation may only be present in stars with [Fe/H] ≳ −1.65. The strong correlation between [La/Eu] and [Fe/H] suggests that the s-process must have dominated the heavy element enrichment at higher metallicities. A small group of stars with low [α/Fe] is identified and may have been accreted from a former surrounding field star population. The cluster's large abundance variations are coupled with a complex, extended, and multimodal blue horizontal branch (HB). The HB morphology and chemical abundances suggest that NGC 6273 may have an origin that is similar to ω Cen and M54.

Connecting in situ measured solar-wind plasma properties with typical regions on the Sun can provide an effective constraint and test to various solar wind models. We examine the statistical characteristics of the solar wind with an origin in different types of source regions. We find that the speed distribution of coronal-hole (CH) wind is bimodal with the slow wind peaking at ∼400 km s⁻¹ and the fast at ∼600 km s⁻¹. An anti-correlation between the solar wind speeds and the O⁷⁺/O⁶⁺ ion ratio remains valid in all three types of solar wind as well during the three studied solar cycle activity phases, i.e., solar maximum, decline, and minimum. The ${N}_{\mathrm{Fe}}/{N}_{{\rm{O}}}$ range and its average values all decrease with the increasing solar wind speed in different types of solar wind. The ${N}_{\mathrm{Fe}}/{N}_{{\rm{O}}}$ range (0.06–0.40, first ionization potential (FIP) bias range 1–7) for active region wind is wider than for CH wind (0.06–0.20, FIP bias range 1–3), while the minimum value of ${N}_{\mathrm{Fe}}/{N}_{{\rm{O}}}$ (∼ 0.06) does not change with the variation of speed, and it is similar for all source regions. The two-peak distribution of CH wind and the anti-correlation between the speed and O⁷⁺/O⁶⁺ in all three types of solar wind can be explained qualitatively by both the wave-turbulence-driven and reconnection-loop-opening (RLO) models, whereas the distribution features of ${N}_{\mathrm{Fe}}/{N}_{{\rm{O}}}$ in different source regions of solar wind can be explained more reasonably by the RLO models.

By exploiting two ACS/HST data sets separated by a temporal baseline of ∼7 years, we have determined the relative stellar proper motions (PMs; providing membership) and the absolute PM of the Galactic globular cluster M71. The absolute PM has been used to reconstruct the cluster orbit within a Galactic, three-component, axisymmetric potential. M71 turns out to be in a low-latitude disk-like orbit inside the Galactic disk, further supporting the scenario in which it lost a significant fraction of its initial mass. Since large differential reddening is known to affect this system, we took advantage of near-infrared, ground-based observations to re-determine the cluster center and density profile from direct star counts. The new structural parameters turn out to be significantly different from the ones quoted in the literature. In particular, M71 has a core and a half-mass radii almost 50% larger than previously thought. Finally, we estimate that the initial mass of M71 was likely one order of magnitude larger than its current value, thus helping to solve the discrepancy with the observed number of X-ray sources.

The James Webb Space Telescope's Medium Resolution Spectrometer (MRS), will offer nearly two orders of magnitude improvement in sensitivity and >3×  improvement in spectral resolution over our previous space-based mid-IR spectrometer, the Spitzer IRS. In this paper, we make predictions for spectroscopic pointed observations and serendipitous detections with the MRS. Specifically, pointed observations of Herschel sources require only a few minutes on source integration for detections of several star-forming and active galactic nucleus lines, out to z = 3 and beyond. But the same data will also include tens of serendipitous 0 ≲ z ≲ 4 galaxies per field with infrared luminosities ranging ∼10⁶–10¹³L_☉. In particular, for the first time and for free we will be able to explore the L_IR < 10⁹L_☉ regime out to z ∼ 3. We estimate that with ∼ 100 such fields, statistics of these detections will be sufficient to constrain the evolution of the low-L end of the infrared luminosity function, and hence the star formation rate function. The above conclusions hold for a wide range in the potential low-L end of the IR luminosity function, and account for the PAH deficit in low-L, low-metallicity galaxies.

The AMS-02 has measured the cosmic-ray electron (plus positron) spectrum up to ∼TeV with unprecedented precision. The spectrum can be well described by a power law without any obvious features above 10 GeV. The satellite instrument Dark Matter Particle Explorer (DAMPE), which was launched a year ago, will measure the electron spectrum up to 10 TeV with high-energy resolution. The cosmic electrons beyond TeV may be attributed to few local cosmic-ray sources, such as supernova remnants. Therefore, spectral features, such as cut-off and bumps, can be expected at high energies. In this work, we provide a careful study on the perspective of the electron spectrum beyond TeV. We first examine our astrophysical source models on the latest leptonic data of AMS-02 to give a self-consistent picture. Then we focus on the discussion about the candidate sources, which could be electron contributors above TeV. Depending on the properties of the local sources (especially on the nature of Vela), DAMPE may detect interesting features in the electron spectrum above TeV in the future.

The 2175 Å UV extinction feature was discovered in the mid-1960s, yet its physical origin remains poorly understood. One suggestion is absorption by polycyclic aromatic hydrocarbon (PAH) molecules, which is supported by theoretical molecular structure computations and by laboratory experiments. PAHs are positively detected by their 3.3, 6.2, 7.7, 8.6, 11.3, and 12.7 μm IR emission bands, which are specified by their modes of vibration. A definitive empirical link between the 2175 Å UV extinction and the IR PAH emission bands, however, is still missing. We present a new sample of hot stars that have both 2175 Å absorption and IR PAH emission. We find significant shifts of the central wavelength of the UV absorption feature, up to 2350 Å, but predominantly in stars that also have IR PAH emission. These UV shifts depend on stellar temperature in a fashion that is similar to the shifts of the 6.2 and 7.7 μm IR PAH bands, that is, the features are increasingly more redshifted as the stellar temperature decreases, but only below ∼15 kK. Above 15 kK both UV and IR features retain their nominal values. Moreover, we find a suggestive correlation between the UV and IR shifts. We hypothesize that these similar dependences of both the UV and IR features on stellar temperature hint at a common origin of the two in PAH molecules and may establish the missing link between the UV and IR observations. We further suggest that the shifts depend on molecular size, and that the critical temperature of ∼15 kK above which no shifts are observed is related to the onset of UV-driven hot-star winds and their associated shocks.

We present a sample of 1483 sources that display spectral peaks between 72 MHz and 1.4 GHz, selected from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. The GLEAM survey is the widest fractional bandwidth all-sky survey to date, ideal for identifying peaked-spectrum sources at low radio frequencies. Our peaked-spectrum sources are the low-frequency analogs of gigahertz-peaked spectrum (GPS) and compact-steep spectrum (CSS) sources, which have been hypothesized to be the precursors to massive radio galaxies. Our sample more than doubles the number of known peaked-spectrum candidates, and 95% of our sample have a newly characterized spectral peak. We highlight that some GPS sources peaking above 5 GHz have had multiple epochs of nuclear activity, and we demonstrate the possibility of identifying high-redshift (z > 2) galaxies via steep optically thin spectral indices and low observed peak frequencies. The distribution of the optically thick spectral indices of our sample is consistent with past GPS/CSS samples but with a large dispersion, suggesting that the spectral peak is a product of an inhomogeneous environment that is individualistic. We find no dependence of observed peak frequency with redshift, consistent with the peaked-spectrum sample comprising both local CSS sources and high-redshift GPS sources. The 5 GHz luminosity distribution lacks the brightest GPS and CSS sources of previous samples, implying that a convolution of source evolution and redshift influences the type of peaked-spectrum sources identified below 1 GHz. Finally, we discuss sources with optically thick spectral indices that exceed the synchrotron self-absorption limit.

Recent radio observations show that giant molecular cloud (GMC) mass functions noticeably vary across galactic disks. High-resolution magnetohydrodynamics simulations show that multiple episodes of compression are required for creating a molecular cloud in the magnetized interstellar medium. In this article, we formulate the evolution equation for the GMC mass function to reproduce the observed profiles, for which multiple compressions are driven by a network of expanding shells due to H ii regions and supernova remnants. We introduce the cloud–cloud collision (CCC) terms in the evolution equation in contrast to previous work (Inutsuka et al.). The computed time evolution suggests that the GMC mass function slope is governed by the ratio of GMC formation timescale to its dispersal timescale, and that the CCC effect is limited only in the massive end of the mass function. In addition, we identify a gas resurrection channel that allows the gas dispersed by massive stars to regenerate GMC populations or to accrete onto pre-existing GMCs. Our results show that almost all of the dispersed gas contributes to the mass growth of pre-existing GMCs in arm regions whereas less than 60% contributes in inter-arm regions. Our results also predict that GMC mass functions have a single power-law exponent in the mass range <10^5.5${M}_{\odot }$ (where ${M}_{\odot }$ represents the solar mass), which is well characterized by GMC self-growth and dispersal timescales. Measurement of the GMC mass function slope provides a powerful method to constrain those GMC timescales and the gas resurrecting factor in various environments across galactic disks.

We study the signatures of reionization and ionizing properties of early galaxies in the cross-correlations between the 21 cm emission from the spin-flip transition of neutral hydrogen (H i) and the underlying galaxy population. In particular, we focus on a sub-population of galaxies visible as Lyα Emitters (LAEs). With both observables simultaneously derived from a $z\simeq 6.6$ hydrodynamical simulation (GADGET-2) snapshot post-processed with a radiative transfer code (pCRASH) and a dust model, we perform a parameter study and aim to constrain both the average intergalactic medium (IGM) ionization state ($1-\langle {\chi }_{{\rm{H}}{\rm{I}}}\rangle$) and the reionization topology (outside-in versus inside-out). We find that, in our model, LAEs occupy the densest and most-ionized regions resulting in a very strong anti-correlation between the LAEs and the 21 cm emission. A 1000 hr Square Kilometer Array (SKA)-LOW1—Subaru Hyper Suprime-Cam experiment can provide constraints on $\langle {\chi }_{{\rm{H}}{\rm{I}}}\rangle$, allowing us to distinguish between IGM ionization levels of 50%, 25%, 10%, and fully ionized at scales $r\lesssim 10$ comoving Mpc (assuming foreground avoidance for SKA). Our results support the inside-out reionization scenario where the densest knots (under-dense voids) are ionized first (last) for $\langle {\chi }_{{\rm{H}}{\rm{I}}}\rangle \gtrsim 0.1$. Further, 1000 hr SKA-LOW1 observations should be able to confirm the inside-out scenario by detecting a lower 21 cm brightness temperature (by about 2–10 mK) in the densest regions (≳2 arcmin scales) hosting LAEs, compared to lower-density regions devoid of them.

We report new spectroscopic and photometric observations of the main-sequence, detached, eccentric, double-lined eclipsing binary V541 Cyg (P = 15.34 days, e = 0.468). Using these observations together with existing measurements, we determine the component masses and radii to better than 1% precision: ${M}_{1}={2.335}_{-0.013}^{+0.017}\,{M}_{\odot }$, ${M}_{2}={2.260}_{-0.013}^{+0.016}\,{M}_{\odot }$, ${R}_{1}={1.859}_{-0.009}^{+0.012}\,{R}_{\odot }$, and ${R}_{2}={1.808}_{-0.013}^{+0.015}\,{R}_{\odot }$. The nearly identical B9.5 stars have estimated effective temperatures of 10650 ± 200 K and 10350 ± 200 K. A comparison of these properties with current stellar evolution models shows excellent agreement at an age of about 190 Myr and [Fe/H] ≈ −0.18. Both components are found to be rotating at the pseudo-synchronous rate. The system displays a slow periastron advance that is dominated by general relativity (GR), and has previously been claimed to be slower than predicted by theory. Our new measurement, $\dot{\omega }={0.859}_{-0.017}^{+0.042}$ deg century⁻¹, has an 88% contribution from GR and agrees with the expected rate within the uncertainties. We also clarify the use of the gravity darkening coefficients in the light-curve fitting Eclipsing Binary Orbit Program (EBOP), a version of which we use here.

Flare-associated coronal shock waves sometimes interact with solar prominences, leading to large-amplitude prominence oscillations (LAPOs). Such prominence activation gives us a unique opportunity to track the time evolution of shock–cloud interaction in cosmic plasmas. Although the dynamics of interstellar shock–cloud interaction has been  extensively studied, coronal shock–solar prominence interaction is rarely studied in the context of shock–cloud interaction. Associated with the X5.4 class solar flare that occurred on 2012 March 7, a globally propagated coronal shock wave interacted with a polar prominence, leading to LAPO. In this paper, we studied bulk acceleration and excitation of the internal flow of the shocked prominence using three-dimensional magnetohydrodynamic (MHD) simulations. We studied eight MHD simulation runs, each with different mass density structure of the prominence, and one hydrodynamic simulation run, and we compared the result. In order to compare the observed motion of activated prominence with the corresponding simulation, we also studied prominence activation by injection of a triangular-shaped coronal shock. We found that the prominence is first accelerated mainly by magnetic tension force as well as direct transmission of the shock, and later decelerated mainly by magnetic tension force. The internal flow, on the other hand, is excited during the shock front sweeps through the prominence and damps almost exponentially. We construct a phenomenological model of bulk momentum transfer from the shock to the prominence, which agreed quantitatively with all the simulation results. Based on the phenomenological prominence activation model, we diagnosed physical parameters of the coronal shock wave. The estimated energy of the coronal shock is several percent of the total energy released during the X5.4 flare.

In light of recent observational results indicating an apparent lack of correlation between the anomalous microwave emission (AME) and mid-infrared emission from polycyclic aromatic hydrocarbons, we assess whether rotational emission from spinning silicate and/or iron nanoparticles could account for the observed AME without violating observational constraints on interstellar abundances, ultraviolet extinction, and infrared emission. By modifying the SpDust code to compute the rotational emission from these grains, we find that nanosilicate grains could account for the entirety of the observed AME, whereas iron grains could be responsible for only a fraction, even for extreme assumptions on the amount of interstellar iron concentrated in ultrasmall iron nanoparticles. Given the added complexity of contributions from multiple grain populations to the total spinning dust emission, as well as existing uncertainties due to the poorly constrained grain size, charge, and dipole moment distributions, we discuss generic, carrier-independent predictions of spinning dust theory and observational tests that could help identify the AME carrier(s).

We report observations of CO(J = 2 → 1) and $\mathrm{CO}(J=3\to 2)$ line emission toward the quadruply-lensed quasar RXS J1131−1231 at z = 0.654 obtained using the Plateau de Bure Interferometer (PdBI) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA). Our lens modeling shows that the asymmetry in the double-horned CO(J = 2 → 1) line profile is mainly a result of differential lensing, where the magnification factor varies from ∼3 to ∼9 across different kinematic components. The intrinsically symmetric line profile and a smooth source-plane velocity gradient suggest that the host galaxy is an extended rotating disk, with a CO size of ${R}_{\mathrm{CO}}\sim 6$ kpc and a dynamical mass of ${M}_{\mathrm{dyn}}\sim 8\times {10}^{10}$M_⊙. We also find a secondary CO-emitting source near RXS J1131−1231, the location of which is consistent with the optically-faint companion reported in previous studies. The lensing-corrected molecular gas masses are M_gas = (1.4 ± 0.3) × 10¹⁰M_⊙ and (2.0 ± 0.1) × 10⁹M_⊙ for RXS J1131−1231 and the companion, respectively. We find a lensing-corrected stellar mass of M_* = (3 ± 1) × 10¹⁰M_⊙ and a star formation rate of SFR_FIR = (120 ± 63) M_⊙ yr⁻¹, corresponding to a specific SFR and star formation efficiency comparable to z ∼ 1 disk galaxies not hosting quasars. The implied gas mass fraction of ∼18 ± 4% is consistent with the previously observed cosmic decline since z ∼ 2. We thus find no evidence for quenching of star formation in RXS J1131−1231. This agrees with our finding of an elevated ${M}_{\mathrm{BH}}/{M}_{\mathrm{bulge}}$ ratio of >0.27${}_{-0.08}^{+0.11}$% compared to the local value, suggesting that the bulk of its black hole mass is largely in place while its stellar bulge is still assembling.

In barred galaxies, the contours of stellar velocity dispersions (σ) are generally expected to be oval and aligned with the orientation of bars. However, many double-barred (S2B) galaxies exhibit distinct σ peaks on the minor axis of the inner bar, which we termed "σ-humps," while two local σ minima are present close to the ends of inner bars, i.e., "σ-hollows." Analysis of numerical simulations shows that ${\sigma }_{z}$-humps or hollows should play an important role in generating the observed σ-humps+hollows in low-inclination galaxies. In order to systematically investigate the properties of ${\sigma }_{z}$ in barred galaxies, we apply the vertical Jeans equation to a group of well-designed three-dimensional bar+disk(+bulge) models. A vertically thin bar can lower ${\sigma }_{z}$ along the bar and enhance it perpendicular to the bar, thus generating ${\sigma }_{z}$-humps+hollows. Such a result suggests that ${\sigma }_{z}$-humps+hollows can be generated by the purely dynamical response of stars in the presence of a sufficiently massive, vertically thin bar, even without an outer bar. Using self-consistent N-body simulations, we verify the existence of vertically thin bars in the nuclear-barred and S2B models that generate prominent σ-humps+hollows. Thus, the ubiquitous presence of σ-humps+hollows in S2Bs implies that inner bars are vertically thin. The addition of a bulge makes the ${\sigma }_{z}$-humps more ambiguous and thus tends to somewhat hide the ${\sigma }_{z}$-humps+hollows. We show that ${\sigma }_{z}$ may be used as a kinematic diagnostic of stellar components that have different thicknesses, providing a direct perspective on the morphology and thickness of nearly face-on bars and bulges with integral field unit spectroscopy.

We present an analysis of the mid-infrared Wide-field Infrared Survey Explorer (WISE) sources seen within the equatorial GAMA G12 field, located in the North Galactic Cap. Our motivation is to study and characterize the behavior of WISE source populations in anticipation of the deep multiwavelength surveys that will define the next decade, with the principal science goal of mapping the 3D large-scale structures and determining the global physical attributes of the host galaxies. In combination with cosmological redshifts, we identify galaxies from their WISE W1 (3.4 μm) resolved emission, and we also perform a star-galaxy separation using apparent magnitude, colors, and statistical modeling of star counts. The resulting galaxy catalog has ≃590,000 sources in 60 deg², reaching a W1 5σ depth of 31 μJy. At the faint end, where redshifts are not available, we employ a luminosity function analysis to show that approximately 27% of all WISE extragalactic sources to a limit of 17.5 mag (31 μJy) are at high redshift, $z\gt 1$. The spatial distribution is investigated using two-point correlation functions and a 3D source density characterization at 5 Mpc and 20 Mpc scales. For angular distributions, we find that brighter and more massive sources are strongly clustered relative to fainter sources with lower mass; likewise, based on WISE colors, spheroidal galaxies have the strongest clustering, while late-type disk galaxies have the lowest clustering amplitudes. In three dimensions, we find a number of distinct groupings, often bridged by filaments and superstructures. Using special visualization tools, we map these structures, exploring how clustering may play a role with stellar mass and galaxy type.

We present a combined morphological and X-ray analysis of Was 49, an isolated, dual-AGN system notable for the presence of a dominant AGN, Was 49b, in the disk of the primary galaxy, Was 49a, at a projected radial distance of 8 kpc from the nucleus. Using X-ray data from Chandra, the Nuclear Spectroscopic Telescope Array, and Swift, we find that this AGN has a bolometric luminosity of L_bol ∼ 10⁴⁵ erg s⁻¹, with a black hole mass of ${M}_{\mathrm{BH}}\,=\,{1.3}_{-0.9}^{+2.9}\times {10}^{8}\,{M}_{\odot }$. Despite the large mass, our analysis of optical data from the Discovery Channel Telescope shows that the supermassive black hole (SMBH) is hosted by a stellar counterpart with a mass of only ${5.6}_{-2.6}^{+4.9}\times {10}^{9}\,{M}_{\odot }$, which makes the SMBH potentially larger than expected from SMBH–galaxy scaling relations, and the stellar counterpart exhibits a morphology that is consistent with dwarf elliptical galaxies. Our analysis of the system in the r and K bands indicates that Was 49 is a minor merger, with the mass ratio of Was 49b to Was 49a between ∼1:7 and ∼1:15. This is in contrast with findings that the most luminous merger-triggered AGNs are found in major mergers and that minor mergers predominantly enhance AGN activity in the primary galaxy.

We use a simple organism lifecycle model to explore the viability of an atmospheric habitable zone (AHZ), with temperatures that could support Earth-centric life, which sits above an environment that does not support life. To illustrate our model, we use a cool Y dwarf atmosphere, such as WISE J085510.83–0714442.5, whose 4.5–5.2 μm spectrum shows absorption features consistent with water vapor and clouds. We allow organisms to adapt to their atmospheric environment (described by temperature, convection, and gravity) by adopting different growth strategies that maximize their chance of survival and proliferation. We assume a constant upward vertical velocity through the AHZ. We found that the organism growth strategy is most sensitive to the magnitude of the atmospheric convection. Stronger convection supports the evolution of more massive organisms. For a purely radiative environment, we find that evolved organisms have a mass that is an order of magnitude smaller than terrestrial microbes, thereby defining a dynamical constraint on the dimensions of life that an AHZ can support. Based on a previously defined statistical approach, we infer that there are of the order of 10⁹ cool Y brown dwarfs in the Milky Way, and likely a few tens of these objects are within 10 pc from Earth. Our work also has implications for exploring life in the atmospheres of temperate gas giants. Consideration of the habitable volumes in planetary atmospheres significantly increases the volume of habitable space in the galaxy.

We study the spectral energy distribution (SED) of the radio continuum (RC) emission from the Key Insight in Nearby Galaxies Emitting in Radio (KINGFISHER) sample of nearby galaxies to understand the energetics and origin of this emission. Effelsberg multi-wavelength observations at 1.4, 4.8, 8.4, and 10.5 GHz combined with archive data allow us, for the first time, to determine the mid-RC (1–10 GHz, MRC) bolometric luminosities and further present calibration relations versus the monochromatic radio luminosities. The 1–10 GHz radio SED is fitted using a Bayesian Markov Chain Monte Carlo technique leading to measurements for the nonthermal spectral index (${S}_{\nu }\sim {\nu }^{-{\alpha }_{\mathrm{nt}}}$) and the thermal fraction (${f}_{\mathrm{th}}$) with mean values of ${\alpha }_{\mathrm{nt}}=0.97\,\pm \,0.16(0.79\,\pm \,0.15$ for the total spectral index) and ${f}_{\mathrm{th}}$ = (10 ± 9)% at 1.4 GHz. The MRC luminosity changes over ∼3 orders of magnitude in the sample, $4.3\times \,{10}^{2}\,{L}_{\odot }\,\lt$ MRC $\,\lt \,3.9\times \,{10}^{5}\,{L}_{\odot }$. The thermal emission is responsible for ∼23% of the MRC on average. We also compare the extinction-corrected diagnostics of the star-formation rate (SFR) with the thermal and nonthermal radio tracers and derive the first star-formation calibration relations using the MRC radio luminosity. The nonthermal spectral index flattens with increasing SFR surface density, indicating the effect of the star-formation feedback on the cosmic-ray electron population in galaxies. Comparing the radio and IR SEDs, we find that the FIR-to-MRC ratio could decrease with SFR, due to the amplification of the magnetic fields in star-forming regions. This particularly implies a decrease in the ratio at high redshifts, where mostly luminous/star-forming galaxies are detected.

We study the time lags between the continuum emission of quasars at different wavelengths, based on more than four years of multi-band (g, r, i, z) light curves in the Pan-STARRS Medium Deep Fields. As photons from different bands emerge from different radial ranges in the accretion disk, the lags constrain the sizes of the accretion disks. We select 240 quasars with redshifts of z ≈ 1 or z ≈ 0.3 that are relatively emission-line free. The light curves are sampled from day to month timescales, which makes it possible to detect lags on the scale of the light crossing time of the accretion disks. With the code JAVELIN, we detect typical lags of several days in the rest frame between the g band and the riz bands. The detected lags are ∼2–3 times larger than the light crossing time estimated from the standard thin disk model, consistent with the recently measured lag in NGC 5548 and microlensing measurements of quasars. The lags in our sample are found to increase with increasing luminosity. Furthermore, the increase in lags going from g − r to g − i and then to g − z is slower than predicted in the thin disk model, particularly for high-luminosity quasars. The radial temperature profile in the disk must be different from what is assumed. We also find evidence that the lags decrease with increasing line ratios between ultraviolet Fe ii lines and Mg ii, which may point to changes in the accretion disk structure at higher metallicity.

Stacks of digital astronomical images are combined in order to increase image depth. The variable seeing conditions, sky background, and transparency of ground-based observations make the coaddition process nontrivial. We present image coaddition methods that maximize the signal-to-noise ratio (S/N) and  optimized for source detection and flux measurement. We show that for these purposes, the best way to combine images is to apply a matched filter to each image using its own point-spread function (PSF) and only then to sum the images with the appropriate weights. Methods that either match the filter after coaddition or perform PSF homogenization prior to coaddition will result in loss of sensitivity. We argue that our method provides an increase of between a few and 25% in the survey speed of deep ground-based imaging surveys compared with weighted coaddition techniques. We demonstrate this claim using simulated data as well as data from the Palomar Transient Factory data release 2. We present a variant of this coaddition method, which is optimal for PSF or aperture photometry. We also provide an analytic formula for calculating the S/N for PSF photometry on single or multiple observations. In the next paper in this series, we present a method for image coaddition in the limit of background-dominated noise, which is optimal for any statistical test or measurement on the constant-in-time image (e.g., source detection, shape or flux measurement, or star–galaxy separation), making the original data redundant. We provide an implementation of these algorithms in MATLAB.

Image coaddition is one of the most basic operations that astronomers perform. In Paper I, we presented the optimal ways to coadd images in order to detect faint sources and to perform flux measurements under the assumption that the noise is approximately Gaussian. Here, we build on these results and derive from first principles a coaddition technique that is optimal for any hypothesis testing and measurement (e.g., source detection, flux or shape measurements, and star/galaxy separation), in the background-noise-dominated case. This method has several important properties. The pixels of the resulting coadded image are uncorrelated. This image preserves all the information (from the original individual images) on all spatial frequencies. Any hypothesis testing or measurement that can be done on all the individual images simultaneously, can be done on the coadded image without any loss of information. The PSF of this image is typically as narrow, or narrower than the PSF of the best image in the ensemble. Moreover, this image is practically indistinguishable from a regular single image, meaning that any code that measures any property on a regular astronomical image can be applied to it unchanged. In particular, the optimal source detection statistic derived in Paper I is reproduced by matched filtering this image with its own PSF. This coaddition process, which we call proper coaddition, can be understood as the maximum signal-to-noise ratio measurement of the Fourier transform of the image, weighted in such a way that the noise in the entire Fourier domain is of equal variance. This method has important implications for multi-epoch seeing-limited deep surveys, weak lensing galaxy shape measurements, and diffraction-limited imaging via speckle observations. The last topic will be covered in depth in future papers. We provide an implementation of this algorithm in MATLAB.

Despite over 40 years of active research, the nature of the white dwarf progenitors of SNe Ia remains unclear. However, in the last decade, various progenitor scenarios have highlighted the need for detonations to be the primary mechanism by which these white dwarfs are consumed, but it is unclear how these detonations are triggered. In this paper we study how detonations are spontaneously initiated due to temperature inhomogeneities, e.g., hotspots, in burning nuclear fuel in a simplified physical scenario. Following the earlier work by Zel'Dovich, we describe the physics of detonation initiation in terms of the comparison between the spontaneous wave speed and the Chapman–Jouguet speed. We develop an analytic expression for the spontaneous wave speed and utilize it to determine a semi-analytic criterion for the minimum size of a hotspot with a linear temperature gradient between a peak and base temperature for which detonations in burning carbon–oxygen material can occur. Our results suggest that spontaneous detonations may easily form under a diverse range of conditions, likely allowing a number of progenitor scenarios to initiate detonations that burn up the star.

We discuss the possibility that gravitational focusing is responsible for the power-law mass function of star clusters $N(\mathrm{log}M)\propto {M}^{-1}$. This power law can be produced asymptotically when the mass accretion rate of an object depends upon the mass of the accreting body, as $\dot{M}\propto {M}^{2}$. Although Bondi–Hoyle–Lyttleton accretion formally produces this dependence on mass in a uniform medium, realistic environments are much more complicated. However, numerical simulations in SPH that allow for sink formation yield such an asymptotic power-law mass function. We perform pure N-body simulations to isolate the effects of gravity from those of gas physics and to show that clusters naturally result with the power-law mass distribution. We also consider the physical conditions necessary to produce clusters on appropriate timescales. Our results help support the idea that gravitationally dominated accretion is the most likely mechanism for producing the cluster mass function.

We describe the structural, stellar population and gas properties of the nearest ultra diffuse galaxy discovered so far: UGC 2162 (z = 0.00392; ${R}_{e,g}=1.7(\pm 0.2)$ kpc; ${\mu }_{g}(0)$ = 24.4 ± 0.1 mag arcsec⁻²; $g-i$ = 0.33 ± 0.02). This galaxy, located at a distance of 12.3(±1.7) Mpc, is a member of the M77 group. UGC 2162 has a stellar mass of $\sim 2{(}_{-1}^{+2})$ × 10⁷${M}_{\odot }$ and is embedded within a cloud of HI gas ∼10 times more massive: ∼1.9(±0.6) × 10⁸${M}_{\odot }$. Using the width of its HI line as a dynamical proxy, the enclosed mass within the inner R ∼ 5 kpc is ∼4.6(±0.8) × 10⁹${M}_{\odot }$ (i.e., M/L ∼ 200). The estimated virial mass from the cumulative mass curve is ∼8(±2)×10¹⁰M_⊙. Ultra-deep imaging from the IAC Stripe82 Legacy Project show that the galaxy is irregular and has many star-forming knots, with a gas-phase metallicity around one-third of the solar value. Its estimated star-formation rate is ∼0.01 ${M}_{\odot }$ yr⁻¹. This SFR would double the stellar mass of the object in ∼2 Gyr. If the object were to stop forming stars at this moment, after a passive evolution, its surface brightness would become extremely faint: ${\mu }_{g}(0)\,\sim$ 27 mag arcsec⁻² and its size would remain large ${R}_{e,g}\,\sim$ 1.8 kpc. Such faintness would make it almost undetectable to most present-day surveys. This suggests that there could be an important population of ${M}_{\star }\,\sim$ 10⁷${M}_{\odot }$ "dark galaxies" in rich environments (depleted of HI gas) waiting to be discovered by current and future ultra-deep surveys.

We seek to characterize how the change of global rotation rate influences the overall dynamics and large-scale flows arising in the convective envelopes of stars covering stellar spectral types from early G to late K. We do so through numerical simulations with the ASH code, where we consider stellar convective envelopes coupled to a radiative interior with various global properties. As solar-like stars spin down over the course of their main sequence evolution, such a change must have a direct impact on their dynamics and rotation state. We indeed find that three main states of rotation may exist for a given star: anti-solar-like (fast poles, slow equator), solar-like (fast equator, slow poles), or a cylindrical rotation profile. Under increasingly strict rotational constraints, the last profile can further evolve into a Jupiter-like profile, with alternating prograde and retrograde zonal jets. We have further assessed how far the convection and meridional flows overshoot into the radiative zone and investigated the morphology of the established tachocline. Using simple mixing length arguments, we are able to construct a scaling of the fluid Rossby number ${R}_{{of}}=\tilde{\omega }/2{{\rm{\Omega }}}_{* }\sim \tilde{v}/2{{\rm{\Omega }}}_{* }{R}_{* }$, which we calibrate based on our 3D ASH simulations. We can use this scaling to map the behavior of differential rotation versus the global parameters of stellar mass and rotation rate. Finally, we isolate a region on this map (R_of ≳ 1.5–2) where we posit that stars with an anti-solar differential rotation may exist in order to encourage observers to hunt for such targets.

The structure of the steady axisymmetric force-free magnetosphere of a Kerr black hole (BH) is governed by a second-order partial differential equation of A_ϕ depending on two "free" functions ${\rm{\Omega }}({A}_{\phi })$ and $I({A}_{\phi })$, where A_ϕ is the ϕ component of the vector potential of the electromagnetic field, Ω is the angular velocity of the magnetic field lines, and I is the poloidal electric current. In this paper, we investigate the solution uniqueness. Taking the asymptotically uniform field as an example, analytic studies imply that there are infinitely many solutions approaching the uniform field at infinity, while only a unique one is found in general relativistic magnetohydrodynamic simulations. To settle the disagreement, we reinvestigate the structure of the governing equation and numerically solve it with given constraint and boundary conditions. We find that the constraint condition (field lines smoothly crossing the light surface) and boundary conditions at the horizon and at infinity are connected via radiation conditions at horizon and at infinity, rather than being independent. With appropriate constraint and boundary conditions, we numerically solve the governing equation and find a unique solution. Contrary to naive expectations, our numerical solution yields a discontinuity in the angular velocity of the field lines and a current sheet along the last field line crossing the event horizon. We also briefly discuss the applicability of the perturbation approach to solving the governing equation.

We present 2.5-square-degree C₂H N = 1–0 and N₂H⁺J = 1–0 maps of the ρ Ophiuchi molecular cloud complex. These are the first large-scale maps of the ρ Ophiuchi molecular cloud complex with these two tracers. The C₂H emission is spatially more extended than the N₂H⁺ emission. One faint N₂H⁺ clump, Oph-M, and one C₂H ring, Oph-RingSW, are identified for the first time. The observed C₂H-to-N₂H⁺ abundance ratio ([C₂H]/[N₂H⁺]) varies between 5 and 110. We modeled the C₂H and N₂H⁺ abundances with 1D chemical models, which show a clear decline of [C₂H]/[N₂H⁺] with chemical age. Such an evolutionary trend is little affected by temperatures when they are below 40 K. At high density (n_H > 10⁵ cm⁻³), however, the time it takes for the abundance ratio to drop at least one order of magnitude becomes less than the dynamical time (e.g., turbulence crossing time of ∼10⁵ yr). The observed [C₂H]/[N₂H⁺] difference between L1688 and L1689 can be explained by L1688 having chemically younger gas in relatively less dense regions. The observed [C₂H]/[N₂H⁺] values are the results of time evolution, accelerated at higher densities. For the relatively low density regions in L1688 where only C₂H emission was detected, the gas should be chemically younger.

The processing of the hydrocarbon ice, ethylene (C₂H₄/C₂D₄), via energetic electrons, thus simulating the processes in the track of galactic cosmic-ray particles, was carried out in an ultrahigh vacuum apparatus. The chemical evolution of the ices was monitored online and in situ utilizing Fourier transform infrared spectroscopy (FTIR) and during temperature programmed desorption, via a quadrupole mass spectrometer utilizing electron impact ionization (EI-QMS) and a reflectron time-of-flight mass spectrometer utilizing a photoionization source (PI-ReTOF-MS). Several previous in situ studies of ethylene ice irradiation using FTIR were substantiated with the detection of six products: [CH₄ (CD₄)], acetylene [C₂H₂ (C₂D₂)], the ethyl radical [C₂H₅ (C₂D₅)], ethane [C₂H₆ (C₂D₆)], 1-butene [C₄H₈ (C₄D₈)], and n-butane [C₄H₁₀ (C₄D₁₀)]. Contrary to previous gas phase studies, the PI-ReTOF-MS detected several groups of hydrocarbon with varying degrees of saturation: C_nH_2n+2 (n = 4–10), C_nH_2n (n = 2–12, 14, 16), C_nH_2n−2 (n = 3–12, 14, 16), C_nH_2n−4 (n = 4–12, 14, 16), C_nH_2n−6 (n = 4–10, 12), C_nH_2n−8 (n = 6–10), and C_nH_2n−10 (n = 6–10). Multiple laboratory studies have shown the facile production of ethylene from methane, which is a known ice constituent in the interstellar medium. Various astrophysically interesting molecules can be associated with the groups detected here, such as allene/methylacetylene (C₃H₄) or 1, 3-butadiene (C₄H₆) and its isomers, which have been shown to lead to polycyclic aromatic hydrocarbons. Finally, several hydrocarbon groups detected here are unique to ethylene ice versus ethane ice and may provide understanding of how complex hydrocarbons form in astrophysical environments.

The high-energy non-thermal universe is dominated by power-law-like spectra. Therefore, results in high-energy astronomy are often reported as parameters of power-law fits, or, in the case of a non-detection, as an upper limit assuming the underlying unseen spectrum behaves as a power law. In this paper, I demonstrate a simple and powerful one-to-one relation of the integral upper limit in the two-dimensional power-law parameter space into the spectrum parameter space and use this method to unravel the so-far convoluted question of the sensitivity of astroparticle telescopes.

In order to study chromospheric magnetosonic wave propagation including, for the first time, the effects of ion–neutral interactions in the partially ionized solar chromosphere, we have developed a new multi-fluid computational model accounting for ionization and recombination reactions in gravitationally stratified magnetized collisional media. The two-fluid model used in our 2D numerical simulations treats neutrals as a separate fluid and considers charged species (electrons and ions) within the resistive MHD approach with Coulomb collisions and anisotropic heat flux determined by Braginskiis transport coefficients. The electromagnetic fields are evolved according to the full Maxwell equations and the solenoidality of the magnetic field is enforced with a hyperbolic divergence-cleaning scheme. The initial density and temperature profiles are similar to VAL III chromospheric model in which dynamical, thermal, and chemical equilibrium are considered to ensure comparison to existing MHD models and avoid artificial numerical heating. In this initial setup we include simple homogeneous flux tube magnetic field configuration and an external photospheric velocity driver to simulate the propagation of MHD waves in the partially ionized reactive chromosphere. In particular, we investigate the loss of chemical equilibrium and the plasma heating related to the steepening of fast magnetosonic wave fronts in the gravitationally stratified medium.

We present 5–20 μm spectral maps of the reflection nebula NGC 2023 obtained with the Infrared Spectrograph SL and SH modes on board the Spitzer Space Telescope, which reveal emission from polycyclic aromatic hydrocarbons (PAHs), C₆₀, and H₂ superposed on a dust continuum. We show that several PAH emission bands correlate with each other and exhibit distinct spatial distributions that reveal a spatial sequence with distance from the illuminating star. We explore the distinct morphology of the 6.2, 7.7, and 8.6 μm PAH bands and find that at least two spatially distinct components contribute to the 7–9 μm PAH emission in NGC 2023. We report that the PAH features behave independently of the underlying plateaus. We present spectra of compact, oval PAHs ranging in size from C₆₆ to C₂₁₀, determined computationally using density functional theory, and we investigate trends in the band positions and relative intensities as a function of PAH size, charge, and geometry. Based on the NASA Ames PAH database, we discuss the 7–9 μm components in terms of band assignments and relative intensities. We assign the plateau emission to very small grains with possible contributions from PAH clusters and identify components in the 7–9 μm emission that likely originate in these structures. Based on the assignments and the observed spatial sequence, we discuss the photochemical evolution of the interstellar PAH family as the PAHs are more and more exposed to the radiation field of the central star in the evaporative flows associated with the Photo-Dissociation Regions in NGC 2023.

The importance of the magnetic (B) field in the formation of infrared dark clouds (IRDCs) and massive stars is an ongoing topic of investigation. We studied the plane-of-sky B field for one IRDC, G028.23-00.19, to understand the interaction between the field and the cloud. We used near-IR background starlight polarimetry to probe the B field and performed several observational tests to assess the field importance. The polarimetric data, taken with the Mimir instrument, consisted of H-band and K-band observations, totaling 17,160 stellar measurements. We traced the plane-of-sky B-field morphology with respect to the sky-projected cloud elongation. We also found the relationship between the estimated B-field strength and gas volume density, and we computed estimates of the normalized mass-to-magnetic flux ratio. The B-field orientation with respect to the cloud did not show a preferred alignment, but it did exhibit a large-scale pattern. The plane-of-sky B-field strengths ranged from 10 to 165 μG, and the B-field strength dependence on density followed a power law with an index consistent with 2/3. The mass-to-magnetic flux ratio also increased as a function of density. The relative orientations and relationship between the B field and density imply that the B field was not dynamically important in the formation of the IRDC. The increase in mass-to-flux ratio as a function of density, though, indicates a dynamically important B field. Therefore, it is unclear whether the B field influenced the formation of G28.23. However, it is likely that the presence of the IRDC changed the local B-field morphology.

Ages and masses of young stars are often estimated by comparing their luminosities and effective temperatures to pre-main-sequence stellar evolution tracks, but magnetic fields and starspots complicate both the observations and evolution. To understand their influence, we study the heavily spotted weak-lined T-Tauri star LkCa 4 by searching for spectral signatures of radiation originating from the starspot or starspot groups. We introduce a new methodology for constraining both the starspot filling factor and the spot temperature by fitting two-temperature stellar atmosphere models constructed from Phoenix synthetic spectra to a high-resolution near-IR IGRINS spectrum. Clearly discernable spectral features arise from both a hot photospheric component ${T}_{\mathrm{hot}}$ ∼ 4100 K and a cool component ${T}_{\mathrm{cool}}$ ∼ 2700–3000 K, which covers ∼80% of the visible surface. This mix of hot and cool emission is supported by analyses of the spectral energy distribution, rotational modulation of colors and of TiO band strengths, and features in low-resolution optical/near-IR spectroscopy. Although the revised effective temperature and luminosity make LkCa 4 appear to be much younger and of much lower mass than previous estimates from unspotted stellar evolution models, appropriate estimates will require the production and adoption of spotted evolutionary models. Biases from starspots likely afflict most fully convective young stars and contribute to uncertainties in ages and age spreads of open clusters. In some spectral regions, starspots act as a featureless "veiling" continuum owing to high rotational broadening and heavy line blanketing in cool star spectra. Some evidence is also found for an anticorrelation between the velocities of the warm and cool components.

We report ALMA Cycle 2 observations of 230 GHz (1.3 mm) dust continuum emission, and ¹²CO, ¹³CO, and C¹⁸O J = 2–1 line emission, from the Upper Scorpius transitional disk [PZ99] J160421.7-213028, with an angular resolution of ∼$0\buildrel{\prime\prime}\over{.} 25$ (35 au). Armed with these data and existing H-band scattered light observations, we measure the size and depth of the disk's central cavity, and the sharpness of its outer edge, in three components: sub-μm-sized "small" dust traced by scattered light, millimeter-sized "big" dust traced by the millimeter continuum, and gas traced by line emission. Both dust populations feature a cavity of radius ∼70 au that is depleted by factors of at least 1000 relative to the dust density just outside. The millimeter continuum data are well explained by a cavity with a sharp edge. Scattered light observations can be fitted with a cavity in small dust that has either a sharp edge at 60 au, or an edge that transitions smoothly over an annular width of 10 au near 60 au. In gas, the data are consistent with a cavity that is smaller, about 15 au in radius, and whose surface density at 15 au is ${10}^{3\pm 1}$ times smaller than the surface density at 70 au; the gas density grades smoothly between these two radii. The CO isotopologue observations rule out a sharp drop in gas surface density at 30 au or a double-drop model, as found by previous modeling. Future observations are needed to assess the nature of these gas and dust cavities (e.g., whether they are opened by multiple as-yet-unseen planets or photoevaporation).

We conducted a large spectroscopic survey of 336 red giants in the direction of the Leo II dwarf galaxy using Hectochelle on the Multiple Mirror Telescope, and we conclude that 175 of them are members based on their radial velocities and surface gravities. Of this set, 40 stars have never before been observed spectroscopically. The systemic velocity of the dwarf is 78.3 ± 0.6 km s⁻¹ with a velocity dispersion of 7.4 ± 0.4 km s⁻¹. We identify one star beyond the tidal radius of Leo II but find no signatures of uniform rotation, kinematic asymmetries, or streams. The stars show a strong metallicity gradient of −1.53 ± 0.10 dex kpc⁻¹ and have a mean metallicity of −1.70 ± 0.02 dex. There is also evidence of two different chemodynamic populations, but the signal is weak. A larger sample of stars would be necessary to verify this feature.

We consider the hydrodynamics of the outer core of a neutron star under conditions when both neutrons and protons are superfluid. Starting from the equation of motion for the phases of the wave functions of the condensates of neutron pairs and proton pairs, we derive the generalization of the Euler equation for a one-component fluid. These equations are supplemented by the conditions for conservation of neutron number and proton number. Of particular interest is the effect of entrainment, the fact that the current of one nucleon species depends on the momenta per nucleon of both condensates. We find that the nonlinear terms in the Euler-like equation contain contributions that have not always been taken into account in previous applications of superfluid hydrodynamics. We apply the formalism to determine the frequency of oscillations about a state with stationary condensates and states with a spatially uniform counterflow of neutrons and protons. The velocities of the coupled sound-like modes of neutrons and protons are calculated from properties of uniform neutron star matter evaluated on the basis of chiral effective field theory. We also derive the condition for the two-stream instability to occur.

Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution-dependent scalings. We present a subresolution model representing the three major phases of supernova blast wave evolution—free expansion, energy-conserving Sedov–Taylor, and momentum-conserving snowplow—with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, ${f}_{\mathrm{hot}}$, as ${(1-{f}_{\mathrm{hot}})}^{-4/5}$. We also include winds from young massive stars and AGB stars, Strömgren sphere gas heating by massive stars, and a mechanism that limits gas cooling that is driven by radiative recombination of dense H ii regions. We present initial tests for isolated Milky Way-like systems simulated with the Gadget-based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass-loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a variation of two orders of magnitude in particle mass in the range (1.3–130) × 10⁴ solar masses.

B2 1215+30 is a BL-Lac-type blazar that was first detected at TeV energies by the MAGIC atmospheric Cherenkov telescopes and subsequently confirmed by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) observatory with data collected between 2009 and 2012. In 2014 February 08, VERITAS detected a large-amplitude flare from B2 1215+30 during routine monitoring observations of the blazar 1ES 1218+304, located in the same field of view. The TeV flux reached 2.4 times the Crab Nebula flux with a variability timescale of $\lt 3.6\,\mathrm{hr}$. Multiwavelength observations with Fermi-LAT, Swift, and the Tuorla Observatory revealed a correlated high GeV flux state and no significant optical counterpart to the flare, with a spectral energy distribution where the gamma-ray luminosity exceeds the synchrotron luminosity. When interpreted in the framework of a one-zone leptonic model, the observed emission implies a high degree of beaming, with Doppler factor $\delta \gt 10$, and an electron population with spectral index $p\lt 2.3$.

We present an X-ray photometric analysis of six gravitationally lensed quasars, with observation campaigns spanning from 5 to 14 years, measuring the total (0.83–21.8 keV restframe), soft- (0.83–3.6 keV), and hard- (3.6–21.8 keV) band image flux ratios for each epoch. Using the ratios of the model-predicted macro-magnifications as baselines, we build differential microlensing light curves and obtain joint likelihood functions for the average X-ray emission region sizes. Our analysis yields a probability distribution function for the average half-light radius of the X-ray emission region in the sample that peaks slightly above 1 gravitational radius and with nearly indistinguishable $68 \%$ confidence (one-sided) upper limits of 17.8 and 18.9 gravitational radii for the soft and hard X-ray emitting regions, assuming a mean stellar mass of 0.3 M_⊙. We see hints of energy dependent microlensing between the soft and hard bands in two of the objects. In a separate analysis on the root-mean-square (rms) of the microlensing variability, we find significant differences between the soft and hard bands, but the sign of the difference is not consistent across the sample. This suggests the existence of some kind of spatial structure to the X-ray emission in an otherwise extremely compact source. We also discover a correlation between the rms microlensing variability and the average microlensing amplitude.

The distribution of heavy elements is anomalously low in the asteroid main belt region compared with elsewhere in the solar system. Observational surveys also indicate a deficit in the number of small (≲50 km size) asteroids, which is two orders of magnitude lower than what is expected from the single power-law distribution that results from a collisional coagulation and fragmentation equilibrium. Here, we consider the possibility that a major fraction of the original asteroid population may have been cleared out by Jupiter's secular resonance, as it swept through the main asteroid belt during the depletion of the solar nebula. This effect leads to the excitation of the asteroids' orbital eccentricities. Concurrently, hydrodynamic drag and planet–disk tidal interaction effectively damp the eccentricities of sub-100 km-size and of super-lunar-size planetesimals, respectively. These combined effects lead to the asteroids' orbital decay and clearing from the present-day main belt region (∼2.1–3.3 au). Eccentricity damping for the intermediate-size (50 to several hundreds of kilometers) planetesimals is less efficient than for small or large planetesimals. These objects therefore preferentially remain as main belt asteroids near their birthplaces, with modest asymptotic eccentricities. The smaller asteroids are the fragments of subsequent disruptive collisions at later times as suggested by the present-day asteroid families. This scenario provides a natural explanation for both the observed low surface density and the size distribution of asteroids in the main belt, without the need to invoke special planetesimal formation mechanisms. It also offers an explanation for the confined spatial extent of the terrestrial planet building blocks without the requirement of extensive migration of Jupiter, which is required in the grand-tack scenario.

The Fermi Large Area Telescope (LAT) has opened the way for comparative studies of cosmic rays (CRs) and high-energy objects in the Milky Way (MW) and in other, external, star-forming galaxies. Using 2 yr of observations with the Fermi LAT, Local Group galaxy M31 was detected as a marginally extended gamma-ray source, while only an upper limit has been derived for the other nearby galaxy M33. We revisited the gamma-ray emission in the direction of M31 and M33 using more than 7 yr of LAT Pass 8 data in the energy range $0.1\mbox{--}100\,\mathrm{GeV}$, presenting detailed morphological and spectral analyses. M33 remains undetected, and we computed an upper limit of $2.0\times {10}^{-12}\,\mathrm{erg}\,{\mathrm{cm}}^{-2}\,{{\rm{s}}}^{-1}\,$ on the $0.1\mbox{--}100\,\mathrm{GeV}$ energy flux (95% confidence level). This revised upper limit remains consistent with the observed correlation between gamma-ray luminosity and star formation rate tracers and implies an average CR density in M33 that is at most half of that of the MW. M31 is detected with a significance of nearly $10\sigma$. Its spectrum is consistent with a power law with photon index ${\rm{\Gamma }}=2.4\pm {0.1}_{\mathrm{stat}+\mathrm{syst}}$ and a $0.1\mbox{--}100\,\mathrm{GeV}$ energy flux of $(5.6\pm {0.6}_{\mathrm{stat}+\mathrm{syst}})\times {10}^{-12}\,\mathrm{erg}\,{\mathrm{cm}}^{-2}\,{{\rm{s}}}^{-1}$. M31 is detected to be extended with a $4\sigma$ significance. The spatial distribution of the emission is consistent with a uniform-brightness disk with a radius of 04 and no offset from the center of the galaxy, but nonuniform intensity distributions cannot be excluded. The flux from M31 appears confined to the inner regions of the galaxy and does not fill the disk of the galaxy or extend far from it. The gamma-ray signal is not correlated with regions rich in gas or star formation activity, which suggests that the emission is not interstellar in origin, unless the energetic particles radiating in gamma rays do not originate in recent star formation. Alternative and nonexclusive interpretations are that the emission results from a population of millisecond pulsars dispersed in the bulge and disk of M31 by disrupted globular clusters or from the decay or annihilation of dark matter particles, similar to what has been proposed to account for the so-called Galactic center excess found in Fermi-LAT observations of the MW.

We present the discovery of two extended ∼0.12 mag dimming events of the weak-lined T Tauri star V1334. The start of the first event was missed but came to an end in late 2003, and the second began in 2009 February, and continues as of 2016 November. Since the egress of the current event has not yet been observed, it suggests a period of >13 years if this event is periodic. Spectroscopic observations suggest the presence of a small inner disk, although the spectral energy distribution shows no infrared excess. We explore the possibility that the dimming events are caused by an orbiting body (e.g., a disk warp or dust trap), enhanced disk winds, hydrodynamical fluctuations of the inner disk, or a significant increase in the magnetic field flux at the surface of the star. We also find a ∼0.32 day periodic photometric signal that persists throughout the 2009 dimming which appears to not be due to ellipsoidal variations from a close stellar companion. High-precision photometric observations of V1334 Tau during K2 campaign 13, combined with simultaneous photometric and spectroscopic observations from the ground, will provide crucial information about the photometric variability and its origin.

We search for high-redshift dropout galaxies behind the Hubble Frontier Fields (HFF) galaxy cluster MACS J1149.5+2223, a powerful cosmic lens that has revealed a number of unique objects in its field. Using the deep images from the Hubble and Spitzer space telescopes, we find 11 galaxies at z > 7 in the MACS J1149.5+2223 cluster field, and 11 in its parallel field. The high-redshift nature of the bright z ≃ 9.6 galaxy MACS1149-JD, previously reported by Zheng et al., is further supported by non-detection in the extremely deep optical images from the HFF campaign. With the new photometry, the best photometric redshift solution for MACS1149-JD reduces slightly to z = 9.44 ± 0.12. The young galaxy has an estimated stellar mass of $(7\pm 2)\times {10}^{8}\,{M}_{\odot }$, and was formed at $z={13.2}_{-1.6}^{+1.9}$ when the universe was ≈300 Myr old. Data available for the first four HFF clusters have already enabled us to find faint galaxies to an intrinsic magnitude of ${M}_{{UV}}\simeq -15.5$, approximately a factor of 10 deeper than the parallel fields.

We present CO observations toward three large supernova remnants (SNRs) in the third Galactic quadrant using the Purple Mountain Observatory Delingha 13.7 m millimeter-wavelength telescope. The observations are part of the high-resolution CO survey of the Galactic plane between Galactic longitudes $l=-10^\circ$ to $250^\circ$ and latitudes $b=-5^\circ$ to $5^\circ$. CO emission was detected toward the three SNRs: G205.5+0.5 (Monoceros Nebula), G206.9+2.3 (PKS 0646+06), and G213.0–0.6. Both SNRs G205.5+0.5 and G213.0–0.6 exhibit the morphological agreement (or spatial correspondences) between the remnant and the surrounding molecular clouds (MCs), as well as kinematic signatures of shock perturbation in the molecular gas. We confirm that the two SNRs are physically associated with their ambient MCs and the shock of SNRs is interacting with the dense, clumpy molecular gas. SNR G206.9+2.3, which is close to the northeastern edge of the Monoceros Nebula, displays the spatial coincidence with molecular partial shell structures at V_LSR ∼ 15 km s⁻¹. While no significant line broadening has been detected within or near the remnant, the strong morphological correspondence between the SNR and the molecular cavity implies that SNR G206.9+2.3 is probably associated with the CO gas and is evolving in the low-density environment. The physical features of individual SNRs, together with the relationship between SNRs and their nearby objects, are also discussed.

A newly developed laboratory, Meteoric Ablation Simulator (MASI), is used to test model predictions of the atmospheric ablation of interplanetary dust particles (IDPs) with experimental Na, Fe, and Ca vaporization profiles. MASI is the first laboratory setup capable of performing time-resolved atmospheric ablation simulations, by means of precision resistive heating and atomic laser-induced fluorescence detection. Experiments using meteoritic IDP analogues show that at least three mineral phases (Na-rich plagioclase, metal sulfide, and Mg-rich silicate) are required to explain the observed appearance temperatures of the vaporized elements. Low melting temperatures of Na-rich plagioclase and metal sulfide, compared to silicate grains, preclude equilibration of all the elemental constituents in a single melt. The phase-change process of distinct mineral components determines the way in which Na and Fe evaporate. Ca evaporation is dependent on particle size and on the initial composition of the molten silicate. Measured vaporized fractions of Na, Fe, and Ca as a function of particle size and speed confirm differential ablation (i.e., the most volatile elements such as Na ablate first, followed by the main constituents Fe, Mg, and Si, and finally the most refractory elements such as Ca). The Chemical Ablation Model (CABMOD) provides a reasonable approximation to this effect based on chemical fractionation of a molten silicate in thermodynamic equilibrium, even though the compositional and geometric description of IDPs is simplistic. Improvements in the model are required in order to better reproduce the specific shape of the elemental ablation profiles.

The formation mechanism of CO clouds observed with the NANTEN2 and Mopra telescopes toward the stellar cluster Westerlund 2 is studied by 3D magnetohydrodynamic simulations, taking into account the interstellar cooling. These molecular clouds show a peculiar shape composed of an arc-shaped cloud on one side of the TeV γ-ray source HESS J1023-575 and a linear distribution of clouds (jet clouds) on the other side. We propose that these clouds are formed by the interaction of a jet with clumps of interstellar neutral hydrogen (H i). By studying the dependence of the shape of dense cold clouds formed by shock compression and cooling on the filling factor of H i clumps, we found that the density distribution of H i clumps determines the shape of molecular clouds formed by the jet–cloud interaction: arc clouds are formed when the filling factor is large. On the other hand, when the filling factor is small, molecular clouds align with the jet. The jet propagates faster in models with small filling factors.

We investigate the physical properties of the inner gaseous disks of three hot Herbig B2e stars, HD 76534, HD 114981, and HD 216629, by modeling CFHT-ESPaDOns spectra using non-LTE radiative transfer codes. We assume that the emission lines are produced in a circumstellar disk heated solely by photospheric radiation from the central star in order to test whether the optical and near-infrared emission lines can be reproduced without invoking magnetospheric accretion. The inner gaseous disk density was assumed to follow a simple power-law in the equatorial plane, and we searched for models that could reproduce observed lines of H i (Hα and Hβ), He i, Ca ii, and Fe ii. For the three stars, good matches were found for all emission line profiles individually; however, no density model based on a single power-law was able to reproduce all of the observed emission lines. Among the single power-law models, the one with the gas density varying as ∼10⁻¹⁰(R_*/R)³ g cm⁻³ in the equatorial plane of a 25 R_* (0.78 au) disk did the best overall job of representing the optical emission lines of the three stars. This model implies a mass for the Hα-emitting portion of the inner gaseous disk of ∼10⁻⁹M_*. We conclude that the optical emission line spectra of these HBe stars can be qualitatively reproduced by a ≈1 au, geometrically thin, circumstellar disk of negligible mass compared to the central star in Keplerian rotation and radiative equilibrium.

Measuring the temperature structure of the solar atmosphere is critical to understanding how it is heated to high temperatures. Unfortunately, the temperature of the upper atmosphere cannot be observed directly, but must be inferred from spectrally resolved observations of individual emission lines that span a wide range of temperatures. Such observations are "inverted" to determine the distribution of plasma temperatures along the line of sight. This inversion is ill posed and, in the absence of regularization, tends to produce wildly oscillatory solutions. We introduce the application of sparse Bayesian inference to the problem of inferring the temperature structure of the solar corona. Within a Bayesian framework a preference for solutions that utilize a minimum number of basis functions can be encoded into the prior and many ad hoc assumptions can be avoided. We demonstrate the efficacy of the Bayesian approach by considering a test library of 40 assumed temperature distributions.

Using cosmological hydrodynamical simulations, we study the effect of supernova (SN) and active galactic nucleus (AGN) feedback on the mass transport (MT) of gas onto galactic nuclei and the black hole (BH) growth down to redshift $z\sim 6$. We study the BH growth in relation to the MT processes associated with gravity and pressure torques and how they are modified by feedback. Cosmological gas funneled through cold flows reaches the galactic outer region close to freefall. Then torques associated with pressure triggered by gas turbulent motions produced in the circumgalactic medium by shocks and explosions from SNe are the main source of MT beyond the central ∼100 pc. Due to high concentrations of mass in the central galactic region, gravitational torques tend to be more important at high redshift. The combined effect of almost freefalling material and both gravity and pressure torques produces a mass accretion rate of order $\sim 1\,{M}_{\odot }$ yr⁻¹ at approximately parsec scales. In the absence of SN feedback, AGN feedback alone does not affect significantly either star formation or BH growth until the BH reaches a sufficiently high mass of $\sim {10}^{6}\,{M}_{\odot }$ to self-regulate. SN feedback alone, instead, decreases both stellar and BH growth. Finally, SN and AGN feedback in tandem efficiently quench the BH growth, while star formation remains at the levels set by SN feedback alone, due to the small final BH mass, ∼few times ${10}^{5}\,{M}_{\odot }$. SNe create a more rarefied and hot environment where energy injection from the central AGN can accelerate the gas further.

An analysis of the physics-rich endgame of reionization at z = 5.7 is performed, jointly utilizing the observations of the Lyα forest, the mean free path (mfp) of ionizing photons, the luminosity function of galaxies, and new physical insight. We find that an upper limit on τ_e provides a constraint on the minimum mfp (of ionizing photons) that is primarily due to dwarf galaxies, which in turn yields a new, yet strongest constraint on the matter power spectrum on 10⁶–10⁹$\,{M}_{\odot }$ scales. With the latest Planck measurements of τ_e = 0.055 ± 0.009, we can place an upper limit of (8.9 × 10⁶, 3.8 × 10⁷, 4.2 × 10⁸)$\,{M}_{\odot }$ on the lower cutoff mass of the halo-mass function, or equivalently, a lower limit on warm dark matter particle mass m_x ≥ (15.1, 9.8, 4.6) keV or on sterile neutrino mass m_s ≥ (161, 90, 33) keV at  the (1, 1.4, 2.2)σ confidence level, respectively.

The stellar density distribution of the bulge is analyzed through one of its tracers. We use oxygen-rich Mira variables from the Catchpole et al. survey and OGLE-III survey as standard candles. The average age of these stars is around 9 Gyr. The population traced by Mira variables matches a boxy bulge prediction, not an X-shaped one, because only one peak is observed in the density along the analyzed lines of sight, whereas the prediction of an X-shape gives two clear peaks.

Magnetohydrodynamic (MHD) waves permeate the solar atmosphere and constitute potential coronal heating agents. Yet, the waves detected so far may be but a small subset of the true existing wave power. Detection is limited by instrumental constraints but also by wave processes that localize the wave power in undetectable spatial scales. In this study, we conduct 3D MHD simulations and forward modeling of standing transverse MHD waves in coronal loops with uniform and non-uniform temperature variation in the perpendicular cross-section. The observed signatures are largely dominated by the combination of the Kelvin–Helmholtz instability (KHI), resonant absorption, and phase mixing. In the presence of a cross-loop temperature gradient, we find that emission lines sensitive to the loop core catch different signatures compared to those that are more sensitive to the loop boundary and the surrounding corona, leading to an out-of-phase intensity and Doppler velocity modulation produced by KHI mixing. In all of the considered models, common signatures include an intensity and loop width modulation at half the kink period, a fine strand-like structure, a characteristic arrow-shaped structure in the Doppler maps, and overall line broadening in time but particularly at the loop edges. For our model, most of these features can be captured with a spatial resolution of 033 and a spectral resolution of 25 km s⁻¹, although we do obtain severe over-estimation of the line width. Resonant absorption leads to a significant decrease of the observed kinetic energy from Doppler motions over time, which is not recovered by a corresponding increase in the line width from phase mixing and KHI motions. We estimate this hidden wave energy to be a factor of 5–10 of the observed value.

We study particle acceleration and radiative processes in blazar jets under recurring conditions set by gravitational perturbations in supermassive binary systems. We consider the action from a companion orbiting a primary black hole of $\sim {10}^{8}\,{M}_{\odot }$, and perturbing its relativistic jet. We discuss how such conditions induce repetitive magneto-hydrodynamic stresses along the jet, and affect its inner electron acceleration and radiative processes. Specifically, we study how macroscopic perturbations related to increased jet "magnetization" end up in higher radiative outputs in the optical, X-ray, and gamma-ray bands. We first find an increase in magnetic field strength, as gauged in the optical band from the synchrotron emission of electrons accelerated in kinetic processes stimulated by reconnecting magnetic lines. The energetic electrons then proceed to up-scatter the synchrotron photons to GeV energies after the canonical synchrotron self-Compton radiation process. Our model implies a specific recurring pattern in the optical to gamma-ray emissions, made of high peaks and wide troughs. Progressing accelerations caused by spreading reconnections will produce an additional synchrotron keV component. Such outbursts provide a diagnostics for enhanced acceleration of electrons, which can up-scatter photons into the TeV range. We discuss how our model applies to the BL Lac object PG 1553+113, arguably the best candidate to now for high amplitude, recurring modulations in its gamma-ray emissions. We also consider other BL Lacs showing correlated keV–TeV radiations such as Mrk 421.

The key aspect determining the postformation luminosity of gas giants has long been considered to be the energetics of the accretion shock at the surface of the planet. We use one-dimensional radiation-hydrodynamical simulations to study the radiative loss efficiency and to obtain postshock temperatures and pressures and thus entropies. The efficiency is defined as the fraction of the total incoming energy flux that escapes the system (roughly the Hill sphere), taking into account the energy recycling that occurs ahead of the shock in a radiative precursor. We focus in this paper on a constant equation of state (EOS) to isolate the shock physics but use constant and tabulated opacities. While robust quantitative results will have to await a self-consistent treatment including hydrogen dissociation and ionization, the results presented here show the correct qualitative behavior and can be understood from semianalytical calculations. The shock is found to be isothermal and supercritical for a range of conditions relevant to the core accretion formation scenario (CA), with Mach numbers ${ \mathcal M }\gtrsim 3$. Across the shock, the entropy decreases significantly by a few times ${k}_{{\rm{B}}}/{\rm{baryon}}$. While nearly 100% of the incoming kinetic energy is converted to radiation locally, the efficiencies are found to be as low as roughly 40%, implying that a significant fraction of the total accretion energy is brought into the planet. However, for realistic parameter combinations in the CA scenario, we find that a nonzero fraction of the luminosity always escapes the Hill sphere. This luminosity could explain, at least in part, recent observations in the young LkCa 15 and HD 100546 systems.

A blue source in pre-explosion Hubble Space Telescope (HST)/Wide-Field Planetary Camera 2 (WFPC2) images falls within the 5σ astrometric error circle (∼024) derived from post-explosion ground-based imaging of SN 2010jl. At the time the ground-based astrometry was published, however, the SN had not faded sufficiently for post-explosion HST follow-up observations to determine a more precise astrometric solution and/or confirm if the pre-explosion source had disappeared, both of which are necessary to ultimately disentangle the possible progenitor scenarios. Here we present HST/WFC3 imaging of the SN 2010jl field obtained in 2014, 2015, and 2016 when the SN had faded sufficiently to allow for new constraints on the progenitor. The SN, which is still detected in the new images, is offset by 0061 ± 0008 (15 ± 2 pc) from the underlying and extended source of emission that contributes at least partially, if not entirely, to the blue source previously suggested as the candidate progenitor in the WFPC2 data. This point alone rules out the possibility that the blue source in the pre-explosion images is the exploding star, but may instead suggest an association with a young (<5–6 Myr) cluster and still argues for a massive (>30 M_⊙) progenitor. We obtain new upper limits on the flux from a single star at the SN position in the pre-explosion WFPC2 and Spitzer/IRAC images that may ultimately be used to constrain the progenitor properties.

Multi-wavelength observations provide a complementary view of the formation of young, directly imaged planet-mass companions. We report the ALMA 1.3 mm and Magellan adaptive optics Hα, $i^{\prime}$, $z^{\prime}$, and Y_S observations of the GQ Lup system, a classical T Tauri star with a $10\mbox{--}40\,{M}_{\mathrm{Jup}}$ substellar companion at ∼110 au projected separation. We estimate the accretion rates for both components from the observed Hα fluxes. In our ∼005 resolution ALMA map, we resolve GQ Lup A's disk in the dust continuum, but no signal is found from the companion. The disk is compact, with a radius of ∼22 au, a dust mass of ∼6 M_⊕, an inclination angle of ∼56°, and a very flat surface density profile indicative of a radial variation in dust grain sizes. No gaps or inner cavity are found in the disk, so there is unlikely a massive inner companion to scatter GQ Lup B outward. Thus, GQ Lup B might have formed in situ via disk fragmentation or prestellar core collapse. We also show that GQ Lup A's disk is misaligned with its spin axis, and possibly with GQ Lup B's orbit. Our analysis on the tidal truncation radius of GQ Lup A's disk suggests that GQ Lup B's orbit might have a low eccentricity.

The phenomenon of subpulse drifting may hold the key to understanding the pulsar emission mechanism. Here, we report on new observations of PSR J0034−0721 (B0031−07) carried out with the Murchison Widefield Array at $185\,\mathrm{MHz}$. We observe three distinct drift modes whose "vertical" drift band separations (P₃) and relative abundances are consistent with previous studies at similar and higher frequencies. The drift bands, however, are observed to change their slopes over the course of individual drift modes, which can be interpreted as a continuously changing drift rate. The implied acceleration of the intrinsic carousel rotation cannot easily be explained by plasma models based on ${\boldsymbol{E}}\times {\boldsymbol{B}}$ drift. Furthermore, we find that methods of classifying the drift modes by means of P₃ measurements can sometimes produce erroneous identifications in the presence of a changing drift rate. The "horizontal" separation between drift bands (P₂) is found to be larger at later rotation phases within the pulse window, which is inconsistent with the established effects of retardation, aberration, and the motion of the visible point. Longer observations spanning at least ∼10,000 pulses are required to determine how the carousel rotation parameters change from one drift sequence to the next.

We found a simultaneous decrease of the Fe–K line and 4.2–6 keV continuum of Cassiopeia A with the monitoring data taken by the Chandra X-ray Observatory in 2000–2013. The flux change rates in the whole remnant are −0.65 ± 0.02% yr⁻¹ in the 4.2–6.0 keV continuum and −0.6 ± 0.1% yr⁻¹ in the Fe–K line. In the eastern region where the thermal emission is considered to dominate, the variations show the largest values: −1.03 ± 0.05% yr⁻¹ (4.2–6 keV band) and −0.6 ± 0.1% yr⁻¹ (Fe–K line). In this region, the time evolution of the emission measure and the temperature have a decreasing trend. This could be interpreted as adiabatic cooling with the expansion of m = 0.66. On the other hand, in the non-thermal emission dominated regions, variations of the 4.2–6 keV continuum show smaller rates: −0.60 ± 0.04% yr⁻¹ in the southwestern region, −0.46 ± 0.05% yr⁻¹ in the inner region, and +0.00 ± 0.07% yr⁻¹ in the forward shock region. In particular, flux does not show significant change in the forward shock region. These results imply that strong braking in shock velocity has not been occurring in Cassiopeia A (<5 km s⁻¹ yr⁻¹). All of our results support the idea that X-ray flux decay in the remnant is mainly caused by thermal components.

FU Orionis-type objects (FUors) are young stellar objects experiencing large optical outbursts due to highly enhanced accretion from the circumstellar disk onto the star. FUors are often surrounded by massive envelopes, which play a significant role in the outburst mechanism. Conversely, the subsequent eruptions might gradually clear up the obscuring envelope material and drive the protostar on its way to become a disk-only T Tauri star. Here we present an APEX ¹²CO and ¹³CO survey of eight southern and equatorial FUors. We measure the mass of the gaseous material surrounding our targets, locate the source of the CO emission, and derive physical parameters for the envelopes and outflows, where detected. Our results support the evolutionary scenario where FUors represent a transition phase from envelope-surrounded protostars to classical T Tauri stars.

The composition of Jupiter and the primordial distribution of the heavy elements are determined by its formation history. As a result, in order to constrain the primordial internal structure of Jupiter, the growth of the core and the deposition and settling of accreted planetesimals must be followed in detail. In this paper we determine the distribution of the heavy elements in proto-Jupiter and determine the mass and composition of the core. We find that while the outer envelope of proto-Jupiter is typically convective and has a homogeneous composition, the innermost regions have compositional gradients. In addition, the existence of heavy elements in the envelope leads to much higher internal temperatures (several times 10⁴ K) than in the case of a hydrogen–helium envelope. The derived core mass depends on the actual definition of the core: if the core is defined as the region in which the heavy-element mass fraction is above some limit (say, 0.5), then it can be much more massive (∼15 ${M}_{\oplus }$) and more extended (10% of the planet's radius) than in the case where the core is just the region with 100% heavy elements. In the former case Jupiter's core also consists of hydrogen and helium. Our results should be taken into account when constructing internal structure models of Jupiter and when interpreting the upcoming data from the Juno (NASA) mission.

We study logamediate inflation in the context of f(T) teleparallel gravity. f(T)-gravity is a generalization of the teleparallel gravity which is formulated on the Weitzenbock spacetime, characterized by the vanishing curvature tensor (absolute parallelism) and the non-vanishing torsion tensor. We consider an f(T)-gravity model which is sourced by a canonical scalar field. Assuming a power-law f(T) function in the action, we investigate an inflationary universe with a logamediate scale factor. Our results show that, although logamediate inflation is completely ruled out by observational data in the standard inflationary scenario based on Einstein gravity, it can be compatible with the 68% confidence limit joint region of Planck 2015 TT,TE,EE+lowP data in the framework of f(T)-gravity.

We analyze the spectrum of the 11.2 μm unidentified infrared band (UIR) from NGC 7027 and identify a small fullerene (C₂₄) as a plausible carrier. The blurring effects of lifetime and vibrational anharmonicity broadening obscure the narrower, intrinsic spectral profiles of the UIR band carriers. We use a spectral deconvolution algorithm to remove the blurring, in order to retrieve the intrinsic profile of the UIR band. The shape of the intrinsic profile—a sharp blue peak and an extended red tail—suggests that the UIR band originates from a molecular vibration–rotation band with a blue band head. The fractional area of the band-head feature indicates a spheroidal molecule, implying a nonpolar molecule and precluding rotational emission. Its rotational temperature should be well approximated by that measured for nonpolar molecular hydrogen, ∼825 K for NGC 7027. Using this temperature, and the inferred spherical symmetry, we perform a spectral fit to the intrinsic profile, which results in a rotational constant implying C₂₄ as the carrier. We show that the spectroscopic parameters derived for NGC 7027 are consistent with the 11.2 μm UIR bands observed for other objects. We present density functional theory (DFT) calculations for the frequencies and infrared intensities of C₂₄. The DFT results are used to predict a spectral energy distribution (SED) originating from absorption of a 5 eV photon, and characterized by an effective vibrational temperature of 930 K. The C₂₄ SED is consistent with the entire UIR spectrum and is the dominant contributor to the 11.2 and 12.7 μm bands.

The role of compact binary mergers as the main production site of r-process elements is investigated by combining stellar abundances of Eu observed in the Milky Way, galactic chemical evolution (GCE) simulations, and binary population synthesis models, and gravitational wave measurements from Advanced LIGO. We compiled and reviewed seven recent GCE studies to extract the frequency of neutron star–neutron star (NS–NS) mergers that is needed in order to reproduce the observed [Eu/Fe] versus [Fe/H] relationship. We used our simple chemical evolution code to explore the impact of different analytical delay-time distribution functions for NS–NS mergers. We then combined our metallicity-dependent population synthesis models with our chemical evolution code to bring their predictions, for both NS–NS mergers and black hole–neutron star mergers, into a GCE context. Finally, we convolved our results with the cosmic star formation history to provide a direct comparison with current and upcoming Advanced LIGO measurements. When assuming that NS–NS mergers are the exclusive r-process sites, and that the ejected r-process mass per merger event is 0.01 M${}_{\odot }$, the number of NS–NS mergers needed in GCE studies is about 10 times larger than what is predicted by standard population synthesis models. These two distinct fields can only be consistent with each other when assuming optimistic rates, massive NS–NS merger ejecta, and low Fe yields for massive stars. For now, population synthesis models and GCE simulations are in agreement with the current upper limit (O1) established by Advanced LIGO during their first run of observations. Upcoming measurements will provide an important constraint on the actual local NS–NS merger rate, will provide valuable insights on the plausibility of the GCE requirement, and will help to define whether or not compact binary mergers can be the dominant source of r-process elements in the universe.

The lensing signal around galaxy clusters can, in principle, be used to test detailed predictions for their average mass profile from numerical simulations. However, the intrinsic shape of the profiles can be smeared out when a sample that spans a wide range of cluster masses is averaged in physical length units. This effect especially conceals rapid changes in gradient such as the steep drop associated with the splashback radius, a sharp edge corresponding to the outermost caustic in accreting halos. We optimize the extraction of such local features by scaling individual halo profiles to a number of spherical overdensity radii, and apply this method to 16 X-ray-selected, high-mass clusters targeted in the Cluster Lensing And Supernova survey with Hubble. By forward-modeling the weak- and strong-lensing data presented by Umetsu et al., we show that, regardless of the scaling overdensity, the projected ensemble density profile is remarkably well described by a Navarro–Frenk–White (NFW) or Einasto profile out to $R\sim 2.5\,{h}^{-1}\,\mathrm{Mpc}$, beyond which the profiles flatten. We constrain the NFW concentration to ${c}_{200{\rm{c}}}=3.66\pm 0.11$ at ${M}_{200{\rm{c}}}\simeq 1.0\times {10}^{15}\,{h}^{-1}\,{M}_{\odot }$, consistent with and improved from previous work that used conventionally stacked lensing profiles, and in excellent agreement with theoretical expectations. Assuming the profile form of Diemer & Kravtsov and generic priors calibrated from numerical simulations, we place a lower limit on the splashback radius of the cluster halos, if it exists, of ${R}_{\mathrm{sp}}^{3{\rm{D}}}/{r}_{200{\rm{m}}}\gt 0.89$ (${R}_{\mathrm{sp}}^{3{\rm{D}}}\gt 1.83\,{h}^{-1}\,\mathrm{Mpc}$) at 68% confidence. The corresponding density feature is most pronounced when the cluster profiles are scaled by ${r}_{200{\rm{m}}}$, and smeared out when scaled to higher overdensities.

The intrinsic colors of Type Ia supernovae (SNe Ia) are important to understanding their use as cosmological standard candles. Understanding the effects of reddening and redshift on the observed colors are complicated and dependent on the intrinsic spectrum, the filter curves, and the wavelength dependence of reddening. We present ultraviolet and optical data of a growing sample of SNe Ia observed with the Ultraviolet/Optical Telescope on the Swift spacecraft and use this sample to re-examine the near-UV (NUV) colors of SNe Ia. We find that a small amount of reddening (E(B − V) = 0.2 mag) could account for the difference between groups designated as NUV-blue and NUV-red, and a moderate amount of reddening (E(B − V) = 0.5 mag) could account for the whole NUV-optical differences. The reddening scenario, however, is inconsistent with the mid-UV colors and color evolution. The effect of redshift alone only accounts for part of the variation. Using a spectral template of SN2011fe, we can forward model the effects of redshift and reddening and directly compare those with the observed colors. We find that some SNe are consistent with reddened versions of SN2011fe, but most SNe Ia are much redder in the uvw1 − v color than SN2011fe reddened to the same b − v color. The absolute magnitudes show that two out of five NUV-blue SNe Ia are blue because their near-UV luminosity is high, and the other three are optically fainter. We also show that SN 2011fe is not a "normal" SN Ia in the UV, but has colors placing it at the blue extreme of our sample.

The multiteraelectronvolt γ-rays from the galactic center (GC) have a cutoff at tens of teraelectronvolts, whereas the diffuse emission has no such cutoff, which is regarded as an indication of petaelectronvolt proton acceleration by the HESS experiment. It is important to understand the inconsistency and study the possibility that petaelectronvolt cosmic-ray acceleration could account for the apparently contradictory point and diffuse γ-ray spectra. In this work, we propose that the cosmic rays are accelerated up to greater than petaelectronvolts in the GC. The interaction between cosmic rays and molecular clouds is responsible for the multiteraelectronvolt γ-ray emissions from both the point and diffuse sources today. Enhanced by the small volume filling factor (VFF) of the clumpy structure, the absorption of the γ-rays leads to a sharp cutoff spectrum at tens of teraelectronvolts produced in the GC. Away from the GC, the VFF grows, and the absorption enhancement becomes negligible. As a result, the spectra of γ-ray emissions for both point and diffuse sources can be successfully reproduced under such a self-consistent picture. In addition, a "surviving tail" at ∼100 TeV is expected from the point source, which can be observed by future projects CTA and LHAASO. Neutrinos are simultaneously produced during proton-proton (PP) collision. With 5–10 years of observations, the KM3Net experiment will be able to detect the petaelectronvolt source according to our calculation.

Stellar streams result from the tidal disruption of satellites and star clusters as they orbit a host galaxy, and can be very sensitive probes of the gravitational potential of the host system. We select and study narrow stellar streams formed in a Milky-Way-like dark matter halo of the Aquarius suite of cosmological simulations, to determine if these streams can be used to constrain the present day characteristic parameters of the halo's gravitational potential. We find that orbits integrated in both spherical and triaxial static Navarro–Frenk–White potentials reproduce the locations and kinematics of the various streams reasonably well. To quantify this further, we determine the best-fit potential parameters by maximizing the amount of clustering of the stream stars in the space of their actions. We show that using our set of Aquarius streams, we recover a mass profile that is consistent with the spherically averaged dark matter profile of the host halo, although we ignored both triaxiality and time evolution in the fit. This gives us confidence that such methods can be applied to the many streams that will be discovered by the Gaia mission to determine the gravitational potential of our Galaxy.

We present the detection of a small eruption and the associated quasi-circular ribbon flare during the emergence of a bipole occurring on 2015 February 3. Under a fan dome, a sigmoid was rooted in a single magnetic bipole, which was encircled by negative polarity. The nonlinear force-free field extrapolation shows the presence of twisted field lines, which can represent a sigmoid structure. The rotation of the magnetic bipole may cause the twisting of magnetic field lines. An initial brightening appeared at one of the footpoints of the sigmoid, where the positive polarity slides toward a nearby negative polarity field region. The sigmoid displayed an ascending motion and then interacted intensively with the spine-like field. This type of null point reconnection in corona led to a violent blowout jet, and a quasi-circular flare ribbon was also produced. The magnetic emergence and rotational motion are the main contributors to the energy buildup for the flare, while the cancellation and collision might act as a trigger.

Exoplanet transit spectroscopy enables the characterization of distant worlds, and will yield key results for NASA's James Webb Space Telescope. However, transit spectra models are often simplified, omitting potentially important processes like refraction and multiple scattering. While the former process has seen recent development, the effects of light multiple scattering on exoplanet transit spectra have received little attention. Here, we develop a detailed theory of exoplanet transit spectroscopy that extends to the full refracting and multiple scattering case. We explore the importance of scattering for planet-wide cloud layers, where the relevant parameters are the slant scattering optical depth, the scattering asymmetry parameter, and the angular size of the host star. The latter determines the size of the "target" for a photon that is back-mapped from an observer. We provide results that straightforwardly indicate the potential importance of multiple scattering for transit spectra. When the orbital distance is smaller than 10–20 times the stellar radius, multiple scattering effects for aerosols with asymmetry parameters larger than 0.8–0.9 can become significant. We provide examples of the impacts of cloud/haze multiple scattering on transit spectra of a hot Jupiter-like exoplanet. For cases with a forward and conservatively scattering cloud/haze, differences due to multiple scattering effects can exceed 200 ppm, but shrink to zero at wavelength ranges corresponding to strong gas absorption or when the slant optical depth of the cloud exceeds several tens. We conclude with a discussion of types of aerosols for which multiple scattering in transit spectra may be important.

We explore the nucleus of the nearby 10⁹${M}_{\odot }$ early-type galaxy, NGC 404, using Hubble Space Telescope (HST)/STIS spectroscopy and WFC3 imaging. We first present evidence for nuclear variability in UV, optical, and infrared filters over a time period of 15 years. This variability adds to the already substantial evidence for an accreting black hole at the center of NGC 404. We then redetermine the dynamical black hole mass in NGC 404 including modeling of the nuclear stellar populations. We combine HST/STIS spectroscopy with WFC3 images to create a local color–M/L relation derived from stellar population modeling of the STIS data. We then use this to create a mass model for the nuclear region. We use Jeans modeling to fit this mass model to adaptive optics stellar kinematic observations from Gemini/NIFS. From our stellar dynamical modeling, we find a 3σ upper limit on the black hole mass of $1.5\times {10}^{5}$${M}_{\odot }$. Given the accretion evidence for a black hole, this upper limit makes NGC 404 the lowest mass central black hole with dynamical mass constraints. We find that the kinematics of H₂ emission line gas show evidence for non-gravitational motions preventing the use of gas dynamical modeling to constrain the black hole mass. Our stellar population modeling also reveals that the central, counter-rotating region of the nuclear cluster is dominated by ∼1 Gyr old populations.

Interstellar Boundary Explorer (IBEX) measurements of energetic neutral atoms (ENAs) from the heliotail show a multi-lobe structure of ENA fluxes as a function of energy between ∼0.71 and 4.29 keV. Below ∼2 keV, there is a single structure of enhanced ENA fluxes centered near the downwind direction. Above ∼2 keV, this structure separates into two lobes, one north and one south of the solar equatorial plane. ENA flux from these two lobes can be interpreted as originating from the fast solar wind (SW) propagating through the inner heliosheath (IHS). Alternatively, a recently published model of the heliosphere suggests that the heliotail may split into a "croissant-like" shape, and that such a geometry could be responsible for the heliotail ENA feature. Here we present results from a time-dependent simulation of the heliosphere that produces a comet-like heliotail, and show that the 11-year solar cycle leads to the formation of ENA lobes with properties remarkably similar to those observed by IBEX. The ENA energy at which the north and south lobes appear suggests that the pickup ion (PUI) temperature in the slow SW of the IHS is ∼10⁷ K. Moreover, we demonstrate that the extinction of PUIs by charge-exchange is an essential process required to create the observed global ENA structure. While the shape and locations of the ENA lobes as a function of energy are well reproduced by PUIs that cross the termination shock, the results appear to be sensitive to the form of the distribution of PUIs injected in the IHS.

The selection of high-redshift sources from broadband photometry using the Lyman-break galaxy (LBG) technique is a well established methodology, but the characterization of its contamination for the faintest sources is still incomplete. We use the optical and near-IR data from four (ultra)deep Hubble Space Telescope legacy fields to investigate the contamination fraction of LBG samples at $z\sim 5\mbox{--}8$ selected using a color–color method. Our approach is based on characterizing the number count distribution of interloper sources, that is, galaxies with colors similar to those of LBGs, but showing detection at wavelengths shorter than the spectral break. Without sufficient sensitivity at bluer wavelengths, a subset of interlopers may not be properly classified, and contaminate the LBG selection. The surface density of interlopers in the sky gets steeper with increasing redshift of LBG selections. Since the intrinsic number of dropouts decreases significantly with increasing redshift, this implies increasing contamination from misclassified interlopers with increasing redshift, primarily by intermediate redshift sources with unremarkable properties (intermediate ages, lack of ongoing star formation and low/moderate dust content). Using Monte-Carlo simulations, we estimate that the CANDELS deep data have contamination induced by photometric scatter increasing from $\sim 2 \%$ at $z\sim 5$ to $\sim 6 \%$ at $z\sim 8$ for a typical dropout color $\geqslant 1$ mag, with contamination naturally decreasing for a more stringent dropout selection. Contaminants are expected to be located preferentially near the detection limit of surveys, ranging from 0.1 to 0.4 contaminants per arcmin² at ${J}_{125}$ = 30, depending on the field considered. This analysis suggests that the impact of contamination in future studies of $z\gt 10$ galaxies needs to be carefully considered.

Chemistry in the interstellar medium (ISM) is capable of producing complex organic molecules (COMs) of great importance to astrobiology. Gas phase and grain surface chemistry almost certainly both contribute to COM formation. Amino acids as building blocks of proteins are some of the most interesting COMs. The simplest one, glycine, has been characterized in meteorites and comets and, its conclusive detection in the ISM seems to be highly plausible. In this work, we analyze the gas phase reaction of glycine and ${{\mathrm{CH}}_{5}}^{+}$ to establish the role of this process in the formation of alanine or other COMs in the ISM. Formation of protonated α- and β-alanine in spite of being exothermic processes is not viable under interstellar conditions because the different paths leading to these isomers present net activation energies. Nevertheless, glycine can evolve to protonated 1-imide-2, 2-propanediol, protonated amino acetone, protonated hydroxyacetone, and protonated propionic acid. However, formation of acetic acid and protonated methylamine is also a favorable process and therefore will be a competitive channel with the evolution of glycine to COMs.

PSR J2032+4127 is a radio-loud gamma-ray-emitting pulsar; it is orbiting around a high-mass Be type star with a very long orbital period of 25–50 years, and is approaching periastron, which will occur in late 2017/early 2018. This system comprises a young pulsar and a Be type star, which is similar to the so-called gamma-ray binary PSR B1259–63/LS2883. It is expected therefore that PSR J2032+4127 shows an enhancement of high-energy emission caused by the interaction between the pulsar wind and Be wind/disk around periastron. Ho et al. recently reported a rapid increase in the X-ray flux from this system. In this paper, we also confirm a rapid increase in the X-ray flux along the orbit, while the GeV flux shows no significant change. We discuss the high-energy emissions from the shock caused by the pulsar wind and stellar wind interaction and examine the properties of the pulsar wind in this binary system. We argue that the rate of increase of the X-ray flux observed by Swift indicates (1) a variation of the momentum ratio of the two-wind interaction region along the orbit, or (2) an evolution of the magnetization parameter of the pulsar wind with the radial distance from the pulsar. We also discuss the pulsar wind/Be disk interaction at the periastron passage, and propose the possibility of formation of an accretion disk around the pulsar. We model high-energy emissions through the inverse-Compton scattering process of the cold-relativistic pulsar wind off soft photons from the accretion disk.

Residual gas in disks around young stars can spin down stars, circularize the orbits of terrestrial planets, and whisk away the dusty debris that is expected to serve as a signpost of terrestrial planet formation. We have carried out a sensitive search for residual gas and dust in the terrestrial planet region surrounding young stars ranging in age from a few to ∼10 Myr. Using high-resolution 4.7 μm spectra of transition objects (TOs) and weak T Tauri stars, we searched for weak continuum excesses and CO fundamental emission, after making a careful correction for the stellar contribution to the observed spectrum. We find that the CO emission from TOs is weaker and located farther from the star than CO emission from nontransition T Tauri stars with similar stellar accretion rates. The difference is possibly the result of chemical and/or dynamical effects (i.e., a low CO abundance or close-in low-mass planets). The weak T Tauri stars show no CO fundamental emission down to low flux levels (5 × 10⁻²⁰ to 10⁻¹⁸ W m⁻²). We illustrate how our results can be used to constrain the residual disk gas content in these systems and discuss their potential implications for star and planet formation.

Hot, million degree gas appears to pervade the Milky Way halo, containing a large fraction of the Galactic missing baryons. This circumgalactic medium (CGM) is probed effectively in X-rays, both in absorption and in emission. The CGM also appears to be anisotropic, so we have started a program to determine CGM properties along several sightlines by combining absorption and emission measurements. Here we present the emission measure close to the Mrk 509 sightline using new Suzaku and XMM-Newton observations. We also present new analysis and modeling of Chandra HETG spectra to constrain the absorption parameters. The emission measure in this sightline is high, EM = 0.0165 ± 0.0008 ± 0.0006 cm⁻⁶ pc, five times larger than the average. The observed O vii column density N(O vii) $=\,2.35\pm 0.4\times {10}^{16}$ cm⁻², however, is close to the average. We find that the temperature of the emitting and absorbing gas is the same: logT(K) = 6.33 ± 0.01 and logT(K) = 6.33 ± 0.16 respectively. We fit the observed column density and emission measure with a β-model density profile. The central density is constrained to be between n₀ = 2.8–6.0 × 10⁻⁴ cm⁻³ and the core radius of the density profile has a lower limit of 40 kpc. This shows that the hot gas is mostly in the CGM of the galaxy, not in the Galactic disk. Our derived density profile is close to the Maller & Bullock profile for adiabatic gas in hydrostatic equilibrium with an NFW dark-matter potential well. Assuming this density profile, the minimum mass of the hot CGM is $3.2\times {10}^{10}\,$M_⊙.

The final evolution of stars in the mass range 70–140 ${\text{}}{M}_{\odot }$ is explored. Depending upon their mass loss history and rotation rates, these stars will end their lives as pulsational pair-instability supernovae (PPISN) producing a great variety of observational transients with total durations ranging from weeks to millennia and luminosities from 10⁴¹ to over 10⁴⁴ erg s⁻¹. No nonrotating model radiates more than $5\times {10}^{50}$ erg of light or has a kinetic energy exceeding $5\times {10}^{51}$ erg, but greater energies are possible, in principle, in magnetar-powered explosions, which are explored. Many events resemble SNe Ibn, SNe Icn, and SNe IIn, and some potential observational counterparts are mentioned. Some PPISN can exist in a dormant state for extended periods, producing explosions millennia after their first violent pulse. These dormant supernovae contain bright Wolf–Rayet stars, possibly embedded in bright X-ray and radio sources. The relevance of PPISN to supernova impostors like Eta Carinae, to superluminous supernovae, and to sources of gravitational radiation is discussed. No black holes between 52 and 133 ${\text{}}{M}_{\odot }$ are expected from stellar evolution in close binaries.

Recent simulations on super-Eddington accretion flows have shown that, apart from the diffusion process, the vertical advection based on magnetic buoyancy can be a more efficient process to release the trapped photons in the optically thick disk. As a consequence, the radiative luminosity from the accretion disk can be far beyond the Eddington value. Following this spirit, we revisit the structure and radiation of hyper-accretion disks with mass accretion rates in the range of ${10}^{-3}\sim 10\,{M}_{\odot }\,{{\rm{s}}}^{-1}$. Our results show that, due to the strong cooling through the vertical advection, the disk temperature becomes lower than that in the classic model without the vertical advection process, and therefore the neutrino luminosity from the disk is lower. On the other hand, the gamma-ray photons released through the vertical advection can be extremely super-Eddington. We argue that the large amount of escaped gamma-ray photons may have more significant contribution to the primordial fireball than the neutrino annihilation, and may hint at a link between gamma-ray bursts and kilonovae in the black hole hyper-accretion scenario.

We investigate the nature of the outer envelope of halo coronal mass ejections (H-CMEs) using multi-viewpoint observations from the Solar Terrestrial Relations Observatory-A, -B, and SOlar and Heliospheric Observatory coronagraphs. The 3D structure and kinematics of the halo envelopes and the driving CMEs are derived separately using a forward modeling method. We analyze three H-CMEs with peak speeds from 1355 to 2157 km s⁻¹; sufficiently fast to drive shocks in the corona. We find that the angular widths of the halos range from 192° to 252°, while those of the flux ropes range between only 58° and 91°, indicating that the halos are waves propagating away from the CMEs. The halo widths are in agreement with widths of Extreme Ultraviolet (EUV) waves in the low corona further demonstrating the common origin of these structures. To further investigate the wave nature of the halos, we model their 3D kinematic properties with a linear fast magnetosonic wave model. The model is able to reproduce the position of the halo flanks with realistic coronal medium assumptions but fails closer to the CME nose. The CME halo envelope seems to arise from a driven wave (or shock) close to the CME nose, but it is gradually becoming a freely propagating fast magnetosonic wave at the flanks. This interpretation provides a simple unifying picture for CME halos, EUV waves, and the large longitudinal spread of solar energetic particles.

We present a statistical study of coherent structures at kinetic scales, using data from the Magnetospheric Multiscale mission in the Earth's magnetosheath. We implemented the multi-spacecraft partial variance of increments (PVI) technique to detect these structures, which are associated with intermittency at kinetic scales. We examine the properties of the electron heating occurring within such structures. We find that, statistically, structures with a high PVI index are regions of significant electron heating. We also focus on one such structure, a current sheet, which shows some signatures consistent with magnetic reconnection. Strong parallel electron heating coincides with whistler emissions at the edges of the current sheet.