The quantum distribution functions are of great interest for understanding the properties of dusty plasma. In statistical mechanics, we can get the thermodynamic functions such as the internal energy, the osmotic pressure, and the excess free energy by using the distribution functions for the plasma particles. The importance of quantum distribution function in statistical physics is to give the particles number density in the phase space at time t.
Many authors have calculated the quantum distribution functions. Hussein and Hassan  have calculated the quantum binary distribution of high temperature plasma. Kraeft et al.  used effective potentials to calculate the binary distribution function. Bogoliubov and Kraeft et al.   defined a set of equations describing the dynamics of a system of a large number of interacting particles by using the quantum BBGKY hierarchy. Kaniadakis  showed that the classical BBGKY hierarchical equation, for the n-particle reduced distribution function. Fesciyan  calculated the quantum time independent BBGKY hierarchy of equations at low densities. Hussein et al.  used the BBGKY hierarchy for calculating the binary and triplet distribution functions for one- and two-component Plasmas in terms of Green’s function technique. A quantum mechanical calculation of the radial distribution function for a Plasma was given by  . Also, many others have studied the binary and triplet distribution functions  -  . Many researchers have been interested in identifying the optimal form of the classical distribution functions  . Hansen  calculated the distribution functions for one component plasma in the classical form and quantum corrections.
The calculation of the distribution functions in general form by using curvilinear coordinates makes it easier for the researchers to find the form of the distribution function in the spherical or cylindrical coordinates or any type of coordinates that facilitates the study of the physical problem. Błaszak and Domanski  calculated the canonical quantization of classical mechanics in curvilinear coordinates. Kjaergaard and Mortensen in 1990  made a simple derivation of the quantum mechanical Hamiltonian in curvilinear coordinates. In 2012 a generalized, curvilinear-coordinate formulation of Poisson’s equations to solve for the electrostatic fields in plasma was given by Fichtl et al.  .
Dusty Plasmas plays an important role in experimental physics and in many astrophysical situations  . To understand many space and astrophysical phenomena, and many industrial and physical applications, the study of plasmas containing heavy dust particles is very important  . It differs from ordinary plasma in the presence of dust particles along with a number of positive and negative charges moving at a high speed if compared with the speed of dust  . Wang and Zhang  employed the quantum hydrodynamic model to study the solution and chaotic structures of dust ion-acoustic waves in quantum dusty plasmas consisting of electrons, ions and charged dust particles.
Dusty plasma was originally important in the field of astrophysics. Examples of astronomical dusty plasmas include planetary ring systems (rings of Saturn). The rings of Saturn are the most important models in space dusty plasma study  . The rings of Saturn are made of billions of particles; these particles mostly range from tiny, dust-sized icy grains to giant chunks. They are made of small chunks of ice and rock coated with another material such as dust. A few particles are as large as mountains.
Saturn’s ring system extends up to 175,000 miles (282,000 kilometres) from the planet, yet the vertical height is typically about 30 feet (10 meters) in the main rings  . Named alphabetically in the order they were discovered, the rings are relatively close to each other, with the exception of a gap measuring 2920 miles (4700 kilometres) wide called the Cassini Division that separates Rings A and B. The main rings are A, B and C. Rings D, E, F and G are fainter and more recently discovered; each ring orbits at a different speed around the planet  . But if we ask a question why plasma molecules move in these rings in circular paths and what is the nature of their components. And how to find the distribution functions of their particles. Does the use of curvilinear coordinates make it easier to find and deal with the functions of predisposition? In this study, we attempted to find the form of the distribution functions of the molecules of planetary rings. A model was found that simulates its shape. How fast the particles of planetary rings and whether it is necessary to use Einstein patches for high speeds. In this research, we tried to answer some of these questions.
Particle-in-cell (PIC) simulations are a useful tool in modelling plasma. The electron velocity distribution function and the plasma potential are found by particle-in-cell (PIC) simulations  . Reinmüller in 1998  determined the plasma potential from PIC simulations. We use particle-in-cell (PIC) simulations to find a visualization of dusty three component plasma quantum phase space in curvilinear coordinates.
This work is aimed to calculate the quantum binary and triplet distribution functions of a dusty plasma in curvilinear coordinates. In geometry, curvilinear coordinates are a coordinate system for Euclidean space in which the coordinate lines may be curved. Commonly used curvilinear coordinate systems include rectangular, spherical, and cylindrical coordinate systems. The calculation is based on the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy  .
The cluster expansion method consists of writing the binary distribution function as a power series in the density. The coefficients of different powers of the density involve integrals of the order of that power. These coefficients are then expressed as a sum of a product of integrals. The power series of the binary distribution function in the density convergence badly at high densities, many attempts have been made to overcome such difficulty.
2. The Basic Equations and Hierarchy
Consider 3 Dimensions space with coordinates . A point p in 3d space can be defined using Cartesian coordinates or it can also be defined by its curvilinear coordinates . The relation between the coordinates is then given by the invertible transformation functions:
The surfaces = constant, = constant, = constant are called the coordinate surfaces. The coordinate axes are determined by the tangents to the coordinate curves at the intersection of three surfaces. They are not in general fixed directions in space, which happens to be the case for simple Cartesian coordinates, and thus there is generally no natural global basis for curvilinear coordinates. The momentum operator in quantum mechanics is the gradiant operator . By defining the Jacobi matrix as:
The Jacobian of the transformation is the determinant of the Jacobi matrix
Define the natural basis vectors:
The reduced s-particle density operators defined by Bogoliubov  in the following form
where is the Hamiltonian of our system given by
where V is the potential of the system and Q is the configuration integral is given by
The solution of N-particle schrödinger (time-dependent) equation with this Hamiltonian are given by and from a complete orthonormal basis
Define the N-particle density operator
where wk are positive real probabilities
The density operator follow the Von Neumann equation
In order to derive the quantum BBGKY-hierarchy, we introduce the reduced s-particle density operator as
The equation of motion for the reduced density operator obeys directly the Von Neuman equation. Now by substituting from Equation (12) into Equation (11) we get
where is the potential between particles i,j and is the s-particle Hamiltonian operator. The above equation constitutes the quantum generalization of the (BBGKY) hierarchy.
3. The Binary Distribution Function
We assume that the momentum of the electron lies between and is . Also the momentum of the positron lies between and is and the momentum of the dust (ion) lies between and is . In this section, we shall find the binary distribution function of a quantum dusty in curvilinear coordinates. Firstly, define the quantum N particle distrbution function as follows:
The one particle distribution function is obtained by reducing by integrating over N − 1 positions and momenta then
If there are no external fields
where the value of can be found by normalization:
By substituting from Equation (16) into Equation (17) we get
Then the one particle distribution function is given by
By putting s = 1 into Equation (13) we can get the first equation of quantum BBGKY as
Then we have the binary distribution function by substituting from Equation (19) into Equation (20):
where is the correlation function.
Let us now study the model of dusty three component plasma i.e. the neutral system of point-like particles of positive and negative charges such as electrons, positrons and dust particles like ions. This model is an important model in both laboratory physics and space physics and has many applications  . By finding the quantum distribution functions of this model we can get the important thermodynamic functions such as the internal energy and the equation of state. Dust particles are heavier and slower in their velocities than electrons and positrons where for dusty plasma the relation between electrons, positrons, and dust (ion) density at equilibrium is , where is the number density of sth species and is the charge state of dust  .
Substituting Equations (19) and (21) into (13) for we obtain
4. The Triplet Distribution Function
The quantum triplet distribution function defined in such a way that is the probability of finding a particle of the type 1th in the volume element surrounding ,with momentum in range , a particle of type 2th in the volume element surrounding , with momentum in range and a particle of the type 3th in the volume element surrounding , with momentum in range respectively at time t.
The quantum triplet distribution function is defined by the calculation of the interaction between three charged particles seems rather involved, because the force on particle 1 at time t would depend on the position and momentum of particles 2 and 3 at a retarded time. For simplification the quantum triplet distribution function can be written as
where is the correlation function between particles 1, 2 and 3.
Substituting Equations (19) and (22) into (13) for we obtain
For three component plasma, we can use the two particle correlation function which is given by
is the Debye-Hückel solution and the three particle correlation function which is given by
Whatever particles 1, 2, 3 are, we can write the quantum triplet distribution function as
5. Quantum Distribution Functions of Saturn’s Rings
An alternative statistical description of planetary rings was formulated in a series of papers by Hämeen-Anttila   . It uses a kinetic equation of Boltzmann type for the description of the evolution of the one-particle phase space distribution function, in a similar manner as gas-kinetics. Hämeen-Anttila gives analytical solutions for the collision integrals, where necessary in terms of appropriate approximations. The effect of self-gravity is taken into account in a self consistent manner in the local vertical gravity field of the disk and in its effect on close particle encounters. In principle the theory can treat the average effect of particle surface irregularities stochastically and it is formulated so that it can be extended to describe particle fragmentation and coagulation.
The biggest advantage is that the balance equations for mass, stress, and scale hight of the ring, are given analytically as partial differential equations. Thus, the theory can be applied to investigate the dynamical evolution of a planetary ring. The disk’s self-gravity potential couples to the surface mass density through Poisson’s equation 
where is the surface mass density, a is gravitation parameter and is the Dirac delta function.
where is the natural basis vectors. The Kinetic theory describes the evolution of the local velocity distribution function of an ensemble of particles in terms of the Boltzmann equation or a suitable generalization of it, like Enskog’s theory of hard sphere gases. The kinetic equation can be derived from Liouville’s theorem, appearing as the leading equation in a hierarchy of equations describing n-particle distribution functions in phase space and neglecting correlations between particle pairs. Then the one particle distribution function in Equation (19) in in the local vertical gravity field can be given by
Kinetic theory allows us to incorporate the full complexity of the dynamics of a planetary ring in a statistical description, such as the effects of the motion of ring particles on curved orbits between inelastic collisions, their finite size, the anisotropy of the velocity dispersion, and in principle also coagulation and fragmentation of the ring particles.
The collision motion of an ensemble of identical particles in the plane, in the frame of reference rotating with angular velocity can be described by the Boltzmann equation  in Curvilinear Coordinates:
where was defined by Equation (1) and is the collision integral
which takes into account effects due to the discrete-point nature of the gravitational charges, or collision effects (including diffusion in space and velocity), and defines the change of the distribution function f(r, v, t) arising from ordinary interparticle collisions (in a plasma this term represents the change of f arising from collisions with particles at distances shorter than a Debye length). The Boltzmann form for the collision integral is based on an assumption that the duration of a collision is much less than the time between collisions instantaneous collisions are considered  . The simple Krook integral in the case of a two-dimensional disk of identical particles has the form
where f is the actual distribution function of particles and f0 is the steady-state equilibrium distribution function (Shu and Stewart, 1985)  . The equilibrium axially symmetric distribution function is the Maxwellian with the surface density :
The gravitational potential of a rotating oblate planet
where r is the distance from the center of the planet to the point at which the potential is sought, is the planetocentric latitude of the point, R is the radius of the planet, and is its mass. For saturn , and for more exact values of the parameters see Ref.  .
In this work we obtained the quantum binary and triplet distribution functions of dusty plasma; the calculation is based on curvilinear coordinates and the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. We consider only the thermal equilibrium plasma. The model under consideration is the three-component dusty plasma i.e. neutral system of point like particles of positive and negative charges (electrons and positrons) interspersed with dust particles (ions)  . Figure 1 and Figure 2 show the effect of dust on the quantum binary and triplet distribution functions. From these figures we note that the existence of dust makes a fluctuation in binary and triplet distribution functions in the range of . It was shown that the local stability criterion obtained from the computer models is in general agreement with the theoretical prediction as outlined in the present paper.
The first derivation of reactive quantum Boltzmann equations by Olmstead and Curtiss  starts with the Wigner transformed version of the standard BBGKY hierarchy for a system of stable atomic and diatomic constituents. Also, they calculated the quantum triplet distribution function of moderately dense gases. Alavi et al.  described the time evolution of a one-particle statistical distribution or density operator, influenced by the interaction of another particle which represents the effect of all other particles in the system.
Figure 1. The quantum binary distribution function from Equation (22) without dust effect (solid line) and with a dust effect (dashed line).
Figure 2. The quantum triplet distribution function from Equation (28) without dust effect (solid line) and with a dust effect (dashed line).
Figure 3 and Figure 4 show the comparison between the quantum binary and triplet distribution function based on our result without dust effect (solid line), Alavi et al.  (dashed line) and from Olmstead and Curtiss  (dot-dashed line). We note from this comparison the convergence between the results obtained and the results of the references  and  .
Figure 5 and Figure 6 show the comparison between the quantum binary and triplet distribution function based on our result with the effect of dust, Alavi et al.  and from Olmstead and Curtiss  . After taking into account the effect of dust, we noticed that our results are closer to the results obtained by Olmstead and Curtiss  in the range of , but with the increase in the value of the variable in the range of the results are closer to results obtained by Alavi et al.  .
Initial quantum curvilinear coordinates phase-space for dusty plasma electrons (gray color), positrons (red color) and dust (orange color) was given in Figure 7. Figure 8 shows the quantum curvilinear coordinates phase-space for dusty plasma based on Particle-In-Cell (PIC) simulation methods.
Figure 3. The comparison between the quantum binary distribution function based on our result without dust effect (solid line), Alavi et al.  (dashed line) and from Olmstead and Curtiss  (dot-dashed line).
Figure 4. The comparison between the quantum triplet distribution function based on our result without dust effect (solid line), Alavi et al.  (dashed line) and from Olmstead and Curtiss  (dot-dashed line).
Figure 5. The comparison between the quantum binary distribution function based on our result with a dust effect (solid line), Alavi et al.  (dashed line) and from Olmstead and Curtiss  (dot-dashed line).
Figure 6. The comparison between the quantum triplet distribution function based on our result without dust effect (solid line), Alavi et al.  (dashed line) and from Olmstead and Curtiss  (dot-dashed line).
Figure 7. Initial quantum curvilinear coordinates phase-space for dusty plasma electrons (gray color), positrons (red color) and dust (orange color).
We observe from the phase space that the speed of dust particles is much lower than the speed of electrons and positrons. The density of dust particles varies due to the large volume of dust particles.
Figure 11 and Figure 12 show the quantum binary and triplet distribution functions for dusty plasma in curvilinear phase space . These solutions are important in different branches of physics and other areas of applied sciences and can provide help for researchers to study and understand the physical interpretation of the dusty plasma model.
Figure 8. The quantum curvilinear coordinates phase-space for dusty plasma electrons (gray color), positrons (red color) and dust (orange color) based on Particle-In-Cell (PIC) simulation methods.
Figure 9. The quantum binary distribution function for dusty plasma with .
Figure 10. The quantum triplet distribution function for dusty plasma with .
Figure 11. The quantum binary distribution function for dusty plasma in curvilinear phase-space .
Figure 12. The quantum triplet distribution function for dusty plasma in curvilinear phase-space .
Figure 13. Model of Saturn planet rings in the radius interval (92,000 km; 139,350 km).
Figure 14. 2D Model of Saturn planet rings in the radius interval (92,000 km; 139,350 km).
Figure 15. Ultraviolet ring’s image; red color refer to dust and blue for water ice.
Figure 16. Simulation of dusty plasma for sturan rings for different values of gravitation parameter in the radius interval (92,000 km; 139,350 km).
Model of Saturn planet rings in the radius interval (92,000 km, 139,350 km) was given in Figure 13 and Figure 14. Figure 15 shows ultraviolet ring’s image; red color refers to dust and blue for water ice. Figure 16 shows the Simulation model of dusty plasma for Saturn rings for different values of gravitation parameter. Also, we note from this figure when gravitation parameter increased the dispersion of the presence of particles becomes more visible than the lower values of the gravitation parameter.