1
|
Wang Z, Zhang A, Hua T, Chen X, Zhu M, Guo Z, Song Y, Yang G, Li S, Feng J, Li M, Yan W. Revealing the interaction forms between Hg(II) and group types (-Cl, -CN, -NH 2, -OH, -COOH) in functionalized Poly(pyrrole methane)s for efficient mercury removal. ENVIRONMENTAL POLLUTION (BARKING, ESSEX : 1987) 2024; 351:124049. [PMID: 38692386 DOI: 10.1016/j.envpol.2024.124049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2024] [Revised: 04/21/2024] [Accepted: 04/23/2024] [Indexed: 05/03/2024]
Abstract
To explore the impact of different functional groups on Hg(II) adsorption, a range of poly(pyrrole methane)s functionalized by -Cl, -CN, -NH2, -OH and -COOH were synthesized and applied to reveal the interaction between different functional groups and mercury ions in water, and the adsorption mechanism was revealed through combined FT-IR, XPS, and DFT calculations. The adsorption performance can be improved to varying degrees by the incorporation of functional groups. Among them, the oxygen-containing functional groups (-OH and -COOH) exhibit stronger affinity for Hg(II) and can increase the adsorption capacity from 180 mg g-1 to more than 1400 mg g-1 at 318 K, with distribution coefficient (Kd) exceeding 105 mL g-1. The variations in the capture and immobilization capabilities of functionalized poly(pyrrole methane)s predominantly stem from the unique interactions between their functional groups and mercury ions. In particular, oxygen-containing -OH and -COOH effectively capture Hg(OH)2 through hydrogen bonding, and further deprotonate to form the -O-Hg-OH and -COO-Hg-OH complexes which are more stable than those obtained from other functionalized groups. Finally, the ecological safety has been fully demonstrated through bactericidal and bacteriostatic experiments to prove the functionalized poly(pyrrole methane)s can be as an environmentally friendly adsorbent for purifying contaminated water.
Collapse
Affiliation(s)
- Zhenyu Wang
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, 710049, China
| | - Aijing Zhang
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, 710049, China
| | - Tingyu Hua
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, 710049, China
| | - Xin Chen
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, 710049, China
| | - Mengyuan Zhu
- International Research Center for Renewable Energy (IRCRE), State Key Laboratory of Multiphase Flow in Power Engineering (MFPE), Xi'an Jiaotong University, Xi'an, Shaanxi, 710049, China
| | - Ziyu Guo
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, 710049, China
| | - Yanna Song
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, 710049, China
| | - Guorui Yang
- School of Chemistry, Engineering Research Center of Energy Storage Materials and Devices, Ministry of Education, "Four Joint Subjects One Union" School-Enterprise Joint Research Center for Power Battery Recycling & Circulation Utilization Technology, Xi'an Jiaotong University, Xi'an, 710049, China
| | - Shanshan Li
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, 710049, China
| | - Jiangtao Feng
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, 710049, China; International Research Center for Renewable Energy (IRCRE), State Key Laboratory of Multiphase Flow in Power Engineering (MFPE), Xi'an Jiaotong University, Xi'an, Shaanxi, 710049, China.
| | - Mingtao Li
- International Research Center for Renewable Energy (IRCRE), State Key Laboratory of Multiphase Flow in Power Engineering (MFPE), Xi'an Jiaotong University, Xi'an, Shaanxi, 710049, China
| | - Wei Yan
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, 710049, China
| |
Collapse
|
2
|
Florez E, Zapata-Escobar AD, Ferraro F, Ibargüen Becerra C, Chamorro Y, Maldonado AF. Coordination of Mercury(II) in Water Promoted over Hydrolysis in Solvated Clusters [Hg(H 2O) 1-6] (aq)2+: Insights from Relativistic Effects and Free Energy Analysis. J Phys Chem A 2023; 127:8032-8049. [PMID: 37672217 DOI: 10.1021/acs.jpca.3c02927] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2023]
Abstract
Understanding the nature of the interaction between mercury(II) ions, Hg2+, and water molecules is crucial to describe the stability and chemical behavior of structures formed during solvation, as well as the conditions that favor the Hg2+ coordination or inducing water hydrolysis. In our study, we explored exhaustively the potential energy surface of Hg2+ with up to six water molecules. We analyzed electronic and Gibbs free energies, binding, and nuclear magnetic resonance parameters. We used the zeroth-order regular approximation Hamiltonian, including scalar and spin-orbit relativistic corrections for free energy calculations and geometry optimizations to explore the interplay between electron correlation and relativistic effects. We analyzed intermolecular interactions with energy decomposition analysis, quantum theory of atoms in molecules, and natural bond orbital. Additionally, we used the four-component Dirac Hamiltonian to compute solvent effect on the magnetic shielding and J-coupling constants. Our results revealed that the water hydrolysis by Hg2+ requires a minimum of three water molecules. We found that the interaction between Hg2+ and water molecules is an orbital interaction due to relativistic effects and the most stable structures are opened-shape clusters, reducing the number of oxygen-mercury contacts and maximizing the formation of hydrogen bonds among water molecules. In these types of clusters, Hg2+ promotes the water hydrolysis over coordination with oxygen atoms. However, when we considered the change associated with the transfer of a cluster from the ideal gas to a solvated system, our solvation free energy analysis revealed that closed-shape clusters are more favorable, maximizing the number of oxygen-mercury contacts and reducing the formation of hydrogen bonds among water molecules. This finding suggests that, under room conditions, the coordination of Hg2+ is more favorable than hydrolysis. Our results have significant implications for understanding Hg2+ behavior in water, helping to develop targeted strategies for mercury remediation and management, and contributing to advancements in the broader field of environmental chemistry.
Collapse
Affiliation(s)
- Edison Florez
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University, 0632 Auckland, New Zealand
| | - Andy D Zapata-Escobar
- Institute for Modeling and Innovative Technology (IMIT), CONICET-UNNE, W3404AAS Corrientes, Argentina
- Physics Department, Natural and Exact Science Faculty, Northeastern University, W3404AAS Corrientes, Argentina
| | - Franklin Ferraro
- Departamento de Ciencias Básicas, Universidad Católica Luis Amigó, 050034 Medellín, Colombia
| | - César Ibargüen Becerra
- Institute of Chemistry, University of Antioquia, 050010 Medellín, Colombia
- Facultad de Arquitectura e Ingeniería, Institución Universitaria Colegio Mayor de Antioquia (IUCMA), 050034 Medellín, Colombia
| | - Yuly Chamorro
- Van Swinderen Institute for Particle Physics and Gravity, University of Groningen, 9747 AG Groningen, The Netherlands
| | - Alejandro F Maldonado
- Institute for Modeling and Innovative Technology (IMIT), CONICET-UNNE, W3404AAS Corrientes, Argentina
| |
Collapse
|
3
|
Wang Z, Liu Y, Zhang W, Wang Y, Xu H, Yang L, Feng J, Hou B, Li M, Yan W. Selective mercury adsorption and enrichment enabled by phenylic carboxyl functionalized poly(pyrrole methane)s chelating polymers. THE SCIENCE OF THE TOTAL ENVIRONMENT 2023; 858:159870. [PMID: 36328257 DOI: 10.1016/j.scitotenv.2022.159870] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/01/2022] [Revised: 10/26/2022] [Accepted: 10/27/2022] [Indexed: 06/16/2023]
Abstract
Mercury decontamination from water requires highly effective and efficient methods for maintaining public health and environmental protection. Herein, based on the coordination theory between functional groups and metal ions, we proposed phenylic carboxyl group-based poly(pyrrole methane)s (PPDCBAs) as highly efficient mercury removal materials for environmental remediation applications. It was found that PPDCBAs can efficiently adsorb and remove mercury(II) from aqueous solutions by functionalizing the molecular structure with phenylic carboxyl groups. Among the as-prepared PPDCBAs, poly[pyrrole-2, 5-diyl (4-carboxybenzylidane)] (PPD4CBA) with the carboxyl group at the para position can not only adsorb mercury over 1400 mg⋅g-1 but also achieve a 92.5 % mercury(II) uptake within 100 min by a very low dosage of 0.1 g⋅L-1. In addition, PPDCBAs exhibited excellent adsorption selectivity for mercury(II) compared with copper(II), cadmium(II), zinc(II) and lead(II). Furthermore, as determined by Fourier transform infrared (FT-IR) spectra, X-ray photoelectron spectroscopy (XPS) and the density functional theory (DFT) calculation, the mercury removal was found to be mainly dependent on the high density of chelating sites, the phenylic carboxyl moieties, which helped us to realize an ultra-trace amount mercury removal (from 10.8 μg⋅L-1 to 0.6-0.8 μg⋅L-1) for meeting drinking water standard requirements (1.0 μg⋅L-1).
Collapse
Affiliation(s)
- Zhenyu Wang
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an 710049, China
| | - Yunpeng Liu
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an 710049, China
| | - Wenlong Zhang
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an 710049, China; School of Eco-Environment, Hebei University, Baoding 071002, China
| | - Yubing Wang
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an 710049, China
| | - Hao Xu
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an 710049, China
| | - Liu Yang
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an 710049, China
| | - Jiangtao Feng
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an 710049, China; Jiangsu Engineering Laboratory of New Materials for Sewage Treatment and Recycling, Soochow University, Suzhou, Jiangsu 215123, China.
| | - Bo Hou
- School of Physics and Astronomy, Cardiff University, The Parade, Cardiff CF24 3AA, UK.
| | - Mingtao Li
- International Research Center for Renewable Energy (IRCRE), State Key Laboratory of Multiphase Flow in Power Engineering (MFPE), Xi'an Jiaotong University, Xi'an, Shaanxi 710049, China
| | - Wei Yan
- Department of Environmental Engineering, Xi'an Key Laboratory of Solid Waste Recycling and Resource Recovery, School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an 710049, China
| |
Collapse
|
4
|
Taamalli S, Pitoňák M, Dibble TS, Černušák I, Louis F. Theoretical Study of the Monohydration of Mercury Compounds of Atmospheric Interest. J Phys Chem A 2021; 125:5819-5828. [PMID: 34180661 DOI: 10.1021/acs.jpca.1c02772] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The structures, vibrational frequencies, and model IR spectra of the monohydrates of oxygenated mercury compounds (BrHgO, BrHgOH, BrHgOOH, BrHgNO2, BrHgONO, and HgOH) have been theoretically studied using the ωB97X-D/aug-cc-pVTZ level of theory. The ground state potential energy surface exhibits several stable structures of these monohydrates. The thermodynamic properties of the hydration reactions have been calculated at different levels of theory including DFT and coupled-cluster calculations DK-CCSD(T) with the ANO-RCC-Large basis sets. Standard enthalpies and Gibbs free energies of hydration were computed. The temperature dependence of ΔrG°(T) was evaluated for the most stable complexes over the temperature range 200-400 K. Thermodynamic data revealed that the highest fraction hydrated at 298 K and 100% relative humidity will be BrHgNO2-H2O at ∼5%.
Collapse
Affiliation(s)
- Sonia Taamalli
- CNRS, UMR 8522 - PC2A - PhysicoChimie des Processus de Combustion et de l'Atmosphère, Université Lille, 59000 Lille, France
| | - Michal Pitoňák
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, 84215 Bratislava, Slovakia.,Computing Center, Centre of Operations of the Slovak Academy of Sciences, Dúbravská cesta 9, 845 35 Bratislava, Slovakia
| | - Theodore S Dibble
- Department of Chemistry, State University of New York-College of Environmental Science and Forestry, 1 Forestry Drive, Syracuse, New York 13210, United States
| | - Ivan Černušák
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, 84215 Bratislava, Slovakia
| | - Florent Louis
- CNRS, UMR 8522 - PC2A - PhysicoChimie des Processus de Combustion et de l'Atmosphère, Université Lille, 59000 Lille, France
| |
Collapse
|
5
|
Al Khalyfeh K, Afaneh AT, Marashdeh A, Almatarneh MH, Al-Mazaideh GM, Mizyed S, Ashram M. Thiacrown Ethers Engaged C 60 through Charge Transfer: Experimental and Theoretical Study. ACS OMEGA 2020; 5:25049-25058. [PMID: 33043183 PMCID: PMC7542602 DOI: 10.1021/acsomega.0c01877] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Accepted: 09/15/2020] [Indexed: 06/11/2023]
Abstract
UV-Vis spectroscopy is used to study the charge transfer complexes of thiacrown ethers 1-6 with fullerene. The size of TCE1-6 and the nature of the heteroatoms (N, O and S) have been systematically changed to examine the effect of these factors on the HOMO/LUMO energy levels, the optical energy gap and the interactions between TCE's and C60. The negative and positive values of ΔS designate the structural forming method and the randomness of the free solvent molecules, respectively. Thermodynamics and stability data show that the complexes have a 1:1 ratio that has been emphasized by density functional theory calculations. Additionally, they show a synergetic interplay of donor-acceptor, π-π, and n-π interactions, which are the basis for the affinity of our novel receptors toward C60. The proposed system of enzyme model suggests a development concept in the future design of enzyme model organic photovoltaic systems.
Collapse
Affiliation(s)
- Khaled Al Khalyfeh
- Department
of Chemistry, Al-Hussein Bin Talal University, Ma’an 71111, Jordan
| | - Akef T. Afaneh
- Department
of Chemistry, Faculty of Science, Al-Balqa
Applied University, P.O. Box 19117, Al-Salt 19117, Jordan
| | - Ali Marashdeh
- Department
of Chemistry, Faculty of Science, Al-Balqa
Applied University, P.O. Box 19117, Al-Salt 19117, Jordan
| | - Mansour H. Almatarneh
- Department
of Chemistry, University of Jordan, Amman 11942, Jordan
- Department
of Chemistry, Memorial University, St. John’s, NL A1B 3X7, Canada
| | - Ghassab M. Al-Mazaideh
- Department
of Chemistry and Chemical Technology, Faculty of Science, Tafila Technical University, P.O. Box 179, Tafila 66110, Jordan
- School of
Pharmaceutical Sciences, Universiti Sains
Malaysia, 11800 Penang , Malaysia
| | - Shehadeh Mizyed
- Department
of Chemistry, Yarmouk University, Irbid 21163, Jordan
| | - Muhammad Ashram
- Department
of Chemistry, Mutah University, Al Karak 61710, Jordan
| |
Collapse
|
6
|
Motta LC, Chien AD, Rask AE, Zimmerman PM. Mercury Magnetic Isotope Effect: A Plausible Photochemical Mechanism. J Phys Chem A 2020; 124:3711-3719. [PMID: 32314920 DOI: 10.1021/acs.jpca.0c00661] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Large mass-independent fractionation signatures in Hg have been observed in the laboratory and the environment, prompting deep questions about the chemical reasons behind these signatures. Since the relative lack of mechanistic information about Hg chemistry in the environment has precluded explanations of these isotope effects, the present study uses high-level electronic structure methods to evaluate the possible photochemical mechanisms of mass-independent isotope effects (MIEs) in HgX2 and CH3HgX (X = Cl, Br, I, and SCH3). The results show that spin-orbit coupling wipes out the potential of MIEs for Hg bound to Br or I, but that complexes involving lighter elements, HgX2 and CH3HgX (X = Cl and SCH3), have relatively small spin-orbit couplings upon photolysis. This unexpected finding shows that magnetic isotope fractionation due to hyperfine coupling is possible, depending on the identity of the Hg complex. By examination of the photolysis potential energy profiles, this study shows that HgX2 complexes can have a positive or a negative MIE (depending on reaction conditions), while CH3HgX complexes exclusively result in a positive MIE. These findings agree with MIE recorded in natural samples, demonstrating a plausible mechanism for the surprising mass-independent fractionation of Hg in the environment.
Collapse
Affiliation(s)
- Laura C Motta
- Department of Chemistry, University of Michigan, 930 N. University Avenue, Ann Arbor, Michigan 48109, United States
| | - Alan D Chien
- Department of Chemistry, University of Michigan, 930 N. University Avenue, Ann Arbor, Michigan 48109, United States
| | - Alan E Rask
- Department of Chemistry, University of Michigan, 930 N. University Avenue, Ann Arbor, Michigan 48109, United States
| | - Paul M Zimmerman
- Department of Chemistry, University of Michigan, 930 N. University Avenue, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
7
|
Afaneh AT, Schreckenbach G, Wang F. Theoretical Study of the Formation of Mercury (Hg2+) Complexes in Solution Using an Explicit Solvation Shell in Implicit Solvent Calculations. J Phys Chem B 2014; 118:11271-83. [DOI: 10.1021/jp5045089] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Akef T. Afaneh
- Department of Chemistry and ‡Centre for Earth Observation Science, Department
of Environment and Geography, University of Manitoba, Winnipeg, MB Canada, R3T 2N2
| | - Georg Schreckenbach
- Department of Chemistry and ‡Centre for Earth Observation Science, Department
of Environment and Geography, University of Manitoba, Winnipeg, MB Canada, R3T 2N2
| | - Feiyue Wang
- Department of Chemistry and ‡Centre for Earth Observation Science, Department
of Environment and Geography, University of Manitoba, Winnipeg, MB Canada, R3T 2N2
| |
Collapse
|
8
|
Sun R, Wang D, Mao W, Zhao S, Zhang C. Roles of chloride ion in photo-reduction/oxidation of mercury. CHINESE SCIENCE BULLETIN 2014. [DOI: 10.1007/s11434-014-0435-y] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
9
|
Šebera J, Burda J, Straka M, Ono A, Kojima C, Tanaka Y, Sychrovský V. Formation of a thymine-Hg(II)-thymine metal-mediated DNA base pair: proposal and theoretical calculation of the reaction pathway. Chemistry 2013; 19:9884-94. [PMID: 23766024 DOI: 10.1002/chem.201300460] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2013] [Revised: 04/24/2013] [Indexed: 11/07/2022]
Abstract
A reaction mechanism that describes the substitution of two imino protons in a thymine:thymine (T:T) mismatched DNA base pair with a Hg(II) ion, which results in the formation of a (T)N3-Hg(II)-N3(T) metal-mediated base pair was proposed and calculated. The mechanism assumes two key steps: The formation of the first Hg(II)-N3(T) bond is triggered by deprotonation of the imino N3 atom in thymine with a hydroxo ligand on the Hg(II) ion. The formation of the second Hg(II)-N3(T) bond proceeds through water-assisted tautomerization of the remaining, metal-nonbonded thymine base or through thymine deprotonation with a hydroxo ligand of the Hg(II) ion already coordinated to the thymine base. The thermodynamic parameters ΔGR =-9.5 kcal mol(-1), ΔHR =-4.7 kcal mol(-1), and ΔSR =16.0 cal mol(-1) K(-1) calculated with the ONIOM (B3LYP:BP86) method for the reaction agreed well with the isothermal titration calorimetric (ITC) measurements by Torigoe et al. [H. Torigoe, A. Ono, T. Kozasa, Chem. Eur. J. 2010, 16, 13218-13225]. The peculiar positive reaction entropy measured previously was due to both dehydration of the metal and the change in chemical bonding. The mercury reactant in the theoretical model contained one hydroxo ligand in accord with the experimental pKa value of 3.6 known for an aqua ligand of a Hg(II) center. The chemical modification of T:T mismatched to the T-Hg(II)-T metal-mediated base pair was modeled for the middle base pair within a trinucleotide B-DNA duplex, which ensured complete dehydration of the Hg(II) ion during the reaction.
Collapse
Affiliation(s)
- Jakub Šebera
- Department of Functional Materials, Institute of Physics, Academy of Sciences of the Czech Republic, v.v.i. Na Slovance 2, 182 21 Prague 8, Czech Republic
| | | | | | | | | | | | | |
Collapse
|
10
|
|