Mirroring within the Fokker-Planck formulation of cosmic ray pitch angle scattering in homogeneous magnetic turbulence

The Fokker-Planck coefficient for pitch angle scattering, appropriate for cosmic rays in homogeneous, stationary, magnetic turbulence, is computed from first principles. No assumptions are made concerning any special statistical symmetries the random field may have. This result can be used to compute the parallel diffusion coefficient for high energy cosmic rays moving in strong turbulence, or low energy cosmic rays moving in weak turbulence. Becuase of the generality of the magnetic turbulence which is allowed in this calculation, special interplanetary magnetic field features such as discontinuities, or particular wave modes, can be included rigorously. The reduction of this results to previously available expressions for the pitch angle scattering coefficient in random field models with special symmetries is discussed. The general existance of a Dirac delta function in the pitch angle scattering coefficient is demonstrated. It is proved that this delta function is the Fokker-Planck prediction for pitch angle scattering due to mirroring in the magnetic field

Test particle propagation in magnetostatic turbulence. 1. Failure of the diffusion approximation

The equation which governs the quasi-linear approximation to the ensemble and gyro-phase averaged one-body probability distribution function is constructed from first principles. This derived equation is subjected to a thorough investigation in order to calculate the possible limitations of the quasi-linear approximation. It is shown that the reduction of this equation to a standard diffusion equation in the Markovian limit can be accomplished through the application of the adiabatic approximation. A numerical solution of the standard diffusion equation in the Markovian limit is obtained for the narrow parallel beam injection. Comparison of the diabatic and adiabatic results explicitly demonstrates the failure of the Markovian description of the probability distribution function. Through the use of a linear time-scale extension the failure of the adiabatic approximation, which leads to the Markovian limit, is shown to be due to mixing of the relaxation and interaction time scales in the presence of the strong mean field

Test particle propagation in magnetostatic turbulence. 3: The approach to equilibrium

The asymptotic behavior, for large time, of the quasi-linear diabatic solutions and their local approximations is considered. A time averaging procedure is introduced which yields the averages of these solutions over time intervals which contain only large time values. A discussion of the quasi-linear diabatic solutions which is limited to those solutions that are bounded from below as functions of time is given. It is shown that as the upper limit of the time averaging interval is allowed to approach infinity the time averaged quasi-linear diabatic solutions must approach isotropy (mu-independence). The first derivative with respect to mu of these solutions is also considered. This discussion is limited to first derivatives which are bounded functions of time. It is shown that as the upper limit of the time averaging interval is allowed to approach infinity, the time averaged first derivative must approach zero everywhere in mu except at mu = 0 where it must approach a large value which is calculated. The impact of this large derivative on the quasi-linear expansion scheme is discussed. An H-theorem for the first local approximation to the quasi-linear diabatic solutions is constructed. Without time averaging, the H-theorem is used to determine sufficient conditions for the first local approximate solutions to asymptote, with increasing time, to exactly the same final state which the time averaged quasi-linear diabatic solutions must approach as discussed above

Test particle propagation in magnetostatic turbulence. 2: The local approximation method

An approximation method for statistical mechanics is presented and applied to a class of problems which contains a test particle propagation problem. All of the available basic equations used in statistical mechanics are cast in the form of a single equation which is integrodifferential in time and which is then used as the starting point for the construction of the local approximation method. Simplification of the integrodifferential equation is achieved through approximation to the Laplace transform of its kernel. The approximation is valid near the origin in the Laplace space and is based on the assumption of small Laplace variable. No other small parameter is necessary for the construction of this approximation method. The n'th level of approximation is constructed formally, and the first five levels of approximation are calculated explicitly. It is shown that each level of approximation is governed by an inhomogeneous partial differential equation in time with time independent operator coefficients. The order in time of these partial differential equations is found to increase as n does. At n = 0 the most local first order partial differential equation which governs the Markovian limit is regained

Methodology for Measuring Distortions to Agricultural Incentives

distorted incentives, agricultural and trade policy reforms, national agricultural development, Agricultural and Food Policy, International Relations/Trade, F13, F14, Q17, Q18,

The Fokker-Planck coefficient for pitch-angle scattering of cosmic rays

For the case of homogeneous, isotropic magnetic field fluctuations, it is shown that most theories which are based on the quasi-linear and adiabatic approximation yield the same integral for the Fokker-Planck coefficient for the pitch angle scattering of cosmic rays. For example, despite apparent differences, the theories due to Jokipii and to Klimas and Sandri yield the same integral. It is also shown, however, that this integral in most cases has been evaluated incorrectly in the past. For large pitch angles these errors become significant, and for pitch angles of 90 deg the actual Fokker-Planck coefficient contains a delta function. The implications for these corrections relating cosmic ray diffusion coefficients to observed properties of the interplanetary magnetic field are discussed

Vorticity Dynamics and Sound Generation in Two-Dimensional Fluid Flow

An approximate solution to the two-dimensional incompressible fluid equations is constructed by expanding the vorticity field in a series of derivatives of a Gaussian vortex. The expansion is used to analyze the motion of a corotating Gaussian vortex pair, and the spatial rotation frequency of the vortex pair is derived directly from the fluid vorticity equation. The resulting rotation frequency includes the effects of finite vortex core size and viscosity and reduces, in the appropriate limit, to the rotation frequency of the Kirchhoff point vortex theory. The expansion is then used in the low Mach number Lighthill equation to derive the far-field acoustic pressure generated by the Gaussian vortex pair. This pressure amplitude is compared with that of a previous fully numerical simulation in which the Reynolds number is large and the vortex core size is significant compared to the vortex separation. The present analytic result for the far-field acoustic pressure is shown to be substantially more accurate than previous theoretical predictions. The given example suggests that the vorticity expansion is a useful tool for the prediction of sound generated by a general distributed vorticity field

Linear Ramps of Interaction in the Fermionic Hubbard Model

We study the out of equilibrium dynamics of the Fermionic Hubbard Model induced by a linear ramp of the repulsive interaction $U$ from the metallic state through the Mott transition. To this extent we use a time dependent Gutzwiller variational method and complement this analysis with the inclusion of quantum fluctuations at the leading order, in the framework of a $Z_2$ slave spin theory. We discuss the dynamics during the ramp and the issue of adiabaticity through the scaling of the excitation energy with the ramp duration $\tau$. In addition, we study the dynamics for times scales longer than the ramp time, when the system is again isolated and the total energy conserved. We establish the existence of a dynamical phase transition analogous to the one present in the sudden quench case and discuss its properties as a function of final interaction and ramp duration. Finally we discuss the role of quantum fluctuations on the mean field dynamics for both long ramps, where spin wave theory is sufficient, and for very short ramps, where a self consistent treatment of quantum fluctuations is required in order to obtain relaxation.Comment: v2: 19 pages, 14 figures, published versio

Probabilistic volcanic hazard assessment and eruption forecasting: the Bayesian Event Tree approach

The purpose of this report is to discuss in detail the importance and prerogatives of quantitative volcanic hazard assessment and to describe the main features of a Bayesian model designed to achieve this goal. Ideas, models and results come out from the work made in the framework of the INGV-DPC V4 project, and partially from the application of the technique to Campi Flegrei (V3_2) and Vesuvius (V3_4). Here, we examine in depth the practical and philosophical implications of the approach, and report only a brief summary of the technical details that can be found on the cited literature

BET VH: a probabilistic tool for long-term volcanic hazard assessment

In this paper we illustrate a Bayesian Event Tree to estimate Volcanic Hazard (BET VH). The procedure enables us to calculate the probability of any kind of long-term hazardous event for which we are interested, accounting for the intrinsic stochastic nature of volcanic eruptions and our limited knowledge regarding related processes. For the input, the code incorporates results from numerical models simulating the impact of hazardous volcanic phenomena on an area, and data from the eruptive history. For the output, the code provides a wide and exhaustive set of spatio-temporal probabilities of different events; these probabilities are estimated by means of a Bayesian approach that allows all uncertainties to be properly accounted for. The code is able to deal with many eruptive settings simultaneously, weighting each with its own probability of occurrence. In a companion paper, we give a detailed example of application of this tool to the Campi Flegrei caldera, in order to estimate the hazard from tephra fall