1
|
Yuk SF, Sargin I, Meyer N, Krogel JT, Beckman SP, Cooper VR. Putting error bars on density functional theory. Sci Rep 2024; 14:20219. [PMID: 39215027 PMCID: PMC11364665 DOI: 10.1038/s41598-024-69194-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2024] [Accepted: 08/01/2024] [Indexed: 09/04/2024] Open
Abstract
Predicting the error in density functional theory (DFT) calculations due to the choice of exchange-correlation (XC) functional is crucial to the success of DFT, but currently, there are limited options to estimate this a priori. This is particularly important for high-throughput screening of new materials. In this work, the structure and elastic properties of binary and ternary oxides are computed using four XC functionals: LDA, PBE-GGA, PBEsol, and vdW-DF with C09 exchange. To analyze the systemic errors inherent to each XC functional, we employed materials informatics methods to predict the expected errors. The predicted errors were also used to better the DFT-predicted lattice parameters. Our results emphasize the link between the computed errors and the electron density and hybridization errors of a functional. In essence, these results provide "error bars" for choosing a functional for the creation of high-accuracy, high-throughput datasets as well as avenues for the development of XC functionals with enhanced performance, thereby enabling the accelerated discovery and design of new materials.
Collapse
Affiliation(s)
- Simuck F Yuk
- Department of Chemistry and Life Science, United States Military Academy, West Point, NY, 10996, USA
| | - Irmak Sargin
- Department of Metallurgical and Materials Engineering, Middle East Technical University, 06800, Ankara, Turkey
| | - Noah Meyer
- Physics Department, Stanford University, Stanford, CA, 94305, USA
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Scott P Beckman
- School of Mechanical and Materials Engineering, Washington State University, Pullman, WA, 99164, USA
| | - Valentino R Cooper
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA.
| |
Collapse
|
2
|
Han S, Lee J, Han S, Moosavi SM, Kim J, Park C. Design of New Inorganic Crystals with the Desired Composition Using Deep Learning. J Chem Inf Model 2023; 63:5755-5763. [PMID: 37683188 DOI: 10.1021/acs.jcim.3c00935] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/10/2023]
Abstract
New solid-state materials have been discovered using various approaches from atom substitution in density functional theory (DFT) to generative models in machine learning. Recently, generative models have shown promising performance in finding new materials. Crystal generation with deep learning has been applied in various methods to discover new crystals. However, most generative models can only be applied to materials with specific elements or generate structures with random compositions. In this work, we developed a model that can generate crystals with desired compositions based on a crystal diffusion variational autoencoder. We generated crystal structures for 14 compositions of three types of materials in different applications. The generated structures were further stabilized using DFT calculations. We found the most stable structures in the existing database for all but one composition, even though eight compositions among them were not in the data set trained in a crystal diffusion variational autoencoder. This substantiates the prospect of the generation of an extensive range of compositions. Finally, 205 unique new crystal materials with energy above hull <100 meV/atom were generated. Moreover, we compared the average formation energy of the crystals generated from five compositions, two of which were hypothetical, with that of traditional methods like atom substitution and a generative model. The generated structures had lower formation energy than those of other models, except for one composition. These results demonstrate that our approach can be applied stably in various fields to design stable inorganic materials based on machine learning.
Collapse
Affiliation(s)
- Seunghee Han
- Department of Chemical and Biomolecular Engineering, Korea Advanced Institute of Science and Technology, Daejeon 34141, Republic of Korea
| | - Jaewan Lee
- LG AI Research, ISC, 30, Magokjungang 10-ro, Gangseogu, Seoul 07796, Republic of Korea
| | - Sehui Han
- LG AI Research, ISC, 30, Magokjungang 10-ro, Gangseogu, Seoul 07796, Republic of Korea
| | - Seyed Mohamad Moosavi
- Department of Chemical Engineering & Applied Chemistry, University of Toronto, Ontario M5S 3E5, Canada
| | - Jihan Kim
- Department of Chemical and Biomolecular Engineering, Korea Advanced Institute of Science and Technology, Daejeon 34141, Republic of Korea
| | - Changyoung Park
- LG AI Research, ISC, 30, Magokjungang 10-ro, Gangseogu, Seoul 07796, Republic of Korea
| |
Collapse
|
3
|
Afify AA, Abd El-Fattah ZM, Esmail MSM, El-Bahnasawy HH. Comparative Mössbauer and first principles calculations for selected iron oxides and ferrites nanoparticles. MATERIALS TODAY COMMUNICATIONS 2023; 35:106193. [DOI: 10.1016/j.mtcomm.2023.106193] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
|
4
|
Saikia D, Alam M, Bera J, Betal A, Gandi AN, Sahu S. A First‐Principles Study on ABBr
3
(A = Cs, Rb, K, Na; B = Ge, Sn) Halide Perovskites for Photovoltaic Applications. ADVANCED THEORY AND SIMULATIONS 2022. [DOI: 10.1002/adts.202200511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Dibyajyoti Saikia
- Department of Physics Indian Institute of Technology Jodhpur Jodhpur 342037 India
| | - Mahfooz Alam
- Department of Metallurgical and Materials Engineering Indian Institute of Technology Jodhpur Jodhpur 342037 India
| | - Jayanta Bera
- Department of Physics Indian Institute of Technology Jodhpur Jodhpur 342037 India
| | - Atanu Betal
- Department of Physics Indian Institute of Technology Jodhpur Jodhpur 342037 India
| | - Appala Naidu Gandi
- Department of Metallurgical and Materials Engineering Indian Institute of Technology Jodhpur Jodhpur 342037 India
| | - Satyajit Sahu
- Department of Physics Indian Institute of Technology Jodhpur Jodhpur 342037 India
| |
Collapse
|
5
|
Rahman MH, Jubair M, Rahaman MZ, Ahasan MS, Ostrikov KK, Roknuzzaman M. RbSnX 3 (X = Cl, Br, I): promising lead-free metal halide perovskites for photovoltaics and optoelectronics. RSC Adv 2022; 12:7497-7505. [PMID: 35424654 PMCID: PMC8982351 DOI: 10.1039/d2ra00414c] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2022] [Accepted: 02/26/2022] [Indexed: 11/21/2022] Open
Abstract
Lead (Pb) free metal halide perovskites by atomistic design are of strong interest to photovoltaics and optoelectronics industries because of the pressing need to resolve Pb-related toxicity and instability challenges. In this study, structural, mechanical, electronic, and optical properties of Pb-free RbSnX3 (X = Cl, Br, I) perovskites have been evaluated by using ab initio density functional theory (DFT) calculations. The computed elastic constants suggest that the Rb-based halide perovskites are mechanically stable and highly ductile, making them suitable as flexible thin films in optoelectronic devices. Besides, the investigated electronic band structures reveal that the RbSnX3 compounds are direct bandgap semiconductors, suitable for photovoltaic and optoelectronic applications. Furthermore, several optical parameters such as dielectric functions, reflectivity, photon absorptions, refractive index, optical conductivity, and loss functions have been investigated and the results predict the excellent optoelectronic efficiency of RbSnX3. Also, the computed mechanical and optical properties of RbSnX3 (X = Cl, Br, I) have been compared with the previously studied CsBX3 (B = Ge, Sn, Pb; X = Cl, Br, I) phases, revealing that the Rb-based perovskites are extremely ductile and possess excellent light absorption and optical conductivity compared to the Cs-based perovskites. Importantly, RbSnI3 shows superior ductility, absorption coefficient, and optical conductivity compared to the CsBX3 (B = Ge, Sn, Pb; X = Cl, Br, I) perovskites. Superior absorption at the ultraviolet region of RbSnI3 holds great promise of this perovskite to be used in next-generation ultraviolet photodetectors.
Collapse
Affiliation(s)
- Md Habibur Rahman
- Department of Mechanical Engineering, Bangladesh University of Engineering and Technology Dhaka 1000 Bangladesh
| | - Md Jubair
- Department of Physics, Rajshahi University of Engineering and Technology Rajshahi 6204 Bangladesh
| | - Md Zahidur Rahaman
- School of Materials Science and Engineering, University of New South Wales Sydney NSW 2052 Australia
| | - Md Shamim Ahasan
- Department of Physics, University of Rajshahi Rajshahi 6205 Bangladesh
| | - Kostya Ken Ostrikov
- School of Chemistry and Physics and QUT Centre for Materials Science, Queensland University of Technology (QUT) Brisbane QLD 4000 Australia
| | - Md Roknuzzaman
- School of Chemistry and Physics and QUT Centre for Materials Science, Queensland University of Technology (QUT) Brisbane QLD 4000 Australia
- School of Physics, University of New South Wales Sydney NSW 2052 Australia
| |
Collapse
|
6
|
How WB, Wang B, Chu W, Tkatchenko A, Prezhdo OV. Significance of the Chemical Environment of an Element in Nonadiabatic Molecular Dynamics: Feature Selection and Dimensionality Reduction with Machine Learning. J Phys Chem Lett 2021; 12:12026-12032. [PMID: 34902248 DOI: 10.1021/acs.jpclett.1c03469] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Using supervised and unsupervised machine learning (ML) on features generated from nonadiabatic (NA) molecular dynamics (MD) trajectories under the classical path approximation, we demonstrate that mutual information with the NA Hamiltonian can be used for feature selection and model simplification. Focusing on CsPbI3, a popular metal halide perovskite, we observe that the chemical environment of a single element is sufficient for predicting the NA Hamiltonian. The conclusion applies even to Cs, although Cs does not contribute to the relevant wave functions. Interatomic distances between Cs and I or Pb and the octahedral tilt angle are the most important features. We reduce a typical 360-parameter ML force-field model to just a 12-parameter NA Hamiltonian model, while maintaining a high NA-MD simulation quality. Because NA-MD is a valuable tool for studying excited state processes, overcoming its high computational cost through simple ML models will streamline NA-MD simulations and expand the ranges of accessible system size and simulation time.
Collapse
Affiliation(s)
- Wei Bin How
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 637371 Singapore
| | - Bipeng Wang
- Department of Chemical Engineering, University of Southern California, Los Angeles, California 90089, United States
| | - Weibin Chu
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Oleg V Prezhdo
- Department of Chemical Engineering, University of Southern California, Los Angeles, California 90089, United States
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
7
|
Khan R, Ur Rahman K, Zhang Q, Rahman AU, Azam S, Dahshan A. The effect of substitutional doping of Yb 2+on structural, electronic, and optical properties of CsCa X3( X: Cl, Br, I) phosphors: a first-principles study. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 34:065502. [PMID: 34727528 DOI: 10.1088/1361-648x/ac3583] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/08/2021] [Accepted: 11/02/2021] [Indexed: 06/13/2023]
Abstract
Using first-principles calculations, the effects of Yb2+substitutional doping on structural, electronic, and optical properties of a series of perovskite compounds CsCaX3(X: Cl, Br, I), have been investigated. We employed generalized gradient approximation (GGA) and HSE hybrid functional to study the electronic and optical properties. A series of pristine CsCaX3(X: Cl, Br, I) is characterized as a non-magnetic insulator with indirect bandgap perovskite materials. These phosphor materials are suitable candidates for doping with lanthanide series elements to tune their electronic bandgaps according to our requirements because of their wide bandgaps. The calculated electronic bandgaps of CsCaX3(X: Cl, Br, I) are 3.7 eV (GGA) and 4.5 eV (HSE) for CsCaI3, 4.5 eV (GGA) and 5.3 eV (HSE) for CsCaBr3, and 5.4 eV (GGA) and 6.4 eV (HSE) for CsCaCl3. According to formation energies, the Yb2+doped at the Ca-site is thermodynamically more stable as compared to all possible atomic sites. The electronic band structures show that the Yb2+doping induces defective states within the bandgaps of pristine CsCaX3(X: Cl, Br, I). As a result, the Yb2+doped CsCaX3(X: Cl, Br, I) become the direct bandgap semiconductors. The defective states above the valence band maximum are produced due to thef-orbital of the Yb atom. The impurity states near the conduction band minimum are induced due to the major contribution ofd-orbital of the Yb atom and the minor contribution ofs-orbital of the Cs atom. The real and imaginary parts of the dielectric function, optical reflectivity, electron energy loss spectrum, extinction coefficient, and refractive index of pristine and Yb2+doped CsCaX3(X: Cl, Br, I) were studied. The optical dispersion results of dielectric susceptibility closely match their relevant electronic structure and align with previously reported theoretical and experimental data. We conclude that the Yb2+doped CsCaX3(X: Cl, Br, I) are appealing candidates for optoelectronic devices.
Collapse
Affiliation(s)
- Rashid Khan
- School of Energy and Power Engineering, Xi'an Jiaotong University (XJTU), 28 Xianning W.Rd., Xi'an 710049, People's Republic of China
| | - Kaleem Ur Rahman
- Department of Physics, Quaid-i-Azam University, Islamabad, Pakistan
| | - Qingmin Zhang
- School of Energy and Power Engineering, Xi'an Jiaotong University (XJTU), 28 Xianning W.Rd., Xi'an 710049, People's Republic of China
| | - Altaf Ur Rahman
- Department of Physics, Riphah International University, Lahore, Pakistan
| | - Sikander Azam
- Department of Physics, Riphah International University, Islamabad, Pakistan
| | - Alaa Dahshan
- Department of Physics, Faculty of Science, King Khalid University, PO Box 9004, Abha, Saudi Arabia
- Department of Physics, Faculty of Science, Port Said University, Port Said, Egypt
| |
Collapse
|
8
|
Chen K, Qi K, Zhou T, Yang T, Zhang Y, Guo Z, Lim CK, Zhang J, Žutic I, Zhang H, Prasad PN. Water-Dispersible CsPbBr 3 Perovskite Nanocrystals with Ultra-Stability and its Application in Electrochemical CO 2 Reduction. NANO-MICRO LETTERS 2021; 13:172. [PMID: 34383132 PMCID: PMC8360258 DOI: 10.1007/s40820-021-00690-8] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/22/2021] [Accepted: 06/24/2021] [Indexed: 05/03/2023]
Abstract
Thanks to the excellent optoelectronic properties, lead halide perovskites (LHPs) have been widely employed in high-performance optoelectronic devices such as solar cells and light-emitting diodes. However, overcoming their poor stability against water has been one of the biggest challenges for most applications. Herein, we report a novel hot-injection method in a Pb-poor environment combined with a well-designed purification process to synthesize water-dispersible CsPbBr3 nanocrystals (NCs). The as-prepared NCs sustain their superior photoluminescence (91% quantum yield in water) for more than 200 days in an aqueous environment, which is attributed to a passivation effect induced by excess CsBr salts. Thanks to the ultra-stability of these LHP NCs, for the first time, we report a new application of LHP NCs, in which they are applied to electrocatalysis of CO2 reduction reaction. Noticeably, they show significant electrocatalytic activity (faradaic yield: 32% for CH4, 40% for CO) and operation stability (> 350 h).
Collapse
Affiliation(s)
- Keqiang Chen
- Key Laboratory of Optoelectronic Devices and Systems of Ministry of Education and Guangdong Province, Institute of Microscale Optoelectronics, Shenzhen University, Shenzhen, 518060, People's Republic of China
- Institute for Lasers, Photonics, and Biophotonics, Department of Chemistry, University at Buffalo, State University of New York , Buffalo, NY, 14260, USA
| | - Kun Qi
- Key Laboratory of Optoelectronic Devices and Systems of Ministry of Education and Guangdong Province, Institute of Microscale Optoelectronics, Shenzhen University, Shenzhen, 518060, People's Republic of China
| | - Tong Zhou
- Department of Physics, University at Buffalo, State University of New York, Buffalo, NY, 14260, USA
| | - Tingqiang Yang
- Key Laboratory of Optoelectronic Devices and Systems of Ministry of Education and Guangdong Province, Institute of Microscale Optoelectronics, Shenzhen University, Shenzhen, 518060, People's Republic of China
| | - Yupeng Zhang
- Key Laboratory of Optoelectronic Devices and Systems of Ministry of Education and Guangdong Province, Institute of Microscale Optoelectronics, Shenzhen University, Shenzhen, 518060, People's Republic of China.
| | - Zhinan Guo
- Key Laboratory of Optoelectronic Devices and Systems of Ministry of Education and Guangdong Province, Institute of Microscale Optoelectronics, Shenzhen University, Shenzhen, 518060, People's Republic of China
| | - Chang-Keun Lim
- Institute for Lasers, Photonics, and Biophotonics, Department of Chemistry, University at Buffalo, State University of New York , Buffalo, NY, 14260, USA
- Department of Chemical and Materials Engineering, School of Engineering, Nazarbayev University, Nur-Sultan City, 010000, Kazakhstan
| | - Jiayong Zhang
- Jiangsu Key Laboratory of Micro and Nano Heat Fluid Flow Technology and Energy Application, School of Physical Science and Technology, Suzhou University of Science and Technology, Suzhou, Jiangsu, 215009, People's Republic of China
| | - Igor Žutic
- Department of Physics, University at Buffalo, State University of New York, Buffalo, NY, 14260, USA
| | - Han Zhang
- Key Laboratory of Optoelectronic Devices and Systems of Ministry of Education and Guangdong Province, Institute of Microscale Optoelectronics, Shenzhen University, Shenzhen, 518060, People's Republic of China.
| | - Paras N Prasad
- Institute for Lasers, Photonics, and Biophotonics, Department of Chemistry, University at Buffalo, State University of New York , Buffalo, NY, 14260, USA.
| |
Collapse
|
9
|
Yu X, Liu Z, Yang X, Wang Y, Zhang J, Duan J, Liu L, Tang Q. Crystal-Plane Controlled Spontaneous Polarization of Inorganic Perovskite toward Boosting Triboelectric Surface Charge Density. ACS APPLIED MATERIALS & INTERFACES 2021; 13:26196-26203. [PMID: 34047548 DOI: 10.1021/acsami.1c05796] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Triboelectric generators (TENGs) have been extensively studied as a new energy for low cost and the universally applicable prospect. Meanwhile, perovskites have been applied in TENG and show a good performance in view of high carrier mobility, long life and dielectric properties. The asymmetry structure of the orthogonal phase CsPbBr3 perovskite endows it with ferroelectric property and induces the misalignment of the positive and negative charge centers. Herein, the surface energy of halogen doped inorganic CsPbX3 (X = Cl-, Br-) perovskites are theoretically investigated by density functional theory (DFT) calculation, the crystal polarizability of pristine CsPbBr3 is improved from 0.47 Ry a.u. to 0.52 Ry a.u. (CsPbCl3), indicating the polarizability of CsPbCl3 is higher than CsPbBr3. In addition, the build-in electric field (Ebuild-in) of perovskite materials can be enhanced by the spontaneous polarization and the aligned dipoles in the Ebuild-in could further improve the tribo-electrostatic electric field by retaining more triboelectric surface charges. In the end, CsPbCl3 achieved a power of 3.06 W m-2 compared to the power of 1.34 W m-2 of CsPbBr3. This work focuses on the regulation of crystal planes using spontaneous polarization of perovskite toward achieving a high built-in electric field for enhancing triboelectric surface charge density.
Collapse
Affiliation(s)
- Xueping Yu
- College of Information Science and Technology, Jinan University, Guangzhou 510632, PR China
- State Key Laboratory of Electronic Thin Films and Integrated Devices, University of Electronic Science & Technology of China, Chengdu, Sichuan 610054, PR China
| | - Zhongzhe Liu
- State Key Laboratory of Electronic Thin Films and Integrated Devices, University of Electronic Science & Technology of China, Chengdu, Sichuan 610054, PR China
| | - Xiya Yang
- College of Information Science and Technology, Jinan University, Guangzhou 510632, PR China
| | - Yudi Wang
- College of Information Science and Technology, Jinan University, Guangzhou 510632, PR China
| | - Jihua Zhang
- State Key Laboratory of Electronic Thin Films and Integrated Devices, University of Electronic Science & Technology of China, Chengdu, Sichuan 610054, PR China
| | - Jialong Duan
- College of Information Science and Technology, Jinan University, Guangzhou 510632, PR China
| | - Liming Liu
- Zhongshan Branch of State Key Laboratory of Electronic Thin Films and Integrated Devices, University of Electronic Science and Technology of China, Zhongshan Institute, Zhongshan, Guangdong 528402, PR China
| | - Qunwei Tang
- College of Information Science and Technology, Jinan University, Guangzhou 510632, PR China
| |
Collapse
|
10
|
Yang Y, Hou C, Liang TX. Energetic and electronic properties of CsPbBr 3 surfaces: a first-principles study. Phys Chem Chem Phys 2021; 23:7145-7152. [PMID: 33876080 DOI: 10.1039/d0cp04893c] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Surface properties of all-inorganic halide perovskites play a crucial role in determining optoelectronic performance of these materials. We investigate the surface energies and electronic structures of cubic CsPbBr3 surfaces systematically using density functional theory (DFT) methods. We calculate the surface phase diagrams of low-index surfaces of CsPbBr3, i.e., (100), (110), (111) surfaces. We found that nonpolar (100) surfaces are more stable than polar (110) and (111) surfaces. The nonpolar CsBr-terminated (100) surface shows the best stability, which is attributed to the effect of surface relaxation and high ionicity of the surface layer. The electronic structures reveal that charge transfer to compensate the polarity raises the energy of polar surfaces, which makes polar surfaces unstable. Furthermore, we found that the modulation of surface chemical composition provides an effective way to compensate polarity and thus make polar surfaces of CsPbBr3 stable. Our results provide physical insights into understanding and further enhancing the surface stability of all-inorganic halide perovskites. This would be helpful in promoting the advancement of all-inorganic halide perovskite-based materials and devices.
Collapse
Affiliation(s)
- Yi Yang
- College of Rare Earths and Faculty of Materials, Metallurgy and Chemistry, Jiangxi University of Science and Technology, Ganzhou, China.
| | | | | |
Collapse
|