Route to Stabilize Cubic Gauche Polynitrogen to Ambient Conditions via Surface Saturation by Hydrogen
-
Abstract
Cubic gauche polynitrogen (cg-N) is an attractive high-energy density material. However, high-pressure synthesized cg-N will decompose at low pressure and cannot exist under ambient conditions. Here, the stabilities of cg-N surfaces with and without saturations at different pressures and temperatures are systematically investigated based on first-principles calculations and molecular dynamics simulations. Pristine surfaces at 0 GPa are very brittle and will decompose at 300 K, especially (110) surface will collapse completely just after structural relaxation, whereas the decompositions of surfaces can be suppressed by applying pressure, indicating that surface instability causes the cg-N decomposition at low pressure. Due to the saturation of dangling bonds and transferring electrons to the surfaces, saturation with H can stabilize surfaces under ambient conditions, while it is impossible for OH saturation to occur solely from obtaining electrons from surfaces. This suggests that polynitrogen is more stable in an acidic environment or when the surface is saturated with less electronegative adsorbates. -
The triple bond in nitrogen molecule (N≡N with bonding energy of 954 kJ/mol) is one of the strongest chemical bonds of all molecules, and the bonding energy will be decrease and total energy will increase significantly if the N≡N is transformed into the double bond (N=N with bonding energy of 418 kJ/mol) and the single bond (N–N with bonding energy of 160 kJ/mol).[1] On the other hand, an enormous amount of energy will be released when the N–N and N=N decompose into the N≡N. Therefore, polymerized nitrogen containing N–N or N=N has been attracted lots of interests for potential applications as high-energy density material (HEDM).[2–26]
To break N≡N and synthesize polymerized nitrogen, extreme synthetic conditions are needed, and high pressure is a popular way.[27] Existing works have shown that within the low pressure range, nitrogen exists in the form of molecule crystals.[28–33] The first polymerized nitrogen with fully single bonding, named as cubic gauche polynitrogen (cg-N), was theoretically predicted to appear at 50 GPa,[29] and its energy density is about five times larger than that of TNT.[34] Furthermore, many kinds of polynitrogens (such as chains,[35,36] layers,[37,38] cage,[39] and polynitrogen[40,41]) were theoretically proposed to be thermodynamically stable at high pressure. In 2004, the cg-N was successfully synthesized under the extreme conditions of 110 GPa and 2000 K.[42] Following cg-N, layered polymeric nitrogen (LP-N),[43] black phosphorus nitrogen (BP-N),[44] and hexagonal layered polymeric nitrogen (HLP-N)[45] have been synthesized at the pressures of 120 GPa, 150 GPa, and 180 GPa, respectively. Nevertheless, the synthetic pressure of cg-N is the smallest among these synthesized polymeric all-nitrogen materials.
Importantly, there are nonnegative frequencies in the phonon spectra of cg-N at 0 GPa,[46] indicating that it potentially can be quenched to the ambient pressure. However, the cg-N becomes unstable and begins to decompose with pressure decreasing to 42 GPa at room temperature or to 25 GPa at temperature of 140 K,[42] and the same phenomenon also occurs to other polymeric all-nitrogen materials. For example, LP-N, HLP-N, and BP-N can only be stabilized up to 52 GPa, 66 GPa, and 48 GPa at room temperature, respectively,[43–45] which hinder their practical applications as HEDMs. In order to preserve cg-N to the ambient pressure, at present, the important tasks are to find out its decomposition mechanism at low pressure and to find a way to stabilize it under ambient conditions.
It is generally believed that decomposition of materials will start from the weakest bonds, which may exist on the surfaces containing energy unfavorable dangling bonds.[47] A typical example is that reconstructions are usually observed on materials’ surfaces because of the instability of free surfaces, such as Au (111) and Si (100).[48,49] The weakest bond of cg-N may also exist on the surface. Actually, the stability is improved when the cg-N is placed in the carbon nanotubes, indicating that surface stability may strongly affect the stability of cg-N.[50] However, so far, the surface stabilities of cg-N have not been well illustrated.
In this work, we investigate the stabilities of cg-N low-index surfaces at different pressures (0 GPa, 40 GPa, 80 GPa, 100 GPa, 120 GPa) under different saturated conditions (H, OH, CH3), using the first-principles method as implemented in the Vienna ab initio approximation package (VASP) combined with molecular dynamics (MD) method. Our results show that the instabilities of pristine surfaces cause the cg-N decomposition at low pressure. The saturation of OH cannot improve the surface stabilities, while all investigated surfaces can be stabilized under ambient conditions by H-saturation due to the electron transform from H to surfaces.
The crystal structure of cg-N is shown in Fig. S1(a) in the Supplementary Materials, where the local coordination environments of N is similar to that of NH3. Each N atom connects with other three neighbors via N–N single bonds, and two electrons form a lone pair of electrons. The N–N single bonds are connected to each other, forming a cubic gauche configuration.[29] All sub-surfaces obtained from the Miller indices of (100), (110), and (111) are shown in Figs. S1(b)–S1(h). By exposing different surfaces to vacuum, two [(100a), (100b)] and four [(110a), (110b), (110c), (110d)] kinds of surfaces are constructed based on the Miller indices of (100) and (110), respectively, and only one surface (111a) is obtained in the case of (111) Miller indices. The corresponding free surfaces constructed are shown in Fig. S2(a).
The configurations of pristine surfaces after structural relaxation at 0 GPa are shown in Fig. 1(a). The (100a) and (111a) surfaces are stable, keeping their original structures, while linear N3 groups (similar to the azide ion) are formed and adsorbed on (110c) surface, giving rise to a stable reconfigured (110c) [i.e., r-(110c)] surface. After releasing one nitrogen molecule, the (100b) surface transfers to the stable (100a) surface. Similar phenomenon occurs in the (110a) surface, and it will transfer to the r-(110c) surface after releasing two nitrogen molecules. Interestingly, the (110b) surface breaks down layer by layer via releasing nitrogen molecules and decomposes completely under the processes of structure relaxation. The (110d) surface is also dissociative completely as well, breaking down into zigzag nitrogen chains and nitrogen molecules. Therefore, the (110b) and (110d) surfaces are the weakest among those of all investigated sub-surfaces. In the following, we investigate the stability of (110b) and (110d) surfaces under high pressure.
Fig. Fig. 1. Stabilities of pristine surfaces. (a) Structures of pristine surfaces after structural relaxations at 0 GPa. (b) Optimized structures of (110d) surfaces at pressures ranging from 0 GPa to 120 GPa, and He atoms are filled in vacuum zone. (c) Structures of pristine (100a), (110c), and (111a) surfaces after MD simulations at 300 K by using the supercell with 192 N atoms. The atoms located in the blue shadow area are fixed. For structural relaxations filled with He at high pressure, all nitrogen atoms and helium atoms are allowed to relax. To clearly show the atomic structure of surfaces, the vacuum between surfaces are not fully presented, and it is indicated by the dotted lines. Blue balls represent nitrogen atoms and orange represents helium atoms.To apply hydrostatic pressure on the surfaces, helium gas employed to fill the vacuum zone. The relaxed structures of (110d) surfaces at high pressure are shown in Fig. 1(b). We can find that at the pressures of 100 GPa and 120 GPa closed to the conditions of experimentally synthesized cg-N, the (110d) surface is stable. Furthermore, the structure is maintained with pressure decreasing to 60 GPa, and slight surface decomposition appears at 40 GPa and 20 GPa, and totally structural decomposition occurs at 0 GPa. As shown in Fig. S3 in the Supplementary Materials, the stability of (110b) surface is also significantly enhanced by applying high pressures larger than 20 GPa. The results agree well with experimental observations, where the cg-N can only be quenched to 25 GPa at low temperature.[42] Based on the above results, we can make the following assumption about the process of cg-N decomposition: Since pressure can enhance the stability of surfaces, the surface decompositions will be suppressed under high-pressure conditions, resulting in a stable cg-N, and surfaces eventually become unstable after depressurization, disintegrating the whole structure.
MD simulations were performed to study the stabilities of the (100a), (111a), and r-(110c) surfaces at 0 GPa and 300 K. To allow for sufficient relaxation of the surface, the supercell for MD simulation is increased 2 × 2 times as much as that used during structure relaxation. All subsequent MD simulations were conducted using this larger supercells. The corresponding z-axis mean square displacements (MSDs) are shown in Fig. S4. We can find that the MSDs of both (100a) and (111a) surfaces at 300 K increase noticeably, and corresponding N–N networks are broken down from surface after MD simulations [Fig. 1(c)], indicating that (100a) and (111a) surfaces cannot exist stably at 300 K. In contrast, the MSD of r-(110c) surface at 300 K is almost a constant, and its structure can exist stably. It is easier to stabilize cg-N at low temperatures, which is consistent with the previous experimental result.[42,51] The above results show that except r-(110c) surface, all investigated pristine surfaces are not stable under ambient conditions.
To stabilize the high-pressure synthesized cg-N under the ambient conditions and to promote its practical applications, it is necessary to find a way to stabilize the unstable surfaces [especially (110b) and (110d) surfaces], and to enhance the finite-temperature stabilities of all investigated surfaces. The high stability of r-(110c) surface, where the dangling bonds on the surface are saturated with the linear N3 groups, indicates that saturation of the dangling bonds is important for enhancing the stabilities of cg-N surfaces. Therefore, saturations of dangling bonds at 0 GPa with H, OH, and CH3 are investigated, and the corresponding saturated surfaces are named as surface-H, surface-OH, and surface-CH3, respectively (see Fig. S2 in the Supplementary Materials).
The stabilities of surfaces with (without) saturations under the static condition and at 300 K are summarized in Fig. 2. The OH saturation does not improve the surface stabilities, since the stabilities of (100b), (110d), and (111a) surfaces are enhanced but (100a) and r-(110c) surfaces changes from stable to unstable, while the others remains unchanged. Furthermore, the CH3 saturation provides finite stabilizing effect that all surfaces are improved except (100a)-CH3 surface, and three surfaces [(100b), (110b) and (110c)] cannot exist stably at 300 K. However, with H saturation, all investigated surfaces are stable at 300 K, indicating that the cg-N can be stabilized to ambient conditions if its surfaces are saturated with hydrogen. Please note that for two stable pristine surfaces [(100a) and (111a)], the H adsorption energies of (100a)-H and (111a)-H are −0.794 eV and −1.278 eV, respectively, and H saturation is an energy-favorable process.
Fig. Fig. 2. The stabilities of surfaces without and with H, OH, and CH3 saturations after structure relaxation and MD simulations at 300 K. The brackets on the y-axis tags show the average number of electrons (in units of e) transferred from adsorbates to surface, where positive (negative) values present electrons transfer from adsorbates to surfaces (from surfaces to adsorbates). The r-(110c) means the (110c) surface with surface reconstruction. SR: structural relaxation.We further run MD simulations for surface-H at higher temperatures of 500 K, 750 K, 1000 K, and 1250 K for illustrating the maximum temperature that can be stabilized by H saturation, which is important for its practical application as HEDM. The corresponding z-axis MSD at temperatures before and after the structure loses stability, as shown in Fig. 3(a). We can find that H saturation can improve the stability of all surfaces to at least 750 K and even more than half of the surfaces [(100a), (100b), (110a) and (110d)] can be stabilized to 1000 K. The stabilities of surface-H are summarized in Fig. 3(b) including that of pristine surfaces used for comparison, indicating that H saturation can enhance the stabilities of cg-N surfaces significantly.
Fig. Fig. 3. Stabilities of surfaces with H saturation. (a) The MSD of surface-H at the temperature before (blue line) and after (red line) structure collapse. (b) Stability of pristine surfaces and surface-H in structural relaxation and MD simulations at finite temperatures up to 1000 K. Different levels of rings represent the different stabilities, the outer the ring is, the more stable it is. The r-(110c) means the (110c) surface with surface reconstruction.It can be observed that at high temperatures, the decomposition products of surface-H include nitrogen gas and nitrogen-hydrogen compounds. Nitrogen-hydrogen compounds observed include N2H3, N3H, N3H2, N4H, N4H2, N4H3, N5H, N6H2, and N8H5. Additionally, N4H chains are observed for (110d)-H, and a small amount of H2 is also detected in (110a)-H.
Interestingly, reconstruction occurring in the pristine (110c) surface with the formation of N3 groups [Fig. 1(a)] is suppressed by H saturation. As shown in Fig. 4, after structural relaxation, the reconstruction does not occur in the H-saturated (110c). With temperature increasing to 300 K and 500 K, the reconstruction is triggered by temperature and occurs in half of (110c)-H surface, and it further extends to the whole surface at 750 K. Note that the reconstruction of (110c) surface cannot be suppressed by the CH3 saturation, and CH3 will be decomposed by the reconstruction (Fig. S5). This phenomenon makes us further convince that H saturation is a better way to stabilize the cg-N surfaces.
Fig. Fig. 4. Structures of (110c)-H surfaces at different temperatures. (a) Structure of (110c)-H after structure relaxation. (b) Structures of (110c)-H after MD simulations at 300 K. (c) Structures of (110c)-H after MD simulations at 500 K. (d) Structures of (110c)-H after MD simulations at 750 K. Reconstruction does not occur after the structural relaxation, while half of surface has reconstruction after MD simulations at 300 K and 500 K. At 750 K, surface reconstruction is completed. To facilitate identification, reconfigured N atoms are marked with yellow circles.Based on the above discussions, we can find that although the dangling bonds are saturated in both surface-OH and surface-H, their stabilities show large difference. To understand this phenomenon, we further investigate the charger transfer between saturated groups and surfaces relying on the Bader charge analysis (Table S1 in the Supplementary Materials), and the average number of electrons transferring from saturated groups into nitrogen atoms is shown in the brackets on the y-axis tag in Fig. 2. Positive (negative) values mean electrons transfer from saturated groups to surfaces (from surfaces to saturated groups). For surface-H, surface-CH3 and surface-OH, the charge transfer is 1.39e, 1.19e, and −0.28e, respectively, and the stability of surfaces increases with the number of electrons on the surfaces (surface-OH < surface-CH3 < surface-H). The electronegativity of OH (3.49)[52] is larger than that of N (3.04),[53] and electrons would transfer from N to OH. However, the opposite process occurs when hydrogen atoms are adsorbed on the surface, because the electronegativity of H (2.20)[52] is notably smaller.
To further rationalize the charge transfer behaviors of surface-OH and surface-H, the electron density difference of (111a)-H and (111a)-OH surfaces are calculated as shown in Fig. 5, which confirms charge transfers from H to N in (111a)-H surface and N to OH in (111a)-OH surface. The results clearly show that electron transfer to nitrogen is beneficial to enhance the surface stabilities, while the opposite process leads to the decreasing stability. Therefore, although the dangling bonds are saturated in surface-OH, the OH saturation takes electrons from the surfaces and cannot promote their stabilities. However, H adsorption not only saturates the dangling bonds but also transfers electrons to the surfaces, giving rise to a very high stability of cg-N surfaces.
Fig. Fig. 5. (a) Electron density difference (EDD) of (111a)-H surface. (b) EDD of (111a)-OH surface. The curves are plotted along the z axis of the structures, which are shown in the insets. The local integral curve shown in red presents the electron gain (positive values) and loss (negative values), and the integral curve shown in blue is the integral of the local integral curve. Isosurface in the structures is set to 0.003 e/Bohr3, and orange (blue) isosurface presents positive (negative) values of electron density difference.In summary, based on the first-principles and molecular dynamics method, stabilities of cg-N surfaces are systematically investigated, and effects of pressure, temperature, and saturations with H, OH, and CH3 on the surface stabilities are presented. For pristine cases, only the reconfigured (110c) surfaces with linear N3 groups adsorption can exist at 300 K, and other surfaces are very brittle, especially the (110b) and (110d) surfaces will collapse totally just after structural relaxation. However, the decompositions of surfaces will be suppressed at high pressure, indicating that surface instability causes the cg-N instability at low pressure. Interestingly, by H saturation, all of investigated surfaces can be stabilized not only under ambient conditions but also at 750 K, since H atoms can saturate dangling bonds combined with transferring electrons to the surfaces. OH saturation cannot improve the stabilities of surfaces, because it gets electrons from the surfaces.
We propose the following plan to realize hydrogen surface-saturated cg-N: The most part of the sample chamber of the diamond anvil cell is filled with nitrogen gas, and small part is filled with a suitable hydrogen-rich material (e.g., ammonia borane) for providing hydrogen source by heating. In the process of synthesizing cg-N, the N2-filled part is heated by laser heating, while the part filled with hydrogen-rich material remains unheated. After cg-N is synthesized, the hydrogen-rich material can be heated to generate hydrogen, which can saturate the surface of cg-N. Besides the high-pressure synthesis, our results can also provide effective guidance for future synthesis of cg-N under atmospheric pressure. For example, an acidic environment is helpful to improve the stabilities of cg-N and its relatives. Furthermore, materials, which can be more conducive to transferring electrons to neighbors, such as graphene and carbon nanotubes with a band gap as small as possible, can be used as assistant to improve the stabilities of polymeric nitrogen.
Acknowledgements: This work was supported by the National Natural Science Foundation of China (Grant No. U2030114), and CASHIPS Director’s Fund (Grant No. YZJJ202207-CX). The calculations were partly performed at Center for Computational Science of CASHIPS, the ScGrid of Supercomputing Center and Computer Network Information Center of Chinese Academy of Sciences, and the Hefei Advanced Computing Center. -
References
[1] Eremets M I, Hemley R J, Mao H, Gregoryanz E 2001 Nature 411 170 doi: 10.1038/35075531[2] Ma Y M, Oganov A R, Li Z W, Xie Y, Kotakoski J 2009 Phys. Rev. Lett. 102 065501 doi: 10.1103/PhysRevLett.102.065501[3] Sun J, Martinez-Canales M, Klug D D, Pickard C J, Needs R J 2013 Phys. Rev. Lett. 111 175502 doi: 10.1103/PhysRevLett.111.175502[4] Nordlund K, Krasheninnikov A, Juslin N, Nord J, Albe K 2004 Europhys. Lett. 65 400 doi: 10.1209/epl/i2003-10087-1[5] Abou-Rachid H, Hu A, Timoshevskii V, Song Y, Lussier L S 2008 Phys. Rev. Lett. 100 196401 doi: 10.1103/PhysRevLett.100.196401[6] Adeleke A A, Greschner M J, Majumdar A, Wan B, Liu H, Li Z, Gou H, Yao Y 2017 Phys. Rev. B 96 224104 doi: 10.1103/PhysRevB.96.224104[7] Alemany M M G, Martins J L 2003 Phys. Rev. B 68 024110 doi: 10.1103/PhysRevB.68.024110[8] Katzke H, Tolédano P 2008 Phys. Rev. B 78 064103 doi: 10.1103/PhysRevB.78.064103[9] Zhang L J, Wang Y C, Lv J, Ma Y M 2017 Nat. Rev. Mater. 2 17005 doi: 10.1038/natrevmats.2017.5[10] Yao Y S, Tse J S, Tanaka K 2008 Phys. Rev. B 77 052103 doi: 10.1103/PhysRevB.77.052103[11] Li Q, Rellán-Piñeiro M, Almora-Barrios N, Garcia-Ratés M, Remediakis I N, López N 2017 Nanoscale 9 13089 doi: 10.1039/C7NR03889E[12] Goncharov A F, Gregoryanz E, Mao H, Liu Z, Hemley R J 2000 Phys. Rev. Lett. 85 1262 doi: 10.1103/PhysRevLett.85.1262[13] Boates B, Bonev S A 2009 Phys. Rev. Lett. 102 015701 doi: 10.1103/PhysRevLett.102.015701[14] Boates B, Bonev S A 2011 Phys. Rev. B 83 174114 doi: 10.1103/PhysRevB.83.174114[15] Gregoryanz E, Goncharov A F, Sanloup C, et al. 2007 J. Chem. Phys. 126 184505 doi: 10.1063/1.2723069[16] Eremets M I, Gavriliuk A G, Serebryanaya N R, Trojan I A, Dzivenko D A, Boehler R, Mao H K, Hemley R J 2004 J. Chem. Phys. 121 11296 doi: 10.1063/1.1814074[17] Zahariev F, Dudiy S V, Hooper J, Zhang F, Woo T K 2006 Phys. Rev. Lett. 97 155503 doi: 10.1103/PhysRevLett.97.155503[18] Zahariev F, Hooper J, Alavi S, Zhang F, Woo T K 2007 Phys. Rev. B 75 140101 doi: 10.1103/PhysRevB.75.140101[19] Hirshberg B, Gerber R B, Krylov A I 2014 Nat. Chem. 6 52 doi: 10.1038/nchem.1818[20] Zhang J, Zeng Z, Lin H, Li Y 2014 Sci. Rep. 4 4358 doi: 10.1038/srep04358[21] Lei L, Tang Q Q, Zhang F, Liu S, B W, Zhou C 2020 Chin. Phys. Lett. 37 068101 doi: 10.1088/0256-307X/37/6/068101[22] Li Y W, Feng X L, Liu H Y, Hao J, Redfern S A T, Lei W W, Liu D, Ma Y M 2018 Nat. Commun. 9 722 doi: 10.1038/s41467-018-03200-4[23] Salke N P, Xia K, Fu S, Zhang Y, Greenberg E, Prakapenka V B, Liu J, Sun J, Lin J F 2021 Phys. Rev. Lett. 126 065702 doi: 10.1103/PhysRevLett.126.065702[24] Liu Y, Su H, Niu C, Wang X, Zhang J, Ge Z, Li Y 2020 Chin. Phys. B 29 106201 doi: 10.1088/1674-1056/aba9bd[25] Liu L L, Zhang S T, Zhang H J 2022 Chin. Phys. Lett. 39 056102 doi: 10.1088/0256-307X/39/5/056102[26] Ding C, Wang J, Han Y, Yuan J, Gao H, Sun J 2022 Chin. Phys. Lett. 39 036101 doi: 10.1088/0256-307X/39/3/036101[27] Lei L, Pu M, Feng L, Qi L, Zhang L 2018 Chin. J. High Press. Phys. 32 020102 in Chinese doi: 10.11858/gywlxb.20170672[28] Yakub L N 2016 Low Temp. Phys. 42 1 doi: 10.1063/1.4940225[29] Mailhiot C, Yang L H, McMahan A K 1992 Phys. Rev. B 46 14419 doi: 10.1103/PhysRevB.46.14419[30] Fan C M, Liu S, Liu J Y, et al. 2022 Chin. Phys. Lett. 39 026401 doi: 10.1088/0256-307X/39/2/026401[31] Frost M, Howie R T, Dalladay-Simpson P, Goncharov A F, Gregoryanz E 2016 Phys. Rev. B 93 024113 doi: 10.1103/PhysRevB.93.024113[32] Liu S, Pu M, Tang Q, Zhang F, Wu B, Lei L 2020 Solid State Commun. 310 113843 doi: 10.1016/j.ssc.2020.113843[33] Pu M F, Liu S, Lei L, Zhang F, Feng L H, Qi L, Zhang L L 2019 Solid State Commun. 298 113645 doi: 10.1016/j.ssc.2019.113645[34] Xia K, Yuan J, Zheng X, Liu C, Gao H, Wu Q, Sun J 2019 J. Phys. Chem. Lett. 10 6166 doi: 10.1021/acs.jpclett.9b02383[35] Mattson W D, Sanchez-Portal D, Chiesa S, Martin R M 2004 Phys. Rev. Lett. 93 125501 doi: 10.1103/PhysRevLett.93.125501[36] Oganov A R, Glass C W 2006 J. Chem. Phys. 124 244704 doi: 10.1063/1.2210932[37] Zahariev F, Hu A, Hooper J, Zhang F, Woo T 2005 Phys. Rev. B 72 214108 doi: 10.1103/PhysRevB.72.214108[38] Pickard C J, Needs R J 2009 Phys. Rev. Lett. 102 125702 doi: 10.1103/PhysRevLett.102.125702[39] Wang X L, Wang Y C, Miao M S, Zhong X, Lv J, Cui T, Li J F, Chen L, Pickard C J, Ma Y M 2012 Phys. Rev. Lett. 109 175502 doi: 10.1103/PhysRevLett.109.175502[40] Zhang J, Wang X, Yang K, Cheng Y, Zeng Z 2018 Sci. Rep. 8 13144 doi: 10.1038/s41598-018-31355-z[41] Zhang J, Niu C, Zhang H, Zhao J, Wang X, Zeng Z 2021 J. Phys. Chem. Lett. 12 5731 doi: 10.1021/acs.jpclett.1c01181[42] Eremets M I, Gavriliuk A G, Trojan I A, Dzivenko D A, Boehler R 2004 Nat. Mater. 3 558 doi: 10.1038/nmat1146[43] Tomasino D, Kim M, Smith J, Yoo C S 2014 Phys. Rev. Lett. 113 205502 doi: 10.1103/PhysRevLett.113.205502[44] Ji C, Adeleke A A, Yang L X, et al. 2020 Sci. Adv. 6 eaba9206 doi: 10.1126/sciadv.aba9206[45] Laniel D, Geneste G, Weck G, Mezouar M, Loubeyre P 2019 Phys. Rev. Lett. 122 066001 doi: 10.1103/PhysRevLett.122.066001[46] Zhang T, Zhang S, Chen Q, Peng L M 2006 Phys. Rev. B 73 094105 doi: 10.1103/PhysRevB.73.094105[47] Christe K O 2007 Prop. Explos. Pyrotech. 32 194 doi: 10.1002/prep.200700020[48] Fan X L, Fang X L, Ran R X, Lau W M 2014 Physica E 59 248 doi: 10.1016/j.physe.2014.01.029[49] Saranin A A, Zotov A V, Ignatovich K V, Lifshits V G, Numata T, Kubo O, Tani H, Katayama M, Oura K 1997 Phys. Rev. B 56 1017 doi: 10.1103/PhysRevB.56.1017[50] Benchafia E M, Yao Z, Yuan G, Chou T, Piao H, Wang X, Iqbal Z 2017 Nat. Commun. 8 930 doi: 10.1038/s41467-017-01083-5[51] Liu S, Pu M, Zhang F, Lei L 2019 J. Light Scattering 31 235 in Chinese doi: 10.13883/j.issn1004-5929.201903004[52] Pauling L 1932 J. Am. Chem. Soc. 54 3570 doi: 10.1021/ja01348a011[53] Inamoto N, Masuda S 1982 Chem. Lett. 11 1003 doi: 10.1246/cl.1982.1003 -
Related Articles
[1] CHEN Jian-Bin, FENG Tai-Fu, GAO Tie-Jun. Kaluza–Klein Corrections to the μ Anomalous Magnetic Moment in the Appelquist–Cheng–Dobrescu Model [J]. Chin. Phys. Lett., 2012, 29(9): 091101. doi: 10.1088/0256-307X/29/9/091101 [2] Eser Olgar. Exact Solution of Klein--Gordon Equation by Asymptotic Iteration Method [J]. Chin. Phys. Lett., 2008, 25(6): 1939-1942. [3] M. Kocak. Bound State Solutions of Klein--Gordon Equation with the Kratzer Potential [J]. Chin. Phys. Lett., 2007, 24(2): 315-317. [4] DUAN Yi-Shi, WU Shao-Feng. Regular Magnetic Monopole from Generalized ’t Hooft Tensor [J]. Chin. Phys. Lett., 2006, 23(11): 2932-2935. [5] B.Gönül. A Search on the Klein--Gordon Equation [J]. Chin. Phys. Lett., 2006, 23(10): 2640-2643. [6] ZHANG Min-Cang, WANG Zhen-Bang. Exact Solutions of the Klein--Gordon Equation with a New Anharmonic Oscillator Potential [J]. Chin. Phys. Lett., 2005, 22(12): 2994-2996. [7] JING Ji-Liang, YAN Mu-Lin. Statistical Entropy of the Kaluza-Klein Black Hole from the Horizon Conformal Field Theory [J]. Chin. Phys. Lett., 2001, 18(3): 331-333. [8] NI Guang-Jiong, GUAN Hong, ZHOU Wei-Min, YAN Jun. Antiparticle in the Light of Einstein-Podolsky-Rosen Paradox and Klein Paradox [J]. Chin. Phys. Lett., 2000, 17(6): 393-395. [9] WEN Wan-xin, ZHONG Ji-quan, JIN Gen-ming. Giant Monopole Resonance with Boltzmann-Nordheim-Vlasov Model [J]. Chin. Phys. Lett., 1996, 13(10): 724-726. [10] SHEN Yougen. COSMICAL WAVE FUNCTION FOR THE COUPLING SCALARS FIELD IN THE THEORY OF KALUZA-KLEIN [J]. Chin. Phys. Lett., 1989, 6(1): 43-46.