A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas
Identifieur interne : 001361 ( Istex/Corpus ); précédent : 001360; suivant : 001362A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas
Auteurs : A L Alexandrov ; I V Schweigert ; F M PeetersSource :
- New Journal of Physics [ 1367-2630 ] ; 2008.
Abstract
Nanoparticle charging in a capacitively coupled radio frequency discharge in argon is studied using a particle in cell Monte Carlo collisions method. The plasma parameters and dust potential were calculated self-consistently for different unmovable dust profiles. A new method for definition of the dust floating potential is proposed, based on the information about electron and ion energy distribution functions, obtained during the kinetic simulations. This approach provides an accurate balance of the electron and ion currents on the dust particle surface and allows us to precisely calculate the dust floating potential. A comparison of the obtained floating potentials with the results of the traditional orbital motion limit (OML) theory shows that in the presence of the ion resonant charge exchange collisions, even when the OML approximation is valid, its results are correct only in the region of a weak electric field, where the ion drift velocity is much smaller than the thermal one. With increasing ion drift velocity, the absolute value of the calculated dust potential becomes significantly smaller than the theory predicts. This is explained by a non-Maxwellian shape of the ion energy distribution function for the case of fast ion drift.
Url:
DOI: 10.1088/1367-2630/10/9/093025
Links to Exploration step
ISTEX:B74E140A61C686AC8CB111B0FC21BA43C88A2403Le document en format XML
<record><TEI wicri:istexFullTextTei="biblStruct"><teiHeader><fileDesc><titleStmt><title xml:lang="en">A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas</title><author><name sortKey="Alexandrov, A L" sort="Alexandrov, A L" uniqKey="Alexandrov A" first="A L" last="Alexandrov">A L Alexandrov</name><affiliation><mods:affiliation>Institute of Theoretical and Applied Mechanics, Russian Academy of Sciences, Novosibirsk 630090, Russia</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: aalex@itam.nsc.ru</mods:affiliation></affiliation></author><author><name sortKey="Schweigert, I V" sort="Schweigert, I V" uniqKey="Schweigert I" first="I V" last="Schweigert">I V Schweigert</name><affiliation><mods:affiliation>Institute of Theoretical and Applied Mechanics, Russian Academy of Sciences, Novosibirsk 630090, Russia</mods:affiliation></affiliation><affiliation><mods:affiliation>Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: ischweig@yahoo.com</mods:affiliation></affiliation></author><author><name sortKey="Peeters, F M" sort="Peeters, F M" uniqKey="Peeters F" first="F M" last="Peeters">F M Peeters</name><affiliation><mods:affiliation>Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</mods:affiliation></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">ISTEX</idno><idno type="RBID">ISTEX:B74E140A61C686AC8CB111B0FC21BA43C88A2403</idno><date when="2008" year="2008">2008</date><idno type="doi">10.1088/1367-2630/10/9/093025</idno><idno type="url">https://api.istex.fr/document/B74E140A61C686AC8CB111B0FC21BA43C88A2403/fulltext/pdf</idno><idno type="wicri:Area/Istex/Corpus">001361</idno></publicationStmt><sourceDesc><biblStruct><analytic><title level="a" type="main" xml:lang="en">A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas</title><author><name sortKey="Alexandrov, A L" sort="Alexandrov, A L" uniqKey="Alexandrov A" first="A L" last="Alexandrov">A L Alexandrov</name><affiliation><mods:affiliation>Institute of Theoretical and Applied Mechanics, Russian Academy of Sciences, Novosibirsk 630090, Russia</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: aalex@itam.nsc.ru</mods:affiliation></affiliation></author><author><name sortKey="Schweigert, I V" sort="Schweigert, I V" uniqKey="Schweigert I" first="I V" last="Schweigert">I V Schweigert</name><affiliation><mods:affiliation>Institute of Theoretical and Applied Mechanics, Russian Academy of Sciences, Novosibirsk 630090, Russia</mods:affiliation></affiliation><affiliation><mods:affiliation>Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</mods:affiliation></affiliation><affiliation><mods:affiliation>E-mail: ischweig@yahoo.com</mods:affiliation></affiliation></author><author><name sortKey="Peeters, F M" sort="Peeters, F M" uniqKey="Peeters F" first="F M" last="Peeters">F M Peeters</name><affiliation><mods:affiliation>Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</mods:affiliation></affiliation></author></analytic><monogr></monogr><series><title level="j">New Journal of Physics</title><title level="j" type="abbrev">New J. Phys.</title><idno type="ISSN">1367-2630</idno><idno type="eISSN">1367-2630</idno><imprint><publisher>Institute of Physics Publishing</publisher><date type="published" when="2008">2008</date><biblScope unit="volume">10</biblScope><biblScope unit="issue">9</biblScope><biblScope unit="page" from="1">1</biblScope><biblScope unit="page" to="12">12</biblScope></imprint><idno type="ISSN">1367-2630</idno></series><idno type="istex">B74E140A61C686AC8CB111B0FC21BA43C88A2403</idno><idno type="DOI">10.1088/1367-2630/10/9/093025</idno><idno type="articleID">286109</idno><idno type="articleNumber">093025</idno></biblStruct></sourceDesc><seriesStmt><idno type="ISSN">1367-2630</idno></seriesStmt></fileDesc><profileDesc><textClass></textClass><langUsage><language ident="en">en</language></langUsage></profileDesc></teiHeader><front><div type="abstract">Nanoparticle charging in a capacitively coupled radio frequency discharge in argon is studied using a particle in cell Monte Carlo collisions method. The plasma parameters and dust potential were calculated self-consistently for different unmovable dust profiles. A new method for definition of the dust floating potential is proposed, based on the information about electron and ion energy distribution functions, obtained during the kinetic simulations. This approach provides an accurate balance of the electron and ion currents on the dust particle surface and allows us to precisely calculate the dust floating potential. A comparison of the obtained floating potentials with the results of the traditional orbital motion limit (OML) theory shows that in the presence of the ion resonant charge exchange collisions, even when the OML approximation is valid, its results are correct only in the region of a weak electric field, where the ion drift velocity is much smaller than the thermal one. With increasing ion drift velocity, the absolute value of the calculated dust potential becomes significantly smaller than the theory predicts. This is explained by a non-Maxwellian shape of the ion energy distribution function for the case of fast ion drift.</div></front></TEI><istex><corpusName>iop</corpusName><author><json:item><name>A L Alexandrov</name><affiliations><json:string>Institute of Theoretical and Applied Mechanics, Russian Academy of Sciences, Novosibirsk 630090, Russia</json:string><json:string>E-mail: aalex@itam.nsc.ru</json:string></affiliations></json:item><json:item><name>I V Schweigert</name><affiliations><json:string>Institute of Theoretical and Applied Mechanics, Russian Academy of Sciences, Novosibirsk 630090, Russia</json:string><json:string>Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</json:string><json:string>E-mail: ischweig@yahoo.com</json:string></affiliations></json:item><json:item><name>F M Peeters</name><affiliations><json:string>Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</json:string></affiliations></json:item></author><language><json:string>eng</json:string></language><originalGenre><json:string>paper</json:string></originalGenre><abstract>Nanoparticle charging in a capacitively coupled radio frequency discharge in argon is studied using a particle in cell Monte Carlo collisions method. The plasma parameters and dust potential were calculated self-consistently for different unmovable dust profiles. A new method for definition of the dust floating potential is proposed, based on the information about electron and ion energy distribution functions, obtained during the kinetic simulations. This approach provides an accurate balance of the electron and ion currents on the dust particle surface and allows us to precisely calculate the dust floating potential. A comparison of the obtained floating potentials with the results of the traditional orbital motion limit (OML) theory shows that in the presence of the ion resonant charge exchange collisions, even when the OML approximation is valid, its results are correct only in the region of a weak electric field, where the ion drift velocity is much smaller than the thermal one. With increasing ion drift velocity, the absolute value of the calculated dust potential becomes significantly smaller than the theory predicts. This is explained by a non-Maxwellian shape of the ion energy distribution function for the case of fast ion drift.</abstract><qualityIndicators><score>7.278</score><pdfVersion>1.4</pdfVersion><pdfPageSize>595.276 x 841.89 pts (A4)</pdfPageSize><refBibsNative>true</refBibsNative><keywordCount>0</keywordCount><abstractCharCount>1258</abstractCharCount><pdfWordCount>4462</pdfWordCount><pdfCharCount>23561</pdfCharCount><pdfPageCount>12</pdfPageCount><abstractWordCount>193</abstractWordCount></qualityIndicators><title>A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas</title><genre><json:string>research-article</json:string></genre><host><volume>10</volume><publisherId><json:string>nj</json:string></publisherId><pages><total>12</total><last>12</last><first>1</first></pages><issn><json:string>1367-2630</json:string></issn><issue>9</issue><genre><json:string>journal</json:string></genre><language><json:string>unknown</json:string></language><eissn><json:string>1367-2630</json:string></eissn><title>New Journal of Physics</title></host><publicationDate>2008</publicationDate><copyrightDate>2008</copyrightDate><doi><json:string>10.1088/1367-2630/10/9/093025</json:string></doi><id>B74E140A61C686AC8CB111B0FC21BA43C88A2403</id><score>0.23106353</score><fulltext><json:item><original>true</original><mimetype>application/pdf</mimetype><extension>pdf</extension><uri>https://api.istex.fr/document/B74E140A61C686AC8CB111B0FC21BA43C88A2403/fulltext/pdf</uri></json:item><json:item><original>false</original><mimetype>application/zip</mimetype><extension>zip</extension><uri>https://api.istex.fr/document/B74E140A61C686AC8CB111B0FC21BA43C88A2403/fulltext/zip</uri></json:item><istex:fulltextTEI uri="https://api.istex.fr/document/B74E140A61C686AC8CB111B0FC21BA43C88A2403/fulltext/tei"><teiHeader><fileDesc><titleStmt><title level="a" type="main" xml:lang="en">A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas</title></titleStmt><publicationStmt><authority>ISTEX</authority><publisher>Institute of Physics Publishing</publisher><availability><p>IOP Publishing and Deutsche Physikalische Gesellschaft</p></availability><date>2008</date></publicationStmt><sourceDesc><biblStruct type="inbook"><analytic><title level="a" type="main" xml:lang="en">A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas</title><author xml:id="author-1"><persName><forename type="first">A L</forename><surname>Alexandrov</surname></persName><email>aalex@itam.nsc.ru</email><affiliation>Institute of Theoretical and Applied Mechanics, Russian Academy of Sciences, Novosibirsk 630090, Russia</affiliation></author><author xml:id="author-2"><persName><forename type="first">I V</forename><surname>Schweigert</surname></persName><email>ischweig@yahoo.com</email><affiliation>Institute of Theoretical and Applied Mechanics, Russian Academy of Sciences, Novosibirsk 630090, Russia</affiliation><affiliation>Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</affiliation></author><author xml:id="author-3"><persName><forename type="first">F M</forename><surname>Peeters</surname></persName><affiliation>Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</affiliation></author></analytic><monogr><title level="j">New Journal of Physics</title><title level="j" type="abbrev">New J. Phys.</title><idno type="pISSN">1367-2630</idno><idno type="eISSN">1367-2630</idno><imprint><publisher>Institute of Physics Publishing</publisher><date type="published" when="2008"></date><biblScope unit="volume">10</biblScope><biblScope unit="issue">9</biblScope><biblScope unit="page" from="1">1</biblScope><biblScope unit="page" to="12">12</biblScope></imprint></monogr><idno type="istex">B74E140A61C686AC8CB111B0FC21BA43C88A2403</idno><idno type="DOI">10.1088/1367-2630/10/9/093025</idno><idno type="articleID">286109</idno><idno type="articleNumber">093025</idno></biblStruct></sourceDesc></fileDesc><profileDesc><creation><date>2008</date></creation><langUsage><language ident="en">en</language></langUsage><abstract><p>Nanoparticle charging in a capacitively coupled radio frequency discharge in argon is studied using a particle in cell Monte Carlo collisions method. The plasma parameters and dust potential were calculated self-consistently for different unmovable dust profiles. A new method for definition of the dust floating potential is proposed, based on the information about electron and ion energy distribution functions, obtained during the kinetic simulations. This approach provides an accurate balance of the electron and ion currents on the dust particle surface and allows us to precisely calculate the dust floating potential. A comparison of the obtained floating potentials with the results of the traditional orbital motion limit (OML) theory shows that in the presence of the ion resonant charge exchange collisions, even when the OML approximation is valid, its results are correct only in the region of a weak electric field, where the ion drift velocity is much smaller than the thermal one. With increasing ion drift velocity, the absolute value of the calculated dust potential becomes significantly smaller than the theory predicts. This is explained by a non-Maxwellian shape of the ion energy distribution function for the case of fast ion drift.</p></abstract></profileDesc><revisionDesc><change when="2008">Published</change></revisionDesc></teiHeader></istex:fulltextTEI><json:item><original>false</original><mimetype>text/plain</mimetype><extension>txt</extension><uri>https://api.istex.fr/document/B74E140A61C686AC8CB111B0FC21BA43C88A2403/fulltext/txt</uri></json:item></fulltext><metadata><istex:metadataXml wicri:clean="corpus iop not found" wicri:toSee="no header"><istex:xmlDeclaration>version="1.0" encoding="ISO-8859-1"</istex:xmlDeclaration><istex:docType SYSTEM="http://ej.iop.org/dtd/iopv1_5_2.dtd" name="istex:docType"></istex:docType><istex:document><article artid="nj286109"><article-metadata><jnl-data jnlid="nj"><jnl-fullname>New Journal of Physics</jnl-fullname><jnl-abbreviation>New J. Phys.</jnl-abbreviation><jnl-shortname>NJP</jnl-shortname><jnl-issn>1367-2630</jnl-issn><jnl-coden>NJOPFM</jnl-coden><jnl-imprint>Institute of Physics Publishing</jnl-imprint><jnl-web-address>stacks.iop.org/NJP</jnl-web-address></jnl-data><volume-data><year-publication>2008</year-publication><volume-number>10</volume-number></volume-data><issue-data><issue-number>9</issue-number><coverdate>September 2008</coverdate></issue-data><article-data><article-type type="paper" sort="regular"></article-type><type-number type="paper" numbering="article" artnum="093025">093025</type-number><article-number>286109</article-number><first-page>1</first-page><last-page>12</last-page><length>12</length><pii></pii><doi>10.1088/1367-2630/10/9/093025</doi><copyright>IOP Publishing and Deutsche Physikalische Gesellschaft</copyright><ccc>1367-2630/08/093025+12$30.00</ccc><printed></printed></article-data></article-metadata><header><title-group><title>A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas</title><short-title>A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas</short-title><ej-title>A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas</ej-title></title-group><author-group><author address="nj286109ad1" email="nj286109ea1"><first-names>A L</first-names><second-name>Alexandrov</second-name></author><author address="nj286109ad1" second-address="nj286109ad2" email="nj286109ea2"><first-names>I V</first-names><second-name>Schweigert</second-name></author><author address="nj286109ad2"><first-names>F M</first-names><second-name>Peeters</second-name></author><short-author-list>A L Alexandrov <italic>et al</italic></short-author-list></author-group><address-group><address id="nj286109ad1" showid="yes">Institute of Theoretical and Applied Mechanics, <orgname>Russian Academy of Sciences</orgname>, Novosibirsk 630090, <country>Russia</country></address><address id="nj286109ad2" showid="yes">Departement Fysica, <orgname>Universiteit Antwerpen</orgname>, Groenenborgerlaan 171, B-2020 Antwerpen, <country>Belgium</country></address><e-address id="nj286109ea1"><email mailto="a_alex@itam.nsc.ru">a_alex@itam.nsc.ru</email></e-address><e-address id="nj286109ea2"><email mailto="ischweig@yahoo.com">ischweig@yahoo.com</email></e-address></address-group><history received="10 July 2008" online="23 September 2008"></history><abstract-group><abstract><heading>Abstract</heading><p indent="no">Nanoparticle charging in a capacitively coupled radio frequency discharge in argon is studied using a particle in cell Monte Carlo collisions method. The plasma parameters and dust potential were calculated self-consistently for different unmovable dust profiles. A new method for definition of the dust floating potential is proposed, based on the information about electron and ion energy distribution functions, obtained during the kinetic simulations. This approach provides an accurate balance of the electron and ion currents on the dust particle surface and allows us to precisely calculate the dust floating potential. A comparison of the obtained floating potentials with the results of the traditional orbital motion limit (OML) theory shows that in the presence of the ion resonant charge exchange collisions, even when the OML approximation is valid, its results are correct only in the region of a weak electric field, where the ion drift velocity is much smaller than the thermal one. With increasing ion drift velocity, the absolute value of the calculated dust potential becomes significantly smaller than the theory predicts. This is explained by a non-Maxwellian shape of the ion energy distribution function for the case of fast ion drift.</p></abstract></abstract-group></header><body refstyle="numeric"><sec-level1 id="nj286109s1" label="1"><heading>Introduction</heading><p indent="no"> The charging of dust particles in a discharge plasma is crucial in the description of dust–plasma interaction. Usually the floating potential of the dust is calculated from the balance of electron and ion currents on the dust particle surface using a Maxwellian distribution function for electrons and ions (see, for example [<cite linkend="nj286109bib1">1</cite>]). However, it is known from experiments that in a capacitively coupled radio frequency (ccrf) discharge without dust the electron energy distribution function (EEDF) can be far from Maxwellian even in a bulk plasma. The EEDFs were measured in [<cite linkend="nj286109bib2">2</cite>] in a ccrf discharge in helium and argon to illustrate the difference in the electron kinetics in non-Ramsauer and Ramsaur gases. It was shown that in helium the shape of the EEDF remains Maxwellian for all pressures. In argon the EEDF shape changes from Druyvesteyn to bi-Maxwellian with decreasing gas pressure. The dust particle charge was investigated in [<cite linkend="nj286109bib3">3</cite>] in a low pressure argon discharge plasma using a Boltzmann equation for the electrons and fluid equations for the ions. It was found that the presence of the dust particles modifies the shape of the EEDF from Druyvesteyn to Maxwellian due to the fact that the dust grains collect fast electrons.</p><p>Unfortunately the modification of the ion energy distribution function (IEDF) with the presence of dust is not completely understood. The shifted Maxwellian distribution is the commonly used assumption for the IEDF in a non-zero external electrical field. Particle in cell Monte Carlo collision (PIC-MCC) simulations of charging of unmovable dust were performed in [<cite linkend="nj286109bib4">4</cite>]. However, the shapes of the EEDF and IEDF, as well as the procedure of equilibration of the electron and ion currents on the dust particle surface were not discussed. In [<cite linkend="nj286109bib4">4</cite>] the unmovable dust was taken to be uniformly distributed over the discharge gap including the electrode sheaths and this assumption makes the results less reliable, because the heavy negatively charged dust is not able to penetrate in the electrode sheath area. Recently the formation of the steady state nanosize dust profile in the rf discharge and the influence of dust size on the plasma parameters were studied in argon in [<cite linkend="nj286109bib5">5</cite>] and in a mixture of argon and acetylene in [<cite linkend="nj286109bib6">6</cite>].</p><p>In this paper, we study the charging of the nanoparticles in a ccrf argon discharge by using the EEDF and IEDF from PIC-MCC simulations. The discharge parameters with the static dust distribution are found from a self-consistent PIC-MCC simulation. Special attention is paid to compare the PIC-MCC EEDF and IEDF with the Maxwellian ones.</p><p>This paper is organized as follows. In section <secref linkend="nj286109s2">2</secref>, we present the model of a discharge plasma with movable nanoparticles. Section <secref linkend="nj286109s3">3</secref> is devoted to the description of our new method of calculation of the dust particle charge. In section <secref linkend="nj286109s4">4</secref>, we present the results of our PIC-MCC calculations and compare them with the analytic model results. Our conclusions are given in section <secref linkend="nj286109s5">5</secref>.</p></sec-level1><sec-level1 id="nj286109s2" label="2"><heading>Physical model</heading><p indent="no"> The PIC-MCC model of a high frequency discharge plasma with nanoparticles is based on the electron and ion kinetic equations, including their interaction with dust and the Poisson equation. In this model, the velocity distribution functions of electrons, <inline-eqn></inline-eqn>, and ions, <inline-eqn></inline-eqn> (three-dimensional over velocity and one-dimensional in space) are obtained by solving the Boltzmann equations numerically,
<display-eqn id="nj286109eqn1" textype="equation" notation="LaTeX" eqnnum="1"></display-eqn><display-eqn id="nj286109eqn2" textype="equation" notation="LaTeX" eqnnum="2"></display-eqn>where <inline-eqn><math-text><italic>v</italic><sub><upright>e</upright></sub></math-text></inline-eqn>, <inline-eqn><math-text><italic>v</italic><sub><upright>i</upright></sub></math-text></inline-eqn>, <inline-eqn><math-text><italic>n</italic><sub><upright>e</upright></sub></math-text></inline-eqn>, <inline-eqn><math-text><italic>n</italic><sub><upright>i</upright></sub></math-text></inline-eqn>, <inline-eqn><math-text><italic>m</italic><sub><upright>e</upright></sub></math-text></inline-eqn> and <inline-eqn><math-text><italic>m</italic><sub><upright>i</upright></sub></math-text></inline-eqn> are the electron and ion velocities, densities and masses, respectively, <inline-eqn></inline-eqn> is the electric field, <inline-eqn><math-text><italic>J</italic><sub><upright>e</upright></sub></math-text></inline-eqn>, <inline-eqn><math-text><italic>J</italic><sub><upright>i</upright></sub></math-text></inline-eqn> are the collisional integrals for electrons and ions which include elastic and inelastic collisions with background gas and negatively charged dust particles.</p><p>The electron and ion kinetics in the PIC-MCC model are simulated with test particles with weight <inline-eqn><math-text><italic>W</italic></math-text></inline-eqn>, which means the number of electrons/ions represented by this test particle. Below we will call the test particles electrons and ions. The equations of motion are integrated for electrons and ions in the discharge electric field accounting for the collisions with neutral atoms and with negatively charged dust. The probabilities of various events (elastic scattering, excitation of neutrals, ionization and absorption on the dust surface) during an electron/ion time step are set by the cross sections depending on the test particle's energy and simulated by using the null-collision method [<cite linkend="nj286109bib7">7</cite>] and a random number generator. If an ionization event takes place, a new electron and ion with the same coordinate and weight are introduced. Electrons/ions leaving the calculation area through the electrodes or captured by dust are excluded from the simulation. The interelectrode gap is covered with a calculation grid (one-dimensional in this case). The macroscopic characteristics of the discharge such as the plasma densities, ion and electron average energy, the ion/electron flux are calculated from the EEDF and IEDF over the calculation grid. The electrical field distribution is calculated at each time step from the Poisson equation:
<display-eqn id="nj286109eqn3" textype="equation" notation="LaTeX" eqnnum="3"></display-eqn>where <inline-eqn><math-text><italic>n</italic><sub><upright>e</upright></sub></math-text></inline-eqn> and <inline-eqn><math-text><italic>n</italic><sub><upright>i</upright></sub></math-text></inline-eqn> are obtained from the PIC-MCC calculation. The boundary conditions are the sinusoidal applied voltage on one electrode and zero voltage on the second grounded electrode. In this model, the space charge of the dust particles with the density <inline-eqn><math-text><italic>n</italic><sub><upright>d</upright></sub>(<italic>x</italic>)</math-text></inline-eqn> and the mean particle charge <inline-eqn><math-text><italic>Z</italic><sub><upright>d</upright></sub>(<italic>x</italic>)</math-text></inline-eqn> is included. The particle charge is found self-consistently during the simulation, as described below. The calculation starts from chosen initial distributions for the electrons, ions and dust corresponding to some analytical approximation for the discharge plasma. The iterative procedure of solving equations (<eqnref linkend="nj286109eqn1">1</eqnref>)–(<eqnref linkend="nj286109eqn3">3</eqnref>) proceeds until a steady state solution is achieved. Usually this takes about 1000 simulated rf cycles for discharge with unmovable dust. The number of particles representing the electrons and ions in our PIC-MCC simulations is usually <inline-eqn><math-text>10<sup>4</sup></math-text></inline-eqn>–<inline-eqn><math-text>10<sup>5</sup></math-text></inline-eqn>. The interelectrode distance is divided into a calculation grid with 101–201 nodes for definition of the macroscopic plasma parameters. The time step is <inline-eqn><math-text>τ<sub><upright>e</upright></sub>=(0.3−1)×10<sup>−11</sup> <upright>s</upright></math-text></inline-eqn> for electrons and <inline-eqn><math-text>τ<sub><upright>i</upright></sub>=(1−3)×10<sup>−10</sup> <upright>s</upright></math-text></inline-eqn> for ions.</p><p>For electron scattering in argon the cross sections are taken from [<cite linkend="nj286109bib8">8</cite>]. For ions we take into account the resonant charge exchange collisions. Additionally our model incorporates the cross sections of electron and ion absorption by negatively charged dust particles, taken from the orbital motion limit (OML) theory [<cite linkend="nj286109bib9">9</cite>]:
<display-eqn textype="equation" notation="LaTeX" number="no"></display-eqn><display-eqn id="nj286109eqn4" textype="equation" notation="LaTeX" eqnnum="4"></display-eqn><display-eqn textype="equation" notation="LaTeX" number="no"></display-eqn>here <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn> is the floating potential of the dust particle surface gained in the discharge plasma, <inline-eqn><math-text>ϵ<sub><upright>e</upright>, <upright>i</upright></sub></math-text></inline-eqn> is the kinetic energy of the electron/ion, <inline-eqn><math-text><italic>r</italic><sub><upright>d</upright></sub></math-text></inline-eqn> is the radius of the dust particle, <inline-eqn><math-text><italic>e</italic></math-text></inline-eqn> is the elementary charge. Since the calculations were performed for low gas pressure and for nanometre-size particles, the conditions of validity of the OML approach <inline-eqn><math-text><italic>l</italic><sub><upright>e</upright>, <upright>i</upright></sub>≫λ<sub><upright>D</upright></sub>≫<italic>r</italic><sub><upright>d</upright></sub></math-text></inline-eqn> (where <inline-eqn><math-text><italic>l</italic><sub><upright>e</upright>, <upright>i</upright></sub></math-text></inline-eqn> is the mean free paths for electrons and ions, <inline-eqn><math-text>λ<sub><upright>D</upright></sub></math-text></inline-eqn> is the Debye screening length) are fulfilled here. Note that our model is not limited by the OML approach, and it is possible to use any other absorption cross sections in analytical or tabular form. In this work, we restrict our study to static dust and concentrate on the determination of the dust floating potential.</p></sec-level1><sec-level1 id="nj286109s3" label="3"><heading>New method for calculation of the dust particle charge</heading><p indent="no">In the steady state of discharge glow, the floating potential of the dust particle <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn> can be derived from the equality of the electron and ion currents onto the dust particle surface. The volume densities of both currents <inline-eqn><math-text><italic>j</italic><sub><upright>e</upright></sub></math-text></inline-eqn> and <inline-eqn><math-text><italic>j</italic><sub><upright>i</upright></sub></math-text></inline-eqn> onto the dust surface (in other words, the charge carrier absorption rate per unit volume) at the position <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn> and time <inline-eqn><math-text><italic>t</italic></math-text></inline-eqn> are given by the following basic equations (see, for example [<cite linkend="nj286109bib10">10</cite>]):
<display-eqn id="nj286109eqn5" textype="equation" notation="LaTeX" eqnnum="5"></display-eqn><display-eqn id="nj286109eqn6" textype="equation" notation="LaTeX" eqnnum="6"></display-eqn>where <inline-eqn><math-text>ϵ</math-text></inline-eqn> is the kinetic energy, <inline-eqn><math-text><italic>V</italic><sub><upright>e, i</upright></sub></math-text></inline-eqn> are the magnitudes of electron and ion velocity vectors, respectively. Averaging of equations (<eqnref linkend="nj286109eqn5">5</eqnref>) and (<eqnref linkend="nj286109eqn6">6</eqnref>) over a discharge cycle leads to an integral equation for <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn> at each position <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>:
<display-eqn textype="equation" notation="LaTeX" number="no"></display-eqn>In the traditional OML approach, using the cross sections (<secref linkend="nj286109s4">4</secref>) and assuming the Maxwellian velocity distribution for electrons and ions, equations (<eqnref linkend="nj286109eqn5">5</eqnref>) and (<eqnref linkend="nj286109eqn6">6</eqnref>) give the analytical expressions for the currents onto a single dust particle [<cite linkend="nj286109bib10">10</cite>]:
<display-eqn id="nj286109eqn7" textype="equation" notation="LaTeX" eqnnum="7"></display-eqn>for electron current, and, taking into account that ions may have drift velocity <inline-eqn><math-text><italic>u</italic></math-text></inline-eqn> of the same order of magnitude as the thermal one, hence they are assumed to have shifted on <inline-eqn><math-text><italic>u</italic></math-text></inline-eqn> Maxwellian distribution:
<display-eqn id="nj286109eqn8" textype="equation" notation="LaTeX" eqnnum="8"></display-eqn>for ion current.</p><p>Here <inline-eqn><math-text><italic>n</italic><sub><upright>e, i</upright></sub></math-text></inline-eqn> and <inline-eqn><math-text><italic>T</italic><sub><upright>e, i</upright></sub></math-text></inline-eqn> are the electron/ion density and temperature, <inline-eqn><math-text><italic>k</italic></math-text></inline-eqn> is the Boltzmann constant, <inline-eqn></inline-eqn> is the ratio of ion drift and thermal velocities. Note also that expression (<eqnref linkend="nj286109eqn8">8</eqnref>) has two well-known limit cases:
<display-eqn textype="equation" notation="LaTeX" number="no"></display-eqn>for the case <inline-eqn><math-text><italic>g</italic>≪1</math-text></inline-eqn>, and
<display-eqn textype="equation" notation="LaTeX" number="no"></display-eqn>for the reverse case.</p><p>The above shown expressions are widely used in OML theory. If the values of <inline-eqn><math-text><italic>n</italic><sub><upright>e, i</upright></sub>(<italic>x</italic>)</math-text></inline-eqn> and <inline-eqn><math-text><italic>T</italic><sub><upright>e, i</upright></sub>(<italic>x</italic>)</math-text></inline-eqn> are known in every position <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>, the floating potential <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn> may be calculated from equations (<eqnref linkend="nj286109eqn7">7</eqnref>) and (<eqnref linkend="nj286109eqn8">8</eqnref>) if <inline-eqn><math-text><italic>I</italic><sub><upright>e</upright></sub>=<italic>I</italic><sub><upright>i</upright></sub></math-text></inline-eqn>. To a first approximation this allows us to complete the model of rf discharge with nanoparticles, assuming the dust particle charge <inline-eqn><math-text><italic>Z</italic><sub><upright>d</upright></sub>(<italic>x</italic>)=<italic>r</italic><sub><upright>d</upright></sub>ϕ<sub><upright>d</upright></sub>(<italic>x</italic>)</math-text></inline-eqn>, which is correct for the particles with size much smaller than the Debye screening length.</p><p>However, since it is known that the energy distribution in a discharge plasma may be far from Maxwellian, the aim of the developed algorithm is to avoid the assumption about the Maxwellian shape of the energy distribution. The PIC-MCC simulation of the complex plasma provides two other possible ways for obtaining of the dust particle charge. The simplest way for the direct definition of the currents is to count the number of events of electron and ion absorption by the dust particles with the cross sections (<secref linkend="nj286109s4">4</secref>) during the test particles motion simulation, with adjustment of the floating potential to minimize the difference in the electron and ion currents. The problem of such an approach is the statistical noise of the currents, that requires averaging of their magnitudes over large numbers of rf cycles. We offer a more accurate and less time consuming method, namely to calculate the currents by the equations (<eqnref linkend="nj286109eqn5">5</eqnref>) and (<eqnref linkend="nj286109eqn6">6</eqnref>), using the EEDF and IEDF <inline-eqn><math-text><italic>f</italic><sub><upright>e</upright></sub>(<italic>t</italic>,<italic>x</italic>, ϵ)</math-text></inline-eqn> and <inline-eqn><math-text><italic>f</italic><sub><upright>i</upright></sub>(<italic>t</italic>,<italic>x</italic>, ϵ)</math-text></inline-eqn>, obtained from the PIC-MCC simulation of the discharge plasma with dust.</p><p>In the PIC-MCC simulation, the integration over energy and the averaging over a discharge cycle in equations (<eqnref linkend="nj286109eqn5">5</eqnref>) and (<eqnref linkend="nj286109eqn6">6</eqnref>) is equivalent to a summation over the test particles ensemble:
<display-eqn id="nj286109eqn9" textype="equation" notation="LaTeX" eqnnum="9"></display-eqn>where <inline-eqn><math-text><italic>k</italic></math-text></inline-eqn> is a summation over electrons, and
<display-eqn id="nj286109eqn10" textype="equation" notation="LaTeX" eqnnum="10"></display-eqn>where the summation is over the ions. The above expressions are the direct analogue of equations (<eqnref linkend="nj286109eqn5">5</eqnref>) and (<eqnref linkend="nj286109eqn6">6</eqnref>) for the case of a discrete ensemble of test particles. Here, <inline-eqn><math-text><italic>W</italic><sub><italic>k</italic>, <italic>l</italic></sub></math-text></inline-eqn> is the weight of the <inline-eqn><math-text><italic>k</italic></math-text></inline-eqn>th electron/<inline-eqn><math-text><italic>l</italic></math-text></inline-eqn>th ion, <inline-eqn><math-text><italic>V</italic><sub><italic>k</italic>, <italic>l</italic></sub></math-text></inline-eqn> and <inline-eqn><math-text>ϵ<sub><italic>k</italic>, <italic>l</italic></sub></math-text></inline-eqn> are its velocity magnitude and kinetic energy, respectively. The summation in equations (<eqnref linkend="nj286109eqn9">9</eqnref>) and (<eqnref linkend="nj286109eqn10">10</eqnref>) is carried out during one rf cycle and only over the particles corresponding to the calculation grid node with coordinate <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn> (the terms <inline-eqn><math-text><italic>W</italic><sub><italic>k</italic>, <italic>l</italic></sub><italic>V</italic><sub><italic>k</italic>, <italic>l</italic></sub><italic>n</italic><sub><upright>d</upright></sub>(<italic>x</italic>)σ<sub><upright>id</upright>, <upright>ed</upright></sub></math-text></inline-eqn> are taken at each grid node using the standard PIC procedure). For the self-consistent definition of the floating potential <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn>, we apply the iterative method, in which after each rf cycle we add a correction <inline-eqn><math-text>△ϕ</math-text></inline-eqn> to the value of <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn>. To accelerate the convergence of the iterations, we take into account that the dependence of <inline-eqn><math-text><italic>j</italic><sub><upright>e</upright></sub></math-text></inline-eqn> on <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn> is exponential, whereas for <inline-eqn><math-text><italic>j</italic><sub><upright>i</upright></sub></math-text></inline-eqn> it is linear. The correction is calculated at each iteration step in all nodes of the calculation grid from the following equation:
<display-eqn id="nj286109eqn11" textype="equation" notation="LaTeX" eqnnum="11"></display-eqn>where <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn> is the value from the previous iteration and <inline-eqn><math-text><italic>T</italic><sub><upright>e</upright></sub></math-text></inline-eqn> is the mean electron temperature taken from the PIC-MCC simulation (here it is used as the iteration parameter). Equation (<eqnref linkend="nj286109eqn11">11</eqnref>) is solved by the Newton method. Then the summation in equations (<eqnref linkend="nj286109eqn9">9</eqnref>) and (<eqnref linkend="nj286109eqn10">10</eqnref>) is performed again at the next rf cycle with the corrected value for the potential <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub>+△ϕ</math-text></inline-eqn> and equation (<eqnref linkend="nj286109eqn11">11</eqnref>) is solved again. After 3–4 iterations, the currents <inline-eqn><math-text><italic>j</italic><sub><upright>e</upright></sub><sup><upright>PIC</upright></sup></math-text></inline-eqn> and <inline-eqn><math-text><italic>j</italic><sub><upright>i</upright></sub><sup><upright>PIC</upright></sup></math-text></inline-eqn> become practically equal for all <inline-eqn><math-text><italic>x</italic></math-text></inline-eqn>.</p><p>The proposed method for the calculation of the dust floating potential <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn> is restricted neither by the OML assumptions (the cross sections <inline-eqn><math-text>σ<sub><upright>ed</upright></sub></math-text></inline-eqn> and <inline-eqn><math-text>σ<sub><upright>id</upright></sub></math-text></inline-eqn> may be given in another form) nor by the assumption about the Maxwellian distribution for the electrons and ions. Our method provides an accurate dust potential determination during the rf discharge simulation and can be used to test the validity of the traditional OML approach (equations (<eqnref linkend="nj286109eqn7">7</eqnref>) and (<eqnref linkend="nj286109eqn8">8</eqnref>)). It should be mentioned that the additional computational time required for the summation in equations (<eqnref linkend="nj286109eqn9">9</eqnref>) and (<eqnref linkend="nj286109eqn10">10</eqnref>) and solving equation (<eqnref linkend="nj286109eqn11">11</eqnref>) is very small in relation to the computation time needed for the solution of equations (<eqnref linkend="nj286109eqn1">1</eqnref>–<eqnref linkend="nj286109eqn3">3</eqnref>).</p></sec-level1><sec-level1 id="nj286109s4" label="4"><heading>Results and discussion</heading><p indent="no">The most typical simulation results are shown for the case of ccrf discharge in argon at 50 mTorr pressure, with the interelectrode distance <inline-eqn><math-text><italic>d</italic>=8 <upright>cm</upright></math-text></inline-eqn>, the voltage frequency 13.6 MHz and the voltage amplitude 120 V. The dust particle radius was chosen as <inline-eqn><math-text><italic>r</italic><sub><upright>d</upright></sub>=100 <upright>nm</upright></math-text></inline-eqn>, the dust density distribution was considered as constant. Two different dust density profiles were taken in our simulation: (a) one with a maximum in the center of the discharge gap, <inline-eqn><math-text><italic>n</italic><sub><upright>d</upright></sub>(<italic>x</italic>)=<italic>n</italic><sub><upright>max</upright></sub><upright>sin</upright>(π(<italic>x</italic>−<italic>l</italic><sub><upright>s</upright></sub>)/(<italic>d</italic>−2<italic>l</italic><sub><upright>s</upright></sub>))</math-text></inline-eqn>, and (b) with two maxima near the electrode sheaths and a void in the gap center, which is often observed in experiments, <inline-eqn><math-text><italic>n</italic><sub><upright>d</upright></sub>(<italic>x</italic>)=<italic>n</italic><sub><upright>max</upright></sub><upright>sin</upright><sup>2</sup>(2π(<italic>x</italic>−<italic>l</italic><sub><upright>s</upright></sub>)/(<italic>d</italic>−2<italic>l</italic><sub><upright>s</upright></sub>))</math-text></inline-eqn>. The maximum of the density <inline-eqn><math-text><italic>n</italic><sub><upright>max</upright></sub>=10<sup>6</sup> <upright>cm</upright><sup>−3</sup></math-text></inline-eqn>, <inline-eqn><math-text><italic>l</italic><sub><upright>s</upright></sub></math-text></inline-eqn> is the width of electrode sheath, where dust is absent (<inline-eqn><math-text><italic>n</italic><sub><upright>d</upright></sub>(<italic>x</italic>)=0</math-text></inline-eqn> for <inline-eqn><math-text><italic>x</italic><<italic>l</italic><sub><upright>s</upright></sub></math-text></inline-eqn> and <inline-eqn><math-text><italic>x</italic>><italic>d</italic>−<italic>l</italic><sub><upright>s</upright></sub></math-text></inline-eqn>). The electrode sheath <inline-eqn><math-text><italic>l</italic><sub><upright>s</upright></sub>=1 <upright>cm</upright></math-text></inline-eqn> means the maximal sheath width in the cathodic phase of the discharge rf cycle. Together with a self-consistent determination of <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub>(<italic>x</italic>)</math-text></inline-eqn> by equations (<eqnref linkend="nj286109eqn9">9</eqnref>)–(<eqnref linkend="nj286109eqn11">11</eqnref>), the traditional analytical estimation by the OML theory, obtained by solving the system of equations (<eqnref linkend="nj286109eqn7">7</eqnref>) and (<eqnref linkend="nj286109eqn8">8</eqnref>) for <inline-eqn><math-text><italic>I</italic><sub><upright>e</upright></sub>=<italic>I</italic><sub><upright>i</upright></sub></math-text></inline-eqn>, and called below <inline-eqn><math-text>ϕ<sub><upright>OML</upright></sub></math-text></inline-eqn>, was also made for comparison.</p><p>Figure <figref linkend="nj286109fig1">1</figref> shows the densities of the charge carriers in the discharge plasma for both types of the dust density profile. It can be seen that the ion density is about equal to the sum of the densities of the electrons and the dust charge everywhere in the gap, except the electrode sheaths. This means that the central region of the discharge plasma remains quasineutral, and the negative charge of the particles is effectively screened by the ions. The shape of the electric potential of the plasma is similar for both dust profiles. The average electron energy in the gap center shown in figure <figref linkend="nj286109fig2">2</figref> is 4.5 eV, and gradually decreases to 2.5 eV in the electrode sheaths. The ion energy in the quasineutral region is close to the thermal one of the neutral gas (<inline-eqn><math-text><italic>kT</italic><sub><upright>i</upright></sub>=0.026 <upright>eV</upright></math-text></inline-eqn>), but significantly increases toward the sheaths (see figure <figref linkend="nj286109fig2">2</figref>). The mean free path for ions <inline-eqn><math-text><italic>l</italic><sub><upright>i</upright></sub></math-text></inline-eqn> under these conditions is near 1 mm, the Debye screening length <inline-eqn><math-text>λ<sub><upright>D</upright></sub></math-text></inline-eqn> in the region containing dust is 0.02–0.06 mm.</p><figure id="nj286109fig1" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="indented"><graphic-file version="print" format="EPS" width="22pc" printcolour="no" filename="images/nj286109fig1.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj286109fig1.jpg"></graphic-file></graphic><caption type="figure" id="nj286109fc1" label="Figure 1"><p indent="no">Density of electrons <inline-eqn><math-text><italic>n</italic><sub><upright>e</upright></sub>(<italic>x</italic>)</math-text></inline-eqn>, ions <inline-eqn><math-text><italic>n</italic><sub><upright>i</upright></sub>(<italic>x</italic>)</math-text></inline-eqn> and dust charge <inline-eqn><math-text><italic>Z</italic><sub><upright>d</upright></sub>(<italic>x</italic>)<italic>n</italic><sub><upright>d</upright></sub>(<italic>x</italic>)</math-text></inline-eqn> in the discharge gap for dust localized in the gap center (a) and near the electrode sheaths (b).</p></caption></figure><figure id="nj286109fig2" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="indented"><graphic-file version="print" format="EPS" width="21.6pc" printcolour="no" filename="images/nj286109fig2.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj286109fig2.jpg"></graphic-file></graphic><caption type="figure" id="nj286109fc2" label="Figure 2"><p indent="no">Average electron energy <inline-eqn><math-text>ϵ<sub><upright>e</upright></sub></math-text></inline-eqn> (solid curve), ion energy <inline-eqn><math-text>ϵ<sub><upright>i</upright></sub></math-text></inline-eqn> (dashed curve) and ion drift velocity <inline-eqn><math-text><italic>u</italic></math-text></inline-eqn> (crosses) in the discharge gap.</p></caption></figure><p>Figure <figref linkend="nj286109fig3">3</figref> shows the volume density of the electron and ion currents onto the dust surface <inline-eqn><math-text><italic>j</italic><sub><upright>e</upright></sub></math-text></inline-eqn> and <inline-eqn><math-text><italic>j</italic><sub><upright>i</upright></sub></math-text></inline-eqn>. These values are the time-averaged amounts of the charge carriers captured by dust in <inline-eqn><math-text>1 <upright>cm</upright><sup>3</sup></math-text></inline-eqn> per 1 s obtained by the direct count of MCCs (as described in the previous section) and from the equations (<eqnref linkend="nj286109eqn9">9</eqnref>) and (<eqnref linkend="nj286109eqn10">10</eqnref>). In the first case, the events of absorption with cross section (<secref linkend="nj286109s4">4</secref>) were counted during the simulation of the test particles motion for 300 rf cycles. Since these events are relatively seldom, the statistical noise is large despite the long averaging time. In the second case, the sums (<eqnref linkend="nj286109eqn9">9</eqnref>) and (<eqnref linkend="nj286109eqn10">10</eqnref>) are plotted, averaged for 10 rf cycles. The values of both sums practically coincide for the self-consistent determination of <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn>. It is seen that the direct count of absorbtion events during the PIC-MCC simulation and equations (<eqnref linkend="nj286109eqn9">9</eqnref>) and (<eqnref linkend="nj286109eqn10">10</eqnref>) give equivalent current values, but the statistical noise for the second case is strongly suppressed.</p><figure id="nj286109fig3" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="indented"><graphic-file version="print" format="EPS" width="21.9pc" printcolour="no" filename="images/nj286109fig3.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj286109fig3.jpg"></graphic-file></graphic><caption type="figure" id="nj286109fc3" label="Figure 3"><p indent="no">Volume density of electron and ion currents onto the dust surface: <inline-eqn><math-text><italic>j</italic><sub><upright>e</upright>, <upright>i</upright></sub><sup><upright>MCC</upright></sup></math-text></inline-eqn>—direct count of absorptions during Monte Carlo simulation; <inline-eqn><math-text><italic>j</italic><sub><upright>e</upright>, <upright>i</upright></sub><sup><upright>PIC</upright></sup></math-text></inline-eqn>—calculation by expressions (<eqnref linkend="nj286109eqn9">9</eqnref>) and (<eqnref linkend="nj286109eqn10">10</eqnref>), for dust localized in the gap center (a) and near the electrode sheaths (b).</p></caption></figure><p>In figure <figref linkend="nj286109fig4">4</figref>, the self-consistent potential <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn> is shown and compared with the analytical evaluation of <inline-eqn><math-text>ϕ<sub><upright>OML</upright></sub></math-text></inline-eqn>. As the particle charge and potential are negative, the absolute value is plotted in figure <figref linkend="nj286109fig4">4</figref> for convenience. The <inline-eqn><math-text>ϕ<sub><upright>OML</upright></sub></math-text></inline-eqn> was obtained for the parameters shown in figures <figref linkend="nj286109fig1">1</figref> and <figref linkend="nj286109fig2">2</figref>, with <inline-eqn><math-text><italic>T</italic><sub><upright>e</upright>, <upright>i</upright></sub>=2ϵ<sub><upright>e</upright>, <upright>i</upright></sub>/3<italic>k</italic></math-text></inline-eqn>. The disagreement between the self-consistent potential and <inline-eqn><math-text>ϕ<sub><upright>OML</upright></sub></math-text></inline-eqn> is mainly due to the fact that the energy distribution functions <inline-eqn><math-text><italic>f</italic><sub><upright>e</upright>, <upright>i</upright></sub></math-text></inline-eqn> are different from Maxwellian. In the central area of the discharge, where ions have thermal energy and Maxwellian energy distribution, the disagreement is minimal and does not exceed 0.3 V. The small local maximum on both curves in the gap center for the dust profile in figure <figref linkend="nj286109fig4">4</figref>(b) is explained by the minimum of the ion density where the dust has a void (see figure <figref linkend="nj286109fig1">1</figref>(b)) and coincides for both potentials. But near the electrode sheaths, where the electric field rises and the ion drift velocity <inline-eqn><math-text><italic>u</italic></math-text></inline-eqn> directed to the electrode increases (see figure <figref linkend="nj286109fig2">2</figref>), the self-consistent <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn> and <inline-eqn><math-text>ϕ<sub><upright>OML</upright></sub></math-text></inline-eqn> begin to diverge and at the edge of dust cloud the difference reaches 1.7 V or 25%.</p><figure id="nj286109fig4" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="indented"><graphic-file version="print" format="EPS" width="21.6pc" printcolour="no" filename="images/nj286109fig4.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj286109fig4.jpg"></graphic-file></graphic><caption type="figure" id="nj286109fc4" label="Figure 4"><p indent="no">Absolute values of the dust floating potentials <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn>, obtained by self-consistent calculation (solid curves) and <inline-eqn><math-text>ϕ<sub><upright>OML</upright></sub></math-text></inline-eqn> from equations (<eqnref linkend="nj286109eqn7">7</eqnref>) and (<eqnref linkend="nj286109eqn8">8</eqnref>) (dashed curves), for dust localized in the gap center (a) and near the electrode sheaths (b).</p></caption></figure><p>The most probable reason for this disagreement is the increase of ion current to the dust in the case of fast ion drift (or large electrical field) in relation to the OML-based expression (<eqnref linkend="nj286109eqn8">8</eqnref>). This may be explained by the shape of the IEDF, which must be different from the shifted Maxwellian in the presence of ion resonant charge exchange. Since the resonant exchange produces ions with thermal energy, the abundance of low-energy ions having <inline-eqn><math-text><italic>v</italic>≪<italic>u</italic></math-text></inline-eqn> is much larger than the shifted Maxwellian IEDF gives. In order to illustrate this fact, in figure <figref linkend="nj286109fig5">5</figref>(a) the PIC-MCC and shifted Maxwellian IEDFs are shown for <inline-eqn><math-text><italic>x</italic>=1.1 <upright>cm</upright></math-text></inline-eqn>. At this point, the ion drift velocity <inline-eqn><math-text><italic>u</italic>=500 <upright>m</upright> <upright>s</upright><sup>−1</sup></math-text></inline-eqn> is two times larger than the thermal one <inline-eqn></inline-eqn>, where <inline-eqn><math-text><italic>T</italic><sub><upright>n</upright></sub>=300 <upright>K</upright></math-text></inline-eqn> is the neutral gas temperature (see figure <figref linkend="nj286109fig2">2</figref>). The shifted Maxwellian IEDF is plotted for the same <inline-eqn><math-text><italic>u</italic></math-text></inline-eqn> and <inline-eqn><math-text><italic>T</italic><sub><italic>n</italic></sub></math-text></inline-eqn>. Both functions are normalized as <inline-eqn><math-text>∫<italic>f</italic>(<italic>E</italic>) <upright>d</upright><italic>E</italic>=1</math-text></inline-eqn>. For the PIC-MCC case the distribution: (i) is shifted to lower energy, and (ii) has a larger population of high energy ions. The combination of fast ion drift and continuous thermalization of the ions, due to the resonant charge exchange collisions, increases the ion current on the dust surface in comparison to the case of ion drift with the shifted Maxwellian distribution. For the equal ion flux toward the dust particle, in the first case the probability of ion absorbtion is much larger. Hence, the absolute value of <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn> reduces toward the electrode sheaths in comparison with the behavior of <inline-eqn><math-text>ϕ<sub><upright>OML</upright></sub></math-text></inline-eqn>.</p><figure id="nj286109fig5" parts="single" width="column" position="float" pageposition="top" printstyle="normal" orientation="port"><graphic position="indented"><graphic-file version="print" format="EPS" width="21.3pc" printcolour="no" filename="images/nj286109fig5.eps"></graphic-file><graphic-file version="ej" format="JPEG" printcolour="no" filename="images/nj286109fig5.jpg"></graphic-file></graphic><caption type="figure" id="nj286109fc5" label="Figure 5"><p indent="no"> Comparison of simulated energy distribution functions (solid curves) and analytic approximation (dashed curves): (a) for ions moving with resonant charge exchange for drift velocity <inline-eqn></inline-eqn>, <inline-eqn><math-text><italic>T</italic><sub><upright>n</upright></sub>=300 <upright>K</upright></math-text></inline-eqn> and the shifted Maxwellian for the same <inline-eqn><math-text><italic>u</italic></math-text></inline-eqn> and <inline-eqn><math-text><italic>T</italic><sub><upright>n</upright></sub></math-text></inline-eqn> and (b) the same for electrons with the Maxwellian for <inline-eqn><math-text><italic>T</italic><sub><upright>e</upright></sub>=2.9 <upright>eV</upright></math-text></inline-eqn>. Both simulated functions and <inline-eqn><math-text><italic>T</italic><sub><upright>e</upright></sub></math-text></inline-eqn> are taken at <inline-eqn><math-text><italic>x</italic>=1.1 <upright>cm</upright></math-text></inline-eqn>.</p></caption></figure><p>It should be mentioned that the EEDF in this simulation was close to Maxwellian (see figure <figref linkend="nj286109fig5">5</figref>(b) and changed weakly over the region containing dust. So the electron current on the dust is adequately described by equation (<eqnref linkend="nj286109eqn7">7</eqnref>).</p><p>The incorrectness of the analytical evaluation of <inline-eqn><math-text>ϕ<sub><upright>OML</upright></sub></math-text></inline-eqn> in the presence of ion resonant charge exchange was noticed in [<cite linkend="nj286109bib11">11</cite>]. There the ion capturing by charged dust particles was simulated in the absence of an external electric field. It was also found that despite the condition of applicability of the OML approach being formally fulfilled (<inline-eqn><math-text><italic>l</italic><sub><upright>e</upright>, <upright>i</upright></sub>≫λ<sub><upright>D</upright></sub>≫<italic>r</italic><sub><upright>d</upright></sub></math-text></inline-eqn>), the analytical evaluation of <inline-eqn><math-text>ϕ<sub><upright>OML</upright></sub></math-text></inline-eqn> does not always provide an accurate value for the floating potential. In other work [<cite linkend="nj286109bib12">12</cite>], where a direct current stratified glow discharge was considered, the dust potential was found using equation (<eqnref linkend="nj286109eqn8">8</eqnref>) for the ion current and equation (<eqnref linkend="nj286109eqn5">5</eqnref>) for the electron component, with the electron distribution function obtained by solving the Boltzmann equation. Like in the present work, the non-Maxwellian approach gives up to 2 V difference from the OML results.</p><p>Note that our model does not consider the detailed ion trajectories near the particles. The presence of dust in the discharge is described by giving the probability of electrons and ions capturing during their motion through the cross section <inline-eqn><math-text>σ<sub><upright>id</upright></sub></math-text></inline-eqn> (<eqnref linkend="nj286109eqn4">4</eqnref>), which depends on the local dust concentration <inline-eqn><math-text><italic>n</italic><sub><upright>d</upright></sub>(<italic>x</italic>)</math-text></inline-eqn>, potential <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub>(<italic>x</italic>)</math-text></inline-eqn> and the instant electron/ion kinetic energy. For small dust particles <inline-eqn><math-text><italic>r</italic><sub><upright>d</upright></sub>≪λ<sub><upright>D</upright></sub></math-text></inline-eqn> the capture cross section <inline-eqn><math-text>σ<sub><upright>id</upright></sub></math-text></inline-eqn> does not depend on the details of the ion motion near the dust particle [<cite linkend="nj286109bib10">10</cite>]. According to equations (<eqnref linkend="nj286109eqn5">5</eqnref>) and (<eqnref linkend="nj286109eqn6">6</eqnref>), the most important thing is the shape of the IEDF (assumed to be Maxwellian in the OML approximation).</p><p>Let us notice that in the central region of the discharge gap, where charge resonant exchange also occurs, the OML approach gives good agreement with our self-consistent calculation (see figure <figref linkend="nj286109fig4">4</figref>). This is explained by the fact that the ion drift in this region is slow (<inline-eqn></inline-eqn>) and the IEDF is close to Maxwellian with <inline-eqn><math-text><italic>T</italic><sub><upright>i</upright></sub>=<italic>T</italic><sub><upright>n</upright></sub></math-text></inline-eqn> (neutral gas temperature), hence the resonant exchange does not influence the ion energy spectrum. But, when the ion drift velocity increases and <inline-eqn><math-text><italic>u</italic></math-text></inline-eqn> becomes comparable to the thermal velocity, the IEDF loses the Maxwellian shape and the disagreement between the PIC-MCC <inline-eqn><math-text>ϕ<sub><upright>d</upright></sub></math-text></inline-eqn> and OML model results becomes visible.</p></sec-level1><sec-level1 id="nj286109s5" label="5"><heading>Conclusion</heading><p indent="no">A new method for calculation of dust charging in a ccrf discharge was developed based on the kinetic treatment of the electron and ion motion using the PIC-MCC method. The floating potential of the dust particles is calculated self-consistently with discharge dynamics through the integral balance of ion and electron currents onto the dust surface, taking into account the energy distribution functions of both the carriers. The PIC-MCC simulations show that the OML-approximation can give incorrect results even when the conditions of its applicability (<inline-eqn><math-text><italic>l</italic><sub><upright>e</upright>, <upright>i</upright></sub>≫λ<sub><upright>D</upright></sub>≫<italic>r</italic><sub><upright>d</upright></sub></math-text></inline-eqn>) are fulfilled. In the presence of ion resonant charge exchange collisions, the OML-approximation gives correct results only in a weak electric field, where the ion drift velocity is much smaller than the thermal one. With increasing ion drift velocity, the OML-approximation gives up to 25% larger absolute value of the dust potential, than that obtained with the PIC-MCC calculations. This disagreement may be explained by the shape of the IEDF, which is different from Maxwellian or shifted Maxwellian. For the case of fast ion drift with resonant charge exchange, the abundance of low energy ions in the energy distribution shape is larger than that expected from a Maxwellian approximation, so the ion current on dust is enhanced in relation to the traditional OML-based expression. Our numerical approach allows to take into account the shape of the real energy distribution and its effect on the charging of the dust particles.</p></sec-level1><acknowledgment><heading>Acknowledgments</heading><p indent="no">This work was supported by the Bilateral grants RFBR-Flanders (05-02-19809-MFa) and RFBR-Ukraine (08-02-90446-Uka).</p></acknowledgment></body><back><references><heading>References</heading><reference-list type="numeric"><journal-ref id="nj286109bib1" num="1"><authors><au><second-name>Akdim</second-name><first-names>M R</first-names></au><au><second-name>Goedheer</second-name><first-names>W J</first-names></au></authors><year>2003</year><jnl-title>Phys. Rev.</jnl-title><part>E</part><volume>67</volume><pages>066407</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevE.67.066407</cr_doi><cr_issn type="print">15393755</cr_issn><cr_issn type="electronic">15502376</cr_issn></crossref></journal-ref><journal-ref id="nj286109bib2" num="2"><authors><au><second-name>Godyak</second-name><first-names>V A</first-names></au><au><second-name>Piejak</second-name><first-names>R B</first-names></au><au><second-name>Alexandrovich</second-name><first-names>B M</first-names></au></authors><year>1992</year><jnl-title>Plasma Sources Sci. Technol.</jnl-title><volume>1</volume><pages>36</pages><crossref><cr_doi>http://dx.doi.org/10.1088/0963-0252/1/1/006</cr_doi><cr_issn type="print">09630252</cr_issn><cr_issn type="electronic">13616595</cr_issn></crossref></journal-ref><journal-ref id="nj286109bib3" num="3"><authors><au><second-name>Denysenko</second-name><first-names>I</first-names></au><au><second-name>Yu</second-name><first-names>M Y</first-names></au><au><second-name>Ostrikov</second-name><first-names>K</first-names></au><au><second-name>Azarenkov</second-name><first-names>N A</first-names></au><au><second-name>Stenflo</second-name><first-names>L</first-names></au></authors><year>2004</year><jnl-title>Phys. Plasmas</jnl-title><volume>11</volume><pages>4959</pages><crossref><cr_doi>http://dx.doi.org/10.1063/1.1793174</cr_doi><cr_issn type="print">1070664X</cr_issn></crossref></journal-ref><journal-ref id="nj286109bib4" num="4"><authors><au><second-name>Chutov</second-name><first-names>Y I</first-names></au><au><second-name>Goedheer</second-name><first-names>W J</first-names></au></authors><year>2003</year><jnl-title>IEEE Trans. Plasma Sci.</jnl-title><volume>31</volume><pages>606</pages><crossref><cr_doi>http://dx.doi.org/10.1109/TPS.2003.815490</cr_doi><cr_issn type="print">00933813</cr_issn></crossref></journal-ref><journal-ref id="nj286109bib5" num="5"><authors><au><second-name>Schweigert</second-name><first-names>I V</first-names></au><au><second-name>Peeters</second-name><first-names>F M</first-names></au></authors><year>2007</year><jnl-title>JETP Lett.</jnl-title><volume>86</volume><pages>572</pages></journal-ref><journal-ref id="nj286109bib6" num="6"><authors><au><second-name>Schweigert</second-name><first-names>I V</first-names></au><au><second-name>Alexandrov</second-name><first-names>A L</first-names></au><au><second-name>Ariskin</second-name><first-names>D A</first-names></au><au><second-name>Peeters</second-name><first-names>F M</first-names></au><au><second-name>Stefanović</second-name><first-names>I</first-names></au><au><second-name>Kovačević</second-name><first-names>E</first-names></au><au><second-name>Berndt</second-name><first-names>J</first-names></au><au><second-name>Winter</second-name><first-names>J</first-names></au></authors><year>2008</year><jnl-title>Phys. Rev.</jnl-title><part>E</part><volume>78</volume><pages>026410</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRevE.78.026410</cr_doi><cr_issn type="print">15393755</cr_issn><cr_issn type="electronic">15502376</cr_issn></crossref></journal-ref><book-ref id="nj286109bib7" num="7"><authors><au><second-name>Hockney</second-name><first-names>R W</first-names></au><au><second-name>Eastwood</second-name><first-names>J W</first-names></au></authors><year>1981</year><book-title>Computer Simulation Using Particles</book-title><publication><place>New York</place><publisher>McGraw-Hill</publisher></publication></book-ref><journal-ref id="nj286109bib8" num="8"><authors><au><second-name>Lagushenko</second-name><first-names>R</first-names></au><au><second-name>Maya</second-name><first-names>J</first-names></au></authors><year>1984</year><jnl-title>J. Appl. Phys.</jnl-title><volume>55</volume><pages>3293</pages><crossref><cr_doi>http://dx.doi.org/10.1063/1.333365</cr_doi><cr_issn type="print">00218979</cr_issn></crossref></journal-ref><journal-ref id="nj286109bib9" num="9"><authors><au><second-name>Mott-Smith</second-name><first-names>H M</first-names></au><au><second-name>Langmuir</second-name><first-names>I</first-names></au></authors><year>1926</year><jnl-title>Phys. Rev.</jnl-title><volume>28</volume><pages>727</pages><crossref><cr_doi>http://dx.doi.org/10.1103/PhysRev.28.727</cr_doi><cr_issn type="print">0031899X</cr_issn><cr_issn type="electronic">15366065</cr_issn></crossref></journal-ref><journal-ref id="nj286109bib10" num="10"><authors><au><second-name>Fortov</second-name><first-names>V E</first-names></au><au><second-name>Khrapak</second-name><first-names>A G</first-names></au><au><second-name>Khrapak</second-name><first-names>S A</first-names></au><au><second-name>Molotkov</second-name><first-names>V I</first-names></au><au><second-name>Petrov</second-name><first-names>O F</first-names></au></authors><year>2004</year><jnl-title>Phys. Usp.</jnl-title><volume>47</volume><pages>447</pages><crossref><cr_doi>http://dx.doi.org/10.1070/PU2004v047n05ABEH001689</cr_doi><cr_issn type="print">10637869</cr_issn><cr_issn type="electronic">14684780</cr_issn></crossref></journal-ref><journal-ref id="nj286109bib11" num="11"><authors><au><second-name>Zobnin</second-name><first-names>A V</first-names></au><au><second-name>Nefedov</second-name><first-names>A P</first-names></au><au><second-name>Sinelshchikov</second-name><first-names>V A</first-names></au><au><second-name>Fortov</second-name><first-names>V E</first-names></au></authors><year>2000</year><jnl-title>J. Exp. Theor. Phys.</jnl-title><volume>91</volume><pages>483</pages><crossref><cr_doi>http://dx.doi.org/10.1134/1.1320081</cr_doi><cr_issn type="print">10637761</cr_issn><cr_issn type="electronic">10906509</cr_issn></crossref></journal-ref><journal-ref id="nj286109bib12" num="12"><authors><au><second-name>Sukhinin</second-name><first-names>G I</first-names></au><au><second-name>Fedoseev</second-name><first-names>A V</first-names></au></authors><year>2007</year><jnl-title>Plasma Phys. Rep.</jnl-title><volume>33</volume><pages>1023</pages><crossref><cr_doi>http://dx.doi.org/10.1134/S1063780X07120070</cr_doi><cr_issn type="print">1063780X</cr_issn><cr_issn type="electronic">15626938</cr_issn></crossref></journal-ref></reference-list></references></back></article></istex:document></istex:metadataXml><mods version="3.6"><titleInfo lang="eng"><title>A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas</title></titleInfo><titleInfo type="abbreviated"><title>A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas</title></titleInfo><titleInfo type="alternative" lang="eng"><title>A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas</title></titleInfo><name type="personal"><namePart type="given">A L</namePart><namePart type="family">Alexandrov</namePart><affiliation>Institute of Theoretical and Applied Mechanics, Russian Academy of Sciences, Novosibirsk 630090, Russia</affiliation><affiliation>E-mail: aalex@itam.nsc.ru</affiliation><role><roleTerm type="text">author</roleTerm></role></name><name type="personal"><namePart type="given">I V</namePart><namePart type="family">Schweigert</namePart><affiliation>Institute of Theoretical and Applied Mechanics, Russian Academy of Sciences, Novosibirsk 630090, Russia</affiliation><affiliation>Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</affiliation><affiliation>E-mail: ischweig@yahoo.com</affiliation><role><roleTerm type="text">author</roleTerm></role></name><name type="personal"><namePart type="given">F M</namePart><namePart type="family">Peeters</namePart><affiliation>Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium</affiliation><role><roleTerm type="text">author</roleTerm></role></name><typeOfResource>text</typeOfResource><genre type="research-article" displayLabel="paper">paper</genre><originInfo><publisher>Institute of Physics Publishing</publisher><dateIssued encoding="w3cdtf">2008</dateIssued><copyrightDate encoding="w3cdtf">2008</copyrightDate></originInfo><language><languageTerm type="code" authority="iso639-2b">eng</languageTerm><languageTerm type="code" authority="rfc3066">en</languageTerm></language><physicalDescription><internetMediaType>text/html</internetMediaType></physicalDescription><abstract>Nanoparticle charging in a capacitively coupled radio frequency discharge in argon is studied using a particle in cell Monte Carlo collisions method. The plasma parameters and dust potential were calculated self-consistently for different unmovable dust profiles. A new method for definition of the dust floating potential is proposed, based on the information about electron and ion energy distribution functions, obtained during the kinetic simulations. This approach provides an accurate balance of the electron and ion currents on the dust particle surface and allows us to precisely calculate the dust floating potential. A comparison of the obtained floating potentials with the results of the traditional orbital motion limit (OML) theory shows that in the presence of the ion resonant charge exchange collisions, even when the OML approximation is valid, its results are correct only in the region of a weak electric field, where the ion drift velocity is much smaller than the thermal one. With increasing ion drift velocity, the absolute value of the calculated dust potential becomes significantly smaller than the theory predicts. This is explained by a non-Maxwellian shape of the ion energy distribution function for the case of fast ion drift.</abstract><relatedItem type="host"><titleInfo><title>New Journal of Physics</title></titleInfo><titleInfo type="abbreviated"><title>New J. Phys.</title></titleInfo><genre type="journal">journal</genre><identifier type="ISSN">1367-2630</identifier><identifier type="eISSN">1367-2630</identifier><identifier type="PublisherID">nj</identifier><identifier type="CODEN">NJOPFM</identifier><identifier type="URL">stacks.iop.org/NJP</identifier><part><date>2008</date><detail type="volume"><caption>vol.</caption><number>10</number></detail><detail type="issue"><caption>no.</caption><number>9</number></detail><extent unit="pages"><start>1</start><end>12</end><total>12</total></extent></part></relatedItem><identifier type="istex">B74E140A61C686AC8CB111B0FC21BA43C88A2403</identifier><identifier type="DOI">10.1088/1367-2630/10/9/093025</identifier><identifier type="articleID">286109</identifier><identifier type="articleNumber">093025</identifier><accessCondition type="use and reproduction" contentType="copyright">IOP Publishing and Deutsche Physikalische Gesellschaft</accessCondition><recordInfo><recordContentSource>IOP</recordContentSource><recordOrigin>IOP Publishing and Deutsche Physikalische Gesellschaft</recordOrigin></recordInfo></mods></metadata><enrichments><json:item><type>multicat</type><uri>https://api.istex.fr/document/B74E140A61C686AC8CB111B0FC21BA43C88A2403/enrichments/multicat</uri></json:item></enrichments><serie></serie></istex></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Belgique/explor/OpenAccessBelV2/Data/Istex/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 001361 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Istex/Corpus/biblio.hfd -nk 001361 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Belgique
|area= OpenAccessBelV2
|flux= Istex
|étape= Corpus
|type= RBID
|clé= ISTEX:B74E140A61C686AC8CB111B0FC21BA43C88A2403
|texte= A non-Maxwellian kinetic approach for charging of dust particles in discharge plasmas
}}
| This area was generated with Dilib version V0.6.25. |  |