Open access peer-reviewed chapter

Numerical Approach to Dynamical Structure Factor of Dusty Plasmas

Written By

Aamir Shahzad, Muhammad Asif Shakoori, Mao-Gang HE and Yan Feng

Submitted: December 5th, 2017 Reviewed: May 4th, 2018 Published: February 27th, 2019

DOI: 10.5772/intechopen.78334

Chapter metrics overview

917 Chapter Downloads

View Full Metrics


The dynamical structure factor [S(k,ω)] gives the information about static and dynamic properties of complex dusty plasma (CDPs). We have used the equilibrium molecular dynamic (EMD) simulations for the investigation of S(k,ω) of strongly coupled CDPs (SCCDPs). In this work, we have computed all possible values of dynamical density with increasing and decreasing sequences of plasma frequency (ωp) and wave number (k) over a wide range of different combinations of the plasma parameters (κ, Γ). Our new simulation results show that the fluctuation of S(k,ω) increases with increasing Г and it decreases with an increase of κ and N. Moreover, investigation shows that the amplitude of S(k,ω) increases by increasing screening (κ) and wave number (k), and it decreases with increasing Г. Our EMD simulation shows that dynamical density of SCCDPs is slightly dependent on N; however, it is nearly independent of other parameters. The presented results obtained through EMD approach are in reasonable agreement with earlier known results based on different numerical methods and plasma states. It is demonstrated that the presented model is the best tool for estimating the density fluctuation in the SCCDPs over a suitable range of parameters.


  • plasma parameter
  • current correlation function
  • dynamical structure factor
  • strongly coupled complex (dusty) plasmas
  • equilibrium molecular dynamics
  • complex dusty plasma

1. Introduction

Dynamical structure factor S(k,ω) is very actively investigated through theoretical, experimental and computer simulation for simple liquids as well as complex systems (dusty plasmas). The S(k,ω) of the dusty plasmas is very significant to understand the dynamic behavior of dust particles in the complex plasmas. The subject of dusty plasma containing micron-size charged condensed particles has recently been actively investigated in the fields of science and technology. In addition, the investigation of dynamical behaviors is also studied in the areas of physics and chemistry of plasmas, ionized gases, and the space environment, environmental sciences, semiconductor plasma processing industries, nuclear energy generation and materials research. Dust in the atmosphere and in the entire universe exists in different shapes and sizes. Mostly dust particles are observed in solid form and sometimes also in liquid and gaseous forms. Current correlations and wave spectra’s in the dusty plasma are generated due to dynamic motion of charged dust particles. The dust particles increase remarkable and unique fundamental physical property of ionized gases and dense plasma. Further dust particle increases future application of dusty plasma in industrial fields including nuclear fusion energy, material modification, and synthesis, environmental remediation, nano, aerospace and medical technologies.

1.1. History of dusty plasma

Current research on the dusty (complex) plasma becomes an interesting field in sciences and technologies. The term “plasma” was first used by the Langmuir in 1924. The most thrilling events in the field of dusty plasmas were occurred in 1980 for planet Saturn mission. Mendis in 1997 observed the bright comet in distant ancestor, which was the excellent cosmic laboratory for the investigation of dusty plasma and their dynamical and physical consequences. The other appearances of dusty plasmas were zodiacal light, the origin nebula, the noctilucent clouds, etc. At the laboratory level, dusty plasma is available in terrestrial laboratory at the remote past. Lyman Spitzer along with Hannes Alfven was recognized that dust in the universe was not merely a hindrance to a visual opinion but it was an essential component of the universe. A dust particle image taken a shape of spokes that rotating around Saturn ring and last surveys designated that these spokes were fine dust material. In 2005, Cassini spacecraft was made a new and improved observation of spokes with a feature that would provide an improved considerate of their source. In 1992, the European spacecraft Ulysses flew by the planet Jupiter and detected the dust particles and measured their masses and impact speed. Again in 1995, NASA spacecraft, Galileo perceived the origin of dust streams around Jupiter. The current enormous interest in complex plasmas started in the mid of 1980s, and started by laboratory investigation of dynamical structure factor of dusty plasmas. Current situation (2000–2017) of dusty plasma is stable in laboratory conditions and it is very significant in industries, science, and technology, medical science, and energy sectors, etc. Different characteristics of dust particles are investigating via theoretical, computer simulation and experimental techniques [1, 2].

1.2. Characteristics and types of dusty plasma

Dusty plasma is the plasma that contains an addition of dust charged particles along with electrons, ions, and neutral particles. The dynamical properties of plasma become more complex when we insert dust particle, so it is also known as complex plasma [3]. Dust particles are much larger in charge as compared to electrons, ions, and neutral particles and its size vary from hundreds of millimeter to 10 of nanometers and having a mass approximately 3 × 10−11 kg. The dynamic behaviors are easily observed experimentally by CCD camera. Dusty plasmas exist in space as well as in the laboratory and these dust particles are negatively charged but sometimes positive charge as well. It is charged through photoionization, electron bombardment, etc. the amount of charge on dust particles depends upon shape and size of dust particles. It has mostly spherical shape but sometimes also having rod type shape and irregular [4]. Dusty plasmas are classified on the basis of density, temperature, potential and thermal energies. The Coulomb coupling parameter describes the classification of complex dusty plasmas. The Coulomb coupling that is defined as “the ratio of average potential energy to the average thermal energy” and mathematically it can be expressed as: Γ=P.EK.E=Q24πε01kBT[5], where kB is the Boltzmann constant, Tis system temperature and ε0 is permittivity of free space.

1.2.1. Weakly coupled dusty (complex) plasmas (WCCDPs)

When dust charged particles have average thermal energy due to neighboring particles much larger than the average potential energy, then that plasma is known as weakly coupled complex dusty plasma (WCDPs). The WCCDPs have a high temperature but low density and value of Coulomb coupling parameter less than 1 (mathematically Γ < 1). In case of WCCDPs, any structural form does not exist. The background of WCCDPs is considered as ionized gases.

1.2.2. Strongly coupled dusty (complex) plasma (SCCDPs)

The SCCDPs is speedily emerging filed from last three decades. It is very significant to astrophysical plasma and quickly progressing in laboratory experiments. In this type of dusty plasma average thermal energy of charged particles due to neighboring particles are much smaller than average potential energy and mathematically Γ ≥ 1. In case of SCCDPs, it has high density and low temperature and can be specified in structural farm [6, 7, 8].

1.3. Dusty plasma in atmosphere and laboratory

There are many systems in the atmosphere where dust particles are established. Spaces between the stars are filled with a large amount of dust and gases. Dust particles in the interstellar region are metallic i.e. graphite, magnetite, and amorphous of carbon, dielectric material i.e. silicates and ices, etc. Comets, planetary space, planetary ring, and earth atmosphere are the region of our solar system. Gossamer ring, halo ring, and main rig are the three systems of Jupiter’s ring. In the Saturn rings systems are mostly ices and its size vary from meter to micron. A Uranian rings system has major rings such as 6, 5, 4, 3, α, β, η, γ. In Neptune’s ring systems appear in curios twisted materials and structure is dirty ice and composition such as iron, nickel, sulfur, earth atmosphere dusty ice, etc. [9]. A simple device for producing dusty plasmas is a dusty plasma device which is a single-ended Q-machine modified to allow the dispersal of dust grains. Dusty plasmas are produced by suspending micron-sized dust particles in a stratum of a dc neon glow discharge. Dusty plasma has been for the first time confined in cylindrical symmetric radio frequency plasma (RF) system also in the semiconductor industry.

1.4. Role of dust particles and applications

Dust particles have charge and chemically active species, it is formed and growth in dusty plasma devices. Sometimes in the form of a mixture of gases such as SiH4, silane, oxygen, O2, and Ar, etc. Secondly dust particles are formed in devices when atoms and molecules are spurted from walls and electrodes into the plasma by electron and ion bombardment. Moreover, the growth of dust particles in the plasma is coagulation, nucleation and surface growth. Thermal fluctuations and Coulomb interaction play the significant role for determination structure of SCCDPs. When the values of Г (>1) increase then system organized from nonideal gases phase to ordered condensed phase. Dust particles are suspended in the gaseous plasma phase with few electrons temperature and charge up to 104 ordered. Interestingly, the dust clouds in a dusty plasma formed into the structural form even at room temperature. Dust particle has large mass as compared to ion and an electron which gives the results slowly downtime scale and it can easily observed the macroscopic structure and its dynamical behavior directly study in space and real time [10]. Dusty plasma used for nanocrystalline silicon particles grow in the silane plasmas used to increase efficiency and lifetime of the silicon solar cells. It is used for thin film coating applied in plasma-enhanced chemical vapor deposition (PECVD) for the improvement of material surface properties. Carbon-based nanostructure growth in the hydrocarbon plasmas or fluorocarbons used to produce thin carbon films. It is used to improve material properties such as chemical inertness, high hardness and wear resistance. Self-lubrication coating and wear resistance using different compound as a dust particle (MoS2). Ar/CH4 plasma used for making the nanocrystalline diamond films fabricated which as exceptional properties such as chemical inertness, high hardness, and extreme smoothness which used to improve the performance of cutting tools. Diamond whiskers fabricated by the etching in RF plasmas for the enhancement of electron field emission. The reactive ion etching (RIE) process are used to a precise efficiently sharpen micro-tips of diamond [11]. Complex (dusty) plasmas (CDPs) have various advantages in a different industries, technologies, and energy sector due to the existence of dust particles. CDPs are stable under the laboratory condition. The CDPs can also be used for the diagnostic purpose because dust particles are trapped at the room temperature and keep their desire dynamical state for hours. Dusty plasma frequency in the range several hertz and easily observed through CCD cameras. It is produced in the gas discharge tube with natural gas pressure range that varies from 1 to 100 pa, which is subject to moderate damping [8]. Moreover, magnetized dustyplasma device used to produce a number of the verity of magnetic fields configuration with the help of four independent superconducting coils [12]. Magnetized glow discharge dusty plasma device, RF plasma device, ISS experiment and DC glow discharge devices used for different applications in industries and diagnostic purpose of dusty plasma. Dust particles are found in a tokamak (fusion plasma) and dusty plasmas depositions techniques devices [11, 13].

1.5. Dynamical structure factor

The dynamical structure factor S(k,ω)] gives the information about static and dynamic properties of the fluid in simple and complex systems. In hydrodynamic condition, the S(k,ω) provides experimental calculable quantities such as the thermal diffusivity, adiabatic sound velocity and the ratio of specific heats [8]. These properties of the fluids are measured through light scattering, x-ray and inelastic neutron experiments on a substance such as dense plasmas, liquids, and glasses. Sound waves are generated through S(k,ω) in strongly coupled CDPs (SCCDPs) and density is more strongly damped at liquid phase [14, 15]. The SCCDP is many body dynamical systems that show different collective excitations and their properties investigated through numerical simulation and theoretical approaches [16]. In condensed matter physics, the pair correlation function usually cannot be determined directly. Rather, the structure factor is determined by scattering of x-rays or neutrons. In dusty plasmas, we are able to measure the pair correlation function directly and to calculate the structure factor in order to compare with condensed matter experiments. The S(k,ω) is just the Fourier transform of the pair correlation function. Fluctuation of the dynamical density of dusty plasma generates current correlation spectra such as longitudinal and transverse currents [17]. Dynamic ion structure factor of warm dense matter and dense plasma consist of the complete information of ions in strongly interacting systems and also influenced by the electrons property. It is closely associated with density fluctuation, thus determines transport properties and many relaxations such as electrons ion temperature equilibrium and stopping power and also the equation of state. It is also used for diagnostics of the extreme states of matter like a warm dense matter of x-ray Thomson scattering [18]. Dynamical scattering function is given through times correlation function, fluctuated density of liquid argon and light scattering function [19].


2. Numerical model and simulation techniques

In this section, we have implemented molecular dynamic (MD) simulation code with Ewald summation for forces and energies which makes it possible to account the long-range Coulomb interparticle interactions. We trace the motion of single charge species and integrated through leapfrog method and assume that the presence of neutralizing homogenous background. In this plasma environment, random fluctuating forces and friction forces are acting on a charged particle in addition to which forces initiating from the interaction of charged particle. Length of simulation cubic box is defined as 4πN31/3 and particles have a random spatial configuration for the beginning of simulation [15, 20]. Fluctuation of microscopic density is observed for different plasma parameters approaching near the equilibrium state [21]. The presented study includes the solution of the equation of motion of a system and particle interacts with each other through Yukawa potential. Provided that an accurate potential can be established for the system of attention under study and equilibrium MD (EMD) can be used irrespective of the phase condition and thermodynamic of the system involved. Yukawa potential is most commonly used potential (screened Coulomb) for SCCDPs including many physical systems such as physics of chemical and polymer, medicine and biology systems, astrophysics, environmental, etc. Major advantage for using this potential is that it reduces the calculation time compared other potentials [22]. The interaction potential energy of a charged particle in Yukawa liquid is given


Here Qis the charge on dust particle, ris the distance between interacting particles, λD is Debye screening length that accounts for the screening of interaction of other plasma species and εo is permittivity of free space. The scaling (dimensionless) parameters, which fully characterized the system, one is known as Coulomb coupling parameter [23],


where, aws is the “Wigner-Seitz” radius and it is defined as πn12with nis the equilibrium dust number density, kBis the Boltzmann constant and Tis absolute system temperature. It is noted that the Γ is measured as the ratio of average potential energy to average kinetic energy per particle. Second scaling parameter is Debye screening parameter and it is given as κ = awsD.

The EMD simulations are performed for a particle number that is chosen between 500 and 1000 particles in a microcanonical ensemble using periodic boundary conditions and minimum image convention of the dust particles. It is to be mentioned that the number density (n) is defined as n = N/V,here Nis the number of particles and Vis the volume of simulation box and it is calculated as V = 4πN/3. On our case, most of simulations are performed for N = 500 and it is observed that the mentioned number of particles is suitable for EMD computations with statistical uncertainty limits. In presented case, the simulation time step is dt = 0.0001 and total simulation time limit was 425,000 step units. The EMD simulations are run between 4.25 106 (1/ωp) to 3.25 106 (1/ωp) time units for each combination of (Γ, κ) in the series of recording dynamical structure factor (DSF), S(k,ω), of SCCDPs. It can be seen that the first patch of S(k,ω) results is obtained after the time limit of (38000) step unit. For our case, 13 patches of S(k,ω) results are obtained and results show that each patch has nearly behavior of S(k,ω) showing the accuracy of numerical algorithm. In this work, the dynamical structure factor computations of SCCDPs are reported for a wide range of Coulomb coupling parameters Γ ≡ (1, 200) and the Debye screening strength κ ≡ (1.4, 4).

2.1. Model of DSF [S(k,ω)]

The number density of single species also known as mass density and dimensionless quantity in molecular dynamic units can be written as


here ris the point at time t. from the practical point of view, it is known as local density. The average occupancy in the small volume of rspace and it calculates over the short interval of time. The definition of density is satisfied matter conserved requirement.


The space and time-dependent density correlation are explain from the van Hove correlation function which is define as G(r, t) = Gs(r, t) + Gs(r, t) and further detail is given in Ref. [23]. Fourier transform of density fluctuation and new expression is given as


where kis the wave number and becomes equal to k = 2π/Land Lis the length of simulation box. The intermediate function can be defined as


The relation between static structure factors with the dynamic structure factor is


This satisfies the sum rule as


The dynamical structure function is related to longitudinal current correlation function, which is expressed as [22].


3. Results and discussion

The EMD simulation has been used for the calculation of S(k,ω) to understand the dynamic phenomenon of particles for 3D SCCDPs. We have analyzed our simulation results of S(k,ω) in term of frequency, amplitude and fluctuation rate with respect to these parameters (κ, Γ, and k).

This section shows an overview of our results obtained through EMD simulations for dynamical structure factor S(k,ω) function of SCCDPs at κ = 1, 2, 3 and 5 with N = 500. Figures 14 shown our results for S(k,ω) over a wide suitable range of plasma states (Г, κ) at four values of wave vectors k = (0, 1, 2 and 3). A sequence of dynamical structure factor in increasing of wave vector (k) is computed to determine the suitable equilibrium values of S(k,ω). We have performed 16 EMD simulation with N = 500 for each screening parameters at different four values of k. There are 64 simulations are carried out for different combination of plasma parameters in order to observe complete behaviors of DSF S(k,ω) at higher varying frequency (ωp) as compared to earlier simulation results [16]. It is observed that the presented results obtained through EMD simulations for varying parameters have suitable signal-to-noise ratio of the DSF. Our results are satisfactory good agreement with earlier EMD estimations and show that the presented EMD results at higher ωp and earlier EMD simulations have comparable performance with small system size, both yielding the close values of the DSF S(k,ω). Moreover, the system temperature (1/ Г), strength of Debye screening (κ), system run time (total time), system size (N), and wave vector (k) are changed to observe how efficiently the presented EMD method computes the DSF S(k,ω) of SCCDPs.

Figure 1.

Variation of dynamical structure factorS(k,ω) as a function of plasma frequency (ω) of strongly coupled complex plasma atκ = 1.4,N = 500 and waves numberk = 0, 1, 2, and 3 for (a) Г = 1, (b) Г = 50, (c) Г = 100, (d) Г = 200.

Figure 2.

Variation of dynamical structure factorS(k,ω) as a function of plasma frequency (ω) of strongly coupled complex plasma atκ = 2,N = 500 and waves numberk = 0, 1, 2, and 3 for (a) Г = 1, (b) Г = 50, (c) Г = 100, (d) Г = 200.

Figure 3.

Variation of dynamical structure factorS(k,ω) as a function of plasma frequency (ω) of strongly coupled complex plasma atκ = 3,N = 500 and waves numberk = 0, 1, 2, and 3 for (a) Г = 1, (b) Г = 50, (c) Г = 100, (d) Г = 200.

Figure 4.

Variation of dynamical structure factorS(k,ω) as a function of plasma frequency (ω) of strongly coupled complex plasma atκ = 4,N = 500 and waves numberk = 0, 1, 2, and 3 for (a)Γ= 1, (b) Г = 50, (c) Г = 100, (d) Г = 200.

In addition, in each panels of Figures 14, we have shown the behaviors of S(k, ω) at four values of k. All these data are excellent due to the good statistics and allow the analysis of the structure–function within a broad dynamical range. It is examined that the peaks of S(k,ω) are decay at lower Г values for different four values of κand k. It can be noted that the peaks survive up to higher kvalues at intermediate to higher Г values and shift toward the higher frequency (ωp). It is further observed that the wave spectra of S(k,ω) shifts toward sinusoidal waves form, at the intermediate values of Г = 50 and wave spectra shifts toward square type waves forms at higher values of (Г = 100,200). Moreover, the panels (a) to (d) of each figures represent the results of S(k,ω) in nonideal gases state to liquid and crystalline order state corresponding to Г (1, 200). In case of crystalline order state, the amplitude of vibrating particle decreases significantly and it converts completely in square waveform. It is interesting to note that the fundamental behavior of dust particles is different and it shows decaying trends in nonideal gases state, sinusoidal form in liquid state and square wave form in crystalline state at fixed screening value. Figures show that the results of S(k,ω) depend on the plasma parameters (κ, Г), as expected. Furthermore, it is investigated that in the ideal form of dusty plasma with high temperatures (low coupling values) the dust particles are exponentially decaying from high to low amplitude. One of justification is the dust particles transferred their energy to the surrounding particles that at high plasma temperature values. Furthermore, in panels (b) represent the EMD results of S(k,ω) in the liquefied state of dusty plasma. In this regime, dust particles comparatively less transfer their energy to the system. At this regime amplitude is maximum for simulation box size (L = max) and frequency of oscillating of dust particles is low and in the sinusoidal waveform. The frequency and amplitude of oscillating dust particles increase with increasing wave numbers k = 2, 3 and exhibit in the periodic wave form. Panels (c) show the results of S(k,ω) of dusty plasma in the nearly crystalline states. We have analyzed that in this regime the particles are tightly bound. In this case, the dust particles have small amplitude as compared to the gaseous and liquefy forms of dusty plasma.

It is observed from each panel of Figures 14, the dynamic of dust particles increases with increasing wave number. It is to be noted that the values of amplitude of the S(k,ω) are 3 = 0.7406, 0.6662, 0.8004 and 0.7902 respectively, for four wave numbers (k = 0, 1, 2, and 3) at the same values of κ, Г and N. It is observed that the dynamical structure factor of SCCDPs depends on plasma parameters (κ, Г). The frequency mode of S(k,ω) is high at low κfor SCCDPs. It is observed that κis equally affecting on the dynamic of dust particles either the dusty plasma in any phase (gaseous, liquid and crystalline).


4. Summary

We have employed EMD simulations for the investigations of S(k,ω) over a wide range of Coulomb and Debye screening parameters (κ, Г). It has been shown that the presented EMD technique and previous numerical methods have equivalent performance for the wide range of plasma state points, both yielding satisfactory data for plasma S(k,ω). New simulations provide more consistent and inclusive results for the plasma S(k,ω) over a complete range of Г (1, 200) and κ(1.4, 4) than the previously known numerical results. Our investigations show that the dynamic of dust particles are exponentially decayed, sinusoidal form and crystalline form, respectively, in the gaseous, liquefy regime and crystalline states of SCCDPs. Moreover, the dynamical spectra of dusty plasma do not observe at very high values of Г and low values of κ. Moreover, our results indicate that the dynamical structure factor in SCCDPs depends on κ, Г and k.It has been shown that the presented simulation has comparable performance with the earlier simulation of S(k,ω) over the wide domain of plasma states. For future work, it is suggested that presented EMD technique based Ewald summation can be used to investigate and explore dynamical structure factor behaviors in other simple liquid, dipolar and ionic materials.



The authors are grateful to the National Advanced Computing, National Center of Physics (NCP), Pakistan, and National High Performance Computing Center of X’ian Jiaotong. The authors thank the University, P.R. China, for allocating computer time to test and run our MD code.


Abbreviation and symbols

EMDequilibrium molecular dynamics
MDmolecular dynamics
SCCDPsstrongly coupled complex dusty plasmas
WCCDPsweakly coupled complex dusty plasmas
CDPscomplex (dusty) plasma
DSFdynamical structure factor
PECVDplasma-enhanced chemical vapor deposition
S(k,ω)dynamical structure factor
kwave number
Nnumber of particles
ωpplasma frequency
κDebye screening strength
ГCoulomb coupling


  1. 1. Shukla PK, Mendis DA, Desai T. Conference Proceeding. Advances in Dusty Plasmas. Singapore: World Scientific; 1997. QB791.153671996
  2. 2. Gregori G. Molecular dynamics simulations of the equilibrium dynamics of non-ideal plasmas. (Doctoral dissertation). University of Oxford;2012
  3. 3. Merlino RL. Dusty plasmas and applications in space and industry. Iowa City USA. Plasma Physics Applied. 2006;81:73-110. ISBN: 81-7895-230-0
  4. 4. Shahzad A, Shakoori MA, He M-G, Bashir S. Sound waves in complex (Dusty) plasmas. In: Computational and Experimental Studies of Acoustic Waves. Rijeka, Croatia: InTech; 2018.
  5. 5. Shahzad A, H M-G, Shakoori MA. Thermal transport and non-newtonian behaviors of 3D complex liquids using molecular simulations. In: Applied Sciences and Technology (IBCAST), 2017 14th International Bhurban Conference on (pp. 472-474). USA: IEEE; 2017 January. DOI: 10.1109/IBCAST.2017.7868096
  6. 6. Bellan PM. Fundamentals of Plasma Physics. UK: Cambridge University Press; 2008. ISBN: 9780521528009
  7. 7. Shahzad A, He MG. Thermal conductivity of three-dimensional Yukawa liquids (dusty plasmas). Contributions to Plasma Physics. 2012;52(8):667. DOI: 10.1002/ctpp.201200002
  8. 8. Killian T, Pattard T, Pohl T, Rost J. Ultracold neutral plasmas. Physics Reports. 2007;449:77. DOI: 10.1016/j.physrep.2007.04.007
  9. 9. Shukla PK, Mamun AA. Introduction to Dusty Plasma Physics. New York: CRC Press; 2015. ISBN 0 7503 0653
  10. 10. Bonitz M, Henning C, Block D. Complex plasmas: A laboratory for strong correlations. Reports on Progress in Physics. 2010;73(6):066501. DOI: 10.1088/0034-4885/73/6/066501
  11. 11. Schoen M, Vogelsang R, Hoheisel C. Computation and analysis of the dynamic structure factor S (k, ω) for small wave vectors: A molecular dynamics study for a Lennard-Jones liquid country. Molecular Physics. 1986;57(3):445-471.
  12. 12. Thomas E, Konopka U, Artis D, Lynch B. The magnetized dusty plasma experiment. Cambridge University Press; 2015. DOI:
  13. 13. Alley WE, Alder BJ, Yip S. The neutron scattering function for hard spheres. Physical Review A. 1983;27(6):3174.
  14. 14. Magyar P, Hartmann P, Kalman GJ, Golden KI, Donkó Z. Factorization of 3-point static structure functions in 3D Yukawa liquids. Contributions to Plasma Physics. 2016;56(9):816-829. DOI: 10.1002/ctpp.201500062
  15. 15. Arkhipov YV, Askaruly A, Davletov AE, Dubovtsev DY, Donkó Z, Hartmann P, Korolov I, Conde L, Tkachenko IM. Direct determination of dynamic properties of coulomb and Yukawa classical one-component plasmas. Physical Review Letters. 2017;119(4). DOI: 045001,
  16. 16. Donkó Z, Hartmann P, Kalman GJ, Golden KI. Simulations of strongly coupled charged particle systems: Static and dynamical properties of classical bilayers. Journal of Physics A: Mathematical and General. 2003;36(22):5877. DOI:
  17. 17. Vorberger J, Donko Z, Tkachenko IM, Gericke DO. Dynamic ion structure factor of warm dense matter. Physical Review Letters. 2012;109(22):225001. DOI:
  18. 18. Enciso E, Almarza NG, del Prado V, Bermejo FJ, Zapata EL, Ujaldón M. Molecular-dynamics simulation on simple fluids: Departure from linearized hydrodynamic behavior of the dynamical structure factor. Physical Review E. 1994;50(2):1336. DOI:
  19. 19. Donkó Z, Hartmann P, Kalman GJ. Molecular dynamics simulations of strongly coupled plasmas: Localization and microscopic dynamics. Physics of Plasmas. 2003;10(5):1563-1568
  20. 20. Piel A, Nosenko V, Goree J. Experiments and molecular-dynamics simulation of elastic waves in a plasma crystal radiated from a small dipole source. Physical Review Letters. 2002;89(8):085004. DOI:
  21. 21. Shahzad A, He M-G. Numerical experiment of thermal conductivity in two-dimensional Yukawa liquids. Physics of Plasmas. 2015;22(12):123707
  22. 22. Rapaport DC, Blumberg RL, McKay SR, Christian W. The art of molecular dynamics simulation. Computers in Physics. 1996;10(5):456-456
  23. 23. Shahzad A, Aslam A, He M-G. Equilibrium molecular dynamics simulation of shear viscosity of two-dimensional complex (dusty) plasmas. Radiation Effects and Defects in Solids;169(11):931-941

Written By

Aamir Shahzad, Muhammad Asif Shakoori, Mao-Gang HE and Yan Feng

Submitted: December 5th, 2017 Reviewed: May 4th, 2018 Published: February 27th, 2019