Context. Solar Energetic Particles (SEPs) with energy in the GeV range can propagate to Earth from their acceleration region near the Sun and produce Ground Level Enhancements (GLEs). The traditional approach to interpreting and modelling GLE observations assumes particle propagation only parallel to the magnetic field lines of interplanetary space, i.e. it is spatially 1D. Recent measurements by PAMELA have characterised SEP properties at 1 AU for the ~100 MeV-1 GeV range at high spectral resolution. Aims. We model the transport of GLE-energy solar protons through the Interplanetary Magnetic Field (IMF) using a 3D approach, to assess the effect of the Heliospheric Current Sheet (HCS) and drifts associated to the gradient and curvature of the Parker spiral. The latter are influenced by the IMF polarity. We derive 1 AU observables and compare the simulation results with data from PAMELA. Methods. We use a 3D test particle model including a HCS. Monoenergetic populations are studied first to obtain a qualitative picture of propagation patterns and numbers of crossings of the 1 AU sphere. Simulations for power law injection are used to derive intensity profiles and fluence spectra at 1 AU. A simulation for a specific event, GLE 71, is used to compare with PAMELA data. Results. Spatial patterns of 1 AU crossings and the average number of crossings are strongly influenced by 3D effects, with significant differences between periods of A+ and A- polarities. The decay time constant of 1 AU intensity profiles varies depending on the polarity and position of the observer, and it is not a simple function of the mean free path as in 1D models. Energy dependent leakage from the injection flux tube is particularly important for GLE energy particles, in many cases resulting in a roll-over in the fluence spectrum.

Axions/axion-like particles (ALPs) are a well motivated extension of the Standard Model and are generic within String Theory. The X-ray transparency of the intracluster medium (ICM) in galaxy clusters is a powerful probe of light ALPs (with mass $<10^{-11}\,{\rm eV}$); as X-ray photons from an embedded or background source propagate through the magnetized ICM, they may undergo energy-dependent quantum mechanical conversion into ALPs (and vice versa), imprinting distortions on the X-ray spectrum. We present Chandra data for the active galactic nucleus NGC1275 at the center of the Perseus cluster. Employing a 490ks High-Energy Transmission Gratings (HETG) exposure, we obtain a high-quality 1-9keV spectrum free from photon pileup and ICM contamination. Apart from iron-band features, the spectrum is described by a power-law continuum, with any spectral distortions at the $<3\%$ level. We compute photon survival probabilities as a function of ALP mass $m_a$ and ALP-photon coupling constant $g_{a\gamma}$ for an ensemble of ICM magnetic field models, and then use the NGC1275 spectrum to constraint the $(m_a, g_{a\gamma})$-plane. Marginalizing over magnetic field realizations, the 99.7% credible region limits the ALP-photon coupling to $g_{a\gamma}<6-8\times 10^{-13}\, {\rm GeV}^{-1}$ (depending upon magnetic field model) for masses $m_a<1\times 10^{-12}\,{\rm eV}$. These are the most stringent limit to date on $g_{a\gamma}$ for these light ALPs, and have already reached the sensitivity limits of next-generation helioscopes and light-shining-through-wall experiments. We highlight the potential of these studies with the next-generation X-ray observatories Athena and Lynx, but note the critical importance of advances in relative calibration of these future X-ray spectrometers.

Inflation can be supported in very steep potentials if it is generated by rapidly turning fields, which can be natural in negatively curved field spaces. The curvature perturbation, $\zeta$, of these models undergoes an exponential, transient amplification around the time of horizon crossing, but can still be compatible with observations at the level of the power spectrum. However, a recent analysis (based on a proposed single-field effective theory with an imaginary speed of sound) found that the trispectrum and other higher-order, non-Gaussian correlators also undergo similar exponential enhancements. This arguably leads to `hyper-large' non-Gaussianities in stark conflict with observations, and even to the loss of perturbative control of the calculations. In this paper, we provide the first analytic solution of the growth of the perturbations in two-field rapid-turn models, and find it in good agreement with previous numerical and single-field EFT estimates. We also show that the nested structure of commutators of the in-in formalism has subtle and crucial consequences: accounting for these commutators, we show analytically that the naively leading-order piece (which indeed is exponentially large) cancels exactly in all relevant correlators. The remaining non-Gaussianities of these models are modest, and there is no problem with perturbative control from the exponential enhancement of $\zeta$. Thus, rapid-turn inflation with negatively curved field spaces remains a viable and interesting class of candidate theories of the early universe.

In negatively curved field spaces, inflation can be realised even in steep potentials. Hyperinflation invokes the `centrifugal force' of a field orbiting the hyperbolic plane to sustain inflation. We generalise hyperinflation by showing that it can be realised in models with any number of fields ($N_f\geq2$), and in broad classes of potentials that, in particular, don't need to be rotationally symmetric. For example, hyperinflation can follow a period of radial slow-roll inflation that undergoes geometric destabilisation, yet this inflationary phase is not identical to the recently proposed scenario of `side-tracked inflation'. We furthermore provide a detailed proof of the attractor mechanism of (the original and generalised) hyperinflation, and provide a novel set of characteristic, explicit models. We close by discussing the compatibility of hyperinflation with observations and the recently much discussed `swampland conjectures'. Observationally viable models can be realised that satisfy either the `de Sitter conjecture' ($V'/V\gtrsim 1$) or the `distance conjecture' ($\Delta \phi \lesssim 1$), but satisfying both simultaneously brings hyperinflation in some tension with successful reheating after inflation. However, hyperinflation can get much closer to satisfying all of these criteria than standard slow-roll inflation. Furthermore, while the original model is in stark tension with the weak gravity conjecture, generalisations can circumvent this issue.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

We consider the robustness of small-field inflation in the presence of scalar field inhomogeneities. Previous numerical work has shown that if the scalar potential is flat only over a narrow interval, such as in commonly considered inflection-point models, even small-amplitude inhomogeneities present at the would-be onset of inflation at $\tau = \tau_i$ can disrupt the accelerated expansion. In this paper, we parametrise and evolve the inhomogeneities from an earlier time $\tau_{IC}$ at which the initial data were imprinted, and show that for a broad range of inflationary and pre-inflationary models, inflection-point inflation withstands initial inhomogeneities. We consider three classes of perturbative pre-inflationary solutions (corresponding to energetic domination by the scalar field kinetic term, a relativistic fluid, and isotropic negative curvature), and two classes of exact solutions to Einstein's equations with large inhomogeneities (corresponding to a stiff fluid with cylindrical symmetry, and anisotropic negative curvature). We derive a stability condition that depends on the Hubble scales $H(\tau_ i)$ and $H(\tau_{IC})$, and a few properties of the pre-inflationary cosmology. For initial data imprinted at the Planck scale, the absence of an inhomogeneous initial data problem for inflection-point inflation leads to a novel, lower limit on the tensor-to-scalar ratio.

It has recently been conjectured that string theory does not admit de Sitter vacua, and that quintessence explains the current epoch of accelerated cosmic expansion. A proposed, key prediction of this scenario is time-varying couplings in the dark sector, induced by the evolving quintessence field. We note that cosmological models with varying couplings suffer from severe problems with quantum corrections, beyond those shared by all quintessence models. The vacuum energy depends on all masses and couplings of the theory, and even small variations of parameters can lead to overwhelmingly large corrections to the effective potential. We find that quintessence models with varying parameters can be realised in consistent quantum theories by either: 1) enforcing exceptional levels of fine-tuning; 2) realising some unknown mechanism that cancels all undesirable contributions to the effective potential with unprecedented accuracy; or 3) ensuring that the quintessence field couples exclusively to very light states, and does not backreact on heavy fields.

Many applications of Gaussian random fields and Gaussian random processes are limited by the computational complexity of evaluating the probability density function, which involves inverting the relevant covariance matrix. In this work, we show how that problem can be completely circumvented for the local Taylor coefficients of a Gaussian random field with a Gaussian (or `square exponential') covariance function. Our results hold for any dimension of the field and to any order in the Taylor expansion. We present two applications. First, we show that this method can be used to explicitly generate non-trivial potential energy landscapes with many fields. This application is particularly useful when one is concerned with the field locally around special points (e.g.~maxima or minima), as we exemplify by the problem of cosmic `manyfield' inflation in the early universe. Second, we show that this method has applications in machine learning, and greatly simplifies the regression problem of determining the hyperparameters of the covariance function given a training data set consisting of local Taylor coefficients at single point. An accompanying Mathematica notebook is available at https://doi.org/10.17863/CAM.22859 .

Forecasting Solar Energetic Particle (SEP) fluence, as integrated over an SEP event, is an important element when estimating the effect of solar eruptions on humans and technology in space. Current real-time estimates are based on SEP measurements at a single location in space. However, the interplanetary magnetic field corotates with the Sun approximately 13$^\circ$ each day with respect to Earth, thus in 4 days a near-Earth spacecraft will have changed their connection about 60$^\circ$ from the original SEP source. We estimate the effect of the corotation on particle fluence using a simple particle transport model, and show that ignoring corotation can cause up to an order of magnitude error in fluence estimations, depending on the interplanetary particle transport conditions. We compare the model predictions with STEREO observations of SEP events.

We construct models of inflation with many randomly interacting fields and use these to study the generation of cosmological observables. We model the potentials as multi-dimensional Gaussian random fields (GRFs) and identify powerful algebraic simplifications that, for the first time, make it possible to access the manyfield limit of inflation in GRF potentials. Focussing on small-field, slow-roll, approximate saddle-point inflation in potentials with structure on sub-Planckian scales, we construct explicit examples involving up to 100 fields and generate statistical ensembles comprising of 164,000 models involving 5 to 50 fields. For the subset of these that support at least sixty e-folds of inflation, we use the 'transport method' and $\delta N$ formalism to determine the predictions for cosmological observables at the end of inflation, including the power spectrum and the local non-Gaussianities of the primordial perturbations. We find three key results: i) Planck compatibility is not rare, but future experiments may rule out this class of models; ii) In the manyfield limit, the predictions from these models agree well with, but are sharper than, previous results derived using potentials constructed through non-equilibrium Random Matrix Theory; iii) Despite substantial multifield effects, non-Gaussianities are typically very small: $f_{\rm nl}^{\rm loc} \ll 1$. We conclude that many of the 'generic predictions' of single-field inflation can be emergent features of complex inflation models.

We study inflation in models with many interacting fields subject to randomly generated scalar potentials. We use methods from non-equilibrium random matrix theory to construct the potentials and an adaption of the 'transport method' to evolve the two-point correlators during inflation. This construction allows, for the first time, for an explicit study of models with up to 100 interacting fields supporting a period of 'approximately saddle-point' inflation. We determine the statistical predictions for observables by generating over 30,000 models with 2-100 fields supporting at least 60 efolds of inflation. These studies lead us to seven lessons: i) Manyfield inflation is not single-field inflation, ii) The larger the number of fields, the simpler and sharper the predictions, iii) Planck compatibility is not rare, but future experiments may rule out this class of models, iv) The smoother the potentials, the sharper the predictions, v) Hyperparameters can transition from stiff to sloppy, vi) Despite tachyons, isocurvature can decay, vii) Eigenvalue repulsion drives the predictions. We conclude that many of the 'generic predictions' of single-field inflation can be emergent features of complex inflation models.

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

Axion-like particles (ALPs) and photons can quantum mechanically interconvert when propagating through magnetic fields, and ALP-photon conversion may induce oscillatory features in the spectra of astrophysical sources. We use deep (370 ks), short frame time Chandra observations of the bright nucleus at the centre of the radio galaxy M87 in the Virgo cluster to search for signatures of light ALPs. The absence of substantial irregularities in the X-ray power-law spectrum leads to a new upper limit on the photon-ALP coupling, $g_{a\gamma}$: using a conservative model of the cluster magnetic field consistent with Faraday rotation measurements from M87 and M84, we find $g_{a \gamma} < 2.6\times10^{-12}$ GeV$^{-1}$ at 95% confidence level for ALP masses $m_a \leq 10^{-13}$ eV. Other consistent magnetic field models lead to stronger limits of $g_{a \gamma} \lesssim 1.1$--$1.5 \times 10^{-12}$ GeV$^{-1}$. These bounds are all stronger than the limit inferred from the absence of a gamma-ray burst from SN1987A, and rule out a substantial fraction of the parameter space accessible to future experiments such as ALPS-II and IAXO.

Solar Energetic Particles (SEPs), a major component of space weather, propagate through the interplanetary medium strongly guided by the Interplanetary Magnetic Field (IMF). In this work, we analyse the implications a flat Heliospheric Current Sheet (HCS) has on proton propagation from SEP release sites to the Earth. We simulate proton propagation by integrating fully 3-D trajectories near an analytically defined flat current sheet, collecting comprehensive statistics into histograms, fluence maps and virtual observer time profiles within an energy range of 1--800 MeV. We show that protons experience significant current sheet drift to distant longitudes, causing time profiles to exhibit multiple components, which are a potential source of confusing interpretation of observations. We find that variation of current sheet thickness within a realistic parameter range has little effect on particle propagation. We show that IMF configuration strongly affects deceleration of protons. We show that in our model, the presence of a flat equatorial HCS in the inner heliosphere limits the crossing of protons into the opposite hemisphere.

Future cosmological surveys will probe the expansion history of the universe and constrain phenomenological models of dark energy. Such models do not address the fine-tuning problem of the vacuum energy, i.e. the cosmological constant problem (c.c.p.), but can make it spectacularly worse. We show that this is the case for 'interacting dark energy' models in which the masses of the dark matter states depend on the dark energy sector. If realised in nature, these models have far-reaching implications for proposed solutions to the c.c.p. that require the number of vacua to exceed the fine-tuning of the vacuum energy density. We show that current estimates of the number of flux vacua in string theory, $N_{\rm vac} \sim {\cal O}(10^{272,000})$, is far too small to realise certain simple models of interacting dark energy \emph{and} solve the cosmological constant problem anthropically. These models admit distinctive observational signatures that can be targeted by future gamma-ray observatories, hence making it possible to observationally rule out the anthropic solution to the cosmological constant problem in theories with a finite number of vacua.

Context. The intensity profiles of iron and oxygen in Solar Energetic Particle (SEP) events often display differences that result in a decreasing Fe/O ratio over time. The physical mechanisms behind this behaviour are not fully understood, but these observational signatures provide important tests of physical modelling efforts. Aims. In this paper we study the propagation of iron and oxygen SEP ions using a 3D model of propagation which includes the effect of guiding centre drift in a Parker spiral magnetic field. We derive time intensity profiles for a variety of observer locations and study the temporal evolution of the Fe/O ratio. Methods. We use a 3D full orbit test particle model which includes scattering. The configuration of the interplanetary magnetic field is a unipolar Parker spiral. Particles are released instantaneously from a compact region at 2 solar radii and allowed to propagate in 3D. Results. Both Fe and O experience significant transport across the magnetic field due to gradient and curvature drifts. We find that Fe ions drift more than O ions due to their larger mass-to-charge ratio, so that an observer that is not magnetically well connected to the source region will observe Fe arriving before O, for particles within the same range in energy per nucleon. As a result, for the majority of observer locations, the Fe/O ratio displays a decrease in time. Conclusions. We conclude that propagation effects associated with drifts produce a decay over time of the Fe/O ratio, qualitatively reproducing that observed in SEP event profiles.

The event-averaged charge state of heavy ion Solar Energetic Particles (SEPs), measured at 1 AU from the Sun, typically increases with the ions' kinetic energy. The origin of this behaviour has been ascribed to processes taking place within the acceleration region. In this paper we study the propagation through interplanetary space of SEP Fe ions, injected near the Sun with a variety of charge states that are uniformly distributed in energy, by means of a 3D test particle model. In our simulations, due to gradient and curvature drifts associated with the Parker spiral magnetic field, ions of different charge propagate with very different efficiencies to an observer that is not magnetically well connected to the source region. As a result we find that, for many observer locations, the 1 AU event-averaged charge state <Q>, as obtained from our model, displays an increase with particle energy E, in qualitative agreement with spacecraft observations. We conclude that drift-associated propagation is a possible explanation for the observed distribution of <Q> versus E in SEP events, and that the distribution measured in interplanetary space cannot be taken to represent that at injection.

We construct ensembles of random scalar potentials for $N_f$ interacting scalar fields using non-equilibrium random matrix theory, and use these to study the generation of observables during small-field inflation. For $N_f={\cal O}({\rm few})$, these heavily featured scalar potentials give rise to power spectra that are highly non-linear, at odds with observations. For $N_f\gg 1$, the superhorizon evolution of the perturbations is generically substantial, yet the power spectra simplify considerably and become more predictive, with most realisations being well approximated by a linear power spectrum. This provides proof of principle that complex inflationary physics can give rise to simple emergent power spectra. We explain how these results can be understood in terms of large $N_f$ universality of random matrix theory.

Context. Current solar energetic particle (SEP) propagation models describe the effects of interplanetary plasma turbulence on SEPs as diffusion, using a Fokker-Planck (FP) equation. However, FP models cannot explain the observed fast access of SEPs across the average magnetic field to regions that are widely separated in longitude within the heliosphere without using unrealistically strong cross-field diffusion. Aims. We study whether the recently suggested early non-diffusive phase of SEP propagation can explain the wide SEP events with realistic particle transport parameters. Methods. We used a novel model that accounts for the SEP propagation along field lines that meander as a result of plasma turbulence. Such a non-diffusive propagation mode has been shown to dominate the SEP cross-field propagation early in the SEP event history. We compare the new model to the traditional approach, and to SEP observations. Results. Using the new model, we reproduce the observed longitudinal extent of SEP peak fluxes that are characterised by a Gaussian profile with $\sigma=30-50^\circ$, while current diffusion theory can only explain extents of 11$^\circ$ with realistic diffusion coefficients. Our model also reproduces the timing of SEP arrival at distant longitudes, which cannot be explained using the diffusion model. Conclusions. The early onset of SEPs over a wide range of longitudes can be understood as a result of the effects of magnetic field-line random walk in the interplanetary medium and requires an SEP transport model that properly describes the non-diffusive early phase of SEP cross-field propagation.

Axion-like particles (ALPs) and photons inter-convert in the presence of a magnetic field. At keV energies in the environment of galaxy clusters, the conversion probability can become unsuppressed for light ALPs. Conversion of thermal X-ray photons into ALPs can introduce a step-like feature into the cluster thermal bremsstrahlung spectrum, and we argue that existing X-ray data on galaxy clusters should be sufficient to extend bounds on ALPs in the low-mass region $m_a \lesssim 1 \times 10^{-12}\,{\rm eV}$ down to $M \sim 7\times 10^{11}\, {\rm GeV}$, and that for $10^{11}\, {\rm GeV} < M \lesssim 10^{12}$ GeV light ALPs give rise to interesting and unique observational signatures that may be probed by existing and upcoming X-ray (and potentially X-ray polarisation) observations of galaxy clusters.

We investigate the deceleration of Solar Energetic Particles (SEPs) during their propagation from the Sun through interplanetary space, in the presence of weak to strong scattering in a Parker spiral configuration, using relativistic full orbit test particle simulations. The calculations retain all three spatial variables describing particles' trajectories, allowing to model any transport across the magnetic field. Large energy change is shown to occur for protons, due to the combined effect of standard adiabatic deceleration and a significant contribution from particle drift in the direction opposite to that of the solar wind electric field. The latter drift-induced deceleration is found to have a stronger effect for SEP energies than for galactic cosmic rays. The kinetic energy of protons injected at 1 MeV is found to be reduced by between 35 and 90% after four days, and for protons injected at 100 MeV by between 20 and 55%. The overall degree of deceleration is a weak function of the scattering mean free path, showing that, although adiabatic deceleration plays a role, a large contribution is due to particle drift. Current SEP transport models are found to account for drift-induced deceleration in an approximate way and their accuracy will need to be assessed in future work.

The capability to predict the parameters of an SEP event such as its onset, peak flux, and duration is critical to assessing any potential space weather impact. We present a new flexible modelling system simulating the propagation of Solar Energetic Particles (SEPs) from locations near the Sun to any given location in the heliosphere to forecast the SEP flux profiles. SPARX uses an innovative methodology, that allows implementation within an operational framework, to overcome the time constraints of test particle modelling of SEP profiles, allowing the production of near real time SEP nowcasts and forecasts, when paired with appropriate near real time triggers. SPARX has the capability to produce SEP forecasts within minutes of being triggered by observations of a solar eruptive event. The model is based on the test particle approach and is spatially 3D, thus allowing for the possibility of transport in the direction perpendicular to the magnetic field. The model naturally includes the effects of perpendicular propagation due to drifts and drift-induced deceleration. The modelling framework and the way in which parameters of relevance for Space Weather forecasting are obtained are described. The first results from the modelling system are presented. These results demonstrate that corotation and drift of SEP streams play an important role in shaping SEP flux profiles.

A statistical analysis of the relationship between solar energetic particles (SEPs) and properties of solar flares and coronal mass ejections (CMEs) is presented. SEP events during solar cycle 23 are selected which are associated with solar flares originating on the visible hemisphere of the Sun and at least of magnitude M1. Taking into account all flares and CMEs that occurred during this period, the probability for the occurrence of an SEP event near Earth is determined. A strong rise of this probability is observed for increasing flare intensities, more western locations, larger CME speeds and halo CMEs. The correlations between the proton peak flux and these solar parameters are derived for a low (>10 MeV) and high (>60 MeV) energy range excluding any flux enhancement due to the passage of fast interplanetary shocks. The obtained correlation coefficients are: 0.55+-0.07 (0.63+-0.06) with flare intensity and 0.56+-0.08 (0.40+-0.09) with the CME speed for E>10 MeV (E>60 MeV). For both energy ranges, the correlations with flare longitude and CME width are very small or non-existent. Furthermore, the occurrence probabilities, correlation coefficients and mean peak fluxes are derived in multi-dimensional bins combining the aforementioned solar parameters. The correlation coefficients are also determined in different proton energy channels ranging from 5 to 200 MeV. The results show that the correlation between the proton peak flux and the CME speed decreases with energy, while the correlation with the flare intensity shows the opposite behavior. Furthermore, the correlation with the CME speed is stronger than the correlation with the flare intensity below 15 MeV and becomes weaker above 20 MeV. Excluding the flux enhancements due to interplanetary shocks, only a small but not very significant change is observed in the correlation between the peak flux below 7 MeV and the CME speed.

We show that the scattering of a general spin-0 sector of dark radiation off the pre-recombination thermal plasma results in undetectably small spectral distortions of the Cosmic Microwave Background.

Relativistic axions are good candidates for the dark radiation for which there are mounting observational hints. The primordial decays of heavy fields produce axions which are ultra-energetic compared to thermalised matter and inelastic axion-matter scattering can occur with $E_{CoM} \gg T_{\gamma}$, thus accessing many interesting processes which are otherwise kinematically forbidden in standard cosmology. Axion-photon scattering into quarks and leptons during BBN affects the light element abundances, and bounds on overproduction of $^4$He constrain a combination of the axion decay constant and the reheating temperature. For supersymmetric models, axion scattering into visible sector superpartners can give direct non-thermal production of dark matter at $T_{\gamma} \ll T_{freezeout}$. Most axions --- or any other dark radiation candidate from modulus decay --- still linger today as a Cosmic Axion Background with $E_{axion} \sim \mathcal{O}(100) eV$, and a flux of $\sim 10^6 cm^{-2} s^{-1}$.

We construct a class of random potentials for N >> 1 scalar fields using non-equilibrium random matrix theory, and then characterize multifield inflation in this setting. By stipulating that the Hessian matrices in adjacent coordinate patches are related by Dyson Brownian motion, we define the potential in the vicinity of a trajectory. This method remains computationally efficient at large N, permitting us to study much larger systems than has been possible with other constructions. We illustrate the utility of our approach with a numerical study of inflation in systems with up to 100 coupled scalar fields. A significant finding is that eigenvalue repulsion sharply reduces the duration of inflation near a critical point of the potential: even if the curvature of the potential is fine-tuned to be small at the critical point, small cross-couplings in the Hessian cause the curvature to grow in the neighborhood of the critical point.

Motivated by the possibility of explaining the 3.5 keV line through dark matter decaying to axion-like particles that subsequently convert to photons, we study ALP-photon conversion for sightlines passing within 50 pc of the galactic centre. Conversion depends on the galactic centre magnetic field which is highly uncertain. For fields at low or mid-range of observational estimates (10--100 $\mu$G), no observable signal is possible. For fields at the high range of observational estimates (a pervasive poloidal mG field over the central 150 pc) it is possible to generate sufficient signal to explain recent observations of a 3.5 keV line in the galactic centre. In this scenario, the galactic centre line signal comes predominantly from the region with $z > 20$ pc, reconciling the results from the Chandra and XMM-Newton X-ray telescopes. The dark matter to ALP to photon scenario also naturally predicts the non-observation of the 3.5 keV line in stacked galaxy spectra. We further explore predictions for the line flux in galaxies and suggest a set of galaxies that is optimised for observing the 3.5 keV line in this model.

Primordial decays of string theory moduli at z \sim 10^{12} naturally generate a dark radiation Cosmic Axion Background (CAB) with 0.1 - 1 keV energies. This CAB can be detected through axion-photon conversion in astrophysical magnetic fields to give quasi-thermal excesses in the extreme ultraviolet and soft X-ray bands. Substantial and observable luminosities may be generated even for axion-photon couplings \ll 10^{-11} GeV^{-1}. We propose that axion-photon conversion may explain the observed excess emission of soft X-rays from galaxy clusters, and may also contribute to the diffuse unresolved cosmic X-ray background. We list a number of correlated predictions of the scenario.

Galaxy clusters can efficiently convert axion-like particles (ALPs) to photons. We propose that the recently claimed detection of a 3.55--3.57 keV line in the stacked spectra of a large number of galaxy clusters and the Andromeda galaxy may originate from the decay of either a scalar or fermionic $7.1$ keV dark matter species into an axion-like particle (ALP) of mass $m_{a} \lesssim 6\cdot 10^{-11}~{\rm eV}$, which subsequently converts to a photon in the cluster magnetic field. In contrast to models in which the photon line arises directly from dark matter decay or annihilation, this can explain the anomalous line strength in the Perseus cluster. As axion-photon conversion scales as $B^2$ and cool core clusters have high central magnetic fields, this model can also explains the observed peaking of the line emission in the cool cores of the Perseus, Ophiuchus and Centaurus clusters, as opposed to the much larger dark matter halos. We describe distinctive predictions of this scenario for future observations.

We show that the soft X-ray excess in the Coma cluster can be explained by a cosmic background of relativistic axions converting into photons in the cluster magnetic field. We provide a detailed self-contained review of the cluster soft X-ray excess, the proposed astrophysical explanations and the problems they face, and explain how a 0.1-1 keV axion background naturally arises at reheating in many string theory models of the early universe. We study the morphology of the soft excess by numerically propagating axions through stochastic, multi-scale magnetic field models that are consistent with observations of Faraday rotation measures from Coma. By comparing to ROSAT observations of the 0.2-0.4 keV soft excess, we find that the overall excess luminosity is easily reproduced for $g_{a\gamma\gamma} \sim 2 \times 10^{-13}$ GeV$^{-1}$. The resulting morphology is highly sensitive to the magnetic field power spectrum. For Gaussian magnetic field models, the observed soft excess morphology prefers magnetic field spectra with most power in coherence lengths on ${\cal O}(3 {\rm ~kpc})$ scales over those with most power on ${\cal O}(12 {\rm ~kpc})$ scales. Within this scenario, we bound the mean energy of the axion background to $50\, {\rm eV}\lesssim \langle E_a \rangle \lesssim 250\, {\rm eV}$, the axion mass to $m_a \lesssim 10^{-12}\,\hbox{eV}$, and derive a lower bound on the axion-photon coupling $g_{a\gamma\gamma} \gtrsim \sqrt{0.5/\Delta N_{\rm eff}}\, 1.4 \times 10^{-13}$ GeV$^{-1}$.

Drifts are known to play a role in galactic cosmic ray transport within the heliosphere and are a standard component of cosmic ray propagation models. However, the current paradigm of Solar Energetic Particle (SEP) propagation holds the effects of drifts to be negligible, and they are not accounted for in most current SEP modelling efforts. We present full-orbit test particle simulations of SEP propagation in a Parker spiral interplanetary magnetic field which demonstrate that high energy particle drifts cause significant asymmetric propagation perpendicular to the interplanetary magnetic field. Thus in many cases the assumption of field aligned propagation of SEPs may not be valid. We show that SEP drifts have dependencies on energy, heliographic latitude, and charge to mass ratio, that are capable of transporting energetic particles perpendicular to the field over significant distances within interplanetary space, e.g. protons of initial energy 100 MeV propagate distances across the field on the order of 1 AU, over timescales typical of a gradual SEP event. Our results demonstrate the need for current models of SEP events to include the effects of particle drift. We show that the drift is considerably stronger for heavy ion SEPs due to their larger mass to charge ratio. This paradigm shift has important consequences for the modelling of SEP events and is crucial to the understanding and interpretation of in-situ observations.

Solar energetic particles (SEPs) have been observed to easily spread across heliographic longitudes, and the mechanisms responsible for this behaviour remain unclear. We use full-orbit simulations of a 10 MeV proton beam in a turbulent magnetic field to study to what extent the spread across the mean field can be described as diffusion early in a particle event. We compare the full-orbit code results to solutions of a Fokker-Planck equation including spatial and pitch angle diffusion, and of one including also propagation of the particles along random-walking magnetic field lines. We find that propagation of the particles along meandering field lines is the key process determining their cross-field spread at 1 AU at the beginning of the simulated event. The mean square displacement of the particles an hour after injection is an order of magnitude larger than that given by the diffusion model, indicating that models employing spatial cross-field diffusion cannot be used to describe early evolution of an SEP event. On the other hand, the diffusion of the particles from their initial field lines is negligible during the first 5 hours, which is consistent with the observations of SEP intensity dropouts. We conclude that modelling SEP events must take into account the particle propagation along meandering field lines for the first 20 hours of the event.

Drifts in the Parker spiral interplanetary magnetic field are known to be an important component in the propagation of galactic cosmic rays, while they are thought to be negligible for Solar Energetic Particles (SEPs). As a result they have so far been ignored in SEP propagation modelling and data analysis. We examine drift velocities in the Parker spiral within single particle first-order adiabatic theory, in a local coordinate system with an axis parallel to the magnetic field. We show that, in the presence of scattering in interplanetary space, protons at the high end of the SEP energy range experience significant gradient and curvature drift. In the scatter-free case, drift due to magnetic field curvature is present. The magnitude of drift velocity increases by more than an order of magnitude at high heliographic latitudes compared to near the ecliptic; it has a strong dependence on radial distance r from the Sun, reaching a maximum at r~1 AU at low heliolatitudes and r~10 AU at high heliolatitudes. Due to the mass over charge dependence of drift velocities, the effect of drift for partially ionised SEP heavy ions is stronger than for protons. Drift is therefore likely to be a considerable source of cross field transport for high energy SEPs.

Observations and modelling suggest that the fluctuations in magnetised plasmas exhibit scale-dependent anisotropy, with more energy in the fluctuations perpendicular to the mean magnetic field than in the parallel fluctuations and the anisotropy increasing at smaller scales. The scale-dependence of the anisotropy has not been studied in full-orbit simulations of particle transport in turbulent plasmas so far. In this paper, we construct a model of critically balanced turbulence, as suggested by \cite{GoSr1995}, and calculate energetic particle spatial diffusion coefficients using full-orbit simulations. The model uses an enveloped turbulence approach, where each 2-dimensional wave mode with wavenumber $k_\perp$ is packed into envelopes of length $L$ following the critical balance condition, $L\propto k_\perp^{-2/3}$, with the wave mode parameters changing between envelopes. Using full-orbit particle simulations, we find that both the parallel and perpendicular diffusion coefficients increase by a factor 2, compared to previous models with scale-independent anisotropy.

Spectroscopic and stereoscopic imaging observations of slow magnetoacoustic wave propagation within a coronal loop are investigated to determine the decay length scale of the slow magnetoacoustic mode in three dimensions and the density profile within the loop system. The slow wave is found to have an e-folding decay length scale of $20,000^{+4000}_{-3000}$km with a uniform density profile along the loop base. These observations place quantitive constraints on the modelling of wave propagation within coronal loops. Theoretical forward modelling suggests that magnetic field line divergence is the dominant damping factor and thermal conduction is insufficient, given the observed parameters of the coronal loop temperature, density and wave mode period.

Recently observed oscillations in the solar atmosphere have been interpreted and modeled as magnetohydrodynamic wave modes. This has allowed the estimation of parameters that are otherwise hard to derive, such as the coronal magnetic-field strength. This work crucially relies on the initial detection of the oscillations, which is commonly done manually. The volume of Solar Dynamics Observatory (SDO) data will make manual detection inefficient for detecting all of the oscillating regions. An algorithm is presented which automates the detection of areas of the solar atmosphere that support spatially extended oscillations. The algorithm identifies areas in the solar atmosphere whose oscillation content is described by a single, dominant oscillation within a user-defined frequency range. The method is based on Bayesian spectral analysis of time-series and image filtering. A Bayesian approach sidesteps the need for an a-priori noise estimate to calculate rejection criteria for the observed signal, and it also provides estimates of oscillation frequency, amplitude and noise, and the error in all these quantities, in a self-consistent way. The algorithm also introduces the notion of quality measures to those regions for which a positive detection is claimed, allowing simple post-detection discrimination by the user. The algorithm is demonstrated on two Transition Region and Coronal Explorer (TRACE) datasets, and comments regarding its suitability for oscillation detection in SDO are made.

The Extreme-Ultraviolet Imaging Spectrometer on board the HINODE satellite is used to examine the loop system described in Marsh et al. (2009) by applying spectroscopic diagnostic methods. A simple isothermal mapping algorithm is applied to determine where the assumption of isothermal plasma may be valid, and the emission measure locii technique is used to determine the temperature profile along the base of the loop system. It is found that, along the base, the loop has a uniform temperature profile with a mean temperature of 0.89 +- 0.09 MK which is in agreement with the temperature determined seismologically in Marsh et al. (2009), using observations interpreted as the slow magnetoacoustic mode. The results further strengthen the slow mode interpretation, propagation at a uniform sound speed, and the analysis method applied in Marsh et al. (2009). It is found that it is not possible to discriminate between the slow mode phase speed and the sound speed within the precision of the present observations.

On 2008 January 10, the twin Solar Terrestrial Relations Observatory (STEREO) A and B spacecraft conducted a high time cadence study of the solar corona with the Extreme UltraViolet Imager (EUVI) instruments with the aim of investigating coronal dynamics. Observations of the three-dimensional propagation of waves within active region coronal loops and a measurement of the true coronal slow mode speed are obtained. Intensity oscillations with a period of approximately 12 minutes are observed to propagate outwards from the base of a loop system, consistent with the slow magnetoacoustic mode. A novel analysis technique is applied to measure the wave phase velocity in the observations of the A and B spacecraft. These stereoscopic observations are used to infer the three-dimensional velocity vector of the wave propagation, with an inclination of 37 +- 6 deg to the local normal and a magnitude of 132 +- 9 and 132 +- 11 km s-1, giving the first measurement of the true coronal longitudinal slow mode speed, and an inferred temperature of 0.84 +- 12 MK and 0.84 +- 15 MK.

A Bayesian probability based approach is applied to the problem of detecting and parameterizing oscillations in the upper solar atmosphere for the first time. Due to its statistical origin, this method provides a mechanism for determining the number of oscillations present, gives precise estimates of the oscillation parameters with a self-consistent statistical error analysis, and allows the oscillatory model signals to be reconstructed within these errors. A highly desirable feature of the Bayesian approach is the ability to resolve oscillations with extremely small frequency separations. The code is applied to SOHO/CDS (Solar and Heliospheric Observatory/Coronal Diagnostic Spectrometer) O V 629A observations and resolves four distinct P4, P5, P6 and P7 p-modes within the same sunspot transition region. This suggests that a spectrum of photospheric p-modes is able to propagate into the upper atmosphere of the Sun and Sun-like stars, and places precise observational constraints on models of umbral eigen modes.