[Paper Review] Effective coupling parameter for 2D Yukawa liquids and non-invasive measurement of plasma parameters

T. Ott, M. Bönitz|arXiv (Cornell University)|Oct 29, 2010

Atomic and Subatomic Physics Research3 citations

TL;DR

This paper introduces an effective coupling parameter Γ* for 2D Yukawa liquids based on the first-peak height of the pair distribution function g(r), enabling accurate structural comparisons across varying screening lengths. It further proposes a non-invasive reference data method using particle trajectories and velocity autocorrelation functions to extract plasma parameters Γ and κ with ≤15% error, requiring only configuration and velocity snapshots without additional physical inputs.

ABSTRACT

We derive an effective coupling parameter for two-dimensional Yukawa systems based on the height of the first maximum of the pair distribution function. Two variants -- one valid in the high-coupling range, the other for arbitrary couplings of the liquid -- are derived. Comparison to previous approaches to Yukawa coupling parameters shows that the present concept is more general and more accurate. Using, in addition, dynamical information contained in the velocity autocorrelation function, we outline a reference data method that can be employed as a non-invasive measurement scheme of the plasma parameters -- the coupling strength and the screening length. This approach requires only input from a time-series of configuration snapshots and particle velocities with no recourse to additional information about the system. Our results should be directly applicable as a simple, yet reliable diagnostic method for a variety of experiments, including dusty plasmas, colloidal suspensions and ions in traps, and can be employed to facilitate comparisons between experiments, theory and simulations.

Motivation & Objective

• To develop a generalized, accurate effective coupling parameter Γ* for 2D Yukawa liquids that accounts for varying screening lengths κ.
• To overcome limitations of prior coupling parameters that are less accurate or less general, especially at high κ.
• To enable non-invasive, data-driven measurement of both Γ and κ from experimental or simulation data without external physical inputs.
• To provide a practical diagnostic tool applicable to dusty plasmas, colloidal suspensions, and trapped ions using only time-series particle snapshots.
• To validate the method’s accuracy across a broad range of physically relevant parameters (κ ≤ 3, Γ* ≈ 18–124).

Proposed method

• Derive Γ* from the height of the first maximum of the radial pair distribution function g(r), with separate formulations for high- and low-coupling regimes.
• Use the first-peak height g_max as a structural reference to define Γ* such that systems with identical g_max are considered structurally equivalent.
• Incorporate dynamical information from the velocity autocorrelation function (VACF) to extract the first oscillation cycle’s amplitude ΔZ.
• Construct a reference data method (RDM) that maps g_max and ΔZ to Γ and κ using precomputed simulation data, avoiding reliance on temperature, mass, or charge.
• Apply the RDM by first computing g(r) and ΔZ from time-series snapshots, then using empirical fits (Eqs. 16 and 18) to infer κ and Γ from Γ*.
• Validate the method by applying it to molecular dynamics simulations, comparing inferred Γ and κ to original simulation parameters.

Experimental results

Research questions

• RQ1 How can an effective coupling parameter Γ* be defined that enables structural comparison of 2D Yukawa systems across different screening lengths κ?
• RQ2 What is the most accurate and general formulation of Γ* that captures both high- and low-coupling regimes?
• RQ3 Can plasma parameters Γ and κ be extracted from equilibrium particle trajectories alone, without additional physical inputs?
• RQ4 How accurate is the proposed reference data method (RDM) in recovering Γ and κ from simulated or experimental data?
• RQ5 To what extent is the method robust across the experimentally relevant range of κ ≤ 3 and Γ* ≈ 18–124?

Key findings

• The proposed Γ* parameter provides a more general and accurate representation of structural features than previous coupling parameters, especially for κ > 3.
• For high coupling (Γ* ≥ 40), Γ* / Γ = w(κ) holds as a simple one-parameter relation, enabling direct comparison across κ values.
• For low coupling, Γ* is defined as Γ* / Γ = w(κ, Γ), capturing the coupling-dependent screening effects with only three free parameters.
• The reference data method (RDM) achieves an average absolute error of Δκ ≤ 0.15 and relative error ΔΓ/Γ ≤ 0.15 for κ ≤ 3.
• The method is robust across a wide range of parameters, with average errors of Δκ ≈ 0.12 and ΔΓ/Γ ≈ 0.12 for κ from 0 to 5.
• The peak height g_max is sensitive to the liquid-solid transition and can serve as a first-order structural diagnostic.