D. Ruffolo

Measurement and simulation of neutron monitor count rate dependence on surrounding structure

Neutron monitors are the premier instruments for precise measurements of time variations (e.g., of solar origin) in the galactic cosmic ray (GCR) flux in the range of ∼1–100u2009GeV. However, it has proven challenging to accurately determine the yield function (effective area) versus rigidity in order to relate a neutron monitors count rate with those of other monitors worldwide and the underlying GCR spectrum. Monte Carlo simulations of the yield function have been developed, but there have been few opportunities to validate these observationally, especially regarding the particular environment surrounding each monitor. Here we have precisely measured the count rate of a calibration neutron monitor near the Princess Sirindhorn Neutron Monitor (PSNM) at Doi Inthanon, Thailand (18.59∘N, 98.49∘E, 2560u2009m altitude), which provides a basis for comparison with count rates of other neutron monitors worldwide that are similarly calibrated. We directly measured the effect of surrounding structure by operating the calibrator outside and inside the building. Using Monte Carlo simulations, we clarify differences in response of the calibrator and PSNM, as well as the calibrator outside and inside the building. The dependence of the calibrator count rate on surrounding structure can be attributed to its sensitivity to neutrons of 0.5–10u2009MeV and a shift of sensitivity to nucleons of higher energy when placed inside the building. Simulated calibrator to PSNM count rate ratios inside and outside agree with observations within a few percent, providing useful validation and improving confidence in our ability to model the yield function for a neutron monitor station.

Latitude Survey Investigation of Galactic Cosmic Ray Solar Modulation during 1994-2007

The Galactic cosmic ray spectrum exhibits subtle variations over the 22 yr solar magnetic cycle in addition to the more dramatic variations over the 11 yr sunspot cycle. Neutron monitors are large ground-based detectors that provide accurate measurements of variations in the cosmic ray flux at the top of the atmosphere above the detector. At any given location the magnetic field of the Earth excludes particles below a well-defined rigidity (momentum per unit charge) known as the cutoff rigidity, which can be accurately calculated using detailed models of the geomagnetic field. By carrying a neutron monitor to different locations, e.g., on a ship, the Earth itself serves as a magnet spectrometer. By repeating such latitude surveys with identical equipment, a sensitive measurement of changes in the spectrum can be made. In this work, we analyze data from the 1994 through 2007 series of latitude surveys conducted by the Bartol Research Institute, the University of Tasmania, and the Australian Antarctic Division. We confirm the curious crossover in spectra measured near solar minima during epochs of opposite solar magnetic polarity, and show that it is directly related to a sudden change in the spectral behavior of solar modulation at the time of the polarity reversal, as revealed from contemporaneous variations in the survey data and a fixed station. We suggest that the spectral change and crossover result from the interaction of effects due to gradient/curvature drifts with a systematic change in the interplanetary diffusion coefficient caused by turbulent magnetic helicity.

EFFECT OF COHERENT STRUCTURES ON ENERGETIC PARTICLE INTENSITY IN THE SOLAR WIND AT 1 AU

Solar energetic particles in the solar wind are accelerated in both solar flares and shocks assocated with fast coronal mass ejections. They follow the interplanetary magnetic field and, upon reaching Earth, have implications for space weather. Space weather affects astronaut health and orbiting equipment through radiation hazard and electrical infrastructure on the ground with ground induced currents. Economic impacts include disruption of GPS and redirection of commercial polar flights due to a dangerous radiation environment over the poles. By studying how these particles interact with the magnetic fields we can better predict onset times and diffusion of these events. We find, using superposed epoch analysis and conditional statisitics from spacecraft observations that there is a strong association between energetic particles in the solar wind and magnetic discontinuities. This may be related to turbulent dissipation mechanisms in which coherent structures in the solar wind seem to be preferred sites of heating, plasma instabilites and dissipation. In the case of energetic particles, magnetic reconnection and transport in flux tubes are likely to play a role. Though we focus on data away from large shocks, trapping can occur in the downstream region of shocks due to the preponderance of compressive turbulence in these areas. This thesis lays the ground work for the results described above with an introduction to solar wind and heliospheric physics in Chapter 1. Chapter 2 is an introduction to the acceleration mechanisms that give rise to observed energetic particle events. Chapter 3 describes various data analysis techniques and statistics that are bread and butter when analyzing spacecraft data for turbulence and energetic particle studies. Chapter 4 is a digression that covers preliminary studies that were done on the side; scale dependent kurtosis, ergodic studies and initial conditions for simulations. Chapter 5 contains that central published results of this thesis, that there is a strong

Monte Carlo simulation of the neutron monitor yield function

Neutron monitors (NMs) are ground-based detectors that measure variations of the Galactic cosmic ray flux at GV range rigidities. Differences in configuration, electronics, surroundings, and location induce systematic effects on the calculation of the yield functions of NMs worldwide. Different estimates of NM yield functions can differ by a factor of 2 or more. In this work, we present new Monte Carlo simulations to calculate NM yield functions and perform an absolute (not relative) comparison with the count rate of the Princess Sirindhorn Neutron Monitor (PSNM) at Doi Inthanon, Thailand, both for the entire monitor and for individual counter tubes. We model the atmosphere using profiles from the Global Data Assimilation System database and the Naval Research Laboratory Mass Spectrometer, Incoherent Scatter Radar Extended model. Using FLUKA software and the detailed geometry of PSNM, we calculated the PSNM yield functions for protons and alpha particles. An agreement better than 9% was achieved between the PSNM observations and the simulated count rate during the solar minimum of December 2009. The systematic effect from the electronic dead time was studied as a function of primary cosmic ray rigidity at the top of the atmosphere up to 1xa0TV. We show that the effect is not negligible and can reach 35% at high rigidity for a dead time >1xa0ms. We analyzed the response function of each counter tube at PSNM using its actual dead time, and we provide normalization coefficients between count rates for various tube configurations in the standard NM64 design that are valid to within ∼1% for such stations worldwide.

Dependence of the neutron monitor count rate and time delay distribution on the rigidity spectrum of primary cosmic rays

Neutron monitors are the premier instruments for precisely tracking time variations in the Galactic cosmic ray flux at GeV-range energies above the geomagnetic cutoff at the location of measurement. Recently, a new capability has been developed to record and analyze the neutron time delay distribution (related to neutron multiplicity) to infer variations in the cosmic ray spectrum as well. In particular, we can determine the leader fraction L, defined as the fraction of neutrons that did not follow a previous neutron detection in the same tube from the same atmospheric secondary particle, from time delay histograms. Using data taken during 2000-2007 by a ship-borne neutron monitor latitude survey, we observe strong dependences of the count rate and L on the geomagnetic cutoff. We have modeled this dependence using Monte Carlo simulations of cosmic ray interactions in the atmosphere and in the neutron monitor. We present new yield functions for the count rate of a neutron monitor at sea level. The simulation results show a variation of L with geomagnetic cutoff as observed by the latitude survey, confirming that these changes in L can be attributed to changes in the cosmic ray spectrum arriving at Earths atmosphere. We also observe a variation in L with time at a fixed cutoff, which reflects the evolution of the cosmic ray spectrum with the sunspot cycle, known as solar modulation.

COMPLEXITY AND DIFFUSION OF MAGNETIC FLUX SURFACES IN ANISOTROPIC TURBULENCE

The complexity of magnetic flux surfaces is investigated analytically and numerically in static homogeneous magnetic turbulence. Magnetic surfaces are computed to large distances in magnetic fields derived from a reduced magnetohydrodynamic model. The question addressed is whether one can define magnetic surfaces over large distances when turbulence is present. Using a flux surface spectral analysis, we show that magnetic surfaces become complex at small scales, experiencing an exponential thinning that is quantified here. The computation of a flux surface is of either exponential or nondeterministic polynomial complexity, which has the conceptual implication that global identification of magnetic flux surfaces and flux exchange, e.g., in magnetic reconnection, can be intractable in three dimensions. The coarse-grained large-scale magnetic flux experiences diffusive behavior. The link between the diffusion of the coarse-grained flux and field-line random walk is established explicitly through multiple scale analysis. The Kubo number controls both large and small scale limits. These results have consequences for interpreting processes such as magnetic reconnection and field-line diffusion in astrophysical plasmas.

Magnetic Field Line Random Walk in Isotropic Turbulence with Zero Mean Field

In astrophysical plasmas, magnetic field lines often guide the motions of thermal and non-thermal particles. The field line random walk (FLRW) is typically considered to depend on the Kubo number R = (b/B 0)(l∥/l⊥) for rms magnetic fluctuation b, large-scale mean field B 0, and parallel and perpendicular coherence scales l∥ and l⊥, respectively. Here we examine the FLRW when R → ∞ by taking B 0 → 0 for finite bz (fluctuation component along B 0), which differs from the well-studied route with bz = 0 or bz B 0 as the turbulence becomes quasi-two-dimensional (quasi-2D). Fluctuations with B 0 = 0 are typically isotropic, which serves as a reasonable model of interstellar turbulence. We use a non-perturbative analytic framework based on Corrsins hypothesis to determine closed-form solutions for the asymptotic field line diffusion coefficient for three versions of the theory, which are directly related to the k –1 or k –2 moment of the power spectrum. We test these theories by performing computer simulations of the FLRW, obtaining the ratio of diffusion coefficients for two different parameterizations of a field line. Comparing this with theoretical ratios, the random ballistic decorrelation version of the theory agrees well with the simulations. All results exhibit an analog to Bohm diffusion. In the quasi-2D limit, previous works have shown that Corrsin-based theories deviate substantially from simulation results, but here we find that as B 0 → 0, they remain in reasonable agreement. We conclude that their applicability is limited not by large R, but rather by quasi-two-dimensionality.

Charged particle diffusion in isotropic random magnetic fields

The investigation of the diffusive transport of charged particles in a turbulent magnetic field remains a subject of considerable interest. Research has most frequently concentrated on determining the diffusion coefficient in the presence of a mean magnetic field. Here we consider diffusion of charged particles in fully three-dimensional isotropic turbulent magnetic fields with no mean field, which may be pertinent to many astrophysical situations. We identify different ranges of particle energy depending upon the ratio of the Larmor radius of the charged particle to the characteristic outer length scale of the turbulence. Two different theoretical models are proposed to calculate the diffusion coefficient, each applicable to a distinct range of particle energies. The theoretical results are compared with those from computer simulations, showing good agreement.

NON-AXISYMMETRIC PERPENDICULAR DIFFUSION OF CHARGED PARTICLES AND THEIR TRANSPORT ACROSS TANGENTIAL MAGNETIC DISCONTINUITIES

We investigate the transport of charged particles across magnetic discontinuities, focusing specifically on stream interfaces associated with co-rotating interaction regions in the solar wind. We argue that the magnetic field fluctuations perpendicular to the magnetic discontinuity, and usually also perpendicular to the mean magnetic field, are strongly damped in the vicinity of such a magnetic structure, leading to anisotropic perpendicular diffusion. Assuming that perpendicular diffusion arises from drifts in a turbulent magnetic field, we adopt a simplified approach to derive the relevant perpendicular diffusion coefficient. This approach, which we believe gives the correct principal dependences as expected from more elaborate calculations, allows us to investigate transport in different turbulent geometries, such as longitudinal compressional turbulence that may be present near the heliopause. Although highly dependent on the (possibly anisotropic) perpendicular length scales and turbulence levels, we generally find perpendicular diffusion to be strongly damped at magnetic discontinuities, which may in turn provide an explanation for the large particle gradients associated with these structures.

Driving reconnection in sheared magnetic configurations with forced fluctuations

We investigate reconnection of magnetic field lines in sheared magnetic field configurations due to fluctuations driven by random forcing by means of numerical simulations. The simulations are performed with an incompressible, pseudo-spectral magnetohydrodynamics code in 2D where we take thick, resistively decaying, current-sheet like sheared magnetic configurations which do not reconnect spontaneously. We describe and test the forcing that is introduced in the momentum equation to drive fluctuations. It is found that the forcing does not change the rate of decay; however, it adds and removes energy faster in the presence of the magnetic shear structure compared to when it has decayed away. We observe that such a forcing can induce magnetic reconnection due to field line wandering leading to the formation of magnetic islands and O-points. These reconnecting field lines spread out as the current sheet decays with time. A semi-empirical formula is derived which reasonably explains the formation and spread o...