Table of contents

High precision proper motion (PM) measurements are available for approximately 20% of all known dwarf satellite galaxies of the Milky Way (MW). Here we extend the Bayesian framework of Patel et al. to include all MW satellites with measured 6D phase-space information and apply it with the Illustris-Dark simulation to constrain the MW's mass. Using the properties of each MW satellite individually, we find that the scatter among mass estimates is reduced when the magnitude of specific orbital angular momentum (j) is adopted, rather than their combined instantaneous positions and velocities. We also find that high j satellites (i.e., Leo II) constrain the upper limits for the MW's mass and low j satellites, rather than the highest speed satellites (i.e., Leo I and Large Magellanic Cloud), set the lower mass limits. When j of all classical satellites is used to simultaneously estimate the MW's mass, we conclude the halo mass is 0.85^+0.23_−0.26 × 10¹²${M}_{\odot }$ (including Sagittarius dSph) and 0.96^+0.29_−0.28 × 10¹²${M}_{\odot }$ (excluding Sagittarius dSph), cautioning that low j satellites on decaying orbits like Sagittarius dSph may bias the distribution. These estimates markedly reduce the current factor of two spread in the mass range of the MW. We also find a well-defined relationship between host halo mass and satellite j distribution, which yields the prediction that upcoming PMs for ultra-faint dwarfs should reveal j within 5 × 10³–10⁴ kpc km s⁻¹. This is a promising method to significantly constrain the cosmologically expected mass range for the MW and eventually M31 as more satellite PMs become available.

We have selected a sample of 11 massive clusters of galaxies observed by the Hubble Space Telescope in order to study the impact of the dynamical state on the intracluster light (ICL) fraction, the ratio of total integrated ICL to the total galaxy member light. With the exception of the Bullet cluster, the sample is drawn from the Cluster Lensing and Supernova Survey and the Frontier Fields program, containing five relaxed and six merging clusters. The ICL fraction is calculated in three optical filters using the CHEFs ICL estimator, a robust and accurate algorithm free of a priori assumptions. We find that the ICL fraction in the three bands is, on average, higher for the merging clusters, ranging between ∼7% and 23%, compared with the ∼2%–11% found for the relaxed systems. We observe a nearly constant value (within the error bars) in the ICL fraction of the regular clusters at the three wavelengths considered, which would indicate that the colors of the ICL and the cluster galaxies are, on average, coincident and, thus, so are their stellar populations. However, we find a higher ICL fraction in the F606W filter for the merging clusters, consistent with an excess of lower-metallicity/younger stars in the ICL, which could have migrated violently from the outskirts of the infalling galaxies during the merger event.

The morphology and optical spectrum of IPHASX J210204.7+471015, a nebula classified as a possible planetary nebula are, however, strikingly similar to those of AT Cnc, a classical nova shell around a dwarf nova. To investigate its true nature, we have obtained high-resolution narrowband [O iii] and [N ii] images and deep optical spectra. The nebula shows an arc of [N ii]-bright knots notably enriched in nitrogen, while an [O iii]-bright bow shock is progressing throughout the ISM. Diagnostic line ratios indicate that shocks are associated with the arc and bow shock. The central star of this nebula has been identified by its photometric variability. Time-resolved photometric and spectroscopic data of this source reveal a period of 4.26 hr, which is attributed to a binary system. The optical spectrum is notably similar to that of RW Sex, a cataclysmic variable star (CV) of the UX UMa nova-like (NL) type. Based on these results, we propose that IPHASX J210204.7 + 471015 is a classical nova shell observed around a CV-NL system in quiescence.

We have recently witnessed the first multi-messenger detection of colliding neutron stars through gravitational waves (GWs) and electromagnetic (EM) waves (GW 170817) thanks to the joint efforts of LIGO/Virgo and Space/Ground-based telescopes. In this paper, we report on the RATIR follow-up observation strategies and show the results for the trigger G194575. This trigger is not of astrophysical interest; however, it is of great interest to the robust design of a follow-up engine to explore large sky-error regions. We discuss the development of an image-subtraction pipeline for the six-color, optical/NIR imaging camera RATIR. Considering a two-band (i and r) campaign in the fall of 2015, we find that the requirement of simultaneous detection in both bands leads to a factor ∼10 reduction in false alarm rate, which can be further reduced using additional bands. We also show that the performance of our proposed algorithm is robust to fluctuating observing conditions, maintaining a low false alarm rate with a modest decrease in system efficiency that can be overcome utilizing repeat visits. Expanding our pipeline to search for either optical or NIR detections (three or more bands), considering separately the optical riZ and NIR YJH bands, should result in a false alarm rate ≈1% and an efficiency ≈90%. RATIR's simultaneous optical/NIR observations are expected to yield about one candidate transient in the vast 100 deg² LIGO error region for prioritized follow-up with larger aperture telescopes.

Type III bursts are generated by fast electron beams originated from magnetic reconnection sites of solar flares. As propagation of radio waves in the interplanetary medium is strongly affected by random electron density fluctuations, type III bursts provide us with a unique diagnostic tool for solar wind remote plasma measurements. Here, we performed a statistical survey of 152 simple and isolated type III bursts observed by the twin-spacecraft Solar TErrestrial RElations Observatory mission. We investigated their time–frequency profiles in order to retrieve decay times as a function of frequency. Next, we performed Monte Carlo simulations to study the role of scattering due to random electron density fluctuations on time–frequency profiles of radio emissions generated in the interplanetary medium. For simplification, we assumed the presence of isotropic electron density fluctuations described by a power law with the Kolmogorov spectral index. Decay times obtained from observations and simulations were compared. We found that the characteristic exponential decay profile of type III bursts can be explained by the scattering of the fundamental component between the source and the observer despite restrictive assumptions included in the Monte Carlo simulation algorithm. Our results suggest that relative electron density fluctuations $\langle \delta {n}_{{\rm{e}}}\rangle /{n}_{{\rm{e}}}$ in the solar wind are 0.06–0.07 over wide range of heliospheric distances.

This study aims at investigating the ability of strong interactions between magnetic field concentrations during their rise through the convection zone to produce complex active regions at the solar surface. To do so, we perform numerical simulations of buoyant magnetic structures evolving and interacting in a model solar convection zone. We first produce a 3D model of rotating convection and then introduce idealized magnetic structures close to the bottom of the computational domain. These structures possess a certain degree of field line twist and they are made buoyant on a particular extension in longitude. The resulting twisted Ω-loops will thus evolve inside a spherical convective shell possessing large-scale mean flows. We present results on the interaction between two such loops with various initial parameters (mainly buoyancy and twist) and on the complexity of the emerging magnetic field. In agreement with analytical predictions, we find that if the loops are introduced with opposite handedness and same axial field direction or the same handedness but opposite axial field, they bounce against each other. The emerging region is then constituted of two separated bipolar structures. On the contrary, if the loops are introduced with the same direction of axial and peripheral magnetic fields and are sufficiently close, they merge while rising. This more interesting case produces complex magnetic structures with a high degree of non-neutralized currents, especially when the convective motions act significantly on the magnetic field. This indicates that those interactions could be good candidates to produce eruptive events like flares or CMEs.

We investigate the role of active galactic nucleus (AGN) feedback in turbulent heating of galaxy clusters. Specifically, we analyze the production of turbulence by g-modes generated by the supersonic expansion and buoyant rise of AGN-driven bubbles. Previous work that neglects magnetic fields has shown that this process is inefficient, with less than 1% of the injected energy ending up in turbulence. This inefficiency primarily arises because the bubbles are shredded apart by hydrodynamic instabilities before they can excite sufficiently strong g-modes. Using a plane-parallel model of the intracluster medium (ICM) and 3D ideal magnetohydrodynamics (MHD) simulations, we examine the role of a large-scale magnetic field that is able to drape around these rising bubbles, preserving them from hydrodynamic instabilities. We find that while magnetic draping appears better able to preserve AGN-driven bubbles, the driving of g-modes and the resulting production of turbulence is still inefficient. The magnetic tension force prevents g-modes from transitioning into the nonlinear regime, suppressing turbulence in our model ICM. Our work highlights the ways in which ideal MHD is an insufficient description for the cluster feedback process, and we discuss future work such as the inclusion of anisotropic viscosity as a means of simulating high β plasma kinetic effects. These results suggest the hypothesis that other mechanisms of heating the ICM plasma such as sound waves or cosmic rays may be responsible for the observed feedback in galaxy clusters.

Why the tenuous solar outer atmosphere, or corona, is much hotter than the underlying layers remains one of the greatest challenges for solar modeling. Detailed diagnostics of the coronal thermal structure come from extreme ultraviolet (EUV) emission. The EUV emission is produced by heavy ions in various ionization states and depends on the amount of these ions and on plasma temperature and density. Any nonuniformity of the elemental distribution in space or variability in time affects thermal diagnostics of the corona. Here we theoretically predict ionized chemical element concentrations in some areas of the solar atmosphere, where the electric current is directed upward. We then detect these areas observationally, by comparing the electric current density with the EUV brightness in an active region. We found a significant excess in EUV brightness in the areas with positive current density rather than negative. Therefore, we report the observational discovery of substantial concentrations of heavy ions in current-carrying magnetic flux tubes, which might have important implications for the elemental fractionation in the solar corona known as the first ionization potential effect. We call such areas of heavy ion concentration the "ion traps." These traps hold enhanced ion levels until they are disrupted by a flare, whether large or small.

NGC 5548 has been intensively monitored by the AGN Space Telescope and Optical Reverberation Mapping collaboration. Approximately after half of the light curves, the correlation between the broad emission lines and the lag-corrected ultraviolet (UV) continua becomes weak. This anomalous behavior is accompanied by an increase of soft X-ray emission. We propose a simple model to understand this anomalous behavior, i.e., the corona might fall down, thereby increasing the covering fraction of the inner disk. Therefore, X-ray and extreme-UV emission suffer from spectral variations. The UV continua variations are driven by both X-ray and extreme-UV variations. Consequently, the spectral variability induced by the falling corona would dilute the correlation between the broad emission lines and the UV continua. Our model can explain many additional observational facts, including the dependence of the anomalous behavior on velocity and ionization energy. We also show that the time lag and correlation between the X-ray and the UV variations change as NGC 5548 displays the anomalous behavior. The time lag is dramatically longer than the expectation from disk reprocessing if the anomalous behavior is properly excluded. During the anomalous state, the time lag approaches the light-travel timescale of disk reprocessing albeit with a much weaker correlation. We speculate that the time lag in the normal state is caused by reprocessing of the broad line region gas. As NGC 5548 enters the abnormal state, the contribution of the broad line region gas is smaller; the time lag reflects disk reprocessing. We also discuss alternative scenarios.

Recent Atacama Large Millimeter and Submillimeter Array (ALMA) observations of the protoplanetary disk around the Herbig Ae star HD 163296 revealed three depleted dust gaps at 60, 100, and 160 au in the 1.3 mm continuum as well as CO depletion in the middle and outer dust gaps. However, no CO depletion was found in the inner dust gap. To examine the planet–disk interaction model, we present results of 2D two fluid (gas + dust) hydrodynamic simulations coupled with 3D radiative transfer simulations. To fit the high gas-to-dust ratio of the first gap, we find that the Shakura–Sunyaev viscosity parameter α must be very small ($\lesssim {10}^{-4}$) in the inner disk. On the other hand, a relatively large α ($\sim 7.5\times {10}^{-3}$) is required to reproduce the dust surface density in the outer disk. We interpret the variation of α as an indicator of the transition from an inner dead zone to the outer magnetorotational instability (MRI) active zone. Within ∼100 au, the HD 163296 disk's ionization level is low, and non-ideal magnetohydrodynamic effects could suppress the MRI, so the disk can be largely laminar. The disk's ionization level gradually increases toward larger radii, and the outermost disk ($r\gt 300$ au) becomes turbulent due to MRI. Under this condition, we find that the observed dust continuum and CO gas line emissions can be reasonably fit by three half-Jovian-mass planets (0.46, 0.46, and $0.58\,{M}_{{\rm{J}}}$) at 59, 105, and 160 au, respectively.

We present Hubble Space Telescope observations and photometric measurements of the Type Ia supernova (SN Ia) SN 2013aa 1500 days after explosion. At this epoch, the luminosity is primarily dictated by the amounts of radioactive ${}^{57}{\rm{Co}}$ and ${}^{55}{\rm{Fe}}$, while at earlier epochs, the luminosity depends on the amount of radioactive ${}^{56}{\rm{Co}}$. The ratio of odd-numbered to even-numbered isotopes depends significantly on the density of the progenitor white dwarf (WD) during the SN explosion, which, in turn, depends on the details of the progenitor system at the time of ignition. From a comprehensive analysis of the entire light curve of SN 2013aa, we measure a $M({}^{57}{\rm{Co}})/M({}^{56}{\rm{Co}})$ ratio of ${0.02}_{-0.02}^{+0.01}$, which indicates a relatively low central density for the progenitor WD at the time of explosion, consistent with DD progenitor channels. We estimate $M({}^{56}{\rm{Ni}})\,=\,0.732\,\pm \,0.151\ {M}_{\odot }$, and place an upper limit on the abundance of ${}^{55}{\rm{Fe}}$. A recent study reported a possible correlation between $M({}^{57}{\rm{Co}})/M({}^{56}{\rm{Co}})$ and stretch for four SNe Ia. SN 2013aa, however, does not fit this trend, indicating either SN 2013aa is an extreme outlier or the correlation does not hold up with a larger sample. The $M({}^{57}{\rm{Co}})/M({}^{56}{\rm{Co}})$ measured for the expanded sample of SNe Ia with photometry at extremely late times has a much larger range than that of explosion models, perhaps limiting conclusions about SN Ia progenitors drawn from extremely late-time photometry.

Infrared spectra of ices containing alcohols that are known or potential interstellar molecules are examined before and after irradiation with 1 MeV protons at ∼20 K. The low-temperature oxidation (hydrogen loss) of six alcohols is followed, and conclusions are drawn based on the results. The formation of reaction products is discussed in terms of the literature on the radiation chemistry of alcohols and a systematic variation in their structures. The results from these new laboratory measurements are then applied to a recent study of propargyl alcohol. Connections are drawn between known interstellar molecules, and several new reaction products in interstellar ices are predicted.

We study the magnetic structure of a successively erupting sigmoid in active region 12371 by modeling the quasi-static coronal field evolution with nonlinear force-free field (NLFFF) equilibria. Helioseismic and Magnetic Imager/Solar Dynamic Observatory vector magnetograms are used as input to the NLFFF model. In all eruption events, the modeled structure resembles the observed pre-eruptive coronal sigmoid and the NLFFF core field is a combination of double inverse-J-shaped and inverse-S field lines with dips touching the photosphere. Such field lines are formed by the flux cancellation reconnection of opposite-J field lines at bald-patch locations, which in turn implies the formation of a weakly twisted flux-rope (FR) from large-scale sheared arcade field lines. Later on, this FR undergoes coronal tether-cutting reconnection until a coronal mass ejection is triggered. The modeled structure captured these major features of sigmoid-to-arcade-to-sigmoid transformation, which is reoccuring under continuous photospheric flux motions. Calculations of the field line twist reveal a fractional increase followed by a decrease of the number of pixels having a range of twist. This traces the buildup process of a twisted core field by slow photospheric motions and the relaxation after eruption, respectively. Our study infers that the large eruptivity of this AR is due to a steep decrease of the background coronal field meeting the torus instability criteria at a low height (≈40 Mm) in contrast to noneruptive ARs.

Efficient numerical approximation of the polarized radiative transfer equation is challenging because this system of ordinary differential equations exhibits stiff behavior, which potentially results in numerical instability. This negatively impacts the accuracy of formal solvers, and small step-sizes are often necessary to retrieve physical solutions. This work presents stability analyses of formal solvers for the radiative transfer equation of polarized light, identifies instability issues, and suggests practical remedies. In particular, the assumptions and the limitations of the stability analysis of Runge–Kutta methods play a crucial role. On this basis, a suitable and pragmatic formal solver is outlined and tested. An insightful comparison to the scalar radiative transfer equation is also presented.

We prove that a second-order logarithmic Hamiltonian leapfrog for the classical general N-body problem (CGNBP) designed by Mikkola and Tanikawa and some higher-order logarithmic Hamiltonian methods based on symmetric multicompositions of the logarithmic algorithm exactly reproduce the orbits of elliptic relative equilibrium solutions in the original CGNBP. These methods are explicit symplectic methods. Before this proof, only some implicit discrete-time CGNBPs proposed by Minesaki had been analytically shown to trace the orbits of elliptic relative equilibrium solutions. The proof is therefore the first existence proof for explicit symplectic methods. Such logarithmic Hamiltonian methods with a variable time step can also precisely retain periodic orbits in the classical general three-body problem, which generic numerical methods with a constant time step cannot do.

In order to assess the relationship between metallicity and exoplanetary systems, we compare the abundances of AF-type main-sequence stars with debris disk properties assessed using Herschel observations of an unbiased survey of nearby stars. Hot stars are not as commonly observed, given their unique constraints in data reduction, lack of metal lines, and "astrophysical noise" from rotation speed. Here, we address that deficiency using new and archival spectra of 83 AF-type stars. We measure the abundances of a few species in addition to Fe in order to classify the stars with Ap/Am or Lambda Boo signatures. Lambda Boo stars have a chemical signature of solar-abundant volatile species and sub-solar refractory abundances that is hypothesized to be altered by the pollution of volatiles. Overall, we see no correlation between debris disks and metallicity, primarily because the sample size is cut significantly when using only reliable fits to the spectroscopic data. The abundance measured from the Mg ii 4481 blend is a useful diagnostic because it can be reliably measured at large v·sin(i) and is found to be lower around stars with bright debris disks. We find that Lambda Boo stars have brighter debris disks compared to a bias-free sample of AF stars. The trend with disk brightness and Mg abundances suggests pollution effects can be significant and used as a marker for the stability of planetary systems. We explore trends with other species, such as with the C/O ratios, but are significantly limited by the low number of reliable detections.

The interstellar abundances of refractory elements indicate a substantial depletion from the gas phase, which increases with gas density. Our recent model of dust evolution, based on hydrodynamic simulations of the life cycle of giant molecular clouds (GMCs), proves that the observed trend for [Si_gas/H] is driven by a combination of dust growth by accretion in the cold diffuse interstellar medium (ISM) and efficient destruction by supernova (SN) shocks. With an analytic model of dust evolution, we demonstrate that even with optimistic assumptions for the dust input from stars and without destruction of grains by SNe it is impossible to match the observed [Si_gas/H]–n_H relation without growth in the ISM. We extend the framework developed in our previous work for silicates to include the evolution of iron grains and address a long-standing conundrum: "Where is the interstellar iron?" Much higher depletion of Fe in the warm neutral medium compared to Si is reproduced by the models, in which a large fraction of interstellar iron (70%) is locked as inclusions in silicate grains, where it is protected from efficient sputtering by SN shocks. The slope of the observed [Fe_gas/H]–n_H relation is reproduced if the remaining depleted iron resides in a population of metallic iron nanoparticles with sizes in the range of 1–10 nm. Enhanced collision rates due to the Coulomb focusing are important for both silicate and iron dust models to match the slopes of the observed depletion–density relations and the magnitudes of depletion at high gas density.

Rapidly spinning, strongly magnetized protoneutron stars ("millisecond protomagnetars") are candidate central engines of long-duration gamma-ray bursts (GRBs), superluminous supernovae (SLSNe), and binary neutron star mergers. Magnetar birth may be accompanied by the fallback of stellar debris, lasting for seconds or longer following the explosion. Accretion alters the magnetar evolution by (1) providing an additional source of rotational energy (or a potential sink, if the propeller mechanism operates), (2) enhancing the spin-down luminosity above the dipole rate by compressing the magnetosphere and expanding the polar cap region of open magnetic field lines, and (3) supplying an additional accretion-powered neutrino luminosity that sustains the wind baryon loading, even after the magnetar's internal neutrino luminosity has subsided. The more complex evolution of the jet power and magnetization of an accreting magnetar more readily accounts for the high ⁵⁶Ni yields of GRB SNe and the irregular time evolution of some GRB light curves (e.g., bursts with precursors followed by a long quiescent interval before the main emission episode). Additional baryon loading from accretion-powered neutrino irradiation of the polar cap lengthens the time frame over which the jet magnetization is in the requisite range σ ≲ 10³ for efficient gamma-ray emission, thereby accommodating GRBs with ultralong durations. Though accretion does not significantly raise the maximum energy budget from the limit of ≲ few × 10⁵² erg for an isolated magnetar, it greatly expands the range of magnetic field strengths and birth spin periods capable of powering GRB jets, reducing the differences between the magnetar properties normally invoked to explain GRBs versus SLSNe.

Recent laboratory efforts have constrained the remanent magnetizations of chondrules and the magnetic field strengths to which the chondrules were exposed as they cooled below their Curie points. An outstanding question is whether the inferred paleofields represent the background magnetic field of the solar nebula or were unique to the chondrule-forming environment. We investigate the amplification of the magnetic field above background values for two proposed chondrule formation mechanisms, large-scale nebular shocks and planetary bow shocks. Behind large-scale shocks, the magnetic field parallel to the shock front is amplified by factors of ∼10–30, regardless of the magnetic diffusivity. Therefore, chondrules melted in these shocks probably recorded an amplified magnetic field. Behind planetary bow shocks, the field amplification is sensitive to the magnetic diffusivity. We compute the gas properties behind a bow shock around a 3000 km radius planetary embryo, with and without atmospheres, using hydrodynamics models. We calculate the ionization state of the hot, shocked gas, including thermionic emission from dust, thermal ionization of gas-phase potassium atoms, and the magnetic diffusivity due to Ohmic dissipation and ambipolar diffusion. We find that the diffusivity is sufficiently large that magnetic fields have already relaxed to background values in the shock downstream where chondrules acquire magnetizations, and that these locations are sufficiently far from the planetary embryos that chondrules should not have recorded a significant putative dynamo field generated on these bodies. We conclude that, if melted in planetary bow shocks, chondrules probably recorded the background nebular field.

A long-standing problem is identifying the elusive progenitors of Type Ia supernovae (SNe Ia), which can roughly be split into Chandraksekhar and sub-Chandrasekhar-mass events. An important difference between these two cases is the nucleosynthetic yield, which is altered by the increased neutron excess in Chandrasekhar progenitors due to their pre-explosion simmering and high central density. Based on these arguments, we show that the chemical composition of the most metal-rich star in the Ursa Minor dwarf galaxy, COS 171, is dominated by nucleosynthesis from a low-metallicity, low-mass, sub-Chandrasekhar-mass SN Ia. Key diagnostic abundance ratios include Mn/Fe and Ni/Fe, which could not have been produced by a Chandrasekhar-mass SN Ia. Large deficiencies of Ni/Fe, Cu/Fe and Zn/Fe also suggest the absence of alpha-rich freeze-out nucleosynthesis, favoring low-mass white dwarf progenitors of SNe Ia, near 0.95 M_⊙, from comparisons to numerical detonation models. We also compare Mn/Fe and Ni/Fe ratios to the recent yields predicted by Shen et al., finding consistent results. To explain the [Fe/H] at −1.35 dex for COS 171 would require dilution of the SN Ia ejecta with ∼10⁴M_⊙ of material, which is expected for an SN remnant expanding into a warm interstellar medium with n ∼ 1 cm⁻³. In the future, finding more stars with the unique chemical signatures we highlight here will be important for constraining the rate and environments of sub-Chandrasekhar SNe Ia.

The advanced rheological models of Andrade and Sundberg & Cooper are compared to the traditional Maxwell model to understand how each affects the tidal dissipation of heat within rocky bodies. We find both Andrade and Sundberg–Cooper rheologies can produce at least 10× the tidal heating compared to a traditional Maxwell model for a warm (1400–1600 K) Io-like satellite. Sundberg–Cooper can cause even larger dissipation around a critical temperature and frequency. These models allow cooler planets to stay tidally active in the face of orbital perturbations—a condition we term "tidal resilience." This has implications for the time evolution of tidally active worlds and the long-term equilibria they fall into. For instance, if Io's interior is better modeled by the Andrade or Sundberg–Cooper rheologies, the number of possible resonance-forming scenarios that still produce a hot, modern Io is expanded, and these scenarios do not require an early formation of the Laplace resonance. The two primary empirical parameters that define the Andrade anelasticity are examined in several phase spaces to provide guidance on how their uncertainties impact tidal outcomes, as laboratory studies continue to constrain their real values. We provide detailed reference tables on the fully general equations required for others to insert the models of Andrade and Sundberg–Cooper into standard tidal formulae. Lastly, we show that advanced rheologies can greatly impact the heating of short-period exoplanets and exomoons, while the properties of tidal resilience could mean a greater number of tidally active worlds among all extrasolar systems.

Recent observations in extreme-ultraviolet (EUV) wavelengths reveal an EUV late phase in some solar flares that is characterized by a second peak in warm coronal emissions (∼3 MK) several tens of minutes to a few hours after the soft X-ray (SXR) peak. Using the model enthalpy-based thermal evolution of loops (EBTEL), we numerically probe the production of EUV late-phase solar flares. Starting from two main mechanisms of producing the EUV late phase, i.e., long-lasting cooling and secondary heating, we carry out two groups of numerical experiments to study the effects of these two processes on the emission characteristics in late-phase loops. In either of the two processes an EUV late-phase solar flare that conforms to the observational criteria can be numerically synthesized. However, the underlying hydrodynamic and thermodynamic evolutions in late-phase loops are different between the two synthetic flare cases. The late-phase peak due to a long-lasting cooling process always occurs during the radiative cooling phase, while that powered by a secondary heating is more likely to take place in the conductive cooling phase. We then propose a new method for diagnosing the two mechanisms based on the shape of EUV late-phase light curves. Moreover, from the partition of energy input, we discuss why most solar flares are not EUV late flares. Finally, by addressing some other factors that may potentially affect the loop emissions, we also discuss why the EUV late phase is mainly observed in warm coronal emissions.

The relationship between dust polarization and extinction was determined for the cold dense starless molecular cloud core FeSt 1-457 based on the background star polarimetry of dichroic extinction at near-infrared wavelengths. Owing to the known (three-dimensional) magnetic field structure, the observed polarizations from the core were corrected by considering (a) the subtraction of the ambient polarization component, (b) the depolarization effect of inclined distorted magnetic fields, and (c) the magnetic inclination angle of the core. After these corrections, a linear relationship between polarization and extinction was obtained for the core in the range up to A_V ≈ 20 mag. The initial polarization versus extinction diagram changed dramatically after the corrections of (a) to (c), with the correlation coefficient being refined from 0.71 to 0.79. These corrections should affect the theoretical interpretation of the observational data. The slope of the finally obtained polarization–extinction relationship is ${P}_{H}/{E}_{H-{K}_{s}}=11.00\pm 0.72 \% \,{\mathrm{mag}}^{-1}$, which is close to the statistically estimated upper limit of the interstellar polarization efficiency. This consistency suggests that the upper limit of interstellar polarization efficiency might be determined by the observational viewing angle toward polarized astronomical objects.

The Magellanic Stream and Leading Arm of H i that stretches from the Large and Small Magellanic Clouds (LMC and SMC) and over 200° of the Southern sky is thought to be formed from multiple encounters between the LMC and SMC. In this scenario, most of the gas in the Stream and Leading Arm is stripped from the SMC, yet recent observations have shown a bifurcation of the Trailing Arm that reveals LMC origins for some of the gas. Absorption measurements in the Stream also reveal an order of magnitude more gas than in current tidal models. We present hydrodynamical simulations of the multiple encounters between the LMC and SMC at their first pass around the Milky Way, assuming that the Clouds were more extended and gas-rich in the past. Our models create filamentary structures of gas in the Trailing Stream from both the LMC and SMC. While the SMC trailing filament matches the observed Stream location, the LMC filament is offset. In addition, the total observed mass of the Stream in these models is underestimated by a factor of four when the ionized component is accounted for. Our results suggest that there should also be gas stripped from both the LMC and SMC in the Leading Arm, mirroring the bifurcation in the Trailing Stream. This prediction is consistent with recent measurements of spatial variation in chemical abundances in the Leading Arm, which show that gas from multiple sources is present, although its nature is still uncertain.

We investigate the properties of bright galaxies of various morphological types in Abell 1139 and Abell 2589, using pixel color–magnitude diagram (pCMD) analysis. The sample contains 32 galaxies brighter than M_r = −21.3 mag with spectroscopic redshifts, which are deeply imaged in the g and r bands using the MegaCam mounted on the Canada–France–Hawaii Telescope. After masking contaminants with two-step procedures, we examine how the detailed properties in the pCMDs depend on galaxy morphology and infrared color. The mean g − r color as a function of surface brightness (μ_r) in the pCMD of a galaxy shows good performance in distinguishing between early- and late-type galaxies, but it is not perfect because of the similarity between elliptical galaxies and bulge-dominated spiral galaxies. On the other hand, the g − r color dispersion as a function of μ_r works better. We find that the best set of parameters for galaxy classification is a combination of the minimum color dispersion at μ_r ≤ 21.2 mag arcsec⁻² and the maximum color dispersion at 20.0 ≤ μ_r ≤ 21.0 mag arcsec⁻²; the latter reflects the complexity of stellar populations at the disk component in a typical spiral galaxy. Finally, the color dispersion measurements of an elliptical galaxy appear to be correlated with the Wide-field Infrared Survey Explorer infrared color ([4.6]–[12]). This indicates that the complexity of stellar populations in an elliptical galaxy is related to its recent star formation activities. From this observational evidence, we infer that gas-rich minor mergers or gas interactions may have usually occurred during the recent growth of massive elliptical galaxies.

Here, we test the weighted horizontal magnetic gradient (WG_M) as a flare precursor, introduced by Korsós et al., by applying it to a magnetohydrodynamic (MHD) simulation of solar-like flares. The preflare evolution of the WG_M and the behavior of the distance parameter between the area-weighted barycenters of opposite-polarity sunspots at various heights is investigated in the simulated δ-type sunspot. Four flares emanated from this sunspot. We found the optimum heights above the photosphere where the flare precursors of the WG_M method are identifiable prior to each flare. These optimum heights agree reasonably well with the heights of the occurrence of flares identified from the analysis of their thermal and ohmic heating signatures in the simulation. We also estimated the expected time of the flare onsets from the duration of the approaching–receding motion of the barycenters of opposite polarities before each single flare. The estimated onset time and the actual time of occurrence of each flare are in good agreement at the corresponding optimum heights. This numerical experiment further supports the use of flare precursors based on the WG_M method.

We present a catalog of extended low surface brightness galaxies (LSBGs) identified in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first ∼200 deg² of the survey, we have uncovered 781 LSBGs, spanning red (g − i ≥ 0.64) and blue (g − i < 0.64) colors and a wide range of morphologies. Since we focus on extended galaxies (r_eff = 25–14''), our sample is likely dominated by low-redshift objects. We define LSBGs to have mean surface brightnesses ${\bar{\mu }}_{\mathrm{eff}}(g)\gt 24.3$ mag arcsec⁻², which allows nucleated galaxies into our sample. As a result, the central surface brightness distribution spans a wide range of μ₀(g) = 18–27.4 mag arcsec⁻², with 50% and 95% of galaxies fainter than 24.3 and 22 mag arcsec⁻², respectively. Furthermore, the surface brightness distribution is a strong function of color, with the red distribution being much broader and generally fainter than that of the blue LSBGs, and this trend shows a clear correlation with galaxy morphology. Red LSBGs typically have smooth light profiles that are well characterized by single-component Sérsic functions. In contrast, blue LSBGs tend to have irregular morphologies and show evidence for ongoing star formation. We cross-match our sample with existing optical, H i, and ultraviolet catalogs to gain insight into the physical nature of the LSBGs. We find that our sample is diverse, ranging from dwarf spheroidals and ultradiffuse galaxies in nearby groups to gas-rich irregulars to giant LSB spirals, demonstrating the potential of the HSC-SSP to provide a truly unprecedented view of the LSBG population.

We investigate nitrogen-isotope fractionation in forming and evolving molecular clouds using gas-ice astrochemical simulations. We find that the bulk gas can become depleted in heavy nitrogen (¹⁵N) due to the formation of ¹⁵N-enriched ices. Around the chemical transition from atomic nitrogen to N₂, N¹⁵N is selectively photodissociated, which results in the enrichment of ¹⁵N in atomic nitrogen. As ¹⁵N-enriched atomic nitrogen is converted to ammonia ice via grain surface reactions, the bulk gas is depleted in ¹⁵N. The level of ¹⁵N depletion in the bulk gas can be up to a factor of two compared to the elemental nitrogen-isotope ratio, depending on the photodesorption yield of ammonia ice. Once the nitrogen isotopes are differentially partitioned between gas and solids in a molecular cloud, the condition should remain in the later stages of star formation (e.g., prestellar core) as long as the sublimation of ammonia ice is inefficient. Our model suggests that all of the N-bearing molecules in the cold gas of star-forming regions can be depleted in ¹⁵N, which is at least qualitatively consistent with the observations toward prestellar core L1544. In our models, icy species show both ¹⁵N and deuterium fractionation. The fractionation pattern within ice mantles is different between ¹⁵N and deuterium, reflecting their fractionation mechanisms; while the concentration of deuterium almost monotonically increases from the lower layers of the ice mantles to the upper layers, the concentration of ¹⁵N reaches the maximum at a certain depth and declines toward the surface.

The stability of Earth's climate on geological timescales is enabled by the carbon–silicate cycle that acts as a negative feedback mechanism stabilizing surface temperatures via the intake and outgassing of atmospheric carbon. On Earth, this thermostat is enabled by plate tectonics that sequesters outgassed CO₂ back into the mantle via weathering and subduction at convergent margins. Here we propose a separate tectonic mechanism—vertical recycling—that can serve as the vehicle for CO₂ outgassing and sequestration over long timescales. The mechanism requires continuous tidal heating, which makes it particularly relevant to planets in the habitable zone of M stars. Dynamical models of this vertical recycling scenario and stability analysis show that temperate climates stable over timescales of billions of years are realized for a variety of initial conditions, even as the M star dims over time. The magnitude of equilibrium surface temperatures depends on the interplay of sea weathering and outgassing, which in turn depends on planetary carbon content, so that planets with lower carbon budgets are favored for temperate conditions. The habitability of planets such as found in the Trappist-1 system may be rooted in tidally driven tectonics.

The ESPERTA (Empirical model for Solar Proton Event Real Time Alert) forecast tool has a Probability of Detection (POD) of 63% for all >10 MeV events with proton peak intensity ≥10 pfu (i.e., ≥S1 events, S1 referring to minor storms on the NOAA Solar Radiation Storms scale), from 1995 to 2014 with a false alarm rate (FAR) of 38% and a median (minimum) warning time (WT) of ∼4.8 (0.4) hr. The NOAA space weather scale includes four additional categories: moderate (S2), strong (S3), severe (S4), and extreme (S5). As S1 events have only minor impacts on HF radio propagation in the polar regions, the effective threshold for significant space radiation effects appears to be the S2 level (100 pfu), above which both biological and space operation impacts are observed along with increased effects on HF propagation in the polar regions. We modified the ESPERTA model to predict ≥S2 events and obtained a POD of 75% (41/55) and an FAR of 24% (13/54) for the 1995–2014 interval with a median (minimum) WT of ∼1.7 (0.2) hr based on predictions made at the time of the S1 threshold crossing. The improved performance of ESPERTA for ≥S2 events is a reflection of the big flare syndrome, which postulates that the measures of the various manifestations of eruptive solar flares increase as one considers increasingly larger events.

Ion cyclotron waves (ICWs) can play important roles in the energization of plasma particles. Charged particle beams are ubiquitous in space, and astrophysical plasmas and can effectively lead to the generation of ICWs. Based on linear kinetic theory, we consider the excitation of ICWs by ion and electron beams in a compensated-current system. We also investigate the competition between reactive and kinetic instabilities. The results show that ion and electron beams both are capable of generating ICWs. For ICWs driven by ion beams, there is a critical beam velocity, v_bi^c, and critical wavenumber, k_z^c, for a fixed beam density; the reactive instability dominates the growth of ICWs when the ion-beam velocity ${v}_{{bi}}\gt {v}_{{bi}}^{c}$ and the wavenumber ${k}_{z}\lt {k}_{z}^{c}$, and the maximal growth rate is reached at ${k}_{z}\simeq 2{k}_{z}^{c}/3$ for a given ${v}_{{bi}}\gt {v}_{{bi}}^{c}$. For the slow ion beams with ${v}_{{bi}}\lt {v}_{{bi}}^{c}$, the kinetic instability can provide important growth rates of ICWs. On the other hand, ICWs driven by electron beams are excited only by the reactive instability, but require a critical velocity, ${v}_{{be}}^{c}\gg {v}_{{\rm{A}}}$ (the Alfvén velocity). In addition, the comparison between the approximate analytical results based on the kinetic theory and the exact numerical calculation based on the fluid model demonstrates that the reactive instabilities can well agree quantitatively with the numerical results by the fluid model. Finally, some possible applications of the present results to ICWs observed in the solar wind are briefly discussed.

We use a simple framework to calculate the time evolution of the composition of the fallback material onto a supermassive black hole arising from the tidal disruption of main-sequence stars. We study stars with masses between 0.8 and 3.0 M_⊙, at evolutionary stages from zero-age main sequence to terminal-age main sequence, built using the Modules for Experiments in Stellar Astrophysics code. We show that most stars develop enhancements in nitrogen (¹⁴N) and depletions in carbon (¹²C) and oxygen (¹⁶O) over their lifetimes, and that these features are more pronounced for higher mass stars. We find that, in an accretion-powered tidal disruption flare, these features become prominent only after the time of peak of the fallback rate and appear at earlier times for stars of increasing mass. We postulate that no severe compositional changes resulting from the fallback material should be expected near peak for a wide range of stellar masses and, as such, are unable to explain the extreme helium-to-hydrogen line ratios observed in some TDEs. On the other hand, the resulting compositional changes could help explain the presence of nitrogen-rich features, which are currently only detected after peak. When combined with the shape of the light curve, the time evolution of the composition of the fallback material provides a clear method to help constrain the nature of the disrupted star. This will enable a better characterization of the event by helping break the degeneracy between the mass of the star and the mass of the black hole when fitting tidal disruption light curves.

With growing data from ongoing and future supernova surveys, it is possible to empirically quantify the shapes of SNIa light curves in more detail, and to quantitatively relate the shape parameters with the intrinsic properties of SNIa. Building such relationships is critical in controlling systematic errors associated with supernova cosmology. Based on a collection of well-observed SNIa samples accumulated in the past years, we construct an empirical SNIa light curve model using a statistical method called the functional principal component analysis (FPCA) for sparse and irregularly sampled functional data. Using this method, the entire light curve of an SNIa is represented by a linear combination of principal component functions, and the SNIa is represented by a few numbers called "principal component scores." These scores are used to establish relations between light curve shapes and physical quantities such as intrinsic color, interstellar dust reddening, spectral line strength, and spectral classes. These relations allow for descriptions of some critical physical quantities based purely on light curve shape parameters. Our study shows that some important spectral feature information is being encoded in the broad band light curves; for instance, we find that the light curve shapes are correlated with the velocity and velocity gradient of the Si II λ6355 line. This is important for supernova surveys (e.g., LSST and WFIRST). Moreover, the FPCA light curve model is used to construct the entire light curve shape, which in turn is used in a functional linear form to adjust intrinsic luminosity when fitting distance models.

Recent theory predicts that first stars are born with a massive initial mass of ≳100 M_⊙. Pair-instability supernova (PISN) is a common fate for such massive stars. Our final goal is to prove the existence of PISNe and thus the high-mass nature of the initial mass function in the early universe by conducting abundance profiling, in which properties of a hypothetical first star is constrained by metal-poor star abundances. In order to determine reliable and useful abundances, we investigate the PISN nucleosynthesis taking both rotating and nonrotating progenitors for the first time. We show that the initial and CO core mass ranges for PISNe depend on the envelope structures: nonmagnetic rotating models developing inflated envelopes have a lower shifted CO mass range of ∼70–125 M_⊙, while nonrotating and magnetic rotating models with deflated envelopes have a range of ∼80–135 M_⊙. However, we find no significant difference in explosive yields from rotating and nonrotating progenitors, except for large nitrogen production in nonmagnetic rotating models. Furthermore, we conduct the first systematic comparison between theoretical yields and a large sample of metal-poor star abundances. We find that the predicted low [Na/Mg] ∼ −1.5 and high [Ca/Mg] ∼0.5–1.3 abundance ratios are the most important to discriminate PISN signatures from normal metal-poor star abundances, and confirm that no currently observed metal-poor star matches with the PISN abundance. An extensive discussion on the nondetection is presented.

High peaks in weak lensing (WL) maps originate dominantly from the lensing effects of single massive halos. Their abundance is therefore closely related to the halo mass function and thus a powerful cosmological probe. However, besides individual massive halos, large-scale structures (LSS) along lines of sight also contribute to the peak signals. In this paper, with ray-tracing simulations, we investigate the LSS projection effects. We show that for current surveys with a large shape noise, the stochastic LSS effects are subdominant. For future WL surveys with source galaxies having a median redshift z_med ∼ 1 or higher, however, they are significant. For the cosmological constraints derived from observed WL high-peak counts, severe biases can occur if the LSS effects are not taken into account properly. We extend the model of Fan et al. by incorporating the LSS projection effects into the theoretical considerations. By comparing with simulation results, we demonstrate the good performance of the improved model and its applicability in cosmological studies.

This study aims to provide further evidence for the potential influence of the global solar magnetic field on localized chromospheric jets, the macrospicules (MS). To find a connection between the long-term variation of properties of MS and other solar activity proxies, including, e.g., the temporal variation of the frequency shift of solar global oscillations, sunspot area, etc., a database overarching seven years of observations was compiled. This database contains 362 MS, based on observations at the 30.4 nm of the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory. Three of the five investigated physical properties of MS show a clear long-term temporal variation after smoothing the raw data. Wavelet analysis of the temporal variation of maximum length, maximum area, and average velocity is carried out. The results reveal a strong pattern of periodicities at around 2 years (also referred to as quasi-biennial oscillations—QBOs). A comparison with solar activity proxies that also possess the properties of QBOs provides some interesting features: the minima and maxima of QBOs of MS properties occur at around the same epoch as the minima and maxima of these activity proxies. For most of the time span investigated, the oscillations are out of phase. This out-of-phase behavior was also corroborated by a cross-correlation analysis. These results suggest that the physical processes that generate and drive the long-term evolution of the global solar activity proxies may be coupled to the short-term local physical processes driving the macrospicules, and, therefore modulate the properties of local dynamics.

Multiple stellar systems are ubiquitous in the Milky Way but are often unresolved and seen as single objects in spectroscopic, photometric, and astrometric surveys. However, modeling them is essential for developing a full understanding of large surveys such as Gaia and connecting them to stellar and Galactic models. In this paper, we address this problem by jointly fitting the Gaia and Two Micron All Sky Survey photometric and astrometric data using a data-driven Bayesian hierarchical model that includes populations of binary and trinary systems. This allows us to classify observations into singles, binaries, and trinaries, in a robust and efficient manner, without resorting to external models. We are able to identify multiple systems and, in some cases, make strong predictions for the properties of their unresolved stars. We will be able to compare such predictions with Gaia Data Release 4, which will contain astrometric identification and analysis of binary systems.

Many solar prominences end their lives in eruptions or abrupt disappearances that are associated with dynamical or thermal instabilities. Such instabilities are important because they may be responsible for energy transport and conversion. We present a clear observation of a streaming kink-mode Kelvin–Helmholtz Instability (KHI) taking place in a solar prominence using the Hα Lyot filter installed at the New Vacuum Solar Telescope, Fuxian-lake Solar Observatory in Yunnan, China. On one side of the prominence, a series of plasma blobs floated up from the chromosphere and streamed parallel to the limb. The plasma stream was accelerated to about 20–60 km s⁻1 and then undulated. We found that 2''- and 5''-size vortices formed, floated along the stream, and then broke up. After the 5''-size vortex, a plasma ejection out of the stream was detected in the Solar Dynamics Observatory/Atmospheric Imaging Assembly images. Just before the formation of the 5''-size vortex, the stream displayed an oscillatory transverse motion with a period of 255 s with the amplitude growing at the rate of 0.001 s⁻¹. We attribute this oscillation of the stream and the subsequent formation of the vortex to the KHI triggered by velocity shear between the stream, guided by the magnetic field and the surrounding media. The plasma ejection suggests the transport of prominence material into the upper layer by the KHI in its nonlinear stage.

Nuclear rings are excellent laboratories for probing diverse phenomena such as the formation and evolution of young massive star clusters and nuclear starbursts, as well as the secular evolution and dynamics of their host galaxies. We have compiled a sample of 17 galaxies with nuclear rings, which are well resolved by high-resolution Hubble and Spitzer Space Telescope imaging. For each nuclear ring, we identified the ring star cluster population, along with their physical properties (ages, masses, and extinction values). We also determined the integrated ring properties, including the average age, total stellar mass, and current star formation rate (SFR). We find that Sb-type galaxies tend to have the highest ring stellar mass fraction with respect to the host galaxy, and this parameter is correlated with the ring's SFR surface density. The ring SFRs are correlated with their stellar masses, which is reminiscent of the main sequence of star-forming galaxies. There are striking correlations between star-forming properties (i.e., SFR and SFR surface density) and nonaxisymmetric bar parameters, appearing to confirm previous inferences that strongly barred galaxies tend to have lower ring SFRs, although the ring star formation histories turn out to be significantly more complicated. Nuclear rings with higher stellar masses tend to be associated with lower cluster mass fractions, but there is no such relation for the ages of the rings. The two youngest nuclear rings in our sample, NGC 1512 and NGC 4314, which have the most extreme physical properties, represent the young extremity of the nuclear ring age distribution.

We present a search for coincidence between IceCube TeV neutrinos and fast radio bursts (FRBs). During the search period from 2010 May 31 to 2016 May 12, a total of 29 FRBs with 13 unique locations have been detected in the whole sky. An unbinned maximum likelihood method was used to search for spatial and temporal coincidence between neutrinos and FRBs in expanding time windows, in both the northern and southern hemispheres. No significant correlation was found in six years of IceCube data. Therefore, we set upper limits on neutrino fluence emitted by FRBs as a function of time window duration. We set the most stringent limit obtained to date on neutrino fluence from FRBs with an E⁻² energy spectrum assumed, which is 0.0021 GeV cm⁻² per burst for emission timescales up to ∼10² s from the northern hemisphere stacking search.

In cold dark-matter (CDM) cosmology, objects in the universe have grown under the effect of gravity of dark matter. The intracluster gas in a galaxy cluster was heated when the dark-matter halo formed through gravitational collapse. The potential energy of the gas was converted to thermal energy through this process. However, this process and the thermodynamic history of the gas have not been clearly characterized in connection with the formation and evolution of the internal structure of dark-matter halos. Here, we show that observational CLASH data of high-mass galaxy clusters lie on a plane in the three-dimensional logarithmic space of their characteristic radius r_s, mass M_s, and X-ray temperature T_X with a very small orthogonal scatter. The tight correlation indicates that the gas temperature was determined at a specific cluster formation time, which is encoded in r_s and M_s. The plane is tilted with respect to T_X ∝ M_s/r_s, which is the plane expected in the case of simplified virial equilibrium. We show that this tilt can be explained by a similarity solution, which indicates that clusters are not isolated but continuously growing through matter accretion from their outer environments. Numerical simulations reproduce the observed plane and its angle. This result holds independently of the gas physics implemented in the code, revealing the fundamental origin of this plane.

Oscillation properties are usually measured by fitting symmetric Lorentzian profiles to the power spectra of Sun-like stars. However, the line profiles of solar oscillations have been observed to be asymmetrical for the Sun. The physical origin of this line asymmetry is not fully understood; though, it should depend on the depth dependence of the source of wave excitation (convective turbulence) and details of the observable (velocity or intensity). For oscillations of the Sun, it has been shown that neglecting the asymmetry leads to systematic errors in the frequency determination. This could subsequently affect the results of seismic inferences of the solar internal structure. Using light curves from the Kepler spacecraft, we have measured mode asymmetries in 43 stars. We confirm that neglecting the asymmetry leads to systematic errors that can exceed the 1σ confidence intervals for seismic observations longer than one year. Therefore, the application of an asymmetric Lorentzian profile should be favored to improve the accuracy of the internal stellar structure and stellar fundamental parameters. We also show that the asymmetry changes sign between cool Sun-like stars and hotter stars. This provides the best constraints to date on the location of the excitation sources across the Hertzsprung–Russel diagram.

For the first time we demonstrate using a robust Bayesian approach to analyze the populations of radio-quiet (RQ) and radio-loud (RL) gamma-ray pulsars. We quantify their differences and obtain their distributions of the radio-cone opening half-angle δ and the magnetic inclination angle α by Bayesian inference. In contrast to the conventional frequentist point estimations that might be non-representative when the distribution is highly skewed or multi-modal, which is often the case when data points are scarce, Bayesian statistics displays the complete posterior distribution that the uncertainties can be readily obtained regardless of the skewness and modality. We found that the spin period, the magnetic field strength at the light cylinder, the spin-down power, the gamma-ray-to-X-ray flux ratio, and the spectral curvature significance of the two groups of pulsars exhibit significant differences at the 99% level. Using Bayesian inference, we are able to infer the values and uncertainties of δ and α from the distribution of RQ and RL pulsars. We found that δ is between 10° and 35° and the distribution of α is skewed toward large values.

We investigate the effects of AGN feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. The inner boundary of the simulation is chosen so that the Bondi radius is resolved. Compared to previous works, the two accretion modes—namely, hot and cold, which correspond to different accretion rates and have different radiation and wind outputs—are carefully discriminated, and the feedback effects by radiation and wind in each mode are taken into account. The most updated AGN physics, including the descriptions of radiation and wind from the hot accretion flows and wind from cold accretion disks, are adopted. Physical processes like star formation and SNe Ia and II are taken into account. We study the AGN light curve, typical AGN lifetime, growth of the black hole mass, AGN duty cycle, star formation, and X-ray surface brightness of the galaxy. We compare our simulation results with observations and find general consistency. Comparisons with previous simulation works find significant differences, indicating the importance of AGN physics. The respective roles of radiation and wind feedback are examined, and it is found that they are different for different problems of interest, such as AGN luminosity and star formation. We find that it is hard to neglect any of them, so we suggest using the names "cold feedback mode" and "hot feedback mode" to replace the currently used ones.

We investigate the evolution of star formation rates (SFRs), stellar masses, and M/L_{3.4 μm} ratios of brightest cluster galaxies (BCGs) in the COSMOS survey since z ∼ 1 to determine the contribution of star formation to the growth-rate of BCG stellar mass over time. Through the spectral energy density (SED) fitting of the GALEX, CFHT, Subaru, Vista, Spitzer, and Herschel photometric data available in the COSMOS2015 catalog, we estimate the stellar mass and SFR of each BCG. We use a modified version of the iSEDfit package to fit the SEDs of our sample with both stellar and dust emission models, as well as constrain the impact of star formation history assumptions on our results. We find that in our sample of COSMOS BCGs, star formation evolves similarly to that in BCGs in samples of more massive galaxy clusters. However, compared to the latter, the magnitude of star formation in our sample is lower by ∼1 dex. Additionally, we find an evolution of BCG baryonic mass-to-light ratio (M/L_{3.4 μm}) with redshift which is consistent with a passively aging stellar population. We use this to build upon Wen et al.'s low-redshift νL_{3.4 μm}–M_Stellar relation, quantifying a correlation between νL_{3.4 μm} and M_Stellar to z ∼ 1. By comparing our results to BCGs in Sunyaev–Zel'dovich and X-ray-selected samples of galaxy clusters, we find evidence that the normalization of star formation evolution in a cluster sample is driven by the mass range of the sample and may be biased upwards by cool cores.

We present the Neil Gehrels Swift Observatory (Swift), Fermi Large Area Telescope (Fermi-LAT), and Karl G. Jansky Very Large Array (VLA) observations of the γ-ray binary PSR J2032+4127/MT91 213, of which the periastron passage recently occurred in 2017 November. In the Swift X-ray light curve, the flux was steadily increasing before 2017 mid-October, however, a sharp X-ray dip on a weekly timescale is seen during the periastron passage, followed by a post-periastron X-ray flare lasting for ∼20 days. We suggest that the X-ray dip is caused by (i) an increase of the magnetization parameter at the shock, and (ii) the suppression due to the Doppler boosting effect. The 20-day post-periastron flare could be a consequence of the Be stellar disk passage by the pulsar. An orbital GeV modulation is also expected in our model; however, no significant variability is seen in the Fermi-LAT light curve. We suspect that the GeV emission from the interaction between the binary's members is hidden behind the bright magnetospheric emission of the pulsar. The pulsar gating technique would be useful to remove the magnetospheric emission and recover the predicted GeV modulation, if an accurate radio timing solution over the periastron passage is provided in the future.

The solar corona is frequently disrupted by coronal mass ejections (CMEs), whose core structure is believed to be a flux rope made of helical magnetic field. This has become a "standard" picture; though, it remains elusive how the flux rope forms and evolves toward eruption. While one-third of the ejecta passing through spacecraft demonstrate a flux-rope structure, the rest have complex magnetic fields. Are they originating from a coherent flux rope, too? Here we investigate the source region of a complex ejecta, focusing on a flare precursor with definitive signatures of magnetic reconnection, i.e., nonthermal electrons, flaring plasma, and bidirectional outflowing blobs. Aided by nonlinear force-free field modeling, we conclude that the reconnection occurs within a system of multiple braided flux ropes with different degrees of coherency. The observation signifies the importance of internal structure and dynamics in understanding CMEs and in predicting their impacts on Earth.

Vortex motions are frequently observed on the solar photosphere. These motions may play a key role in the transport of energy and momentum from the lower atmosphere into the upper solar atmosphere, contributing to coronal heating. The lower solar atmosphere also consists of complex networks of flux tubes that expand and merge throughout the chromosphere and upper atmosphere. We perform numerical simulations to investigate the behavior of vortex-driven waves propagating in a pair of such flux tubes in a non-force-free equilibrium with a realistically modeled solar atmosphere. The two flux tubes are independently perturbed at their footpoints by counter-rotating vortex motions. When the flux tubes merge, the vortex motions interact both linearly and nonlinearly. The linear interactions generate many small-scale transient magnetic substructures due to the magnetic stress imposed by the vortex motions. Thus, an initially monolithic tube is separated into a complex multithreaded tube due to the photospheric vortex motions. The wave interactions also drive a superposition that increases in amplitude until it exceeds the local Mach number and produces shocks that propagate upward with speeds of approximately 50 km s⁻¹. The shocks act as conduits transporting momentum and energy upward, and heating the local plasma by more than an order of magnitude, with a peak temperature of approximately 60,000 K. Therefore, we present a new mechanism for the generation of magnetic waveguides from the lower solar atmosphere to the solar corona. This wave guide appears as the result of interacting perturbations in neighboring flux tubes. Thus, the interactions of photospheric vortex motions is a potentially significant mechanism for energy transfer from the lower to upper solar atmosphere.

We examine the kinematics of the gas in the environments of galaxies hosting quasars at z ∼ 2. We employ 148 projected quasar pairs to study the circumgalactic gas of the foreground quasars in absorption. The sample selects foreground quasars with precise redshift measurements, using emission lines with precision ≲300 km s⁻¹ and average offsets from the systemic redshift $\lesssim | 100\,\mathrm{km}\,{{\rm{s}}}^{-1}|$. We stack the background quasar spectra at the foreground quasar's systemic redshift to study the mean absorption in C ii, C iv, and Mg ii. We find that the mean absorptions exhibit large velocity widths σ_v ≈ 300 km s⁻¹. Further, the mean absorptions appear to be asymmetric about the systemic redshifts. The mean absorption centroids exhibit small redshift relative to the systemic δv ≈ +200 km s⁻¹, with large intrinsic scatter in the centroid velocities of the individual absorption systems. We find the observed widths are consistent with gas in gravitational motion and Hubble flow. However, while the observation of large widths alone does not require galactic-scale outflows, the observed offsets suggest that the gas is on average outflowing from the galaxy. The observed offsets also suggest that the ionizing radiation from the foreground quasars is anisotropic and/or intermittent.

We use a 200 $\,{h}^{-1}\,\mathrm{Mpc}$ a-side N-body simulation to study the mass accretion history (MAH) of dark matter halos to be accreted by larger halos, which we call infall halos. We define a quantity ${a}_{\mathrm{nf}}\equiv (1+{z}_{{\rm{f}}})/(1+{z}_{\mathrm{peak}})$ to characterize the MAH of infall halos, where ${z}_{\mathrm{peak}}$ and ${z}_{{\rm{f}}}$ are the accretion and formation redshifts, respectively. We find that, at given ${z}_{\mathrm{peak}}$, their MAH is bimodal. Infall halos are dominated by a young population at high redshift and by an old population at low redshift. For the young population, the ${a}_{\mathrm{nf}}$ distribution is narrow and peaks at about 1.2, independent of ${z}_{\mathrm{peak}}$, while for the old population, the peak position and width of the ${a}_{\mathrm{nf}}$ distribution both increase with decreasing ${z}_{\mathrm{peak}}$ and are both larger than those of the young population. This bimodal distribution is found to be closely connected to the two phases in the MAHs of halos. While members of the young population are still in the fast accretion phase at z_peak, those of the old population have already entered the slow accretion phase at ${z}_{\mathrm{peak}}$. This bimodal distribution is not found for the whole halo population, nor is it seen in halo merger trees generated with the extended Press–Schechter formalism. The infall halo population at ${z}_{\mathrm{peak}}$ are, on average, younger than the whole halo population of similar masses identified at the same redshift. We discuss the implications of our findings in connection to the bimodal color distribution of observed galaxies and to the link between central and satellite galaxies.

The jet breaks in the afterglow light curves of short gamma-ray bursts (SGRBs), rarely detected so far, are crucial for estimating the half-opening angles of the ejecta (θ_j) and hence the neutron star merger rate. In this work, we report the detection of jet decline behaviors in GRB 150424A and GRB 160821B, and find θ_j ∼ 0.1 rad. Together with five events reported before 2015 and three others "identified" recently (GRB 050709, GRB 060614, and GRB 140903A), we have a sample consisting of nine SGRBs and one long-short GRB with reasonably estimated θ_j. In particular, three Swift bursts in the sample have redshifts z ≤ 0.2, with which we estimate the local neutron star merger rate density to be $\sim {1109}_{-657}^{+1432}\,{\mathrm{Gpc}}^{-3}\,{\mathrm{yr}}^{-1}$ or ${162}_{-83}^{+140}\,{\mathrm{Gpc}}^{-3}\,{\mathrm{yr}}^{-1}$ if the narrowly beamed GRB 061201 is excluded. Inspired by the typical θ_j ∼ 0.1 rad found currently, we further investigate whether the off-beam GRBs (in the uniform jet model) or the off-axis events (in the structured jet model) can significantly enhance the GRB/GW association. For the former, the enhancement is at most moderate, while for the latter the enhancement can be much greater and a high GRB/GW association probability of ∼10% is possible. We also show that the data of GRB 160821B may contain a macronova/kilonova emission component with a temperature of ∼3100 K at ∼3.6 days after the burst and more data are needed to ultimately clarify.

Intense heating in the quiet-Sun chromosphere raises the temperature from 4000 to 6500 K but, despite decades of study, the underlying mechanism remains a mystery. This study continues to explore the possibility that the Farley–Buneman instability contributes to chromospheric heating. This instability occurs in weakly ionized collisional plasmas in which electrons are magnetized, but ions are not. A mixture of metal ions generate the plasma density in the coolest parts of the chromosphere; while some ions are weakly magnetized, others are demagnetized by neutral collisions. This paper incorporates the effects of multiple, arbitrarily magnetized species of ions to the theory of the Farley–Buneman instability and examines the ramifications on instability in the chromosphere. The inclusion of magnetized ions introduces new restrictions on the regions in which the instability can occur in the chromosphere—in fact, it confines the instability to the regions in which heating is observed. For a magnetic field of 30 G, the minimum ambient electric field capable of driving the instability is 13.5 V/m at the temperature minimum.

We used our newly developed magnetohydrodynamic (MHD) code to perform 2.5D simulations of a fast-mode MHD wave interacting with coronal holes (CHs) of varying Alfvén speed that result from assuming different CH densities. We find that this interaction leads to effects like reflection, transmission, stationary fronts at the CH boundary, and the formation of a density depletion that moves in the opposite direction to the incoming wave. We compare these effects with regard to the different CH densities and present a comprehensive analysis of morphology and kinematics of the associated secondary waves. We find that the density value inside the CH influences the phase speed and the amplitude values of density and magnetic field for all different secondary waves. Moreover, we observe a correlation between the CH density and the peak values of the stationary fronts at the CH boundary. The findings of reflection and transmission on the one hand and the formation of stationary fronts caused by the interaction of MHD waves with CHs on the other hand strongly support the theory that large-scale disturbances in the corona are fast-mode MHD waves.

We present timing solutions for 10 pulsars discovered in 350 MHz searches with the Green Bank Telescope. Nine of these were discovered in the Green Bank Northern Celestial Cap survey and one was discovered by students in the Pulsar Search Collaboratory program during an analysis of drift-scan data. Following the discovery and confirmation with the Green Bank Telescope, timing has yielded phase-connected solutions with high-precision measurements of rotational and astrometric parameters. Eight of the pulsars are slow and isolated, including PSR J0930−2301, a pulsar with a nulling fraction lower limit of ∼30% and a nulling timescale of seconds to minutes. This pulsar also shows evidence of mode changing. The remaining two pulsars have undergone recycling, accreting material from binary companions, resulting in higher spin frequencies. PSR J0557−2948 is an isolated, 44 ms pulsar that has been partially recycled and is likely a former member of a binary system that was disrupted by a second supernova. The paucity of such so-called "disrupted binary pulsars" (DRPs) compared to double neutron star (DNS) binaries can be used to test current evolutionary scenarios, especially the kicks imparted on the neutron stars in the second supernova. There is some evidence that DRPs have larger space velocities, which could explain their small numbers. PSR J1806+2819 is a 15 ms pulsar in a 44-day orbit with a low-mass white dwarf companion. We did not detect the companion in archival optical data, indicating that it must be older than 1200 Myr.

NGC 796 is a massive young cluster located 59 kpc from us in the diffuse intergalactic medium of the 1/5–1/10 Z_☉ Magellanic Bridge, allowing us to probe variations in star formation and stellar evolution processes as a function of metallicity in a resolved fashion, and providing a link between resolved studies of nearby solar-metallicity and unresolved distant metal-poor clusters located in high-redshift galaxies. In this paper, we present adaptive optics griHα imaging of NGC 796 (at 05, which is ∼0.14 pc at the cluster distance) along with optical spectroscopy of two bright members to quantify the cluster properties. Our aim is to explore whether star formation and stellar evolution vary as a function of metallicity by comparing the properties of NGC 796 to higher-metallicity clusters. We find an age of ${20}_{-5}^{+12}$ Myr from isochronal fitting of the cluster main sequence in the color–magnitude diagram. Based on the cluster luminosity function, we derive a top-heavy stellar initial mass function (IMF) with a slope α = 1.99 ± 0.2, hinting at a metallicity and/or environmental dependence of the IMF, which may lead to a top-heavy IMF in the early universe. Study of the Hα emission-line stars reveals that classical Be stars constitute a higher fraction of the total B-type stars when compared with similar clusters at greater metallicity, providing some support to the chemically homogeneous theory of stellar evolution. Overall, NGC 796 has a total estimated mass of 990 ± 200 M_⊙, and a core radius of 1.4 ± 0.3 pc, which classifies it as a massive young open cluster, unique in the diffuse interstellar medium of the Magellanic Bridge.

As part of a wider search for radio emission from nearby systems known or suspected to contain extrasolar planets, Eridani was observed by the Jansky Very Large Array (VLA) in the 2–4 GHz and 4–8 GHz frequency bands. In addition, as part of a separate survey of thermal emission from solar-like stars, Eri was observed in the 8–12 GHz and the 12–18 GHz bands of the VLA. Quasi-steady continuum radio emission from Eri was detected in the three high-frequency bands at levels ranging from 67 to 83 μJy. No significant variability is seen in the quasi-steady emission. The emission in the 2–4 GHz emission, however, is shown to be the result of a circularly polarized (up to 50%) radio pulse or flare of a few minutes in duration that occurred at the beginning of the observation. We consider the astrometric position of the radio source in each frequency band relative to the expected position of the K2V star and the purported planet. The quasi-steady radio emission at frequencies ≥8 GHz is consistent with a stellar origin. The quality of the 4–8 GHz astrometry provides no meaningful constraint on the origin of the emission. The location of the 2–4 GHz radio pulse is >2.5σ from the star; however, based on the ephemeris of Benedict et al., it is not consistent with the expected location of the planet either. If the radio pulse has a planetary origin, then either the planetary ephemeris is incorrect or the emission originates from another planet.

The evolution of the remnant of the merger of two white dwarfs is still an open problem. Furthermore, few studies have addressed the case in which the remnant is a magnetic white dwarf with a mass larger than the Chandrasekhar limiting mass. Angular momentum losses might bring the remnant of the merger to the physical conditions suitable for developing a thermonuclear explosion. Alternatively, the remnant may be prone to gravitational or rotational instabilities, depending on the initial conditions reached after the coalescence. Dipole magnetic braking is one of the mechanisms that can drive such losses of angular momentum. However, the timescale on which these losses occur depends on several parameters, like the strength of the magnetic field. In addition, the coalescence leaves a surrounding Keplerian disk that can be accreted by the newly formed white dwarf. Here we compute the post-merger evolution of a super-Chandrasekhar magnetized white dwarf taking into account all the relevant physical processes. These include magnetic torques acting on the star, accretion from the Keplerian disk, the threading of the magnetic field lines through the disk, and the thermal evolution of the white dwarf core. We find that the central remnant can reach the conditions suitable to develop a thermonuclear explosion before other instabilities (such as the inverse beta-decay instability or the secular axisymmetric instability) are reached, which would instead lead to gravitational collapse of the magnetized remnant.

Spectroscopic observations of some metal-rich white dwarfs (WDs), believed to be polluted by planetary material, reveal the presence of compact gaseous metallic disks orbiting them. The observed variability of asymmetric, double-peaked emission-line profiles in about half of such systems could be interpreted as the signature of precession of an eccentric gaseous debris disk. The variability timescales—from decades down to 1.4 year (recently inferred for the debris disk around HE 1349–2305)—are in rough agreement with the rate of general relativistic (GR) precession in the test-particle limit. However, it has not been demonstrated that this mechanism can drive such a fast, coherent precession of a radially extended (out to $1\,{R}_{\odot }$) gaseous disk mediated by internal stresses (pressure). Here, we use the linear theory of eccentricity evolution in hydrodynamic disks to determine several key properties of eccentric modes in gaseous debris disks around WDs. We find a critical dependence of both the precession period and radial eccentricity distribution of the modes on the inner disk radius, r_in. For small inner radii, ${r}_{\mathrm{in}}\lesssim (0.2\mbox{--}0.4)\,{R}_{\odot }$, the modes are GR-driven, with periods of ≈1–10 year. For ${r}_{\mathrm{in}}\gtrsim (0.2\mbox{--}0.4)\,{R}_{\odot }$, the modes are pressure dominated, with periods of ≈3–20 year. Correspondence between the variability periods and inferred inner radii of the observed disks is in general agreement with this trend. In particular, the short period of HE 1349–2305 is consistent with its small r_in. Circum-WD debris disks may thus serve as natural laboratories for studying the evolution of eccentric gaseous disks.

The power spectrum of magnetic fluctuations has a break at the high-frequency end of the inertial range. Beyond this break, the spectrum becomes steeper than the Kolmogorov law f^−5/3. The break frequency was found to be associated with plasma beta (β). However, the full-range β dependence of the ion-scale spectral break has not been presented before in observational studies. Here we show the continuous variation of the break frequency on full-range β in the solar wind turbulence. By using measurements from the WIND and Ulysses spacecraft, we show the break frequency (f_b) normalized, respectively, by the frequencies corresponding to ion inertial length (f_di), ion gyroradius (${f}_{\rho i}$), and cyclotron resonance scale (f_ri) as a function of β for 1306 intervals. Their β values spread from 0.005 to 20, which nearly covers the full β range of the observed solar wind turbulence. It is found that ${f}_{b}/{f}_{{di}}$ (${f}_{b}/{f}_{\rho i}$) generally decreases (increases) with β, while ${f}_{b}/{f}_{{ri}}$ is nearly a constant. We perform a linear fit on the statistical result, and obtain the empirical formulas ${f}_{b}/{f}_{{di}}\sim {\beta }^{-1/4}$, ${f}_{b}/{f}_{\rho i}\sim {\beta }^{1/4}$, and ${f}_{b}/{f}_{{ri}}\sim 0.90$ to describe the relation between f_b and β. We also compare our observations with a numerical simulation and the prediction by ion cyclotron resonance theory. Our result favors the idea that the cyclotron resonance is an important mechanism for energy dissipation at the spectral break. When β ≪ 1 and β ≫ 1, the break at f_di and ${f}_{\rho i}$ may also be associated with other processes.

Using the observations recorded by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory and the Interface Region Imaging Spectrograph (IRIS) and the Extreme-ultraviolet Imaging Spectrometer and X-Ray Telescope both on board Hinode, we present evidence of chromospheric evaporation in a coronal loop after the occurrence of two active region transient brightenings (ARTBs) at the two footpoints. The chromospheric evaporation started nearly simultaneously in all of the three hot channels of AIA 131, 94, and 335 Å and was observed to be temperature dependent, being fastest in the highest temperature channel. The whole loop became fully brightened following the ARTBs after ≈25 s in 131 Å, ≈40 s in 94 Å, and ≈6.5 minutes in 335 Å. The differential emission measurements at the two footpoints (i.e., of two ARTBs) and at the loop top suggest that the plasma attained a maximum temperature of ∼10 MK at all these locations. The spectroscopic observations from IRIS revealed the presence of redshifted emission of ∼20 km s⁻¹ in cooler lines like C ii and Si iv during the ARTBs that was cotemporal with the evaporation flow at the footpoint of the loop. During the ARTBs, the line width of C ii and Si iv increased nearly by a factor of two during the peak emission. Moreover, enhancement in the line width preceded that in the Doppler shift, which again preceded enhancement in the intensity. The observed results were qualitatively reproduced by 1D hydrodynamic simulations, where energy was deposited at both of the footpoints of a monolithic coronal loop that mimicked the ARTBs identified in the observations.

We present a model for the evolution of supermassive protostars from their formation at ${M}_{\star }\simeq 0.1\,{M}_{\odot }$ until their growth to ${M}_{\star }\simeq {10}^{5}\,{M}_{\odot }$. To calculate the initial properties of the object in the optically thick regime, we follow two approaches: one based on idealized thermodynamic considerations, and another based on a more detailed one-zone model. Both methods derive a similar value of ${n}_{{\rm{F}}}\simeq 2\times {10}^{17}\,{\mathrm{cm}}^{-3}$ for the density of the object when opacity becomes important, i.e., the opacity limit. The subsequent evolution of the growing protostar is determined by the accretion of gas onto the object and can be described by a mass–radius relation of the form ${R}_{\star }\propto {M}_{\star }^{1/3}$ during the early stages, and of the form ${R}_{\star }\propto {M}_{\star }^{1/2}$ when internal luminosity becomes important. For the case of a supermassive protostar, this implies that the radius of the star grows from ${R}_{\star }\simeq 0.65\,\mathrm{au}$ to ${R}_{\star }\simeq 250\,\mathrm{au}$ during its evolution. Finally, we use this model to construct a subgrid recipe for accreting sink particles in numerical simulations. A prime ingredient thereof is a physically motivated prescription for the accretion radius and the effective temperature of the growing protostar embedded inside it. From the latter, we can conclude that photoionization feedback can be neglected until very late in the assembly process of the supermassive object.

The Fourier technique of extracting rotation rates and macroturbulence-velocity dispersions from the shapes and broadening of stellar spectral lines is applied to the solar-flux spectrum. Lines with equivalent widths less than ∼0.055 Å are shown to have the advantage over stronger lines by allowing the residual transform to be followed to higher frequencies. The standard radial-tangential macroturbulence formulation fits the observations well and yields an equatorial velocity that is within a few percent of the correct rate.

Electrons accelerated in relativistic collisionless shocks are usually assumed to follow a power-law energy distribution with an index of p. Observationally, although most gamma-ray bursts (GRBs) have afterglows that are consistent with p > 2, there are still a few GRBs suggestive of a hard (p < 2) electron energy spectrum. Our previous work showed that GRB 091127 gave strong evidence for a double power-law hard electron energy (DPLH) spectrum with 1 < p₁ < 2, p₂ > 2 and an "injection break" assumed as γ_b ∝ γ^q in the highly relativistic regime, where γ is the bulk Lorentz factor of the jet. In this paper, we show that GRB 060614 and GRB 060908 provide further evidence for such a DPLH spectrum. We interpret the multiband afterglow of GRB 060614 with the DPLH model in a homogeneous interstellar medium by taking into account a continuous energy injection process, while, for GRB 060908, a wind-like circumburst density profile is used. The two bursts, along with GRB 091127, suggest a similar behavior in the evolution of the injection break, with q ∼ 0.5. Whether this represents a universal law of the injection break remains uncertain and more afterglow observations such as these are needed to test this conjecture.

We present a comprehensive analysis of 21 light curves of Type 1 active galactic nuclei (AGN) from the Kepler spacecraft. First, we describe the necessity and development of a customized pipeline for treating Kepler data of stochastically variable sources like AGN. We then present the light curves, power spectral density functions (PSDs), and flux histograms. The light curves display an astonishing variety of behaviors, many of which would not be detected in ground-based studies, including switching between distinct flux levels. Six objects exhibit PSD flattening at characteristic timescales that roughly correlate with black hole mass. These timescales are consistent with orbital timescales or free-fall accretion timescales. We check for correlations of variability and high-frequency PSD slope with accretion rate, black hole mass, redshift, and luminosity. We find that bolometric luminosity is anticorrelated with both variability and steepness of the PSD slope. We do not find evidence of the linear rms–flux relationships or lognormal flux distributions found in X-ray AGN light curves, indicating that reprocessing is not a significant contributor to optical variability at the 0.1%–10% level.

We perform numerical simulations of the TRAPPIST-1 system of seven exoplanets orbiting a nearby M dwarf, starting with a previously suggested stable configuration. The long-term stability of this configuration is confirmed, but the motion of planets is found to be chaotic. The eccentricity values are found to vary within finite ranges. The rates of tidal dissipation and tidal evolution of orbits are estimated, assuming an Earth-like rheology for the planets. We find that under this assumption, the planets b, d, and e were captured in the 3:2 or higher spin–orbit resonances during the initial spin-down, but slipped further down into the 1:1 resonance. Depending on its rheology, the innermost planet b may be captured in a stable pseudosynchronous rotation. Nonsynchronous rotation ensures higher levels of tidal dissipation and internal heating. The positive feedback between the viscosity and the dissipation rate—and the ensuing runaway heating—are terminated by a few self-regulation processes. When the temperature is high and the viscosity is low enough, the planet spontaneously leaves the 3:2 resonance. Further heating is stopped either by passing the peak dissipation or by the emergence of partial melt in the mantle. In the post-solidus state, the tidal dissipation is limited to the levels supported by the heat transfer efficiency. The tides on the host star are unlikely to have had a significant dynamical impact. The tides on the synchronized inner planets tend to reduce these planets' orbital eccentricity, possibly contributing thereby to the system's stability.

We apply a newly developed on-the-fly mosaicing technique on the Jansky Very Large Array (VLA) at 3 GHz in order to carry out a sensitive search for an afterglow from the Advanced LIGO binary black hole merger event GW151226. In three epochs between 1.5 and 6 months post-merger, we observed a 100 deg² region, with more than 80% of the survey region having an rms sensitivity of better than 150 μJy/beam, in the northern hemisphere with a merger containment probability of 10%. The data were processed in near real time and analyzed to search for transients and variables. No transients were found but we have demonstrated the ability to conduct blind searches in a time-frequency phase space where the predicted afterglow signals are strongest. If the gravitational wave event is contained within our survey region, the upper limit on any late-time radio afterglow from the merger event at an assumed mean distance of 440 Mpc is about 10²⁹ erg s⁻¹ Hz⁻¹. Approximately 1.5% of the radio sources in the field showed variability at a level of 30%, and can be attributed to normal activity from active galactic nuclei. The low rate of false positives in the radio sky suggests that wide-field imaging searches at a few Gigahertz can be an efficient and competitive search strategy. We discuss our search method in the context of the recent afterglow detection from GW170817 and radio follow-up in future gravitational wave observing runs.

We study tidal features around galaxies in the REsolved Spectroscopy Of a Local VolumE (RESOLVE) survey. Our sample consists of 1048 RESOLVE galaxies that overlap with the DECam Legacy Survey, which reaches an r-band 3σ depth of ∼27.9 mag arcsec⁻² for a 100 arcsec² feature. Images were masked, smoothed, and inspected for tidal features such as streams, shells, or tails/arms. We find tidal features in 17^±2% of our galaxies, setting a lower limit on the true frequency. The frequency of tidal features in the gas-poor (gas-to-stellar mass ratio <0.1) subsample is lower than in the gas-rich subsample (13^±3% versus 19^±2%). Within the gas-poor subsample, galaxies with tidal features have higher stellar and halo masses, ∼3× closer distances to nearest neighbors (in the same group), and possibly fewer group members at fixed halo mass than galaxies without tidal features, but similar specific star formation rates. These results suggest tidal features in gas-poor galaxies are typically streams/shells from dry mergers or satellite disruption. In contrast, the presence of tidal features around gas-rich galaxies does not correlate with stellar or halo mass, suggesting these tidal features are often tails/arms from resonant interactions. Similar to tidal features in gas-poor galaxies, tidal features in gas-rich galaxies imply 1.7× closer nearest neighbors in the same group; however, they are associated with diskier morphologies, higher star formation rates, and higher gas content. In addition to interactions with known neighbors, we suggest that tidal features in gas-rich galaxies may arise from accretion of cosmic gas and/or gas-rich satellites below the survey limit.

We present Magellan/IMACS, Anglo-Australian Telescope/AAOmega+2dF, and Very Large Telescope/GIRAFFE+FLAMES spectroscopy of the Carina II (Car II) and Carina III (Car III) dwarf galaxy candidates, recently discovered in the Magellanic Satellites Survey (MagLiteS). We identify 18 member stars in Car II, including two binaries with variable radial velocities and two RR Lyrae stars. The other 14 members have a mean heliocentric velocity ${v}_{\mathrm{hel}}=477.2\pm 1.2$$\mathrm{km}\,{{\rm{s}}}^{-1}$ and a velocity dispersion of ${\sigma }_{v}={3.4}_{-0.8}^{+1.2}$$\mathrm{km}\,{{\rm{s}}}^{-1}$. Assuming Car II is in dynamical equilibrium, we derive a total mass within the half-light radius of ${1.0}_{-0.4}^{+0.8}\times {10}^{6}$${M}_{\odot }$, indicating a mass-to-light ratio of ${369}_{-161}^{+309}$${M}_{\odot }$/${L}_{\odot }$. From equivalent width measurements of the calcium triplet lines of nine red giant branch (RGB) stars, we derive a mean metallicity of ${\rm{[Fe/H]}}=-2.44\pm 0.09$ with dispersion ${\sigma }_{{\rm{[Fe/H]}}}={0.22}_{-0.07}^{+0.10}$. Considering both the kinematic and chemical properties, we conclude that Car II is a dark-matter-dominated dwarf galaxy. For Car III, we identify four member stars, from which we calculate a systemic velocity of ${v}_{\mathrm{hel}}={284.6}_{-3.1}^{+3.4}$$\mathrm{km}\,{{\rm{s}}}^{-1}$. The brightest RGB member of Car III has a metallicity of ${\rm{[Fe/H]}}\,=-1.97\pm 0.12$. Due to the small size of the Car III spectroscopic sample, we cannot conclusively determine its nature. Although these two systems have the smallest known physical separation (${\rm{\Delta }}d\sim 10\,\mathrm{kpc}$) among Local Group satellites, the large difference in their systemic velocities, $\sim 200\,\mathrm{km}\,{{\rm{s}}}^{-1}$, indicates that they are unlikely to be a bound pair. One or both systems are likely associated with the Large Magellanic Cloud (LMC), and may remain LMC satellites today. No statistically significant excess of γ-ray emission is found at the locations of Car II and Car III in eight years of Fermi-LAT data.