Engineering topics

Parker, Scott E.

Publications and source records attributed to Parker, Scott E..

7 records

DOE OSTI2025

Numerical simulations of three-dimensional ion crystal dynamics in a Penning trap using the fast multipole method

We simulate the dynamics, including laser cooling, of three-dimensional (3-D) ion crystals confined in a Penning trap using a newly developed molecular dynamics-like code. The numerical integration of the ions’ equations of motion is accelerated using the fast multipole method to calculate the Coulomb interaction between ions, which allows us to efficiently study large ion crystals with thousands of ions. In particular, we show that the simulation time scales linearly with ion number, rather than with the square of the ion number. By treating the ions’ absorption of photons as a Poisson process, we simulate individual photon scattering events to study laser cooling of 3-D ellipsoidal ion crystals. Initial simulations suggest that these crystals can be efficiently cooled to ultracold temperatures, aided by the mixing of the easily cooled axial motional modes with the low frequency planar modes. In our simulations of a spherical crystal of 1000 ions, the planar kinetic energy is cooled to several millikelvin in a few milliseconds while the axial kinetic energy and total potential energy are cooled even further. This suggests that 3-D ion crystals could be well suited as platforms for future quantum science experiments.

Zaris, John (ORCID:0009000196476323)↗

DOE OSTI2024

Rapid cooling of the in-plane motion of two-dimensional ion crystals in a Penning trap to millikelvin temperatures

Here, we propose a highly feasible technique with no experimental overhead to rapidly cool the in-plane degrees of freedom of large two-dimensional ion crystals in Penning traps. Through simulations, we demonstrate that our approach enables the in-plane modes to cool down to a temperature of around 1 mK in less than 10 ms. Our technique relies on near-resonant coupling of the poorly cooled in-plane motions and the efficiently cooled out-of-plane motions, and is achieved without introducing additional potentials. The rapid cooling enabled by our approach is in contrast to typical operating conditions, where our simulations of the laser cooling dynamics suggest that the ion crystal's in-plane motion cools very slowly on a timescale of several hundreds of milliseconds, a rate likely slower than experimental heating rates. Our work sets the stage for sub-Doppler laser cooling of the planar motion, and more robust and versatile quantum simulation and quantum sensing experiments with two-dimensional crystals in Penning traps.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Correspondence between open bosonic systems and stochastic differential equations

Bosonic mean-field theories can approximate the dynamics of systems of $n$ bosons provided that $n \gg 1$. Here, we show that there can also be an exact correspondence at finite $n$ when the bosonic system is generalized to include interactions with the environment and the mean-field theory is replaced by a stochastic differential equation. When the $n \to \infty$ limit is taken, the stochastic terms in this differential equation vanish, and a mean-field theory is recovered. Besides providing insight into the differences between the behavior of finite quantum systems and their classical limits given by $n \to \infty$, the developed mathematics can provide a basis for quantum algorithms that solve some stochastic nonlinear differential equations. We discuss conditions on the efficiency of these quantum algorithms, with a focus on the possibility for the complexity to be polynomial in the log of the stochastic system size. A particular system with the form of a stochastic discrete nonlinear Schrödinger equation is analyzed in more detail.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Effects of the heat source on the steady-state transport in gradient-driven global gyrokinetic simulations

Heat sources are routinely used in gradient-driven global gyrokinetic simulations to prevent profile relaxation and ensure steady-state turbulence and transport. To have any predictive capability, gradient-driven simulations need to be insensitive to the particular choice of a source. Several forms of the heat source are investigated using gyrokinetic δf particle-in-cell simulations. The ion heat flux is found to be sensitive to the source strength when the source rate is below a fraction (∼1/10) of the linear growth rate.

Physics↗

DOE OSTI2022

Evolution of the marker distribution in gyrokinetic $δf$ particle-in-cell simulations

The evolution of the particle weight in a δf particle-in-cell simulation depends on the marker distribution that can evolve in a turbulent field due to turbulent diffusion. When Monte Carlo methods are used to implement the test particle collision operator, or when the particle motion is not strictly Hamiltonian in a collisionless simulation, the marker distribution will evolve along the particle trajectory and, in general, cannot be known exactly. A two-dimensional numerical marker distribution is proposed as an approximation. It is shown to be advantageous over other common methods for evaluating the marker distribution in long-time turbulence simulations. A generalized two-weight δf-method is proposed to mitigate the marker evolution problem.

, Monte Carlo methods↗

DOE OSTI2021

Equilibration of the planar modes of ultracold two-dimensional ion crystals in a Penning trap

In this work, planar thermal equilibration is studied using direct numerical simulations of ultracold two-dimensional ion crystals in a Penning trap with a rotating wall. The large magnetic field of the trap splits the modes that describe in-plane motion of the ions into two branches: high-frequency cyclotron modes dominated by kinetic energy and low-frequency E × B modes dominated by potential energy associated with thermal position displacements. Using an eigenmode analysis we extract the equilibration rate between these two branches as a function of the ratio of the frequencies that characterize the two branches and observe this equilibration rate to be exponentially suppressed as the ratio increases. Under experimental conditions relevant for current work at NIST, the predicted equilibration time is orders of magnitude longer than any relevant experimental timescales. We also study the coupling rate dependence on the initial energy of the planar modes and the number of ions. In addition, we show how increasing the rotating wall strength improves crystal stability. These details of in-plane mode dynamics help set the stage for developing strategies to efficiently cool the in-plane modes and improve the performance of single plane ion crystals for quantum information processing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Linear embedding of nonlinear dynamical systems and prospects for efficient quantum algorithms

The simulation of large nonlinear dynamical systems, including systems generated by discretization of hyperbolic partial differential equations, can be computationally demanding. Such systems are important in both fluid and kinetic computational plasma physics. This motivates exploring whether a future error-corrected quantum computer could perform these simulations more efficiently than any classical computer. In this work, we describe a method for mapping any finite nonlinear dynamical system to an infinite linear dynamical system (embedding) and detail three specific cases of this method that correspond to previously studied mappings. Then we explore an approach for approximating the resulting infinite linear system with finite linear systems (truncation). Using a number of qubits only logarithmic in the number of variables of the nonlinear system, a quantum computer could simulate truncated systems to approximate output quantities if the nonlinearity is sufficiently weak. Other aspects of the computational efficiency of the three detailed embedding strategies are also discussed.

97 MATHEMATICS AND COMPUTING↗