Landau Collision Integral Solver with Adaptive Mesh Refinement on Emerging Architectures

The Landau collision integral is an accurate model for the small-angle dominated Coulomb collisions in fusion plasmas. We investigate a high order accurate, fully conservative, finite element discretization of the nonlinear multi-species Landau integral with adaptive mesh refinement using the PETSc library (www.mcs.anl.gov/petsc). We develop algorithms and techniques to efficiently utilize emerging architectures with an approach that minimizes memory usage and movement and is suitable for vector processing. The Landau collision integral is vectorized with Intel AVX-512 intrinsics and the solver sustains as much as 22% of the theoretical peak flop rate of the Second Generation Intel Xeon Phi, Knights Landing, processor

Intraspecific trait variation and coordination: Root and leaf economics spectra in coffee across environmental gradients

Hypotheses on the existence of a universal “Root Economics Spectrum” (RES) have received arguably the least attention of all trait spectra, despite the key role root trait variation plays in resource acquisition potential. There is growing interest in quantifying intraspecific trait variation (ITV) in plants, but there are few studies evaluating (i) the existence of an intraspecific RES within a plant species, or (ii) how a RES may be coordinated with other trait spectra within species, such as a leaf economics spectrum (LES). Using Coffea arabica (Rubiaceae) as a model species, we measured seven morphological and chemical traits of intact lateral roots, which were paired with information on four key LES traits. Field collections were completed across four nested levels of biological organization. The intraspecific trait coefficient of variation (cv) ranged from 25 to 87% with root diameter and specific root tip density showing the lowest and highest cv, respectively. Between 27 and 68% of root ITV was explained by site identity alone for five of the seven traits measured. A single principal component explained 56.2% of root trait covariation, with plants falling along a RES from resource acquiring to conserving traits. Multiple factor analysis revealed significant orthogonal relationships between root and leaf spectra. RES traits were strongly orthogonal with respect to LES traits, suggesting these traits vary independently from one another in response to environmental cues. This study provides among the first evidence that plants from the same species differentiate from one another along an intraspecific RES. We find that in one of the world's most widely cultivated crops, an intraspecific RES is orthogonal to an intraspecific LES, indicating that above and belowground responses of plants to managed (or natural) environmental gradients are likely to occur independently from one another. (Résumé d'auteur

Limitations of Quantum Simulation Examined by Simulating a Pairing Hamiltonian using Nuclear Magnetic Resonance

Quantum simulation uses a well-known quantum system to predict the behavior of another quantum system. Certain limitations in this technique arise, however, when applied to specific problems, as we demonstrate with a theoretical and experimental study of an algorithm to find the low-lying spectrum of a Hamiltonian. While the number of elementary quantum gates does scale polynomially with the size of the system, it increases inversely to the desired error bound $\epsilon$. Making such simulations robust to decoherence using fault-tolerance constructs requires an additional factor of $1/ \epsilon$ gates. These constraints are illustrated by using a three qubit nuclear magnetic resonance system to simulate a pairing Hamiltonian, following the algorithm proposed by Wu, Byrd, and Lidar.Comment: 6 pages, 2 eps figure

Universal Baxterization for Z\mathbb{Z}-graded Hopf algebras

We present a method for Baxterizing solutions of the constant Yang-Baxter equation associated with $\mathbb{Z}$-graded Hopf algebras. To demonstrate the approach, we provide examples for the Taft algebras and the quantum group $U_q[sl(2)]$.Comment: 8 page

The effects of bell vibrations on the sound of the modern trumpet

The acoustic spectrum of a modern trumpet with the bell section heavily damped has been compared to the spectrum of the same instrument with the bell section left free to vibrate. Measurements of the amplitude of vibration indicate that the damping significantly reduces the movement of the metal, and a corresponding change in the acoustic spectrum between the two cases is found. It is shown that the relative power in the fundamental may change by more than 3 dB when the vibrations in the bell section are damped. Two possible causes for the effects are considered: a change in input impedance, and the transfer of mechanical vibrations through the instrument to the lips. Results of modelling and experiments are presented that indicate the latter is the more plausible explanation; however, the etiology of the effect is still unknown

Density Waves Inside Inner Lindblad Resonance: Nuclear Spirals in Disk Galaxies

We analyze formation of grand-design two-arm spiral structure in the nuclear regions of disk galaxies. Such morphology has been recently detected in a number of objects using high-resolution near-infrared observations. Motivated by the observed (1) continuity between the nuclear and kpc-scale spiral structures, and by (2) low arm-interarm contrast, we apply the density wave theory to explain the basic properties of the spiral nuclear morphology. In particular, we address the mechanism for the formation, maintenance and the detailed shape of nuclear spirals. We find, that the latter depends mostly on the shape of the underlying gravitational potential and the sound speed in the gas. Detection of nuclear spiral arms provides diagnostics of mass distribution within the central kpc of disk galaxies. Our results are supported by 2D numerical simulations of gas response to the background gravitational potential of a barred stellar disk. We investigate the parameter space allowed for the formation of nuclear spirals using a new method for constructing a gravitational potential in a barred galaxy, where positions of resonances are prescribed.Comment: 18 pages, 9 figures, higher resolution available at http://www.pa.uky.edu/~ppe/papers/nucsp.ps.g

Randomized benchmarking of single and multi-qubit control in liquid-state NMR quantum information processing

Being able to quantify the level of coherent control in a proposed device implementing a quantum information processor (QIP) is an important task for both comparing different devices and assessing a device's prospects with regards to achieving fault-tolerant quantum control. We implement in a liquid-state nuclear magnetic resonance QIP the randomized benchmarking protocol presented by Knill et al (PRA 77: 012307 (2008)). We report an error per randomized $\frac{\pi}{2}$ pulse of $1.3 \pm 0.1 \times 10^{-4}$ with a single qubit QIP and show an experimentally relevant error model where the randomized benchmarking gives a signature fidelity decay which is not possible to interpret as a single error per gate. We explore and experimentally investigate multi-qubit extensions of this protocol and report an average error rate for one and two qubit gates of $4.7 \pm 0.3 \times 10^{-3}$ for a three qubit QIP. We estimate that these error rates are still not decoherence limited and thus can be improved with modifications to the control hardware and software.Comment: 10 pages, 6 figures, submitted versio

Why Buckling Stellar Bars Weaken in Disk Galaxies

Young stellar bars in disk galaxies experience a vertical buckling instability which terminates their growth and thickens them, resulting in a characteristic peanut/boxy shape when viewed edge on. Using N-body simulations of galactic disks embedded in live halos, we have analyzed the bar structure throughout this instability and found that the outer third of the bar dissolves completely while the inner part (within the vertical inner Lindblad resonance) becomes less oval. The bar acquires the frequently observed peanut/boxy-shaped isophotes. We also find that the bar buckling is responsible for a mass injection above the plane, which is subsequently trapped by specific 3-D families of periodic orbits of particular shapes explaining the observed isophotes, in line with previous work. Using a 3-D orbit analysis and surfaces of sections, we infer that the outer part of the bar is dissolved by a rapidly widening stochastic region around its corotation radius -- a process related to the bar growth. This leads to a dramatic decrease in the bar size, decrease in the overall bar strength and a mild increase in its pattern speed, but is not expected to lead to a complete bar dissolution. The buckling instability appears primarily responsible for shortening the secular diffusion timescale to a dynamical one when building the boxy isophotes. The sufficiently long timescale of described evolution, ~1 Gyr, can affect the observed bar fraction in local universe and at higher redshifts, both through reduced bar strength and the absence of dust offset lanes in the bar.Comment: 7 pages, 4 figures, ApJ Letters, in pres