We describe a new algorithm for reconstruction of Differential Emission Measures (DEMs) in the solar corona. Although a number of such algorithms currently exist, they can have difficulty converging for some cases, and can be complex, slow, or idiosyncratic in their output (i.e., their inversions can have features that are a result of the inversion code and instrument response, not of the solar source); we will document some of these issues in this paper. The new algorithm described here significantly reduces these drawbacks and is particularly notable for its simplicity; it is reproduced here, in full, on a single page. After we describe the algorithm, we compare its performance and fidelity with some prevalent methods. Although presented here for extreme ultraviolet (EUV) data, the algorithm is robust and extensible to any other wavelengths (e.g., X-rays) where the DEM treatment is valid.

This is the last of three papers describing an `absolute' calibration of the GONG magnetograph using and end-to-end simulation of its measurement process. The simulation begins with a MURaM 3D MHD datacube and ends with a `synthetic magnetogram' of the corresponding magnetic field values as they would be observed by GONG. We determine a calibration by comparing the synthetic magnetic field measurements with the MURaM magnetic field values that produced them. The previous two papers have described the GONG measurement process (both instrument and data processing), our simulation of it, and the theory of magnetogram comparison and calibration. In this paper, we address some final points on calibration, combine all of this work into a set of calibration curves, and consider the results. We also review the results of the previous two papers for locality of reference. Our calibration indicates that GONG magnetograms underestimate weak flux by a factor of over 2 near disk center, but that factor decreases to about 1 as the line-of-sight approaches the limb. A preliminary investigation of the generalizability of these results suggests other instruments will be affected in a similar way. We also find that some differences in previous magnetograph comparisons are artifacts of instrumental resolution which do not reflect an intrinsic calibration difference, and the measurements are more similar than sometimes thought. These results are directly applicable to question of solar wind prediction model accuracies, particularly in the search for the cause of the common discrepancy between predicted solar wind magnetic flux at 1 AU and values measured in situ by current satellite missions.

This is the second of three papers describing an `absolute' calibration of the GONG magnetograph using an end-to-end simulation of its measurement process. In the first paper, we described the GONG instrument and our `end-to-end' simulation of its measurement process. In this paper, we consider the theory of calibration, and magnetograph comparison in general, identifying some of the significant issues and pitfalls. The calibration of a magnetograph is a function of whether or not it preserves flux, independent of its spatial resolution. However, we find that the one-dimensional comparison methods most often used for magnetograph calibration and comparison will show dramatic differences between two magnetograms with differing spatial resolution, even if they both preserve flux. Some of the apparent disagreement between magnetograms found in the literature are likely a result of these instrumental resolution differences rather than any intrinsic calibration differences. To avoid them, spatial resolution must be carefully matched prior to comparing magnetograms or making calibration curves. In the third paper, we apply the lessons learned here to absolute calibration of GONG using our `end-to-end' measurement simulation.

This is the first of three papers describing an `absolute' calibration of the GONG magnetograph using an end-to-end simulation of its measurement process. The input to this simulation is a MURaM 3D MHD photospheric simulation and the output is the corresponding set of simulated data numbers which would be recorded by the GONG detectors. These simulated data numbers are then used to produce `synthetic magnetograms' which can be compared with the simulation inputs. This paper describes the GONG instrument, the MURaM datacube, our instrument simulator, and calculation of synthetic magnetograms, setting the stage for the subsequent two papers. These will first lay groundwork for calibration (and magnetogram comparison in general), then apply them to calibration of GONG using the simulation results.

Events in the solar corona are often widely separated in their timescales, which can allow them to be identified when they would otherwise be confused with emission from other sources in the corona. Methods for cleanly separating such events based on their timescales are thus desirable for research in the field. This paper develops a technique for identifying time-varying signals in solar coronal image sequences which is based on a per-pixel running median filter and an understanding of photon-counting statistics. Example applications to 'EIT Waves' and small-scale dynamics are shown, both using data from the 193 Angstrom channel on AIA. The technique is found to discriminate EIT Waves more cleanly than the running and base difference techniques most commonly used. It is also demonstrated that there is more signal in the data than is commonly appreciated, finding that the waves can be traced to the edge of the AIA field of view when the data are rebinned to increase the signal-to-noise ratio.

We present a fast method for reconstructing Differential Emission Measures (DEMs) using solar coronal data. On average, the method computes over 1000 DEMs per second for a sample active region observed by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO), and achieves reduced chi-squared of order unity with no negative emission in all but a few test cases. The high performance of this method is especially relevant in the context of AIA, which images of order one million solar pixels per second. This paper describes the method, analyzes its fidelity, compares its performance and results with other DEM methods, and applies it to an active region and loop observed by AIA and by the Extreme-ultraviolet Imaging Spectrometer (EIS) on Hinode.

A number of scenarios have been proposed for the origin of the supermassive black holes (SMBHs) that are found in the centres of most galaxies. Many such scenarios predict a high-redshift population of massive black holes (MBHs), with masses in the range 100 to 100000 times that of the Sun. When the Laser Interferometer Space Antenna (LISA) is finally operational, it is likely that it will detect on the order of 100 of these MBH binaries as they merge. The differences between proposed population models produce appreciable effects in the portion of the population which is detectable by LISA, so it is likely that the LISA observations will allow us to place constraints on them. However, gravitational wave detectors such as LISA will not be able to detect all such mergers nor assign precise black hole parameters to the merger, due to weak gravitational wave signal strengths. This paper explores LISA's ability to distinguish between several MBH population models. In this way, we go beyond predicting a LISA observed population and consider the extent to which LISA observations could inform astrophysical modellers. The errors in LISA parameter estimation are applied with a direct method which generates random sample parameters for each source in a population realisation. We consider how the distinguishability varies depending on the choice of source parameters (1 or 2 parameters chosen from masses, redshift or spins) used to characterise the model distributions, with confidence levels determined by 1 and 2-dimensional tests based on the Kolmogorov-Smirnov test.

Many scenarios have been proposed for the origin of the supermassive black holes (SMBHs) that are found in the centres of most galaxies. Many of these formation scenarios predict a high-redshift population of intermediate-mass black holes (IMBHs), with masses in the range 100 to 100000 times that of the Sun. A powerful way to observe these IMBHs is via gravitational waves the black holes emit as they merge. The statistics of the observed black hole population should, in principle, allow us to discriminate between competing astrophysical scenarios for the origin and formation of SMBHs. However, gravitational wave detectors such as LISA will not be able to detect all such mergers nor assign precise black hole parameters to the merger, due to weak gravitational wave signal strengths. In order to use LISA observations to infer the statistics of the underlying population, these errors must be taken into account. We describe here a method for folding the LISA gravitational wave parameter error estimates into an `error kernel' designed for use at the population model level. The effects of this error function are demonstrated by applying it to several recent models of black hole mergers, and some tentative conclusions are made about LISA's ability to test scenarios of the origin and formation of supermassive black holes.

One puzzling observed property of coronal loops is that they are of roughly constant thickness along their length. Various studies have found no consistent pattern of width variation along the length of loops observed by TRACE and SOHO. This is at odds with expectations of magnetic flux tube expansion properties, which suggests that loops are widest at their tops, and significantly narrower at their footpoints. Coronal loops correspond to areas of the solar corona which have been preferentially heated by some process, so this observed property might be connected to the mechanisms that heat the corona. One means of energy deposition is magnetic reconnection, which occurs along field lines called separators. These field lines begin and end on magnetic null points, and loops forming near them can therefore be relatively wide at their bases. Thus, coronal energization by magnetic reconnection may replicate the puzzling expansion properties observed in coronal loops. We present results of a Monte Carlo survey of separator field line expansion properties, comparing them to the observed properties of coronal loops.