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

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!

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

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

Randomised comparison of 5 years of adjuvant tamoxifen with continuous therapy for operable breast cancer. The Scottish Cancer Trials Breast Group.

In 1985 a second randomisation was initiated for women in the treatment arm of the Scottish Tamoxifen Trial either to stop tamoxifen at 5 years or to continue indefinitely. A preliminary analysis of outcome in 342 patients at a median follow-up of 6 years suggests that a worthwhile gain in disease control from continuing adjuvant tamoxifen beyond 5 years is unlikely. [Hazard ratio for events (relapse or death without relapse) is 1.27, 95% CI = 0.87 - 1.85.] There is a suggestion that therapy for longer than 5 years may increase the risk of endometrial carcinoma (P = 0.064)

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

Nonlinear saturation of electrostatic waves: mobile ions modify trapping scaling

The amplitude equation for an unstable electrostatic wave in a multi-species Vlasov plasma has been derived. The dynamics of the mode amplitude $\rho(t)$ is studied using an expansion in $\rho$; in particular, in the limit $\gamma\rightarrow0^+$, the singularities in the expansion coefficients are analyzed to predict the asymptotic dependence of the electric field on the linear growth rate $\gamma$. Generically $|E_k|\sim \gamma^{5/2}$, as $\gamma\rightarrow0^+$, but in the limit of infinite ion mass or for instabilities in reflection-symmetric systems due to real eigenvalues the more familiar trapping scaling $|E_k|\sim \gamma^{2}$ is predicted.Comment: 13 pages (Latex/RevTex), 4 postscript encapsulated figures which are included using the utility "uufiles". They should be automatically included with the text when it is downloaded. Figures also available in hard copy from the authors ([email protected]

Predictive use of the Maximum Entropy Production principle for Past and Present Climates

In this paper, we show how the MEP hypothesis may be used to build simple climate models without representing explicitly the energy transport by the atmosphere. The purpose is twofold. First, we assess the performance of the MEP hypothesis by comparing a simple model with minimal input data to a complex, state-of-the-art General Circulation Model. Next, we show how to improve the realism of MEP climate models by including climate feedbacks, focusing on the case of the water-vapour feedback. We also discuss the dependence of the entropy production rate and predicted surface temperature on the resolution of the model

Universal trapping scaling on the unstable manifold for a collisionless electrostatic mode

An amplitude equation for an unstable mode in a collisionless plasma is derived from the dynamics on the two-dimensional unstable manifold of the equilibrium. The mode amplitude $\rho(t)$ decouples from the phase due to the spatial homogeneity of the equilibrium, and the resulting one-dimensional dynamics is analyzed using an expansion in $\rho$. As the linear growth rate $\gamma$ vanishes, the expansion coefficients diverge; a rescaling $\rho(t)\equiv\gamma^2\,r(\gamma t)$ of the mode amplitude absorbs these singularities and reveals that the mode electric field exhibits trapping scaling $|E_1|\sim\gamma^2$ as $\gamma\rightarrow0$. The dynamics for $r(\tau)$ depends only on the phase $e^{i\xi}$ where $d\epsilon_{{k}} /dz=|{\epsilon_{{k}}}|e^{-i\xi/2}$ is the derivative of the dielectric as $\gamma\rightarrow0$.Comment: 11 pages (Latex/RevTex), 2 figures available in hard copy from the Author ([email protected]); paper accepted by Physical Review Letter