Strong "quantum" chaos in the global ballooning mode spectrum of three-dimensional plasmas

The spectrum of ideal magnetohydrodynamic (MHD) pressure-driven (ballooning) modes in strongly nonaxisymmetric toroidal systems is difficult to analyze numerically owing to the singular nature of ideal MHD caused by lack of an inherent scale length. In this paper, ideal MHD is regularized by using a $k$-space cutoff, making the ray tracing for the WKB ballooning formalism a chaotic Hamiltonian billiard problem. The minimum width of the toroidal Fourier spectrum needed for resolving toroidally localized ballooning modes with a global eigenvalue code is estimated from the Weyl formula. This phase-space-volume estimation method is applied to two stellarator cases.Comment: 4 pages typeset, including 2 figures. Paper accepted for publication in Phys. Rev. Letter

A comparison of incompressible limits for resistive plasmas

The constraint of incompressibility is often used to simplify the magnetohydrodynamic (MHD) description of linearized plasma dynamics because it does not affect the ideal MHD marginal stability point. In this paper two methods for introducing incompressibility are compared in a cylindrical plasma model: In the first method, the limit $\gamma \to \infty$ is taken, where $\gamma$ is the ratio of specific heats; in the second, an anisotropic mass tensor $\mathbf{\rho}$ is used, with the component parallel to the magnetic field taken to vanish, $\rho_{\parallel} \to 0$. Use of resistive MHD reveals the nature of these two limits because the Alfv\'en and slow magnetosonic continua of ideal MHD are converted to point spectra and moved into the complex plane. Both limits profoundly change the slow-magnetosonic spectrum, but only the second limit faithfully reproduces the resistive Alfv\'en spectrum and its wavemodes. In ideal MHD, the slow magnetosonic continuum degenerates to the Alfv\'en continuum in the first method, while it is moved to infinity by the second. The degeneracy in the first is broken by finite resistivity. For numerical and semi-analytical study of these models, we choose plasma equilibria which cast light on puzzling aspects of results found in earlier literature.Comment: 14 pages, 10 figure

Singularity theory study of overdetermination in models for L-H transitions

Two dynamical models that have been proposed to describe transitions between low and high confinement states (L-H transitions) in confined plasmas are analysed using singularity theory and stability theory. It is shown that the stationary-state bifurcation sets have qualitative properties identical to standard normal forms for the pitchfork and transcritical bifurcations. The analysis yields the codimension of the highest-order singularities, from which we find that the unperturbed systems are overdetermined bifurcation problems and derive appropriate universal unfoldings. Questions of mutual equivalence and the character of the state transitions are addressed.Comment: Latex (Revtex) source + 13 small postscript figures. Revised versio

Nonequilibrium statistical mechanics of shear flow: invariant quantities and current relations

In modeling nonequilibrium systems one usually starts with a definition of the microscopic dynamics, e.g., in terms of transition rates, and then derives the resulting macroscopic behavior. We address the inverse question for a class of steady state systems, namely complex fluids under continuous shear flow: how does an externally imposed shear current affect the microscopic dynamics of the fluid? The answer can be formulated in the form of invariant quantities, exact relations for the transition rates in the nonequilibrium steady state, as discussed in a recent letter [A. Baule and R. M. L. Evans, Phys. Rev. Lett. 101, 240601 (2008)]. Here, we present a more pedagogical account of the invariant quantities and the theory underlying them, known as the nonequilibrium counterpart to detailed balance (NCDB). Furthermore, we investigate the relationship between the transition rates and the shear current in the steady state. We show that a fluctuation relation of the Gallavotti-Cohen type holds for systems satisfying NCDB.Comment: 24 pages, 11 figure

Turbulent edge structure formation in complex configurations

Formation of nonlinear structures in drift-Alfvén turbulence is investigated in the often complex edge geometries of stellarator and tokamak configurations, by analysis of drift waveturbulence simulations using a model in which three-dimensional magnetic geometries are approximated. The structures of parallel mode extension, radially sheared zonal flows and perpendicular mode spectra are highlighted in particular for three-dimensional stellaratormagnetic fields and shaped tokamaks. Specific characteristics of advanced stellarators in comparison to (lower aspect ratio) circular tokamaks are a less pronounced ballooning structure of the modes, a strong influence of local magnetic shear on amplitude structure and average, and stronger level of zonal flows due to lower geodesic curvature.This work was partly funded by grants within the ‘‘Australian-German Joint Research Co-operation scheme’’ ~PPP project no. D/0205403!

Lattice-gas model for alkali-metal fullerides: face-centered-cubic structure

A lattice-gas model is suggested for describing the ordering phenomena in alkali-metal fullerides of face-centered-cubic structure assuming the electric charge of alkali ions residing in either octahedral or tetrahedral interstitial sites is completely screened by the first-neighbor C_60 molecules. This approximation allows us to derive an effective ion-ion interaction. The van der Waals interaction between the ion and C_60 molecule is characterized by introducing an additional energy at the tetrahedral sites. This model is investigated by using a three-sublattice mean-field approximation and a simple cluster-variation method. The analysis shows a large variety of phase diagrams when changing the site energy parameter.Comment: 10 twocolumn pages (REVTEX) including 12 PS figure

Model Data Fusion: developing Bayesian inversion to constrain equilibrium and mode structure

Recently, a new probabilistic "data fusion" framework based on Bayesian principles has been developed on JET and W7-AS. The Bayesian analysis framework folds in uncertainties and inter-dependencies in the diagnostic data and signal forward-models, together with prior knowledge of the state of the plasma, to yield predictions of internal magnetic structure. A feature of the framework, known as MINERVA (J. Svensson, A. Werner, Plasma Physics and Controlled Fusion 50, 085022, 2008), is the inference of magnetic flux surfaces without the use of a force balance model. We discuss results from a new project to develop Bayesian inversion tools that aim to (1) distinguish between competing equilibrium theories, which capture different physics, using the MAST spherical tokamak; and (2) test the predictions of MHD theory, particularly mode structure, using the H-1 Heliac.Comment: submitted to Journal of Plasma Fusion Research 10/11/200

Anderson localization of ballooning modes, quantum chaos and the stability of compact quasiaxially symmetric stellarators

The radially local magnetohydrodynamic(MHD) ballooning stability of a compact, quasiaxially symmetric stellarator (QAS), is examined just above the ballooning beta limit with a method that can lead to estimates of global stability. Here MHDstability is analyzed through the calculation and examination of the ballooning modeeigenvalue isosurfaces in the 3-space (s,α,θk); s is the edge normalized toroidal flux, α is the field linevariable, and θk is the perpendicular wave vector or ballooning parameter. Broken symmetry, i.e., deviations from axisymmetry, in the stellarator magnetic field geometry causes localization of the ballooning mode eigenfunction, and gives rise to new types of nonsymmetric eigenvalue isosurfaces in both the stable and unstable spectrum. For eigenvalues far above the marginal point, isosurfaces are topologically spherical, indicative of strong “quantum chaos.” The complexity of QAS marginal isosurfaces suggests that finite Larmor radius stabilization estimates will be difficult and that fully three-dimensional, high-nMHD computations are required to predict the beta limit.Research supported by U.S. DOE Contract No. DEAC02-76CH0373. John Canik held a U.S. DOE National Undergraduate Fellowship at Princeton Plasma Physics Laboratory, during the summer of 2000

Automated detection of atrial fibrillation using RR intervals and multivariate-based classification

Automated detection of AF from the electrocardiogram (ECG) still remains a challenge. In this study, we investigated two multivariate-based classification techniques, Random Forests (RF) and k-nearest neighbor (k-nn), for improved automated detection of AF from the ECG. We have compiled a new database from ECG data taken from existing sources. R-R intervals were then analyzed using four previously described R-R irregularity measurements: (1) the coefficient of sample entropy (CoSEn), (2) the coefficient of variance (CV), (3) root mean square of the successive differences (RMSSD), and (4) median absolute deviation (MAD). Using outputs from all four R-R irregularity measurements, RF and k-nn models were trained. RF classification improved AF detection over CoSEn with overall specificity of 80.1% vs. 98.3% and positive predictive value of 51.8% vs. 92.1% with a reduction in sensitivity, 97.6% vs. 92.8%. k-nn also improved specificity and PPV over CoSEn; however, the sensitivity of this approach was considerably reduced (68.0%)