<?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">JAMP</journal-id><journal-title-group><journal-title>Journal of Applied Mathematics and Physics</journal-title></journal-title-group><issn pub-type="epub">2327-4352</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jamp.2016.42027</article-id><article-id pub-id-type="publisher-id">JAMP-63510</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>
 
 
  Quest for Potentials in the Quintessence Scenario
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Tetsuya</surname><given-names>Hara</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Physics, Kyoto Sangyo University, Kyoto 603-8555, Japan</addr-line></aff><author-notes><corresp id="cor1">* E-mail:</corresp></author-notes><pub-date pub-type="epub"><day>17</day><month>02</month><year>2016</year></pub-date><volume>04</volume><issue>02</issue><fpage>211</fpage><lpage>214</lpage><history><date date-type="received"><day>3</day>	<month>January</month>	<year>2016</year></date><date date-type="rev-recd"><day>accepted</day>	<month>13</month>	<year>February</year>	</date><date date-type="accepted"><day>17</day>	<month>February</month>	<year>2016</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><html>
 <head></head>
 
   The time evolution of the equation of state w for quintessence scenario with a scalar field as dark energy is studied up to the third derivative (d3w/da3) with respect to the scale factor a, in order to predict the future observations and specify the scalar potential parameters with the observables. The third derivative of w for general potential V is derived and applied to several types of potentials. They are the inverse power-law (V = M4 + α/Qα), the exponential  <img src="Edit_d936fade-ff15-4362-8cd9-3813511368f6.bmp" alt="" />, the mixed <img src="Edit_d7329c2f-0573-41da-a4d2-690125603200.bmp" alt="" />
   
    <!--[if gte mso 9]><xml>
 <o:oleobject type="Embed" progid="Equation.DSMT4" shapeid="_x0000_i1026" drawaspect="Content" objectid="_1517233582">
 </o:oleobject>
</xml><![endif]-->, the cosine <img src="Edit_e975299d-36bf-4076-97f2-e5021e7bfa59.bmp" alt="" />
   
    <!--[if gte mso 9]><xml>
 <o:oleobject type="Embed" progid="Equation.DSMT4" shapeid="_x0000_i1027" drawaspect="Content" objectid="_1517233583">
 </o:oleobject>
</xml><![endif]--> and the Gaussian types <img src="Edit_f3c7eaa9-7a68-46db-8b07-0f9919e3ff3c.bmp" alt="" />
   
    <!--[if gte mso 9]><xml>
 <o:oleobject type="Embed" progid="Equation.DSMT4" shapeid="_x0000_i1028" drawaspect="Content" objectid="_1517233584">
 </o:oleobject>
</xml><![endif]-->, which are prototypical potentials for the freezing and thawing models. If the parameter number for a potential form is n, it is necessary to find at least for n + 2 independent observations to identify the potential for0m and the evolution of the scalar field (Q and <img src="Edit_d7ce7053-34ac-40ac-a7ef-e069f8e11949.bmp" alt="" />
   
    <!--[if gte mso 9]><xml>
 <o:oleobject type="Embed" progid="Equation.DSMT4" shapeid="_x0000_i1029" drawaspect="Content" objectid="_1517233585">
 </o:oleobject>
</xml><![endif]-->). Such observations would be the values of ΩQ, w, dw/da, ???, and dwn/dan. From these specific potentials, we can predict the n + 1 and higher derivative of w; dwn + 1/dan + 1, ???. Since four of the above mentioned potentials have two parameters, it is necessary to calculate the third derivative of w for them to estimate the predict values. If they are tested observationally, it will be under- stood whether the dark energy could be described by the scalar field with this potential. At least it will satisfy the necessary conditions. Numerical analysis for d3w/da3 is made under some specified parameters in the investigated potentials, except the mixed one. It becomes possible to distinguish the potentials by the accurate observing dw/da and d<sup>2</sup>w/da<sup>2</sup> in some parameters.  
 
</html></p></abstract><kwd-group><kwd>Dark Energy</kwd><kwd> Quintessence</kwd><kwd> Thawing Model</kwd><kwd> Freezing Model</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>There are mainly two theoretical viewpoints to explain the accelerated universe. One is related to modification of gravity and the other is associated with vacuum energy and/or matter field theories. Taking the latter viewpoint, we investigate the scalar fields in quintessence scenario how relevant it to the dark energy.</p><p>In this scenario, the potential of the scalar field has n independent parameters, so we recognize that in principle n time derivatives of the equation of state with observable Ω<sub>Q</sub> and w are enough to specify the scalar potentials and to predict the higher derivatives. In the paper [<xref ref-type="bibr" rid="scirp.63510-ref1">1</xref>], we have calculated the third derivative of the equation of state for five scalar potentials to identify the models and to predict the future observations. The first and second derivatives have been investigated in the paper [<xref ref-type="bibr" rid="scirp.63510-ref2">2</xref>].</p><p>Usually, the variation of the equation of state w for the dark energy is described by [<xref ref-type="bibr" rid="scirp.63510-ref4">4</xref>]</p><disp-formula id="scirp.63510-formula18"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/63510x9.png"  xlink:type="simple"/></disp-formula><p>where a, w<sub>0</sub>, and w<sub>a</sub> are the scale factor (a = 1 at current), the current value of w(a) and the first derivative of w(a) by w<sub>a</sub> = −dw/da, respectively.</p><p>We have extended the parameters pace, in the paper [<xref ref-type="bibr" rid="scirp.63510-ref1">1</xref>],</p><disp-formula id="scirp.63510-formula19"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/63510x10.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x11.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x12.png" xlink:type="simple"/></inline-formula>.</p><p>Recent Planck and other observations for w(z) are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref> [<xref ref-type="bibr" rid="scirp.63510-ref3">3</xref>]. The typical values for each z are adopted and the derived values of dw/da and d<sup>2</sup>w/da<sup>2</sup> are estimated in <xref ref-type="table" rid="table1">Table 1</xref>.</p><p>We follow the single scalar field formalism of Steinhardt et al. (1999) [<xref ref-type="bibr" rid="scirp.63510-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.63510-ref6">6</xref>] and investigate three potentials for so-called freezing model [<xref ref-type="bibr" rid="scirp.63510-ref7">7</xref>], in which the field is rolling towards down its potential minimum, as V = M<sup>4+</sup><sup>α</sup>/Q<sup>α</sup> (inverse power law) [<xref ref-type="bibr" rid="scirp.63510-ref8">8</xref>], <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x13.png" xlink:type="simple"/></inline-formula>(exponential), and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x14.png" xlink:type="simple"/></inline-formula> (mixed) [<xref ref-type="bibr" rid="scirp.63510-ref9">9</xref>]. In this freezing model, w(z) approaches to −1.</p><p>We study other two potentials for so-called thawing model, in which the field is nearly constant at first and then starts to evolve slowly down the potential; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x15.png" xlink:type="simple"/></inline-formula>(cosine) and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x16.png" xlink:type="simple"/></inline-formula> (Gaussian). In this thawing model, w(z) starts from −1 and increases later.</p><p>Because four of the above mentioned potentials have two parameters, it is necessary to calculate the third derivative of w for them to estimate the predict values. If they are the predicted one, it will be understood that the dark energy could be described by the scalar field with this potential. At least it will satisfy the necessary</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> From the observations, the reconstructed equation of state w(z) as a function of red shift z where 1 + z = 1/a [<xref ref-type="bibr" rid="scirp.63510-ref3">3</xref>]. The simulated curves in <xref ref-type="table" rid="table1">Table 1</xref> are designated by signs</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/63510x17.png"/></fig><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> The values of w(z) for each z are adopted and the values dw/da and dw<sup>2</sup>/da<sup>2</sup> are estimated. The typical values are designated by signs, which are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref> and plotted in <xref ref-type="fig" rid="fig2">Figure 2</xref></title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >z = 0.1</th><th align="center" valign="middle" >z = 0.3</th><th align="center" valign="middle" >z = 0.5</th><th align="center" valign="middle" >dw/da</th><th align="center" valign="middle" >d<sup>2</sup>w/da<sup>2</sup></th><th align="center" valign="middle" >sign</th></tr></thead><tr><td align="center" valign="middle" >w(z)</td><td align="center" valign="middle" >−0.90</td><td align="center" valign="middle" >−0.91</td><td align="center" valign="middle" >−0.85</td><td align="center" valign="middle" >0.45</td><td align="center" valign="middle" >5.42</td><td align="center" valign="middle" >a)</td></tr><tr><td align="center" valign="middle" >w(z)</td><td align="center" valign="middle" >−0.90</td><td align="center" valign="middle" >−0.91</td><td align="center" valign="middle" >−0.90</td><td align="center" valign="middle" >0.17</td><td align="center" valign="middle" >1.39</td><td align="center" valign="middle" >b)</td></tr><tr><td align="center" valign="middle" >w(z)</td><td align="center" valign="middle" >−0.90</td><td align="center" valign="middle" >−0.93</td><td align="center" valign="middle" >−0.95</td><td align="center" valign="middle" >0.23</td><td align="center" valign="middle" >0.16</td><td align="center" valign="middle" >c)</td></tr><tr><td align="center" valign="middle" >w(z)</td><td align="center" valign="middle" >−0.90</td><td align="center" valign="middle" >−0.93</td><td align="center" valign="middle" >−0.958</td><td align="center" valign="middle" >0.18</td><td align="center" valign="middle" >−0.48</td><td align="center" valign="middle" >d)</td></tr></tbody></table></table-wrap><p>conditions. Numerical analysis are made for d<sup>3</sup>w/da<sup>3</sup> under some specified parameters in the investigated potentials except mixed one which has three parameters [<xref ref-type="bibr" rid="scirp.63510-ref1">1</xref>].</p></sec><sec id="s2"><title>2. Equation of state w<sub>Q</sub> by a Scalar Field</title><p>For the dark energy, we consider a scalar field Q(x,t), where the action for this field in the gravitational field is described in [<xref ref-type="bibr" rid="scirp.63510-ref2">2</xref>].</p><p>Neglecting the coordinate dependence, the equations for Q(t) becomes</p><disp-formula id="scirp.63510-formula20"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/63510x18.png"  xlink:type="simple"/></disp-formula><p>where H is the Hubble parameter, over-dot is the derivative with time, and V' is the derivative with Q. The equation of state w<sub>Q</sub> due to the scalar field is described by</p><disp-formula id="scirp.63510-formula21"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/63510x19.png"  xlink:type="simple"/></disp-formula><p>We put w<sub>Q</sub> = −1 + ∆ for the later convenience (0 &lt; ∆ &lt; 0.2).</p></sec><sec id="s3"><title>3. Second, and Third Derivative of w<sub>Q</sub></title><p>The detailed calculations of the second, and third derivatives of w<sub>Q</sub> for potentials are displayed in the paper [<xref ref-type="bibr" rid="scirp.63510-ref1">1</xref>]. The numerical calculations for the freezing and thawing models under limited parameters are analyzed there.</p><p>In <xref ref-type="fig" rid="fig2">Figure 2</xref>, the curve for α = 0 is presented for the case of V = M<sup>4+</sup><sup>α</sup>/Q<sup>α</sup> with ∆ = 0.1 by the red solid curve in the dw<sub>Q</sub>/da and d<sup>2</sup>w<sub>Q</sub>/da<sup>2</sup> coordinates. The signature of α will change beyond the parabolic curve. We assume α &gt; 0, so that the upper part of the red curve is forbidden for this potential and the freezing type potentials as well. The green (inner) dotted curve is the case of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x20.png" xlink:type="simple"/></inline-formula> with ∆ = 0.1. Upper part of the green dotted curve is allowed region for this potential. The allowed region of the other thawing potential (Gaussian) is the upper part of the red curve.</p><p>The interesting point is that the forbidden regions for the freezing type potentials are allowed region for the thawing type potentials and the reverse is also true. It is possible to distinguish the potentials among each type due to the different predicted values of d<sup>3</sup>w<sub>Q</sub>/da<sup>3</sup> [<xref ref-type="bibr" rid="scirp.63510-ref1">1</xref>], however it is necessary to make accurate observations for the values of dw<sub>Q</sub>/da, d<sup>2</sup>w<sub>Q</sub>/da<sup>2</sup>, d<sup>3</sup>w<sub>Q</sub>/da<sup>3</sup> and other parameters such as ∆.</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> The curve for α = 0 is presented for the case of V =M<sup>4+</sup><sup>α</sup>/Q<sup>α</sup> with ∆ = 0.1 by the red solid curve in the dwQ/da and d<sup>2</sup>wQ/da<sup>2</sup> coordinates. The signature of α will change beyond the parabolic curve and the upper part is forbidden for the freezing model. The lower part of this curve is forbidden for the thawing model (Gausssian type). The upper part of the green (inner) dotted curve is allowed region for thawing model of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x22.png" xlink:type="simple"/></inline-formula> with ∆ = 0.1. The typical values adopted in <xref ref-type="table" rid="table1">Table 1</xref> are plotted. Notice that the values of d<sup>2</sup>wQ/da<sup>2</sup> for a) and b) are out of frame</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/63510x21.png"/></fig></sec><sec id="s4"><title>4. Conclusions</title><p>At present, backward observations, such as Planck, baryon acoustic oscillation, Supernova Ia, Hubble constant, weak lensing, and red shift distortion, have been undertaken to estimate w<sub>Q</sub> at the age (1 + z) as in <xref ref-type="fig" rid="fig1">Figure 1</xref> [<xref ref-type="bibr" rid="scirp.63510-ref3">3</xref>]. From <xref ref-type="fig" rid="fig1">Figure 1</xref>, the rough values of w(a = a<sub>0</sub>), dw/da and d<sup>2</sup>w/da<sup>2</sup> have been estimated which are presented in <xref ref-type="table" rid="table1">Table 1</xref>. They are pointed in the dw<sub>Q</sub>/da and d<sup>2</sup>w<sub>Q</sub>/da<sup>2</sup> plane in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><p>The adopted values from observation show −0.5 &lt; d<sup>2</sup>w<sub>Q</sub>/da<sup>2</sup>&lt; 6 within the region 0.1 &lt; dw<sub>Q</sub>/da &lt; 0.5. Although there is a lot of uncertainty, at the moment, it seems to be preferable for the thawing model against the freezing model under the comparison with the numerical results and the observations [<xref ref-type="bibr" rid="scirp.63510-ref3">3</xref>].</p><p>About observations in <xref ref-type="fig" rid="fig1">Figure 1</xref>, it seems to be difficult to accept that the equation of state w = p/ρ is almost smaller than −1 in the region z −1 ~ 2.</p><p>Usually matter density increases as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x23.png" xlink:type="simple"/></inline-formula>, then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x24.png" xlink:type="simple"/></inline-formula> must increase with z where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x25.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/63510x26.png" xlink:type="simple"/></inline-formula>, taking that p<sub>Q</sub> and p<sub>m</sub> are pressure for scalar field and matter. There seems to be no such features that w increases with z in the observations in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>If w &lt; −1 which means ∆ &lt; 0 is correct in <xref ref-type="fig" rid="fig1">Figure 1</xref>, we must consider fully different models such as phantom, k-essence, chameleon, tachyon, dilaton, quintom, modified gravity, and so on [<xref ref-type="bibr" rid="scirp.63510-ref10">10</xref>].</p></sec><sec id="s5"><title>Cite this paper</title><p>Tetsuya Hara, (2016) Quest for Potentials in the Quintessence Scenario. Journal of Applied Mathematics and Physics,04,211-214. doi: 10.4236/jamp.2016.42027</p></sec></body><back><ref-list><title>References</title><ref id="scirp.63510-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Muromachi, Y., Okabayashi, A., Okada, D., Hara, T. and Itoh, Y. (2015) Search for Dark Energy Potentials in Quintessence. Progress of Theoretical &amp; Experimental Physics, 2015, 093E01. [arXiv:1503.03678 [astro-ph]]</mixed-citation></ref><ref id="scirp.63510-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Hara, T., Sakata, R., Muromachi, Y. and Itoh, Y. (2014) Time Variation of Equation of State for Dark Energy. Progress of Theoretical &amp; Experimental Physics, 2014, 113E01. [arXiv:1409.2726[astro-ph]]</mixed-citation></ref><ref id="scirp.63510-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Planck Collaborations: Ade, P.A.R., et al. (2015) Planck 2015 Results. XIV. Dark energy and Modified Gravity. [ar-Xiv:1502.01590[astro-ph]]</mixed-citation></ref><ref id="scirp.63510-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Linder, E.V. (2003) Exploring the Expansion History of the Universe. Physical Review Letters, 90, 091301. [arXiv:0208512v1[astro-ph]]</mixed-citation></ref><ref id="scirp.63510-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Steinhardt, P.J., Wang, L. and Zlatev, I. (1999) Cosmological Tracking Solutions. Physical Review D, 59, 123504.  
http://dx.doi.org/10.1103/PhysRevD.59.123504</mixed-citation></ref><ref id="scirp.63510-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Zlatev, I., Wang, L. and Steinhardt, P.J. (1999) Quintessence, Cosmic Coincidence, and the Cosmological Constant. Physical Review Letters, 82, 896. http://dx.doi.org/10.1103/PhysRevLett.82.896</mixed-citation></ref><ref id="scirp.63510-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Caldwell, R.R. and Linder, E.V. (2005) Limits of Quintessence. Physical Review Letters, 95, 141301.  
http://dx.doi.org/10.1103/PhysRevLett.95.141301</mixed-citation></ref><ref id="scirp.63510-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Peebles, P.J.E. and Ratra, B. (1988) The Cosmology with a Time-Variable Cosmological “Constant”. Astrophysical Journal, 325, L17-L20. http://dx.doi.org/10.1086/185100</mixed-citation></ref><ref id="scirp.63510-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Brax, Ph. and Martin, J. (1999) Quintessence and Supergravity. Physics Letters B, 468, 40-45.  
http://dx.doi.org/10.1016/S0370-2693(99)01209-5</mixed-citation></ref><ref id="scirp.63510-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Amendola, L. and Tsujikawa, S. (2010) Dark Energy. Cambridge University Press, Cambridge.  
http://dx.doi.org/10.1017/CBO9780511750823</mixed-citation></ref></ref-list></back></article>