Reheating and entropy perturbations in fibre inflation

  • We study reheating in some one and two field realizations of fibre inflation. We find that reheating begins with a phase of preheating in which long wavelength fluctuation modes are excited. In two field models there is a danger that the parametric amplification of infrared fluctuations in the second scalar field - associated with an entropy mode - might induce an instability of the curvature fluctuations. We show that, at least in the models we consider, the entropy mode has a sufficiently large mass to prevent this instability. Hence, from the point of view of reheating the models we consider are well-behaved.
  • [1] R. Allahverdi, R. Brandenberger, F. Y. Cyr-Racine et al, Ann. Rev. Nucl. Part. Sci., 60: 27 (2010), arXiv:1001.2600 [hep-th doi: 10.1146/annurev.nucl.012809.104511
    [2] M. A. Amin, M. P. Hertzberg, D. I. Kaiser et al, Int. J. Mod. Phys. D, 24: 1530003 (2014), arXiv:1410.3808 [hep-ph
    [3] A. Berera, Phys. Rev. Lett., 75: 3218 (1995), arXiv:astro-ph/9509049 doi: 10.1103/PhysRevLett.75.3218
    [4] L. F. Abbott, E. Farhi, and M. B. Wise, Phys. Lett., 117B: 29 (1982)
    [5] A. D. Dolgov and A. D. Linde, Phys. Lett., 116B: 329 (1982)
    [6] A. Albrecht, P. J. Steinhardt, M. S. Turner et al, Phys. Rev. Lett., 48: 1437 (1982) doi: 10.1103/PhysRevLett.48.1437
    [7] J. H. Traschen and R. H. Brandenberger, Phys. Rev. D, 42: 2491 (1990) doi: 10.1103/PhysRevD.42.2491
    [8] A. D. Dolgov and D. P. Kirilova, Sov. J. Nucl. Phys., 51: 172 (1990)
    [9] L. Kofman, A. D. Linde, and A. A. Starobinsky, Phys. Rev. Lett., 73: 3195 (1994), arXiv:hepth/9405187 doi: 10.1103/PhysRevLett.73.3195
    [10] Y. Shtanov, J. H. Traschen, and R. H. Brandenberger, Phys. Rev. D, 51: 5438 (1995), arXiv:hep-ph/9407247 doi: 10.1103/PhysRevD.51.5438
    [11] L. Kofman, A. D. Linde, and A. A. Starobinsky, Phys. Rev. D, 56: 3258 (1997), arXiv:hep-ph/9704452 doi: 10.1103/PhysRevD.56.3258
    [12] B. A. Bassett, D. I. Kaiser, and R. Maartens, Phys. Lett. B, 455: 84 (1999) [arXiv:hep-ph/9808404]
    [13] B. A. Bassett, F. Tamburini, D. I. Kaiser et al, Nucl. Phys. B, 561: 188 (1999) [arXiv:hep-ph/9901319]
    [14] B. A. Bassett, C. Gordon, R. Maartens et al, Phys. Rev. D, 61: 061302 (2000) [hep-ph/9909482]
    [15] F. Finelli and R. H. Brandenberger, Phys. Rev. Lett., 82: 1362 (1999), arXiv:hep-ph/9809490 doi: 10.1103/PhysRevLett.82.1362
    [16] D. Langlois and F. Vernizzi, Phys. Rev. Lett., 95: 091303 (2005) [astro-ph/0503416]; D. Langlois and F. Vernizzi, Phys. Rev. D, 72: 103501 (2005) [astroph/0509078]
    [17] N. Afshordi and R. H. Brandenberger, Phys. Rev. D, 63: 123505 (2001), arXiv:gr-qc/0011075 doi: 10.1103/PhysRevD.63.123505
    [18] F. Finelli and R. H. Brandenberger, Phys. Rev. D, 62: 083502 (2000), arXiv:hep-ph/0003172 doi: 10.1103/PhysRevD.62.083502
    [19] A. Taruya and Y. Nambu, Phys. Lett. B, 428: 37 (1998) [arXiv:gr-qc/9709035]; B. A. Bassett and F. Viniegra, Phys. Rev. D, 62: 043507 (2000) [hep-ph/9909353]
    [20] H. B. Moghaddam, R. H. Brandenberger, Y. F. Cai et al, arXiv:1409.1784 [astro-ph.CO]
    [21] R. H. Brandenberger, K. Dasgupta, and A. C. Davis, Phys. Rev. D, 78: 083502 (2008), arXiv:0801.3674 [hep-th doi: 10.1103/PhysRevD.78.083502
    [22] E. McDonough, H. Bazrafshan Moghaddam, and R. H. Brandenberger, JCAP, 1605(05): 012 (2016), arXiv:1601.07749 [hep-th
    [23] Y. F. Cai, R. Brandenberger, and X. Zhang, JCAP, 1103: 003 (2011), arXiv:1101.0822 [hep-th
    [24] M. Cicoli, C. P. Burgess, and F. Quevedo, JCAP, 0903: 013 (2009), arXiv:0808.0691 [hep-th
    [25] C. Vafa, hep-th/0509212; T. D. Brennan, F. Carta, and C. Vafa, arXiv:1711.00864 [hep-th]
    [26] G. Obied, H. Ooguri, L. Spodyneiko et al, arXiv:1806.08362 [hep-th]; P. Agrawal, G. Obied, P. J. Steinhardt et al, arXiv:1806.09718 [hep-th]
    [27] S. Antusch, F. Cefala, S. Krippendorf et al, JHEP, 1801: 083 (2018), arXiv:1708.08922 [hep-th
    [28] P. Cabella, A. Di Marco, and G. Pradisi, Phys. Rev. D, 95: 123528 (2017), arXiv:1704.03209 [astro-ph.CO doi: 10.1103/PhysRevD.95.123528
    [29] D. Baumann, arXiv:0907.5424 [hep-th]
    [30] M. B. Green, J. H. Schwarz, and E. Witten, Superstring Theory. Vol. 2: Loop Amplitudes, Anomalies And Phenomenology, (Cambridge, Uk: Univ. Pr., 1987) p. 596 (Cambridge Monographs On Mathematical Physics)
    [31] D. Baumann and L. McAllister, arXiv:1404.2601 [hep-th]
    [32] S. H. S. Alexander, Phys. Rev. D, 65: 023507 (2002), arXiv:hep-th/0105032
    [33] G. R. Dvali and S. H. H. Tye, Phys. Lett. B, 450: 72 (1999), [hep-ph/9812483]
    [34] C. P. Burgess, M. Majumdar, D. Nolte et al, JHEP, 0107: 047 (2001), [hep-th/0105204]
    [35] G. R. Dvali, Q. Shafi, and S. Solganik, hep-th/0105203
    [36] J. J. Blanco-Pillado, C. P. Burgess, J. M. Cline et al, JHEP, 0411: 063 (2004), arXiv:hep-th/0406230
    [37] J. P. Conlon and F. Quevedo, JHEP, 0601: 146 (2006), arXiv:hep-th/0509012
    [38] P. Svrcek and E. Witten, JHEP, 0606: 051 (2006), arXiv:hep-th/0605206
    [39] K. Dasgupta, H. Firouzjahi, and R. Gwyn, JHEP, 0806: 056 (2008), arXiv:0803.3828 [hep-th
    [40] R. Brandenberger, Int. J. Mod. Phys. D, 26(01): 1740002 (2016), arXiv:1601.01918 [hep-th
    [41] H. Kurki-Suonio, R. A. Matzner, J. Centrella et al, Phys. Rev. D, 35: 435 (1987)
    [42] R. H. Brandenberger and J. H. Kung, Phys. Rev. D, 42: 1008 (1990)
    [43] H. Kurki-Suonio, P. Laguna, and R. A. Matzner, Phys. Rev. D, 48: 3611 (1993), [astro-ph/9306009]
    [44] E. Silverstein and A. Westphal, Phys. Rev. D, 78: 106003 (2008), arXiv:0803.3085 [hep-th doi: 10.1103/PhysRevD.78.106003
    [45] V. Balasubramanian, P. Berglund, J. P. Conlon et al, JHEP, 0503: 007 (2005), arXiv:hep-th/0502058
    [46] S. B. Giddings, S. Kachru, and J. Polchinski, Phys. Rev. D, 66: 106006 (2002), arXiv:hep-th/0105097 doi: 10.1103/PhysRevD.66.106006
    [47] M. R. Douglas and S. Kachru, Rev. Mod. Phys, 79: 733 (2007), arXiv:hep-th/0610102 doi: 10.1103/RevModPhys.79.733
    [48] R. H. Brandenberger and C. Vafa, Nucl. Phys. B, 316: 391 (1989) doi: 10.1016/0550-3213(89)90037-0
    [49] R. H. Brandenberger, Class. Quant. Grav., 28: 204005 (2011)
    [50] R. H. Brandenberger, String Gas Cosmology, (String Cosmology, J.Erdmenger (Editor). Wiley, 2009) p.193-230, [arXiv:0808.0746 [hep-th]]
    [51] T. Battefeld and S. Watson, Rev. Mod. Phys., 78: 435 (2006), [arXiv:hep-th/0510022]
    [52] S. P. Patil and R. Brandenberger, Phys. Rev. D, 71: 103522 (2005), [hep-th/0401037]
    [53] S. P. Patil and R. H. Brandenberger, JCAP, 0601: 005 (2006), [hep-th/0502069]
    [54] S. Watson and R. Brandenberger, JCAP, 0311: 008 (2003), [hepth/0307044]
    [55] S. Watson, Phys. Rev. D, 70: 066005 (2004), [hep-th/0404177]
    [56] R. Brandenberger, Y. K. E. Cheung, and S. Watson, JHEP, 0605: 025 (2006), [hep-th/0501032]
    [57] V.F. Mukhanov, H.A. Feldman, and R.H. Brandenberger, Physics Reports, 215: 203 (1992) doi: 10.1016/0370-1573(92)90044-Z
    [58] R. H. Brandenberger, Lect. Notes Phys., 646: 127 (2004), arXiv:hepth/0306071
    [59] C. Gordon, D. Wands, B. A. Bassett et al, Phys. Rev. D, 63: 023506 (2001), arXiv:astro-ph/0009131
    [60] J. Jiang, Q. Liang, Y. F. Cai et al, Astrophys. J., 876: 136 (2019), arXiv:1812.08220 [astroph.CO doi: 10.3847/1538-4357/ab189e
  • [1] R. Allahverdi, R. Brandenberger, F. Y. Cyr-Racine et al, Ann. Rev. Nucl. Part. Sci., 60: 27 (2010), arXiv:1001.2600 [hep-th doi: 10.1146/annurev.nucl.012809.104511
    [2] M. A. Amin, M. P. Hertzberg, D. I. Kaiser et al, Int. J. Mod. Phys. D, 24: 1530003 (2014), arXiv:1410.3808 [hep-ph
    [3] A. Berera, Phys. Rev. Lett., 75: 3218 (1995), arXiv:astro-ph/9509049 doi: 10.1103/PhysRevLett.75.3218
    [4] L. F. Abbott, E. Farhi, and M. B. Wise, Phys. Lett., 117B: 29 (1982)
    [5] A. D. Dolgov and A. D. Linde, Phys. Lett., 116B: 329 (1982)
    [6] A. Albrecht, P. J. Steinhardt, M. S. Turner et al, Phys. Rev. Lett., 48: 1437 (1982) doi: 10.1103/PhysRevLett.48.1437
    [7] J. H. Traschen and R. H. Brandenberger, Phys. Rev. D, 42: 2491 (1990) doi: 10.1103/PhysRevD.42.2491
    [8] A. D. Dolgov and D. P. Kirilova, Sov. J. Nucl. Phys., 51: 172 (1990)
    [9] L. Kofman, A. D. Linde, and A. A. Starobinsky, Phys. Rev. Lett., 73: 3195 (1994), arXiv:hepth/9405187 doi: 10.1103/PhysRevLett.73.3195
    [10] Y. Shtanov, J. H. Traschen, and R. H. Brandenberger, Phys. Rev. D, 51: 5438 (1995), arXiv:hep-ph/9407247 doi: 10.1103/PhysRevD.51.5438
    [11] L. Kofman, A. D. Linde, and A. A. Starobinsky, Phys. Rev. D, 56: 3258 (1997), arXiv:hep-ph/9704452 doi: 10.1103/PhysRevD.56.3258
    [12] B. A. Bassett, D. I. Kaiser, and R. Maartens, Phys. Lett. B, 455: 84 (1999) [arXiv:hep-ph/9808404]
    [13] B. A. Bassett, F. Tamburini, D. I. Kaiser et al, Nucl. Phys. B, 561: 188 (1999) [arXiv:hep-ph/9901319]
    [14] B. A. Bassett, C. Gordon, R. Maartens et al, Phys. Rev. D, 61: 061302 (2000) [hep-ph/9909482]
    [15] F. Finelli and R. H. Brandenberger, Phys. Rev. Lett., 82: 1362 (1999), arXiv:hep-ph/9809490 doi: 10.1103/PhysRevLett.82.1362
    [16] D. Langlois and F. Vernizzi, Phys. Rev. Lett., 95: 091303 (2005) [astro-ph/0503416]; D. Langlois and F. Vernizzi, Phys. Rev. D, 72: 103501 (2005) [astroph/0509078]
    [17] N. Afshordi and R. H. Brandenberger, Phys. Rev. D, 63: 123505 (2001), arXiv:gr-qc/0011075 doi: 10.1103/PhysRevD.63.123505
    [18] F. Finelli and R. H. Brandenberger, Phys. Rev. D, 62: 083502 (2000), arXiv:hep-ph/0003172 doi: 10.1103/PhysRevD.62.083502
    [19] A. Taruya and Y. Nambu, Phys. Lett. B, 428: 37 (1998) [arXiv:gr-qc/9709035]; B. A. Bassett and F. Viniegra, Phys. Rev. D, 62: 043507 (2000) [hep-ph/9909353]
    [20] H. B. Moghaddam, R. H. Brandenberger, Y. F. Cai et al, arXiv:1409.1784 [astro-ph.CO]
    [21] R. H. Brandenberger, K. Dasgupta, and A. C. Davis, Phys. Rev. D, 78: 083502 (2008), arXiv:0801.3674 [hep-th doi: 10.1103/PhysRevD.78.083502
    [22] E. McDonough, H. Bazrafshan Moghaddam, and R. H. Brandenberger, JCAP, 1605(05): 012 (2016), arXiv:1601.07749 [hep-th
    [23] Y. F. Cai, R. Brandenberger, and X. Zhang, JCAP, 1103: 003 (2011), arXiv:1101.0822 [hep-th
    [24] M. Cicoli, C. P. Burgess, and F. Quevedo, JCAP, 0903: 013 (2009), arXiv:0808.0691 [hep-th
    [25] C. Vafa, hep-th/0509212; T. D. Brennan, F. Carta, and C. Vafa, arXiv:1711.00864 [hep-th]
    [26] G. Obied, H. Ooguri, L. Spodyneiko et al, arXiv:1806.08362 [hep-th]; P. Agrawal, G. Obied, P. J. Steinhardt et al, arXiv:1806.09718 [hep-th]
    [27] S. Antusch, F. Cefala, S. Krippendorf et al, JHEP, 1801: 083 (2018), arXiv:1708.08922 [hep-th
    [28] P. Cabella, A. Di Marco, and G. Pradisi, Phys. Rev. D, 95: 123528 (2017), arXiv:1704.03209 [astro-ph.CO doi: 10.1103/PhysRevD.95.123528
    [29] D. Baumann, arXiv:0907.5424 [hep-th]
    [30] M. B. Green, J. H. Schwarz, and E. Witten, Superstring Theory. Vol. 2: Loop Amplitudes, Anomalies And Phenomenology, (Cambridge, Uk: Univ. Pr., 1987) p. 596 (Cambridge Monographs On Mathematical Physics)
    [31] D. Baumann and L. McAllister, arXiv:1404.2601 [hep-th]
    [32] S. H. S. Alexander, Phys. Rev. D, 65: 023507 (2002), arXiv:hep-th/0105032
    [33] G. R. Dvali and S. H. H. Tye, Phys. Lett. B, 450: 72 (1999), [hep-ph/9812483]
    [34] C. P. Burgess, M. Majumdar, D. Nolte et al, JHEP, 0107: 047 (2001), [hep-th/0105204]
    [35] G. R. Dvali, Q. Shafi, and S. Solganik, hep-th/0105203
    [36] J. J. Blanco-Pillado, C. P. Burgess, J. M. Cline et al, JHEP, 0411: 063 (2004), arXiv:hep-th/0406230
    [37] J. P. Conlon and F. Quevedo, JHEP, 0601: 146 (2006), arXiv:hep-th/0509012
    [38] P. Svrcek and E. Witten, JHEP, 0606: 051 (2006), arXiv:hep-th/0605206
    [39] K. Dasgupta, H. Firouzjahi, and R. Gwyn, JHEP, 0806: 056 (2008), arXiv:0803.3828 [hep-th
    [40] R. Brandenberger, Int. J. Mod. Phys. D, 26(01): 1740002 (2016), arXiv:1601.01918 [hep-th
    [41] H. Kurki-Suonio, R. A. Matzner, J. Centrella et al, Phys. Rev. D, 35: 435 (1987)
    [42] R. H. Brandenberger and J. H. Kung, Phys. Rev. D, 42: 1008 (1990)
    [43] H. Kurki-Suonio, P. Laguna, and R. A. Matzner, Phys. Rev. D, 48: 3611 (1993), [astro-ph/9306009]
    [44] E. Silverstein and A. Westphal, Phys. Rev. D, 78: 106003 (2008), arXiv:0803.3085 [hep-th doi: 10.1103/PhysRevD.78.106003
    [45] V. Balasubramanian, P. Berglund, J. P. Conlon et al, JHEP, 0503: 007 (2005), arXiv:hep-th/0502058
    [46] S. B. Giddings, S. Kachru, and J. Polchinski, Phys. Rev. D, 66: 106006 (2002), arXiv:hep-th/0105097 doi: 10.1103/PhysRevD.66.106006
    [47] M. R. Douglas and S. Kachru, Rev. Mod. Phys, 79: 733 (2007), arXiv:hep-th/0610102 doi: 10.1103/RevModPhys.79.733
    [48] R. H. Brandenberger and C. Vafa, Nucl. Phys. B, 316: 391 (1989) doi: 10.1016/0550-3213(89)90037-0
    [49] R. H. Brandenberger, Class. Quant. Grav., 28: 204005 (2011)
    [50] R. H. Brandenberger, String Gas Cosmology, (String Cosmology, J.Erdmenger (Editor). Wiley, 2009) p.193-230, [arXiv:0808.0746 [hep-th]]
    [51] T. Battefeld and S. Watson, Rev. Mod. Phys., 78: 435 (2006), [arXiv:hep-th/0510022]
    [52] S. P. Patil and R. Brandenberger, Phys. Rev. D, 71: 103522 (2005), [hep-th/0401037]
    [53] S. P. Patil and R. H. Brandenberger, JCAP, 0601: 005 (2006), [hep-th/0502069]
    [54] S. Watson and R. Brandenberger, JCAP, 0311: 008 (2003), [hepth/0307044]
    [55] S. Watson, Phys. Rev. D, 70: 066005 (2004), [hep-th/0404177]
    [56] R. Brandenberger, Y. K. E. Cheung, and S. Watson, JHEP, 0605: 025 (2006), [hep-th/0501032]
    [57] V.F. Mukhanov, H.A. Feldman, and R.H. Brandenberger, Physics Reports, 215: 203 (1992) doi: 10.1016/0370-1573(92)90044-Z
    [58] R. H. Brandenberger, Lect. Notes Phys., 646: 127 (2004), arXiv:hepth/0306071
    [59] C. Gordon, D. Wands, B. A. Bassett et al, Phys. Rev. D, 63: 023506 (2001), arXiv:astro-ph/0009131
    [60] J. Jiang, Q. Liang, Y. F. Cai et al, Astrophys. J., 876: 136 (2019), arXiv:1812.08220 [astroph.CO doi: 10.3847/1538-4357/ab189e
  • 加载中

Cited by

1. Cicoli, M., Conlon, J.P., Maharana, A. et al. String cosmology: From the early universe to today[J]. Physics Reports, 2024. doi: 10.1016/j.physrep.2024.01.002
2. Bhattacharya, S., Dutta, K., Gangopadhyay, M.R. et al. Fibre inflation and precision CMB data[J]. Physical Review D, 2020, 102(12): 123531. doi: 10.1103/PhysRevD.102.123531
3. Bayarsaikhan, B., Koh, S., Tsedenbaljir, E. et al. Constraints on dark energy models from the Horndeski theory[J]. Journal of Cosmology and Astroparticle Physics, 2020, 2020(11): 057. doi: 10.1088/1475-7516/2020/11/057
4. Aragam, V., Paban, S., Rosati, R. Multi-field inflation in high-slope potentials[J]. Journal of Cosmology and Astroparticle Physics, 2020, 2020(4): 022. doi: 10.1088/1475-7516/2020/04/022

Figures(7)

Get Citation
Bao-Min Gu and Robert Brandenberger. Reheating and entropy perturbations in fibre inflation[J]. Chinese Physics C, 2020, 44(1): 015103. doi: 10.1088/1674-1137/44/1/015103
Bao-Min Gu and Robert Brandenberger. Reheating and entropy perturbations in fibre inflation[J]. Chinese Physics C, 2020, 44(1): 015103.  doi: 10.1088/1674-1137/44/1/015103 shu
Milestone
Received: 2019-09-05
Article Metric

Article Views(2696)
PDF Downloads(48)
Cited by(4)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Reheating and entropy perturbations in fibre inflation

    Corresponding author: Bao-Min Gu, gubm@ncu.edu.cn
    Corresponding author: Robert Brandenberger, rhb@physics.mcgill.ca
  • 1. Department of Physics, Nanchang University, Nanchang 330031, China
  • 2. Center for Relativistic Astrophysics and High Energy Physics, Nanchang University, Nanchang 330031, China
  • 3. Department of Physics, McGill University, Montr′eal, QC, H3A 2T8, Canada

Abstract: We study reheating in some one and two field realizations of fibre inflation. We find that reheating begins with a phase of preheating in which long wavelength fluctuation modes are excited. In two field models there is a danger that the parametric amplification of infrared fluctuations in the second scalar field - associated with an entropy mode - might induce an instability of the curvature fluctuations. We show that, at least in the models we consider, the entropy mode has a sufficiently large mass to prevent this instability. Hence, from the point of view of reheating the models we consider are well-behaved.

    HTML

    1.   Introduction
    • The reheating phase (see e.g. [1, 2] for recent reviews) is an integral part of any successfull inflationary model. The inflationary phase leaves behind a state with an exponentially suppressed density of matter particles [3], and a mechanism is required which transforms the energy trapped in the homogeneous inflaton condensate to quanta of the regular matter fields.

      Initially, reheating was studied perturbatively [4-6], but this neglects the coherent nature of the inflaton condensate. It was then realized [7, 8] that the reheating phase may begin with a period of parametric resonance instability called preheating [9-11] during which the oscillations of the inflaton condensate lead to exponential increase in the number density of long wavelength fluctuations of either the inflaton field itself (“self-resonance”) or of fields which are coupled to the inflaton. The resulting state of matter after preheating has a non-thermal distribution, and hence a second stage of reheating is required in order to obtain full kinetic and chemical equilibrium.

      The equation of state of matter becomes radiation-dominated even after an efficient preheating period. Since the detailed predictions of any inflationary model for observations depends on the length of the period after inflation before the onset of radiation-domination it is important in any inflationary model to study the possible presence of parametric resonance.

      One might fear that the exponential increase in infrared matter fluctuations during preheating might lead to an exponential increase of the cosmological fluctuations [12-15]. In single matter field models, however, this does not happen since the curvature fluctuations are conserved on super-Hubble scales, as can be shown even beyond perturbation theory (see e.g. [16, 17]). However, in models in which there are extra light fields in addition to the inflaton, there is the danger than fluctuations in the extra fields will also experience parametric resonance of infrared modes [18, 19]. If the fluctuations of these entropy modes are coupled to the curvature fluctuations, then the exponentially amplified entropy fluctuations could induce exponentially enhanced curvature fluctuations - on scales which are measured today, hence destroying the agreement between theory and observations. Toy models in which such a dangerous amplification of entropy modes occurs were studied in [20]. On the other hand, it was shown that a number of string-motivated inflationary models such as D3-D7 brane inflation [21] and axion monodromy inflation [22] are safe from this potential problem. Although it is not expected to have amplification of primordial curvature perturbation, it is interesting to note that this amplification may suppress the tensor-to-scalar ratio hence can be applied in some particular theories to better fit current observations. A typical example is the two-field matter bounce curvaton scenario [23], in which the scalar modes are amplified through a “kinetic amplification” while the tensor modes are not, providing a mechanism to suppress the tensor-to-scalar ratio and new sources of non-Gaussianity.

      Fibre inflation [24] is a popular model of inflation motivated by ideas from string theory [25, 26]. In this paper we show that in the versions of the scenario which we consider here, reheating begins with a period of preheating (the possibility of preheating in fibre inflation was also considered in [27], and the macro reheating properties in [28]). On the other hand, we show that the entropy modes are sufficiently heavy such that no efficient resonance of these modes occurs during the phase of preheating. Hence, it appears that from the point of view of reheating constraints, the fibre inflation model is safe.

      A word on our notation: we use natural units in which the speed of light and Planck's constant are set to 1. Unless otherwise indicated we work in units in which the Planck mass is also set to 1. We work in the context of a spatially flat homogeneous and isotropic background cosmology given by the metric

      ds2=dt2+a(t)2dx2,

      (1)

      where t is physical time, x are the comoving spatial coordinates, and a(t) is the cosmological scale factor. The Hubble expansion rate is

      H(t)˙aa.

      (2)

      There are two coupling constants which appear in string theory. The first is the string coupling constant gs which measures the strength of the string interactions. The value of gs is set by the expectation value of the string theory dilaton field. The second coupling constant α measures the strength of quantum effects. It is given by the square of the string length.

    2.   Review of fibre inflation
    • Inflation driven by a scalar field has become the paradigm of early universe cosmology (see e.g. [29] for a review). Unless the potential of the inflaton field is finely tuned, the energy scale at which inflation takes place corresponds that of particle physics Grand Unification, which is close to the expected string scale in many setups [30]. Hence, it is reasonable to search for realizations of inflation in the context of string theory.

      String theory is anomaly-free only in ten space-time dimensions. Hence, to make contact with our world, six of the spatial dimensions need to be compactified to a small size. In this context, the effective field theory in our four space-time dimensions will contain many scalar fields, some of which could be candidates for the inflaton (see e.g. [31] for a review of the connection between inflation and string theory).

      The scalar field candidates to be the inflaton include the separation between branes [32-35], Kähler moduli (roughly speaking the radii of the compactified dimensions) [36, 37], and axions [38, 39].

      Scalar field-driven inflation usually requires the field to be slowly rolling (technically this means that the acceleration term in the Klein-Gordon equation for the field is negligible, and the kinetic energy is also negligible). These conditions are in general hard to realize. In particular, if we want the slow-roll trajectory of the scalar field to be a local attractor in initial condition space, then field values greater than the Planck scale are required (see e.g. [40] for a recent review and [41-43] for some original works). As reviewed in [31], such large field values are hard to obtain in a controlled way from string theory. There are, however, a couple of interesting possibilities, the first being axion monodromy inflation [44], the second the Large Volume Scenario [45]. Both of these scenarios may, in fact, lie in the string theory swampland [26], but this is not the topic of our work.

      Fibre inflation [24] is a particular realization of inflation in the context of the Large Volume Scenario. It is based on a flux compactification of Type IIB superstring theory, with the compact manifold being a Calabi-Yau three fold X [46, 47] with fibres. The resulting low energy theory is four space-time dimensional N=1 supergravity which in turn is given by a Kähler potential K (which determines the kinetic part of the action) and a superpotential W which determines the potential energy term.

      The Kähler potential depends on the total volume V of the threefold, on the axio-dilaton S, on a holomorphic (3,0) form Ω, on complex structure moduli Uα;α= 1,,h2,1(X) and on Kähler moduli Ti;i=1,,h1,1(X). Here the numbers h2,1 and h1,1 characterize the topology of the compact space X.

      The axio-dilaton field S and the complex structure moduli are fixed by turning on fluxes (both NS-NS and RR fluxes). At tree level the Kähler moduli Ti are not fixed and neither is the volume V, but the volume is determined from the Ti (Note that in the string gas cosmology [48] approach to merging string theory with cosmology (see [49-51] for a review) the Kähler moduli are automatically fixed [52-56] by the interplay of winding and momentum modes. The winding modes are frozen out in the effective field theory approach). There are no corrections to the superpotential at any finite order in gs and α. However, there are non-perturbative corrections (e.g. via gaugino condensation), and they give exponential terms of the form

      ΔW=iAieaiTi,

      (3)

      where Ai are coefficients. The Kähler potential obtains perturbative corrections of the form:

      ΔK=2ln(V+ξ2g3/2s),

      (4)

      where ξ is a constant proportional to the Euler number of X.

      With these corrections it is possible to generate a minimum of the potential corresponding to large volume. But this depends on the geometry of X. The Calabi-Yau manifold X must have a “blow-up” mode (a direction in which the radius is large in string units) to get the minimum of the potential to be at a large value of V. In this case, the induced potential stabilizes the Kähler moduli in the directions which are not large, but not the value of the blowup mode. Note that extra structures are needed to uplift the AdS (anti-de-Sitter) minimum to a dS (de Sitter) (and this might not be consistent from the point of view of string theory).

      The specific example considered in [24] has two blowup directions which yield two light moduli fields τ1 and τ3. The overall volume V is determined by the values of the blowup fields. Without string loop corrections, the potential is flat in a direction which in field space is roughly τ1. If we fix the total volume and τ3, then in fact the flat direction is exactly the τ1 direction. In τ3 direction there is a minimum of the potential, yielding a valley in the potential viewed as a function of τ1 and τ3 which is the inflaton direction. For values of τ3 larger than the value at the minimum of the potential, the potential increases slowly. How slowly will be the key question we ask in Section 4. String loop effects will slightly lift the potential in τ1 direction, thus allowing for inflation.

      As done in [24], we here consider two versions of the fibre inflation model. In the first we fix τ3 and V at the minima of their potentials. We then obtain a single field inflation model given by the Lagrangian

      L=38(μτ1μτ1τ21)V(τ1)

      (5)

      with potential

      V(τ1)=V0+(Aτ21BVτ1+Cτ1V2)W20V2,

      (6)

      where V0,A,B,C and W0 are constants given by the three-fold X, and the volume V is taken to be constant.

      In the second version we keep τ3 fixed at the minimum but allow V to fluctuate away from its minimum. The kinetic piece of the Lagrangian is then both non-canonical and non-diagonal

      Lkin=38(μτ1μτ1τ21)+12(μτ1μVτ1V)12(μVμVV2),

      (7)

      and then potential is

      V=[μ4(ln(cV))3/2+μ3]W20V3+δupV4/3+(Aτ21BVτ1+Cτ1V2)W20V2

      (8)

      where c,μ4,μ3,W0,δup,A,B and C are constants which are determined by the specific manifold X.

    3.   Self-resonance in the single-field model
    • Here we consider the model given by (5) and (6). In order to obtain the Lagrangian of a canonically normalized field we make the field redefinition (recall that we are using units in which the Planck mass is set to 1)

      ϕ=32lnτ1.

      (9)

      Then, the Lagrangian has the canonical form

      L=12μϕμϕV(ϕ),

      (10)

      where the potential of the inflaton is

      V(ϕ)=C2V10/3[(3R)4(1+R6)eκϕ/2+(1+2R3)e2κϕ+Reκϕ],

      (11)

      where R is a small parameter, R1, and κ23Mp.

      The homogeneous mode of the inflaton satisfies the Klein-Gordon equation

      ¨ϕ+3H˙ϕ+Vϕ=0,

      (12)

      where the subscript indicates the argument with respect to which the derivative is taken. For large field values κϕ>1 (but not Rexp(κϕ)>1) the inflationary slow-roll conditions on the potential are satisfied. Hence, for such values of ϕ inflation can take place. Once κϕ<1 inflation ends and ϕ begins to oscillate about the minimum of the potential. In the absence of cosmic expansion (i.e. setting H=0 in the above equations), the frequency of the oscillations is mϕ. Through the nonlinear term in the evolution equation, self-resonance (i.e. resonance of fluctuations of the field ϕ itself) is possible. However, to check whether this indeed takes place we need to take into account the expansion of space.

    • 3.1.   Evolution without cosmic expansion

    • First we consider in more detail the evolution for H=0. Once κϕ<1, the slow-roll conditions break down, and the inflaton rolls to the minimum of its potential and begins to oscillate about it.

      For small field amplitudes ϕ0<1 (in Planck units, Mp=1), the potential can be approximated as

      Vapp(ϕ)12m2ϕϕ2+16gϕ3,

      (13)

      with

      m2ϕ=C2(6+7R)κ22V10/3

      (14)

      g=C2(17R+30)κ34V10/3.

      (15)

      Once the amplitude becomes small enough such that the effects of the nonlinear terms in the equation of motion are negligible (this always happens since the inflaton loses energy), the oscillations become harmonic:

      ϕϕ0cos(ωt),

      (16)

      with ω=mϕ. We have denoted the amplitude of the homogeneous mode by ϕ0.

      The potential and a typical example of the evolution of the inflaton are shown in Fig. 1. The right panel shows the potentian and its approximated version, the left panel shows the time evolution of ϕ.

      Figure 1.  (color online) The evolution (with H=0) of the homogeneous mode of the inflaton (vertical axis) as a function of time (horizontal axis) during reheating is shown in the left panel. The right panel shows both the full and the approximated potential, the dashed red line being Vapp and the blue line the full potential V.

      To study the self-resonance, let us consider the perturbation of the inflaton,

      δϕ(x,t)ϕ(x,t)ϕ(t),

      (17)

      where the brackets indicate spatial averaging. To linear order, the perturbation mode (in Fourier space) satisfies the equation of motion

      ¨δϕk+(k2+Vϕϕ)δϕk=0,

      (18)

      where we omitted the Hubble damping and set a=1. For a large amplitude of oscillation, one should calculate Vϕϕ with the complete potential (11) since the higher order interactions are important in this case. As the amplitude decays to be small, the approximation (13) works well. We will consider this case in what follows. With the solution (16) we have

      d2Vappdϕ2=m2ϕ+gϕ0cos(ωt).

      (19)

      Thus we have a Mathieu equation for the perturbation mode δϕk,

      ¨δϕk+(k2+m2ϕ+gϕ0cos(ωt))δϕk=0.

      (20)

      Using the rescaled time coordinate τ=12ωt, we get a standard Mathieu equation,

      δϕk+[Ak2qcos(2τ)]δϕk=0,

      (21)

      where

      Ak=4(k2+m2ϕ)ω2q=2gϕ0ω2.

      (22)

      This equation has the following solution

      δϕk=e˜μkτP1(k,τ)+e˜μkτP2(k,τ),

      (23)

      where P1,2 are periodic functions, and ˜μk is the famous Floquet exponent. If the Floquet exponent has nonvanishing real part, then there is parametric resonance of the fluctuations [7]. In terms of physical time we can write

      μkt˜μτ.

      (24)

      Note that μk has dimensions of energy, whereas ˜μ is dimensionless.

      To know how large the Floquet exponent is, let us analyze the parameters involved. After the end of inflation, the field oscillates about the minimum of the potential, and the amplitude is small compared to the Planck scale. Hence, we can regard the background system as a harmonic oscillator. Combining the equations (16) and (13) we have ωmϕ. According the reference [24], we have

      κ=23Mpg=5m2ϕ3Mp.

      (25)

      We immediately get

      Ak4k2m2ϕ+4,q10ϕ03Mp.

      (26)

      Since ϕ0<Mp (as discussed earlier), we can conclude that

      Ak>4,q<103.

      (27)

      From the Fig. 2, we see that the rescaled Floquet exponent is smaller than 0.5. The first band is the most important one. This implies that the long wavelength modes dominate particle production. Since the phase space of modes grows as k3, the dominant modes are in fact the modes with a value of k2 close to the upper end of the resonance band. Note that these modes are sub-Hubble (as can be seen by making use of the Friedmann equation), hence justifying the particle interpretation. The Floquet exponent for some typical values of q is shown in Fig. 2.

      Figure 2.  (color online) The band structure for the Floquet exponent of equations (26) and (36). On the left, the horizontal axis represents the parameter q, and the vertical axis is the parameter k2/m2ϕ, which represents the rescaled wavelength. The colored region is the parameter space where Re[μk]0, with the color coding indicating the value of the Floquet exponent. The right panel shows the Floquet exponent for values of k in the first parametric resonance band. The three curves are for q=3 (blue line), q=2.4 (red line), and q=1.8 (black line), which correspond to the inflaton oscillation with amplitude Φ=0.5, Φ=0.4, and Φ=0.3 respectively.

      To have an exact feeling about the efficiency of the self-resonance, we need to compare the Floquet exponent with the Hubble parameter H. If the former is larger than the latter, the exponential increase in the particle number is significant. We have

      Re[μk]H=mϕ2H2Re[μk]mϕ,

      (28)

      which is larger than one. For the value of mϕ given by the parameters of the model (see 14) we have H101mϕ at the beginning of reheating, and with the maximal value of the rescaled Floquet exponent being about 0.5, the ratio (28) is about 5. On the other hand, since the ratio is not significantly larger than one, it is important to extend our analysis and include the effects of the expanding background.

    • 3.2.   Including the effects of the expanding background

    • Now let us include the cosmic expansion. In this case the evolution should be described by the equation (1). For some given initial conditions, it has a typical evolution as in the Fig. 3 (which is obtained by numericallly solving the evolution equation). As expected, the inflaton starts to oscillate about the minimum of the potential with an amplitude smaller than the Planck scale. Hence, it can be viewed as small perturbation about the minimum ϕ, and we can write

      Figure 3.  (color online) The evolution of the inflaton (left) and the comoving number density (right), calculated by nk=ωk2(|˙Xk|2ω2k+|Xk|2)12, with k=103mϕ. The inflaton (vertical axis) oscillates after inflation as a function of time (horizontal axis), with the amplitude decaying and the periodicity being almost a constant. The field is in Planck unit and time in m1ϕ unit.

      δϕ=ϕϕ.

      (29)

      The perturbation mode (in Fourier space) satisfies the equation of motion

      ¨δϕk+3H˙δϕk+(k2a2+Vϕϕ)δϕk=0.

      (30)

      It is convenient to eliminate the fraction term by making use of a field transformation

      Xk=a3/2δϕk.

      (31)

      In terms of Xk, the equation (30) becomes

      ¨Xk+(k2a2+Vϕϕ94H232˙H)Xk=0,

      (32)

      where

      ω2k(t)=k2a2+Vϕϕ94H232˙H.

      (33)

      This equation was studied in detail in [11], where it was shown that particle production occurs when the adiabaticity condition on the time dependence of the frequency is violated. This happens when the background field ϕ takes on a local maximum (since then the Floquet exponent takes on a local minimum), as long as the amplitude of the background field remains sufficiently large.

      We will give an approximate treatment (see e.g. [10] where it was shown that this analysis gives qualitatively the same results as the rigorous treatment of [11]) and set a(t)=1 in the fluctuation equation. The equation for the perturbation mode then becomes

      ¨Xk+(k2+Vϕϕ)Xk=0.

      (34)

      Now we compute Vϕϕ, using the approximation (13). During the oscillations, the background field goes like ϕ(t)=Φ(t)cos(ωt) with Φ(t)λt, which is shown in Fig. 3. Note that the inflaton starts to oscillate when the amplitude Φ(t) is small, Φ(t)<1 (in Planck units), and the amplitude continues to decrease. The equation (34) then becomes

      ¨Xk+(k2+m2ϕ+gΦ(t)cos(ωt))Xk=0.

      (35)

      It is convenient to write it in terms of the new coordinate τ=12ωt, since this yields the standard form of the Mathieu equation,

      Xk+(Ak2qcos(2τ))Xk=0,Ak=4k2ω2+4m2ϕω2,q=2gΦω2

      (36)

      where a prime represents the derivative with respect to τ. As a good approximation, we can use ωmϕ. Then we find

      Ak4k2m2ϕ+4,q10Φ(t)3Mp.

      (37)

      Typically, the inflaton oscillates for Φ(t)<1, hence q<6. We see that in this case we have narrow resonance. Since the amplitude Φ(t) decays, the parameter q becomes smaller and smaller, and this makes the resonance band narrower. From the Fig. 2, we conclude that the cosmic expansion reduces but does not cut off the efficiency of particle production.

    • 3.3.   Termination of preheating

    • Due to the nonlinearities in the equation of motion for the background field, particle production will lead to correction terms. Once these correction terms become comparable to the mass term driving the oscillations, preheating will shut off. The leading nonlinearity comes from the cubic term in the approximate potential and yield the corrected equation of motion

      ¨ϕ+3H˙ϕ=m2ϕϕ12gϕ2.

      (38)

      Each Fourier fluctuation mode δϕk contributes to the correction term. The criteria for termination of preheating thus becomes

      m2ϕΦ0<g(δϕ)2,

      (39)

      where, as before, Φ0 is the amplitude of the background field. The quantity (δϕ)2 is the contribution of the fluctuation modes to the background. It is given by

      (δϕ)2=d3k|δϕk|2,

      (40)

      where δϕk are the fluctuations computed earlier.

      To compute this expectation value, consider the Fourier expansion of the fluctuation field

      δϕ=d3kfkϵkeikx,

      (41)

      where fk is the amplitude of the fluctuation mode and ϵk is a random variable with normalization

      ϵkϵk=δ3(kk).

      (42)

      Since particle production is dominated by the first resonance band, we can estimate (40) by restricting the integral to this first band and using

      fk=1keμkt,

      (43)

      where we have assumed vacuum initial conditions for the amplitude of the Fourier modes. Inserting (42) and (43) into (40), we can then estimate the resulting integral by k2 (the upper cutoff value of the first resonance band) multiplied by the factor of increase, using the maximal value of the Floquet exponent in the first resonance band. This yields

      (δϕ)2k2e2μkt,

      (44)

      where μk can be taken to be the maximal value of the Floquet exponent.

      The termination criterium (39) then becomes

      2μkt<ln(m2ϕΦ0gk2).

      (45)

      Since the oscillation period is given by m1ϕ and since the maximal value of the Floquet exponent is about 0.5 for q=0.7, we find that the number of oscillation periods before the termination criterion is satisfied is given by the right hand side of the above equation. Making use of k101mϕ and that g3m2ϕ/Mp, we find that the number of oscillations is of the order 10.

      A better termination criterium is to check until which time adiabaticity violation takes place. There is adiabaticity violation for a few oscillation times. This second analysis confirms the conclusion that the expansion of the universe and back-reaction do not prevent the preheating instability.

      On the other hand, the preheating phase does not convert all of the energy of the initial inflaton condensate to particles. To estimate the fraction of the energy which is transferred during the initial stage of parametric resonance we can compute the energy density in particles, using

      ρp4πk0dkk3e2μkt.

      (46)

      Inserting the result from (45) and using k101mϕ we find

      ρpρ0102m2ϕgΦ0

      (47)

      which is of the order of a percent.

      In conclusion, in this section we have shown that parametric self-resonance is effective in the fibre inflation model with the model parameters suggested in [24]. Note that in [27], a much larger value of the Hubble constant was chosen. This choice plus small differences in the shape of the potential lead to significant differences in the value of the Floquet exponent. The authors of [27] obtain very similar results as we do if they apply their analysis to the same form of the potential and parameter values as we use (private communication from F. Cefala). It is also important to stress, as done in [27], that the initial presence of the preheating instability does not imply that this instability will drain most of the inflaton energy.

    4.   Two-field model
    • Now that we have shown that the initial stages of preheating occur via parametric self-resonance of inflaton fluctuations, we must study whether the parametric resonance might lead to entropy fluctuations of other light fields. In the fibre inflation model of [24], the Lagrangian including the second light degree of freedom is given by (see Section 2)

      L=38τ21μτ1μτ1+12τ1Vμτ1μV12V2μVμVV(τ1,V),

      (48)

      where the potential is

      V(τ1,V)=[μ4(ln(cV))3/2+μ3]W20V3+δupV4/3+(Aτ21BVτ1+Cτ1V2)W20V2.

      (49)

      We transform to canonically normalized fields ϕ1 and ϕ2 which are defined by

      τ1=eaϕ1+bϕ2,V=ebϕ1+cϕ2.

      (50)

      The dimensionless parameters above are b=2/(7+42), a=2b2, and c=32b2. In terms of ϕ1 and ϕ2 one has the standard two-field theory

      L=12(μϕ1μϕ1+μϕ2μϕ2)V(ϕ1,ϕ2).

      (51)

      Here the potential V(ϕ1,ϕ2) is

      V(ϕ1,ϕ2)=μ4W20(ln(c)+bϕ1+cϕ2)3/2+μ3W20e3bϕ1+3cϕ2+δupe4(bϕ1+cϕ2)/3+AW20e2(a+b)ϕ12(b+c)ϕ2BW20e(12a+3b)ϕ112(b+3c)ϕ2+CW20e(a4b)ϕ1+(b4c)ϕ2.

      (52)

      The evolution of the background fields is described by the equations

      ¨ϕi+3H˙ϕi+Vϕi=0.

      (53)

      The evolution of the two fields is shown in the Fig. 4.

      Figure 4.  (color online) The inflation potential V(ϕ1,ϕ2) and the evolution of the scalar fields, with the blue line being ϕ1(t) and the red line ϕ2(t). The fields are in Planck unit and t is in 103M1p unit. We see that inflation mostly happens along the ϕ1. Note that we used the parameters of SV2 in [24].

      In the following, we will show that in spite of the presence of two relevant field directions, entropy perturbations produced on super-Hubble scales during the preheating phase are too small to have an important effect of the spectrum of curvature perturbations. It turns ouit that the basic reason for this conclusion is that the second field is not sufficiently light compared to the Hubble expansion rate.

      We will work in longitudinal gauge (see e.g. [57] for a review of the theory of cosmological perturbations and [58] for an introductory overview) in which the metric, including scalar metric fluctuations, takes the form

      ds2=(1+2Ψ)dt2+a(t)2(12Φ)δijdxidxj,

      (54)

      where Ψ(x,t) and Φ(x,t) are the fluctuating degrees of freedom of the metric. If matter has no anisotropic stress (which is the case for scalar field matter), then one of the Einstein equations leads to the conclusion that the two potentials are equal. In this case, the curvature perturbation is

      R=ΦH(HΦ+˙Φ)4πG(ρ+P).

      (55)

      In the absence of entropy fluctuations, the curvature perturbation R is conserved on super-Hubble scales. More generallly, the equation of motion for R is

      ˙R=Hρ+PδPnadHρ+P˙ρ˙P2Φ4πGa2,

      (56)

      where the non-adiabatic pressure perturbation δPnad is defined by

      δPnad=δP˙P˙ρδρ.

      (57)

      Recalling the definition of the total entropy perturbation

      SH(δP˙Pδρ˙ρ),

      (58)

      one immediately gets

      ˙R=˙Pρ+PSHρ+P˙ρ˙P2Φ4πGa2.

      (59)

      On super-horizon scales, kaH, the second term on the right-hand side can be ignored. Hence, on these large scales the evolution of the curvature perturbation is dominated by the entropy perturbation,

      ˙R˙Pρ+PS2˙θ˙P3˙σ(ρ+P)δs.

      (60)

      where δs=δϕ2cosθδϕ1sinθ, and the variables θ and σ are defined in the Appendix.

      The application of the above general formulas to two scalar field models has been studied in many references. We follow the treatment of [59] which is summarized in the Appendix. The two fields are expanded into adiabatic field σ and entropy field s. It is found that the entropy field fluctuation δs obeys the equation

      ¨δs+3H˙δs+(k2a2+Vss+3˙θ2)δs=˙θ˙σk22πGa2Φ,

      (61)

      where the subscripts on V indicate that the second derivative with respect to the entropy field is taken. The variable θ describes the angle between the original field basis and the adiabatic/entropy field basis. We see that the effective mass of the entropy mode is given by

      m2eff=Vss+3˙θ2.

      (62)

      For the two-field theory (51), we have used the numerically determined background trajectory to determine the angle θ and the second derivative of the potential with respect to the entropy field. In this way, we have numerically evaluated the above effective mass. The results for the evolution of the effective mass square m2eff are shown in Fig. 5.

      Figure 5.  (color online) The evolution of the effective mass of the entropy mode and its comparison to H2. Here m2eff is in M2p unit and t is in 103M1p unit.

      We see that the effective mass is much larger than 3H22, so the entropy mode is heavy. Hence, we expect the entropy mode to undergo damped oscillations. The top two panels show the evolution towards the end of the period of slow-roll inflation, the bottom two during the period of oscillations of the inflaton field about the minimum of its potential. The top left panel shows the evolution of the square of the mass, the top right the ratio of this square mass to the square of the Hubble expansion rate. The two bottom panels show a different range on the horizontal axis. With the scaling on the left we can track the overall amplitude, As is apparent, the oscillations of the inflaton lead to large amplitude peaks in the effective mass. However, the amplitudes vary dramatically from peak to peak, although the distance in time between the peaks appears to be nearly periodic. Basically, this is because of the structure of the potential (53). When expanded around the minimum for small amplitudes of oscillations, the potential has a crossed mass term for ϕ1 and ϕ2 and this breaks the harmonic nature of the oscillations. Since the time dependence of the effective square mass is far from periodic, there is no Floquet theory (parametric resonance) instability.

      Indeed, given the numerically computed mass term we can numerically solve the entropy mode equation. The result is shown in Fig. 6. To extract the effects of the expansion of space, we first make a transformation of variables, δ˜s=a3/2δs. Then the equation (79) for the entropy mode becomes

      Figure 6.  (color online) The evolution of a3/2δs and the scale factor a(t) during inflation and reheating, with t in 103M1p unit. Note that t=0 is just the starting point of our numerical evolution but not of inflation, so we do not show the full 60 e-foldings.

      ¨δ˜s+(Vss+3˙θ294H232˙H)δ˜s=0

      (63)

      on large scales where gradient terms are negligible. Since the effective square mass is large compared to H2, the H2 and ˙H terms in the equation are negligible. We have neglected them in the numerical solution for δ˜s shown in Fig. 6. We see that δ˜s oscillates with almost constant amplitude during inflation. In the reheating process, the amplitude has both amplification and decay. This can be understood by looking at the evolution of the effective mass in Fig. 5. The effective mass of the entropy mode does not change much during inflation and this results in a nearly constant amplitude of δ˜s. As mentioned above, during reheating, the effective mass has periodic narrow peaks. Except during the narrow peak intervals, the effective mass is almost constant, and no growth of δ˜s. Since the peak structures are not periodic, they do not allow for the self-resonance for the entropy mode. Since δs=a3/2δ˜s, considering the expansion of the universe, we conclude that the entropy mode is suppressed.

      Note that in this analysis we have assumed that the background scalar field dominates the cosmology during reheating. This is a good approximation at the beginning of the reheating phase, but will no longer be good once backreaction of the produced particles becomes important.

      Since the entropy fluctuations do not grow, the induced effect they have on curvature fluctuations is expected to be negligible. We can check this by solving the equation (80) for the adiabatic mode (which gives the curvature perturbation via R=H˙σQ) given the entropy source function which we have determined. Like in the case of the entropy mode, we first extract the effects of the expansion of space via the transformation Q=a3/2˜Q. This rescaled field obeys the equation

      ¨˜Q+(k2a2+Vσσ˙θ28πG(3˙σ2+2˙σ¨σH˙σ2˙HH2)94H232˙H)˜Q=2(¨θ32˙θHVσ˙σ˙HH)δ˜s+2˙θδ˙˜s.

      (64)

      The Evolution of the adiabatic mode is shown in Fig. 7. The left plot shows the evolution of ˜Q. Combining the evolution of the scale factor in Fig. 6, we see that the adiabatic mode Q=a3/2˜Q does not grow. Actually, since

      Figure 7.  (color online) The evolution of the mode ˜Q and ln(H˙σ˜Q) during inflation and reheating, with t in 103M1p unit.

      lnR=ln(H˙σ˜Q)32lna+const,

      (65)

      from the plot of ln(H˙σ˜Q) in Fig. 7 (right) we conclude that the curvature perturbation R is a constant during inflation and reheating.

    5.   Conclusions
    • We have studied the reheating phase in fibre inflation, taking the parameters used in [24]. In this model, there are two relevant “light” degrees of freedom. They come from the moduli fields of the model. In a first approximation, we (following [24]) have focused on the modulus field which taken by itself can lead to inflation. We have shown that the reheating phase in this model begins with a phase of parametric self-resonance during which a not negligible fraction of the inflaton energy is transferred to low momenta quanta of the inflaton.

      We then considered the effects of the second lightest field. In order for the model to be safe, it is important to check that there is no strong parametric resonance of the fluctuations in this second field. Otherwise, the induced entropy fluctuations might lead to a contribution to curvature fluctuations which would destroy the successful predictions of the model. We have shown that, although the second field is light compared to the string scale, it is sufficiently heavy compared to the Hubble scale at the end of inflation such that super-Hubble entropy fluctuations cannot grow. A similar approach on this issue can be found in the nice work [60]. The authors numerically studied the preheating phase in chaotic inflation model, with metric perturbations included. They showed that the curvature perturbation could be amplified due to the conversion from the entropy mode, thus the theory can be a realization of curvaton mechanism. There is no conflict with our results since the entropy mode considered in [60] is light enough.

      Whether the entropy perturbations violate the reheating process or not depends on the effective mass of the entropy mode. The effective mass is related to the parameters in equation (52) and the evolution of the inflaton, as can be seen from Eq. (82) in appendix. They are determined by string compactifications and moduli stabilization hence are not free parameters. It is not clear how the entropy mode mass is affected by these parameters. We have chosen the parameters of fibre inflation suggested in [24]. It would be interesting to scan a wider parameter space of fibre inflation models and search for regions where there is a parametric instability of the entropy mode.

      We thank Ryo Namba, Ziwei Wang, and Daisuke Yoshida for helpful discussions. In particular, we thank F. Cefala for useful communications concerning the relation of our work to that of [27].

    Appendix A
    • Here we review the analysis of entropy perturbation S and induced curvature fluctuations in the case of two scalar field matter models, following the analysis of [59].

      For scalar field matter (fields with canonical kinetic terms), the total energy density ρ and pressure P are given by

      ρ=I12˙ϕ2I+V

      and

      P=I12˙ϕ2IV,

      where V is the potential energy function and the sum runs of the different fields. The time derivatives and fluctuations of these functions are given by

      ˙ρ=I(˙ϕI¨ϕI+VϕI˙ϕI),

      ˙P=I(˙ϕI¨ϕIVϕI˙ϕI),

      δρ=I(˙ϕIδ˙ϕI˙ϕ2IΦ)+δV,

      δP=I(˙ϕIδ˙ϕI˙ϕ2IΦ)δV.

      Using these quantities, the total entropy perturbation (58) turns to be

      S=2(˙V+3HI˙ϕ2I)δV+2˙VI(˙ϕIδ˙ϕI˙ϕ2IΦ)3I˙ϕ2IJ(2˙V+3HJ˙ϕ2J).

      If there is only one scalar field, then

      S=2Vϕ3˙ϕ(2˙V+3H˙ϕ2)(˙ϕδ˙ϕ˙ϕ2Φ¨ϕδϕ)=2Vϕ3˙ϕ(2˙V+3H˙ϕ2)(δρ+3H˙ϕδϕ).

      In this case, the total entropy perturbation can be reduced by using the perturbed Poisson equation (the 00 component of the Einstein equation) and the momentum constraint (the 0i component of the perturbed Einstein equation),

      δρ+3H˙ϕδϕ=2Φ4πGa2.

      The total entropy perturbation hence scales as Sk2a2Φ and is hence suppressed on super-Hubble scales. Hence, it is clear that the entropy perturbations vanish on large scales for single scalar field inflation models.

      For two field models there is a convenient field space decomposition. The original scalar fields ϕ1 and ϕ2 can be described by the adiabatic field along the inflationary trajectory, and the entropy field orthogonal to the trajectory. The adiabatic field and its perturbation are

      ˙σ=˙ϕ1cosθ+˙ϕ2sinθ,δσ=δϕ1cosθ+δϕ2sinθ,

      where the angle θ is defined by

      cosθ=˙ϕ1˙ϕ21+˙ϕ22,sinθ=˙ϕ2˙ϕ21+˙ϕ22.

      The entropy field simply vanishes along the background trajectory.

      The entropy perturbation is

      7 δs=δϕ2cosθδϕ1sinθ.

      Using the density perturbation (A5) and the momentum constraint (0i component of the Einstein equation) one can show that

      R=Φ+H˙σδσ,S=Vσ6πG˙σ2(3H˙σ+2Vσ)(k2a2Φ)+2˙θ3˙σδs.

      To know how large the entropy perturbation is, and its contribution to the curvature perturbation, it is necessary to know the evolution of the entropy perturbation mode δs and the adiabatic mode δσ. The evolution equations of the entropy mode and the adiabatic mode are

      ..δs+3H.δs+(k2a2+Vss+3˙θ2)δs=˙θ˙σk22πGa2Φ,

      ¨Q+3H˙Q+(k2a2+Vσσ˙θ28πGa3t(a3˙σ2H))Q=2t(˙θδs)2(Vσ˙σ+˙HH)˙θδs,

      where

      Qδσ+˙σHΦ.

      Now we see that the effective mass of the entropy mode is

      m2eff=Vss+3˙θ2.

      The effective mass can be expressed in terms of the slow-roll parameters,

      εIJ=M2pl2VϕIVϕJV2,ηIJ=M2plVϕIϕJV,

      where

      VϕI=VϕIandVϕIϕJ=2VϕIϕJ.

      These are the slow-roll parameters in the original field space. In the adiabatic and entropy field space, one has

      Vσ=Vϕ1cosθ+Vϕ2sinθ,

      Vs=Vϕ2cosθVϕ1sinθ,

      Vσσ=Vϕ1ϕ1cos2θ+2Vϕ1ϕ2sinθcosθ+Vϕ2ϕ2sin2θ,

      Vss=Vϕ2ϕ2cos2θ2Vϕ1ϕ2sinθcosθ+Vϕ1ϕ1sin2θ,

      Vσs=Vϕ1ϕ1sinθcosθ+Vϕ1ϕ2(cos2θsin2θ)+Vϕ2ϕ2sinθcosθ.

      One can show that

      (Vϕ1ϕ1Vϕ1ϕ2Vϕ2ϕ1Vϕ2ϕ2)=(cosθsinθsinθcosθ)(VσσVσsVsσVss)(cosθsinθsinθcosθ)

      Similarly, one has

      (η11η12η21η22)=(cosθsinθsinθcosθ)(ησσησsηsσηss)(cosθsinθsinθcosθ)

      Hence the result is Vss=3H2ηss. For ˙θ, we have

      ˙θ=Vs˙σ=˙ϕ1¨ϕ2˙ϕ2¨ϕ1˙σ2.

      Using the slow-roll equations it can be shown that

      ˙θ=Hησs.

      Hence, we finally have

      m2eff=3H2(ηss+η2σs).

      If the effective mass of the entropy mode is too heavy (larger than 32H), then the fluctuations coming from the entropy mode do not grow on large scales. Instead, they perform damped oscillations. Now let us compare this with the quantity 32H2,

      m2eff32H2=3(ηss+η2σs12)H2.

      This implies that if ηss+η2σs>12 then the entropy perturbation are suppressed on large scales, Otherwise, the entropy perturbation can grow and give a significant contribution to the curvature perturbations. However, one should note that this analysis is only applicable during inflation since the equation (A30) is obtained by using the slow-roll condition. For the evolution after slow-rolling, we should use the expression (A17) rather than (A29).

Reference (60)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return