We present the first reconstruction of dark matter maps from weak lensing observational data using deep learning. We train a convolution neural network (CNN) with a Unet based architecture on over $3.6\times10^5$ simulated data realizations with non-Gaussian shape noise and with cosmological parameters varying over a broad prior distribution. We interpret our newly created DES SV map as an approximation of the posterior mean $P(\kappa | \gamma)$ of the convergence given observed shear. Our DeepMass method is substantially more accurate than existing mass-mapping methods. With a validation set of 8000 simulated DES SV data realizations, compared to Wiener filtering with a fixed power spectrum, the DeepMass method improved the mean-square-error (MSE) by 11 per cent. With N-body simulated MICE mock data, we show that Wiener filtering with the optimal known power spectrum still gives a worse MSE than our generalized method with no input cosmological parameters; we show that the improvement is driven by the non-linear structures in the convergence. With higher galaxy density in future weak lensing data unveiling more non-linear scales, it is likely that deep learning will be a leading approach for mass mapping with Euclid and LSST.

When inferring unknown parameters or comparing different models, data must be compared to underlying theory. Even if a model has no closed-form solution to derive summary statistics, it is often still possible to simulate mock data in order to generate theoretical predictions. For realistic simulations of noisy data, this is identical to drawing realizations of the data from a likelihood distribution. Though the estimated summary statistic from simulated data vectors may be unbiased, the estimator has variance which should be accounted for. We show how to correct the likelihood in the presence of an estimated summary statistic by marginalizing over the true summary statistic in the framework of a Bayesian hierarchical model. For Gaussian likelihoods where the covariance must also be estimated from simulations, we present an alteration to the Sellentin-Heavens corrected likelihood. We show that excluding the proposed correction leads to an incorrect estimate of the Bayesian evidence with JLA data. The correction is highly relevant for cosmological inference that relies on simulated data for theory (e.g. weak lensing peak statistics and simulated power spectra) and can reduce the number of simulations required.

We present a simulated cosmology analysis using the second and third moments of the weak lensing mass (convergence) maps. The analysis is geared towards the third year (Y3) data from the Dark Energy Survey (DES), but the methodology can be applied to other weak lensing data sets. The second moment, or variances, of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. We present the formalism for obtaining the convergence maps from the measured shear and for obtaining the second and third moments of these maps given partial sky coverage. We estimate the covariance matrix from a large suite of numerical simulations. We test our pipeline through a simulated likelihood analyses varying 5 cosmological parameters and 10 nuisance parameters. Our forecast shows that the combination of second and third moments provides a 1.5 percent constraint on $S_8 \equiv \sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ for DES Y3 data. This is 20 percent better than an analysis using a simulated DES Y3 shear 2-point statistics, owing to the non-Gaussian information captured by the inclusion of higher-order statistics. The methodology developed here can be applied to current and future weak lensing datasets to make use of information from non-Gaussianity in the cosmic density field and to improve constraints on cosmological parameters.

The observed properties (i.e., source size, source position, time duration, decay time) of solar radio emission produced through plasma processes near the local plasma frequency, and hence the interpretation of solar radio bursts, are strongly influenced by propagation effects in the inhomogeneous turbulent solar corona. In this work, a 3D stochastic description of the propagation process is presented, based on the Fokker-Planck and Langevin equations of radio-wave transport in a medium containing anisotropic electron density fluctuations. Using a numerical treatment based on this model, we investigate the characteristic source sizes and burst decay times for Type III solar radio bursts. Comparison of the simulations with the observations of solar radio bursts shows that predominantly perpendicular density fluctuations in the solar corona are required, with an anisotropy factor $\sim 0.3$ for sources observed at around 30~MHz. The simulations also demonstrate that the photons are isotropized near the region of primary emission, but the waves are then focused by large-scale refraction, leading to plasma radio emission directivity that is characterized by a half-width-half-maximum of about 40~degrees near 30~MHz. The results are applicable to various solar radio bursts produced via plasma emission.

One of the key problems in solar flare physics is the determination of the low-energy cut-off; the value that determines the energy of nonthermal electrons and hence flare energetics. We discuss different approaches to determine the low-energy cut-off in the spectrum of accelerated electrons: (i) the total electron number model, (ii) the time-of-flight model (based on the equivalence of the time-of-flight and the collisional deflection time); (iii) the warm target model of Kontar et al.~(2015), and (iv) the model of the spectral cross-over between thermal and nonthermal components. We find that the first three models are consistent with a low-energy cutoff with a mean value of $\approx 10$ keV, while the cross-over model provides an upper limit for the low-energy cutoff with a mean value of $ \approx 21$ keV. Combining the first three models we find that the ratio of the nonthermal energy to the dissipated magnetic energy in solar flares has a mean value of $q_E=0.57\pm0.08$, which is consistent with an earlier study based on the simplified approximation of the warm target model alone ($q_E=0.51\pm0.17$). This study corroborates the self-consistency between three different low-energy cutoff models in the calculation of nonthermal flare energies.

During solar flares, a large fraction of the released magnetic energy is carried by energetic electrons that transfer and deposit energy in the Sun's atmosphere. Electron transport is often approximated by a cold thick-target model (CTTM), assuming that electron energy is much larger than the temperature of the ambient plasma, and electron energy evolution is modeled as a systematic loss. Using kinetic modeling of electrons, we re-evaluate the transport and deposition of flare energy. Using a full collisional warm-target model (WTM), we account for electron thermalization and for the properties of the ambient coronal plasma such as its number density, temperature and spatial extent. We show that the deposition of non-thermal electron energy in the lower atmosphere is highly dependent on the properties of the flaring coronal plasma. In general, thermalization and a reduced WTM energy loss rate leads to an increase of non-thermal energy transferred to the chromosphere, and the deposition of non-thermal energy at greater depths. The simulations show that energy is deposited in the lower atmosphere initially by high energy non-thermal electrons, and later by lower energy non-thermal electrons that partially or fully thermalize in the corona, over timescales of seconds, unaccounted for in previous studies. This delayed heating may act as a diagnostic of both the injected non-thermal electron distribution and the coronal plasma, vital for constraining flare energetics.

We present new wide-field weak lensing mass maps for the Year 1 Dark Energy Survey data, generated via a forward fitting approach. This method of producing maps does not impose any prior constraints on the mass distribution to be reconstructed. The technique is found to improve the map reconstruction on the edges of the field compared to the conventional Kaiser-Squires method, which applies a direct inversion on the data; our approach is in good agreement with the previous direct approach in the central regions of the footprint. The mapping technique is assessed and verified with tests on simulations; together with the Kaiser-Squires method, the technique is then applied to data from the Dark Energy Survey Year 1 data and the differences between the two methods are compared. We also produce the first DES measurements of the convergence Minkowski functionals and compare them to those measured in simulations.

Solar flare hard X-ray spectroscopy serves as a key diagnostic of the accelerated electron spectrum. However, the standard approach using the collisional cold thick-target model poorly constrains the lower-energy part of the accelerated electron spectrum, and hence the overall energetics of the accelerated electrons are typically constrained only to within one or two orders of magnitude. Here we develop and apply a physically self-consistent warm-target approach which involves the use of both hard X-ray spectroscopy and imaging data. The approach allows an accurate determination of the electron distribution low-energy cutoff, and hence the electron acceleration rate and the contribution of accelerated electrons to the total energy released, by constraining the coronal plasma parameters. Using a solar flare observed in X-rays by the {\em RHESSI} spacecraft, we demonstrate that using the standard cold-target methodology, the low-energy cutoff (and hence the energy content in electrons) is essentially undetermined. However, the warm-target methodology can determine the low-energy electron cutoff with $\sim$7\% uncertainty at the $3\sigma$ level and hence permits an accurate quantitative study of the importance of accelerated electrons in solar flare energetics.

We present the first observational study of the onset and evolution of solar flare turbulence in the lower solar atmosphere on an unprecedented time scale of 1.7 s using the Interface Region Imaging Spectrograph observing plasma at a temperature of 80,000 K. At this time resolution, nonthermal spectral line broadening, indicating turbulent velocity fluctuations, precedes the flare onset at this temperature and is coincident with net blue-shifts. The broadening decreases as the flare brightens and then oscillates with a period of ~10 s. These observations are consistent with turbulence in the lower solar atmosphere at the flare onset, heating that region as it dissipates. This challenges the current view of energy release and transport in the standard solar flare model, suggesting that turbulence partly heats the lower atmosphere.

Observations of solar flares at sub-THz frequencies (mm and sub-mm wavelengths) over the last two decades often show a spectral component rising with frequency. Unlike a typical gyrosynchrotron spectrum decreasing with frequency, or a weak thermal component from hot coronal plasma, the observations can demonstrate a high flux level (up to ~10^4 s.f.u. at 0.4 THz) and fast variability on sub-second time scales. Although, many models has been put forward to explain the puzzling observations, none of them have clear observational support. Here we propose a scenario to explain the intriguing sub-THz observations. We show that the model, based on free-free emission from the plasma of flare ribbons at temperatures 10^4-10^6 K, is consistent with all existing observations of frequency-rising sub-THz flare emission. The model provides a temperature diagnostic of the flaring chromosphere and suggests fast heating and cooling of the dense transition region plasma.

We use Dataflow Engines (DFE) to construct an efficient Wiener filter of noisy and incomplete image data, and to quickly draw probabilistic samples of the compatible true underlying images from the Wiener posterior. Dataflow computing is a powerful approach using reconfigurable hardware, which can be deeply pipelined and is intrinsically parallel. The unique Wiener-filtered image is the minimum-variance linear estimate of the true image (if the signal and noise covariances are known) and the most probable true image (if the signal and noise are Gaussian distributed). However, many images are compatible with the data with different probabilities, given by the analytic posterior probability distribution referred to as the Wiener posterior. The DFE code also draws large numbers of samples of true images from this posterior, which allows for further statistical analysis. Naive computation of the Wiener-filtered image is impractical for large datasets, as it scales as $n^3$, where $n$ is the number of pixels. We use a messenger field algorithm, which is well suited to a DFE implementation, to draw samples from the Wiener posterior, that is, with the correct probability we draw samples of noiseless images that are compatible with the observed noisy image. The Wiener-filtered image can be obtained by a trivial modification of the algorithm. We demonstrate a lower bound on the speed-up, from drawing 10$^5$ samples of a 128$^2$ image, of 11.3 ${\pm}$ 0.8 with 8 DFEs in a 1U MPC-X box when compared with a 1U server presenting 32 CPU threads. We also discuss a potential application in astronomy, to provide better dark matter maps and improved determination of the parameters of the Universe.

Mapping the underlying density field, including non-visible dark matter, using weak gravitational lensing measurements is now a standard tool in cosmology. Due to its importance to the science results of current and upcoming surveys, the quality of the convergence reconstruction methods should be well understood. We compare three methods: Kaiser-Squires (KS), Wiener filter, and GLIMPSE. KS is a direct inversion, not accounting for survey masks or noise. The Wiener filter is well-motivated for Gaussian density fields in a Bayesian framework. GLIMPSE uses sparsity, aiming to reconstruct non-linearities in the density field. We compare these methods with several tests using public Dark Energy Survey (DES) Science Verification (SV) data and realistic DES simulations. The Wiener filter and GLIMPSE offer substantial improvements over smoothed KS with a range of metrics. Both the Wiener filter and GLIMPSE convergence reconstructions show a 12 per cent improvement in Pearson correlation with the underlying truth from simulations. To compare the mapping methods' abilities to find mass peaks, we measure the difference between peak counts from simulated {\Lambda}CDM shear catalogues and catalogues with no mass fluctuations (a standard data vector when inferring cosmology from peak statistics); the maximum signal-to-noise of these peak statistics is increased by a factor of 3.5 for the Wiener filter and 9 for GLIMPSE. With simulations we measure the reconstruction of the harmonic phases; the phase residuals' concentration is improved 17 per cent by GLIMPSE and 18 per cent by the Wiener filter. The correlation between reconstructions from data and foreground redMaPPer clusters is increased 18 per cent by the Wiener filter and 32 per cent by GLIMPSE.

Grain surface chemistry and its treatment in gas-grain chemical models is an area of large uncertainty. Whilst laboratory experiments are making progress, there is still much that is unknown about grain surface chemistry. Further, the results and parameters produced by experiment are often not easily translated to the rate equation approach most commonly used in astrochemical modelling. It is possible that statistical methods can reduce the uncertainty in grain surface chemical networks. In this work, a simple model of grain surface chemistry in a molecular cloud is developed and a Bayesian inference of the reactions rates is performed through MCMC sampling. Using observational data of the solid state abundances of major chemical species in molecular clouds, the posterior distributions for the rates of seven reactions producing CO, CO$_2$, CH$_3$OH and H$_2$O are calculated, in a form that is suitable for rate equation models. This represents a vital first step in the development of a method to infer reaction rates from observations of chemical abundances in astrophysical environments.

Solar radio observations provide a unique diagnostic of the outer solar atmosphere. However, the inhomogeneous turbulent corona strongly affects the propagation of the emitted radio waves, so decoupling the intrinsic properties of the emitting source from the effects of radio-wave propagation has long been a major challenge in solar physics. Here we report quantitative spatial and frequency characterization of solar radio burst fine structures observed with the LOw Frequency Array (LOFAR), an instrument with high time resolution that also permits imaging at scales much shorter than those corresponding to radio-wave propagation in the corona. The observations demonstrate that radio-wave propagation effects, and not the properties of the intrinsic emission source, dominate the observed spatial characteristics of radio burst images. These results permit more accurate estimates of source brightness temperatures, and open opportunities for quantitative study of the mechanisms that create the turbulent coronal medium through which the emitted radiation propagates.

In situ measurements of the fast solar wind reveal non-thermal distributions of electrons, protons and, minor ions extending from $0.3$ AU to the heliopause. The physical mechanisms responsible for these non-thermal properties and the location where these properties originate remain open questions. Here we present spectroscopic evidence, from extreme ultraviolet spectroscopy, that the velocity distribution functions (VDFs) of minor ions are already non-Gaussian at the base of the fast solar wind in a coronal hole, at altitudes of $< 1.1 R_{\odot}$. Analysis of Fe, Si, and Mg spectral lines reveal a peaked line-shape core and broad wings that can be characteristed by a kappa VDF. A kappa distribution fit gives very small kappa indices off-limb of $\kappa\approx1.9-2.5$, indicating either (a) ion populations far from thermal equilibrium, (b) fluid motions such as non-Gaussian turbulent fluctuations or non-uniform wave motions, or (c) some combination of both. These observations provide important empirical constraints for the source region of the fast solar wind and for the theoretical models of the different acceleration, heating, and energy deposition processes therein. To the best of our knowledge, this is the first time that the ion VDF in the fast solar wind has been probed so close to its source region. The findings are also a timely precursor to the upcoming 2018 launch of the Parker Solar Probe, which will provide the closest in situ measurements of the solar wind at approximately $0.04$ AU ($8.5$ solar radii).

We construct the largest curved-sky galaxy weak lensing mass map to date from the DES first-year (DES Y1) data. The map, about 10 times larger than previous work, is constructed over a contiguous $\approx1,500 $deg$^2$, covering a comoving volume of $\approx10 $Gpc$^3$. The effects of masking, sampling, and noise are tested using simulations. We generate weak lensing maps from two DES Y1 shear catalogs, Metacalibration and Im3shape, with sources at redshift $0.2<z<1.3,$ and in each of four bins in this range. In the highest signal-to-noise map, the ratio between the mean signal-to-noise in the E-mode and the B-mode map is $\sim$1.5 ($\sim$2) when smoothed with a Gaussian filter of $\sigma_{G}=30$ (80) arcminutes. The second and third moments of the convergence $\kappa$ in the maps are in agreement with simulations. We also find no significant correlation of $\kappa$ with maps of potential systematic contaminants. Finally, we demonstrate two applications of the mass maps: (1) cross-correlation with different foreground tracers of mass and (2) exploration of the largest peaks and voids in the maps.

Magnetic reconnection is a well-accepted part of the theory of solar eruptive events, though the evidence is still circumstantial. Intrinsic to the reconnection picture of a solar eruptive event, particularly in the standard model for two-ribbon flares ("CSHKP" model), are an advective flow of magnetized plasma into the reconnection region, expansion of field above the reconnection region as a flux rope erupts, retraction of heated post-reconnection loops, and downflows of cooling plasma along those loops. We report on a unique set of SDO/AIA imaging and Hinode/EIS spectroscopic observations of the disk flare SOL2016-03-23T03:54 in which all four flows are present simultaneously. This includes spectroscopic evidence for a plasma upflow in association with large-scale expanding closed inflow field. The reconnection inflows are symmetric, and consistent with fast reconnection, and the post-reconnection loops show a clear cooling and deceleration as they retract. Observations of coronal reconnection flows are still rare, and most events are observed at the solar limb, obscured by complex foregrounds, making their relationship to the flare ribbons, cusp field and arcades formed in the lower atmosphere difficult to interpret. The disk location and favorable perspective of this event have removed these ambiguities giving a clear picture of the reconnection dynamics.

During the next solar maximum, two upcoming space-borne X-ray missions, STIX on board Solar Orbiter and MiSolFA, will perform stereoscopic X-ray observations of solar flares at two different locations: STIX at 0.28 AU (at perihelion) and up to inclinations of $\sim25^{\circ}$, and MiSolFA in a low-Earth orbit. The combined observations from these cross-calibrated detectors will allow us to infer the electron anisotropy of individual flares confidently for the first time. We simulated both instrumental and physical effects for STIX and MiSolFA including thermal shielding, background and X-ray Compton backscattering (albedo effect) in the solar photosphere. We predict the expected number of observable flares available for stereoscopic measurements during the next solar maximum. We also discuss the range of useful spacecraft observation angles for the challenging case of close-to-isotropic flare anisotropy. The simulated results show that STIX and MiSolFA will be capable of detecting low levels of flare anisotropy, for M1-class or stronger flares, even with a relatively small spacecraft angular separation of 20-30{\deg}. Both instruments will directly measure the flare X-ray anisotropy of about 40 M- and X-class solar flares during the next solar maximum. Near-future stereoscopic observations with Solar Orbiter/STIX and MiSolFA will help distinguishing between competing flare-acceleration mechanisms, and provide essential constraints regarding collisional and non-collisional transport processes occurring in the flaring atmosphere for individual solar flares.

The energy released in solar flares derives from a reconfiguration of magnetic fields to a lower energy state, and is manifested in several forms, including bulk kinetic energy of the coronal mass ejection, acceleration of electrons and ions, and enhanced thermal energy that is ultimately radiated away across the electromagnetic spectrum from optical to X-rays. Using an unprecedented set of coordinated observations, from a suite of instruments, we here report on a hitherto largely overlooked energy component -- the kinetic energy associated with small-scale turbulent mass motions. We show that the spatial location of, and timing of the peak in, turbulent kinetic energy together provide persuasive evidence that turbulent energy may play a key role in the transfer of energy in solar flares. Although the kinetic energy of turbulent motions accounts, at any given time, for only $\sim (0.5-1)$\% of the energy released, its relatively rapid ($\sim$$1-10$~s) energization and dissipation causes the associated throughput of energy (i.e., power) to rival that of major components of the released energy in solar flares, and thus presumably in other astrophysical acceleration sites.

In a solar flare, a large fraction of the magnetic energy released is converted rapidly to the kinetic energy of non-thermal particles and bulk plasma motion. This will likely result in non-equilibrium particle distributions and turbulent plasma conditions. We investigate this by analysing the profiles of high-temperature extreme ultraviolet emission lines from a major flare (SOL2014-03-29T17:44) observed by the EUV Imaging Spectrometer (EIS) on Hinode. We find that in many locations the line profiles are non-Gaussian, consistent with a kappa-distribution of emitting ions with properties that vary in space and time. At the flare footpoints, close to sites of hard X-ray emission from non-thermal electrons, the kappa-index for the Fe XVI 262.976 angstrom line at 3 MK takes values of 3-5. In the corona, close to a low-energy HXR source, the Fe XXIII 263.760 angstrom line at 15 MK shows kappa values of typically 4-7. The observed trends in the kappa parameter show that we are most likely detecting the properties of the ion population rather than any instrumental effects. We calculate that a non-thermal ion population could exist if locally accelerated on timescales < 0.1 s. However, observations of net redshifts in the lines also imply the presence of plasma downflows which could lead to bulk turbulence, with increased non-Gaussianity in cooler regions. Both interpretations have important implications for theories of solar flare particle acceleration.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

The properties of solar flare plasma can be determined from the observation of optically thin lines. The emitting ion distribution determines the shape of the line, with an isothermal Maxwellian ion distribution producing a Gaussian profile. Non-Gaussian lines may indicate more complex ion distributions. We investigate the possibility of determining flare-accelerated non-thermal ion or plasma velocity distributions. We study EUV spectral lines during a flare SOL2013-05-15T01:45 using the Hinode EUV Imaging Spectrometer (EIS). The flare is located close to the eastern solar limb with an extended loop structure, allowing the different flare features: ribbons, hard X-ray (HXR) footpoints and the loop-top source to be clearly observed in UV, EUV and X-rays. EUV line spectroscopy is performed in seven different regions. We study the line profiles of isolated and unblended Fe XVI lines (262.976 A) formed at temperatures of 2-4 MK. Fe XVI line profiles at one time close to the peak soft X-ray emission and free of directed mass motions are examined using: 1. a higher moments analysis, 2. Gaussian fitting, and 3. by fitting a kappa distribution profile convolved with a Gaussian to account for the EIS instrumental profile. Fe XVI line profiles in the loop-top, HXR footpoint and ribbon regions can be confidently fitted with a kappa line profile, giving low, non-thermal kappa values between 2-3.3. A higher moments analysis also finds that many of the line kurtosis values are higher than the Gaussian value of 3, even with the presence of a broad Gaussian instrumental profile. A flare-accelerated non-thermal ion population could account for both the observed non-Gaussian line profiles, and for the Fe XVI excess broadening found from Gaussian fitting, if the emitting ions are interacting with a thermalised 4 MK electron population, and the instrumental profile is well-approximated by a Gaussian.

The detailed knowledge of plasma heating and acceleration region properties presents a major observational challenge in solar flare physics. Using the Ramaty High Energy Solar Spectroscopic Imager (RHESSI), the high temperature differential emission measure, DEM(T), and the energy-dependent spatial structure of solar flare coronal sources are studied quantitatively. The altitude of the coronal X-ray source is observed to increase with energy by ~+0.2 arcsec/keV between 10 and 25 keV. Although an isothermal model can fit the thermal X-ray spectrum observed by RHESSI, such a model cannot account for the changes in altitude, and multi-thermal coronal sources are required where the temperature increases with altitude. For the first time, we show how RHESSI imaging information can be used to constrain the DEM(T) of a flaring plasma. We develop a thermal bremsstrahlung X-ray emission model with inhomogeneous temperature and density distributions to simultaneously reproduce: i) DEM(T), ii) altitude as a function of energy, and iii) vertical extent of the flaring coronal source versus energy. We find that the temperature-altitude gradient in the region is ~+0.08 keV/arcsec (~1.3 MK/Mm). Similar altitude-energy trends in other flares suggest that the majority of coronal X-ray sources are multi-thermal and have strong vertical temperature and density gradients with a broad DEM(T).

X-ray observations are a direct diagnostic of fast electrons produced in solar flares, energized during the energy release process and directed towards the Sun. Since the properties of accelerated electrons can be substantially changed during their transport and interaction with the background plasma, a model must ultimately be applied to X-ray observations in order to understand the mechanism responsible for their acceleration. A cold thick target model is ubiquitously used for this task, since it provides a simple analytic relationship between the accelerated electron spectrum and the emitting electron spectrum in the X-ray source, with the latter quantity readily obtained from X-ray observations. However, such a model is inappropriate for the majority of solar flares in which the electrons propagate in a hot megaKelvin plasma, because it does not take into account the physics of thermalization of fast electrons. The use of a more realistic model, properly accounting for the properties of the background plasma, and the collisional diffusion and thermalization of electrons, can alleviate or even remove many of the traditional problems associated with the cold thick target model and the deduction of the accelerated electron spectrum from X-ray spectroscopy, such as the number problem and the need to impose an ad hoc low energy cut-off.

Extending previous studies of nonthermal electron transport in solar flares which include the effects of collisional energy diffusion and thermalization of fast electrons, we present an analytic method to infer more accurate estimates of the accelerated electron spectrum in solar flares from observations of the hard X-ray spectrum. Unlike for the standard cold-target model, the spatial characteristics of the flaring region, especially the necessity to consider a finite volume of hot plasma in the source, need to be taken into account in order to correctly obtain the injected electron spectrum from the source-integrated electron flux spectrum (a quantity straightforwardly obtained from hard X-ray observations). We show that the effect of electron thermalization can be significant enough to nullify the need to introduce an {\it ad hoc} low-energy cutoff to the injected electron spectrum in order to keep the injected power in non-thermal electrons at a reasonable value. Rather the suppression of the inferred low-energy end of the injected spectrum compared to that deduced from a cold-target analysis allows the inference from hard X-ray observations of a more realistic energy in injected non-thermal electrons in solar flares.

X-rays are a valuable diagnostic tool for the study of high energy accelerated electrons. Bremsstrahlung X-rays produced by, and directly related to, high energy electrons accelerated during a flare, provide a powerful diagnostic tool for determining both the properties of the accelerated electron distribution, and of the flaring coronal and chromospheric plasmas. This thesis is specifically concerned with the study of spatial, spectral and polarization properties of solar flare X-ray sources via both modelling and X-ray observations using the Ramaty High Energy Solar Spectroscopic Imager (RHESSI). Firstly, a new model is presented, accounting for finite temperature, pitch angle scattering and initial pitch angle injection. This is developed to accurately infer the properties of the acceleration region from the observations of dense coronal X-ray sources. Moreover, examining how the spatial properties of dense coronal X-ray sources change in time, interesting trends in length, width, position, number density and thermal pressure are found and the possible causes for such changes are discussed. Further analysis of data in combination with the modelling of X-ray transport in the photosphere, allows changes in X-ray source positions and sizes due to the X-ray albedo effect to be deduced. Finally, it is shown, for the first time, how the presence of a photospheric X-ray albedo component produces a spatially resolvable polarization pattern across a hard X-ray (HXR) source. It is demonstrated how changes in the degree and direction of polarization across a single HXR source can be used to determine the anisotropy of the radiating electron distribution.

Observations with {\em RHESSI} have enabled the detailed study of the structure of dense hard X-ray coronal sources in solar flares. The variation of source extent with electron energy has been discussed in the context of streaming of non-thermal particles in a one-dimensional cold-target model, and the results used to constrain both the physical extent of, and density within, the electron acceleration region. Here we extend this investigation to a more physically realistic model of electron transport that takes into account the finite temperature of the ambient plasma, the initial pitch-angle distribution of the accelerated electrons, and the effects of collisional pitch-angle scattering. The finite temperature results in the thermal diffusion of electrons, that leads to the observationally-inferred value of the acceleration region volume being an overestimate of its true value. The different directions of the electron trajectories, a consequence of both the non-zero injection pitch-angle and scattering within the target, cause the projected propagation distance parallel to the guiding magnetic field to be reduced, so that a one-dimensional interpretation can overestimate the actual density by a factor of up to $\sim 6$. The implications of these results for the determination of acceleration region properties (specific acceleration rate, filling factor, etc.) are discussed.

The spatial and spectral properties of three solar flare coronal X-ray loops are studied before, during and after the peak X-ray emission. Using observations from the Ramaty High Energy Solar Spectroscopic Imager (RHESSI), we deduce the temporal changes in emitting X-ray length, corpulence, volume, position, number density and thermal pressure. We observe a decrease in the loop length, width and volume before the X-ray peak, and an increasing number density and thermal pressure. After the X-ray peak, volume increases and loop corpulence grows due to an increasing width. The volume variations are more pronounced than the position variations, often known as magnetic line contraction. We believe this is the first dedicated study of the temporal evolution of X-ray loop lengths and widths. Collectively, the observations also show for the first time three temporal phases given by peaks in temperature, X-ray emission and thermal pressure, with minimum volume coinciding with the X-ray peak. Although the volume of the flaring plasma decreases before the peak in X-ray emission, the relationship between temperature and volume does not support simple compressive heating in a collapsing magnetic trap model. Within a low beta plasma, shrinking loop widths perpendicular to the guiding field can be explained by squeezing the magnetic field threading the region. Plasma heating leads to chromospheric evaporation and growing number density, producing increasing thermal pressure and decreasing loop lengths as electrons interact at shorter distances and we believe after the X-ray peak, the increasing loop corpulence.

This paper aims to study the polarization of hard X-ray (HXR) sources in the solar atmosphere, including Compton backscattering of photons in the photosphere (the albedo effect) and the spatial distribution of polarization across the source. HXR photon polarization and spectra produced via electron-ion bremsstrahlung are calculated from electron distributions typical for solar flares. Compton scattering and photoelectric absorption are then modelled using Monte Carlo simulations of photon transport in the photosphere. Polarization maps across HXR sources (primary and albedo components) for each of the modelled electron distributions are calculated at various source locations from the solar centre to the limb. We show that Compton scattering produces a distinct polarization variation across the albedo patch at peak albedo energies of 20-50 keV for all anisotropies modelled. The results show that there are distinct spatial polarization changes in both the radial and perpendicular to radial directions across the extent of the HXR source at a given disk location. In the radial direction, the polarization magnitude and direction at specific positions along the HXR source will either increase or decrease with increased photon distribution directivity towards the photosphere. We also show how high electron cutoff energies influence the direction of polarization at above ~100 keV. Spatially resolved HXR polarization measurements can provide important information about the directivity and energetics of the electron distribution. Our results indicate the preferred angular resolution of polarization measurements required to distinguish between the scattered and primary components. We also show how spatially resolved polarization measurements could be used to probe the emission pattern of an HXR source, using both the magnitude and the direction of polarization.

The newly developed X-ray visibility forward fitting technique is applied to Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) data of a limb flare to investigate the energy and height dependence on sizes, shapes, and position of hard X-ray chromospheric footpoint sources. This provides information about the electron transport and chromospheric density structure. The spatial distribution of two footpoint X-ray sources is analyzed using PIXON, Maximum Entropy Method, CLEAN and visibility forward fit algorithms at nonthermal energies from $\sim 20$ to $\sim 200$ keV. We report, for the first time, the vertical extents and widths of hard X-ray chromospheric sources measured as a function of energy for a limb event. Our observations suggest that both the vertical and horizontal sizes of footpoints are decreasing with energy. Higher energy emission originates progressively deeper in the chromosphere consistent with downward flare accelerated streaming electrons. The ellipticity of the footpoints grows with energy from $\sim 0.5$ at $ \sim 20$ keV to $\sim 0.9$ at $\sim 150$ keV. The positions of X-ray emission are in agreement with an exponential density profile of scale height $\sim 150$~km. The characteristic size of the hard X-ray footpoint source along the limb is decreasing with energy suggesting a converging magnetic field in the footpoint. The vertical sizes of X-ray sources are inconsistent with simple collisional transport in a single density scale height but can be explained using a multi-threaded density structure in the chromosphere.

We investigate the positions and source sizes of X-ray solar flare sources taking into account Compton backscattering (albedo). Using a Monte Carlo simulation of X-ray photon transport including photo-electric absorption and Compton scattering, we calculate the apparent source sizes and positions of X-ray sources at the solar disk for various source sizes, spectral indices and directivities of the primary source. We show that the albedo effect will alter the true source positions and substantially increase the measured source sizes. The source positions are shifted up to $\sim 0.5"$ radially towards the disk centre and 5 arcsecond source sizes can be two times larger even for an isotropic source (minimum albedo effect) at 1 Mm above the photosphere. X-ray sources therefore should have minimum observed sizes, thus FWHM source size (2.35 times second-moment) will be as large as $\sim 7"$ in the 20-50 keV range for a disk-centered point source at a height of 1 Mm ($\sim 1.4"$) above the photosphere. The source size and position change is the largest for flatter primary X-ray spectra, stronger downward anisotropy, for sources closer to the solar disk centre, and between the energies of 30 and 50 keV. Albedo should be taken into account when X-ray footpoint positions, footpoint motions or source sizes from e.g. RHESSI or Yohkoh data are interpreted, and suggest that footpoint sources should be larger in X-rays than in optical or EUV ranges.