<?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">AM</journal-id><journal-title-group><journal-title>Applied Mathematics</journal-title></journal-title-group><issn pub-type="epub">2152-7385</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/am.2014.53035</article-id><article-id pub-id-type="publisher-id">AM-42643</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>
 
 
  &lt;i&gt;G&lt;/i&gt;-Function Solutions for Schr&#246;dinger Equation in Cylindrical Coordinates System
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>mir</surname><given-names>Pishkoo</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>Maslina</surname><given-names>Darus</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>School of Mathematical Sciences, Faculty of Science and Technology,
Universiti Kebangsaan Malaysia, Selangor, Malaysia</addr-line></aff><aff id="aff1"><addr-line>Nuclear Science Research School, NSTRI, Tehran, Iran;School of Mathematical Sciences, Faculty of Science and Technology,
Universiti Kebangsaan Malaysia, Selangor, Malaysia</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>apishkoo@gmail.com(MP)</email>;<email>maslina@ukm.my(MD)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>07</day><month>02</month><year>2014</year></pub-date><volume>05</volume><issue>03</issue><fpage>342</fpage><lpage>346</lpage><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>
 
 
   In this paper, the Schr?dinger equation is solved by Modified separation of variables (MSV) method suggested by Pishkoo and Darus. Using this method, Meijer’s G-function solutions are derived in cylindrical coordinate system for quantum particle in cylindrical can. All elementary functions and most of the special functions which are the solution of extensive problems in physics and engineering are special cases of Meijer’s G-functions. 
 
</p></abstract><kwd-group><kwd>Meijer’s &lt;i&gt;G&lt;/i&gt;-Function; Partial Differential Equation; Modified Separation of Variables; Schr&#246;dinger Equation</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Perhaps, Cornelis Simon Meijer (1936) was the first to introduce the G-function in mathematics. This very general function intended to include most of the known special functions as particular cases; and for a long time, many studies have been done involving this type of functions. These functions have a lot of practical applications in the fields of mathematical physics, theoretical physics, mathematical analysis, etc.</p><p>Meijer’s G-functions are defined as Mellin-Barnes contour integrals which have been in existence for over 60 years [1-5]. Meijer’s G-function satisfies the linear ordinary differential equation (LODE) of the generalised hypergeometric type whose order is equal to <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\096bd9d7-4c2c-4cda-9159-8f845df0eac6.png" xlink:type="simple"/></inline-formula> [6-8]. This fact triggered us to verify the equality conditions between Meijer’s G-function’s LODE and some partial differential equations governing physical phenomena. In physics, we have many ordinary and partial differential equations, in which their solutions are elementary functions, special functions or a combination of both of them. Thus, Meijer’s G-functions can be the solution for many physical problems if the equality requirement between Meijer’s G-function’s LODE and those differential equations are verified. As such, we seek to deduce the solution of physical problems explicitly in terms of Meijer’s G-functions.</p><p>Our previous works had focused on the introduction of the Modified separation of variables method (MSV), and applying it to solve partial differential equation related to the Reaction-Diffusion process [<xref ref-type="bibr" rid="scirp.42643-ref9">9</xref>], Laplace’s diffusion and Schr&#246;dinger equations [10-12] which led to representing its solution in terms of Meijer’s Gfunctions. The Cartesian coordinate system is used to derive their solutions. However, in this paper, we obtain G-function solutions for the same problem related to Schr&#246;dinger equation, solved by “separation of variables (SV)” in [<xref ref-type="bibr" rid="scirp.42643-ref13">13</xref>], by using modified separation of variables (MSV) method and cylindrical coordinates system as follows.</p></sec><sec id="s2"><title>2. Meijer’s G-Function</title><p>We begin with the definition of Meijer’s G-function as the following:</p><p>Definition 1 A definition of the Meijer’s G-function is given by the following path integral in the complex plane, called Mellin-Barnes type integral:</p><disp-formula id="scirp.42643-formula84343"><label>(1.1)</label><graphic position="anchor" xlink:href="htmlimages\3-7401949x\1cb66f01-3a52-45c3-a4a8-2351284d6768.png"  xlink:type="simple"/></disp-formula><p>Here, the integers <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\4cfdaf62-a0d5-45f5-bff9-3a9d139b8bc8.png" xlink:type="simple"/></inline-formula> are called “orders” of the G-function, or the components of the order<inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\c27ba49b-5a3a-4247-a3d7-e08a848f97fd.png" xlink:type="simple"/></inline-formula>; <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\a5bce936-2fe5-4c2f-be98-0e5a8cdce320.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\a3a78034-7b78-4aae-bdd4-2348e8fa7fab.png" xlink:type="simple"/></inline-formula> are called “parameters” and in general, they are complex numbers. The definition holds under the following assumptions: <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\97fcbc2c-513b-4800-b298-7844e542bd06.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\a8a18ec2-8082-4e05-b1a4-2f7e0d6c32d8.png" xlink:type="simple"/></inline-formula>, where <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\4244c503-f7c6-4b90-9e4d-21d4558394d0.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\81e860a1-f000-4351-ac24-5a2427a99870.png" xlink:type="simple"/></inline-formula> are integer numbers. Subtracting parameters <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\252be12b-4523-42cd-966b-a88c760b9f12.png" xlink:type="simple"/></inline-formula> for <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\0d85d1ae-ac7c-4adc-b266-7dc517ecda22.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\0853e97d-8d0f-4804-8b90-0e839ee2dc2b.png" xlink:type="simple"/></inline-formula> imply that no pole of any <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\e70efda2-524a-4189-9a3c-128ce5637b54.png" xlink:type="simple"/></inline-formula> coincides with any pole of any<inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\b8881e89-4aee-4a44-9814-a94374d8d447.png" xlink:type="simple"/></inline-formula>.</p><p>Exercise 1 Using (1.1), wo obtain the follows <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\7ab22048-b11e-49bf-bab7-37b1be12f0e9.png" xlink:type="simple"/></inline-formula> see [<xref ref-type="bibr" rid="scirp.42643-ref14">14</xref>]</p><disp-formula id="scirp.42643-formula84344"><label>(1.2)</label><graphic position="anchor" xlink:href="htmlimages\3-7401949x\4a77409e-4ed4-4415-9f69-31942688d4a8.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.42643-formula84345"><label>(1.3)</label><graphic position="anchor" xlink:href="htmlimages\3-7401949x\fecf1ed2-137e-4ce9-a22f-04a134feb52f.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.42643-formula84346"><label>(1.4)</label><graphic position="anchor" xlink:href="htmlimages\3-7401949x\2a124691-575d-436d-a7f4-ba356fe508ba.png"  xlink:type="simple"/></disp-formula><p>The Meijer’s G-function <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\9cc9150e-8044-43de-82c1-42dec9761770.png" xlink:type="simple"/></inline-formula> satisfies the linear ordinary differential equation of the generalised hypergeometric type</p><disp-formula id="scirp.42643-formula84347"><label>(1.5)</label><graphic position="anchor" xlink:href="htmlimages\3-7401949x\d94ff23f-7c83-4223-aa40-8ae4a018e64e.png"  xlink:type="simple"/></disp-formula><p>whose order is equal to<inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\03294654-bdbe-4f77-846e-782a5d640b86.png" xlink:type="simple"/></inline-formula>.</p><p>Choosing appropriate values for <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\fb76e155-8c77-4fc9-b7d3-0b0b5d13f437.png" xlink:type="simple"/></inline-formula> orders of G-functions, Equation (1.5) can be changed to complex first and second order linear differential equations. The following section discusses the properties of the solutions of complex first and second order differential equations, and then studies the properties of coefficient functions of these differential equations in the complex plane.</p><p>We start with using Modified separation of variables method (MSV) in cylindrical coordinates system as follows:</p></sec><sec id="s3"><title>3. The G-Function Solutions for the Schr&#246;dinger Equation</title><p>When the geometry of the boundaries is cylindrical, the appropriate coordinate system is the cylindrical one. Separation of variables leads to ODEs in which certain constants (eigenvalues) appear. Different choices of signs for these constants can lead to different functional forms of the general solution. Thus general form of the solution is indeterminate. However, once the boundary conditions are imposed, the unique solutions will emerge regardless of the initial functional form of the solutions. Writing <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\1fd0549d-2886-4359-88b3-3e2550caf8bf.png" xlink:type="simple"/></inline-formula> as a product of three functions, <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\2a569808-7ef4-4226-bdf0-b3a4a33d5eae.png" xlink:type="simple"/></inline-formula>and a separation of variables transforms Laplace’s equation, <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\f42fe82f-9846-4f94-bd3c-1b172e5cd198.png" xlink:type="simple"/></inline-formula>, into three ODEs:</p><p><img src="htmlimages\3-7401949x\d30759d2-edec-4bc0-94e8-962952ac8cf8.png" /></p><p><img src="htmlimages\3-7401949x\fb75f9be-4369-4dfe-a433-b2c53c5618fb.png" /></p><p>where in anticipation of the correct BCs, we have written the constants as <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\fc38fa37-d5f9-409e-a1f6-eed171949e38.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\980feaa0-7ee4-4a8f-b2eb-f1eb58230f87.png" xlink:type="simple"/></inline-formula> with <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\fcabb91e-4988-465c-84ef-8ee1768a4bb7.png" xlink:type="simple"/></inline-formula> an integer.</p><p><img src="htmlimages\3-7401949x\ee67c9ae-48cc-4cd0-a6f0-29fb835ad099.png" /></p><p>is a combination of two independent solutions deduced by “separation of variables method” or</p><p><img src="htmlimages\3-7401949x\6d4c96ba-3941-4436-a5e0-8da0f259e2fb.png" /></p><p>is a combination of two independent solutions deduced by “modified separation of variables method”. Similarly</p><p><img src="htmlimages\3-7401949x\d70d3540-ee19-4852-aa93-fd1a30cc5eac.png" /></p><p>is a combination of two independent solutions deduced by “separation of variables method” or</p><p><img src="htmlimages\3-7401949x\41d6f5d3-8f7d-47cb-a0d3-fe535c80c622.png" /></p><p>is a combination of two independent solutions deduced by “modified separation of variables method”.</p><p>1) For<inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\5b6a4718-87b9-435b-8cbe-426be473082c.png" xlink:type="simple"/></inline-formula>, (1.5) reduces to</p><p><img src="htmlimages\3-7401949x\8c6251b8-40ca-47df-9c56-a6862a8855dc.png" /></p><p>By changing<inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\6329bb0b-cb4c-42e3-9a5b-36e6144be631.png" xlink:type="simple"/></inline-formula>, we have</p><p><img src="htmlimages\3-7401949x\bf75eb00-7645-4025-8b38-ca4e34dd2b97.png" /></p><p>By multiplying both sides of the equation by −4, we have</p><p><img src="htmlimages\3-7401949x\d6de9d57-47ac-44cc-9877-69b8bc22342b.png" /></p><p>On the other hand, let Bessel equation</p><p><img src="htmlimages\3-7401949x\e9f6b4d8-0e55-4cc4-9093-590250f2bf10.png" /></p><p>The conditions for equivalence or these two differential equations are</p><p><img src="htmlimages\3-7401949x\c64a3cb8-3b74-4223-ac49-60d7412f2185.png" /></p><p>and its solutions is</p><p><img src="htmlimages\3-7401949x\1f01e0bc-f168-411d-b2ed-bd438c5ec02e.png" /></p></sec><sec id="s4"><title>4. Quantum Particle in a Cylindrical Can</title><p>Exercise 2 Considering a quantum particle in a cylindrical can, for an atomic particle of mass μ confined in a cylindrical can of length L and radius a, the relevant Schr&#246;dinger equation is</p><p><img src="htmlimages\3-7401949x\31ebd302-c333-4b4c-96f1-5d91806c7b37.png" /></p><p>Let us solve this equation subject to the BCs that <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\3b16a286-2f2f-46db-be60-d3c7c78b6464.png" xlink:type="simple"/></inline-formula> vanishes at the sides of the can by using MSV method.</p><p>Here is the steps: A separation of variables,</p><p><inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\7a1d04b4-3cc5-4968-bc69-9830b2247e03.png" xlink:type="simple"/></inline-formula>leads to the following ODEs:</p><disp-formula id="scirp.42643-formula84348"><label>(1.6)</label><graphic position="anchor" xlink:href="htmlimages\3-7401949x\ad59d95a-fcd7-42a7-ba21-8bca751a5805.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.42643-formula84349"><label>(1.7)</label><graphic position="anchor" xlink:href="htmlimages\3-7401949x\d11bffac-b02b-4b12-a382-727f7b7a756a.png"  xlink:type="simple"/></disp-formula><p>The second DE, the <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\484abb35-02a6-4af7-8558-467567d2f848.png" xlink:type="simple"/></inline-formula> equation, along with its BCs, constitutes an S-L system with G-function solutions.</p><p><img src="htmlimages\3-7401949x\2906b9f8-f343-4fdc-8c33-28d1dad91e74.png" /></p><p>Since the extra condition of periodicity is usually imposed on the potential for variable<inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\d108e7a8-981f-4df2-8f22-205abc57e20f.png" xlink:type="simple"/></inline-formula>, the third DE, in terms of Meijer’s G-functions, has the general solution.</p><p><img src="htmlimages\3-7401949x\b503e100-32bd-4bd7-b44a-fe4045e753e5.png" /></p><p>If we let<inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\0f657075-e640-4162-8c9a-6d7b29765b8c.png" xlink:type="simple"/></inline-formula>, then the solution of the last equation that is well-behaved is at <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\caa5cf9c-0aa5-4ead-96e7-a7155099a62d.png" xlink:type="simple"/></inline-formula> and is <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\c58d5bb7-3d2a-4294-8a84-6aed7e815919.png" xlink:type="simple"/></inline-formula> From the last BC, R(a) = 0, we obtain the quantisation condition <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\5179a4b1-d186-4b46-add9-5d6906e32c30.png" xlink:type="simple"/></inline-formula> for <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\8ef0c1b7-7a84-41bd-88c9-891a09f8bba3.png" xlink:type="simple"/></inline-formula></p><p>Then, the energy eigenvalues are</p><p><img src="htmlimages\3-7401949x\1e3990d3-4f05-4b06-b773-ab50ef584dea.png" /></p><p><inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\d9702628-4517-4d1e-b901-09d879455726.png" xlink:type="simple"/></inline-formula>component: for<inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\a7c9ba6d-3714-4687-9d04-cd52b313ea8a.png" xlink:type="simple"/></inline-formula>, (1.5) is reduced to:</p><p><img src="htmlimages\3-7401949x\01f6684c-a7d2-43e7-9e0b-e64368923834.png" /></p><p>On the other hand, <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\70af2ce7-b23d-4bb7-8e55-9be0e4d87c1f.png" xlink:type="simple"/></inline-formula>with changing variable<inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\c2a3937d-1014-4943-9db9-aa6e301337e1.png" xlink:type="simple"/></inline-formula>.</p><p>The condition for the equivalence of these two differential equations is given by the solution of t-component. Thus, the general solution can be written as</p><disp-formula id="scirp.42643-formula84350"><label>(1.8)</label><graphic position="anchor" xlink:href="htmlimages\3-7401949x\50abea60-b9c0-494d-8e39-6ef1e2e52a2b.png"  xlink:type="simple"/></disp-formula><p>&amp; <inline-formula><inline-graphic xlink:href="tmlimages\3-7401949x\01c458cf-357a-4143-9009-95453c61dfae.png" xlink:type="simple"/></inline-formula></p></sec><sec id="s5"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.42643-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">V. Kiryakova, “Generalized Fractional Calculus and Applications, Vol. 301,” Chapman &amp; Hall/CRC, London, 1994.</mixed-citation></ref><ref id="scirp.42643-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">L. Andrews, “Special Functions for Engineers and Applied Mathematicians,” Macmillan, 1985.</mixed-citation></ref><ref id="scirp.42643-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Y. Luck, “The Special Functions and Their Approximations, Complex Variables,” Academic Press, New York, 1969.</mixed-citation></ref><ref id="scirp.42643-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">R. Askey, “Meijer G-Function,” Cambridge University Press, Cambridge, 2010.</mixed-citation></ref><ref id="scirp.42643-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">A. Klimyik, “Meijer G-Function,” Springer, Berlin, 2001.</mixed-citation></ref><ref id="scirp.42643-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">A. Erdelyi, W. Magnus, F. Oberhettinger and F. Tricomi, “Higher Transcendental Functions,” McGraw-Hill, New York, 1953.</mixed-citation></ref><ref id="scirp.42643-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">A. Prudnikov, I. Brychkov and O. Marichev, “Integrals and Series: Some More Special Functions,” Gordon &amp; Breach, Philadelphia, 1992.</mixed-citation></ref><ref id="scirp.42643-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">A. Kilbas, H. Srivastava and J. Trujillo, “Theory and Applications of Fractional Differential Equations,” Elsevier, Amsterdam, 2006.</mixed-citation></ref><ref id="scirp.42643-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">A. Pishkoo and M. Darus, “Meijer’s G-Functions (MGFs) in Microand Nanostructures,” Journal of Computational and Theoretical Nanoscience, Vol. 10, No. 10, 2013, pp. 2478-2483. http://dx.doi.org/10.1166/jctn.2013.3234</mixed-citation></ref><ref id="scirp.42643-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">A. Pishkoo and M. Darus, “Some Applications of Meijer G-Functions as Solutions of Differential Equations in Physical Models,” Journal of Mathematical Physics, Analysis, Geometry, Vol. 3, No. 9, 2013, pp. 379-391.</mixed-citation></ref><ref id="scirp.42643-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">A. Pishkoo and M. Darus, “Meijer’s G-Functions as the Solution of Schrodinger Equation,” Proceeding of the 43rd Iranian Mathematics Conference, University of Tabriz, Tabriz, 27-30 August 2012.</mixed-citation></ref><ref id="scirp.42643-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">A. Pishkoo and M. Darus, “G-Function Solutions for Diffusion and Laplace’s Equations,” Journal of Advances in Mathematics, Vol. 4, No. 1, 2013, pp. 359-365.</mixed-citation></ref><ref id="scirp.42643-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">S. Hassani, “Mathematical Physics: A Modern Introduction to Its Foundations,” Springer, New York, 1998.</mixed-citation></ref><ref id="scirp.42643-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">A. Pishkoo and M. Darus, “Fractional Differintegral Transformations of Univalent Meijer’s G-Functions,” Journal of Inequalities and Applications, Vol. 2012, 2012, p. 36. http://dx.doi.org/10.1186/1029-242X-2012-36</mixed-citation></ref></ref-list></back></article>