LOCAL TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF THE COLLISIONLESS MAGNETOROTATIONAL INSTABILITY

The magnetorotational instability (MRI) is a crucial mechanism of angular momentum transport in a variety of astrophysical accretion disks. In systems accreting at well below the Eddington rate, such as the central black hole in the Milky Way (Sgr A*), the plasma in the disk is essentially collisionless. We present a nonlinear study of the collisionless MRI using first-principles particle-in-cell plasma simulations. We focus on local two-dimensional (axisymmetric) simulations, deferring more realistic three-dimensional simulations to future work. For simulations with net vertical magnetic flux, the MRI continuously amplifies the magnetic field, B, until the Alfven velocity, v(A), is comparable to the speed of light, c (independent of the initial value of v(A)/c). This is consistent with the lack of saturation of MRI channel modes in analogous axisymmetric MHD simulations. The amplification of the magnetic field by the MRI generates a significant pressure anisotropy in the plasma (with the pressure perpendicular to B being larger than the parallel pressure). We find that this pressure anisotropy in turn excites mirror modes and that the volume-averaged pressure anisotropy remains near the threshold for mirror mode excitation. Particle energization is due to both reconnection and viscous heating associated with the pressure anisotropy. Reconnection produces a distinctive power-law component in the energy distribution function of the particles, indicating the likelihood of non-thermal ion and electron acceleration in collisionless accretion disks. This has important implications for interpreting the observed emission-from the radio to the gamma-rays-of systems such as Sgr A*.

Laser wakefield and self-modulation of driving pulse

A self-consistent theory of plasma response to a single laser beam is proposed. The driving pump is not viewed as invariant during its interaction with the plasmas. Its modulation by the plasmas has an obvious influence on the strength of the wakefield behind the pulse. This suggests that the compression of the low-intensity pulse by the plasmas might be a possible way to excite largae-amplitude wakefield. (C) 2003 American Institute of Physics.

紧聚焦的超短超强激光脉冲在真空中加速斜入射的相对论电子

研究了紧聚焦的线偏振飞秒超强高斯激光脉冲俘获并剧烈加速斜入射低能相对论电子的效应，发现被俘获的电子在激光脉冲纵向有质动力的强大加速作用下，可以获得GeV量级的能量，并详细研究了入射电子的初能量、斜入射角、电子与激光脉冲的相对延迟时间和激光脉冲宽度等条件对电子能量增益的影响，发现当激光脉宽超过10λ时，脉宽对电子能量增益影响不大．

相对论飞秒激光脉冲在真空中对预加速电子的加速

通过求解电子运动的相对论方程，发现预加速电子在超强超短激光脉冲的作用下可以获得很高的能量增益．飞秒激光脉冲的上升沿在焦点附近的区域有效加速电子后，电子和光脉冲一起传播一段距离（远大于瑞利长度）后，激光强度变得很弱，从而使脉冲下降沿对电子的减速作用可以忽略不计，因此电子只经历加速过程而没有被减速，当电子和光脉冲分离时，电子获得了很高的能量增益．当光强为10^19W／cm^2，电子的初始能量为MeV量级时，电子的能量增益可以达到0．1GeV．进一步讨论了电子的能量增益与电子的初始条件与激光脉冲的参数之间的关系

激光脉冲宽度对有质动力加速电子的影响

基于真空中单电子运动模型，编程计算得到了高斯激光脉冲与初始位于激光传播轴上电子的相互作用结果。不同激光参鼍条件下，得到了电子的能量增益与激光强度、焦斑大小和脉冲宽度关系。结果表明，高斯激光脉冲焦斑较大时，电子没有明显的能量增益，高斯激光脉冲焦斑太小时，电子也没有明显的能量增益。电子的能量增益有一个最佳焦斑大小。在相同激光强度下，电子能量增益的最佳焦斑大小随脉冲宽度的增大而增大，但最佳焦斑大小与脉冲宽度的比值基本上是不变的。

Theory of terahertz-photocurrent resonances in miniband superlattices

The currents of de and ac components and their phase-angle cosines for a superlattice under a direct bias and alternating field are calculated with the balance equations. It is found that the de currents as functions of the direct field show resonance peaks at the fields corresponding to the Bloch frequency equal to n omega. With increasing alternating field intensity the resonance peaks of higher harmonic increase, and simultaneously the first peak caused by the de field decreases. The results are in good agreement with the experimental results, indicating that this resonance can be understood in terms of electron acceleration within the miniband, i.e., it is a bulk superlattice effect, rather than caused by the electric-field localization mechanism (Wannier Stark ladder). The phase-angle cosine for the first harmonic cos phi(1) becomes negative when the Bloch frequency increases to be larger than the frequency of the ac field omega, and it also shows resonance peaks at the resonance frequencies n omega. The negative cos phi(1) may cause the energy transferred to the alternating field, i.e., oscillation of the system.

Self-consistent wave-particle interactions in dispersive scale long-period field-line-resonances

Using 1D Vlasov drift-kinetic computer simulations, it is shown that electron trapping in long period standing shear Alfven waves (SAWs) provides an efficient energy sink for wave energy that is much more effective than Landau damping. It is also suggested that the plasma environment of low altitude auroral-zone geomagnetic field lines is more suited to electron acceleration by inertial or kinetic scale Alfven waves. This is due to the self-consistent response of the electron distribution function to SAWs, which must accommodate the low altitude large-scale current system in standing waves. We characterize these effects in terms of the relative magnitude of the wave phase and electron thermal velocities. While particle trapping is shown to be significant across a wide range of plasma temperatures and wave frequencies, we find that electron beam formation in long period waves is more effective in relatively cold plasma.

Multiwavelength study of quiescent states of MrK 421 with unprecedented hard x-Ray coverage provided by NuStar in 2013

We present coordinated multiwavelength observations of the bright, nearby BL Lacertae object Mrk 421 taken in 2013 January–March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very high energy (VHE) γ-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3–79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep G » 3 power law, with no evidence for an exponential cutoff or additional hard components up "aprox" 80keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure that relates to the two bumps of the broadband SED. In each bump, the variability increases with energy, which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi band variability, the significant X-ray-toVHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in situ electron acceleration and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

Green's function for paraxial equation

The theory and experimental applications of optical Airy beams are in active development recently. The Airy beams are characterised by very special properties: they are non-diffractive and propagate along parabolic trajectories. Among the striking applications of the optical Airy beams are optical micro-manipulation implemented as the transport of small particles along the parabolic trajectory, Airy-Bessel linear light bullets, electron acceleration by the Airy beams, plasmonic energy routing. The detailed analysis of the mathematical aspects as well as physical interpretation of the electromagnetic Airy beams was done by considering the wave as a function of spatial coordinates only, related by the parabolic dependence between the transverse and the longitudinal coordinates. Their time dependence is assumed to be harmonic. Only a few papers consider a more general temporal dependence where such a relationship exists between the temporal and the spatial variables. This relationship is derived mostly by applying the Fourier transform to the expressions obtained for the harmonic time dependence or by a Fourier synthesis using the specific modulated spectrum near some central frequency. Spatial-temporal Airy pulses in the form of contour integrals is analysed near the caustic and the numerical solution of the nonlinear paraxial equation in time domain shows soliton shedding from the Airy pulse in Kerr medium. In this paper the explicitly time dependent solutions of the electromagnetic problem in the form of time-spatial pulses are derived in paraxial approximation through the Green's function for the paraxial equation. It is shown that a Gaussian and an Airy pulse can be obtained by applying the Green's function to a proper source current. We emphasize that the processes in time domain are directional, which leads to unexpected conclusions especially for the paraxial approximation.

L'impact du glissement en fréquence lors de l'accélération directe d'électrons par le faisceau laser

L’accélération directe d’électrons par des impulsions ultrabrèves de polarisation radiale fortement focalisées démontre un grand potentiel, notamment, pour la production de paquets d’électrons ultrabrefs. Plusieurs aspects de ce schéma d’accélération restent toutefois à être explorés pour en permettre une maîtrise approfondie. Dans le cadre du présent mémoire, on s’intéresse à l’ajout d’une dérive de fréquence au champ de l’impulsion TM01 utilisée. Les expressions exactes des composantes du champ électromagnétique de l’impulsion TM01 sont établies à partir d’une généralisation du spectre de Poisson. Il s’agit, à notre connaissance, du premier modèle analytique exact pour la description d’une impulsion avec une dérive de fréquence. Ce modèle est utilisé pour étudier l’impact du glissement en fréquence sur le schéma d’accélération, grâce à des simulations “particule test” unidimensionnelles, considérant en premier lieu une énergie constante par impulsion, puis un champ maximum constant. Les résultats révèlent que le glissement en fréquence diminue le gain en énergie maximum atteignable dans le cadre du schéma d’accélération à l’étude ; une baisse d’efficacité de plusieurs dizaines de pourcents peut survenir. De plus, les simulations mettent en évidence certaines différences reliées à l’utilisation d’impulsions avec une dérive vers les basses fréquences ou avec une dérive vers les hautes fréquences : il se trouve que, pour un glissement en fréquence de même grandeur, l’impulsion avec une dérive vers les basses fréquences conduit à un gain en énergie cinétique maximum plus élevé pour l’électron que l’impulsion avec une dérive vers les hautes fréquences.

Electron self-injection and acceleration driven by a tightly focused intense laser beam in an underdense plasma

A scheme for electron self-injection in the laser wakefield acceleration is proposed. In this scheme, the transverse wave breaking of the wakefield and the tightly focused geometry of the laser beam play important roles. A large number of the background electrons are self-injected into the acceleration phase of the wakefield during the defocusing of the tightly focused laser beam as it propagates through an underdense plasma. Particle-in-cell simulations performed using a 2D3V code have shown generation of a collimated electron bunch with a total number of 1.4 x 109 and energies up to 8 MeV. (C) 2005 American Institute of Physics.

Acceleration of ultra-thin electron layer.Analytical treatment compared with 1D-PIC simulation

In this paper, we apply an analytical model [V.V. Kulagin et al., Phys. Plasmas 14, 113101 (2007)] to describe the acceleration of an ultra-thin electron layer by a schematic single-cycle laser pulse and compare with one-dimensional particle-in-cell (1D-PIC) simulations. This is in the context of creating a relativistic mirror for coherent backscattering and supplements two related papers in this EPJD volume. The model is shown to reproduce the 1D-PIC results almost quantitatively for the short time of a few laser periods sufficient for the backscattering of ultra-short probe pulses.

Effects of Ultrafast Solvation on the Rate of Adiabatic Outer-Sphere Electron Transfer Reactions

Recent studies have demonstrated that solvation dynamics in many common dipolar liquids contain an initial, ultrafast Gaussian component which may contribute even more than 60% to the total solvation energy. It is also known that adiabatic electron transfer reactions often probe the high-frequency components of the relevant solvent friction (Hynes, J. T. J. Phys. Chem. 1986, 90, 3701). In this paper, we present a theoretical study of the effects of the ultrafast solvent polar modes on the adiabatic electron transfer reactions by using the formalism of Hynes. Calculations have been carried out for a model system and also for water and acetonitrile. It is found that, in general, the ultrafast modes can greatly enhance the rate of electron transfer, even by more than an order of magnitude, over the rate obtained by using only the slow overdamped modes usually considered. For water, this acceleration of the rate can be attributed to the high-frequency intermolecular vibrational and librational modes. For a weakly adiabatic reaction, the rate is virtually indistinguishable from the rate predicted by the Marcus transition state theory. Another important result is that even in this case of ultrafast underdamped solvation, energy diffusion appears to be efficient so that electron transfer reaction in water is controlled essentially by the barrier crossing dynamics. This is because the reactant well frequency is-directly proportional to the rate of the initial Gaussian decay of the solvation time correlation function. As a result, the value of the friction at the reactant well frequency rarely falls below the value required for the Kramers turnover except when the polarizability of the water molecules may be neglected. On the other hand, in acetonitrile, the rate of electron transfer reaction is found to be controlled by the energy diffusion dynamics, although a significant contribution to the rate comes also from the barrier crossing rate. Therefore, the present study calls for a need to understand the relaxation of the high-frequency modes in dipolar liquids.

Parallel Fully-Implicit Computation of MHD Acceleration Experiments

A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETSc package. The results of convergence tests are plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.