Dipole Oscillations of a Bose-Einstein Condensate in Presence of Defects and Disorder

We consider dipole oscillations of a trapped dilute Bose-Einstein condensate in the presence of a scattering potential consisting either in a localized defect or in an extended disordered potential. In both cases the breaking of superfluidity and the damping of the oscillations are shown to be related to the appearance of a nonlinear dissipative flow. At supersonic velocities the flow becomes asymptotically dissipationless.Comment: 4 pages, 4 figure

Band Distributions for Quantum Chaos on the Torus

Band distributions (BDs) are introduced describing quantization in a toral phase space. A BD is the uniform average of an eigenstate phase-space probability distribution over a band of toral boundary conditions. A general explicit expression for the Wigner BD is obtained. It is shown that the Wigner functions for {\em all} of the band eigenstates can be reproduced from the Wigner BD. Also, BDs are shown to be closer to classical distributions than eigenstate distributions. Generalized BDs, associated with sets of adjacent bands, are used to extend in a natural way the Chern-index characterization of the classical-quantum correspondence on the torus to arbitrary rational values of the scaled Planck constant.Comment: 12 REVTEX page

Average ground-state energy of finite Fermi systems

Semiclassical theories like the Thomas-Fermi and Wigner-Kirkwood methods give a good description of the smooth average part of the total energy of a Fermi gas in some external potential when the chemical potential is varied. However, in systems with a fixed number of particles N, these methods overbind the actual average of the quantum energy as N is varied. We describe a theory that accounts for this effect. Numerical illustrations are discussed for fermions trapped in a harmonic oscillator potential and in a hard wall cavity, and for self-consistent calculations of atomic nuclei. In the latter case, the influence of deformations on the average behavior of the energy is also considered.Comment: 10 pages, 8 figure

Dual equilibrium in a finite aspect ratio tokamak

A new approach to high pressure magnetically-confined plasmas is necessary to design efficient fusion devices. This paper presents an equilibrium combining two solutions of the Grad-Shafranov equation, which describes the magnetohydrodynamic equilibrium in toroidal geometry. The outer equilibrium is paramagnetic and confines the inner equilibrium, whose strong diamagnetism permits to balance large pressure gradients. The existence of both equilibria in the same volume yields a dual equilibrium structure. Their combination also improves free-boundary mode stability

Fermi-surface reconstruction and two-carrier model for the Hall effect in YBa2Cu4O8

Pulsed field measurements of the Hall resistivity and magnetoresistance of underdoped YBa2Cu4O8 are analyzed self-consistently using a simple model based on coexisting electron and hole carriers. The resultant mobilities and Hall numbers are found to vary markedly with temperature. The conductivity of the hole carriers drops by one order of magnitude below 30 K, explaining the absence of quantum oscillations from these particular pockets. Meanwhile the Hall coefficient of the electron carriers becomes strongly negative below 50 K. The overall quality of the fits not only provides strong evidence for Fermi-surface reconstruction in Y-based cuprates, it also strongly constrains the type of reconstruction that might be occurring.Comment: 5 pages, 4 figures, updated after publication in Physical Review B (Rapid Communication

Quantum Chaotic Dynamics and Random Polynomials

We investigate the distribution of roots of polynomials of high degree with random coefficients which, among others, appear naturally in the context of "quantum chaotic dynamics". It is shown that under quite general conditions their roots tend to concentrate near the unit circle in the complex plane. In order to further increase this tendency, we study in detail the particular case of self-inversive random polynomials and show that for them a finite portion of all roots lies exactly on the unit circle. Correlation functions of these roots are also computed analytically, and compared to the correlations of eigenvalues of random matrices. The problem of ergodicity of chaotic wave-functions is also considered. For that purpose we introduce a family of random polynomials whose roots spread uniformly over phase space. While these results are consistent with random matrix theory predictions, they provide a new and different insight into the problem of quantum ergodicity. Special attention is devoted all over the paper to the role of symmetries in the distribution of roots of random polynomials.Comment: 33 pages, Latex, 6 Figures not included (a copy of them can be requested at [email protected]); to appear in Journal of Statistical Physic

Implementation of 2D Domain Decomposition in the UCAN Gyrokinetic Particle-in-Cell Code and Resulting Performance of UCAN2

The massively parallel, nonlinear, three-dimensional (3D), toroidal, electrostatic, gyrokinetic, particle-in-cell (PIC), Cartesian geometry UCAN code, with particle ions and adiabatic electrons, has been successfully exercised to identify non-diffusive transport characteristics in present day tokamak discharges. The limitation in applying UCAN to larger scale discharges is the 1D domain decomposition in the toroidal (or z-) direction for massively parallel implementation using MPI which has restricted the calculations to a few hundred ion Larmor radii or gyroradii per plasma minor radius. To exceed these sizes, we have implemented 2D domain decomposition in UCAN with the addition of the y-direction to the processor mix. This has been facilitated by use of relevant components in the P2LIB library of field and particle management routines developed for UCLA's UPIC Framework of conventional PIC codes. The gyro-averaging specific to gyrokinetic codes is simplified by the use of replicated arrays for efficient charge accumulation and force deposition. The 2D domain-decomposed UCAN2 code reproduces the original 1D domain nonlinear results within round-off. Benchmarks of UCAN2 on the Cray XC30 Edison at NERSC demonstrate ideal scaling when problem size is increased along with processor number up to the largest power of 2 available, namely 131,072 processors. These particle weak scaling benchmarks also indicate that the 1 nanosecond per particle per time step and 1 TFlops barriers are easily broken by UCAN2 with 1 billion particles or more and 2000 or more processors.This work was supported in part in the USA by Grant No. DE-FG02-04ER54741 to the University of Alaska, Fairbanks, AK, from the Office of Fusion Energy Sciences, Office of Science, United States Department of Energy. It was also supported in part at Universidad Carlos III, Madrid, Spain, by Spanish National Project No. ENE2009-12213-C03-03. This research used resources of the National Energy Research Scientific Computing Center (NERSC), which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. It also took advantage of resources at the Barcelona Supercomputing Center (BSC), Centro Nacional de Supercomputación, Barcelona, Spain. One of us (Leboeuf) would particularly like to thank David Vicente from BSC and Zhengji Zhao from NERSC for their help in the porting, debugging, and optimization of UCAN2 on the mainframes at their respective centers

Nature of Transport across Sheared Zonal Flows in Electrostatic Ion-Temperature-Gradient Gyrokinetic Plasma Turbulence

4 pages, 4 figures.-- PACS nrs.: 52.35.Ra, 05.40.Fb, 52.55.Fa, 52.65.Tt.It is shown that the usual picture for the suppression of turbulent transport across a stable sheared flow based on a reduction of diffusive transport coefficients is, by itself, incomplete. By means of toroidal gyrokinetic simulations of electrostatic, collisionless ion-temperature-gradient turbulence, it is found that the nature of the transport is altered fundamentally, changing from diffusive to anticorrelated and subdiffusive. Additionally, whenever the flows are self-consistently driven by turbulence, the transport gains an additional non-Gaussian character. These results suggest that a description of transport across sheared flows using effective diffusivities is oversimplified.Research carried out at ORNL, managed by UT-Battelle LLC, for US DOE under Contract No. DE-AC05-00OR22725. Research funded by DOE Office of Science Grants No. DE-FG02-04ER54741 at University of Alaska and No. DE-FG02-04ER54740 at UCLA.Publicad

Enhancement of the Nernst effect by stripe order in a high-Tc superconductor

The Nernst effect in metals is highly sensitive to two kinds of phase transition: superconductivity and density-wave order. The large positive Nernst signal observed in hole-doped high-Tc superconductors above their transition temperature Tc has so far been attributed to fluctuating superconductivity. Here we show that in some of these materials the large Nernst signal is in fact caused by stripe order, a form of spin / charge modulation which causes a reconstruction of the Fermi surface. In LSCO doped with Nd or Eu, the onset of stripe order causes the Nernst signal to go from small and negative to large and positive, as revealed either by lowering the hole concentration across the quantum critical point in Nd-LSCO, or lowering the temperature across the ordering temperature in Eu-LSCO. In the latter case, two separate peaks are resolved, respectively associated with the onset of stripe order at high temperature and superconductivity near Tc. This sensitivity to Fermi-surface reconstruction makes the Nernst effect a promising probe of broken symmetry in high-Tc superconductors

Shubnikov-de Haas oscillations in YBa_2Cu_4O_8

We report the observation of Shubnikov-de Haas oscillations in the underdoped cuprate superconductor YBa$_2$Cu$_4$O$_8$ (Y124). For field aligned along the c-axis, the frequency of the oscillations is $660\pm 30$ T, which corresponds to $\sim 2.4$ % of the total area of the first Brillouin zone. The effective mass of the quasiparticles on this orbit is measured to be $2.7\pm0.3$ times the free electron mass. Both the frequency and mass are comparable to those recently observed for ortho-II YBa$_2$Cu$_3$O$_{6.5}$ (Y123-II). We show that although small Fermi surface pockets may be expected from band structure calculations in Y123-II, no such pockets are predicted for Y124. Our results therefore imply that these small pockets are a generic feature of the copper oxide plane in underdoped cuprates.Comment: v2: Version of paper accepted for publication in Physical Review Letters. Only minor changes to the text and reference