An extended hybrid magnetohydrodynamics gyrokinetic model for numerical simulation of shear Alfv\'en waves in burning plasmas

Adopting the theoretical framework for the generalized fishbonelike dispersion relation, an extended hybrid magnetohydrodynamics gyrokinetic simulation model has been derived analytically by taking into account both thermal ion compressibility and diamagnetic effects in addition to energetic particle kinetic behaviors. The extended model has been used for implementing an eXtended version of Hybrid Magnetohydrodynamics Gyrokinetic Code (XHMGC) to study thermal ion kinetic effects on Alfv\'enic modes driven by energetic particles, such as kinetic beta induced Alfv\'en eigenmodes in tokamak fusion plasmas

Weakly collisional Landau damping and three-dimensional Bernstein-Greene-Kruskal modes: New results on old problems

Landau damping and Bernstein-Greene-Kruskal (BGK) modes are among the most fundamental concepts in plasma physics. While the former describes the surprising damping of linear plasma waves in a collisionless plasma, the latter describes exact undamped nonlinear solutions of the Vlasov equation. There does exist a relationship between the two: Landau damping can be described as the phase-mixing of undamped eigenmodes, the so-called Case-Van Kampen modes, which can be viewed as BGK modes in the linear limit. While these concepts have been around for a long time, unexpected new results are still being discovered. For Landau damping, we show that the textbook picture of phase-mixing is altered profoundly in the presence of collision. In particular, the continuous spectrum of Case-Van Kampen modes is eliminated and replaced by a discrete spectrum, even in the limit of zero collision. Furthermore, we show that these discrete eigenmodes form a complete set of solutions. Landau-damped solutions are then recovered as true eigenmodes (which they are not in the collisionless theory). For BGK modes, our interest is motivated by recent discoveries of electrostatic solitary waves in magnetospheric plasmas. While one-dimensional BGK theory is quite mature, there appear to be no exact three-dimensional solutions in the literature (except for the limiting case when the magnetic field is sufficiently strong so that one can apply the guiding-center approximation). We show, in fact, that two- and three-dimensional solutions that depend only on energy do not exist. However, if solutions depend on both energy and angular momentum, we can construct exact three-dimensional solutions for the unmagnetized case, and two-dimensional solutions for the case with a finite magnetic field. The latter are shown to be exact, fully electromagnetic solutions of the steady-state Vlasov-Poisson-Amp\`ere system

Free streaming in mixed dark matter

Free streaming in a \emph{mixture} of collisionless non-relativistic dark matter (DM) particles is studied by implementing methods from the theory of multicomponent plasmas. The mixture includes Fermionic, condensed and non condensed Bosonic particles decoupling in equilibrium while relativistic, heavy non-relativistic thermal relics (WIMPs), and sterile neutrinos that decouple \emph{out of equilibrium} when they are relativistic. The free-streaming length $\lambda_{fs}$ is obtained from the marginal zero of the gravitational polarization function, which separates short wavelength Landau-damped from long wavelength Jeans-unstable \emph{collective} modes. At redshift $z$ we find $\frac{1}{\lambda^2_{fs}(z)}= \frac{1}{(1+z)} \big[\frac{0.071}{\textrm{kpc}} \big]^2 \sum_{a}\nu_a g^{2/3}_{d,a}({m_a}/{\mathrm{keV}})^2 I_a$,where $0\leq \nu_a \leq 1$ are the \emph{fractions} of the respective DM components of mass $m_a$ that decouple when the effective number of ultrarelativistic degrees of freedom is $g_{d,a}$, and $I_a$ only depend on the distribution functions at decoupling, given explicitly in all cases. If sterile neutrinos produced either resonantly or non-resonantly that decouple near the QCD scale are the \emph{only} DM component,we find $\lambda_{fs}(0) \simeq 7 \mathrm{kpc} (\mathrm{keV}/m)$ (non-resonant), $\lambda_{fs}(0) \simeq 1.73 \mathrm{kpc} (\mathrm{keV}/m)$ (resonant).If WIMPs with $m_{wimp} \gtrsim 100 \mathrm{GeV}$ decoupling at $T_d \gtrsim 10 \mathrm{MeV}$ are present in the mixture with $\nu_{wimp} \gg 10^{-12}$,$\lambda_{fs}(0) \lesssim 6.5 \times 10^{-3} \mathrm{pc}$ is \emph{dominated} by CDM. If a Bose Einstein condensate is a DM component its free streaming length is consistent with CDM because of the infrared enhancement of the distribution function.Comment: 19 pages, 2 figures. More discussions same conclusions and results. Version to appear in Phys. Rev.

Understanding the effect of sheared flow on microinstabilities

The competition between the drive and stabilization of plasma microinstabilities by sheared flow is investigated, focusing on the ion temperature gradient mode. Using a twisting mode representation in sheared slab geometry, the characteristic equations have been formulated for a dissipative fluid model, developed rigorously from the gyrokinetic equation. They clearly show that perpendicular flow shear convects perturbations along the field at a speed we denote by $Mc_s$ (where $c_s$ is the sound speed), whilst parallel flow shear enters as an instability driving term analogous to the usual temperature and density gradient effects. For sufficiently strong perpendicular flow shear, $M >1$, the propagation of the system characteristics is unidirectional and no unstable eigenmodes may form. Perturbations are swept along the field, to be ultimately dissipated as they are sheared ever more strongly. Numerical studies of the equations also reveal the existence of stable regions when $M < 1$, where the driving terms conflict. However, in both cases transitory perturbations exist, which could attain substantial amplitudes before decaying. Indeed, for $M \gg 1$, they are shown to exponentiate $\sqrt{M}$ times. This may provide a subcritical route to turbulence in tokamaks.Comment: minor revisions; accepted to PPC

Local and global Fokker-Planck neoclassical calculations showing flow and bootstrap current modification in a pedestal

In transport barriers, particularly H-mode edge pedestals, radial scale lengths can become comparable to the ion orbit width, causing neoclassical physics to become radially nonlocal. In this work, the resulting changes to neoclassical flow and current are examined both analytically and numerically. Steep density gradients are considered, with scale lengths comparable to the poloidal ion gyroradius, together with strong radial electric fields sufficient to electrostatically confine the ions. Attention is restricted to relatively weak ion temperature gradients (but permitting arbitrary electron temperature gradients), since in this limit a delta-f (small departures from a Maxwellian distribution) rather than full-f approach is justified. This assumption is in fact consistent with measured inter-ELM H-Mode edge pedestal density and ion temperature profiles in many present experiments, and is expected to be increasingly valid in future lower collisionality experiments. In the numerical analysis, the distribution function and Rosenbluth potentials are solved for simultaneously, allowing use of the exact field term in the linearized Fokker-Planck collision operator. In the pedestal, the parallel and poloidal flows are found to deviate strongly from the best available conventional neoclassical prediction, with large poloidal variation of a different form than in the local theory. These predicted effects may be observable experimentally. In the local limit, the Sauter bootstrap current formulae appear accurate at low collisionality, but they can overestimate the bootstrap current near the plateau regime. In the pedestal ordering, ion contributions to the bootstrap and Pfirsch-Schluter currents are also modified

Markov Properties of Electrical Discharge Current Fluctuations in Plasma

Using the Markovian method, we study the stochastic nature of electrical discharge current fluctuations in the Helium plasma. Sinusoidal trends are extracted from the data set by the Fourier-Detrended Fluctuation analysis and consequently cleaned data is retrieved. We determine the Markov time scale of the detrended data set by using likelihood analysis. We also estimate the Kramers-Moyal's coefficients of the discharge current fluctuations and derive the corresponding Fokker-Planck equation. In addition, the obtained Langevin equation enables us to reconstruct discharge time series with similar statistical properties compared with the observed in the experiment. We also provide an exact decomposition of temporal correlation function by using Kramers-Moyal's coefficients. We show that for the stationary time series, the two point temporal correlation function has an exponential decaying behavior with a characteristic correlation time scale. Our results confirm that, there is no definite relation between correlation and Markov time scales. However both of them behave as monotonic increasing function of discharge current intensity. Finally to complete our analysis, the multifractal behavior of reconstructed time series using its Keramers-Moyal's coefficients and original data set are investigated. Extended self similarity analysis demonstrates that fluctuations in our experimental setup deviates from Kolmogorov (K41) theory for fully developed turbulence regime.Comment: 25 pages, 9 figures and 4 tables. V3: Added comments, references, figures and major correction

Analytic fluid theory of beam spiraling in high-intensity cyclotrons

Using a two-dimensional fluid description, we investigate the nonlinear radial-longitudinal dynamics of intense beams in isochronous cyclotrons in the nonrelativistic limit. With a multiscale analysis separating the time scale associated with the betatron motion and the slower time scale associated with space-charge effects, we show that the longitudinal-radial vortex motion can be understood in the frame moving with the charged beam as the nonlinear advection of the beam by the E×B velocity field, where E is the electric field due to the space charge and B is the external magnetic field. This interpretation provides simple explanations for the stability of round beams and for the development of spiral halos in elongated beams. By numerically solving the nonlinear advection equation for the beam density, we find that it is also in quantitative agreement with results obtained in particle-in-cell simulations

The Spatial and Temporal Construction of Confidence in the Visual Scene

Human subjects can report many items of a cluttered field a few hundred milliseconds after stimulus presentation. This memory decays rapidly and after a second only 3 or 4 items can be stored in working memory. Here we compared the dynamics of objective performance with a measure of subjective report and we observed that 1) Objective performance beyond explicit subjective reports (blindsight) was significantly more pronounced within a short temporal interval and within specific locations of the visual field which were robust across sessions 2) High confidence errors (false beliefs) were largely confined to a small spatial window neighboring the cue. The size of this window did not change in time 3) Subjective confidence showed a moderate but consistent decrease with time, independent of all other experimental factors. Our study allowed us to asses quantitatively the temporal and spatial access to an objective response and to subjective reports

Carbon-Nanotube-Embedded Hydrogel Sheets for Engineering Cardiac Constructs and Bioactuators

We engineered functional cardiac patches by seeding neonatal rat cardiomyocytes onto carbon nanotube (CNT)-incorporated photo-cross-linkable gelatin methacrylate (GelMA) hydrogels. The resulting cardiac constructs showed excellent mechanical integrity and advanced electrophysiological functions. Specifically, myocardial tissues cultured on 50 μm thick CNT-GelMA showed 3 times higher spontaneous synchronous beating rates and 85% lower excitation threshold, compared to those cultured on pristine GelMA hydrogels. Our results indicate that the electrically conductive and nanofibrous networks formed by CNTs within a porous gelatin framework are the key characteristics of CNT-GelMA leading to improved cardiac cell adhesion, organization, and cell–cell coupling. Centimeter-scale patches were released from glass substrates to form 3D biohybrid actuators, which showed controllable linear cyclic contraction/extension, pumping, and swimming actuations. In addition, we demonstrate for the first time that cardiac tissues cultured on CNT-GelMA resist damage by a model cardiac inhibitor as well as a cytotoxic compound. Therefore, incorporation of CNTs into gelatin, and potentially other biomaterials, could be useful in creating multifunctional cardiac scaffolds for both therapeutic purposes and in vitro studies. These hybrid materials could also be used for neuron and other muscle cells to create tissue constructs with improved organization, electroactivity, and mechanical integrity.United States. Army Research Office. Institute for Soldier NanotechnologiesNational Institutes of Health (U.S.) (HL092836)National Institutes of Health (U.S.) (EB02597)National Institutes of Health (U.S.) (AR057837)National Institutes of Health (U.S.) (HL099073)National Science Foundation (U.S.) (DMR0847287)United States. Office of Naval Research (ONR PECASE Award)United States. Office of Naval Research (Young Investigator award)National Research Foundation of Korea (grant (NRF-2010-220-D00014)