<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JMP</journal-id><journal-title-group><journal-title>Journal of Modern Physics</journal-title></journal-title-group><issn pub-type="epub">2153-1196</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jmp.2013.411187</article-id><article-id pub-id-type="publisher-id">JMP-40187</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Ion-Acoustic Higher Order Non-Linear Structures in Quantum Dusty Plasma
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>d</surname><given-names>Manir Hossain</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Monirul</surname><given-names>Hasan</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Md</surname><given-names>A saduzzaman</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Md</surname><given-names>Mahfuzul Haque</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Institute of Natural Science, United International University, Dhaka, Bangladesh</addr-line></aff><aff id="aff1"><addr-line>Department of Electrical and Electronic Engineering, Bangladesh University of Business and Technology, 
Dhaka, Bangladesh</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>monirul.h@bubt.edu.bd(MH)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>12</day><month>11</month><year>2013</year></pub-date><volume>04</volume><issue>11</issue><fpage>1530</fpage><lpage>1535</lpage><history><date date-type="received"><day>June</day>	<month>3,</month>	<year>2013</year></date><date date-type="rev-recd"><day>July</day>	<month>5,</month>	<year>2013</year>	</date><date date-type="accepted"><day>August</day>	<month>1,</month>	<year>2013</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   The prominent features of higher order nonlinear ion-acoustic waves involving quantum corrections in an unmagnetized quantum dusty plasma are revisited with the theoretical framework of Hossain et al. [1]. The fluid model is demonstrated here by its constituent inertial ions, Fermi electrons with quantum effect, and immovable dust grain with negative charge. We have used the ideology of Gardner equation. The well-known RPM method is employed to derive the equation. Indeed, the basic features of quantum dust ion-acoustic Gardner solitons (GSs) are pronounced here. GSs are shown to exist for the value of dust to ion ratio around 2/3 which is valid for space plasma [2], and are different from those of K-dV (Korteweg-de Vries) solitons, which do not exist for the value around 2/3. The implications of our results are suitable for cosmological and astrophysical environments. 
 
</p></abstract><kwd-group><kwd>Negative Dust; Quantum Plasma; Fermi Electron; Modified K-dV Equation; Gardner Solitons; RPM</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Quantum plasmas have attracted a great deal of attention because of their potential applications in dense plasma particularly in different astrophysical and cosmological systems [3-5] (e.g. interstellar or molecular clouds, planetary rings, comets, interior of white dwarf stars, etc.), in nanostructures [<xref ref-type="bibr" rid="scirp.40187-ref6">6</xref>], in microelectronic device [<xref ref-type="bibr" rid="scirp.40187-ref7">7</xref>] as well as in the next-generation intense laser [<xref ref-type="bibr" rid="scirp.40187-ref8">8</xref>]. Many authors have proposed some theories including the quantum corrections to the quantum plasma echoes [<xref ref-type="bibr" rid="scirp.40187-ref9">9</xref>], the self-consistent dynamics of Fermi gases [<xref ref-type="bibr" rid="scirp.40187-ref10">10</xref>], quantum beam instabilities [<xref ref-type="bibr" rid="scirp.40187-ref11">11</xref>], wave interactions in quantum magnetoplasmas [<xref ref-type="bibr" rid="scirp.40187-ref12">12</xref>], classical and quantum kinetics of the Zakharov system [<xref ref-type="bibr" rid="scirp.40187-ref13">13</xref>], quantum corrections to the Zakharov equations [<xref ref-type="bibr" rid="scirp.40187-ref14">14</xref>], expansion of quantum electron gas into vacuum [<xref ref-type="bibr" rid="scirp.40187-ref15">15</xref>], quantum ion acoustic waves [<xref ref-type="bibr" rid="scirp.40187-ref16">16</xref>], quantum Landau damping [<xref ref-type="bibr" rid="scirp.40187-ref17">17</xref>], magnetohydrodynamics of quantum plasmas [<xref ref-type="bibr" rid="scirp.40187-ref18">18</xref>], etc. Quantum plasmas have extremely high plasma number densities and low temperatures. At extremely low temperatures, the thermal de Broglie wavelength becomes comparable to the interelectron distance and the electron temperature becomes comparable to the electron Fermi temperature <img src="12-7501410\b62592b9-9b8d-4d0d-b76e-7b0a1adc6549.jpg" /> and the electrons follow Fermi Dirac distribution law. In this condition, quantum mechanical effects are expected to play a significant role in the behavior of charged particles [19-21]. As electrons are lighter than ions, the quantum behavior of electron is reached faster than ions. The dust particles are quite common in various plasma systems. The inclusion of immobile charged dust in electron-ion plasmas leads to introduce a new mode. Shukla and Silin [<xref ref-type="bibr" rid="scirp.40187-ref22">22</xref>] have first theoretically shown the existence of low-frequency dust ion-acoustic (DIA) waves in a dusty plasma, which was latter observed in laboratory experiments [23,24]. The phase speed of the DIA waves is much smaller (larger) than electron (ion) thermal speed. The inertia is provided by the ion mass while the restoring force comes from the electron thermal pressure. These waves differ from usual ion-acoustic waves [<xref ref-type="bibr" rid="scirp.40187-ref25">25</xref>] due to the conservation of equilibrium charge density <img src="12-7501410\bd47c356-cda1-4088-a14f-19040ef02f3d.jpg" /> and the strong inequality<img src="12-7501410\83b96fde-60b5-4a8c-b412-6cbba273c927.jpg" />, where <img src="12-7501410\eff8cdb5-3678-4d78-97a0-5c11bdcda1eb.jpg" /> is the particle number density of the species <img src="12-7501410\4556f53a-efed-4479-86c1-c93ef260ee77.jpg" /> with <img src="12-7501410\74daf1d1-f305-4389-b334-b7f6950b2103.jpg" /> for electrons, <img src="12-7501410\c696e29a-8715-4d98-b541-fa8f179d2239.jpg" />for ions and <img src="12-7501410\e4367ac8-44a6-4ff7-bc75-4946d5bdad14.jpg" /> for dust, <img src="12-7501410\b0a82f0f-ccda-4929-a968-cfefee80cbe4.jpg" />is the number of electrons residing onto the dust grain surface, and <img src="12-7501410\b257d389-9fe8-43be-9ada-ed6cc2d22c75.jpg" /> is the magnitude of an electronic charge. Therefore, a dusty plasma can not support the usual ion-acoustic waves, but can do the DIA waves of Shukla and Silin [<xref ref-type="bibr" rid="scirp.40187-ref22">22</xref>]. The nonlinear waves associated with the DIA and QDIA waves particularly solitary waves (SWs) [26,27] and shock waves [<xref ref-type="bibr" rid="scirp.40187-ref28">28</xref>] have received a great deal of interest in understanding the basic properties of localized electrostatic perturbation in space [29,30] and laboratory dusty plasmas [31-34]. A number of investigations have been made on QDIA SWs [<xref ref-type="bibr" rid="scirp.40187-ref35">35</xref>] and shocks [<xref ref-type="bibr" rid="scirp.40187-ref28">28</xref>] by using K-dV equation. For plasmas with more than two species, it can arise cases where the K-dV equation is not valid near a critical value of a certain parameter (say<img src="12-7501410\40078c55-4587-44ad-9186-61c7a629e6c1.jpg" />). The nonlinear term vanishes at this critical value (at<img src="12-7501410\c25a4f7d-1f69-4c67-8675-f8d8a987ae1d.jpg" />) [<xref ref-type="bibr" rid="scirp.40187-ref36">36</xref>] which makes soliton amplitude large enough to break down the validity of K-dV equation. The mmK-dV (mixed modified K-dV) equation, also known as Gardner equation, can give soliton solution around this critical value [<xref ref-type="bibr" rid="scirp.40187-ref37">37</xref>]. The technique of analyzing SWs is Gardner approach which leads to a standard Gardner equation. From the analysis of standard Gardner equation, SW of permanent profile is found, which is known as Gardner soliton (GS) [1,38,39]. In our present manuscript, we attempt to study the basic features of QDIA GSs by deriving modified Gardner equation, which is valid around<img src="12-7501410\c379a8d0-0ae3-4db9-85e0-7efd59e3426c.jpg" />, in a quantum dusty plasma containing inertial ions, Fermi electrons with quantum effect, and negatively charged immobile dust. The manuscript is organized as follows. The model equations are provided in Section 2. The Gardner equation is derived by using the reductive perturbation method in Section 3. The analytical solutions are presented in Section 4. A brief discussion is finally given in Section 5.</p></sec><sec id="s2"><title>2. Model Equations</title><p>We consider a one-dimensional, collisionless, unmagnetized quantum dusty plasma system composed of inertial ions, massless Fermi electrons with quantum effect, and negatively charged immobile dust. Thus, at equilibrium we have<img src="12-7501410\8121f7e0-0de0-4fd9-be06-73d71d047e8a.jpg" />. The nonlinear dynamics of these low-frequency (purely electrostatic) QDIA waves in such a plasma system is described by the normalized equations of the form</p><disp-formula id="scirp.40187-formula29746"><label>(1)</label><graphic position="anchor" xlink:href="12-7501410\55aefe1b-d6ff-47a7-a816-4482d68059cc.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29747"><label>(2)</label><graphic position="anchor" xlink:href="12-7501410\cde422a9-7ea9-4e76-a0c7-0d9650068e4c.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29748"><label>(3)</label><graphic position="anchor" xlink:href="12-7501410\0f616d50-acff-4f34-ba45-1d02e8af9b5d.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29749"><label>(4)</label><graphic position="anchor" xlink:href="12-7501410\37ca1396-90a3-461c-a699-b47df3f30689.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29750"><label>(5)</label><graphic position="anchor" xlink:href="12-7501410\6f33c6d3-8a0a-4fa9-b3db-20f4af6bd97d.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="12-7501410\01ada8d2-e594-4cbc-b716-b612f3505a28.jpg" /> is the ion (electron) number density normalized by its equilibrium value<img src="12-7501410\49d4bad3-130d-480c-bba8-4e27cef76c65.jpg" />, <img src="12-7501410\e30bae23-8938-4313-bf58-6d7c10de5768.jpg" />is the ion fluid speed normalized by quantum ion-acoustic speed <img src="12-7501410\00c963bc-5252-439a-84c3-29296ae69553.jpg" /> with <img src="12-7501410\54b3393a-ebb3-4dcc-9174-56603fe2b29c.jpg" /> being the ion rest mass, <img src="12-7501410\54b1a98b-79a3-424d-995d-05105f216a4d.jpg" />is the electron Fermi energy, <img src="12-7501410\c8beadde-442a-49f2-8c8a-19457550354d.jpg" />is the Boltzmann constant, and T<sub>Fe</sub> is the Fermi temperature of electron, <img src="12-7501410\80d5d0b4-252c-4926-a3e1-8d2b9b2e87e5.jpg" />is the electrostatic wave potential normalized by <img src="12-7501410\c36c509e-d57f-4ba6-8061-2ced635eeb3f.jpg" /> with e being the magnitude of the charge of an electron, ρ is the normalized surface charge density, and<img src="12-7501410\ca4906b8-ba7c-44d0-849f-45537e57a0b7.jpg" />. The time variable <img src="12-7501410\59195e77-43ab-4d8b-985c-ccf9503ff0e9.jpg" /> is normalized by <img src="12-7501410\61df6d24-9c69-4037-8170-b7b2d892c8ab.jpg" /> and the space variable is normalized by<img src="12-7501410\9a2d56d7-bdb5-4671-aa89-84e381caf31c.jpg" />. In Equation (3) we have used the following Fermi pressure law for the electron species [40,41]</p><disp-formula id="scirp.40187-formula29751"><label>(6)</label><graphic position="anchor" xlink:href="12-7501410\4afd067b-831e-4feb-88b6-1b4cff14b194.jpg"  xlink:type="simple"/></disp-formula><p>Also <img src="12-7501410\27758132-212e-4aa3-adf8-e19812722a99.jpg" /> with <img src="12-7501410\ae93b65e-db93-4e5a-b3cb-c7cae51b3e1e.jpg" /> and <img src="12-7501410\33db8047-3270-40b5-9d47-2d34731f2d47.jpg" /> is the ratio between the plasmon energy and the electron Fermi energy where <img src="12-7501410\0b3493c0-9d95-4ca7-b969-d6f13e4f73d3.jpg" /> is the electron Fermi speed at temperature<img src="12-7501410\fba02cd7-a9a0-4ba2-9d24-898ece6e018d.jpg" />.</p></sec><sec id="s3"><title>3. Gardner Equation</title><p>We first obtain the well known K-dV equation and see why Gardner equation is needed to find SW solution.</p><sec id="s3_1"><title>3.1. Derivation of the K-dV Equation</title><p>To obtain the QDIA K-dV equation, we introduce the stretched coordinates</p><disp-formula id="scirp.40187-formula29752"><label>(7)</label><graphic position="anchor" xlink:href="12-7501410\0d7a08da-edf8-4369-bed2-86e37d2d63d6.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29753"><label>(8)</label><graphic position="anchor" xlink:href="12-7501410\9db631ec-e830-4a16-9c2b-56e75972a3c1.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="12-7501410\6edb9518-ed85-4ce3-96ab-2bdc86d7a7f3.jpg" /> is the QDIA wave phase speed <img src="12-7501410\97957ee9-b0df-48d1-ae2e-4f9f7968944c.jpg" /> and <img src="12-7501410\73cecf9c-c47b-453f-a551-8534165d7df0.jpg" /> is a smallness parameter measuring the weakness of the dispersion<img src="12-7501410\da4ad6f8-2438-4696-b9d5-d9a0aa93a367.jpg" />. We then expand<img src="12-7501410\53cbe1a2-5dbb-41c6-b53c-446f2785457c.jpg" />, <img src="12-7501410\d75d7422-f855-494f-a745-7493ad55ff02.jpg" />, <img src="12-7501410\396d2266-9193-4ae8-b97c-fdb564a1a749.jpg" />, <img src="12-7501410\5c78564f-8488-49bb-979f-8d62cfdaa9b2.jpg" />, and <img src="12-7501410\7985101d-72ab-45f8-8f1e-24f1786613cc.jpg" /> in power series of <img src="12-7501410\3cdb7659-6e59-4328-aada-936045bf7b42.jpg" /></p><disp-formula id="scirp.40187-formula29754"><label>(9)</label><graphic position="anchor" xlink:href="12-7501410\59919680-4859-4ef0-af3b-11c2b7293c4e.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29755"><label>(10)</label><graphic position="anchor" xlink:href="12-7501410\9e2b47d1-6a1a-4445-9ff5-454f8606196c.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29756"><label>(11)</label><graphic position="anchor" xlink:href="12-7501410\679cf168-d391-4bb6-bb79-7721e171d1f2.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29757"><label>(12)</label><graphic position="anchor" xlink:href="12-7501410\94b6bf08-bff3-4e01-a318-45561ef17f56.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29758"><label>(13)</label><graphic position="anchor" xlink:href="12-7501410\49962e1d-a43b-4b48-961c-9b694f74120a.jpg"  xlink:type="simple"/></disp-formula><p>and develop equations in various powers of<img src="12-7501410\91273c8a-ae70-4f17-98da-b0b5f247d04a.jpg" />. To the lowest order in<img src="12-7501410\b7f6a62b-4d18-44b7-81a6-0dbfcaf3dcdb.jpg" />, Equations (1)-(13) give</p><disp-formula id="scirp.40187-formula29759"><label>(14)</label><graphic position="anchor" xlink:href="12-7501410\c41b0b3b-9885-4300-8141-641d367c6bb5.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29760"><label>(15)</label><graphic position="anchor" xlink:href="12-7501410\d8ac61e7-0834-4efe-9839-fcbefaa22a79.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29761"><label>(16)</label><graphic position="anchor" xlink:href="12-7501410\0c0d2148-5f60-4c9b-a4d0-066ed3ba8447.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29762"><label>(17)</label><graphic position="anchor" xlink:href="12-7501410\5fbc3ad0-e9fe-4868-b84a-7d62baf2c4c6.jpg"  xlink:type="simple"/></disp-formula><p>where<img src="12-7501410\0fac30e2-814d-4f80-a985-76469a0f6fce.jpg" />. Equation (17) represents the linear dispersion relation for the QDIA waves. This clearly indicates that the QDIA wave phase speed <img src="12-7501410\091200e6-095d-43eb-8785-62e726e32ae4.jpg" /> increases with the increase of the dust charge density<img src="12-7501410\cb023e9e-6934-4dbe-9173-9295a3af2de7.jpg" />.</p><p>To the next higher order of<img src="12-7501410\371e00d9-e5f9-4af8-847a-c4cf0391a5ea.jpg" />, one can obtain another set of coupled equations for<img src="12-7501410\62d83e2e-09b7-4f27-8c36-6c6aeaf56611.jpg" />, <img src="12-7501410\66262f22-c23a-441c-9b7c-d75c90c622b7.jpg" />, and<img src="12-7501410\44528180-36b8-4b25-bc86-65e8fcaa1024.jpg" />, which -along with the first set of coupled linear equations for<img src="12-7501410\5bc2ca59-f0b4-47e9-b354-4702ddcb8bc9.jpg" />, <img src="12-7501410\01201cb8-e9e2-4d7c-bd1d-1d94a47346ef.jpg" />, and <img src="12-7501410\625b16c4-57ef-445b-aba8-cde7b5d4722d.jpg" />-reduce to a nonlinear dynamical equation of the form</p><disp-formula id="scirp.40187-formula29763"><label>(18)</label><graphic position="anchor" xlink:href="12-7501410\3fff7adc-9d35-4e84-994d-74e2b61d19b7.jpg"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.40187-formula29764"><label>(19)</label><graphic position="anchor" xlink:href="12-7501410\dfa5aadf-a187-498e-a4b0-a39bd6563b4b.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29765"><label>(20)</label><graphic position="anchor" xlink:href="12-7501410\53f4181e-f089-4d13-98a9-e17fc53bd682.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29766"><label>(21)</label><graphic position="anchor" xlink:href="12-7501410\01788ffc-3a24-430e-8be4-5eea6e3b895e.jpg"  xlink:type="simple"/></disp-formula><p>Equation (18) is known as K-dV equation. The stationary localized solution of Equation (18) is given by</p><disp-formula id="scirp.40187-formula29767"><label>(22)</label><graphic position="anchor" xlink:href="12-7501410\e7187d36-6364-4a10-afce-a79f3ec0e0b5.jpg"  xlink:type="simple"/></disp-formula><p>where the amplitude <img src="12-7501410\b81bb8c3-bdcb-4f79-aeb3-81b3a45a7e4c.jpg" /> and the width <img src="12-7501410\fb48f48b-0c06-4721-9cb2-72f7c571def5.jpg" /> are given by <img src="12-7501410\87117e5d-f6e3-43e1-a6b4-e9d23a28a48a.jpg" /> and<img src="12-7501410\05ccac0e-23c6-4bb3-8fe2-246a8ec0cf6a.jpg" />, respectively. <img src="12-7501410\b8334270-6b22-4885-891f-521f709af6ad.jpg" />is the mach number. As <img src="12-7501410\6c06b42b-1fa6-41ce-9b5c-54274bf1bd65.jpg" /> and<img src="12-7501410\32d729f8-8845-4d6f-8b24-e4f066e56f1f.jpg" /><img src="12-7501410\0a2bef58-ac10-4150-9b85-e302b2b8ed0f.jpg" />, (22) clearly indicates that 1) small amplitude solitary waves with<img src="12-7501410\28af80e4-5cb4-444d-a1b2-69851eb06128.jpg" />, i.e. positive soliton exists if<img src="12-7501410\75ac3c0b-f1c7-43b3-9ae0-909c6fa0973c.jpg" />, 2) small amplitude solitary waves with<img src="12-7501410\b2a8e5e2-52db-463e-b2ac-9eda29b87c22.jpg" />, i.e. negative soliton exists if<img src="12-7501410\395c84fb-78be-46f1-89c5-2d5c219b4c42.jpg" />, and 3) <img src="12-7501410\3b098695-42c4-446e-aacc-28c177baec47.jpg" />for <img src="12-7501410\9f1685c6-1ed7-42f8-9468-a27a1fb311a0.jpg" /> i.e. the nonlinear term vanishes at μ = μ<sub>c</sub> and is not valid near μ = μ<sub>c</sub> which makes soliton amplitude large enough to break down its validity. To find soliton solution around<img src="12-7501410\49a12aaa-cfd4-4d80-b865-3844f430afac.jpg" />, we now obtain gardner equation.</p></sec><sec id="s3_2"><title>3.2. Derivation of the Gardner Equation</title><p>To study QDIA GSs by analyzing the ingoing solutions of Equations (1)-(5), we first introduce the stretched coordinates [<xref ref-type="bibr" rid="scirp.40187-ref37">37</xref>]</p><disp-formula id="scirp.40187-formula29768"><label>(23)</label><graphic position="anchor" xlink:href="12-7501410\ad665706-b7bf-4ac3-b8bf-90a21e9f715c.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29769"><label>(24)</label><graphic position="anchor" xlink:href="12-7501410\29a6c51e-2056-40b5-82cc-8a1c2b980880.jpg"  xlink:type="simple"/></disp-formula><p>By using Equations (23) and (24) in Equations (1)-(6), and Equations (9)-(13) and to the lowest order in ε, we find the same values of<img src="12-7501410\231f09db-b024-46f8-8ae7-6f459cac15ae.jpg" />, <img src="12-7501410\731049ac-34a3-490e-9405-4d5142e6d50a.jpg" />, <img src="12-7501410\6df5bc28-45ad-4371-aba1-ab3d89fc9471.jpg" />, and <img src="12-7501410\3cd749ee-c485-452f-8d81-1a6dadee1547.jpg" /> as like as that of the K-dV. To the next higher order in<img src="12-7501410\1852cd54-65e2-4dfa-ae9b-a2cca11e3eb0.jpg" />, we obtain a set of equations, which, after using Equations (14)-(17), can be simplified as</p><disp-formula id="scirp.40187-formula29770"><label>(25)</label><graphic position="anchor" xlink:href="12-7501410\bdd75f6e-cb0b-4302-bb47-7086816f584e.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29771"><label>(26)</label><graphic position="anchor" xlink:href="12-7501410\58eaa6bc-9d77-4a9a-a3cb-d7bfaab8d9b6.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29772"><label>(27)</label><graphic position="anchor" xlink:href="12-7501410\cdd9b8ae-fbe7-4b86-a006-2cbaaf784ed3.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29773"><label>(28)</label><graphic position="anchor" xlink:href="12-7501410\fad6629c-0c9f-444b-8421-4c20d333ddcb.jpg"  xlink:type="simple"/></disp-formula><p>It is obvious from Equation (28) that <img src="12-7501410\59de5a1a-5f5d-4a6a-9286-492ca7b5d806.jpg" /> since<img src="12-7501410\9a882fec-d229-4589-b7d2-00b2aa1d2ffd.jpg" />. One can find that <img src="12-7501410\57a91700-f53c-42b8-947e-846bb7712946.jpg" /> at its critical value <img src="12-7501410\0397523c-8794-4399-8f93-9e64e9d4130e.jpg" /> (which is a solution of<img src="12-7501410\8ac2d6a4-e416-4e23-bdf8-726240703439.jpg" />). So, for <img src="12-7501410\ee66004f-b2b9-41ba-8f06-80ec765c1094.jpg" /> around its critical value<img src="12-7501410\d433a18e-dcb4-4fda-854a-e878519196ab.jpg" />, <img src="12-7501410\98a2a815-d719-4d81-97f5-66b81292a395.jpg" />can be expressed as</p><disp-formula id="scirp.40187-formula29774"><label>(29)</label><graphic position="anchor" xlink:href="12-7501410\047bcdcf-deae-4419-a745-33f44e8e19dc.jpg"  xlink:type="simple"/></disp-formula><p>where<img src="12-7501410\5275b09f-ae7e-431e-b0a8-4810c4d17ff9.jpg" />, <img src="12-7501410\c2977cf9-1817-433f-be5c-120ec4ed77e0.jpg" />is a small and dimensionless parameter, and can be taken as the expansion parameter<img src="12-7501410\e9f04af8-4f3f-445c-a91f-469beef9d8b1.jpg" />, i.e.<img src="12-7501410\f0f4c31b-5255-4b0a-86d6-e772ab6a1428.jpg" />, and <img src="12-7501410\38cd5a69-8dfb-4b83-8126-701b0eadf48d.jpg" /> for <img src="12-7501410\f8878b68-47ca-4b49-bcd6-2d10f4a30ad1.jpg" /> and <img src="12-7501410\864a4d7f-790b-4aae-891d-e7375c374b78.jpg" /> for<img src="12-7501410\ae1c3659-65ae-429d-aa94-c0366e4ba298.jpg" />. So, <img src="12-7501410\71e342b4-b5c7-47a5-996a-fabc1fbc2e6f.jpg" />can be expressed as</p><disp-formula id="scirp.40187-formula29775"><label>(30)</label><graphic position="anchor" xlink:href="12-7501410\3a28e334-807d-4e7a-b498-665c28f14a34.jpg"  xlink:type="simple"/></disp-formula><p>which, therefore, must be included in the third order Poissons equation. To the next higher order in<img src="12-7501410\f9fe4368-4180-4fb3-a341-35f7f0ad7585.jpg" />, and after some mathematical calculations we obtain a set of equations</p><disp-formula id="scirp.40187-formula29776"><label>(31)</label><graphic position="anchor" xlink:href="12-7501410\a9dc0b3e-ff70-42a9-bba7-3d205979d29e.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29777"><label>(32)</label><graphic position="anchor" xlink:href="12-7501410\0ec8c306-c191-4326-b387-6affdae5e7db.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29778"><label>(33)</label><graphic position="anchor" xlink:href="12-7501410\4eaa185b-6780-47cb-b34c-1a646246b3df.jpg"  xlink:type="simple"/></disp-formula><p>Now, combining Equations (31)-(33), we obtain a equation of the form</p><disp-formula id="scirp.40187-formula29779"><label>(34)</label><graphic position="anchor" xlink:href="12-7501410\2b0e6140-0583-4e29-bd58-3aed22056f29.jpg"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.40187-formula29780"><label>(35)</label><graphic position="anchor" xlink:href="12-7501410\87cf7672-f96c-435b-a97f-cb5b8bd0a8f9.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29781"><label>(36)</label><graphic position="anchor" xlink:href="12-7501410\773e1f34-32ae-4ac2-b72e-25a0d2f2913a.jpg"  xlink:type="simple"/></disp-formula><p>And <img src="12-7501410\d41d5344-6758-41f0-95ce-2cc15ce0c59f.jpg" /> is given in Equation (21). Equation (34) is known as Gardner equation. It is important to note that if we neglect <img src="12-7501410\8c6ae5f2-ec0b-42e0-b195-3a3c00471b8c.jpg" /> term and put<img src="12-7501410\ee10ab48-70f1-4e2b-85eb-8dcae6895dcf.jpg" />, the Gardner equation reduces to K-dV equation which has derived in Equation (18). However, in this K-dV equation the nonlinear term vanishes at<img src="12-7501410\acca91f2-2968-4ce4-a9c2-2abbb7f7c13f.jpg" />, and is not valid near <img src="12-7501410\dd71684b-c1b2-4b42-852a-b3b4f390fccd.jpg" /> which makes soliton amplitude large enough to break down its validity. But the Gardner equation derived here is valid for <img src="12-7501410\84fab0e3-7f34-417a-b660-1e917fd98710.jpg" /> near its critical value.</p></sec></sec><sec id="s4"><title>4. SW Solution of the Gardner Equation</title><p>To analyze stationary GSs, we first introduce a transformation <img src="12-7501410\16a11b3a-04b8-4880-a571-1ef0860da455.jpg" /> which allows us to write Equation (34), under the steady state condition, as</p><disp-formula id="scirp.40187-formula29782"><label>(37)</label><graphic position="anchor" xlink:href="12-7501410\b1ec0d05-d9b1-4086-8ffc-a2cdc9313aeb.jpg"  xlink:type="simple"/></disp-formula><p>where the pseudo-potential <img src="12-7501410\62952f0d-5ad9-49dc-b251-4e811dc4b894.jpg" /> is</p><disp-formula id="scirp.40187-formula29783"><label>(38)</label><graphic position="anchor" xlink:href="12-7501410\7596c9ad-da79-4c34-aadc-9b5630a8adb7.jpg"  xlink:type="simple"/></disp-formula><p>It is obvious from Equation (38) that</p><disp-formula id="scirp.40187-formula29784"><label>(39)</label><graphic position="anchor" xlink:href="12-7501410\a34fa70e-a845-4df8-8e3a-2eff446ff2a5.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29785"><label>(40)</label><graphic position="anchor" xlink:href="12-7501410\5d79e299-80f6-46ae-8c9b-5c780bfb6c10.jpg"  xlink:type="simple"/></disp-formula><p>The conditions of Equations (39) and (40) imply that SW solutions of (37) exist if</p><disp-formula id="scirp.40187-formula29786"><label>(41)</label><graphic position="anchor" xlink:href="12-7501410\15dab1f5-ce88-4778-aee4-a5821ddb87b2.jpg"  xlink:type="simple"/></disp-formula><p>The latter can be solved as</p><disp-formula id="scirp.40187-formula29787"><label>(42)</label><graphic position="anchor" xlink:href="12-7501410\3759e820-146d-4aa2-954d-65c4cc4e6452.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.40187-formula29788"><label>(43)</label><graphic position="anchor" xlink:href="12-7501410\7e0eae63-782b-43eb-b9b2-8dbfa8d8f083.jpg"  xlink:type="simple"/></disp-formula><p>where<img src="12-7501410\0b14bb6a-37ea-40a4-8d25-31e10f32b318.jpg" />, and<img src="12-7501410\e00bb94e-1a42-45ce-aa72-994ca6015e47.jpg" />. Now, using Equations (38) and (43) in Equation (37) we have</p><disp-formula id="scirp.40187-formula29789"><label>(44)</label><graphic position="anchor" xlink:href="12-7501410\b79a00bf-105a-45ac-a9fe-0ef2262faa44.jpg"  xlink:type="simple"/></disp-formula><p>where<img src="12-7501410\a846be98-662c-4d43-a715-08ac89bfd2e0.jpg" />. The SW solution of Equations (37) or (44) is, therefore, directly given by</p><disp-formula id="scirp.40187-formula29790"><label>(45)</label><graphic position="anchor" xlink:href="12-7501410\bffddbad-ef8a-419e-a7df-742d8be374af.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="12-7501410\4dfd3778-6cdd-42d0-85c8-310c563736c5.jpg" /> are given in Equation (43) and SWs width <img src="12-7501410\1f8b9697-39f5-419e-b1e6-51956d8a2865.jpg" /> is</p><disp-formula id="scirp.40187-formula29791"><label>(46)</label><graphic position="anchor" xlink:href="12-7501410\014c6ec7-37f4-49c4-8a87-d58606af9e83.jpg"  xlink:type="simple"/></disp-formula><p>Figures 1-4 show the variation of amplitude of positive (negative) GSs with μ for U<sub>0</sub> = 0.5 and H = 0.3. These figures clearly indicate that both positive and negative GSs exist around the crical value,<img src="12-7501410\f0bbb1a3-23d6-4584-9ed1-9a77fc854224.jpg" />. It has been found that the amplitude (magnitude of the amplitude) of both positive and negative GSs decrease with the increase of μ. Figures 2-5 represent the variation of amplitude of positive (negative) GSs with U<sub>0</sub> for μ = 0.66 (μ = 0.67) and H = 0.3. These figures indicate that the amplitude of both positive and negative GSs increase with the increase of U<sub>0</sub>. We have found that the amplitude of positive and negative GSs does not vary with the quantum diffraction parameter, H but the width of the both positive and negative GSs vary with it. Figures 3-6 imply that the width of both positive and negative GSs decrease with the increase of H and increase with the increase of μ. We have also noticed that in our present system the GSs exist when the quantum effect of electron is neglected.</p></sec><sec id="s5"><title>5. Discussion</title><p>We have investigated QDIA GSs in quantum dusty plasma by deriving Gardner equation. The K-dV solitons are not valid for <img src="12-7501410\6bba6fd1-8709-4b2a-893b-953e8a7fa392.jpg" /> and<img src="12-7501410\ff054436-9b37-4f6b-aa94-66e0df795ee1.jpg" />, which vanish</p><p>at the nonlinear coefficients of the K-dV equation. However, the QDIA GSs investigated in our present work are valid for<img src="12-7501410\7a732447-03ad-40b3-a292-d34102a44122.jpg" />. The results, which have been obtained from this investigation, can be summarized as follows:</p><p>1) The quantum dusty plasma system under consideration supports finite amplitude GSs, whose basic features (polarity, amplitude, width, etc.) depend on the ion and dust number densities and quantum diffraction (tunneling) parameter, H.</p><p>2) GSs are shown to exist for<img src="12-7501410\201505e2-34c3-4c2d-9624-4eb47de268de.jpg" />, and are found to be different from K-dV solitons, which do not exist for<img src="12-7501410\3daffa27-f9db-4950-900e-2b84ad38e95a.jpg" />.</p><p>3) It has been found that at<img src="12-7501410\f056d128-6f3e-4899-a59b-881e334afa83.jpg" />, positive GSs exist, whereas at<img src="12-7501410\94f32907-e9fe-4218-88ac-294bcc20d734.jpg" />, negative GSs exist.</p><p>4) We have seen that the amplitude of positive and negative GSs decreases with μ, increases with<img src="12-7501410\65ccca4e-cf2a-4c68-af92-d45862124271.jpg" />, and does not depend on H.</p><p>5) We have also observed that the width of the GSs increases with μ but decreases with the increase of H.</p><p>It should be mentioned here that in our present investigation, we have neglected the quantum effect of ions since ions are heavier than electrons. However, QDIA solitary waves in quantum dusty plasma with or without the effects of obliqueness and external magnetic field are also problems of recent interest for many space and laboratory dusty plasma situations, but beyond the scope of our present investigation. In conclusion, we propose that a new experiment may be designed based on our results to observe such waves in both laboratory and space quantum dusty plasma system.</p></sec><sec id="s6"><title>6. Acknowledgements</title><p>The Third World Academy of Science (TWAS) Research Grant for research equipment is gratefully acknowledged.</p></sec><sec id="s7"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.40187-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">M. M. Hossain, A. A. Mamun and K. S. Ashrafi, Physics of Plasmas, Vol. 18, 2011, Article ID: 103704. http://dx.doi.org/10.1063/1.3646738</mixed-citation></ref><ref id="scirp.40187-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">A. A. Mamun and P. K. Shukla, Physics Letters A, Vol. 290, 2001, pp. 173-175. http://dx.doi.org/10.1016/S0375-9601(01)00669-7</mixed-citation></ref><ref id="scirp.40187-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Y. D. Jung, Physics of Plasmas, Vol. 8, 2001, p. 3842. http://dx.doi.org/10.1063/1.1386430</mixed-citation></ref><ref id="scirp.40187-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">G. Chabier, F. Dquchin and Y. Potekhin, Journal of Physics: Condensed Matter, Vol. 14, 2002, p. 133.</mixed-citation></ref><ref id="scirp.40187-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">M. Opher, L. O. Silvia, D. E. Danger, V. K. Decyk and J. M. Dawson, Physics of Plasmas, Vol. 8, 2001, p. 2454. http://dx.doi.org/10.1063/1.1362533</mixed-citation></ref><ref id="scirp.40187-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">G. Manfredi, Fields Institute Communications, Plasmas, Vol. 46, 2005, pp. 263-287.</mixed-citation></ref><ref id="scirp.40187-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">A. Markowich, C. Ringhofer and C. Schmeiser, “Semiconductor Equations,” Springer, Vienna, 1990. http://dx.doi.org/10.1007/978-3-7091-6961-2</mixed-citation></ref><ref id="scirp.40187-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">D. Kremp, Th. Bornath, M. Bonitz and M. Shlanges, Physical Review E, Vol. 60, 1999, pp. 4725-4732. http://dx.doi.org/10.1103/PhysRevE.60.4725</mixed-citation></ref><ref id="scirp.40187-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">G. Manfredi and M. Feix, Physical Review E, Vol. 53, 1996, pp. 6460-6470. http://dx.doi.org/10.1103/PhysRevE.53.6460</mixed-citation></ref><ref id="scirp.40187-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">G. Manfredi and F. Haas, Physical Review B, Vol. 64, 2001, Article ID: 075316. http://dx.doi.org/10.1103/PhysRevB.64.075316</mixed-citation></ref><ref id="scirp.40187-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">D. Anderson, B. Hall, M. Lisak and M. Marklund, Physical Review E, Vol. 65, 2002, Article ID: 046417. http://dx.doi.org/10.1103/PhysRevE.65.046417</mixed-citation></ref><ref id="scirp.40187-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">P. K. Shukla and S. Ali, Physics of Plasmas, Vol. 13, 2006, Article ID: 112111. http://dx.doi.org/10.1063/1.2390688</mixed-citation></ref><ref id="scirp.40187-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">M. Marklund, Physics of Plasmas, Vol. 12, 2005, Article ID: 082110. http://dx.doi.org/10.1063/1.2012147</mixed-citation></ref><ref id="scirp.40187-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">L. G. Garcia, F. Haas, L. P. L. de Oliveira and J. Goedert, Physics of Plasmas, Vol. 12, 2005, Article ID: 012302. http://dx.doi.org/10.1063/1.1819935</mixed-citation></ref><ref id="scirp.40187-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">S. Mola, G. Manfredi and M. R. Feix, Journal of Physics of Plasmas, Vol. 50, 1993, pp. 145-162. http://dx.doi.org/10.1017/S0022377800026969</mixed-citation></ref><ref id="scirp.40187-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">F. Haas, L. G. Garcia, J. Goedert and G. Manfredi, Physics of Plasmas, Vol. 10, 2003, p. 3858. http://dx.doi.org/10.1063/1.1609446</mixed-citation></ref><ref id="scirp.40187-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">N. Suh, M. R. Feix and P. Bertrand, Journal of Computational Physics, Vol. 94, 1991, pp. 403-418. http://dx.doi.org/10.1016/0021-9991(91)90227-C</mixed-citation></ref><ref id="scirp.40187-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">F. Haas, Physics of Plasmas, Vol. 12, 1993, Article ID: 062117. http://dx.doi.org/10.1063/1.1939947</mixed-citation></ref><ref id="scirp.40187-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">P. K. Shukla and S. Ali, Physics of Plasmas, Vol. 12, 2005, Article ID: 114502. http://dx.doi.org/10.1063/1.2136376</mixed-citation></ref><ref id="scirp.40187-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">P. K. Shukla, Physics Letters A, Vol. 352, 2006, pp. 242-243. http://dx.doi.org/10.1016/j.physleta.2005.11.065</mixed-citation></ref><ref id="scirp.40187-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">L. Stenflo, P. K. Shukla and M. Marklund, Europhysics Letters, Vol. 74, 2006, p. 844.http://dx.doi.org/10.1209/epl/i2006-10032-x</mixed-citation></ref><ref id="scirp.40187-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">P. K. Shukla and V. P. Silin, Physica Scripta, Vol. 45, 2002, p. 508. http://dx.doi.org/10.1088/0031-8949/45/5/015</mixed-citation></ref><ref id="scirp.40187-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">A. Barkan, N. D’Angelo and R. L. Merlino, Planetary and Space Science, Vol. 44, 1996, pp. 239-242. http://dx.doi.org/10.1016/0032-0633(95)00109-3</mixed-citation></ref><ref id="scirp.40187-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">R. L. Merlino, A. Barkan, C. Thompson and N. D’Angelo, Physics of Plasmas, Vol. 5, 1998, p. 1607. http://dx.doi.org/10.1063/1.872828</mixed-citation></ref><ref id="scirp.40187-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">K. E. Lonngren, Physics of Plasmas, Vol. 25, 1983, p. 943. http://dx.doi.org/10.1088/0032-1028/25/9/001</mixed-citation></ref><ref id="scirp.40187-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">W. Masood, A. Mushtaq and R. Khan, Physics of Plasmas, Vol. 14, 2007, Article ID: 123702. http://dx.doi.org/10.1063/1.2803775</mixed-citation></ref><ref id="scirp.40187-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">S. A. Khan and A. Mushtaq, Physics of Plasmas, Vol. 14, 2007, Article ID: 083703. http://dx.doi.org/10.1063/1.2756752</mixed-citation></ref><ref id="scirp.40187-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">W. Masood, M. Siddiq, S. Nargis and A. M. Mirza, Physics of Plasmas, Vol. 16, 2009, Article ID: 013705. http://dx.doi.org/10.1063/1.3068171</mixed-citation></ref><ref id="scirp.40187-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">C. K. Goertz, Reviews of Geophysics, Vol. 27, 1989, pp. 271-292. http://dx.doi.org/10.1029/RG027i002p00271</mixed-citation></ref><ref id="scirp.40187-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">G. E. Morfill and A. V. Ivlev, Reviews of Modern Physics, Vol. 81, 2009, pp. 1353-1404. http://dx.doi.org/10.1103/RevModPhys.81.1353</mixed-citation></ref><ref id="scirp.40187-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">R. L. Merlino and J. Goree, Physics Today, Vol. 57, 2004, p. 32. http://dx.doi.org/10.1063/1.1784300</mixed-citation></ref><ref id="scirp.40187-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">P. K. Shukla and A. A. Mamun, Introduction to Dusty Plasma Physics, IoP Publishing Ltd., Bristol, 2002, 8-27. http://dx.doi.org/10.1887/075030653X</mixed-citation></ref><ref id="scirp.40187-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">M. Asaduzzaman, “Electro-Acoustic Waves in Presence of Polarization Force,” LAP LAMBERT Academic Publishing, Saarbrücken, 2012.</mixed-citation></ref><ref id="scirp.40187-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">M. Asaduzzaman and A. A. Mamun, Physical Review E, Vol. 86, 2012, Article ID: 016409. http://dx.doi.org/10.1103/PhysRevE.86.016409</mixed-citation></ref><ref id="scirp.40187-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">S. A. Khan, S. Mahmood and A. M. Mirza, Physics Letters A, Vol. 372, 2009, pp. 148-153. http://dx.doi.org/10.1016/j.physleta.2007.10.062</mixed-citation></ref><ref id="scirp.40187-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">M. Asaduzzaman and A. A. Mamun, Astrophysics and Space Science, Vol. 341, 2012, pp. 535-542. http://dx.doi.org/10.1007/s10509-012-1102-x</mixed-citation></ref><ref id="scirp.40187-ref37"><label>37</label><mixed-citation publication-type="other" xlink:type="simple">N. C. Lee, Physics of Plasmas, Vol. 16, 2009, Article ID: 042316. http://dx.doi.org/10.1063/1.3121242</mixed-citation></ref><ref id="scirp.40187-ref38"><label>38</label><mixed-citation publication-type="other" xlink:type="simple">M. M. Hossain and A. A. Mamun, Journal of Physics A: Mathematical and Theoretical, Vol. 45, 2012, Article ID: 125501. http://dx.doi.org/10.1088/1751-8113/45/12/125501</mixed-citation></ref><ref id="scirp.40187-ref39"><label>39</label><mixed-citation publication-type="other" xlink:type="simple">M. Hasan, M. M. Hossain and A. A. Mamun, Astrophysics and Space Science, Vol. 345, 2013, pp. 113-118. http://dx.doi.org/10.1007/s10509-013-1384-7</mixed-citation></ref><ref id="scirp.40187-ref40"><label>40</label><mixed-citation publication-type="other" xlink:type="simple">A. Mushtaq and S. A. Khan, Physics of Plasmas, Vol. 14, 2007, Article ID: 052307. http://dx.doi.org/10.1063/1.2727474</mixed-citation></ref><ref id="scirp.40187-ref41"><label>41</label><mixed-citation publication-type="other" xlink:type="simple">L. D. Landau and E. M. Lifshitz, “Statistical Physics, Pt. 1,” Butterworth-Heinemann, Oxford, 1980.</mixed-citation></ref></ref-list></back></article>