Lepton mass matrix from double covering of A4 modular flavor symmetry

  • We study a double covering of modular $ A_4 $ flavor symmetry. To this end, we construct lepton models for canonical and radiative seesaw scenarios. Using irreducible doublet representations, heavier Majorana fermion masses are characterized by one free parameter that would differentiate from $ A_4 $ symmetry. symmetry. Through χ square numerical analysis, we demonstrate that both scenarios produce some predictions in case of normal hierarchy reproducing neutrino oscillation data. However, no solution satisfies the neutrino oscillation data in case of radiative seesaw of inverted hierarchy.
  • 加载中
  • [1] F. Feruglio, doi: 10.1142/9789813238053_0012,arXiv:1706.08749[hep-ph].
    [2] R. de Adelhart Toorop, F. Feruglio, and C. Hagedorn, Nucl. Phys. B 858, 437-467 (2012), arXiv:1112.1340[hep-ph doi: 10.1016/j.nuclphysb.2012.01.017
    [3] J. C. Criado and F. Feruglio, SciPost Phys. 5(5), 042 (2018), arXiv:1807.01125[hep-ph doi: 10.21468/SciPostPhys.5.5.042
    [4] T. Kobayashi, N. Omoto, Y. Shimizu et al., JHEP 11, 196 (2018), arXiv:1808.03012[hep-ph
    [5] H. Okada and M. Tanimoto, Phys. Lett. B 791, 54-61 (2019), arXiv:1812.09677[hep-ph doi: 10.1016/j.physletb.2019.02.028
    [6] T. Kobayashi, H. Okada, and Y. Orikasa, arXiv: 2111.05674 [hep-ph]
    [7] T. Nomura and H. Okada, Phys. Lett. B 797, 134799 (2019), arXiv:1904.03937[hep-ph doi: 10.1016/j.physletb.2019.134799
    [8] H. Okada and M. Tanimoto, Eur. Phys. J. C 81(1), 52 (2021), arXiv:1905.13421[hep-ph doi: 10.1140/epjc/s10052-021-08845-y
    [9] F. J. de Anda, S. F. King, and E. Perdomo, Phys. Rev. D 101(1), 015028 (2020), arXiv:1812.05620[hep-ph doi: 10.1103/PhysRevD.101.015028
    [10] P. P. Novichkov, S. T. Petcov, and M. Tanimoto, Phys. Lett. B 793, 247-258 (2019), arXiv:1812.11289[hep-ph doi: 10.1016/j.physletb.2019.04.043
    [11] T. Nomura and H. Okada, Nucl. Phys. B 966, 115372 (2021), arXiv:1906.03927[hep-ph doi: 10.1016/j.nuclphysb.2021.115372
    [12] H. Okada and Y. Orikasa, arXiv: 1907.13520 [hep-ph]
    [13] G. J. Ding, S. F. King, and X. G. Liu, JHEP 09, 074 (2019), arXiv:1907.11714[hep-ph
    [14] T. Nomura, H. Okada, and O. Popov, Phys. Lett. B 803, 135294 (2020), arXiv:1908.07457[hep-ph doi: 10.1016/j.physletb.2020.135294
    [15] T. Kobayashi, Y. Shimizu, K. Takagi et al., Phys. Rev. D 100(11), 115045 (2019) [Erratum: Phys. Rev. D 101(3), 039904 (2020)], arXiv: 1909.05139 [hep-ph]
    [16] T. Asaka, Y. Heo, T. H. Tatsuishi et al., JHEP 01, 144 (2020), arXiv:1909.06520[hep-ph
    [17] D. Zhang, Nucl. Phys. B 952, 114935 (2020), arXiv:1910.07869[hep-ph doi: 10.1016/j.nuclphysb.2020.114935
    [18] G. J. Ding, S. F. King, X. G. Liu et al., JHEP 12, 030 (2019), arXiv:1910.03460[hep-ph
    [19] T. Kobayashi, T. Nomura, and T. Shimomura, Phys. Rev. D 102(3), 035019 (2020), arXiv:1912.00637[hep-ph doi: 10.1103/PhysRevD.102.035019
    [20] T. Nomura, H. Okada, and S. Patra, Nucl. Phys. B 967, 115395 (2021), arXiv:1912.00379[hep-ph doi: 10.1016/j.nuclphysb.2021.115395
    [21] X. Wang, Nucl. Phys. B 957, 115105 (2020), arXiv:1912.13284[hep-ph doi: 10.1016/j.nuclphysb.2020.115105
    [22] H. Okada and Y. Shoji, Nucl. Phys. B 961, 115216 (2020), arXiv:2003.13219[hep-ph doi: 10.1016/j.nuclphysb.2020.115216
    [23] H. Okada and M. Tanimoto, arXiv: 2005.00775 [hep-ph]
    [24] M. K. Behera, S. Singirala, S. Mishra et al., arXiv: 2009.01806 [hep-ph]
    [25] M. K. Behera, S. Mishra, S. Singirala et al., arXiv: 2007.00545 [hep-ph]
    [26] T. Nomura and H. Okada, arXiv: 2007.04801 [hep-ph]
    [27] T. Nomura and H. Okada, arXiv: 2007.15459 [hep-ph]
    [28] T. Asaka, Y. Heo, and T. Yoshida, Phys. Lett. B 811, 135956 (2020), arXiv:2009.12120[hep-ph doi: 10.1016/j.physletb.2020.135956
    [29] H. Okada and M. Tanimoto, Phys. Rev. D 103(1), 015005 (2021), arXiv:2009.14242[hep-ph doi: 10.1103/PhysRevD.103.015005
    [30] K. I. Nagao and H. Okada, arXiv: 2010.03348 [hep-ph]
    [31] H. Okada and M. Tanimoto, JHEP 03, 010 (2021), arXiv:2012.01688[hep-ph
    [32] C. Y. Yao, J. N. Lu, and G. J. Ding, JHEP 05, 102 (2021), arXiv:2012.13390 [hep-ph
    [33] P. Chen, G. J. Ding, and S. F. King, JHEP 04, 239 (2021), arXiv:2101.12724 [hep-ph
    [34] M. Kashav and S. Verma, arXiv: 2103.07207 [hep-ph]
    [35] H. Okada, Y. Shimizu, M. Tanimoto et al., arXiv: 2105.14292 [hep-ph]
    [36] I. de Medeiros Varzielas, and J. Lourenço, arXiv: 2107.04042 [hep-ph]
    [37] T. Nomura, H. Okada, and Y. Orikasa, arXiv: 2106.12375 [hep-ph]
    [38] P. T. P. Hutauruk, D. W. Kang, J. Kim et al., arXiv: 2012.11156 [hep-ph]
    [39] G. J. Ding, S. F. King, and J. N. Lu, arXiv: 2108.09655 [hep-ph]
    [40] K. I. Nagao and H. Okada, arXiv: 2108.09984 [hep-ph]
    [41] Georgianna Charalampous, Stephen F. King, George K. Leontaris et al., arXiv: 2109.11379 [hep-ph]
    [42] H. Okada and Y. h. Qi, arXiv: 2109.13779 [hep-ph]
    [43] T. Nomura, H. Okada, and Y. h. Qi, arXiv: 2111.10944 [hep-ph]
    [44] T. Kobayashi, H. Otsuka, M. Tanimoto et al., arXiv: 2112.00493 [hep-ph]
    [45] A. Dasgupta, T. Nomura, H. Okada et al., arXiv: 2111.06898 [hep-ph]
    [46] X. G. Liu and G. J. Ding, arXiv: 2112.14761 [hep-ph]
    [47] T. Nomura and H. Okada, arXiv: 2201.10244 [hep-ph]
    [48] H. Otsuka and H. Okada, arXiv: 2202.10089 [hep-ph]
    [49] D. W. Kang, J. Kim, T. Nomura et al., arXiv: 2205.08269 [hep-ph]
    [50] K. Ishiguro, H. Okada, and H. Otsuka, arXiv: 2206.04313 [hep-ph]
    [51] T. Nomura, H. Okada, and Y. Shoji, arXiv: 2206.04466 [hep-ph]
    [52] T. Kobayashi, H. Otsuka, M. Tanimoto et al., arXiv: 2204.12325 [hep-ph]
    [53] T. Kobayashi, K. Tanaka, and T. H. Tatsuishi, Phys. Rev. D 98(1), 016004 (2018), arXiv:1803.10391[hep-ph doi: 10.1103/PhysRevD.98.016004
    [54] T. Kobayashi, Y. Shimizu, K. Takagi et al., Phys. Lett. B 794, 114-121 (2019), arXiv:1812.11072[hep-ph doi: 10.1016/j.physletb.2019.05.034
    [55] T. Kobayashi, Y. Shimizu, K. Takagi et al., PTEP 2020(5), 053B05 (2020), arXiv:1906.10341[hep-ph
    [56] H. Okada and Y. Orikasa, Phys. Rev. D 100(11), 115037 (2019), arXiv:1907.04716[hep-ph doi: 10.1103/PhysRevD.100.115037
    [57] S. Mishra, arXiv: 2008.02095 [hep-ph]
    [58] X. Du and F. Wang, JHEP 02, 221 (2021), arXiv:2012.01397[hep-ph
    [59] J. T. Penedo and S. T. Petcov, Nucl. Phys. B 939, 292-307 (2019), arXiv:1806.11040[hep-ph doi: 10.1016/j.nuclphysb.2018.12.016
    [60] P. P. Novichkov, J. T. Penedo, S. T. Petcov et al., JHEP 04, 005 (2019), arXiv:1811.04933[hep-ph
    [61] T. Kobayashi, Y. Shimizu, K. Takagi et al., JHEP 02, 097 (2020), arXiv:1907.09141[hep-ph
    [62] S. F. King and Y. L. Zhou, Phys. Rev. D 101(1), 015001 (2020), arXiv:1908.02770[hep-ph doi: 10.1103/PhysRevD.101.015001
    [63] H. Okada and Y. Orikasa, arXiv: 1908.08409 [hep-ph]
    [64] J. C. Criado, F. Feruglio, and S. J. D. King, JHEP 02, 001 (2020), arXiv:1908.11867[hep-ph
    [65] X. Wang and S. Zhou, JHEP 05, 017 (2020), arXiv:1910.09473[hep-ph
    [66] Y. Zhao and H. H. Zhang, JHEP 03, 002 (2021), arXiv:2101.02266 [hep-ph
    [67] S. F. King and Y. L. Zhou, JHEP 04, 291 (2021), arXiv:2103.02633 [hep-ph
    [68] G. J. Ding, S. F. King, and C. Y. Yao, arXiv: 2103.16311 [hep-ph]
    [69] X. Zhang and S. Zhou, arXiv: 2106.03433 [hep-ph]
    [70] Bu-Yao Qu, Xiang-Gan Liu, Ping-Tao Chen et al., arXiv: 2106.11659 [hep-ph]
    [71] T. Nomura and H. Okada, arXiv: 2109.04157 [hep-ph]
    [72] P. P. Novichkov, J. T. Penedo, S. T. Petcov et al., JHEP 04, 174 (2019), arXiv:1812.02158[hep-ph
    [73] G. J. Ding, S. F. King, and X. G. Liu, Phys. Rev. D 100(11), 115005 (2019), arXiv:1903.12588[hep-ph doi: 10.1103/PhysRevD.100.115005
    [74] X. G. Liu and G. J. Ding, JHEP 08, 134 (2019), arXiv:1907.01488[hep-ph
    [75] P. Chen, G. J. Ding, J. N. Lu et al., Phys. Rev. D 102(9), 095014 (2020), arXiv:2003.02734[hep-ph doi: 10.1103/PhysRevD.102.095014
    [76] C. C. Li, X. G. Liu, and G. J. Ding, arXiv: 2108.02181 [hep-ph]
    [77] P. P. Novichkov, J. T. Penedo, and S. T. Petcov, Nucl. Phys. B 963, 115301 (2021), arXiv:2006.03058[hep-ph doi: 10.1016/j.nuclphysb.2020.115301
    [78] X. G. Liu, C. Y. Yao, and G. J. Ding, Phys. Rev. D 103(5), 056013 (2021), arXiv:2006.10722[hep-ph doi: 10.1103/PhysRevD.103.056013
    [79] X. Wang, B. Yu, and S. Zhou, Phys. Rev. D 103(7), 076005 (2021), arXiv:2010.10159[hep-ph doi: 10.1103/PhysRevD.103.076005
    [80] C. Y. Yao, X. G. Liu, and G. J. Ding, Phys. Rev. D 103(9), 095013 (2021), arXiv:2011.03501[hep-ph doi: 10.1103/PhysRevD.103.095013
    [81] X. Wang and S. Zhou, arXiv: 2102.04358 [hep-ph]
    [82] M. K. Behera and R. Mohanta, arXiv: 2108.01059 [hep-ph]
    [83] S. Kikuchi, T. Kobayashi, M. Tanimoto et al., arXiv: 2206.08538 [hep-ph]
    [84] S. Kikuchi, T. Kobayashi, K. Nasu et al., arXiv: 2203.14667 [hep-ph]
    [85] I. de Medeiros Varzielas, S. F. King, and Y. L. Zhou, Phys. Rev. D 101(5), 055033 (2020), arXiv:1906.02208[hep-ph doi: 10.1103/PhysRevD.101.055033
    [86] T. Kobayashi and S. Tamba, Phys. Rev. D 99(4), 046001 (2019), arXiv:1811.11384[hep-th doi: 10.1103/PhysRevD.99.046001
    [87] S. Kikuchi, T. Kobayashi, H. Otsuka et al., JHEP 11, 101 (2020), arXiv:2007.06188[hep-th
    [88] Y. Almumin, M. C. Chen, V. Knapp-Pérez et al., JHEP 05, 078 (2021), arXiv:2102.11286 [hep-th
    [89] G. J. Ding, F. Feruglio, and X. G. Liu, SciPost Phys. 10, 133 (2021), arXiv:2102.06716 [hep-ph doi: 10.21468/SciPostPhys.10.6.133
    [90] F. Feruglio, V. Gherardi, A. Romanino et al., SciPost Phys. 05, 242 (2021), arXiv:2101.08718 [hep-ph
    [91] S. Kikuchi, T. Kobayashi, and H. Uchida, arXiv: 2101.00826 [hep-th]
    [92] P. P. Novichkov, J. T. Penedo, and S. T. Petcov, SciPost Phys. 04, 206 (2021), arXiv:2102.07488 [hep-ph
    [93] S. Kikuchi, T. Kobayashi, Y. Ogawa et al., arXiv: 2112.01680 [hep-ph]
    [94] P. P. Novichkov, J. T. Penedo and S. T. Petcov, arXiv: 2201.02020 [hep-ph]
    [95] G. Altarelli and F. Feruglio, Rev. Mod. Phys. 82, 2701-2729 (2010), arXiv:1002.0211[hep-ph doi: 10.1103/RevModPhys.82.2701
    [96] H. Ishimori, T. Kobayashi, H. Ohki et al., Prog. Theor. Phys. Suppl. 183, 1-163 (2010), arXiv:1003.3552[hep-th doi: 10.1143/PTPS.183.1
    [97] H. Ishimori, T. Kobayashi, H. Ohki et al., Lect. Notes Phys. 858, 1-227 (2012)
    [98] D. Hernandez and A. Y. Smirnov, Phys. Rev. D 86, 053014 (2012), arXiv:1204.0445[hep-ph doi: 10.1103/PhysRevD.86.053014
    [99] S. F. King and C. Luhn, Rept. Prog. Phys. 76, 056201 (2013), arXiv:1301.1340[hep-ph doi: 10.1088/0034-4885/76/5/056201
    [100] S. F. King, A. Merle, S. Morisi et al., New J. Phys. 16, 045018 (2014), arXiv:1402.4271[hep-ph doi: 10.1088/1367-2630/16/4/045018
    [101] S. F. King, Prog. Part. Nucl. Phys. 94, 217-256 (2017), arXiv:1701.04413[hep-ph doi: 10.1016/j.ppnp.2017.01.003
    [102] S. T. Petcov, Eur. Phys. J. C 78(9), 709 (2018), arXiv:1711.10806[hep-ph doi: 10.1140/epjc/s10052-018-6158-5
    [103] T. Kobayashi, H. Ohki, H. Okada et al. doi:10.1007/978-3-662-64679-3
    [104] A. Baur, H. P. Nilles, A. Trautner et al., Nucl. Phys. B 947, 114737 (2019), arXiv:1908.00805[hep-th doi: 10.1016/j.nuclphysb.2019.114737
    [105] T. Kobayashi, Y. Shimizu, K. Takagi et al., Phys. Rev. D 101(5), 055046 (2020), arXiv:1910.11553[hep-ph doi: 10.1103/PhysRevD.101.055046
    [106] P. P. Novichkov, J. T. Penedo, S. T. Petcov et al., JHEP 07, 165 (2019), arXiv:1905.11970[hep-ph
    [107] A. Baur, H. P. Nilles, A. Trautner et al., Phys. Lett. B 795, 7-14 (2019), arXiv:1901.03251[hep-th doi: 10.1016/j.physletb.2019.03.066
    [108] T. Kobayashi and H. Otsuka, Phys. Rev. D 101(10), 106017 (2020), arXiv:2001.07972[hep-th doi: 10.1103/PhysRevD.101.106017
    [109] T. Kobayashi and H. Otsuka, Phys. Rev. D 102(2), 026004 (2020), arXiv:2004.04518[hep-th doi: 10.1103/PhysRevD.102.026004
    [110] K. Ishiguro, T. Kobayashi, and H. Otsuka, Nucl. Phys. B 973, 115598 (2021), arXiv:2010.10782[hep-th doi: 10.1016/j.nuclphysb.2021.115598
    [111] M. Tanimoto and K. Yamamoto, arXiv: 2106.10919 [hep-ph]
    [112] K. Ishiguro, T. Kobayashi, and H. Otsuka, JHEP 01, 020 (2022), arXiv:2107.00487[hep-th
    [113] M. C. Chen, S. Ramos-Sánchez, and M. Ratz, Phys. Lett. B 801, 135153 (2020), arXiv:1909.06910[hep-ph doi: 10.1016/j.physletb.2019.135153
    [114] I. de Medeiros Varzielas, M. Levy, and Y. L. Zhou, JHEP 11, 085 (2020), arXiv:2008.05329[hep-ph
    [115] K. Ishiguro, T. Kobayashi, and H. Otsuka, JHEP 03, 161 (2021), arXiv:2011.09154[hep-ph
    [116] H. Abe, T. Kobayashi, S. Uemura, and J. Yamamoto, Phys. Rev. D 102(4), 045005 (2020), arXiv:2003.03512[hep-th doi: 10.1103/PhysRevD.102.045005
    [117] S. Kikuchi, T. Kobayashi, H. Otsuka et al., arXiv: 2201.04505 [hep-ph]
    [118] T. Kobayashi and H. Otsuka, Eur. Phys. J. C 82(1), 25 (2022), arXiv:2108.02700[hep-ph doi: 10.1140/epjc/s10052-022-09986-4
    [119] S. Kikuchi, T. Kobayashi, K. Nasu et al., arXiv: 2202.05425 [hep-th]
    [120] E. Ma, Phys. Rev. D 73, 077301 (2006), arXiv:hep-ph/0601225 doi: 10.1103/PhysRevD.73.077301
    [121] Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28, 870-880 (1962) doi: 10.1143/PTP.28.870
    [122] A. Gando et al. (KamLAND-Zen Collaboration), Phys. Rev. Lett. 117(8), 082503 (2016) [Erratum: Phys. Rev. Lett. 117(10), 109903 (2016)], arXiv: 1605.02889 [hep-ex]
    [123] I. Esteban, M. C. Gonzalez-Garcia, M. Maltoni et al., JHEP 09, 178 (2020), arXiv:2007.14792 [hep-ph doi: 10.1007/JHEP09(2020)178
  • 加载中

Figures(3) / Tables(4)

Get Citation
Hiroshi Okada and Yuta Orikasa. Lepton mass matrix from double covering of A4 modular flavor symmetry[J]. Chinese Physics C. doi: 10.1088/1674-1137/ac92d8
Hiroshi Okada and Yuta Orikasa. Lepton mass matrix from double covering of A4 modular flavor symmetry[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ac92d8 shu
Milestone
Received: 2022-06-11
Article Metric

Article Views(528)
PDF Downloads(22)
Cited by(0)
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:

Lepton mass matrix from double covering of A4 modular flavor symmetry

  • 1. Asia Pacific Center for Theoretical Physics (APCTP) - Headquarters San 31, Hyoja-dong, Nam-gu, Pohang 790-784, Korea
  • 2. Department of Physics, Pohang University of Science and Technology, Pohang 37673, Republic of Korea
  • 3. Institute of Experimental and Applied Physics, Czech Technical University in Prague, Husova 240/5, 110 00 Prague 1, Czech Republic

Abstract: We study a double covering of modular $ A_4 $ flavor symmetry. To this end, we construct lepton models for canonical and radiative seesaw scenarios. Using irreducible doublet representations, heavier Majorana fermion masses are characterized by one free parameter that would differentiate from $ A_4 $ symmetry. symmetry. Through χ square numerical analysis, we demonstrate that both scenarios produce some predictions in case of normal hierarchy reproducing neutrino oscillation data. However, no solution satisfies the neutrino oscillation data in case of radiative seesaw of inverted hierarchy.

    HTML

    I.   INTRODUCTION
    • Even after the discovery of the standard model (SM) Higgs, we wonder whether there is a way to accommodate neutral particles such as the active tiny neutrinos and dark matter (DM) candidates. To describe the observed neutrino sector, we might need a prescription about how to determine three mixings and two mass square differences in addition to CP phases, namely Majorana and Dirac phases, which are not accurately measured yet. Modular flavor symmetries are one of the most promising candidates to obtain predictive scenarios in the neutrino sector, given that these symmetries do not require many neutral bosons owing to a new degree of freedom: "modular weight." Moreover, DM can be considered stable by applying this degree of freedom. In fact, a great number of scientific papers on this research line have been published after the original paper [1]. For example, the modular $ A_4 $ flavor symmetry has been discussed in Refs. [1,352], in addition to $ S_3 $ in Refs. [5358], $ S_4 $ in Refs. [5971], $ A_5 $ in Refs. [64, 72, 73], double covering of $ A_4 $ in Refs. [7476], double covering of $ S_4 $ in Refs. [77, 78], and double covering of $ A_5 $ in Refs. [7982]. Reference [83] discusses the CP phase of quark mass matrices in modular flavor symmetric models at the fixed point of τ. Soft-breaking terms on modular symmetry are discussed in Ref. [84]. Other types of modular symmetries have also been proposed to understand masses, mixings, and phases of the SM in Refs. [8594]. Different applications to various physics fields such as dark matter and the origin of CP are found in Refs. [6, 7, 11, 14, 63, 104112]. Mathematical studies such as possible correction from Kähler potential, systematic analysis of the fixed points, and moduli stabilization are discussed in Refs. [94, 113116]. Recently, the authors of Ref. [117] proposed a scenario to derive four-dimensional modular flavor symmetric models from higher-dimensional theory on extra-dimensional spaces with modular symmetry. It constrains modular weights and representations of fields and modular couplings in the four-dimensional effective field theory. Higher-dimensional operators in the SM effective field theory are also constrained in the higher-dimensional theory, in particular, the string theory [118]. Non-perturbative effects relevant to neutrino masses have been studied in the context of modular symmetry anomaly [119].

      In this study, we apply double covering of modular $ A_4 $ symmetry $ T' $ to the lepton sector and show several predictions in cases of canonical seesaw and radiative seesaw scenarios [120]. Interestingly, $ T' $ has three has three irreducible doublet representations. Using these representations, heavier Majorana fermion masses are described by one free parameter (except τ) that would differentiate from $ A_4 $ symmetry. Note symmetry. Note that the mathematical part of ath $ T' $ has been thoroughly s has been thoroughly studied in Ref. [74], where the authors demonstrate a prediction in case of inverted hierarchy based on the canonical seesaw model.

      This paper is organized as follows. In Sec. II, we review our model setup in the lepton sector, for deriving the superpotential, charged-lepton mass matrix, Dirac Yukawa matrix, and Majorana mass matrix. In Sec. III, we formulate the neutrino mass matrix and its observables in case of canonical seesaw. Then, we address the radiative seesaw model, showing the soft breaking terms that play a crucial role in generating the neutrino mass matrix at one-loop level. In Sec. IV, we perform the numerical χ square analysis and show predictive figures for normal and inverted hierarchies of the canonical and radiative seesaws. Finally, we conclude and summarize our model in Sec. V. In Appendix A, we summarize formulas on the double covering of modular ormula $ A_4 $ symmetry.

    II.   MODEL
    • Next, we review our model in order to obtain the neutrino mass matrix. In addition to the minimal supersymmetric SM (MSSM), we introduce matter superfields including two right-handed neutral fermions $ N^c_{1,2} $ that belong to doublet under the modular $ T^\prime $ group with modular weight $ -1 $. We also add three chiral superfields $ \{\hat\chi, \hat{\eta}_1, \hat{\eta}_2 \} $ including two bosons $ \{\chi, \eta_1, \eta_2\} $ where there are superfields that are true singlets under the $ T^\prime $ group with $ \{-1,-1,-3\} $ modular weight. χ only plays a role in generating the neutrino mass matrix at one-loop level; therefore, $ \eta_{1,2} $ are inert bosons in addition to χ. Left-handed lepton doublets $ \{L_e,L_\mu,L_\tau\} $ are assigned to be triplet with $ -1 $ modular weight, while the right-handed ones $ \{e^c,\mu^c,\tau^c\} $ are set to be $ \{1,1'',1'\} $ with $ -1 $ modular weight. Two Higgs doublet $ H_{1,2} $ are invariant under the modular $ T^\prime $ symmetry. All the fields and their assignments are summarized in Table 1. Under these symmetries, the renormalizable superpotential is expressed as follows :

      Chiral superfields
      $\{{\hat{L}_{e}},{\hat{L}_{\mu}},{\hat{L}_{\tau}}\}$ $\{\hat{e}^c,\hat{\mu}^c,\hat{\tau}^c\}$ $\{\hat{N}^c_{1},\hat{N}^c_{2}\}$ $\hat{H}_1$ $\hat{H}_2$ $\hat{\eta}_1$ $\hat{\eta}_2$ $\hat{\chi}$
      $SU(2)_L$${\bf{2}}$${\bf{1}}$${\bf{1}}$${\bf{2}}$${\bf{2}}$${\bf{2}}$${\bf{2}}$${\bf{1}}$
      $U(1)_Y$$-\frac12$$1$$0$$\frac12$$-\frac12$$\frac12$$-\frac12$$0$
      $T'$$3$$\{1,1'',1'\}$$2$$1$$1$$1$$1$$1$
      $-k$${-1}$$-1$$-1$$0$$0$$-1$$-1$$-3$

      Table 1.  Field contents of matter chiral superfields and their charge assignments under $S U(2)_L\times U(1)_Y\times A_{4}$ in the lepton and boson sectors; $-k_I$ is the number of modular weight and the quark sector is the same as that of the SM.

      $ \begin{aligned}[b] {\cal W} =& \alpha_e [Y^{(2)}_3 \otimes\hat{e}^c\otimes \hat{L} \otimes\hat{H}_2] +\beta_e [Y^{(2)}_3 \otimes \hat{\mu}^c\otimes \hat{L} \otimes\hat{H}_2] +\gamma_e [Y^{(2)}_3 \otimes \hat{\tau}^c\otimes \hat{L}\otimes \hat{H}_2] \\ &+\alpha_\eta [Y^{(3)}_2 \otimes\hat{N}^c\otimes \hat{L} \otimes\hat{\eta}_1] +\beta_\eta [Y^{(3)}_{2''} \otimes\hat{N}^c\otimes \hat{L} \otimes\hat{\eta}_1] + M_{0} [Y^{(2)}_3 \otimes\hat{N}^c\otimes\hat{N}^c] \\ &+\mu_H \hat{H}_1 \hat{H}_2+\mu_\chi Y^{(6)}_1 \hat{\chi} \hat{\chi} + a Y^{(4)}_1 \hat{H}_1 \hat{\eta}_2 \hat{\chi} + b Y^{(4)}_1 \hat{H}_2 \hat{\eta}_1 \hat{\chi}, \end{aligned} $

      (1)

      where R-parity is implicitly imposed in the above superpotential, $ Y^{(2)}_3\equiv(f_1,f_2,f_3)^T $ is $ T' $ triplet with modular weight $ 2 $, and $ Y^{(3)}_{2^{(\prime\prime)}} \equiv(y^{(\prime\prime)}_1,y^{(\prime\prime)}_2)^T $ is $ T^\prime $ doublet with modular weight $ 3 $. The first line in Eq. (1) corresponds to the charged-lepton sector, while the second and third lines are related to the neutrino sector. The third line is particularly important if the neutrino mass matrix is induced at one-loop level as dominant contribution.

      After the electroweak spontaneous symmetry breaking, the charged-lepton mass matrix is given by

      $ \begin{align} m_\ell&= \frac {v_2}{\sqrt{2}} \left[\begin{array}{ccc} \alpha_e & 0 & 0 \\ 0 &\beta_e & 0 \\ 0 & 0 &\gamma_e \\ \end{array}\right] \left[\begin{array}{ccc} f_1 &f_3 & f_2 \\ f_2 &f_1 & f_3 \\ f_3 &f_2 & f_1 \\ \end{array}\right], \end{align} $

      (2)

      where $ \langle H_2\rangle\equiv [v_2/\sqrt2,0]^T $. Then the charged-lepton mass eigenstate is found as $ {\rm diag}( |m_e|^2, |m_\mu|^2, |m_\tau|^2)\equiv V_{e_L}^\dagger m^\dagger_\ell m_\ell V_{e_L} $. In our numerical analysis, we fix the free parameters $ \alpha_e,\beta_e,\gamma_e $ inserting the observed three charged-lepton masses by applying the following relations:

      $ \begin{align} &{\rm Tr}[m_\ell {m_\ell}^\dagger] = |m_e|^2 + |m_\mu|^2 + |m_\tau|^2, \end{align} $

      (3)

      $ \begin{align} &{\rm Det}[m_\ell {m_\ell}^\dagger] = |m_e|^2 |m_\mu|^2 |m_\tau|^2, \end{align} $

      (4)

      $ \begin{aligned}[b] &({\rm Tr}[m_\ell {m_\ell}^\dagger])^2 -{\rm Tr}[(m_\ell {m_\ell}^\dagger)^2]\\ =&2( |m_e|^2 |m_\mu|^2 + |m_\mu|^2 |m_\tau|^2+ |m_e|^2 |m_\tau|^2 ). \end{aligned} $

      (5)

      The Dirac matrix consists of $ \alpha_\eta $ and $ \beta_\eta $; $ N^c y_\eta L\eta_1 $ is given by

      $ \begin{align} y_\eta &= \left[\begin{array}{ccc} \dfrac{\beta_\eta}{\sqrt2} {\rm e}^{\frac{7\pi}{12}i}y^{\prime\prime}_2 & \alpha_\eta y_1 & \dfrac{\alpha_\eta}{\sqrt2} {\rm e}^{\frac{7\pi}{12}i}y_2 + \beta_\eta y^{\prime\prime}_1 \\ \dfrac{\beta_\eta}{\sqrt2} {\rm e}^{\frac{7\pi}{12}i}y^{\prime\prime}_1 + \alpha_\eta {\rm e}^{\frac{\pi}{6}i} y_2 & \beta_\eta {\rm e}^{\frac{\pi}{6}i} y^{\prime\prime}_2 & \dfrac{\alpha_\eta}{\sqrt2} {\rm e}^{\frac{7\pi}{12}i}y_1 \\ \end{array}\right]. \end{align} $

      (6)

      The heavier Majorana mass matrix is given by

      $ \begin{align} M_N &= M_0 \left[\begin{array}{ccc} f_2 & \dfrac1{\sqrt2} {\rm e}^{\frac{7\pi}{12}i} f_3 \\ \dfrac1{\sqrt2} {\rm e}^{\frac{7\pi}{12}i} f_3 & {\rm e}^{\frac{\pi}{6}i} f_1 \\ \end{array}\right] = M_0 {\tilde M}. \end{align} $

      (7)

      The heavy Majorana mass matrix is diagonalized by a unitary matrix $ V_N $ as follows: $ D_N\equiv V_N M_N V_N^T $, where $ N^c\equiv \psi^cV_N^T $, $ \psi^c $ is the mass eigenstate.

    III.   ACTIVE NEUTRINO MASS MATRIX

      A.   Canonical seesaw

      1.   Mass matrix of $ m_\nu $
    • If all the bosons have nonzero VEVs, the neutrino mass matrix is generated via tree-level as follows:

      $ \begin{align} m_\nu = \frac{v_{\eta_1}^2}{2 M_0} y_\eta^T \tilde M^{-1} y_\eta \equiv \kappa \tilde m_\nu, \end{align} $

      (8)

      where $ \kappa\equiv \dfrac{v_{\eta_1}^2}{2 M_0} $ and $ \langle\eta_1\rangle \equiv [0,v_{\eta_1}/\sqrt2]^T $. $ m_\nu $ is diagonalized by a unitary matrix $ V_{\nu} $; $ D_\nu=|\kappa| \tilde D_\nu= V_{\nu}^T m_\nu V_{\nu}= |\kappa| V_{\nu}^T \tilde m_\nu V_{\nu} $. Then, $ |\kappa| $ is determined by

      $ \begin{align} (\mathrm{NH}):\ |\kappa|^2= \frac{|\Delta m_{\rm atm}^2|}{\tilde D_{\nu_3}^2}, \quad (\mathrm{IH}):\ |\kappa|^2= \frac{|\Delta m_{\rm atm}^2|}{\tilde D_{\nu_2}^2}, \end{align} $

      (9)

      where $ \Delta m_{\rm atm}^2 $ denotes atmospheric neutrino mass difference squares, and NH and IH represent the normal hierarchy and inverted hierarchy, respectively. Subsequently, the solar mass different squares can be written in terms of $ |\kappa| $ as follows:

      $ \begin{align} (\mathrm{NH}):\ \Delta m_{\rm sol}^2= |\kappa|^2 \tilde D_{\nu_2}^2, \quad(\mathrm{IH}):\ \Delta m_{\rm sol}^2= |\kappa|^2 ({\tilde D_{\nu_2}^2-\tilde D_{\nu_1}^2}), \end{align} $

      (10)

      which can be compared to the observed value. The observed mixing matrix is defined by $ U=V^\dagger_L V_\nu $ [121], where it is parametrized by three mixing angles, i.e., i.e., $ \theta_{ij} (i,j=1,2,3; i < j) $, one CP violating Dirac phase $ \delta_{CP} $, and one Majorana phase $ \alpha_{21} $ as follows:

      $ \begin{equation} U = \begin{pmatrix} c_{12} c_{13} & s_{12} c_{13} & s_{13} {\rm e}^{-{\rm i} \delta_{CP}} \\ -s_{12} c_{23} - c_{12} s_{23} s_{13} {\rm e}^{{\rm i} \delta_{CP}} & c_{12} c_{23} - s_{12} s_{23} s_{13} {\rm e}^{{\rm i} \delta_{CP}} & s_{23} c_{13} \\ s_{12} s_{23} - c_{12} c_{23} s_{13} {\rm e}^{{\rm i} \delta_{CP}} & -c_{12} s_{23} - s_{12} c_{23} s_{13} {\rm e}^{{\rm i} \delta_{CP}} & c_{23} c_{13} \end{pmatrix} \begin{pmatrix} 1 & 0 & 0 \\ 0 & {\rm e}^{{\rm i} \frac{\alpha_{21}}{2}} & 0 \\ 0 & 0 & 1 \end{pmatrix}, \end{equation} $

      (11)

      where $ c_{ij} $ and $ s_{ij} $ stand for $ \cos \theta_{ij} $ and $ \sin \theta_{ij} $, respectively. Then, each mixing is expressed in terms of the component of U as follows:

      $ \begin{aligned}[b]& \sin^2\theta_{13}=|U_{e3}|^2,\quad \sin^2\theta_{23}=\frac{|U_{\mu3}|^2}{1-|U_{e3}|^2},\\& \sin^2\theta_{12}=\frac{|U_{e2}|^2}{1-|U_{e3}|^2}, \end{aligned} $

      (12)

      and the Majorana phase $ \alpha_{21} $ and Dirac phase $ \delta_{CP} $ are expressed in terms of the following relations:

      $ \begin{aligned}[b] & {\rm{Im}}[U^*_{e1} U_{e2}] = c_{12} s_{12} c_{13}^2 \sin \left( \frac{\alpha_{21}}{2} \right), \\& {\rm{Im}}[U^*_{e1} U_{e3}] = - c_{12} s_{13} c_{13} \sin \delta_{CP} , \end{aligned} $

      (13)

      $ \begin{aligned}[b] & {\rm{Re}}[U^*_{e1} U_{e2}] = c_{12} s_{12} c_{13}^2 \cos \left( \frac{\alpha_{21}}{2} \right), \\& {\rm{Re}}[U^*_{e1} U_{e3}] = c_{12} s_{13} c_{13} \cos \delta_{CP} , \end{aligned} $

      (14)

      where $ \alpha_{21}/2,\ \delta_{CP} $ are subtracted from π when $ \cos(\alpha_{21}/2),\ \cos\delta_{CP} $ are negative. In addition, the effective mass for the neutrinoless double beta decay is given by

      $ \begin{aligned}[b] \langle m_{ee}\rangle=&|\kappa||\tilde D_{\nu_1} \cos^2\theta_{12} \cos^2\theta_{13}+\tilde D_{\nu_2} \sin^2\theta_{12} \cos^2\theta_{13}{\rm e}^{{\rm i}\alpha_{21}}\\&+\tilde D_{\nu_3} \sin^2\theta_{13}{\rm e}^{-2{\rm i}\delta_{CP}}|, \end{aligned} $

      (15)

      where its observed value could be measured by KamLAND-Zen in future [122].

    • 2.   Neutrino masses at the fixed points
    • In general, it is difficult to show analytical predictions in arbitrary τ. However, we might be able to conduct demonstrations in specific points of τ such as fixed points. There are three important points $ \tau=i,\omega,i\infty $ according to the string theory. At these fixed points, modular forms are obtained by simple forms. Table 2 shows modular forms at $ \tau=i \infty $. In this case, a triplet modular form becomes $ (1, 0, 0)^T $ and there is a massless right-handed neutrino. Table 3 shows modular forms at $ \tau=\omega $. In this case, doublet modular forms become $ (0, 0)^T $ and all neutrinos are massless. These two cases are not suitable in our model. Table 4 shows modular forms at $ \tau=i $. Using Eqs. (9) and (10), we obtain the following relations:

      k$\boldsymbol r$$ \tau=i \infty $
      $ 2n-1 $$ \bf 2 $$ (0, 1)^T $
      $ \bf 2' $$ (0, 0)^T $
      $ \bf 2'' $$ (1, 0)^T $
      $ 2n $$ \bf 1 $$ 1 $
      $ \bf 1' $$ 0 $
      $ \bf 1'' $$ 0 $
      $ \bf 3 $$ (1, 0, 0)^T $

      Table 2.  Modular forms at $ \tau=i \infty $; note that we ignore overall factors, and n is a positive integer.

      k$\boldsymbol r$$ \tau=\omega $
      $ 6n-5 $$ \bf 2 $, $ \bf 2' $, $ \bf 2'' $$ (0, 0)^T $
      $ 6n-4 $$ \bf 1 $$ 0 $
      $ \bf 1' $$ 0 $
      $ \bf 1'' $$ 1 $
      $ \bf 3 $$ (1, \omega, -\frac12 \omega^2)^T $
      $ 6n-3 $$ \bf 2 $, $ \bf 2' $, $ \bf 2'' $$ (0, 0)^T $
      $ 6n-2 $$ \bf 1 $$ 0 $
      $ \bf 1' $$ 1 $
      $ \bf 1'' $$ 0 $
      $ \bf 3 $$ (1, -\frac12 \omega, \omega^2)^T $
      $ 6n-1 $$ \bf 2 $, $ \bf 2' $, $ \bf 2'' $$ (0, 0)^T $
      $ 6n $$ \bf 1 $$ 1 $
      $ \bf 1' $$ 0 $
      $ \bf 1'' $$ 0 $
      $ \bf 3 $$ (1, -2 \omega, -2\omega^2)^T $

      Table 3.  Modular forms at $ \tau=\omega $; note that we ignore overall factors, and n is a positive integer.

      k$\boldsymbol r$$ \tau=i $
      $ 4n-3 $$ \bf 2 $, $ \bf 2' $, $ \bf 2'' $$ ( (-1)^{\frac{7}{12}} (1 + \sqrt{3}), -\sqrt{2})^T $
      $ 4n-2 $$ \bf 1 $, $ \bf 1' $, $ \bf 1'' $$ 0 $
      $ \bf 3 $$ (1, 1 + \sqrt{3}, -2 - \sqrt{3})^T $, $ (1, -2 + \sqrt{3}, 1 - \sqrt{3})^T $
      $ 4n-1 $$ \bf 2 $, $ \bf 2' $, $ \bf 2'' $$ (-1 + (-1)^{1/6}, 1)^T $
      $ 4n $$ \bf 1 $, $ \bf 1' $, $ \bf 1'' $$ 1 $
      $ \bf 3 $$ (1, 1, 1)^T $

      Table 4.  Modular forms at $ \tau=i $; note that we ignore overall factors, and n is a positive integer.

      $ \begin{align} (\mathrm{NH}):\ \frac{\Delta m_{\rm sol}^2}{|\Delta m_{\rm atm}^2|}= \frac{{\tilde D_{\nu_2}^2}}{\tilde D_{\nu_3}^2}, \quad (\mathrm{IH}):\ \frac{\Delta m_{\rm sol}^2}{|\Delta m_{\rm atm}^2|}= \frac{{\tilde D_{\nu_2}^2-\tilde D_{\nu_1}^2}}{\tilde D_{\nu_2}^2}. \end{align} $

      (16)

      These relations are functions of $ \beta_\eta/\alpha_\eta $ and the equations have solutions if and only if the neutrino mass ordering is NH case. In the next section, we numerically check whether these analytical estimations are reasonable or not.

    • B.   Radiative seesaw

    • When $ \eta_{1,2},\ \chi $ are inert bosons, the neutrino mass matrix is induced at one-loop level via mixings among neutral components of inert bosons. Before discussing the neutrino sector, we formulate the Higgs sector. The valid soft SUSY-breaking terms to construct the neutrino mass matrix are found as follows:

      $ \begin{aligned}[b] -{\cal L}_{\rm soft} =& \mu_{BH}^2 H_1 H_2 + \mu_{B\chi}^2 Y^{(6)}_1 \chi\chi + A_a Y^{(4)}_1 H_1\eta_2 \chi\\&+ A_b Y^{(4)}_1 H_2\eta_1 \chi +m^2_{H_1}|H_1|^2+m^2_{H_2}|H_2|^2 \\&+m^2_{\eta_1}|\eta_1|^2+m^2_{\eta_2}|\eta_2|^2+m^2_{\chi}|\chi|^2 + {\rm h.c.},\end{aligned} $

      where $ m^2_{\eta_{1,2}} $ and $ m^2_{\chi} $ include the invariant coefficients $ 1/(\tau^*-\tau)^{k_{\eta_{1,2},\chi}} $.

    • 1.   Inert boson
    • Inert bosons χ, $ \eta_1 $, and $ \eta_2 $ mix each other through the soft SUSY-breaking terms of $ A_{a,b} $ and $ \mu_{B\eta} $ after the spontaneous electroweak symmetry breaking. Here, we suppose that $ \mu_{B\eta},\ A_a<<A_b $ for simplicity. Then, the mixing dominantly comes from χ and $ \eta_1 $ only. This assumption does not affect the structure of the neutrino mass matrix. Thus, the mass eigenstate is defined by

      $ \begin{align} \left[\begin{array}{c} \chi_{R,I} \\ \eta_{1_{R,I}} \\ \end{array}\right]= \left[\begin{array}{cc} c_{\theta_{R,I}} & -s_{\theta_{R,I}} \\ s_{\theta_{R,I}} & c_{\theta_{R,I}} \\ \end{array}\right] \left[\begin{array}{c} \xi_{1_{R,I}} \\ \xi_{2_{R,I}} \\ \end{array}\right], \end{align} $

      (17)

      where $ c_{\theta_{R,I}}, s_{\theta_{R,I}} $ are the shorthand notations of $ \sin\theta_{R,I} $ and $ \cos\theta_{R,I} $, respectively; $ \xi_{1,2} $ denotes the mass eigenstates for $ \chi,\eta_1 $, and their mass eigenvalues are denoted by ed by $m_{i_{R,I}}\ (i=1,2)$. Note that the mixing angle θ simultaneously diagonalizes the mass matrix of real and imaginary parts.

    • 2.   Mass matrix of $ m_\nu $
    • The active neutrino mass matrix $ m_\nu $ is induced at one-loop level as follows:

      $ \begin{aligned}[b] m_\nu =- \frac{1}{2(4\pi)^2} (y^T_\eta)_{i\alpha} (V_N)_{\alpha a} D_{N_a} (V_N^T)_{a\beta} (y_\eta)_{\beta j} \end{aligned} $

      $ \begin{aligned}[b]\quad\quad &\times \Big[ s^2_{\theta_R} f(m_{\xi_{1_R}},D_{N_a}) +c^2_{\theta_R} f(m_{\xi_{2_R}},D_{N_a}) \\& -s^2_{\theta_I} f(m_{\xi_{1_I}},D_{N_a}) -c^2_{\theta_I} f(m_{\xi_{2_I}},D_{N_a}) \Big] , \end{aligned} $

      (18)

      $ \begin{align} f(m_1,m_2)&=\int_0^1\ln\left[ x\left(\frac{m_1^2}{m_2^2}-1\right)+1 \right]. \end{align} $

      (19)

      In order to fit the atmospheric mass square difference, we extract $ \alpha_\eta $ from $ y_\eta $ and redefine $ m_\nu\equiv \alpha_\eta^2 \tilde m_\nu $. Then, we can proceed with the discussion by following the same approach as in the case of canonical seesaw by regarding $ \alpha_\eta^2 $ as κ, where κ is a parameter in the canonical seesaw model. In case of radiative seesaw, one might find that two fixed points at $ \tau=\omega,i\infty $ are not favorable owing to absence of enough degrees of freedom of non-vanishing right-handed neutrino masses, analogous to the observation in the canonical seesaw case. However, it is not possible to realize analytical predictions because of the highly complicated loop function even for a fixed point of $ \tau=i $. Thus, we cannot help relying on numerical analysis only.

    IV.   NUMERICAL ANALYSIS
    • In this section, we show numerical $ \Delta \chi^2 $ analysis for each of the cases, fitting the four reliable experimental data; $ \Delta m_{\rm sol}^2, \sin^2\theta_{13},\sin^2\theta_{23},\sin^2\theta_{12} $ in Ref. [123], where $ \Delta m^2_{\rm atm} $ is supposed to be the input value. In case of IH for the radiative seesaw model, we would not find any allowed region within $ 5\sigma $. Thus, we do not discuss this case hereafter. The dimensionful input parameters are randomly selected in the range of [$ 10^{2}-10^7 $] GeV, while the dimensionless ones are selected in the range of [$ 10^{-10}-10^{-1} $] except for τ.

    • A.   NH for the canonical seesaw model

    • Figure 1 represents NH for the canonical seesaw model. We show an allowed region of τ in the top left panel, $ \langle m_{ee}\rangle $ in terms of sum of neutrino masses $ \sum m_i $ in the top right one, Majorana phase $ \alpha_{21} $ and Dirac CP phase $\delta_{CP}$ in the bottom left one, and Dirac CP phase $\delta_{CP}$ versus sum of neutrino masses $ \sum m_i $ in the bottom right one, respectively. Each color corresponds to the range of $ \Delta\chi^2 $ values such that blue represents$ \Delta\chi^2 \leq 1 $, green represents $ 1< \Delta\chi^2 \le 4 $, yellow represents $ 4< \Delta\chi^2\le 9 $, and red represents $ 9< \Delta\chi^2 \le 25 $. These figures suggest that, within $ 5\sigma $ , 0.058 eV $\lesssim \sum{ m_i }\lesssim $ 0.06 eV, 0.001 ${\rm {eV}}\lesssim\langle m_{ee}\rangle\lesssim 0.004$ eV, any value is possible for $ \alpha_{21} $, and $ \delta_{CP} $ tends to be localized near 120 $ ^\circ $ and 240 $ ^\circ$. As analytically estimated, we have found several solutions near $ \tau=i $.

      Figure 1.  (color online) NH for the canonical seesaw model: an allowed region of τ is shown in the top left panel, $ \langle m_{ee}\rangle $ in terms of sum of neutrino masses $ \sum m_i $ in the top right one, Majorana phase $ \alpha_{21} $and Dirac CP phase $ \delta_{CP} $ in the bottom left one, and Dirac CP phase $ \delta_{CP} $ versus sum of neutrino masses $ \sum m_i $ in the bo in the bottom right one, respectively. Each color corresponds to the range of $ \Delta\chi^2 $ values such that blue represents $ \Delta\chi^2 \leq 1 $, green represents $ 1< \Delta\chi^2 \le 4 $, yellow represents $ 4< \Delta\chi^2 \le 9 $, and red represents $ 9< \Delta\chi^2 \le 25 $.

    • B.   IH for the canonical seesaw model

    • Figure 2 represents IH for the canonical seesaw model, where the legends and colors are the same as those of NH for the canonical seesaw case. These figures suggest that $ 0.098\ {\rm eV}\lesssim\sum m_i\lesssim0.102\ {\rm eV} $, $ \langle m_{ee}\rangle \simeq 0.05 $ eV, $ \alpha_{21} \simeq 0^\circ $, $ 50^\circ\lesssim\delta_{CP}\lesssim130^\circ $, and $ 230^\circ\lesssim\delta_{CP}\lesssim310^\circ $. In our analytical estimation, we have no solutions near $ \tau=i $. However, according to this figure, there would exist solutions near $ \tau=i $. Thus, we investigated the behavior of modular Yukawa functions in terms of τ and found that this function is rather sensitive to deviations from $ \tau=i $.

      Figure 2.  (color online) IH for the canonical seesaw model; the legends and colors are the same as those of NH for the canonical seesaw case.

    • C.   NH for the radiative seesaw model

    • Figure 3 represents NH for the radiative seesaw model; the legends and colors are the same as those of NH for the canonical seesaw case. These figures suggest that 0.057 ${\rm eV}\le(\sum m_i,\langle m_{ee}\rangle)\le0.06$ eV, and any values are allowed for phases. In the radiative case, there are solutions near $ \tau=i $ only in case of NH. This is what we expect from our analytical estimation.

      Figure 3.  (color online) NH for the radiative seesaw model; the legends and colors are the same as those of NH for the canonical seesaw case.

    V.   CONCLUSION AND DISCUSSION
    • We studied a double covering of modular $ A_4 $ flavor symmetry in which we constructed lepton models in cases of canonical seesaw and radiative seesaw applying irreducible doublet representations to heavier Majorana fermions that do not have $ A_4 $ symmetry. Then, we have some predictions for both cases except for the IH of radiative seesaw. From our numerical analyses, we conclude the following:

      1. In case of NH in the canonical seesaw scenario, within hin $ 5\sigma $, 0.058 ${\rm eV}\lesssim\sum m_i\lesssim0.06$ eV, 0.001 ${\rm eV}\lesssim\langle m_{ee}\rangle\lesssim0.004$ eV, any value is possible for $ \alpha_{21} $, and $ \delta_{CP} $ tends to be localized near 120 $ ^\circ $ and 240 $ ^\circ $.

      2. In case of IH in the canonical seesaw scenario, within hin $ 5\sigma $, $ 0.098\ {\rm eV}\lesssim\sum m_i\lesssim0.102\ {\rm eV} $, $ \langle m_{ee}\rangle \simeq 0.05 $ eV, $ \alpha_{21} \simeq 0^\circ $, $ 50^\circ\lesssim\delta_{CP}\lesssim130^\circ $, and $ 230^\circ\lesssim\delta_{CP}\lesssim310^\circ $.

      3. In case of NH in the radiative seesaw scenario, within $ 5\sigma $, 0.057 ${\rm eV}\le(\sum m_i,\langle m_{ee}\rangle)\le0.06$ eV, and any values are allowed for phases.

      As a future work, it would be interesting to apply this modular symmetry to both the quark and lepton sectors with the common modulus τ. We expect it to result in different predictions from $ A_4 $ modular symmetry due to especially irreducible doublet representations.

    APPENDIX A: FORMULAS IN MODULAR $ T^\prime $ FRAMEWORK
    • In this appendix, we summarize some formulas in the framework of $ T^\prime $ modular symmetry belonging to the $ SL(2,\mathbb{Z}) $ modular symmetry. The $ SL(2,Z_3) $ modular symmetry corresponds to the $ T^\prime $ modular symmetry. The modulus τ transforms as

      $ \begin{align} & \tau \longrightarrow \gamma\tau= \frac{a\tau + b}{c \tau + d}, \end{align}\tag{A1} $

      with $ \{a,b,c,d\} \in Z_3 $ satisfying $ ad-bc=1 $ and $ {\rm Im} [\tau]>0 $. The transformation of modular forms $ f(\tau) $ are given by

      $ \begin{align} & f(\gamma\tau)= (c\tau+d)^k f(\tau)\; , \; \; \gamma \in SL(2,Z_3)\; , \end{align}\tag{A2} $

      where $ f(\tau) $ denotes holomorphic functions of τ with the modular weight k.

      In a similar way, the modular transformation of a matter chiral superfield $ \phi^{(I)} $ with the modular weight $ -k_I $ is given by

      $ \begin{equation} \phi^{(I)} \to (c\tau+d)^{-k_I}\rho^{(I)}(\gamma)\phi^{(I)}, \end{equation} \tag{A3}$

      where $ \rho^{(I)}(\gamma) $ stands for a unitary matrix corresponding to $ T^\prime $ transformation. Note that the superpotential is invariant when the sum of modular weight from fields and modular form is zero and the term is a singlet under the $ T^\prime $ symmetry. It restricts a form of the superpotential, as expressed in Eq. (1).

      Modular forms are constructed on the basis of weight 1 modular form, $ Y^{(1)}_2=(Y_1, Y_2)^T $, transforming as a doublet of $ T^\prime $. Their explicit forms are written by the Dedekind eta-function $ \eta(\tau) $ with respect to τ [1, 74]:

      $ \begin{aligned}[b] Y_{1}(\tau) =& \sqrt{2}{\rm e}^{{\rm i}\frac{7\pi}{12}} \frac{\eta^3(3\tau)}{\eta(\tau)}, \\ Y_{1}(\tau) =& \sqrt{2}{\rm e}^{{\rm i}\frac{7\pi}{12}} \frac{\eta^3(3\tau)}{\eta(\tau)}. \end{aligned} $

      Modular forms of higher weight can be obtained from tensor products of $ Y^{(1)}_2 $. We enumerate some modular forms used in our analysis:

      $ \begin{align} Y_1^{(4)} = -4 Y_1^3 Y_2 - (1-i) Y_2^4, \end{align} \tag{A4}$

      $ \begin{align} Y^{(6)}_{\bf 1} &= (1-i){\rm e}^{{\rm i}\pi/6}Y_2^6 -(1+i){\rm e}^{{\rm i}\pi/6}Y_1^6 - 10 {\rm e}^{{\rm i}\pi/6}Y_1^3Y_2, \end{align}\tag{A5} $

      $ \begin{align} Y^{(2)}_3 &\equiv(f_1,f_2,f_3)^T = ( {\rm e}^{{\rm i}\pi/6}Y_2^2, \sqrt{2} {\rm e}^{{\rm i}7\pi/12} Y_1 Y_2, Y_1^2 )^T, \end{align} \tag{A6}$

      $ \begin{align} Y^{(3)}_{2} &\equiv(y_1,y_2)^T = ( 3{\rm e}^{{\rm i}\pi/6}Y_1 Y_2^2, \sqrt{2}{\rm e}^{{\rm i}5\pi/12}Y_1^3 - {\rm e}^{{\rm i}\pi/6}Y_2^3 )^T, \end{align}\tag{A7} $

      $ \begin{align} Y^{(3)}_{2^{\prime\prime}} &\equiv (y^{\prime\prime}_1, y^{\prime\prime}_2)^T = ( Y_1^3 + (1-i) Y_2^3, -3 Y_1^2 Y_2 )^T. \end{align}\tag{A8} $

Reference (123)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return