Table of contents

We investigate galaxy formation in models with dark matter (DM) constituted by sterile neutrinos. Given their large parameter space, defined by the combinations of sterile neutrino mass ${m}_{\nu }$ and mixing parameter ${\sin }^{2}(2\theta )$ with active neutrinos, we focus on models with ${m}_{\nu }=7\,\mathrm{keV}$, consistent with the tentative 3.5 keV line detected in several X-ray spectra of clusters and galaxies. We consider (1) two resonant production models with ${\sin }^{2}(2\theta )=5\,\times \,{10}^{-11}$ and ${\sin }^{2}(2\theta )=2\,\times \,{10}^{-10}$, to cover the range of mixing parameters consistent with the 3.5 keV line; (2) two scalar-decay models, representative of the two possible cases characterizing such a scenario: a freeze-in and a freeze-out case. We also consider thermal warm DM with particle mass ${m}_{X}=3\,\mathrm{keV}$. Using a semianalytic model, we compare the predictions for the different DM scenarios with a wide set of observables. We find that comparing the predicted evolution of the stellar mass function, the abundance of satellites of Milky Way–like galaxies, and the global star formation history of galaxies with observations does not allow us to disentangle the effects of the baryonic physics from those related to the different DM models. On the other hand, the distribution of the stellar-to-halo mass ratios, the abundance of faint galaxies in the UV luminosity function at $z\gtrsim 6$, and the specific star formation and age distribution of local, low-mass galaxies constitute potential probes for the DM scenarios considered. We discuss how future observations with upcoming facilities will enable us to rule out or to strongly support DM models based on sterile neutrinos.

The leading models for launching astrophysical jets rely on strong poloidal magnetic fields threading the central parts of their host accretion disks. Numerical simulations of magneto-rotationally turbulent disks suggest that such fields are actually advected from the environment by the accreting matter rather than generated by internal dynamos. This is puzzling from a theoretical point of view, since the reconnection of the radial field across the midplane should cause an outward drift on timescales much shorter than the accretion time. We suggest that a combination of effects are responsible for reducing the radial field near the midplane, causing efficient inward advection of the poloidal field. Magnetic buoyancy in subsonic turbulence pushes the field lines away from the midplane, decreasing the large-scale radial field in the main body of the disk. In magneto-rotationally driven turbulence, magnetic buoyancy dominates over the effects of turbulent pumping, which works against it, and turbulent diamagnetism, which works with it, in determining the vertical drift of the magnetic field. Balancing buoyancy with diffusion implies that the bending angle of the large-scale poloidal field can be very large near the surface, as required for outflows, but vanishes near the midplane, which impedes turbulent reconnection and outward diffusion. This effect becomes less efficient as the poloidal flux increases. This suggests that accretion disks are less likely to form jets if they have a modest ratio of outer to inner radii or if the ambient field is very weak. The former effect is probably responsible for the scarcity of jets in cataclysmic variable systems.

The ice-albedo feedback on rapidly rotating terrestrial planets in the habitable zone can lead to abrupt transitions (bifurcations) between a warm and a snowball (ice-covered) state, bistability between these states, and hysteresis in planetary climate. This is important for planetary habitability because snowball events may trigger rises in the complexity of life, but could also endanger complex life that already exists. Recent work has shown that planets tidally locked in synchronous rotation states will transition smoothly into the snowball state rather than experiencing bifurcations. Here we investigate the structure of snowball bifurcations on planets that are tidally influenced, but not synchronously rotating, so that they experience long solar days. We use PlaSIM, an intermediate-complexity global climate model, with a thermodynamic mixed layer ocean and the Sun's spectrum. We find that the amount of hysteresis (the range in stellar flux for which there is bistability in climate) is significantly reduced for solar days with lengths of tens of Earth days, and disappears for solar days of hundreds of Earth days. These results suggest that tidally influenced planets orbiting M and K stars that are not synchronously rotating could have much less hysteresis associated with the snowball bifurcations than they would if they were rapidly rotating. This implies that the amount of time it takes them to escape a snowball state via CO₂ outgassing would be greatly reduced, as would the period of cycling between the warm and snowball state if they have low CO₂ outgassing rates.

We explore the expected X-ray (0.5–2 keV) signatures from supermassive black holes (SMBHs) at high redshifts ($z\sim 5\mbox{--}12$) assuming various models for their seeding mechanism and evolution. Seeding models are approximated through deviations from the ${M}_{\mathrm{BH}}\mbox{--}\sigma$ relation observed in the local universe, while N-body simulations of the large-scale structure are used to estimate the density of observable SMBHs. We focus on two seeding model families: (i) light seed BHs from remnants of Pop-III stars and (ii) heavy seeds from the direct collapse of gas clouds. We investigate several models for the accretion history, such as sub-Eddington accretion, slim disk models, and torque-limited growth models. We consider observations with two instruments: (i) the Chandra X-ray Observatory and (ii) the proposed Lynx. We find that all of the simulated models are in agreement with the current results from the Chandra Deep Field South, i.e., consistent with zero SMBHs in the field of view. In deep Lynx exposures, the number of observed objects is expected to become statistically significant. We demonstrate the capability to limit the phase space of plausible scenarios of the birth and evolution of SMBHs by performing deep observations at a flux limit of $1\times {10}^{-19}\,\mathrm{erg}\,{\mathrm{cm}}^{-2}\,{{\rm{s}}}^{-1}$. Finally, we show that our models are in agreement with current limits on the cosmic X-ray background (CXRB) and the expected contribution from unresolved quasars. We find that an analysis of CXRB contributions down to the Lynx confusion limit yields valuable information that can help identify the correct scenario for the birth and evolution of SMBHs.

Numerical solutions of the cosmic-ray (CR) magnetohydrodynamic equations are dogged by a powerful numerical instability, which arises from the constraint that CRs can only stream down their gradient. The standard cure is to regularize by adding artificial diffusion. Besides introducing ad hoc smoothing, this has a significant negative impact on either computational cost or complexity and parallel scalings. We describe a new numerical algorithm for CR transport, with close parallels to two-moment methods for radiative transfer under the reduced speed of light approximation. It stably and robustly handles CR streaming without any artificial diffusion. It allows for both isotropic and field-aligned CR streaming and diffusion, with arbitrary streaming and diffusion coefficients. CR transport is handled explicitly, while source terms are handled implicitly. The overall time step scales linearly with resolution (even when computing CR diffusion) and has a perfect parallel scaling. It is given by the standard Courant condition with respect to a constant maximum velocity over the entire simulation domain. The computational cost is comparable to that of solving the ideal MHD equation. We demonstrate the accuracy and stability of this new scheme with a wide variety of tests, including anisotropic streaming and diffusion tests, CR-modified shocks, CR-driven blast waves, and CR transport in multiphase media. The new algorithm opens doors to much more ambitious and hitherto intractable calculations of CR physics in galaxies and galaxy clusters. It can also be applied to other physical processes with similar mathematical structure, such as saturated, anisotropic heat conduction.

We study a truncated accretion disk using a well-resolved, semi-global magnetohydrodynamic simulation that is evolved for many dynamical times (6096 inner disk orbits). The spectral properties of hard-state black hole binary systems and low-luminosity active galactic nuclei are regularly attributed to truncated accretion disks, but a detailed understanding of the flow dynamics is lacking. In these systems the truncation is expected to arise through thermal instability driven by sharp changes in the radiative efficiency. We emulate this behavior using a simple bistable cooling function with efficient and inefficient branches. The accretion flow takes on an arrangement where a "transition zone" exists in between hot gas in the innermost regions and a cold, Shakura & Sunyaev thin disk at larger radii. The thin disk is embedded in an atmosphere of hot gas that is fed by a gentle outflow originating from the transition zone. Despite the presence of hot gas in the inner disk, accretion is efficient. Our analysis focuses on the details of the angular momentum transport, energetics, and magnetic field properties. We find that the magnetic dynamo is suppressed in the hot, truncated inner region of the disk which lowers the effective α-parameter by 65%.

Understanding the origins of stellar radio emission can provide invaluable insight into the strength and geometry of stellar magnetic fields and the resultant space weather environment experienced by exoplanets. Here, we present the first model capable of predicting radio emission through the electron cyclotron maser instability using observed stellar magnetic maps of low-mass stars. We determine the structure of the coronal magnetic field and plasma using spectropolarimetric observations of the surface magnetic fields and the X-ray emission measure. We then model the emission of photons from the locations within the corona that satisfy the conditions for electron cyclotron maser emission. Our model predicts the frequency and intensity of radio photons from within the stellar corona. We have benchmarked our model against the low-mass star V374 Peg. This star has both radio observations from the Very Large Array and a nearly simultaneous magnetic map. Using our model we are able to fit the radio observations of V374 Peg, providing additional evidence that the radio emission observed from low-mass stars may originate from the electron cyclotron maser instability. Our model can now be extended to all stars with observed magnetic maps to predict the expected frequency and variability of stellar radio emission in an effort to understand and guide future radio observations of low-mass stars.

Disequilibrium chemical processes significantly affect the spectra of substellar objects. To study these effects, dynamical disequilibrium has been parameterized using the quench and eddy diffusion approximations, but little work has been done to explore how these approximations perform under realistic planetary conditions in different dynamical regimes. As a first step toward addressing this problem, we study the localized, small-scale convective dynamics of planetary atmospheres by direct numerical simulation of fully compressible hydrodynamics with reactive tracers using the Dedalus code. Using polytropically stratified, plane-parallel atmospheres in 2D and 3D, we explore the quenching behavior of different abstract chemical species as a function of the dynamical conditions of the atmosphere as parameterized by the Rayleigh number. We find that in both 2D and 3D, chemical species quench deeper than would be predicted based on simple mixing-length arguments. Instead, it is necessary to employ length scales based on the chemical equilibrium profile of the reacting species in order to predict quench points and perform chemical kinetics modeling in 1D. Based on the results of our simulations, we provide a new length scale, derived from the chemical scale height, that can be used to perform these calculations. This length scale is simple to calculate from known chemical data and makes reasonable predictions for our dynamical simulations.

Spicule oscillations involve high-frequency components with a typical period approximately corresponding to 40–50 s. The typical timescale of the photospheric oscillation is a few minutes, and thus, the origin of this high-frequency component is not trivial. In this study, a one-dimensional numerical simulation is performed to demonstrate that the observed spicule oscillations originate from longitudinal-to-transverse mode conversion that occurs around the equipartition layer in the chromosphere. Calculations are conducted in a self-consistent manner with the exception of additional heating to maintain coronal temperature. The analyses indicate the following features: (1) mode conversion efficiently excites high-frequency transverse waves; (2) the typical period of the high-frequency waves corresponds to the sound-crossing time of the mode conversion region; and (3) simulated root-mean-square velocity of the high-frequency component is consistent with the observed value. These results indicate that the observation of spicule oscillation provides direct evidence of mode conversion in the chromosphere.

We have analyzed the multipolar magnetic field structure variation at neutron star surface by means of the catastrophic eruption model and find that the variation of the geometry of multipolar fields on the magnetar surface could result in the catastrophic rearrangement of the magnetosphere, which provides certain physical mechanism for the outburst of giant flares. The magnetospheric model we adopted consists of two assumptions: (1) a helically twisted flux rope is suspended in an ideal force-free magnetosphere around the magnetar, and (2) a current sheet emerges during the flux rope evolution. Magnetic energy accumulates during the flux rope's gradual evolution along with the variation of magnetar surface magnetic structure before the eruption. The two typical behaviors, either state transition or catastrophic escape, would take place once the flux rope loses equilibrium; thus, tremendous accumulated energy is radiated. We have investigated the equilibrium state of the flux rope and the energy release affected by different multipolar structures and find structures that could trigger violent eruption and provide the radiation approximately 0.5% of the total magnetic energy during the giant flare outburst. Our results provide certain multipolar structures of the neutron star's magnetic field with an energy release percentage 0.42% in the state transition and 0.51% in the catastrophic escape case, which are sufficient for the previously reported energy release from SGR 1806–20 giant flares.

The BL Lac PG 1553+113 has been continuously monitored in gamma-rays with Fermi-LAT for over 9 years. Its updated light curve now includes five iterations of a main pattern comprising a high peak and a longer trough, with a period $P\simeq 2.2\,\mathrm{year}$. Our analysis of 2015–2017 data confirms the occurrence in 2017 January of a new peak fitting in with the previous trend. In addition, we identify secondary peaks ("twin peaks") that occur in closely symmetric pairs on both sides of most main peaks, including the last one; their occurrence is supported by correlated X-ray outbursts. We stress that the above features strongly point to binary dynamics in a system of two black holes (BHs) of some 10⁸ and ${10}^{7}\,{M}_{\odot }$. At periastron the smaller BH periodically stresses the jet j₁ launched by the heavier companion, and triggers MHD–kinetic tearing instabilities. These lead to magnetic reconnections and to acceleration of electrons that produce synchrotron emission from the optical to X-ray bands, and inverse Compton scattering into the GeV range. We discuss two possible origins of the twin peaks : a single-jet model, based on added instabilities induced in j₁ by the smaller companion BH on its inner orbital arc; and a two-jet model with the smaller BH supporting its own, precessing jet j₂ that contributes lower, specific GeV emissions. Such behaviors combining time stability with amplitude variations betray plasma instabilities driven in either jet by binary dynamics, and can provide a double signature of the long-sought supermassive BH binaries.

The star S0-2, which orbits the supermassive black hole (SMBH) in our Galaxy with a period of 16 years, provides the strongest constraint on both the mass of the SMBH and the distance to the Galactic center. S0-2 will soon provide the first measurement of relativistic effects near a SMBH. We report the first limits on the binarity of S0-2 from radial velocity (RV) monitoring, which has implications for both understanding its origin and robustness as a probe of the central gravitational field. With 87 RV measurements, which include 12 new observations that we present, we have the requisite data set to look for RV variations from S0-2's orbital model. Using a Lomb–Scargle analysis and orbit-fitting for potential binaries, we detect no RV variation beyond S0-2's orbital motion and do not find any significant periodic signal. The lack of a binary companion does not currently distinguish different formation scenarios for S0-2. The upper limit on the mass of a companion star (${M}_{\mathrm{comp}}$) still allowed by our results has a median upper limit of ${M}_{\mathrm{comp}}$ sin i ≤ 1.6 M_⊙ for periods between 1 and 150 days, the longest period to avoid tidal break-up of the binary. We also investigate the impact of the remaining allowed binary system on the measurement of the relativistic redshift at S0-2's closest approach in 2018. While binary star systems are important to consider for this experiment, we find that plausible binaries for S0-2 will not alter a 5σ detection of the relativistic redshift.

Massive young stellar objects (MYSOs) in the Magellanic Clouds show infrared absorption features corresponding to significant abundances of CO, CO₂, and H₂O ice along the line of sight, with the relative abundances of these ices differing between the Magellanic Clouds and the Milky Way. CO ice is not detected toward sources in the Small Magellanic Cloud, and upper limits put its relative abundance well below sources in the Large Magellanic Cloud and the Milky Way. We use our gas-grain chemical code MAGICKAL, with multiple grain sizes and grain temperatures, and further expand it with a treatment for increased interstellar radiation field intensity to model the elevated dust temperatures observed in the MCs. We also adjust the elemental abundances used in the chemical models, guided by observations of H ii regions in these metal-poor satellite galaxies. With a grid of models, we are able to reproduce the relative ice fractions observed in MC MYSOs, indicating that metal depletion and elevated grain temperature are important drivers of the MYSO envelope ice composition. Magellanic Cloud elemental abundances have a subgalactic C/O ratio, increasing H₂O ice abundances relative to the other ices; elevated grain temperatures favor CO₂ production over H₂O and CO. The observed shortfall in CO in the Small Magellanic Cloud can be explained by a combination of reduced carbon abundance and increased grain temperatures. The models indicate that a large variation in radiation field strength is required to match the range of observed LMC abundances. CH₃OH abundance is found to be enhanced in low-metallicity models, providing seed material for complex organic molecule formation in the Magellanic Clouds.

Observations of solar and stellar flares have revealed the presence of power-law dependences between the flare energy and the time interval between flares. Various models have been proposed to explain these dependences and the numerical value of the power-law indices. Here, we propose a model in which convective flows in granules force the footpoints of coronal magnetic loops, which are frozen-in to photospheric gas, to undergo a random walk. In certain conditions, this can lead to a twist in the loop, which drives the loop unstable if the twist exceeds a critical value. The possibility that a solar flare is caused by such a twist-induced instability in a loop has been in the literature for decades. Here, we quantify the process in an approximate way with a view to replicating the power-law index. We find that, for relatively small flares, the random walk twisting model leads to a rather steep power-law slope that agrees very well with the index derived from a sample of 56,000+ solar X-ray flares reported by the GOES satellites. For relatively large flares, we find that the slope of the power law is shallower. The empirical power-law slopes reported for flare stars also have a range that overlaps with the slopes obtained here. We suggest that in the coolest stars, a significant change in slope should occur when the frozen-flux assumption breaks down due to low electrical conductivity.

Episodic accretion has been proposed as a solution to the long-standing luminosity problem in star formation; however, the process remains poorly understood. We present observations of line emission from N₂H⁺ and CO isotopologues using the Atacama Large Millimeter/submillimeter Array (ALMA) in the envelopes of eight very low luminosity objects (VeLLOs). In five of the sources the spatial distribution of emission from N₂H⁺ and CO isotopologues shows a clear anticorrelation. It is proposed that this is tracing the CO snow line in the envelopes: N₂H⁺ emission is depleted toward the center of these sources, in contrast to the CO isotopologue emission, which exhibits a peak. The positions of the CO snow lines traced by the N₂H⁺ emission are located at much larger radii than those calculated using the current luminosities of the central sources. This implies that these five sources have experienced a recent accretion burst because the CO snow line would have been pushed outward during the burst because of the increased luminosity of the central star. The N₂H⁺ and CO isotopologue emission from DCE161, one of the other three sources, is most likely tracing a transition disk at a later evolutionary stage. Excluding DCE161, five out of seven sources (i.e., ∼70%) show signatures of a recent accretion burst. This fraction is larger than that of the Class 0/I sources studied by Jørgensen et al. and Frimann et al., suggesting that the interval between accretion episodes in VeLLOs is shorter than that in Class 0/I sources.

The Kennicutt–Schmidt (KS) relationship between the surface density of the star formation rate (SFR) and the gas surface density has three distinct power laws that may result from one model in which gas collapses at a fixed fraction of the dynamical rate. The power-law slope is 1 when the observed gas has a characteristic density for detection, 1.5 for total gas when the thickness is about constant as in the main disks of galaxies, and 2 for total gas when the thickness is regulated by self-gravity and the velocity dispersion is about constant, as in the outer parts of spirals, dwarf irregulars, and giant molecular clouds. The observed scaling of the star formation efficiency (SFR per unit CO) with the dense gas fraction (HCN/CO) is derived from the KS relationship when one tracer (HCN) is on the linear part and the other (CO) is on the 1.5 part. Observations of a threshold density or column density with a constant SFR per unit gas mass above the threshold are proposed to be selection effects, as are observations of star formation in only the dense parts of clouds. The model allows a derivation of all three KS relations using the probability distribution function of density with no thresholds for star formation. Failed galaxies and systems with sub-KS SFRs are predicted to have gas that is dominated by an equilibrium warm phase where the thermal Jeans length exceeds the Toomre length. A squared relation is predicted for molecular gas-dominated young galaxies.

The γ-ray-bright blazar CTA 102 is studied using imaging (new 15 GHz and archival 43 GHz Very Long Baseline Array, VLBA data) and time variable optical flux density, polarization degree, and electric vector position angle (EVPA) spanning between 2015 June 1 and 2016 October 1, covering a prominent γ-ray flare during 2016 January. The pc-scale jet indicates expansion with oscillatory features up to 17 mas. Component proper motions are in the range 0.04–0.33 mas yr⁻¹ with acceleration up to 1.2 mas followed by a slowing down beyond 1.5 mas. A jet bulk Lorentz factor ≥17.5, position angle of 1283, inclination angle ≤66 and intrinsic half opening angle ≤18 are derived from the VLBA data. These inferences are employed in a helical jet model to infer long-term variability in flux density, polarization degree, EVPA, and a rotation of the Stokes Q and U parameters. A core distance of r_{core,43 GHz} = 22.9 pc, and a magnetic field strength at 1 pc and the core location of 1.57 G and 0.07 G, respectively, are inferred using the core-shift method. The study is useful in the context of estimating jet parameters and in offering clues to distinguish mechanisms responsible for variability over different timescales.

The production of the heavy stable proton-rich isotopes between ⁷⁴Se and ¹⁹⁶Hg—the p nuclides—is due to the contribution from different nucleosynthesis processes, activated in different types of stars. Whereas these processes have been subject to various studies, their relative contributions to Galactic chemical evolution (GCE) are still a matter of debate. Here we investigate for the first time the nucleosynthesis of p nuclides in GCE by including metallicity and progenitor mass-dependent yields of core-collapse supernovae (ccSNe) into a chemical evolution model. We used a grid of metallicities and progenitor masses from two different sets of stellar yields and followed the contribution of ccSNe to the Galactic abundances as a function of time. In combination with previous studies on p-nucleus production in thermonuclear supernovae (SNIa), and using the same GCE description, this allows us to compare the respective roles of SNeIa and ccSNe in the production of p-nuclei in the Galaxy. The γ process in ccSN is very efficient for a wide range of progenitor masses (13 M_⊙–25 M_⊙) at solar metallicity. Since it is a secondary process with its efficiency depending on the initial abundance of heavy elements, its contribution is strongly reduced below solar metallicity. This makes it challenging to explain the inventory of the p nuclides in the solar system by the contribution from ccSNe alone. In particular, we find that ccSNe contribute less than 10% of the solar p nuclide abundances, with only a few exceptions. Due to the uncertain contribution from other nucleosynthesis sites in ccSNe, such as neutrino winds or α-rich freeze out, we conclude that the light p-nuclides ⁷⁴Se, ⁷⁸Kr, ⁸⁴Sr, and ⁹²Mo may either still be completely or only partially produced in ccSNe. The γ-process accounts for up to twice the relative solar abundances for ⁷⁴Se in one set of stellar models and ¹⁹⁶Hg in the other set. The solar abundance of the heaviest p nucleus ¹⁹⁶Hg is reproduced within uncertainties in one set of our models due to photodisintegration of the Pb isotopes ^208,207,206Pb. For all other p nuclides, abundances as low as 2% of the solar level were obtained.

Future observations of terrestrial exoplanet atmospheres will occur for planets at different stages of geological evolution. We expect to observe a wide variety of atmospheres and planets with alternative evolutionary paths, with some planets resembling Earth at different epochs. For an Earth-like atmospheric time trajectory, we simulate planets from the prebiotic to the current atmosphere based on geological data. We use a stellar grid F0V to M8V (${T}_{\mathrm{eff}}=7000$–2400 K) to model four geological epochs of Earth's history corresponding to a prebiotic world (3.9 Ga), the rise of oxygen at 2.0 Ga and at 0.8 Ga, and the modern Earth. We show the VIS–IR spectral features, with a focus on biosignatures through geological time for this grid of Sun-like host stars and the effect of clouds on their spectra. We find that the observability of biosignature gases reduces with increasing cloud cover and increases with planetary age. The observability of the visible O₂ feature for lower concentrations will partly depend on clouds, which, while slightly reducing the feature, increase the overall reflectivity, and thus the detectable flux of a planet. The depth of the IR ozone feature contributes substantially to the opacity at lower oxygen concentrations, especially for the high near-UV stellar environments around F stars. Our results are a grid of model spectra for atmospheres representative of Earth's geological history to inform future observations and instrument design and are available online at http://carlsaganinstitute.org/data/.

We present in situ observations of magnetic turbulence in the draped interstellar magnetic field ${\boldsymbol{B}}$ measured by Voyager 1 during an undisturbed interval from 2015.3987 to 2016.6759 confirming the existence of the turbulence observed previously from 2013.3593 to 2014.6373. The power spectral density of the turbulence was the same in both cases. The turbulence had a Kolmogorov k^−5/3 spectrum in the range from k = 1.3 × 10⁻¹³ cm⁻¹ to 4 × 10⁻¹² cm⁻¹. The ratio of the turbulent fluctuations to the average magnetic field strength was only 0.02, indicating that the turbulence was very weak. Extrapolating the power-law slope to lower frequencies yields an upper limit on the turbulence outer scale of 0.01 pc = 2000 au, which may be regarded as the distance at which Voyager 1 will enter the undisturbed local interstellar medium, beyond the outer heliosheath or bow wave in the upstream direction. The maximum variance of the fluctuations was in the two directions transverse to the average magnetic field in the recent interval, whereas it was parallel to the average magnetic field in the earlier interval, suggesting a transformation from turbulence with a dominant compressive component to turbulence dominated by transverse fluctuations. As the magnitude of the fluctuations was approaching that of the uncertainties of the measurements, the latter result requires confirmation by further observations.

Low-mass low-density planets discovered by Kepler in the super-Earth mass regime typically have large radii for their inferred masses, implying the presence of H₂–He atmospheres. These planets are vulnerable to atmospheric mass loss due to heating by the parent star's XUV flux. Models coupling atmospheric mass loss with thermal evolution predicted a bimodal distribution of planetary radii, which has gained observational support. However, a key component that has been ignored in previous studies is the dissolution of these gases into the molten core of rock and iron that constitute most of their mass. Such planets have high temperatures (>2000 K) and pressures (∼kbars) at the core-envelope boundary, ensuring a molten surface and a subsurface reservoir of hydrogen that can be 5–10 times larger than the atmosphere. This study bridges this gap by coupling the thermal evolution of the planet and the mass loss of the atmosphere with the thermodynamic equilibrium between the dissolved H₂ and the atmospheric H₂ (Henry's law). Dissolution in the interior allows a planet to build a larger hydrogen repository during the planet formation stage. We show that the dissolved hydrogen outgasses to buffer atmospheric mass loss. The slow cooling of the planet also leads to outgassing because solubility decreases with decreasing temperature. Dissolution of hydrogen in the interior therefore increases the atmosphere retention ability of super-Earths. The study highlights the importance of including the temperature- and pressure-dependent solubility of gases in magma oceans and coupling outgassing to planetary evolution models.

The stellar disk size of a galaxy depends on the ratio of the disk stellar mass to the halo mass, m_⋆ ≡ M_⋆/M_dh, and the fraction of the dark halo angular momentum transferred to the stellar disk, j_⋆ ≡ J_⋆/J_dh. Since m_⋆ and j_⋆ are determined by many star formation–related processes, measuring j_⋆ and m_⋆ at various redshifts is essential to understanding the formation history of disk galaxies. We use the 3D-HST GOODS-S, COSMOS, and AEGIS imaging data and photo-z catalog to examine j_⋆ and m_⋆ for star-forming galaxies at z ∼ 2–4, when disks are actively forming. We find that the j_⋆/m_⋆ ratio is ≃0.77 ± 0.06 for all three redshifts over the entire mass range examined, 8 × 10¹⁰ < M_dh/h⁻¹M_⊙ < 2 × 10¹², with a possible (<30%) decrease with mass. This high ratio is close to those of local disk galaxies, descendants of our galaxies in terms of M_dh growth, implying a nearly constant j_⋆/m_⋆ over the past 12 Gyr. These results are remarkable because mechanisms controlling angular momentum transfer to disks such as inflows and feedback depend on both cosmic time and halo mass, and, indeed, theoretical studies tend to predict j_⋆/m_⋆ changing with redshift and mass. It is found that some of the recent theoretical galaxy formation simulations predict a smaller j_⋆/m_⋆ than our values. We also find that a significant fraction of our galaxies appears to be unstable against bar formation.

Understanding the transport of solar energetic particles (SEPs) from acceleration sites at the Sun into interplanetary space and to the Earth is an important question for forecasting space weather. The interplanetary magnetic field (IMF), with two distinct polarities and a complex structure, governs energetic particle transport and drifts. We analyze for the first time the effect of a wavy heliospheric current sheet (HCS) on the propagation of SEPs. We inject protons close to the Sun and propagate them by integrating fully 3D trajectories within the inner heliosphere in the presence of weak scattering. We model the HCS position using fits based on neutral lines of magnetic field source surface maps (SSMs). We map 1 au proton crossings, which show efficient transport in longitude via HCS, depending on the location of the injection region with respect to the HCS. For HCS tilt angles around 30°–40°, we find significant qualitative differences between A+ and A− configurations of the IMF, with stronger fluences along the HCS in the former case but with a distribution of particles across a wider range of longitudes and latitudes in the latter. We show how a wavy current sheet leads to longitudinally periodic enhancements in particle fluence. We show that for an A+ IMF configuration, a wavy HCS allows for more proton deceleration than a flat HCS. We find that A− IMF configurations result in larger average fluences than A+ IMF configurations, due to a radial drift component at the current sheet.

It is established that there is a dependence of the luminosity of type Ia supernovae (SNe Ia) on environment: SNe Ia in young, star-forming, metal-poor stellar populations appear fainter after light-curve shape corrections than those in older, passive, metal-rich environments. This is accounted for in cosmological studies using a global property of the SN host galaxy, typically the host galaxy stellar mass. However, recent low-redshift studies suggest that this effect manifests itself most strongly when using the local star formation rate (SFR) at the SN location, rather than the global SFR or the stellar mass of the host galaxy. At high-redshift, such local SFRs are difficult to determine; here, we show that an equivalent local correction can be made by restricting the SN Ia sample in globally star-forming host galaxies to a low-mass host galaxy subset (≤10¹⁰M_☉). Comparing this sample of SNe Ia (in locally star-forming environments) to those in locally passive host galaxies, we find that SNe Ia in locally star-forming environments are 0.081 ± 0.018 mag fainter (4.5σ), consistent with the result reported by Rigault et al., but our conclusion is based on a sample ∼5 times larger over a wider redshift range. This is a larger difference than when splitting the SN Ia sample based on global host galaxy SFR or host galaxy stellar mass. This method can be used in ongoing and future high-redshift SN surveys, where local SN Ia environments are difficult to determine.

${\mathrm{OH}}^{+}$ and ${{\rm{H}}}_{2}{{\rm{O}}}^{+}$ cations play a significant role in the chemistry of the cold interstellar medium and hence their hydrogen abstraction reactions with ${{\rm{H}}}_{2}$ have to be included in ion chemical models. The reactions lead directly or indirectly to ${{\rm{H}}}_{3}{{\rm{O}}}^{+}$ ions that subsequently recombine with electrons and dissociate into H atoms and ${{\rm{H}}}_{2}{\rm{O}}$. The experiments described in this paper provide rate coefficients (${k}_{{\mathrm{OH}}^{+}}$ and ${k}_{{{\rm{H}}}_{2}{{\rm{O}}}^{+}}$) for the reactions of ${\mathrm{OH}}^{+}$ and ${{\rm{H}}}_{2}{{\rm{O}}}^{+}$ with ${{\rm{H}}}_{2}$ over a wide temperature range (from 15 to 300 K). A cryogenic 22-pole RF ion trap instrument was employed for this purpose. It was found that ${k}_{{\mathrm{OH}}^{+}}$ increases from $(0.76\pm 0.30)\times {10}^{-9}\,{\mathrm{cm}}^{3}\,{{\rm{s}}}^{-1}$ at 17 K to $(1.24\pm 0.25)\times {10}^{-9}\,{\mathrm{cm}}^{3}\,{{\rm{s}}}^{-1}$ at 263 K while ${k}_{{{\rm{H}}}_{2}{{\rm{O}}}^{+}}$ is nearly constant, varying from $(0.93\pm 0.35)\,\times {10}^{-9}\,{\mathrm{cm}}^{3}\,{{\rm{s}}}^{-1}$ at 17 K to $(1.00\pm 0.25)\times {10}^{-9}\,{\mathrm{cm}}^{3}\,{{\rm{s}}}^{-1}$ at 218 K.

We investigate the relationship between the blueshifts of a hot emission line and the nonthermal emissions in microwave and hard X-ray (HXR) wavelengths in the precursor of a solar flare on 2014 October 27. The flare precursor is identified as a small but well-developed peak in the soft X-ray and extreme-ultraviolet passbands before the GOES flare onset, which is accompanied by a pronounced burst in microwave 17 and 34 GHz and in HXR 25–50 keV. The slit of the Interface Region Imaging Spectrograph (IRIS) stays on one ribbon-like transient during the flare precursor phase, which shows visible nonthermal emissions in Nobeyama Radioheliograph and RHESSI images. The IRIS spectroscopic observations show that the hot line of Fe xxi 1354.09 Å (log T ∼ 7.05) displays blueshifts, while the cool line of Si iv 1402.77 Å (log T ∼ 4.8) exhibits redshifts. The blueshifts and redshifts are well correlated with each other, indicative of an explosive chromospheric evaporation during the flare precursor phase combining a high nonthermal energy flux with a short characteristic timescale. In addition, the blueshifts of Fe xxi 1354.09 Å are well correlated with the microwave and HXR emissions, implying that the explosive chromospheric evaporation during the flare precursor phase is driven by nonthermal electrons.

The first gas-phase infrared spectra of two isolated astronomically relevant and large polycyclic aromatic hydrocarbon (PAH) cations—diindenoperylene (DIP) and dicoronylene (DC)—in the 530–1800 cm⁻¹ (18.9−5.6 μm) range—are presented. Vibrational band positions are determined for comparison to the aromatic infrared bands. The spectra are obtained via infrared multiphoton dissociation spectroscopy of ions stored in a quadrupole ion trap using the intense and tunable radiation of the free electron laser for infrared experiments (FELIX). DIP⁺ shows its main absorption peaks at 737 (13.57), 800 (12.50), 1001 (9.99), 1070 (9.35), 1115 (8.97), 1152 (8.68), 1278 (7.83), 1420 (7.04), and 1550 (6.45) cm⁻¹(μm), in good agreement with density functional theory (DFT) calculations that are uniformly scaled to take anharmonicities into account. DC⁺ has its main absorption peaks at 853 (11.72), 876 (11.42), 1032 (9.69), 1168 (8.56), 1300 (7.69), 1427 (7.01), and 1566 (6.39) cm⁻¹(μm), which also agree well with the scaled DFT results presented here. The DIP⁺ and DC⁺ spectra are compared with the prominent infrared features observed toward NGC 7023. This results both in matches and clear deviations. Moreover, in the 11.0–14.0 μm region, specific bands can be linked to CH out-of-plane (oop) bending modes of different CH edge structures in large PAHs. The molecular origin of these findings and their astronomical relevance are discussed.

The analysis of a series of seven observations of the nearby (z = 0.0809) QSO PG 1211+143, taken with the Reflection Grating Spectrometer (RGS) onboard XMM-Newton in 2014, are presented. The high-resolution soft X-ray spectrum, with a total exposure exceeding 600 ks, shows a series of blueshifted absorption lines from the He and H-like transitions of N, O, and Ne, as well as from L-shell Fe. The strongest absorption lines are all systematically blueshifted by −0.06c, originating in two absorption zones from low- and high-ionization gas. Both zones are variable on timescales of days, with the variations in absorber opacity effectively explained by either column density changes or the absorber ionization responding directly to the continuum flux. We find that the soft X-ray absorbers probably exist in a two-phase wind at a radial distance of ∼10¹⁷–10¹⁸ cm from the black hole with the lower-ionization gas as denser clumps embedded within a higher-ionization outflow. The overall mass outflow rate of the soft X-ray wind may be as high as $2{M}_{\odot }$ yr⁻¹, close to the Eddington rate for PG 1211+143 and similar to that previously deduced from the Fe K absorption.

The multiplexing capability of slitless spectroscopy is a powerful asset in creating large spectroscopic data sets, but issues such as spectral confusion make the interpretation of the data challenging. Here we present a new method to search for emission lines in the slitless spectroscopic data from the 3D-HST survey utilizing the Wide-Field Camera 3 on board the Hubble Space Telescope. Using a novel statistical technique, we can detect compact (extended) emission lines at 90% completeness down to fluxes of $1.5(3.0)\times {10}^{-17}\,\mathrm{erg}\,{{\rm{s}}}^{-1}\,{\mathrm{cm}}^{-2}$, close to the noise level of the grism exposures, for objects detected in the deep ancillary photometric data. Unlike previous methods, the Bayesian nature allows for probabilistic line identifications, namely redshift estimates, based on secondary emission line detections and/or photometric redshift priors. As a first application, we measure the comoving number density of Extreme Emission Line Galaxies (restframe [O iii] λ5007 equivalent widths in excess of 500 Å). We find that these galaxies are nearly 10× more common above z ∼ 1.5 than at z ≲ 0.5. With upcoming large grism surveys such as Euclid and WFIRST, as well as grisms featured prominently on the NIRISS and NIRCam instruments on the James Webb Space Telescope, methods like the one presented here will be crucial for constructing emission line redshift catalogs in an automated and well-understood manner.

We study the effects of galaxy environment on the evolution of the stellar mass function (SMF) over 0.2 < z < 2.0 using the FourStar Galaxy Evolution (ZFOURGE) Survey and NEWFIRM Medium-Band Survey (NMBS) down to the stellar mass completeness limit, $\mathrm{log}{M}_{* }/{M}_{\odot }\gt 9.0$ (9.5) at z = 1.0 (2.0). We compare the SMFs for quiescent and star-forming galaxies in the highest and lowest environments using a density estimator based on the distance to the galaxies' third-nearest neighbors. For star-forming galaxies, at all redshifts there are only minor differences with environment in the shape of the SMF. For quiescent galaxies, the SMF in the lowest densities shows no evolution with redshift other than an overall increase in number density (ϕ*) with time. This suggests that the stellar mass dependence of quenching in relatively isolated galaxies both is universal and does not evolve strongly. While at $z\gtrsim 1.5$, the SMF of quiescent galaxies is indistinguishable in the highest and lowest densities, at lower redshifts, it shows a rapidly increasing number density of lower-mass galaxies, $\mathrm{log}{M}_{* }/{M}_{\odot }\simeq 9\mbox{--}10$, in the highest-density environments. We argue that this evolution can account for all the redshift evolution in the shape of the total quiescent galaxy SMF. This evolution in the quiescent galaxy SMF at higher redshift (z > 1) requires an environmental quenching efficiency that decreases with decreasing stellar mass at 0.5 < z < 1.5 or it would overproduce the number of lower-mass quiescent galaxies in denser environments. This requires a dominant environmental process such as starvation combined with rapid gas depletion and ejection at z > 0.5–1.0 for galaxies in our mass range. The efficiency of this process decreases with redshift, allowing other processes (such as galaxy interactions and ram-pressure stripping) to become more important at later times, z < 0.5.

We analyze results from the first 18 months of monthly submillimeter monitoring of eight star-forming regions in the JCMT Transient Survey. In our search for stochastic variability in 1643 bright peaks, only the previously identified source, EC 53, shows behavior well above the expected measurement uncertainty. Another four sources—two disks and two protostars—show moderately enhanced standard deviations in brightness, as expected for stochastic variables. For the two protostars, this apparent variability is the result of single epochs that are much brighter than the mean. In our search for secular brightness variations that are linear in time, we measure the fractional brightness change per year for 150 bright peaks, 50 of which are protostellar. The ensemble distribution of slopes is well fit by a normal distribution with σ ∼ 0.023. Most sources are not rapidly brightening or fading at submillimeter wavelengths. Comparison against time-randomized realizations shows that the width of the distribution is dominated by the uncertainty in the individual brightness measurements of the sources. A toy model for secular variability reveals that an underlying Gaussian distribution of linear fractional brightness change σ = 0.005 would be unobservable in the present sample, whereas an underlying distribution with σ = 0.02 is ruled out. Five protostellar sources, 10% of the protostellar sample, are found to have robust secular measures deviating from a constant flux. The sensitivity to secular brightness variations will improve significantly with a sample over a longer time duration, with an improvement by factor of two expected by the conclusion of our 36 month survey.

A new model describing the transport and evolution of turbulence in the quiet solar corona is presented. In the low plasma beta environment, transverse photospheric convective fluid motions drive predominantly quasi-2D (nonpropagating) turbulence in the mixed-polarity "magnetic carpet," together with a minority slab (Alfvénic) component. We use a simplified sub-Alfvénic flow velocity profile to solve transport equations describing the evolution and dissipation of turbulence from $1\hspace{0.5em}{\rm{t}}{\rm{o}}\,15\,{R}_{\odot }$ (including the Alfvén surface). Typical coronal base parameters are used, although one model uses correlation lengths derived observationally by Abramenko et al., and the other assumes values 10 times larger. The model predicts that (1) the majority quasi-2D turbulence evolves from a balanced state at the coronal base to an imbalanced state, with outward fluctuations dominating, at and beyond the Alfvén surface, i.e., inward turbulent fluctuations are dissipated preferentially; (2) the initially imbalanced slab component remains imbalanced throughout the solar corona, being dominated by outwardly propagating Alfvén waves, and wave reflection is weak; (3) quasi-2D turbulence becomes increasingly magnetized, and beyond $\sim 6\,{R}_{\odot }$, the kinetic energy is mainly in slab fluctuations; (4) there is no accumulation of inward energy at the Alfvén surface; (5) inertial range quasi-2D rather than slab fluctuations are preferentially dissipated within $\sim 3\,{R}_{\odot }$; and (6) turbulent dissipation of quasi-2D fluctuations is sufficient to heat the corona to temperatures $\sim 2\times {10}^{6}$ K within $2\,{R}_{\odot }$, consistent with observations that suggest that the fast solar wind is accelerated most efficiently between $\sim 2\hspace{0.5em}{\rm{a}}{\rm{n}}{\rm{d}}\,4\,{R}_{\odot }$.

We discuss the spectral analysis of a sample of 63 active galactic nuclei (AGN) detected above a limiting flux of $S(8\mbox{--}24\,\mathrm{keV})=7\times {10}^{-14}$$\mathrm{erg}\,{{\rm{s}}}^{-1}\,{\mathrm{cm}}^{-2}$ in the multi-tiered NuSTAR extragalactic survey program. The sources span a redshift range $z=0\mbox{--}2.1$ (median $\langle z\rangle =0.58$). The spectral analysis is performed over the broad 0.5–24 keV energy range, combining NuSTAR with Chandra and/or XMM-Newton data and employing empirical and physically motivated models. This constitutes the largest sample of AGN selected at $\gt 10\,\mathrm{keV}$ to be homogeneously spectrally analyzed at these flux levels. We study the distribution of spectral parameters such as photon index, column density (${N}_{{\rm{H}}}$), reflection parameter (${\boldsymbol{R}}$), and 10–40 keV luminosity (${L}_{{\rm{X}}}$). Heavily obscured ($\mathrm{log}[{N}_{{\rm{H}}}/{\mathrm{cm}}^{-2}]\geqslant 23$) and Compton-thick (CT; $\mathrm{log}[{N}_{{\rm{H}}}/{\mathrm{cm}}^{-2}]\geqslant 24$) AGN constitute ∼25% (15–17 sources) and ∼2–3% (1–2 sources) of the sample, respectively. The observed ${N}_{{\rm{H}}}$ distribution agrees fairly well with predictions of cosmic X-ray background population-synthesis models (CXBPSM). We estimate the intrinsic fraction of AGN as a function of ${N}_{{\rm{H}}}$, accounting for the bias against obscured AGN in a flux-selected sample. The fraction of CT AGN relative to $\mathrm{log}[{N}_{{\rm{H}}}/{\mathrm{cm}}^{-2}]=20\mbox{--}24$ AGN is poorly constrained, formally in the range 2–56% (90% upper limit of 66%). We derived a fraction (f_abs) of obscured AGN ($\mathrm{log}[{N}_{{\rm{H}}}/{\mathrm{cm}}^{-2}]=22\mbox{--}24$) as a function of ${L}_{{\rm{X}}}$ in agreement with CXBPSM and previous $z\lt 1$ X-ray determinations. Furthermore, f_abs at $z=0.1\mbox{--}0.5$ and $\mathrm{log}({L}_{{\rm{x}}}/\mathrm{erg}\,{{\rm{s}}}^{-1})\approx 43.6\mbox{--}44.3$ agrees with observational measurements/trends obtained over larger redshift intervals. We report a significant anti-correlation of R with ${L}_{{\rm{X}}}$ (confirmed by our companion paper on stacked spectra) with considerable scatter around the median R values.

By numerically solving the equations of rotating magnetohydrodynamics, we study the magnetic effect on dynamical tide. We find that a magnetic field has a significant impact not only on the flow structure, i.e., the internal shear layers in a rotating flow can be destroyed in the presence of a moderate or stronger magnetic field (in the sense that the Alfvén velocity is at least of the order of 0.1 of the surface rotational velocity), but also on the dispersion relation of waves excited by tidal force such that the range of tidal resonance is broadened by a magnetic field. A major result is that the total tidal dissipation scales as a square of the field strength, which can be used to estimate the strength of the internal magnetic field in the astronomical object of a binary system. Moreover, with a moderate or stronger field, the ratio of magnetic dissipation to viscous dissipation is almost inversely proportional to the magnetic Prandtl number (i.e., the ratio of viscosity to magnetic diffusivity); thus, in the astrophysical situation at a small magnetic Prandtl number magnetic dissipation dominates over viscous dissipation with a moderate or stronger field.

We address the turbulent fragmentation scenario for the origin of the stellar initial mass function (IMF), using a large set of numerical simulations of randomly driven supersonic MHD turbulence. The turbulent fragmentation model successfully predicts the main features of the observed stellar IMF assuming an isothermal equation of state without any stellar feedback. As a test of the model, we focus on the case of a magnetized isothermal gas, neglecting stellar feedback, while pursuing a large dynamic range in both space and timescales covering the full spectrum of stellar masses from brown dwarfs to massive stars. Our simulations represent a generic 4 pc region within a typical Galactic molecular cloud, with a mass of 3000 M_⊙ and an rms velocity 10 times the isothermal sound speed and 5 times the average Alfvén velocity, in agreement with observations. We achieve a maximum resolution of 50 au and a maximum duration of star formation of 4.0 Myr, forming up to a thousand sink particles whose mass distribution closely matches the observed stellar IMF. A large set of medium-size simulations is used to test the sink particle algorithm, while larger simulations are used to test the numerical convergence of the IMF and the dependence of the IMF turnover on physical parameters predicted by the turbulent fragmentation model. We find a clear trend toward numerical convergence and strong support for the model predictions, including the initial time evolution of the IMF. We conclude that the physics of isothermal MHD turbulence is sufficient to explain the origin of the IMF.

We investigate the evolution of dust content in galaxies from redshifts z = 0 to z = 9.5. Using empirically motivated prescriptions, we model galactic-scale properties—including halo mass, stellar mass, star formation rate, gas mass, and metallicity—to make predictions for the galactic evolution of dust mass and dust temperature in main-sequence galaxies. Our simple analytic model, which predicts that galaxies in the early universe had greater quantities of dust than their low-redshift counterparts, does a good job of reproducing observed trends between galaxy dust and stellar mass out to z ≈ 6. We find that for fixed galaxy stellar mass, the dust temperature increases from z = 0 to z = 6. Our model forecasts a population of low-mass, high-redshift galaxies with interstellar dust as hot as, or hotter than, their more massive counterparts; but this prediction needs to be constrained by observations. Finally, we make predictions for observing 1.1 mm flux density arising from interstellar dust emission with the Atacama Large Millimeter Array.

We present observations of DES16C2nm, the first spectroscopically confirmed hydrogen-free superluminous supernova (SLSN-I) at redshift $z\approx 2$. DES16C2nm was discovered by the Dark Energy Survey (DES) Supernova Program, with follow-up photometric data from the Hubble Space Telescope, Gemini, and the European Southern Observatory Very Large Telescope supplementing the DES data. Spectroscopic observations confirm DES16C2nm to be at z = 1.998, and spectroscopically similar to Gaia16apd (a SLSN-I at z = 0.102), with a peak absolute magnitude of $U=-22.26\pm 0.06$. The high redshift of DES16C2nm provides a unique opportunity to study the ultraviolet (UV) properties of SLSNe-I. Combining DES16C2nm with 10 similar events from the literature, we show that there exists a homogeneous class of SLSNe-I in the UV (${\lambda }_{\mathrm{rest}}\approx 2500$ Å), with peak luminosities in the (rest-frame) U band, and increasing absorption to shorter wavelengths. There is no evidence that the mean photometric and spectroscopic properties of SLSNe-I differ between low ($z\lt 1$) and high redshift ($z\gt 1$), but there is clear evidence of diversity in the spectrum at ${\lambda }_{\mathrm{rest}}\lt 2000\,\mathring{\rm A}$, possibly caused by the variations in temperature between events. No significant correlations are observed between spectral line velocities and photometric luminosity. Using these data, we estimate that SLSNe-I can be discovered to z = 3.8 by DES. While SLSNe-I are typically identified from their blue observed colors at low redshift ($z\lt 1$), we highlight that at $z\gt 2$ these events appear optically red, peaking in the observer-frame z-band. Such characteristics are critical to identify these objects with future facilities such as the Large Synoptic Survey Telescope, Euclid, and the Wide-field Infrared Survey Telescope, which should detect such SLSNe-I to z = 3.5, 3.7, and 6.6, respectively.

We study the stellar bar growth in high-resolution numerical galaxy models with and without dark matter halos. In all models, the galactic disk is exponential, and the halos are rigid or live Plummer spheres. More specifically, when there is no dark matter halo, we modify the gravitational force between point particles. To do so, we use the weak field limit of an alternative theory of dark matter known as MOG in the literature. The galaxy model in MOG has the same initial conditions as galaxy models with a dark matter halo. On the other hand, the initial random velocities and Toomre's local stability parameter are the same for all of the models. We show that the evolution and growth of the bar in MOG is substantially different from the standard cases including dark matter halo. More importantly, we find that the bar growth rate and its final magnitude are smaller in MOG. On the other hand, the maximum value of the bar in MOG is smaller than that in the Newtonian models. It is shown that although the live dark matter halo may support bar instability, MOG has stabilizing effects. Furthermore, we show that MOG supports fast pattern speeds, and unlike in the dark matter halo models, the pattern speed does not decrease with time. These differences, combined with the relevant observations, may help to distinguish between dark matter and modified gravity in galactic scales.

We present deep spectroscopic observations of a Lyman break galaxy (LBG) candidate (hereafter MACS1149-JD) at z ∼ 9.5 with the Hubble Space Telescope (HST) WFC3/IR grisms. The grism observations were taken at four distinct position angles, totaling 34 orbits with the G141 grism, although only 19 of the orbits are relatively uncontaminated along the trace of MACS1149-JD. We fit a three-parameter (z, F160W mag, and Lyα equivalent width [EW]) LBG template to the three least contaminated grism position angles using a Markov chain Monte Carlo approach. The grism data alone are best fit with a redshift of ${z}_{\mathrm{grism}}={9.53}_{-0.60}^{+0.39}$ (68% confidence), in good agreement with our photometric estimate of ${z}_{\mathrm{phot}}={9.51}_{-0.12}^{+0.06}$ (68% confidence). Our analysis rules out Lyα emission from MACS1149-JD above a 3σ EW of 21 Å, consistent with a highly neutral IGM. We explore a scenario where the red Spitzer/IRAC [3.6]–[4.5] color of the galaxy previously pointed out in the literature is due to strong rest-frame optical emission lines from a very young stellar population rather than a 4000 Å break. We find that while this can provide an explanation for the observed IRAC color, it requires a lower redshift (z ≲ 9.1), which is less preferred by the HST imaging data. The grism data are consistent with both scenarios, indicating that the red IRAC color can still be explained by a 4000 Å break, characteristic of a relatively evolved stellar population. In this interpretation, the photometry indicates that a ${340}_{-35}^{+29}$ Myr stellar population is already present in this galaxy only ∼500 Myr after the big bang.

Observations of heavy metal pollution in white dwarf stars indicate that metal-rich planetesimals are frequently scattered into star-grazing orbits, tidally disrupted, and accreted onto the white dwarf surface, offering direct insight into the dynamical evolution of post-main-sequence exoplanetary systems. Emission lines from the gaseous debris in the accretion disks of some of these systems show variations on timescales of decades, and have been interpreted as the general relativistic precession of a recently formed, elliptical disk. Here we present a comprehensive spectroscopic monitoring campaign of the calcium infrared triplet emission in one system, HE 1349–2305, which shows morphological emission profile variations suggestive of a precessing, asymmetric intensity pattern. The emission profiles are shown to vary on a timescale of one to two years, which is an order of magnitude shorter than what has been observed in other similar systems. We demonstrate that this timescale is likely incompatible with general relativistic precession, and consider alternative explanations for the rapid evolution, including the propagation of density waves within the gaseous debris. We conclude with recommendations for follow-up observations, and discuss how the rapid evolution of the gaseous debris in HE 1349–2305 could be leveraged to test theories of exoplanetary debris disk evolution around white dwarf stars.

It has been proposed that primordial black holes (PBHs) form binaries in the radiation dominated era. Once formed, some fraction of them may merge within the age of the universe by gravitational radiation reaction. We investigate the merger rate of the PBH binaries when the PBHs have a distribution of masses around ${ \mathcal O }(10){M}_{\odot }$, which is a generalization of the previous studies where the PBHs are assumed to have the same mass. After deriving a formula for the merger time probability distribution in the PBH mass plane, we evaluate it under two different approximations. We identify a quantity constructed from the mass distribution of the merger rate density per unit cosmic time and comoving volume ${ \mathcal R }({m}_{1},{m}_{2})$, $\alpha =-{({m}_{1}+{m}_{2})}^{2}{\partial }^{2}\mathrm{ln}{ \mathcal R }/\partial {m}_{1}\partial {m}_{2}$, which universally satisfies 0.97 ≲ α ≲ 1.05 for all binary masses independently of the PBH mass function. This result suggests that the measurement of this quantity is useful for testing the PBH scenario.

The basic unified model of active galactic nuclei (AGNs) invokes an anisotropic obscuring structure, usually referred to as a torus, to explain AGN obscuration as an angle-dependent effect. We present a new grid of X-ray spectral templates based on radiative transfer calculations in neutral gas in an approximately toroidal geometry, appropriate for CCD-resolution X-ray spectra (FWHM ≥ 130 eV). Fitting the templates to broadband X-ray spectra of AGNs provides constraints on two important geometrical parameters of the gas distribution around the supermassive black hole: the average column density and the covering factor. Compared to the currently available spectral templates, our model is more flexible, and capable of providing constraints on the main torus parameters in a wider range of AGNs. We demonstrate the application of this model using hard X-ray spectra from NuSTAR (3–79 keV) for four AGNs covering a variety of classifications: 3C 390.3, NGC 2110, IC 5063, and NGC 7582. This small set of examples was chosen to illustrate the range of possible torus configurations, from disk-like to sphere-like geometries with column densities below, as well as above, the Compton-thick threshold. This diversity of torus properties challenges the simple assumption of a standard geometrically and optically thick toroidal structure commonly invoked in the basic form of the unified model of AGNs. Finding broad consistency between our constraints and those from infrared modeling, we discuss how the approach from the X-ray band complements similar measurements of AGN structures at other wavelengths.

We construct a simple but self-consistent collapsar model for gamma-ray bursts (GRBs) and SNe associated with GRBs (GRB-SNe). Our model includes a black hole, an accretion disk, and the envelope surrounding the central system. The evolutions of the different components are connected by the transfer of the mass and angular momentum. To address properties of the jet and the wind-driven SNe, we consider competition of the ram pressure from the infalling envelope and those from the jet and wind. The expected properties of the GRB jet and the wind-driven SN are investigated as a function of the progenitor mass and angular momentum. We find two conditions that should be satisfied if the wind-driven explosion is to explain the properties of the observed GRB-SNe: (1) the wind should be collimated at its base, and (2) it should not prevent further accretion even after the launch of the SN explosion. Under these conditions, some relations seen in the properties of the GRB-SNe could be reproduced by a sequence of different angular momentum in the progenitors. Only the model with the largest angular momentum could explain the observed (energetic) GRB-SNe, and we expect that the collapsar model can result in a wide variety of observational counterparts, mainly depending on the angular momentum of the progenitor star.

Solar-type binaries with short orbital periods (${P}_{\mathrm{close}}\equiv 1\mbox{--}10$ days; a ≲ 0.1 au) cannot form directly via fragmentation of molecular clouds or protostellar disks, yet their component masses are highly correlated, suggesting interaction during the pre-main-sequence (pre-MS) phase. Moreover, the close binary fraction of pre-MS stars is consistent with that of their MS counterparts in the field (${F}_{\mathrm{close}}=2.1 \%$). Thus, we can infer that some migration mechanism operates during the early pre-MS phase (τ ≲ 5 Myr) that reshapes the primordial separation distribution. We test the feasibility of this hypothesis by carrying out a population synthesis calculation which accounts for two formation channels: Kozai–Lidov (KL) oscillations and dynamical instability in triple systems. Our models incorporate (1) more realistic initial conditions compared to previous studies, (2) octupole-level effects in the secular evolution, (3) tidal energy dissipation via weak-friction equilibrium tides at small eccentricities and via non-radial dynamical oscillations at large eccentricities, and (4) the larger tidal radius of a pre-MS primary. Given a 15% triple-star fraction, we simulate a close binary fraction from KL oscillations alone of ${F}_{\mathrm{close}}\approx 0.4 \%$ after τ = 5 Myr, which increases to ${F}_{\mathrm{close}}\approx 0.8 \%$ by τ = 5 Gyr. Dynamical ejections and disruptions of unstable coplanar triples in the disk produce solitary binaries with slightly longer periods P ≈ 10–100 days. The remaining ≈60% of close binaries with outer tertiaries, particularly those in compact coplanar configurations with log ${P}_{\mathrm{out}}$ (days) ≈ 2–5 (${a}_{\mathrm{out}}\lt 50$ au), can be explained only with substantial extra energy dissipation due to interactions with primordial gas.

In this paper, we report a new estimate of the absolute proper motion (PM) of the globular cluster NGC 5139 (ω Cen) as part of the HST large program GO-14118+14662. We analyzed a field 17 arcmin southwest of the center of ω Cen and computed PMs with epoch spans of ∼15.1 years. We employed 45 background galaxies to link our relative PMs to an absolute reference-frame system. The absolute PM of the cluster in our field is $({\mu }_{\alpha }\cos \delta ,{\mu }_{\delta })=(-3.341\pm 0.028,-6.557\pm 0.043)$ mas yr⁻¹. Upon correction for the effects of viewing perspective and the known cluster rotation, this implies that for the cluster center of mass $({\mu }_{\alpha }\cos \delta ,{\mu }_{\delta })=(-3.238\pm 0.028,-6.716\pm 0.043)$ mas yr⁻¹. This measurement is direct and independent, has the highest random and systematic accuracy to date, and will provide an external verification for the upcoming Gaia Data Release 2. It also differs from most reported PMs for ω Cen in the literature by more than 5σ, but consistency checks compared to other recent catalogs yield excellent agreement. We computed the corresponding Galactocentric velocity, calculated the implied orbit of ω Cen in two different Galactic potentials, and compared these orbits to the orbits implied by one of the PM measurements available in the literature. We find a larger (by about 500 pc) perigalactic distance for ω Cen with our new PM measurement, suggesting a larger survival expectancy for the cluster in the Galaxy.

There is substantial and growing observational evidence from the normalized luminosity density in the near-infrared that the local universe is underdense on scales of several hundred megaparsecs. We test whether our parameterization of the observational data of such a "void" is compatible with the latest supernovae type Ia data and with constraints from line-of-sight peculiar-velocity motions of galaxy clusters with respect to the cosmic microwave background rest-frame, known as the linear kinematic Sunyaev–Zel'dovich (kSZ) effect. Our study is based on the large local void (LLV) radial profile observed by Keenan, Barger, and Cowie (KBC) and a theoretical void description based on the Lemaître–Tolman–Bondi model with a nonzero cosmological constant (ΛLTB). We find consistency with the measured luminosity distance–redshift relation on radial scales relevant to the KBC LLV through a comparison with 217 low-redshift supernovae type Ia over the redshift range $0.0233\lt z\lt 0.15$. We assess the implications of the KBC LLV in light of the tension between "local" and "cosmic" measurements of the Hubble constant, H₀. We find that when the existence of the KBC LLV is fully accounted for, this tension is reduced from $3.4\sigma$ to $2.75\sigma$. We find that previous linear kSZ constraints, as well as new ones from the South Pole Telescope and the Atacama Cosmology Telescope, are fully compatible with the existence of the KBC LLV.

The Monoceros Ring (also known as the Galactic Anticenter Stellar Structure) and A13 are stellar overdensities at estimated heliocentric distances of d ∼ 11 kpc and 15 kpc observed at low Galactic latitudes toward the anticenter of our Galaxy. While these overdensities were initially thought to be remnants of a tidally disrupted satellite galaxy, an alternate scenario is that they are composed of stars from the Milky Way (MW) disk kicked out to their current location due to interactions between a satellite galaxy and the disk. To test this scenario, we study the stellar populations of the Monoceros Ring and A13 by measuring the number of RR Lyrae and M giant stars associated with these overdensities. We obtain low-resolution spectroscopy for RR Lyrae stars in the two structures and measure radial velocities to compare with previously measured velocities for M giant stars in the regions of the Monoceros Ring and A13, to assess the fraction of RR Lyrae to M giant stars (f_RR:MG) in A13 and Mon/GASS. We perform velocity modeling on 153 RR Lyrae stars (116 in the Monoceros Ring and 37 in A13) and find that both structures have very low f_RR:MG. The results support a scenario in which stars in A13 and Mon/GASS formed in the MW disk. We discuss a possible association between Mon/GASS, A13, and the Triangulum-Andromeda overdensity based on their similar velocity distributions and f_RR:MG.

Based on extensive air shower simulations, it is shown that electron distributions with respect to two angles determining the electron direction at a given shower age, for a fixed electron energy and lateral distance, are universal. This means that the distributions do not depend on the primary particle energy or mass (thus, neither on the interaction model), shower zenith angle, or shower to shower fluctuations, if they are taken at the same shower age. Together with previous work showing the universality of the distributions of the electron energy, lateral distance (integrated over angles), and angle (integrated over lateral distance) for fixed electron energy, this paper completes a full universal description of the electron states at various shower ages. Analytical parametrizations of the full electron states are given. It is also shown that some distributions can be described by a number of variables smaller than five, with the new ones being products of old ones raised to some power. The accuracy of the present parametrization is sufficiently good to apply to showers with a primary energy uncertainty of 14% (as is the case at the Pierre Auger Observatory). The shower fluctuations in the chosen bins of the multidimensional variable space are about 6%, determining the minimum uncertainty needed for the parametrization of the universal distributions. An analytical way of estimating the effect of the geomagnetic field is given. Thanks to the universality of the electron distribution in any shower, a new method of shower reconstruction can be worked out from the data from observatories using the fluorescence technique. The light fluxes (both fluorescence and Cherenkov) for any shower age can be exactly predicted for a shower with any primary energy and shower maximum depth, so that the two quantities can be obtained by best fitting the predictions to the measurements.

We present the 2–100 keV spectral analysis of 30 candidate Compton-thick-(CT-)active galactic nuclei (AGNs) selected in the Swift-Burst Alert Telescope (BAT) 100 month survey. The average redshift of these objects is $\langle z\rangle \sim 0.03$, and they all lie within ∼500 Mpc. We used the MyTorus model to perform X-ray spectral fittings both without and with the contribution of the Nuclear Spectroscopic Telescope Array (NuSTAR) data in the 3–50 keV energy range. When the NuSTAR data are added to the fit, 13 out of 30 of these objects (43% of the whole sample) have intrinsic absorption N_H < 10²⁴ cm⁻² at the >3σ confidence level, i.e., they are reclassified from Compton thick to Compton thin. Consequently, we infer an overall observed fraction of the CT-AGN, with respect to the whole AGN population, lower than the one reported in previous works, as low as ∼4%. We find evidence that this overestimation of N_H is likely due to the low quality of a subsample of spectra, either in the 2–10 keV band or in the Swift-BAT one.

While monitoring a sample of apparently slowly rotating superficially normal early-A stars, we have discovered that HR 8844 (A0 V) is actually a new chemically peculiar star. We first compared the high-resolution spectrum of HR 8844 with that of four slow rotators near A0V (ν Cap, ν Cnc, Sirius A, and HD 72660) to highlight similarities and differences. The lines of Ti ii, Cr ii, Sr ii, and Ba ii are conspicuous features in the high-resolution high signal-to-noise SOPHIE spectra of HR 8844 and much stronger than in the spectra of the normal star ν Cap. The Hg ii line at 3983.93 Å is also present in a 3.5% blend. Selected unblended lines of 31 chemical elements from He up to Hg have been synthesized using model atmospheres computed with ATLAS9 and the spectrum synthesis code SYNSPEC48 including hyperfine structure of various isotopes when relevant. These synthetic spectra have been adjusted to the mean SOPHIE spectrum of HR 8844, and high-resolution spectra of the comparison stars. Chi-squares were minimized to derive abundances or upper limits to the abundances of these elements for HR 8844 and the comparison stars. HR 8844 is found to have underabundances of He, C, O, Mg, Ca, and Sc, mild enhancements of Ti, V, Cr, Mn, and distinct enhancements of the heavy elements Sr, Y, Zr, Ba, La, Pr, Sm, Eu, and Hg, the overabundances increasing steadily with atomic number. This chemical pattern suggests that HR 8844 may actually be a new transition object between the coolest HgMn stars and the Am stars.

Relativistic electrons accelerated by both the first-order and the second-order Fermi accelerations in some synchrotron sources have a hybrid shape of thermal and nonthermal energy distribution. This particle acceleration result is supported by some recent numerical simulations. We calculate the synchrotron polarization by applying this electron energy distribution. The polarization degrees in the cases of active galactic nucleus jets and gamma-ray bursts are given as examples. The possible application for the polarization study of Sgr A* is also mentioned. We finally suggest high-energy polarization measurements for these synchrotron sources to test our results.

The detonation of a sub-Chandrasekhar-mass white dwarf (WD) has emerged as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. Recent studies have suggested that the rapid transfer of a very small amount of helium from one WD to another is sufficient to ignite a helium shell detonation that subsequently triggers a carbon core detonation, yielding a "dynamically driven double-degenerate double-detonation" SN Ia. Because the helium shell that surrounds the core explosion is so minimal, this scenario approaches the limiting case of a bare C/O WD detonation. Motivated by discrepancies in previous literature and by a recent need for detailed nucleosynthetic data, we revisit simulations of naked C/O WD detonations in this paper. We disagree to some extent with the nucleosynthetic results of previous work on sub-Chandrasekhar-mass bare C/O WD detonations; for example, we find that a median-brightness SN Ia is produced by the detonation of a 1.0 ${M}_{\odot }$ WD instead of a more massive and rarer 1.1 ${M}_{\odot }$ WD. The neutron-rich nucleosynthesis in our simulations agrees broadly with some observational constraints, although tensions remain with others. There are also discrepancies related to the velocities of the outer ejecta and light curve shapes, but overall our synthetic light curves and spectra are roughly consistent with observations. We are hopeful that future multidimensional simulations will resolve these issues and further bolster the dynamically driven double-degenerate double-detonation scenario's potential to explain most SNe Ia.

We present the first scattered-light images of two debris disks around the F8 star HD 104860 and the F0V star HD 192758, respectively ∼45 and ∼67 pc away. We detected these systems in the F110W and F160W filters through our reanalysis of archival Hubble Space Telescope (HST) NICMOS data with modern starlight-subtraction techniques. Our image of HD 104860 confirms the morphology previously observed by Herschel in thermal emission with a well-defined ring at a radius of ∼114 au inclined by ∼58°. Although the outer edge profile is consistent with dynamical evolution models, the sharp inner edge suggests sculpting by unseen perturbers. Our images of HD 192758 reveal a disk at radius ∼95 au inclined by ∼59°, never resolved so far. These disks have low scattering albedos of 10% and 13%, respectively, inconsistent with water ice grain compositions. They are reminiscent of several other disks with similar inclination and scattering albedos: Fomalhaut, HD 92945, HD 202628, and HD 207129. They are also very distinct from brighter disks in the same inclination bin, which point to different compositions between these two populations. Varying scattering albedo values can be explained by different grain porosities, chemical compositions, or grain size distributions, which may indicate distinct formation mechanisms or dynamical processes at work in these systems. Finally, these faint disks with large infrared excesses may be representative of an underlying population of systems with low albedo values. Searches with more sensitive instruments on HST or on the James Webb Space Telescope and using state-of-the art starlight-subtraction methods may help discover more of such faint systems.

We discuss the production of ultra-high-energy cosmic-ray (UHECR) nuclei and neutrinos from blazars. We compute the nuclear cascade in the jet for both BL Lac objects and flat-spectrum radio quasars (FSRQs), and in the ambient radiation zones for FSRQs as well. By modeling representative spectral energy distributions along the blazar sequence, two distinct regimes are identified, which we call "nuclear survival" (typically found in low-luminosity BL Lacs) and "nuclear cascade" (typically found in high-luminosity FSRQs). We quantify how the neutrino and cosmic-ray (CR) emission efficiencies evolve over the blazar sequence, and we demonstrate that neutrinos and CRs come from very different object classes. For example, high-frequency-peaked BL Lacs (HBLs) tend to produce CRs, and high-luminosity FSRQs are the more efficient neutrino emitters. This conclusion does not depend on the CR escape mechanism, for which we discuss two alternatives (diffusive and advective escape). Finally, the neutrino spectrum from blazars is shown to significantly depend on the injection composition into the jet, especially in the nuclear cascade case: Injection compositions heavier than protons lead to reduced neutrino production at the peak, which moves at the same time to lower energies. Thus, these sources will exhibit better compatibility with the observed IceCube and UHECR data.

Optical polarimetry is an effective way of probing the environment of a supernova for dust. We acquired linear HST ACS/WFC polarimetry in bands $F475W$, $F606W$, and $F775W$ of the supernova (SN) 2014J in M82 at six epochs from ∼277 days to ∼1181 days after the B-band maximum. The polarization measured at day 277 shows conspicuous deviations from other epochs. These differences can be attributed to at least ∼${10}^{-6}\,{M}_{\odot }$ of circumstellar dust located at a distance of $\sim 5\times {10}^{17}\,\mathrm{cm}$ from the SN. The scattering dust grains revealed by these observations seem to be aligned with the dust in the interstellar medium that is responsible for the large reddening toward the supernova. The presence of this circumstellar dust sets strong constraints on the progenitor system that led to the explosion of SN 2014J; however, it cannot discriminate between single- and double-degenerate models.

Polarized emission is detected in two young nearly edge-on protostellar disks in 343 GHz continuum at ∼50 au (∼012) resolution with Atacama Large Millimeter/submillimeter Array. One disk is in HH 212 (Class 0) and the other in the HH 111 (early Class I) protostellar system. The polarization fraction is ∼1%. The disk in HH 212 has a radius of ∼60 au. The emission is mainly detected from the nearside of the disk. The polarization orientations are almost perpendicular to the disk major axis, consistent with either self-scattering or emission by grains aligned with a poloidal field around the outer edge of the disk because of the optical depth effect and temperature gradient; the presence of a poloidal field would facilitate the launching of a disk wind, for which there is already tentative evidence in the same source. The disk of HH 111 VLA 1 has a larger radius of ∼220 au and is thus more resolved. The polarization orientations are almost perpendicular to the disk major axis in the nearside, but more along the major axis in the farside, forming roughly half of an elliptical pattern there. It appears that toroidal and poloidal magnetic field may explain the polarization on the near and far sides of the disk, respectively. However, it is also possible that the polarization is due to self-scattering. In addition, alignment of dust grains by radiation flux may play a role in the farside. Our observations reveal a diversity of disk polarization patterns that should be taken into account in future modeling efforts.

The DArk Matter Particle Explorer, a space-based high precision cosmic-ray detector, has just reported the new measurement of the total electron plus positron energy spectrum up to 4.6 TeV. A notable feature in the spectrum is the spectral break at ∼0.9 TeV, with the spectral index softening from −3.1 to −3.9. Such a feature is very similar to the knee at the cosmic nuclei energy spectrum. In this work, we propose that the knee-like feature can be explained naturally by assuming that the electrons are accelerated at the supernova remnants (SNRs) and released when the SNRs die out with lifetimes around 10⁵ years. The cut-off energy of those electrons have already decreased to several TeV due to radiative cooling, which may induce the observed TeV spectral break. Another possibility is that the break is induced by a single nearby old SNR. Such a scenario may bring a large electron flux anisotropy that may be observable by the future detectors. We also show that a minor part of electrons escaping during the acceleration in young and nearby SNRs is able to contribute to a several TeV or higher energy region of the spectrum.

We have observed the neutron star low-mass X-ray binary SAX J1810.8−2609 in quiescence with XMM-Newton. SAX J1810.8−2609 is one of the faintest non-pulsing neutron star low-mass X-ray binaries in quiescence and previously only had upper limits on its quiescent thermal emission. We found SAX J1810.8−2609 at the same 0.5–10 keV, unabsorbed luminosity as the previous quiescent observation in 2003, ${L}_{{\rm{X}}}=1.5\times {10}^{32}$ erg s⁻¹. We show that the spectrum requires both thermal and nonthermal components, each contributing approximately half the total emission. The low neutron star luminosity suggests a time-averaged outburst accretion rate of $\dot{M}\approx {10}^{-12}$M_⊙ yr⁻¹, in conflict with its observed outburst activity corresponding to a mass accretion rate that is an order of magnitude larger ($\dot{M}\approx {10}^{-11}$M_⊙ yr⁻¹). Our observation designates SAX J1810.8−2609 more firmly as a member of a population of faint quiescent neutron star LMXBs whose quiescent thermal luminosity is not aligned with standard cooling models.

We present in this article a new method to derive the observed properties of outbursts in relativistic jets. We use the VLBI MOJAVE maps to obtain the light curves, based on the principle that the variability of extragalactic sources, in this case 3C 279 and 4C +29.45, should appear in high resolution observations since 1996 until 2016. The use of the cross entropy method (CE) can accurately determine the ranges of parameters for a sequence of outbursts based on the shock-wave model, where the decay/rise timescale ratio has a small spread and the use of a unique index 1.3 generates a good fit modeled by functions of outbursts and by the model of the three stages. By the CE method, one can automatically get the start epochs as well as the occurrence of rise and decline times of the outbursts in the light curves. The values found are in agreement with the power-law distribution of energy, which shows that the cooling of electrons is a predominant process during the initial phase of the shock model evolution. The results of the decomposition show that the outbursts match the VLBI components observed in jets in addition to showing strong evidence of the peaks in the frequencies of 15.3 GHz. With this, we can model the shock waves with reference to the distance at the core of AGN to obtain the Doppler factor and the Lorentz factor.

The recent detection of gravitational waves and electromagnetic counterparts from the double neutron star merger event GW+EM170817 supports the standard paradigm of short gamma-ray bursts (SGRBs) and kilonovae/macronovae. It is important to reveal the nature of the compact remnant left after the merger, either a black hole or neutron star, and their physical link to the origin of the long-lasting emission observed in SGRBs. The diversity of the merger remnants may also lead to different kinds of transients that can be detected in future. Here we study the high-energy emission from the long-lasting central engine left after the coalescence, under certain assumptions. In particular, we consider the X-ray emission from a remnant disk and the nonthermal nebular emission from disk-driven outflows or pulsar winds. We demonstrate that late-time X-ray and high-frequency radio emission can provide useful constraints on properties of the hidden compact remnants and their connections to long-lasting SGRB emission, and we discuss the detectability of nearby merger events through late-time observations at ∼30–100 days after the coalescence. We also investigate the GeV–TeV gamma-ray emission that occurs in the presence of long-lasting central engines and show the importance of external inverse Compton radiation due to upscattering of X-ray photons by relativistic electrons in the jet. We also search for high-energy gamma rays from GW170817 in the Fermi-LAT data and report upper limits on such long-lasting emission. Finally, we consider the implications of GW+EM170817 and discuss the constraints placed by X-ray and high-frequency radio observations.

Polycyclic aromatic hydrocarbons (PAHs), comprised of fused benzene (C₆H₆) rings, emit infrared radiation (3–12 μm) due to the vibrational transitions of the C–H bonds of the aromatic rings. The 3.3 μm aromatic band is generally accompanied by the band at 3.4 μm assigned to the vibration of aliphatic C–H bonds of compounds such as PAHs with an excess of peripheral H atoms (H_n–PAHs). Herein we study the stability of fully hydrogenated benzene (or cyclohexane, C₆H₁₂) under the impact of stellar radiation in the photodissociation region (PDR) of NGC 7027. Using synchrotron radiation and time-of-flight mass spectrometry, we investigated the ionization and dissociation processes at energy ranges of UV (10–200 eV) and soft X-rays (280–310 eV). Density Functional Theory (DFT) calculations were used to determine the most stable structures and the relevant low-lying isomers of singly charged C₆H₁₂ ions. Partial Ion Yield (PIY) analysis gives evidence of the higher tendency toward dissociation of cyclohexane in comparison to benzene. However, because of the high photoabsorption cross-section of benzene at the C1s resonance edge, its photodissociation and photoionization cross-sections are enhanced, leading to a higher efficiency of dissociation of benzene in the PDR of NGC 7027. We suggest that a similar effect is experienced by PAHs in X-ray photon-rich environments, which ultimately acts as an auxiliary protection mechanism of super-hydrogenated polycyclic hydrocarbons. Finally, we propose that the single photoionization of cyclohexane could enhance the abundance of branched molecules in interstellar and circumstellar media.

Forward modeling of the full galaxy SED is a powerful technique, providing self-consistent constraints on stellar ages, dust properties, and metallicities. However, the accuracy of these results is contingent on the accuracy of the model. One significant source of uncertainty is the contribution of obscured AGN, as they are relatively common and can produce substantial mid-IR (MIR) emission. Here we include emission from dusty AGN torii in the Prospector SED-fitting framework, and fit the UV–IR broadband photometry of 129 nearby galaxies. We find that 10% of the fitted galaxies host an AGN contributing >10% of the observed galaxy MIR luminosity. We demonstrate the necessity of this AGN component in the following ways. First, we compare observed spectral features to spectral features predicted from our model fit to the photometry. We find that the AGN component greatly improves predictions for observed Hα and Hβ luminosities, as well as mid-infrared Akari and Spitzer/IRS spectra. Second, we show that inclusion of the AGN component changes stellar ages and SFRs by up to a factor of 10, and dust attenuations by up to a factor of 2.5. Finally, we show that the strength of our model AGN component correlates with independent AGN indicators, suggesting that these galaxies truly host AGN. Notably, only 46% of the SED-detected AGN would be detected with a simple MIR color selection. Based on these results, we conclude that SED models which fit MIR data without AGN components are vulnerable to substantial bias in their derived parameters.

We present results from simulations of core-collapse supernovae in FLASH using a newly implemented multidimensional neutrino transport scheme and a newly implemented general relativistic (GR) treatment of gravity. We use a two-moment method with an analytic closure (so-called M1 transport) for the neutrino transport. This transport is multienergy, multispecies, velocity dependent, and truly multidimensional, i.e., we do not assume the commonly used "ray-by-ray" approximation. Our GR gravity is implemented in our Newtonian hydrodynamics simulations via an effective relativistic potential that closely reproduces the GR structure of neutron stars and has been shown to match GR simulations of core collapse quite well. In axisymmetry, we simulate core-collapse supernovae with four different progenitor models in both Newtonian and GR gravity. We find that the more compact proto–neutron star structure realized in simulations with GR gravity gives higher neutrino luminosities and higher neutrino energies. These differences in turn give higher neutrino heating rates (upward of ∼20%–30% over the corresponding Newtonian gravity simulations) that increase the efficacy of the neutrino mechanism. Three of the four models successfully explode in the simulations assuming GREP gravity. In our Newtonian gravity simulations, two of the four models explode, but at times much later than observed in our GR gravity simulations. Our results, in both Newtonian and GR gravity, compare well with several other studies in the literature. These results conclusively show that the approximation of Newtonian gravity for simulating the core-collapse supernova central engine is not acceptable. We also simulate four additional models in GR gravity to highlight the growing disparity between parameterized 1D models of core-collapse supernovae and the current generation of 2D models.

We present observations of an M5.7 white-light flare (WLF) associated with a small filament eruption in NOAA active region 11476 on 2012 May 10. During this flare, a circular flare ribbon appeared in the east and a remote brightening occurred in the northwest of the active region. Multi-wavelength data are employed to analyze the WLF, including white light (WL), ultraviolet, extreme ultraviolet, hard X-ray (HXR), and microwave. A close spatial and temporal relationship between the WL, HXR, and microwave emissions is found in this WLF. However, the peak time of the WL emission lagged that of the HXR and microwave emissions by about 1–2 minutes. Such a result tends to support the backwarming mechanism for the WL emission. Interestingly, the enhanced WL emission occurred at the two footpoints of the filament. Through forced and potential field extrapolations, we find that the 3D magnetic field in the flare region has a fan-spine feature and that a flux rope lies under the dome-like field structure. We describe the entire process of flare evolution into several steps, each one producing the sequent brightening below the filament, the circular flare ribbons, and the WL enhancement, respectively. We suggest that a reconnection between the magnetic field of the filament and the overlying magnetic field or reconnection within the flux rope leads to the WL enhancement.

The shape of line-of-sight velocity distributions (LOSVDs) carries important information about the internal dynamics of galaxies. The skewness of LOSVDs represents their asymmetric deviation from a Gaussian profile. Correlations between the skewness parameter (h₃) and the mean velocity ($\overline{V}$) of a Gauss–Hermite series reflect the underlying stellar orbital configurations of different morphological components. Using two self-consistent N-body simulations of disk galaxies with different bar strengths, we investigate ${h}_{3}-\overline{V}$ correlations at different inclination angles. Similar to previous studies, we find anticorrelations in the disk area, and positive correlations in the bar area when viewed edge-on. However, at intermediate inclinations, the outer parts of bars exhibit anticorrelations, while the core areas dominated by the boxy/peanut-shaped (B/PS) bulges still maintain weak positive correlations. When viewed edge-on, particles in the foreground/background disk (the wing region) in the bar area constitute the main velocity peak, whereas the particles in the bar contribute to the high-velocity tail, generating the ${h}_{3}-\overline{V}$ correlation. If we remove the wing particles, the LOSVDs of the particles in the outer part of the bar only exhibit a low-velocity tail, resulting in a negative ${h}_{3}-\overline{V}$ correlation, whereas the core areas in the central region still show weakly positive correlations. We discuss implications for IFU observations on bars, and show that the variation of the ${h}_{3}-\overline{V}$ correlation in the disk galaxy may be used as a kinematic indicator of the bar and the B/PS bulge.

We present detailed results of Swift observations of the nearby TeV-detected blazar Mrk 421, based on the rich archival data obtained during 2005 March–2008 June. The best fits of the 0.3–10 keV spectra were mainly obtained using the log-parabolic model, yielding low spectral curvatures expected in the case of the efficient stochastic acceleration of particles. During strong X-ray flares, the position of the synchrotron spectral energy distribution peak ${E}_{{\rm{p}}}$ was beyond 8 keV for 41 spectra, while it sometimes was situated at the UV frequencies in quiescent states. The photon index at 1 keV exhibited a broad range, and the values $a\lt 1.70$ were observed during the strong flares, hinting at the possible presence of a jet hadronic component. The spectral parameters were correlated in some periods, expected in the framework of the first- and second-order Fermi accelerations of X-ray emitting particles, as well as in the case of turbulence spectrum. The 0.3–10 keV flux and spectral parameters sometimes showed very fast variability down to the fluctuations by 6–20% in 180–960 s, possibly related to the small-scale turbulent areas containing strongest magnetic fields. X-ray and very high-energy fluxes often showed correlated variability, although several occurrences of more complicated variability patterns are also revealed, indicating that the multifrequency emission of Mrk 421 could not be generated in a single zone.

Recent advances in helioseismology, numerical simulations and mean-field theory of solar differential rotation have shown that the meridional circulation pattern may consist of two or more cells in each hemisphere of the convection zone. According to the mean-field theory the double-cell circulation pattern can result from the sign inversion of a nondiffusive part of the radial angular momentum transport (the so-called Λ-effect) in the lower part of the solar convection zone. Here, we show that this phenomenon can result from the radial inhomogeneity of the Coriolis number, which depends on the convective turnover time. We demonstrate that if this effect is taken into account then the solar-like differential rotation and the double-cell meridional circulation are both reproduced by the mean-field model. The model is consistent with the distribution of turbulent velocity correlations determined from observations by tracing motions of sunspots and large-scale magnetic fields, indicating that these tracers are rooted just below the shear layer.

By using the Hectospec 6.5 m Multiple Mirror Telescope and the 2.16 m telescope of the National Astronomical Observatories, of the Chinese Academy of Sciences, we obtained 188 high signal-to-noise ratio spectra of ${\rm{H}}\,{\rm{II}}$ regions in the nearby galaxy M101, which is the largest spectroscopic sample of ${\rm{H}}\,{\rm{II}}$ regions for this galaxy so far. These spectra cover a wide range of regions on M101, which enables us to analyze two-dimensional distributions of its physical properties. The physical parameters are derived from emission lines or stellar continua, including stellar population age, electron temperature, oxygen abundance, etc. The oxygen abundances are derived using two empirical methods based on O3N2 and R₂₃ indicators, as well as the direct ${T}_{e}$ method when $[{\rm{O}}\,{\rm{III}}]\,\lambda 4363$ is available. By applying the harmonic decomposition analysis to the velocity field, we obtained a line-of-sight rotation velocity of 71 $\mathrm{km}\,{{\rm{s}}}^{-1}$ and a position angle of 36°. The stellar age profile shows an old stellar population in the galaxy center and a relatively young stellar population in outer regions, suggesting an old bulge and a young disk. The oxygen abundance profile exhibits a clear break at ∼18 kpc, with a gradient of −0.0364 dex kpc⁻¹ in the inner region and −0.00686 dex kpc⁻¹ in the outer region. Our results agree with the "inside-out" disk growth scenario of M101.

We performed spectroscopic observations of the small-separation lensed quasar SDSS J1001+5027, whose images have an angular separation $\theta =2\buildrel{\prime\prime}\over{.} 86$, and placed constraints on the physical properties of gas clouds in the vicinity of the quasar (i.e., in the outflowing wind launched from the accretion disk). The two cylinders of sight to the two lensed images go through the same region of the outflowing wind and they become fully separated with no overlap at a very large distance from the source (∼330 pc). We discovered a clear difference in the profile of the C iv broad absorption line (BAL) detected in the two lensed images in two observing epochs. Because the kinematic components in the BAL profile do not vary in concert, the observed variations cannot be reproduced by a simple change of ionization state. If the variability is due to gas motion around the background source (i.e., the continuum source), the corresponding rotational velocity is ${v}_{\mathrm{rot}}$ ≥ 18,000 km s⁻¹, and their distance from the source is $r\leqslant 0.06$ pc assuming Keplerian motion. Among three Mg ii and three C iv NAL systems that we detected in the spectra, only the Mg ii system at ${z}_{\mathrm{abs}}$ = 0.8716 shows a hint of variability in its Mg i profile on a rest-frame timescale of ${\rm{\Delta }}{t}_{\mathrm{rest}}$$\leqslant$ 191 days and an obvious velocity shear between the sightlines whose physical separation is ∼7 kpc. We interpret this as the result of motion of a cosmologically intervening absorber, perhaps located in a foreground galaxy.

We have measured the radial profiles of isophotal ellipticity (ε) and disky/boxy parameter A₄ out to radii of about three times the semimajor axes for ∼4600 star-forming galaxies (SFGs) between redshift 0.5 and 1.8 in the CANDELS/GOODS-S and UDS fields. Based on the average size–mass relation in each redshift bin, we divide our galaxies at a given mass into Small SFGs (SSFGs; smaller than the average) and Large SFGs (LSFGs; larger than the average). We show that, at low masses (${M}_{* }\lt {10}^{10}{M}_{\odot }$), the SSFGs generally have nearly flat ε and A₄ profiles in both edge-on and face-on views, especially at $z\gt 1$. Moreover, the median A₄ values at all radii are almost zero. In contrast, the highly inclined low-mass LSFGs in the same mass-redshift bins generally have monotonically increasing ε profiles with radius and disky feature dominated in the intermediate regions. These findings imply that at these redshifts, the low-mass SSFGs are not disk-like, whereas the low-mass LSFGs likely harbour disk-like components flattened by significant rotations. At high masses (${M}_{* }\gt {10}^{10}{M}_{\odot }$), both highly inclined SSFGs and LSFGs generally exhibit distinct trends in both ε and A₄ profiles, which increase at lower radii, reach maxima, then decrease at larger radii. Such the feature is more prevalent for more massive (${M}_{* }\gt {10}^{10.5}{M}_{\odot }$) galaxies or at lower redshifts ($z\lt 1.4$). This feature can be simply explained if galaxies possess all three components: central bulges, disks in the intermediate regions, and halo-like stellar components in the outskirts.

Supernova remnants (SNRs) are the primary candidate of Galactic cosmic-ray accelerators. It is still an open issue when and how young SNRs, which typically exhibit strong synchrotron X-rays and GeV and TeV gamma rays, undergo the state transition to middle-aged SNRs dominated by thermal X-rays and GeV gamma rays. SNR N132D in the Large Magellanic Cloud is an ideal target to study such a transition, exhibiting bright X-rays and gamma rays, and with an expected age of ∼2500 years. In this paper we present results of NuSTAR and Suzaku spectroscopy. We reveal that N132D has a nearly equilibrium plasma with a temperature of >5 keV or a recombining plasma with a lower temperature (∼1.5 keV) and a recombining timescale (${n}_{e}t$) of $8.8\,(7.0\mbox{--}10.0)\times {10}^{12}$ cm⁻³s. Together with the center-filled morphology observed in the iron K line image, our results suggest that N132D is now at the transition stage from being a young SNR to being middle-aged. We have constrained the tight upper limit of nonthermal X-rays. Bright gamma rays compared to faint nonthermal X-rays suggest that the gamma rays are hadronic in origin. The spectral energy distribution from radio to gamma rays shows a proton cutoff energy of ∼30 TeV. These facts confirm that N132D is undergoing the transition from a young to a middle-aged SNR. The large thermal energy of $\gt {10}^{51}$ erg and accelerated proton energy of $\sim {10}^{50}$ erg suggest the supernova explosion might have been very energetic.

We consider the magnetic interaction of exoplanets orbiting M dwarfs, calculating the expected Poynting flux carried upstream along Alfvén wings to the central star. A region of emission analogous to the Io footprint observed in Jupiter's aurora is produced, and we calculate the radio flux density generated near the surface of the star via the electron-cyclotron maser instability. We apply the model to produce individual case studies for the TRAPPIST-1, Proxima Centauri, and dwarf NGTS-1 systems. We predict steady-state flux densities of up to ∼10 μJy and sporadic bursts of emission of up to ∼1 mJy from each case study, suggesting these systems may be detectable with the Very Large Array and the Giant Metrewave Radio Telescope, and perhaps the Square Kilometre Array in the future. Finally, we present a survey of 85 exoplanets orbiting M dwarfs, identifying 11 such objects capable of generating radio emission above 10 μJy.

We present UV luminosity functions of dropout galaxies at $z\sim 6\mbox{--}10$ with the complete Hubble Frontier Fields data. We obtain a catalog of ∼450 dropout-galaxy candidates (350, 66, and 40 at $z\sim 6\mbox{--}7$, 8, and 9, respectively), with UV absolute magnitudes that reach $\sim -14$ mag, ∼2 mag deeper than the Hubble Ultra Deep Field detection limits. We carefully evaluate number densities of the dropout galaxies by Monte Carlo simulations, including all lensing effects such as magnification, distortion, and multiplication of images as well as detection completeness and contamination effects in a self-consistent manner. We find that UV luminosity functions at $z\sim 6\mbox{--}8$ have steep faint-end slopes, $\alpha \sim -2$, and likely steeper slopes, $\alpha \lesssim -2$ at $z\sim 9\mbox{--}10$. We also find that the evolution of UV luminosity densities shows a non-accelerated decline beyond $z\sim 8$ in the case of ${M}_{\mathrm{trunc}}=-15$, but an accelerated one in the case of ${M}_{\mathrm{trunc}}=-17$. We examine whether our results are consistent with the Thomson scattering optical depth from the Planck satellite and the ionized hydrogen fraction Q_{H ii} at $z\lesssim 7$ based on the standard analytic reionization model. We find that reionization scenarios exist that consistently explain all of the observational measurements with the allowed parameters of ${f}_{\mathrm{esc}}={0.17}_{-0.03}^{+0.07}$ and ${M}_{\mathrm{trunc}}\gt -14.0$ for $\mathrm{log}{\xi }_{\mathrm{ion}}/[{\mathrm{erg}}^{-1}\ \mathrm{Hz}]=25.34$, where ${f}_{\mathrm{esc}}$ is the escape fraction, M_trunc is the faint limit of the UV luminosity function, and ${\xi }_{\mathrm{ion}}$ is the conversion factor of the UV luminosity to the ionizing photon emission rate. The length of the reionization period is estimated to be ${\rm{\Delta }}z={3.9}_{-1.6}^{+2.0}$ (for $0.1\lt {Q}_{{\rm{H}}{\rm{II}}}\lt 0.99$), consistent with the recent estimate from Planck.

Solar oscillation frequencies change with the level of magnetic activity. Localizing subsurface magnetic field concentrations in the Sun with helioseismology will help us to understand the solar dynamo. Because the magnetic fields are not considered in standard solar models, adding them to the basic equations of stellar structure changes the eigenfunctions and eigenfrequencies. We use quasi-degenerate perturbation theory to calculate the effect of toroidal magnetic fields on solar oscillation mean multiplet frequencies for six field configurations. In our calculations, we consider both the direct effect of the magnetic field, which describes the coupling of modes, and the indirect effect, which accounts for changes in stellar structure due to the magnetic field. We limit our calculations to self-coupling of modes. We find that the magnetic field affects the multiplet frequencies in a way that depends on the location and the geometry of the field inside the Sun. Comparing our theoretical results with observed shifts, we find that strong tachocline fields cannot be responsible for the observed frequency shifts of p modes over the solar cycle. We also find that part of the surface effect in helioseismic oscillation frequencies might be attributed to magnetic fields in the outer layers of the Sun. The theory presented here is also applicable to models of solar-like stars and their oscillation frequencies.

Direct observational searches for Population III (Pop III) stars at high redshift are faced with the question of how to select the most promising targets for spectroscopic follow-up. To help answer this, we use a large-scale cosmological simulation, augmented with a new subgrid model that tracks the fraction of pristine gas, to follow the evolution of high-redshift galaxies and the Pop III stars they contain. We generate rest-frame ultraviolet (UV) luminosity functions for our galaxies and find that they are consistent with current $z\geqslant 7$ observations. Throughout the redshift range $7\leqslant z\leqslant 15$, we identify "Pop III–bright" galaxies as those with at least 75% of their flux coming from Pop III stars. While less than 1% of galaxies brighter than ${m}_{\mathrm{UV},\mathrm{AB}}=31.4$ mag are Pop III–bright in the range $7\leqslant z\leqslant 8$, roughly 17% of such galaxies are Pop III–bright at z = 9, immediately before reionization occurs in our simulation. Moving to z = 10, ${m}_{\mathrm{UV},\mathrm{AB}}=31.4$ mag corresponds to larger, more luminous galaxies, and the Pop III–bright fraction falls off to 5%. Finally, at the highest redshifts, a large fraction (29% at z = 14 and 41% at z = 15) of all galaxies are Pop III–bright regardless of magnitude. While ${m}_{\mathrm{UV},\mathrm{AB}}=31.4$ mag galaxies are extremely rare during this epoch, we find that 13% of galaxies at z = 14 are Pop III–bright with ${m}_{\mathrm{UV},\mathrm{AB}}\leqslant 33$ mag, a intrinsic magnitude within reach of the James Webb Space Telescope using lensing. Thus, we predict that the best redshift to search for luminous Pop III–bright galaxies is just before reionization, while lensing surveys for fainter galaxies should push to the highest redshifts possible.