1
|
Liu Y, Li X, Qian C, Müller S, Skrzyńska K, Gallego-Parra S, Garbarino G, Chumakov A, Kupenko I, Wu X. Pressure-Induced Hydrogen Bond Symmetrization in Babingtonite. Inorg Chem 2025; 64:7570-7579. [PMID: 40179372 DOI: 10.1021/acs.inorgchem.5c00522] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/05/2025]
Abstract
Hydrogen bond symmetrization and electronic spin transition of iron under pressure are key factors in understanding the physical properties and behaviors of Fe- and H-bearing minerals. However, the coexistence of these two effects may often result in inconsistent interpretations. Babingtonite is characterized by prominent hydrogen bonds and diverse iron oxidation states within its crystal structure. These features make it an ideal subject for studying the interactions between hydrogen bond symmetrization and spin transitions of iron. Here, we conducted high-pressure infrared absorption spectroscopy, single-crystal X-ray diffraction, and Mössbauer spectroscopy experiments on babingtonite to investigate the behaviors of hydrogen bonds and spin states of iron. Pressure-induced hydrogen bond symmetrization occurs at approximately 22 GPa in babingtonite, while no spin transition of iron is detected up to 40 GPa. Our results show that the notable deformation of the unit cell and the kink in the trend of the hydrogen-bonded O-O distance with increasing pressure provide insightful evidence for hydrogen bond symmetrization. The symmetrical rigid configuration of the hydrogen bond leads to increased distortion in the neighboring Fe3+-O octahedron. These findings will contribute to the determination of hydrogen bond symmetrization through crystallographic structural analysis, partly overcoming the challenge of accurately locating hydrogen positions using X-ray diffraction.
Collapse
Affiliation(s)
- Yungui Liu
- College of Gems and Materials Technology, Hebei GEO University, Shijiazhuang 050031, China
- State key laboratory of geological processes and mineral resources, China University of Geosciences, Wuhan 430074, China
- Hebei Key Laboratory of Green Development of Rock and Mineral Materials, Hebei GEO University, Shijiazhuang 050031, China
| | - Xiang Li
- European Synchrotron Radiation Facility, The European Synchrotron, Grenoble 38043, France
- Institut für Mineralogie, Universität Münster, Münster D-48149, Germany
| | - Cheng Qian
- State key laboratory of geological processes and mineral resources, China University of Geosciences, Wuhan 430074, China
| | - Susanne Müller
- European Synchrotron Radiation Facility, The European Synchrotron, Grenoble 38043, France
| | - Katarzyna Skrzyńska
- European Synchrotron Radiation Facility, The European Synchrotron, Grenoble 38043, France
- Faculty of Natural Sciences, Institute of Earth Sciences, University of Silesia, Sosnowiec 41-200, Poland
| | - Samuel Gallego-Parra
- European Synchrotron Radiation Facility, The European Synchrotron, Grenoble 38043, France
| | - Gaston Garbarino
- European Synchrotron Radiation Facility, The European Synchrotron, Grenoble 38043, France
| | - Aleksandr Chumakov
- European Synchrotron Radiation Facility, The European Synchrotron, Grenoble 38043, France
| | - Ilya Kupenko
- European Synchrotron Radiation Facility, The European Synchrotron, Grenoble 38043, France
| | - Xiang Wu
- State key laboratory of geological processes and mineral resources, China University of Geosciences, Wuhan 430074, China
| |
Collapse
|
2
|
Wang Y, Luo R, Chen J, Zhou X, Wang S, Wu J, Kang F, Yu K, Sun B. Proton Collective Quantum Tunneling Induces Anomalous Thermal Conductivity of Ice under Pressure. PHYSICAL REVIEW LETTERS 2024; 132:264101. [PMID: 38996295 DOI: 10.1103/physrevlett.132.264101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2023] [Revised: 03/18/2024] [Accepted: 05/20/2024] [Indexed: 07/14/2024]
Abstract
Proton tunneling is believed to be nonlocal in ice, but its range has been shown to be limited to only a few molecules. Here, we measured the thermal conductivity of ice under pressure up to 50 GPa and found it increases with pressure until 20 GPa but decreases at higher pressures. We attribute this nonmonotonic thermal conductivity to the collective tunneling of protons at high pressures, supported by large-scale quantum molecular dynamics simulations. The collective tunneling loops span several picoseconds in time and are as large as nanometers in space, which match the phonon periods and wavelengths, leading to strong phonon scattering at high pressures. Our results show direct evidence of global quantum motion existing in high-pressure ice and provide a new perspective to understanding the coupling between phonon propagation and atomic tunneling.
Collapse
|
3
|
Kim T, Wei X, Chariton S, Prakapenka VB, Ryu YJ, Yang S, Shim SH. Stability of hydrides in sub-Neptune exoplanets with thick hydrogen-rich atmospheres. Proc Natl Acad Sci U S A 2023; 120:e2309786120. [PMID: 38109550 PMCID: PMC10756278 DOI: 10.1073/pnas.2309786120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2023] [Accepted: 10/27/2023] [Indexed: 12/20/2023] Open
Abstract
Many sub-Neptune exoplanets have been believed to be composed of a thick hydrogen-dominated atmosphere and a high-temperature heavier-element-dominant core. From an assumption that there is no chemical reaction between hydrogen and silicates/metals at the atmosphere-interior boundary, the cores of sub-Neptunes have been modeled with molten silicates and metals (magma) in previous studies. In large sub-Neptunes, pressure at the atmosphere-magma boundary can reach tens of gigapascals where hydrogen is a dense liquid. A recent experiment showed that hydrogen can induce the reduction of Fe[Formula: see text] in (Mg,Fe)O to Fe[Formula: see text] metal at the pressure-temperature conditions relevant to the atmosphere-interior boundary. However, it is unclear whether Mg, one of the abundant heavy elements in the planetary interiors, remains oxidized or can be reduced by H. Our experiments in the laser-heated diamond-anvil cell found that heating of MgO + Fe to 3,500 to 4,900 K (close to or above their melting temperatures) in an H medium leads to the formation of Mg[Formula: see text]FeH[Formula: see text] and H[Formula: see text]O at 8 to 13 GPa. At 26 to 29 GPa, the behavior of the system changes, and Mg-H in an H fluid and H[Formula: see text]O were detected with separate FeH[Formula: see text]. The observations indicate the dissociation of the Mg-O bond by H and subsequent production of hydride and water. Therefore, the atmosphere-magma interaction can lead to a fundamentally different mineralogy for sub-Neptune exoplanets compared with rocky planets. The change in the chemical reaction at the higher pressures can also affect the size demographics (i.e., "radius cliff") and the atmosphere chemistry of sub-Neptune exoplanets.
Collapse
Affiliation(s)
- Taehyun Kim
- School of Earth and Space Exploration, Arizona State University, Tempe, AZ85287
| | - Xuehui Wei
- School of Earth and Space Exploration, Arizona State University, Tempe, AZ85287
| | - Stella Chariton
- GeoSoilEnviroCARS, Center for Advanced Radiation Sources, University of Chicago, Argonne, IL60439
| | - Vitali B. Prakapenka
- GeoSoilEnviroCARS, Center for Advanced Radiation Sources, University of Chicago, Argonne, IL60439
| | - Young-Jay Ryu
- GeoSoilEnviroCARS, Center for Advanced Radiation Sources, University of Chicago, Argonne, IL60439
| | - Shize Yang
- Eyring Materials Center, Arizona State University, Tempe, AZ85287
| | - Sang-Heon Shim
- School of Earth and Space Exploration, Arizona State University, Tempe, AZ85287
| |
Collapse
|
4
|
Chuai M, Yang J, Tan R, Liu Z, Yuan Y, Xu Y, Sun J, Wang M, Zheng X, Chen N, Chen W. Theory-Driven Design of a Cationic Accelerator for High-Performance Electrolytic MnO 2 -Zn Batteries. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2022; 34:e2203249. [PMID: 35766725 DOI: 10.1002/adma.202203249] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2022] [Revised: 06/03/2022] [Indexed: 06/15/2023]
Abstract
Aqueous electrolytic MnO2 -Zn batteries are considered as one of the most promising energy-storage devices for their cost effectiveness, high output voltage, and safety, but their electrochemical performance is limited by the sluggish kinetics of cathodic MnO2 /Mn2+ and anodic Zn/Zn2+ reactions. To overcome this critical challenge, herein, a cationic accelerator (CA) strategy is proposed based on the prediction of first-principles calculations. Poly(vinylpyrrolidone) is utilized as a model to testify the rational design of the CA strategy. It manifests that the CA effectively facilitates rapid cations migration in electrolyte and adequate charge transfer at electrode-electrolyte interface, benefiting the deposition/dissolution processes of both Mn2+ and Zn2+ cations to simultaneously improve kinetics of cathodic MnO2 /Mn2+ and anodic Zn/Zn2+ reactions. The resulting MnO2 -Zn battery regulated by CA exhibits large reversible capacities of 455 mAh g-1 and 3.64 mAh cm-2 at 20 C, as well as a long lifespan of 2000 cycles with energy density retention of 90%, achieving one of the best overall performances in the electrolytic MnO2 -Zn batteries. This comprehensive work integrating theoretical prediction with experimental studies provides opportunities to the development of high-performance energy-storage devices.
Collapse
Affiliation(s)
- Mingyan Chuai
- Guangdong Research Center for Interfacial Engineering of Functional Materials, College of Materials Science and Engineering, Shenzhen University, Shenzhen, 518060, China
- Department of Applied Chemistry, School of Chemistry and Materials Science, Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Jinlong Yang
- Guangdong Research Center for Interfacial Engineering of Functional Materials, College of Materials Science and Engineering, Shenzhen University, Shenzhen, 518060, China
| | - Rui Tan
- Department of Chemical Engineering, Imperial College London, London, SW72AZ, UK
| | - Zaichun Liu
- Guangdong Research Center for Interfacial Engineering of Functional Materials, College of Materials Science and Engineering, Shenzhen University, Shenzhen, 518060, China
- Department of Applied Chemistry, School of Chemistry and Materials Science, Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Yuan Yuan
- Guangdong Research Center for Interfacial Engineering of Functional Materials, College of Materials Science and Engineering, Shenzhen University, Shenzhen, 518060, China
- Department of Applied Chemistry, School of Chemistry and Materials Science, Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Yan Xu
- Department of Applied Chemistry, School of Chemistry and Materials Science, Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Jifei Sun
- Department of Applied Chemistry, School of Chemistry and Materials Science, Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Mingming Wang
- Department of Applied Chemistry, School of Chemistry and Materials Science, Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Xinhua Zheng
- Department of Applied Chemistry, School of Chemistry and Materials Science, Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Na Chen
- Department of Applied Chemistry, School of Chemistry and Materials Science, Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, China
| | - Wei Chen
- Department of Applied Chemistry, School of Chemistry and Materials Science, Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, China
| |
Collapse
|
5
|
Adachi Y, Koga K. Structure and phase behavior of high-density ice from molecular-dynamics simulations with the ReaxFF potential. J Chem Phys 2020; 153:114501. [PMID: 32962394 DOI: 10.1063/5.0016565] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
We report a molecular dynamics simulation study of dense ice modeled by the reactive force field (ReaxFF) potential, focusing on the possibility of phase changes between crystalline and plastic phases as observed in earlier simulation studies with rigid water models. It is demonstrated that the present model system exhibits phase transitions, or crossovers, among ice VII and two plastic ices with face-centered cubic (fcc) and body-centered cubic (bcc) lattice structures. The phase diagram derived from the ReaxFF potential is different from those of the rigid water models in that the bcc plastic phase lies on the high-pressure side of ice VII and does the fcc plastic phase on the low-pressure side of ice VII. The phase boundary between the fcc and bcc plastic phases on the pressure, temperature plane extends to the high-temperature region from the triple point of ice VII, fcc plastic, and bcc plastic phases. Proton hopping, i.e., delocalization of a proton, along between two neighboring oxygen atoms in dense ice is observed for the ReaxFF potential but only at pressures and temperatures both much higher than those at which ice VII-plastic ice transitions are observed.
Collapse
Affiliation(s)
- Yuji Adachi
- Graduate School of Natural Sciences, Okayama University, Okayama 700-8530, Japan
| | - Kenichiro Koga
- Department of Chemistry, Okayama University, Okayama 700-8530, Japan
| |
Collapse
|
6
|
Abstract
Helmholtz energy of ice VII–X is determined in a pressure regime extending to 450 GPa at 300 K using local-basis-functions in the form of b-splines. The new representation for the equation of state is embedded in a physics-based inverse theory framework of parameter estimation. Selected pressures as a function of volume from 14 prior experimental studies and two theoretical studies constrain the behavior of Helmholtz energy. Separately measured bulk moduli, not used to construct the representation, are accurately replicated below about 20 GPa and above 60 GPa. In the intermediate range of pressure, the experimentally determined moduli are larger and have greater scatter than values predicted using the Helmholtz representation. Although systematic error in the determination of elastic moduli is possible and likely, the alternative hypothesis is a slow relaxation time associated with changes in proton mobility or the ice VII to X transition. A correlation is observed between anomalies in the pressure derivative of the predicted bulk modulus and previously suggested higher-order phase transitions. Improved determinations of elastic properties at high pressure would allow refinement of the current equation of state. More generally, the current method of data assimilation is broadly applicable to other materials in high-pressure studies and for investigations of planetary interiors.
Collapse
|
7
|
Li F, Wang Y, Li Z, Men Z, Sun C. Enhanced Stimulated Raman Scattering by a Pressure-Controlled Shock Wave in Liquid Water. J Phys Chem Lett 2019; 10:4812-4816. [PMID: 31390212 DOI: 10.1021/acs.jpclett.9b01956] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Stimulated Raman scattering (SRS) is observed using a Nd:YAG laser in liquid water at both forward and backward directions under different pressures. The spectra at atmospheric pressure and high pressure exhibit different characteristic features. For high pressure, the main SRS peak (about 3400 cm-1) of liquid water shifts to low frequency. Interestingly, a new peak is observed in both directions. The position of the new peak is lower than that at atmospheric pressure, which belongs to strong hydrogen bonds. Especially, a low peak is obtained at around 3140 cm-1 in the backward direction at 400 MPa, indicating the formation of an ice-like structure. In addition, the normalized SRS intensity of high pressure is higher than that of atmospheric pressure. These results indicate that high pressure can significantly enhance the SRS of water molecules. The enhancement mechanism is attributed to the long duration and slightly slow velocity of the shock wave induced by high pressure.
Collapse
Affiliation(s)
- Fabing Li
- Coherent Light and Atomic and Molecular Spectroscopy Laboratory, College of Physics, Jilin University, Changchun 130012, China
| | - Ying Wang
- Coherent Light and Atomic and Molecular Spectroscopy Laboratory, College of Physics, Jilin University, Changchun 130012, China
| | - Zhanlong Li
- Coherent Light and Atomic and Molecular Spectroscopy Laboratory, College of Physics, Jilin University, Changchun 130012, China
| | - Zhiwei Men
- Coherent Light and Atomic and Molecular Spectroscopy Laboratory, College of Physics, Jilin University, Changchun 130012, China
| | - Chenglin Sun
- Coherent Light and Atomic and Molecular Spectroscopy Laboratory, College of Physics, Jilin University, Changchun 130012, China
- Key Laboratory of Physics and Technology for Advanced Batteries, College of Physics, Jilin University, Changchun 130012, China
| |
Collapse
|
8
|
Ikeda T. First principles isothermal-isobaric centroid molecular dynamics simulation of high pressure ices. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.01.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
9
|
Shin D, Hwang J, Jhe W. Ice-VII-like molecular structure of ambient water nanomeniscus. Nat Commun 2019; 10:286. [PMID: 30655538 PMCID: PMC6336866 DOI: 10.1038/s41467-019-08292-0] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2018] [Accepted: 12/25/2018] [Indexed: 11/09/2022] Open
Abstract
Structural transformations originating from diverse rearrangements of the hydrogen bonding in water create various phases. Although most phases have been well investigated down to the molecular level, the molecular structure of the nanomeniscus, a ubiquitous form of nanoscale water in nature, still remains unexplored. Here, we demonstrate that the water nanomeniscus exhibits the stable, ice-VII-like molecular structure in ambient condition. Surface-enhanced Raman spectroscopy on trace amounts of water, confined in inter-nanoparticle gaps, shows a narrowed tetrahedral peak at 3340 cm-1 in the OH-stretching band as well as a lattice-vibrational mode at 230 cm-1. In particular, the ice-VII-like characteristics are evidenced by the spectral independence with respect to temperature variations and differing surface types including the material, size and shape of nanoparticles. Our results provide un unambiguous identification of the molecular structure of nanoconfined water, which is useful for understanding the molecular aspects of water in various nanoscale, including biological, environments.
Collapse
Affiliation(s)
- Dongha Shin
- Center for 0D Nanofluidics, Institute of Applied Physics, Department of Physics and Astronomy, Seoul National University, Seoul, 08826, Republic of Korea
| | - Jonggeun Hwang
- Center for 0D Nanofluidics, Institute of Applied Physics, Department of Physics and Astronomy, Seoul National University, Seoul, 08826, Republic of Korea
| | - Wonho Jhe
- Center for 0D Nanofluidics, Institute of Applied Physics, Department of Physics and Astronomy, Seoul National University, Seoul, 08826, Republic of Korea.
| |
Collapse
|
10
|
Direct observation of symmetrization of hydrogen bond in δ-AlOOH under mantle conditions using neutron diffraction. Sci Rep 2018; 8:15520. [PMID: 30341340 PMCID: PMC6195538 DOI: 10.1038/s41598-018-33598-2] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2018] [Accepted: 10/01/2018] [Indexed: 11/09/2022] Open
Abstract
At ambient pressure, the hydrogen bond in materials such as ice, hydrates, and hydrous minerals that compose the Earth and icy planets generally takes an asymmetric O-H···O configuration. Pressure significantly affects this configuration, and it is predicted to become symmetric, such that the hydrogen is centered between the two oxygen atoms at high pressure. Changes of physical properties of minerals relevant to this symmetrization have been found; however, the atomic configuration around this symmetrization has remained elusive so far. Here we observed the pressure response of the hydrogen bonds in the aluminous hydrous minerals δ-AlOOH and δ-AlOOD by means of a neutron diffraction experiment. We find that the transition from P21nm to Pnnm at 9.0 GPa, accompanied by a change in the axial ratios of δ-AlOOH, corresponds to the disorder of hydrogen bond between two equivalent sites across the center of the O···O line. Symmetrization of the hydrogen bond is observed at 18.1 GPa, which is considerably higher than the disorder pressure. Moreover, there is a significant isotope effect on hydrogen bond geometry and transition pressure. This study indicates that disorder of the hydrogen bond as a precursor of symmetrization may also play an important role in determining the physical properties of minerals such as bulk modulus and seismic wave velocities in the Earth's mantle.
Collapse
|
11
|
Meier T, Petitgirard S, Khandarkhaeva S, Dubrovinsky L. Observation of nuclear quantum effects and hydrogen bond symmetrisation in high pressure ice. Nat Commun 2018; 9:2766. [PMID: 30018359 PMCID: PMC6050302 DOI: 10.1038/s41467-018-05164-x] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2018] [Accepted: 06/15/2018] [Indexed: 11/19/2022] Open
Abstract
Hydrogen bond symmetrisations in H-bonded systems triggered by pressure-induced nuclear quantum effects (NQEs) is a long-known concept but experimental evidence in high-pressure ices has remained elusive with conventional methods. Theoretical works predicted quantum-mechanical tunneling of protons within water ices to occur at pressures above 30 GPa, and the H-bond symmetrisation transition to occur above 60 GPa. Here we used 1H-NMR on high-pressure ice up to 97 GPa, and demonstrate that NQEs govern the behavior of the hydrogen bonded protons in ice VII already at significantly lower pressures than previously expected. A pronounced tunneling mode was found to be present up to the highest pressures of 97 GPa, well into the stability field of ice X, where NQEs are not anticipated in a fully symmetrised H-bond network. We found two distinct transitions in the NMR shift data at about 20 GPa and 75 GPa attributed to the step-wise symmetrisation of the H-bond.
Collapse
Affiliation(s)
- Thomas Meier
- Bayerisches Geoinstitut, Bayreuth University, Universitätsstraße 30, 95447, Bayreuth, Germany.
| | - Sylvain Petitgirard
- Bayerisches Geoinstitut, Bayreuth University, Universitätsstraße 30, 95447, Bayreuth, Germany
| | - Saiana Khandarkhaeva
- Bayerisches Geoinstitut, Bayreuth University, Universitätsstraße 30, 95447, Bayreuth, Germany
| | - Leonid Dubrovinsky
- Bayerisches Geoinstitut, Bayreuth University, Universitätsstraße 30, 95447, Bayreuth, Germany
| |
Collapse
|
12
|
Ikeda T. First principles centroid molecular dynamics simulation of high pressure ices. J Chem Phys 2018; 148:102332. [DOI: 10.1063/1.5003055] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Takashi Ikeda
- Synchrotron Radiation Research Center, Quantum Beam Science Research Directorate (QuBS), National Institutes for Quantum and Radiological Science and Technology (QST), 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan
| |
Collapse
|