In this paper we present the differential emission measures (DEMs) of two sub-A class microflares observed in hard X-rays (HXRs) by the FOXSI-2 sounding rocket experiment, on 2014 December 11. The second FOXSI (Focusing Optics X-ray Solar Imager) flight was coordinated with instruments Hinode/XRT and SDO/AIA, which provided observations in soft X-rays (SXR) and Extreme Ultraviolet (EUV). This unique dataset offers an unprecedented temperature coverage useful for characterizing the plasma temperature distribution of microflares. By combining data from FOXSI-2, XRT, and AIA, we determined a well-constrained DEM for the microflares. The resulting DEMs peak around 3MK and extend beyond 10MK. The emission measures determined from FOXSI-2 were lower than 10 26cm-5 for temperatures higher than 5MK; faint emission in this range is best measured in HXRs. The coordinated FOXSI-2 observations produce one of the few definitive measurements of the distribution and the amount of plasma above 5MK in microflares. We utilize the multi-thermal DEMs to calculate the amount of thermal energy released during both the microflares as ~ 5.0 x 10 28 ergs for Microflare 1 and ~ 1.6 x 10 28 ergs for Microflare 2. We also show the multi-thermal DEMs provide a more comprehensive thermal energy estimates than isothermal approximation, which systematically underestimates the amount of thermal energy released.

We have obtained Hubble Space Telescope STIS and NICMOS, and Gemini/GPI scattered light images of the HD 191089 debris disk. We identify two spatial components: a ring resembling Kuiper Belt in radial extent (FWHM: ${\sim}$25 au, centered at ${\sim}$46 au), and a halo extending to ${\sim}$640 au. We find that the halo is significantly bluer than the ring, consistent with the scenario that the ring serves as the "birth ring" for the smaller dust in the halo. We measure the scattering phase functions in the 30{\deg}-150{\deg} scattering angle range and find the halo dust is both more forward- and backward-scattering than the ring dust. We measure a surface density power law index of -0.68${\pm}$0.04 for the halo, which indicates the slow-down of the radial outward motion of the dust. Using radiative transfer modeling, we attempt to simultaneously reproduce the (visible) total and (near-infrared) polarized intensity images of the birth ring. Our modeling leads to mutually inconsistent results, indicating that more complex models, such as the inclusion of more realistic aggregate particles, are needed.

We present new high resolution imaging of a light-scattering dust ring and halo around the young star HD 35841. Using spectroscopic and polarimetric data from the Gemini Planet Imager in H-band (1.6 microns), we detect the highly inclined (i=85 deg) ring of debris down to a projected separation of ~12 au (~0.12") for the first time. Optical imaging from HST/STIS shows a smooth dust halo extending outward from the ring to >140 au (>1.4"). We measure the ring's scattering phase function and polarization fraction over scattering angles of 22-125 deg, showing a preference for forward scattering and a polarization fraction that peaks at ~30% near the ansae. Modeling of the scattered-light disk indicates that the ring spans radii of ~60-220 au, has a vertical thickness similar to that of other resolved dust rings, and contains grains as small as 1.5 microns in diameter. These models also suggest the grains have a low porosity, are more likely to consist of carbon than astrosilicates, and contain significant water ice. The halo has a surface brightness profile consistent with that expected from grains pushed by radiation pressure from the main ring onto highly eccentric but still bound orbits. We also briefly investigate arrangements of a possible inner disk component implied by our spectral energy distribution models, and speculate about the limitations of Mie theory for doing detailed analyses of debris disk dust populations.

At a distance of 2~pc, our nearest brown dwarf neighbor, Luhman 16 AB, has been extensively studied since its discovery 3 years ago, yet its most fundamental parameter -- the masses of the individual dwarfs -- has not been constrained with precision. In this work we present the full astrometric orbit and barycentric motion of Luhman 16 AB and the first precision measurements of the individual component masses. We draw upon archival observations spanning 31 years from the European Southern Observatory (ESO) Schmidt Telescope, the Deep Near-Infrared Survey of the Southern Sky (DENIS), public FORS2 data on the Very Large Telescope (VLT), and new astrometry from the Gemini South Multiconjugate Adaptive Optics System (GeMS). Finally, we include three radial velocity measurements of the two components from VLT/CRIRES, spanning one year. With this new data sampling a full period of the orbit, we use a Markov Chain Monte Carlo algorithm to fit a 16-parameter model incorporating mutual orbit and barycentric motion parameters and constrain the individual masses to be~$27.9^{+1.1}_{-1.0}$~$M_{J}$ for the T dwarf and~$34.2^{+1.3}_{-1.1}$~$M_{J}$ for the L dwarf. Our measurements of Luhman 16 AB's mass ratio and barycentric motion parameters are consistent with previous estimates in the literature utilizing recent astrometry only. The GeMS-derived measurements of the Luhman 16 AB separation in 2014-2015 agree closely with Hubble Space Telescope (HST) measurements made during the same epoch Bedin et al. 2017, and the derived mutual orbit agrees with those measurements to within the HST uncertainties of $0.3 - 0.4$ milliarcseconds.

We describe a Bayesian rejection sampling algorithm designed to efficiently compute posterior distributions of orbital elements for data covering short fractions of long-period exoplanet orbits. Our implementation of this method, Orbits for the Impatient (OFTI), converges up to several orders of magnitude faster than two implementations of MCMC in this regime. We illustrate the efficiency of our approach by showing that OFTI calculates accurate posteriors for all existing astrometry of the exoplanet 51 Eri b up to 100 times faster than a Metropolis-Hastings MCMC. We demonstrate the accuracy of OFTI by comparing our results for several orbiting systems with those of various MCMC implementations, finding the output posteriors to be identical within shot noise. We also describe how our algorithm was used to successfully predict the location of 51 Eri b six months in the future based on less than three months of astrometry. Finally, we apply OFTI to ten long-period exoplanets and brown dwarfs, all but one of which have been monitored over less than 3% of their orbits, producing fits to their orbits from astrometric records in the literature.

In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles ($E_{\mathrm{nt}}$), the energy of direct heating ($E_{\mathrm{dir}}$), and the energy in coronal mass ejections ($E_{\mathrm{CME}}$), which are the primary energy dissipation processes in a flare, is found to have a ratio of $(E_{\mathrm{nt}}+E_{\mathrm{dir}}+ E_{\mathrm{CME}})/E_{\mathrm{mag}} = 0.87 \pm 0.18$, compared with the dissipated magnetic free energy $E_{\mathrm{mag}}$, which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: $0.51\pm0.17$ in nonthermal energy of $\ge 6$ keV electrons, $0.17\pm0.17$ in nonthermal $\ge 1$ MeV ions, $0.07\pm0.14$ in CMEs, and $0.07\pm0.17$ in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction $\approx 0.03$ of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at $e_c \approx 6$ keV, based on the mean differential emission measure (DEM) peak temperature of $T_e=8.6$ MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events.

How impulsive magnetic energy release leads to solar eruptions and how those eruptions are energized and evolve are vital unsolved problems in Heliophysics. The standard model for solar eruptions summarizes our current understanding of these events. Magnetic energy in the corona is released through drastic restructuring of the magnetic field via reconnection. Electrons and ions are then accelerated by poorly understood processes. Theories include contracting loops, merging magnetic islands, stochastic acceleration, and turbulence at shocks, among others. Although this basic model is well established, the fundamental physics is poorly understood. HXR observations using grazing-incidence focusing optics can now probe all of the key regions of the standard model. These include two above-the-looptop (ALT) sources which bookend the reconnection region and are likely the sites of particle acceleration and direct heating. The science achievable by a direct HXR imaging instrument can be summarized by the following science questions and objectives which are some of the most outstanding issues in solar physics (1) How are particles accelerated at the Sun? (1a) Where are electrons accelerated and on what time scales? (1b) What fraction of electrons is accelerated out of the ambient medium? (2) How does magnetic energy release on the Sun lead to flares and eruptions? A Focusing Optics X-ray Solar Imager (FOXSI) instrument, which can be built now using proven technology and at modest cost, would enable revolutionary advancements in our understanding of impulsive magnetic energy release and particle acceleration, a process which is known to occur at the Sun but also throughout the Universe.

Quasi-periodic pulsations (QPP) are often observed in X-ray emission from solar flares. To date, it is unclear what their physical origins are. Here, we present a multi-instrument investigation of the nature of QPP during the impulsive and decay phases of the X1.0 flare of 28 October 2013. We focus on the character of the fine structure pulsations evident in the soft X-ray time derivatives and compare this variability with structure across multiple wavelengths including hard X-ray and microwave emission. We find that during the impulsive phase of the flare, high correlations between pulsations in the thermal and non-thermal emissions are seen. A characteristic timescale of ~20s is observed in all channels and a second timescale of ~55s is observed in the non-thermal emissions. Soft X-ray pulsations are seen to persist into the decay phase of this flare, up to 20 minutes after the non-thermal emission has ceased. We find that these decay phase thermal pulsations have very small amplitude and show an increase in characteristic timescale from ~40s up to ~70s. We interpret the bursty nature of the co-existing multi-wavelength QPP during the impulsive phase in terms of episodic particle acceleration and plasma heating. The persistent thermal decay phase QPP are most likely connected with compressive MHD processes in the post-flare loops such as the fast sausage mode or the vertical kink mode.

We model the possible afterglow of the Fermi GBM event associated with LIGO detection GW150914, under the assumption that the gamma-ray are produced by a short GRB-like relativistic outflow. We model GW150914-GBM as both a weak, on-axis short GRB and normal short GRB seen far off axis. Given the large uncertainty in the position of GW150914, we determine that the best chance of finding the afterglow is with ASKAP or possibly the MWA, with the flux from an off-axis short GRB reaching 0.2 - 4 mJy (0.12 - 16 mJy) at 150 MHz (863.5 MHz) by 1 - 12 months after the initial event. At low frequencies, the source would evolve from a hard to soft spectrum over several months. The radio afterglow would be detectable for several months to years after it peaks, meaning the afterglow may still be detectable and increasing in brightness NOW (mid-July 2016). With a localization from the MWA or ASKAP, the afterglow would be detectable at higher radio frequencies with the ATCA and in X-rays with Chandra or XMM.

The statistical examination of solar flares is crucial to revealing their global characteristics and behaviour. Such examinations can tackle large-scale science questions or give context to detailed single-event studies. However, they are often performed using standard but basic flare detection algorithms relying on arbitrary thresholds. This arbitrariness may lead to important scientific conclusions being drawn from results caused by subjective choices in algorithms rather than the true nature of the Sun. In this paper, we explore the effect of the arbitrary thresholds used in the GOES (Geostationary Operational Environmental Satellite) event list and LYRA (Large Yield RAdiometer) Flare Finder algorithms. We find that there is a small but significant relationship between the power law exponent of the GOES flare peak flux frequency distribution and the flare start thresholds of the algorithms. We also find that the power law exponents of these distributions are not stable, but appear to steepen with increasing peak flux. This implies that the observed flare size distribution may not be a power law at all. We show that depending on the true value of the exponent of the flare size distribution, this deviation from a power law may be due to flares missed by the flare detection algorithms. However, it is not possible determine the true exponent from GOES/XRS observations. Additionally we find that the PROBA2/LYRA flare size distributions are artificially steep and clearly non-power law. We show that this is consistent with an insufficient degradation correction. This means that PROBA2/LYRA should not be used for flare statistics or energetics unless degradation is adequately accounted for. However, it can be used to study variations over shorter timescales and for space weather monitoring.

We present new Gemini Planet Imager observations of the young exoplanet 51 Eridani b which provide further evidence that the companion is physically associated with 51 Eridani. Combining this new astrometric measurement with those reported in the literature, we significantly reduce the posterior probability that 51 Eridani b is an unbound foreground or background T-dwarf in a chance alignment with 51 Eridani to $2\times10^{-7}$, an order of magnitude lower than previously reported. If 51 Eridani b is indeed a bound object, then we have detected orbital motion of the planet between the discovery epoch and the latest epoch. By implementing a computationally efficient Monte Carlo technique, preliminary constraints are placed on the orbital parameters of the system. The current set of astrometric measurements suggest an orbital semimajor axis of $14^{+7}_{-3}$ AU, corresponding to a period of $41^{+35}_{-12}$ years (assuming a mass of $1.75$ $M_{\odot}$ for the central star), and an inclination of $138^{+15}_{-13}$ deg. The remaining orbital elements are only marginally constrained by the current measurements. These preliminary values suggest an orbit which does not share the same inclination as the orbit of the distant M-dwarf binary, GJ 3305, which is a wide physically bound companion to 51 Eridani.

We compare the ability of 11 Differential Emission Measure (DEM) forward-fitting and inversion methods to constrain the properties of active regions and solar flares by simulating synthetic data using the instrumental response functions of SDO/AIA, SDO/EVE, RHESSI, and GOES/XRS. The codes include the single-Gaussian DEM, a bi-Gaussian DEM, a fixed-Gaussian DEM, a linear spline DEM, the spatial synthesis DEM, the Monte-Carlo Markov chain DEM, the regularized DEM inversion, the Hinode/XRT method, a polynomial spline DEM, an EVE+GOES, and an EVE+RHESSI method. Averaging the results from all 11 DEM methods, we find the following accuracies in the inversion of physical parameters: the EM-weighted temperature $T_w^{fit}/T_w^{sim}=0.9\pm0.1$, the peak emission measure $EM_p^{fit}/EM_p^{sim}=0.6\pm0.2$, the total emission measure $EM_t^{fit}/EM_t^{sim}=0.8\pm0.3$, and the multi-thermal energies $E_{th}^{fit}/EM_{th}^{sim}=1.2\pm0.4$. We find that the AIA spatial synthesis, the EVE+GOES, and the EVE+RHESSI method yield the most accurate results.

We present the second part of a project on the global energetics of solar flares and CMEs that includes about 400 M- and X-class flares observed with AIA/SDO during the first 3.5 years of its mission. In this Paper II we compute the differential emission measure (DEM) distribution functions and associated multi-thermal energies, using a spatially-synthesized Gaussian DEM forward-fitting method. The multi-thermal DEM function yields a significantly higher (by an average factor of $\approx 14$), but more comprehensive (multi-)thermal energy than an isothermal energy estimate from the same AIA data. We find a statistical energy ratio of $E_{th}/E_{diss} \approx 2\%-40\%$ between the multi-thermal energy $E_{th}$ and the magnetically dissipated energy $E_{diss}$, which is an order of magnitude higher than the estimates of Emslie et al.~2012. For the analyzed set of M and X-class flares we find the following physical parameter ranges: $L=10^{8.2}-10^{9.7}$ cm for the length scale of the flare areas, $T_p=10^{5.7}-10^{7.4}$ K for the DEM peak temperature, $T_w=10^{6.8}-10^{7.6}$ K for the emission measure-weighted temperature, $n_p=10^{10.3}-10^{11.8}$ cm$^{-3}$ for the average electron density, $EM_p=10^{47.3}-10^{50.3}$ cm$^{-3}$ for the DEM peak emission measure, and $E_{th}=10^{26.8}-10^{32.0}$ erg for the multi-thermal energies. The deduced multi-thermal energies are consistent with the RTV scaling law $E_{th,RTV} = 7.3 \times 10^{-10} \ T_p^3 L_p^2$, which predicts extremal values of $E_{th,max} \approx 1.5 \times 10^{33}$ erg for the largest flare and $E_{th,min} \approx 1 \times 10^{24}$ erg for the smallest coronal nanoflare. The size distributions of the spatial parameters exhibit powerlaw tails that are consistent with the predictions of the fractal-diffusive self-organized criticality model combined with the RTV scaling law.

In this paper, the cooling of 72 M- and X-class flares is examined using GOES/XRS and SDO/EVE. The observed cooling rates are quantified and the observed total cooling times are compared to the predictions of an analytical 0-D hydrodynamic model. It is found that the model does not fit the observations well, but does provide a well defined lower limit on a flare's total cooling time. The discrepancy between observations and the model is then assumed to be primarily due to heating during the decay phase. The decay phase heating necessary to account for the discrepancy is quantified and found be ~50% of the total thermally radiated energy as calculated with GOES. This decay phase heating is found to scale with the observed peak thermal energy. It is predicted that approximating the total thermal energy from the peak is minimally affected by the decay phase heating in small flares. However, in the most energetic flares the decay phase heating inferred from the model can be several times greater than the peak thermal energy.

We test the compatibility and biases of multi-thermal flare DEM (differential emission measure) peak temperatures determined with AIA with those determined by GOES and RHESSI using the isothermal assumption. In a set of 149 M- and X-class flares observed during the first two years of the SDO mission, AIA finds DEM peak temperatures at the time of the peak GOES 1-8 A flux to have an average of Tp = 12.0+/-2.9 MK and Gaussian DEM widths of log10(sigma_T) = 0.50+/-0.13. From GOES observations of the same 149 events, a mean temperature of Tp = 15.6+/-2.4 MK is inferred, which is systematically higher by a factor of TGOES/TAIA = 1.4+/-0.4. We demonstrate that this discrepancy results from the isothermal assumption in the inversion of the GOES filter ratio. From isothermal fits to photon spectra at energies of E ~ 6-12 keV of 61 of these events, RHESSI finds the temperature to be higher still by a factor of TRHESSI/TAIA = 1.9+/-1.0. We find that this is partly a consequence of the isothermal assumption. However, RHESSI is not sensitive to the low-temperature range of the DEM peak, and thus RHESSI samples only the high-temperature tail of the DEM function. This can also contribute to the discrepancy between AIA and RHESSI temperatures. The higher flare temperatures found by GOES and RHESSI imply correspondingly lower emission measures. We conclude that self-consistent flare DEM temperatures and emission measures require simultaneous fitting of EUV (AIA) and soft X-ray (GOES and RHESSI) fluxes.

Solar flare X-ray emission results from rapidly increasing temperatures and emission measures in flaring active region loops. To date, observations from the X-Ray Sensor (XRS) onboard the Geostationary Operational Environmental Satellite (GOES) have been used to derive these properties, but have been limited by a number of factors, including the lack of a consistent background subtraction method capable of being automatically applied to large numbers of flares. In this paper, we describe an automated temperature and emission measure-based background subtraction method (TEBBS), which builds on the methods of Bornmann (1990). Our algorithm ensures that the derived temperature is always greater than the instrumental limit and the pre-flare background temperature, and that the temperature and emission measure are increasing during the flare rise phase. Additionally, TEBBS utilizes the improved estimates of GOES temperatures and emission measures from White et al. (2005). TEBBS was successfully applied to over 50,000 solar flares occurring over nearly three solar cycles (1980-2007), and used to create an extensive catalog of the solar flare thermal properties. We confirm that the peak emission measure and total radiative losses scale with background subtracted GOES X-ray flux as power-laws, while the peak temperature scales logarithmically. As expected, the peak emission measure shows an increasing trend with peak temperature, although the total radiative losses do not. While these results are comparable to previous studies, we find that flares of a given GOES class have lower peak temperatures and higher peak emission measures than previously reported. The resulting TEBBS database of thermal flare plasma properties is publicly available on Solar Monitor (www.solarmonitor.org/TEBBS/) and will be available on Heliophysics Integrated Observatory (www.helio-vo.eu).