Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “charge density waves”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42

Possible role of grain-boundary and dislocation structure for the magnetic-flux trapping behavior of niobium: A first-principles study

In this work, first-principles methods were used to understand magnetic flux trapping at vacancies, dislocations, and grain boundaries in high-purity superconducting niobium. Full-potential linear augmented plane-wave methods were applied in progressively greater complexity, starting at simple vacancies and extending to screw dislocations and tilt grain boundaries to analyze the effects of magnetic field on the superconducting state surrounding these defects. Density-functional theory calculations identified changes in electronic structure at the dislocation core and different types of symmetric tilt grain boundaries relative to bulk niobium. Electron redistribution enhanced nonparamagnetic effects and thus perturb superconductivity, resulting in local conditions suitable for premature flux penetration and subsequently flux pinning. Since the coherence length of superconducting niobium at 0 K is significantly larger than the lattice parameter, the effects of line and planar defects in niobium are predicted to be stronger for defect clusters than single defects in isolation, which is consistent with recent experimental observations. Controlling accumulation or depletion of charge at the defects, e.g., by segregation of an impurity atom, can mitigate these tendencies thus increasing the quality of superconducting niobium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cascades between Light and Heavy Fermions in the Normal State of Magic-Angle Twisted Bilayer Graphene

We present a framework for understanding the cascade transitions and the Landau level degeneracies of twisted bilayer graphene. The Coulomb interaction projected onto narrow bands causes the charged excitations at an integer filling to disperse, forming new bands. If the excitation moves the filling away from the charge neutrality point, then it has a band minimum at the moiré Brillouin zone center with a small mass that compares well with the experiment; if towards the charge neutrality point, then it has a much larger mass and a higher degeneracy. At a nonzero density away from an integer filling the excitations interact. The system on the small mass side has a large bandwidth and forms a Fermi liquid. On the large mass side the bandwidth is narrow, the compressibility is negative and the Fermi liquid is likely unstable. This explains the observed sawtooth features in compressibility, the Landau fans pointing away from charge neutrality and their degeneracies. The framework sets the stage for superconductivity at lower temperatures.

36 MATERIALS SCIENCE↗

Interaction of ring current and radiation belt protons with ducted plasmaspheric hiss. 1: Diffusion coefficients and timescales

Protons that are convected into the inner magnetosphere in response to enhanced magnetic activity can resonate with ducted plasmaspheric hiss in the outer plasmasphere via an anomalous Doppler-shifted cyclotron resonance. Plasmaspheric hiss is a right-hand-polarized electromagnetic emission that is observed to fill the plasmasphere on a routine basis. When plasmaspheric hiss is confined within field-aligned ducts or guided along density gradients, wave normal angles remain largely below 45 deg. This allows resonant interactions with ions at typical ring current and radiation belt energies to take place. Such field-aligned ducts have been observed both within the plasmasphere and in regions outside of the plasmasphere. Wave intensities are estimated using statistical information from studies of detached plasma regions. Diffusion coefficients are presented for a range of L shells and proton energies for a fixed wave distribution. Harmonic resonances in the range N = +/-100 are considered in order to include interactions between hiss at 100 Hz to 2 kHz frequencies, and protons in the energy range between approximately 10 keV and 1000 keV. Diffusion timescales are estimated to be of the order of tens of days and comparable to or shorter than lifetimes for Coulomb decay and charge exchange losses over most of the energy and spatial ranges of interest.

Kozyra, J. U.↗

Voyager observations of ion phase space densities in the Jovian magnetosphere

Data from the Voyager low-energy charged particle experiment (LECP) were used to calculate ion phase space densities in Jupiter's magnetosphere. The calculation of F at constant mu and J(2) requires the determination of particle fluxes at specific pitch angles and energies. It is shown that the greatest uncertainties in the determination of F from LECP data arise from the fits to the measured pitch angle distributions and differential energy spectra. An estimate is provided of this uncertainty and of others arising from model radial diffusion coefficients and magnetic fields. The general nature of the curves is consistent with inward diffusion of these energetic ions from the outer magnetosphere combined with losses near but starting beyond the orbit of Io. These losses are not consistent with simple satellite sweeping by Io alone and are probably due to an Io-torus wave-particle interaction. The lifetime against loss deduced from the data is approximately 20,000 s near Io and is a value consistent with strong diffusion losses. The measured ion loss rate declines much more rapidly than the strong diffusion loss rate as L increases from 7 to 9, suggesting that the loss rate is well below the strong-diffusion rate beyond 7 Jupiter radii for ions observed by the LECP.

Paonessa, M.↗

Holographic thermal correlators: a tale of Fuchsian ODEs and integration contours

We analyze real-time thermal correlation functions of conserved currents in holographic field theories using the grSK geometry, which provides a contour prescription for their evaluation. We demonstrate its efficacy, arguing that there are situations involving components of conserved currents, or derivative interactions, where such a prescription is, in fact, essential. To this end, we first undertake a careful analysis of the linearized wave equations in AdS black hole backgrounds and identify the branch points of the solutions as a function of (complexified) frequency and momentum. All the equations we study are Fuchsian with only regular singular points that for the most part are associated with the geometric features of the background. Special features, e.g., the appearance of apparent singular points at the horizon, whence outgoing solutions end up being analytic, arise at higher codimension loci in parameter space. Using the grSK geometry, we demonstrate that these apparent singularities do not correspond to any interesting physical features in higher-point functions. We also argue that the Schwinger-Keldysh collapse and KMS conditions, implemented by the grSK geometry, continue to hold even in the presence of such singularities. For charged black holes above a critical charge, we furthermore demonstrate that the energy density operator does not possess an exponentially growing mode, associated with ‘pole-skipping’, from one such apparent singularity. Our analysis suggests that the connection between the scrambling physics of black holes and energy transport has, at best, a limited domain of validity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energetic particles at Uranus

The energetic particle measurements by the low-energy charged-particle and cosmic-ray instruments on the Voyager 2 spacecraft in the magnetosphere of Uranus are reviewed. Upstream events were observed outside the Uranian bow shock, probably produced by ion escape from the magnetosphere. Evidence of earthlike substorm activity was discovered within the Uranian magnetosphere. A proton injection event was observed within the orbit of Umbriel and proton events were observed in the magnetotail plasma-sheet boundary layer that are diagnostic of earthlike substorms. The magnetospheric composition is totally dominated by protons, with only a trace abundance of H(2+) and no evidence for He or heavy ions; the Uranian atmophere is argued to be the principal plasma source. Phase-space densities of medium energy protons show inward radial diffusion and are quantitatively similar to those observed at the earth, Jupiter, and Saturn. These findings and plasma wave data suggest the existence of structures analogous to the earth's plasmasphere and plasmapause.

Cheng, Andrew F.↗

Facile ab initio approach for self-localized polarons from canonical transformations

Electronic states in a crystal can localize due to strong electron-phonon (e-ph) interactions, forming so-called small polarons. Methods to predict the formation and energetics of small polarons are either computationally costly or not geared toward quantitative predictions. Here we show a formalism based on canonical transformations to compute the polaron formation energy and wave function using ab initio e-ph interactions. Comparison of the calculated polaron and band-edge energies allows us to determine whether charge carriers in a material favor a localized small polaron over a delocalized Bloch state. Due to its low computational cost, our approach enables efficient studies of the formation and energetics of small polarons, as we demonstrate by investigating electron and hole polaron formation in alkali halides and metal oxides and peroxides. We outline refinements of our scheme and extensions to compute transport in the polaron hopping regime.

36 MATERIALS SCIENCE↗

A theoretical model for sheath dynamics in the two-plasma mode of a cylindrical filament discharge

We present a refined understanding of cylindrical plasma-facing emitter cathodes by examining the limitations of the conventional space charge limited model, which overlooks ion trapping dynamics within the virtual cathode. While conventional theory distinguishes between temperature-limited and space charge limited regimes based on electron current constraints, our findings reveal that trapped ions can form a quasi-neutral “second plasma,” expanding upstream and interacting with the primary plasma. In cylindrical and other non-planar geometries, an expanding trapped-ions layer produces current enhancement by converting more of the emitted electron flux to passing flux. The current enhancement leads to a complex coupling between the trapped-ions plasma and upstream plasma through ionizing collisions in the upstream region and charge exchange collisions in the virtual cathode. The coupling results in mutual density enhancement and at the same time, antagonistic expansion dynamics between the two plasma layers. The process is understood using a new “aid-compete” model. Numerical solutions, validated against simulations, reveal the model's capacity to predict system evolution and complex discharge behaviors, albeit with qualitative approximations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of hole polaron formation on excitonic transitions in MgO from first principles

Here, we present a first-principles investigation of the excitonic properties of magnesia (MgO), an ionic insulator known to host hole polarons. We combine a density functional theory-based approach for structural relaxation in the presence of the hole and many-body perturbation theory to describe the excitonic properties. We determine that the hole polaron introduces new in-gap occupied states 0.6–0.8 eV above the valence band maximum that lead to two low-energy peaks in the optical spectrum. The predicted redshift of the lowest-energy transition due to polaron formation of 0.8 eV agrees well with the experimental Stokes shift of 0.8–0.9 eV. Analysis of the exciton wave function indicates that the electron-hole pair consists of a localized hole and delocalized electron, but that the wave function retains its Wannier-Mott character even in the presence of the hole polaron. Our study demonstrates that combining these previously established methods allows for a relatively computationally inexpensive approach to studying the exciton polaron in materials where only one charge carrier forms a polaron.

electronic structure↗

Valley magnetism, nematicity, and density wave orders in twisted bilayer graphene

We analyze density wave and Pomeranchuk orders in twisted bilayer graphene. This complements our earlier analysis of the pairing instabilities. We assume that near half filling of either conduction or valence band, the Fermi level is close to Van Hove points, where the density of states diverges, and study potential instabilities in the particle-hole channel within a patch model with two valley degrees of freedom. The hexagonal symmetry of twisted bilayer graphene allows for either six or twelve Van Hove points. We consider both cases and find the same two leading candidates for particle-hole order. One is an SU(2)-breaking spin state with ferromagnetism within a valley. A subleading intervalley hopping induces antiferromagnetism between the valleys. The same state has also been obtained in strong-coupling approaches, indicating that this order is robust. Here, the other is a mixed state with 120° complex spin order and orthogonal complex charge order. In addition, we find a weaker but still attractive interaction in nematic channels, and discuss the type of a nematic order.

36 MATERIALS SCIENCE↗

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗

Characterizing IHE Response to Multiple Shock Loading

The response of high explosives to shock loading is traditionally measured with a steady loading pressure. In many accident scenarios involving fragment impact, however, a loading duration that is shorter than the build up to detonation may occur. Fragments passing through multiple materials before reaching a high explosive charge may produce loading that is comprised of more than one shock wave. Additionally, the build up to detonation in high explosive corner turning loads the explosive a short duration pressure pulse, since rarefactions can often rapidly overtake the reactive wave. For these reasons, we have studied the response of the insensitive high explosive (IHE) materials PBX 9502 and LX-17 to complex loadings of varied intensity and duration. We refer to a single loading of limited duration as a “thin pulse”, whereas more complex scenarios were studied with an impactor that produces a double shock in the explosive. The following report presents experimental data and analyses of thin pulse shock initiation and double shock experiments designed to guide development of models of Insensitive High Explosives (IHEs) under controlled one-dimensional conditions relevant to accident scenarios and corner turning. Thin pulse shock initiation data on PBX 9502 and LX-17 were obtained under varied pulse duration, pressed density, and temperature conditions in order to probe various parameters essential for the development of a physics-based Cheetah reactive flow hotspot model. In situ pressure gauges provide insight into the degree of reaction in the explosive that are not obtainable with optical PDV measurements or distance measurements such as run to detonation. Double shock data was obtained to inform a Composition Aware Cheetah model which can be applied to any TATB-based IHEs. This model supports efforts to find a new IHE formulation and potentially incorporate new binders into IHE formulations. Simulations of each experiment are included to demonstrate the utility of these focused experiments to developing models of HE behavior. One-dimensional gas gun experiments are essential for characterizing shocked HE behavior and informing HE models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pb-Apatite Framework as a Generator of Novel Flat-Band CuO-Based Physics

Based upon density functional theory (DFT) calculations, we present the basic electronic structure of CuPb 9 (PO 4 ) 6 O (Cu-doped lead apatite, aka LK-99), in two scenarios: (1) where the structure is constrained to the P3 symmetry and (2) where no symmetry is imposed. At the DFT level, the former is predicted to be metallic while the latter is found to be a charge-transfer insulator. In both cases the filling of these states is nominally d 9 , consistent with the standard Cu 2+ valence state, and Cu with a local magnetic moment of order 0.7 μ B . In the metallic case we find these states to be unusually flat (∼0.2 eV dispersion), giving a very high density of electronic states (DOS) at the Fermi level that we argue can be a host for novel electronic physics. The flatness of the bands is the likely origin of symmetry-lowering gapping possibilities that would remove the spectral weight from E F . Motivated by some initial experimental observations of metallic or semiconducting behavior, we propose that disorder (likely structural) is responsible for closing the gap. Here, we consider a variety of possibilities that could possibly close the charge-transfer gap but limit consideration to kinds of disorder that preserve electron count. Of the possible kinds we considered (spin disorder, O populating vacancy sites, and Cu on less energetically favorable Pb sites), the local Cu moment, and consequently the charge-transfer gap, remains robust. We conclude that disorder responsible for metallic behavior entails some kind of doping where the electron count changes. Further, we claim that the emergence of the flat bands should be due to weak wave function overlap between the orbitals on Cu and O sites, owing to the directional character of the constituent orbitals. Therefore, finding an appropriate host structure for minimizing hybridization between Cu and O while allowing them to still weakly interact should be a promising route for generating flat bands at E F which can lead to interesting electronic phenomena, regardless of whether LK-99 is a superconductor.

36 MATERIALS SCIENCE↗

Particle-in-cell simulations of the critical ionization velocity effect in finite size clouds

The critical ionization velocity (CIV) mechanism in a finite size cloud is studied with a series of electrostatic particle-in-cell simulations. It is observed that an initial seed ionization, produced by non-CIV mechanisms, generates a cross-field ion beam which excites a modified beam-plasma instability (MBPI) with frequency in the range of the lower hybrid frequency. The excited waves accelerate electrons along the magnetic field up to the ion drift energy that exceeds the ionization energy of the neutral atoms. The heated electrons in turn enhance the ion beam by electron-neutral impact ionization, which establishes a positive feedback loop in maintaining the CIV process. It is also found that the efficiency of the CIV mechanism depends on the finite size of the gas cloud in the following ways: (1) Along the ambient magnetic field the finite size of the cloud, L (sub parallel), restricts the growth of the fastest growing mode, with a wavelength lambda (sub m parallel), of the MBPI. The parallel electron heating at wave saturation scales approximately as (L (sub parallel)/lambda (sub m parallel)) (exp 1/2); (2) Momentum coupling between the cloud and the ambient plasma via the Alfven waves occurs as a result of the finite size of the cloud in the direction perpendicular to both the ambient magnetic field and the neutral drift. This reduces exponentially with time the relative drift between the ambient plasma and the neutrals. The timescale is inversely proportional to the Alfven velocity. (3) The transvers e charge separation field across the cloud was found to result in the modulation of the beam velocity which reduces the parallel heating of electrons and increases the transverse acceleration of electrons. (4) Some energetic electrons are lost from the cloud along the magnetic field at a rate characterized by the acoustic velocity, instead of the electron thermal velocity. The loss of energetic electrons from the cloud seems to be larger in the direction of plasma drift relative to the neutrals, where the loss rate is characterized by the neutral drift velocity. It is also shown that a factor of 4 increase in the ambient plasma density, increases the CIV ionization yield by almost 2 orders of magnitude at the end of a typical run. It is concluded that a larger ambient plasma density can result in a larger CIV yield because of (1) larger seed ion production by non-CIV mechanisms, (2) smaller Alfven velocity and hence weak momentum coupling, and (3) smaller ratio of the ion beam density to the ambient ion density, and therefore a weaker modulation of the beam velocity. The simulation results are used to interpret various chemical release experiments in space.

Moghaddam-Taaheri, E.↗

Femtosecond control of phonon dynamics near a magnetic order critical point

The spin-phonon interaction in spin density wave (SDW) systems often determines the free energy landscape that drives the evolution of the system. When a passing energy flux, such as photoexcitation, drives a crystalline system far from equilibrium, the resulting lattice displacement generates transient vibrational states. Manipulating intermediate vibrational states in the vicinity of the critical point, where the SDW order parameter changes dramatically, would then allow dynamical control over functional properties. Here we combine double photoexcitation with an X-ray free-electron laser (XFEL) probe to control and detect the lifetime and magnitude of the intermediate vibrational state near the critical point of the SDW in chromium. We apply Landau theory to identify the mechanism of control as a repeated partial quench and sub picosecond recovery of the SDW. Our results showcase the capabilities to influence and monitor quantum states by combining multiple optical photoexcitations with an XFEL probe. They open new avenues for manipulating and researching the behaviour of photoexcited states in charge and spin order systems near the critical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Velocity Structure and Plasma Properties in Halo CMEs

We have identified a set of 23 Halo CMEs through July 2002 and 21 Partial Halo CMEs from the LASCO Halo CME Mail Archive for which Ultraviolet Coronagraph Spectrometer (UVCS) spectra exist. For each event we have collected basic information such as the event speed, whether or not UVCS caught the bright front, lines detected, Doppler shift and associated flare class. We have also obtained excellent observations of some of the spectacular events in November 2003, and we have made theoretical calculations pertaining to CME expansion at the heights observed by UVCS. We first analyzed the halo CMEs on 21 April and 24 August 2002 and the partial halo on 23 July 2002, because the X-class flares associated with these CMEs were extensively observed by RHESSI and other instruments as part of the MAX MILLENIUM campaign. These very fast CMEs showed extremely violent disruption of the pre-CME streamers, little or no cool prominence material, and the unusual (for UVCS heights) hot emission line [Fe XVIII]. Results, including a discussion of the current sheet interpretation for the [Fe XVIII] emission, are published in Raymond et al. and presented at the Fall 2002 AGU meeting and the solar physics summer school in L'Aquila, Italy. We are currently preparing two papers on the Dec. 28, 2000 partial halo event. This event was chosen to take advantage of the SEP event measured by WIND and ACE, and because a Type II radio burst coincides with the time that broad, blue-shifted O VI emission appeared in the UVCS spectra. One paper deals with a new density and velocity diagnostic for very fast CMEs; pumping of O VI lambda 1032 by Ly beta and pumping of O VI lambda 1038 by O VI lambda 1032. The other discusses physics of the shock wave and association with the SEP event. In the coming year we plan to expand the list of Halo and Partial Halo events observed by UVCS through the end of 2003. We will look at those events as a class to search for correlation between UV spectral characteristics and other CME and flare parameters. We will also choose several events for more detailed study, probably including the November 2003 events. We expect to support extended visits to CfA by S. Mancuso and A. Ciaravella. In the past year the grant covered some salary support for members of the SAO UVCS team and a 1 month visit to CfA by Angela Ciaravella, along with trips to meetings by J. Lin and J. Raymond and page charges for two papers.

Wagner, William↗

Optical Imaging and Radiometric Modeling and Simulation

OPTOOL software is a general-purpose optical systems analysis tool that was developed to offer a solution to problems associated with computational programs written for the James Webb Space Telescope optical system. It integrates existing routines into coherent processes, and provides a structure with reusable capabilities that allow additional processes to be quickly developed and integrated. It has an extensive graphical user interface, which makes the tool more intuitive and friendly. OPTOOL is implemented using MATLAB with a Fourier optics-based approach for point spread function (PSF) calculations. It features parametric and Monte Carlo simulation capabilities, and uses a direct integration calculation to permit high spatial sampling of the PSF. Exit pupil optical path difference (OPD) maps can be generated using combinations of Zernike polynomials or shaped power spectral densities. The graphical user interface allows rapid creation of arbitrary pupil geometries, and entry of all other modeling parameters to support basic imaging and radiometric analyses. OPTOOL provides the capability to generate wavefront-error (WFE) maps for arbitrary grid sizes. These maps are 2D arrays containing digital sampled versions of functions ranging from Zernike polynomials to combination of sinusoidal wave functions in 2D, to functions generated from a spatial frequency power spectral distribution (PSD). It also can generate optical transfer functions (OTFs), which are incorporated into the PSF calculation. The user can specify radiometrics for the target and sky background, and key performance parameters for the instrument s focal plane array (FPA). This radiometric and detector model setup is fairly extensive, and includes parameters such as zodiacal background, thermal emission noise, read noise, and dark current. The setup also includes target spectral energy distribution as a function of wavelength for polychromatic sources, detector pixel size, and the FPA s charge diffusion modulation transfer function (MTF).

Ha, Kong Q.↗

Particle Acceleration in Active Galactic Nuclei

The high efficiency of energy generation inferred from radio observations of quasars and X-ray observations of Seyfert active galactic nuclei (AGNs) is apparently achieved only by the gravitational conversion of the rest mass energy of accreting matter onto supermassive black holes. Evidence for the acceleration of particles to high energies by a central engine is also inferred from observations of apparent superluminal motion in flat spectrum, core-dominated radio sources. This phenomenon is widely attributed to the ejection of relativistic bulk plasma from the nuclei of active galaxies, and accounts for the existence of large scale radio jets and lobes at large distances from the central regions of radio galaxies. Reports of radio jets and superluminal motion from galactic black hole candidate X-ray sources indicate that similar processes are operating in these sources. Observations of luminous, rapidly variable high-energy radiation from active galactic nuclei (AGNs) with the Compton Gamma Ray Observatory show directly that particles are accelerated to high energies in a compact environment. The mechanisms which transform the gravitational potential energy of the infalling matter into nonthermal particle energy in galactic black hole candidates and AGNs are not conclusively identified, although several have been proposed. These include direct acceleration by static electric fields (resulting from, for example, magnetic reconnection), shock acceleration, and energy extraction from the rotational energy of Kerr black holes. The dominant acceleration mechanism(s) operating in the black hole environment can only be determined, of course, by a comparison of model predictions with observations. The purpose of the work proposed for this grant was to investigate stochastic particle acceleration through resonant interactions with plasma waves that populate the magnetosphere surrounding an accreting black hole. Stochastic acceleration has been successfully applied to the problem of ion and electron energization in solar flares, and is capable of accounting for a wide range of both neutral and charged particle emissions. It is also a component in diffusive shock acceleration, since pitch-angle scattering (which is necessary for multiple shock crossings) is accompanied by diffusion in momentum space, which in turn yields a net systematic energy gain; however, stochastic energization will dominate the first-order shock process only in certain parameter regimes. Although stochastic acceleration has been applied to particle energization in the lobes of radio galaxies, its application to the central regions of AGNs has only recently been considered, but not in detail. We proposed to systematically investigate the plasma processes responsible for stochastic particle acceleration in black hole magnetospheres along with the energy-loss processes which impede particle energization. To this end we calculated acceleration rates and escape time scales for protons and electrons resonating with Alfven waves, and for electrons resonating with whistlers. Assuming either a Kolmogorov or Kraichnan wave spectrum, accretion at the Eddington limit, magnetic field strengths near equipartition, and turbulence energy densities approx. 10% of the total magnetic field energy density, we find that Alfven waves accelerate protons to Lorentz factors approx, equals 10(exp 4) - 10(exp 6) before they escape from the system. Acceleration of electrons by fast mode and whistler waves can produce a nonthermal population of relativistic electrons whose maximum energy is determined by a competition with radiation losses.

Miller, James A.↗