Particle acceleration by turbulent magnetohydro-dynamic reconnection

Test particles in a two dimensional, turbulent MHD simulation are found to undergo significant acceleration. The magnetic field configuration is a periodic sheet pinch which undergoes reconnection. The test particles are trapped in the reconnection region for times of order an Alfven transit time in the large electric fields that characterize the turbulent reconnection process at the relatively large magnetic Reynolds number used in the simulation. The maximum speed attained by these particles is consistent with an analytic estimate which depends on the reconnection electric field, the Alfven speed, and the ratio of Larmor period to the Alfven transit time

Charged-Particle Motion in Electromagnetic Fields Having at Least One Ignorable Spatial Coordinate

We give a rigorous derivation of a theorem showing that charged particles in an arbitrary electromagnetic field with at least one ignorable spatial coordinate remain forever tied to a given magnetic-field line. Such a situation contrasts the significant motions normal to the magnetic field that are expected in most real three-dimensional systems. It is pointed out that, while the significance of the theorem has not been widely appreciated, it has important consequences for a number of problems and is of particular relevance for the acceleration of cosmic rays by shocks.Comment: 7 pages, emulateapj format, including 1 eps figure, to appear in The Astrophysical Journal, Dec. 10 1998 issu

Particle Acceleration in Multiple Dissipation Regions

The sharp magnetic discontinuities which naturally appear in solar magnetic flux tubes driven by turbulent photospheric motions are associated with intense currents. \citet{Par83} proposed that these currents can become unstable to a variety of microscopic processes, with the net result of dramatically enhanced resistivity and heating (nanoflares). The electric fields associated with such ``hot spots'' are also expected to enhance particle acceleration. We test this hypothesis by exact relativistic orbit simulations in strong random phase magnetohydrodynamic (MHD) turbulence which is forming localized super-Dreicer Ohm electric fields ($E_\Omega/E_D$ = $10^2 ... 10^5$) occurring in 2..15 % of the volume. It is found that these fields indeed yield a large amplification of acceleration of electrons and ions, and can effectively overcome the injection problem. We suggest in this article that nanoflare heating will be associated with sporadic particle acceleration.Comment: 12 pages, 5 figures, to appear in ApJ

Fourier Acceleration of Langevin Molecular Dynamics

Fourier acceleration has been successfully applied to the simulation of lattice field theories for more than a decade. In this paper, we extend the method to the dynamics of discrete particles moving in continuum. Although our method is based on a mapping of the particles' dynamics to a regular grid so that discrete Fourier transforms may be taken, it should be emphasized that the introduction of the grid is a purely algorithmic device and that no smoothing, coarse-graining or mean-field approximations are made. The method thus can be applied to the equations of motion of molecular dynamics (MD), or its Langevin or Brownian variants. For example, in Langevin MD simulations our acceleration technique permits a straightforward spectral decomposition of forces so that the long-wavelength modes are integrated with a longer time step, thereby reducing the time required to reach equilibrium or to decorrelate the system in equilibrium. Speedup factors of up to 30 are observed relative to pure (unaccelerated) Langevin MD. As with acceleration of critical lattice models, even further gains relative to the unaccelerated method are expected for larger systems. Preliminary results for Fourier-accelerated molecular dynamics are presented in order to illustrate the basic concepts. Possible extensions of the method and further lines of research are discussed.Comment: 11 pages, two illustrations included using graphic

Development of a common data model for scientific simulations

The problem of sharing data among scientific simulation models is a difficult and persistent one. Computational scientists employ an enormous variety of discrete approximations in modeling physical processes on computers. Problems occur when models based on different representations are required to exchange data with one another, or with some other software package. Within the DOE`s Accelerated Strategic Computing Initiative (ASCI), a cross-disciplinary group called the Data Models and Formats (DMF) group, has been working to develop a common data model. The current model is comprised of several layers of increasing semantic complexity. One of these layers is an abstract model based on set theory and topology called the fiber bundle kernel (FBK). This layer provides the flexibility needed to describe a wide range of mesh-approximated functions as well as other entities. This paper briefly describes the ASCI common data model, its mathematical basis, and ASCI prototype development. These prototypes include an object-oriented data management library developed at Los Alamos called the Common Data Model Library or CDMlib, the Vector Bundle API from the Lawrence Livermore Laboratory, and the DMF API from Sandia National Laboratory

Particle Acceleration in an Evolving Network of Unstable Current Sheets

We study the acceleration of electrons and protons interacting with localized, multiple, small-scale dissipation regions inside an evolving, turbulent active region. The dissipation regions are Unstable Current Sheets (UCS), and in their ensemble they form a complex, fractal, evolving network of acceleration centers. Acceleration and energy dissipation are thus assumed to be fragmented. A large-scale magnetic topology provides the connectivity between the UCS and determines in this way the degree of possible multiple acceleration. The particles travel along the magnetic field freely without loosing or gaining energy, till they reach a UCS. In a UCS, a variety of acceleration mechanisms are active, with the end-result that the particles depart with a new momentum. The stochastic acceleration process is represented in the form of Continuous Time Random Walk (CTRW), which allows to estimate the evolution of the energy distribution of the particles. It is found that under certain conditions electrons are heated and accelerated to energies above 1 MeV in much less than a second. Hard X-ray (HXR) and microwave spectra are calculated from the electrons' energy distributions, and they are found to be compatible with the observations. Ions (protons) are also heated and accelerated, reaching energies up to 10 MeV almost simultaneously with the electrons. The diffusion of the particles inside the active region is extremely fast (anomalous super-diffusion). Although our approach does not provide insight into the details of the specific acceleration mechanisms involved, its benefits are that it relates acceleration to the energy release, and it well describes the stochastic nature of the acceleration process.Comment: 37 pages, 10 figures, one of them in color; in press at ApJ (2004

Impulsive electron acceleration by Gravitational Waves

We investigate the non-linear interaction of a strong Gravitational Wave with the plasma during the collapse of a massive magnetized star to form a black hole, or during the merging of neutron star binaries (central engine). We found that under certain conditions this coupling may result in an efficient energy space diffusion of particles. We suggest that the atmosphere created around the central engine is filled with 3-D magnetic neutral sheets (magnetic nulls). We demonstrate that the passage of strong pulses of Gravitational Waves through the magnetic neutral sheets accelerates electrons to very high energies. Superposition of many such short lived accelerators, embedded inside a turbulent plasma, may be the source for the observed impulsive short lived bursts. We conclude that in several astrophysical events, gravitational pulses may accelerate the tail of the ambient plasma to very high energies and become the driver for many types of astrophysical bursts.Comment: 13 pages, 8 figures, accepted to The Astrophysical Journa

Quantum Algorithm Implementations for Beginners

As quantum computers become available to the general public, the need has arisen to train a cohort of quantum programmers, many of whom have been developing classical computer programs for most of their careers. While currently available quantum computers have less than 100 qubits, quantum computing hardware is widely expected to grow in terms of qubit count, quality, and connectivity. This review aims to explain the principles of quantum programming, which are quite different from classical programming, with straightforward algebra that makes understanding of the underlying fascinating quantum mechanical principles optional. We give an introduction to quantum computing algorithms and their implementation on real quantum hardware. We survey 20 different quantum algorithms, attempting to describe each in a succinct and self-contained fashion. We show how these algorithms can be implemented on IBM's quantum computer, and in each case, we discuss the results of the implementation with respect to differences between the simulator and the actual hardware runs. This article introduces computer scientists, physicists, and engineers to quantum algorithms and provides a blueprint for their implementations

Plasmoids in Reconnecting Current Sheets: Solar and Terrestrial Contexts Compared

Magnetic reconnection plays a crucial role in violent energy conversion occurring in the environments of high electrical conductivity, such as the solar atmosphere, magnetosphere, and fusion devices. We focus on the morphological features of the process in two different environments, the solar atmosphere and the geomagnetic tail. In addition to indirect evidence that indicates reconnection in progress or having just taken place, such as auroral manifestations in the magnetosphere and the flare loop system in the solar atmosphere, more direct evidence of reconnection in the solar and terrestrial environments is being collected. Such evidence includes the reconnection inflow near the reconnecting current sheet, and the outflow along the sheet characterized by a sequence of plasmoids. Both turbulent and unsteady Petschek-type reconnection processes could account for the observations. We also discuss other relevant observational consequences of both mechanisms in these two settings. While on face value, these are two completely different physical environments, there emerge many commonalities, for example, an Alfven speed of the same order of magnitude, a key parameter determining the reconnection rate. This comparative study is meant as a contribution to current efforts aimed at isolating similarities in processes occurring in very different contexts in the heliosphere, and even in the universe.Comment: 21 pages, 9 figures, in press at J. Geophys. Res. (Space Physics), for the special NESSC section on Comparative Aspects of Magnetic Reconnectio