1
|
Riefer A, Hackert-Oschätzchen M, Plänitz P, Meichsner G. Characterization of iron(III) in aqueous and alkaline environments with ab initio and ReaxFF potentials. J Chem Phys 2024; 160:082501. [PMID: 38411229 DOI: 10.1063/5.0182460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2023] [Accepted: 02/05/2024] [Indexed: 02/28/2024] Open
Abstract
The iron(III) complexes [Fe(H2O)n(OH)m]3-m (n + m = 5, 6, m ≤ 3) and corresponding proton transfer reactions are studied with total energy calculations, the nudged elastic band (NEB) method, and molecular dynamics (MD) simulations using ab initio and a modification of reactive force field potentials, the ReaxFF-AQ potentials, based on the implementation according to Böhm et al. [J. Phys. Chem. C 120, 10849-10856 (2016)]. Applying ab initio potentials, the energies for the reactions [Fe(H2O)n(OH)m]3-m + H2O → [Fe(H2O)n-1(OH)m+1]2-m + H3O+ in a gaseous environment are in good agreement with comparable theoretical results. In an aqueous (aq) or alkaline environment, with the aid of NEB computations, respective minimum energy paths with energy barriers of up to 14.6 kcal/mol and a collective transfer of protons are modeled. Within MD simulations at room temperature, a permanent transfer of protons around the iron(III) ion is observed. The information gained concerning the geometrical and energetic properties of water and the [Fe(H2O)n(OH)m]3-m complexes from the ab initio computations has been used as reference data to optimize parameters for the O-H-Fe interaction within the ReaxFF-AQ approach. For the optimized ReaxFF-AQ parameter set, the statistical properties of the basic water model, such as the radial distribution functions and the proton hopping functions, are evaluated. For the [Fe(H2O)n(OH)m]3-m complexes, it was found that while geometrical and energetic properties are in good agreement with the ab initio data for gaseous environment, the statistical properties as obtained from the MD simulations are only partly in accordance with the ab initio results for the iron(III) complexes in aqueous or alkaline environments.
Collapse
Affiliation(s)
- Arthur Riefer
- Chair of Manufacturing Technology with Focus Machining, Institute of Manufacturing Technology and Quality Management (IFQ), Faculty of Mechanical Engineering, Otto von Guericke University Magdeburg, 39106 Magdeburg, Germany
| | - Matthias Hackert-Oschätzchen
- Chair of Manufacturing Technology with Focus Machining, Institute of Manufacturing Technology and Quality Management (IFQ), Faculty of Mechanical Engineering, Otto von Guericke University Magdeburg, 39106 Magdeburg, Germany
| | - Philipp Plänitz
- AQcomputare Gesellschaft für Materialberechnung mbH, 09125 Chemnitz, Germany
| | - Gunnar Meichsner
- Chair of Manufacturing Technology with Focus Machining, Institute of Manufacturing Technology and Quality Management (IFQ), Faculty of Mechanical Engineering, Otto von Guericke University Magdeburg, 39106 Magdeburg, Germany
| |
Collapse
|
2
|
Shi H, Gong LD, Liu C, Lu LN, Yang ZZ. ABEEM/MM OH - Models for OH -(H 2O) n Clusters and Aqueous OH -: Structures, Charge Distributions, and Binding Energies. J Phys Chem A 2020; 124:5963-5978. [PMID: 32520555 DOI: 10.1021/acs.jpca.0c03941] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Based on the atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM), two fluctuating charge models of OH--water system were proposed. The difference between these two models is whether there is charge transfer between OH- and its first-shell water molecules. The structures, charge distributions, charge transfer, and binding energies of the OH-(H2O)n (n = 1-8, 10, 15, 23) clusters were studied by these two ABEEM/MM models, the OPLS/AA force field, the OPLS-SMOOTH/AA force field, and the QM methods. The results demonstrate that two ABEEM/MM models can search out all stable structures just as the QM methods, and the structures and charge distributions agree well with those from the QM calculations. The structures, the charge transfer, and the strength of hydrogen bonds in the first hydration shell are closely related to the coordination number of OH-. Molecular dynamics simulations on the aqueous OH- solution are performed at 298 and 278 K using ABEEM/MM-I model. The MD results show that the populations of three-, four-, and five-coordinated OH- are 29.6%, 67.1%, and 3.4% at 298 K, respectively, and those of two-, three-, four-, and five-coordinated OH- are 10.8%, 44.9%, 39.2%, and 4.9% at 278 K, respectively; the average hydrogen bond lengths and the hydrogen bond angle in the first shell increase with the temperature decreasing.
Collapse
Affiliation(s)
- Hua Shi
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.,School of Marine Science and Environment, Dalian Ocean University, Dalian 116023, People's Republic of China
| | - Li-Dong Gong
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| | - Cui Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| | - Li-Nan Lu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| |
Collapse
|
3
|
Weitkamp RF, Neumann B, Stammler HG, Hoge B. Generation and Applications of the Hydroxide Trihydrate Anion, [OH(OH 2 ) 3 ] - , Stabilized by a Weakly Coordinating Cation. Angew Chem Int Ed Engl 2019; 58:14633-14638. [PMID: 31373109 PMCID: PMC6790940 DOI: 10.1002/anie.201908589] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2019] [Indexed: 11/17/2022]
Abstract
The reaction of a strongly basic phosphazene (Schwesinger base) with water afforded the corresponding metastable hydroxide trihydrate [OH(OH2)3]− salt. This is the first hydroxide solvate that is not in contact with a cation and furthermore one of rare known water‐stabilized hydroxide anions. Thermolysis in vacuum results in the decomposition of the hydroxide salt and quantitative liberation of the free phosphazene base. This approach was used to synthesize the Schwesinger base from its hydrochloride salt after anion exchange in excellent yields of over 97 %. This deprotonation method can also be used for the phosphazene‐base‐catalyzed preparation of the Ruppert–Prakash reagent Me3SiCF3 using fluoroform (HCF3) as the trifluoromethyl building block and sodium hydroxide as the formal deprotonation agent.
Collapse
Affiliation(s)
- Robin F Weitkamp
- Centrum für Molekulare Materialien, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, 33615, Bielefeld, Germany
| | - Beate Neumann
- Centrum für Molekulare Materialien, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, 33615, Bielefeld, Germany
| | - Hans-Georg Stammler
- Centrum für Molekulare Materialien, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, 33615, Bielefeld, Germany
| | - Berthold Hoge
- Centrum für Molekulare Materialien, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, 33615, Bielefeld, Germany
| |
Collapse
|
4
|
Weitkamp RF, Neumann B, Stammler H, Hoge B. Generierung und Anwendung des Hydroxid‐Trihydrat‐Anions, [OH(OH
2
)
3
]
−
, stabilisiert durch ein schwach koordinierendes Kation. Angew Chem Int Ed Engl 2019. [DOI: 10.1002/ange.201908589] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Robin F. Weitkamp
- Centrum für Molekulare Materialien Fakultät für Chemie Universität Bielefeld Universitätsstraße 25 33615 Bielefeld Deutschland
| | - Beate Neumann
- Centrum für Molekulare Materialien Fakultät für Chemie Universität Bielefeld Universitätsstraße 25 33615 Bielefeld Deutschland
| | - Hans‐Georg Stammler
- Centrum für Molekulare Materialien Fakultät für Chemie Universität Bielefeld Universitätsstraße 25 33615 Bielefeld Deutschland
| | - Berthold Hoge
- Centrum für Molekulare Materialien Fakultät für Chemie Universität Bielefeld Universitätsstraße 25 33615 Bielefeld Deutschland
| |
Collapse
|
5
|
Ghosh S, Verma P, Cramer CJ, Gagliardi L, Truhlar DG. Combining Wave Function Methods with Density Functional Theory for Excited States. Chem Rev 2018; 118:7249-7292. [PMID: 30044618 DOI: 10.1021/acs.chemrev.8b00193] [Citation(s) in RCA: 128] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.
Collapse
Affiliation(s)
- Soumen Ghosh
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Christopher J Cramer
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| |
Collapse
|
6
|
Gadre SR, Yeole SD, Sahu N. Quantum chemical investigations on molecular clusters. Chem Rev 2014; 114:12132-73. [PMID: 25341561 DOI: 10.1021/cr4006632] [Citation(s) in RCA: 147] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Shridhar R Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur , Kanpur 208 016, India
| | | | | |
Collapse
|
7
|
Accuracy of density functionals in the description of dispersion interactions and IR spectra of phosphates and phosphorylated compounds. J Mol Model 2014; 20:2426. [DOI: 10.1007/s00894-014-2426-y] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2014] [Accepted: 08/07/2014] [Indexed: 10/24/2022]
|
8
|
Goerigk L, Grimme S. Double-hybrid density functionals. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2014. [DOI: 10.1002/wcms.1193] [Citation(s) in RCA: 209] [Impact Index Per Article: 20.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Lars Goerigk
- School of Chemistry; The University of Sydney; Sydney New South Wales Australia
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie; Universität Bonn; Bonn Germany
| |
Collapse
|
9
|
Yu HZ, Li C, Chen BH, Yang CT, Wang D, Fu Y, Hu S, Dang Z. Promising density functional theory methods for predicting the structures of uranyl complexes. RSC Adv 2014. [DOI: 10.1039/c4ra08264h] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
By examining the overall accuracy of different theoretical methods in predicting the U–X bond distances (of a series uranyl complexes), we found that both the global-hybrid meta-GGA functional of BB1K and the range-seperated LC-BLYP functional are fairly good (even better than the popular B3LYP method).
Collapse
Affiliation(s)
- Hai-Zhu Yu
- Department of Polymer Science and Engineering
- University of Science and Technology Beijing
- Beijing 100083, China
| | - Can Li
- Department of Polymer Science and Engineering
- University of Science and Technology Beijing
- Beijing 100083, China
| | - Bai-Hua Chen
- Institute of Nuclear Physics and Chemistry
- CAEP
- Mianyang, China
| | - Chu-Ting Yang
- Institute of Nuclear Physics and Chemistry
- CAEP
- Mianyang, China
| | - Dongrui Wang
- Department of Polymer Science and Engineering
- University of Science and Technology Beijing
- Beijing 100083, China
| | - Yao Fu
- Department of Chemistry
- University of Science and Technology of China
- Hefei 230026, China
| | - Sheng Hu
- Institute of Nuclear Physics and Chemistry
- CAEP
- Mianyang, China
| | - Zhimin Dang
- Department of Polymer Science and Engineering
- University of Science and Technology Beijing
- Beijing 100083, China
| |
Collapse
|
10
|
Li R, Peverati R, Isegawa M, Truhlar DG. Assessment and Validation of Density Functional Approximations for Iron Carbide and Iron Carbide Cation. J Phys Chem A 2012; 117:169-73. [DOI: 10.1021/jp3079106] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ruifang Li
- Department of Chemistry, Nankai University, Tianjin, P. R. China
300071
- Department
of Chemistry, Supercomputing Institute, and Chemical
Theory Center, University of Minnesota,
Minneapolis, Minnesota 55455,
United States
| | - Roberto Peverati
- Department
of Chemistry, Supercomputing Institute, and Chemical
Theory Center, University of Minnesota,
Minneapolis, Minnesota 55455,
United States
| | - Miho Isegawa
- Department
of Chemistry, Supercomputing Institute, and Chemical
Theory Center, University of Minnesota,
Minneapolis, Minnesota 55455,
United States
| | - Donald G. Truhlar
- Department
of Chemistry, Supercomputing Institute, and Chemical
Theory Center, University of Minnesota,
Minneapolis, Minnesota 55455,
United States
| |
Collapse
|
11
|
Theoretical study on reaction mechanism of sulfuric acid and ammonia and hydration of (NH4)2SO4. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1103-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
12
|
Zhang P, Fiedler L, Leverentz HR, Truhlar DG, Gao J. Polarized Molecular Orbital Model Chemistry. II. The PMO Method. J Chem Theory Comput 2011; 7:857-867. [PMID: 23378824 PMCID: PMC3560573 DOI: 10.1021/ct100638g] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
We present a new semiempirical molecular orbital method based on neglect of diatomic differential overlap. This method differs from previous NDDO-based methods in that we include p orbitals on hydrogen atoms to provide a more realistic modeling of polarizability. As in AM1-D and PM3-D, we also include damped dispersion. The formalism is based on the original MNDO one, but in the process of parameterization we make some specific changes to some of the functional forms. The present article is a demonstration of the capability of the new approach, and it presents a successful parametrization for compounds composed only of hydrogen and oxygen atoms, including the important case of water clusters.
Collapse
Affiliation(s)
- Peng Zhang
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, MN 55455-0431
| | - Luke Fiedler
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, MN 55455-0431
| | - Hannah R. Leverentz
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, MN 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, MN 55455-0431
| | - Jiali Gao
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, MN 55455-0431
| |
Collapse
|
13
|
Affiliation(s)
- Sławomir Janusz Grabowski
- Kimika Fakultatea, Euskal Herriko Unibertsitatea and Donostia International Physics Center (DIPC) P.K. 1072, 20080 Donostia, Euskadi, Spain
- IKERBASQUE, Basque Foundation for Science, 48011, Bilbao, Spain
| |
Collapse
|
14
|
Keil FJ. Multiscale Modelling in Computational Heterogeneous Catalysis. MULTISCALE MOLECULAR METHODS IN APPLIED CHEMISTRY 2011; 307:69-107. [DOI: 10.1007/128_2011_128] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
15
|
Šakić D, Zipse H, Vrček V. Base-catalyzed reactions of environmentally relevant N-chloro-piperidines. A quantum-chemical study. Org Biomol Chem 2011; 9:4336-46. [DOI: 10.1039/c1ob05077j] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
16
|
Neela YI, Mahadevi AS, Sastry GN. Hydrogen Bonding in Water Clusters and Their Ionized Counterparts. J Phys Chem B 2010; 114:17162-71. [DOI: 10.1021/jp108634z] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Y. Indra Neela
- Molecular Modeling Group, Organic Chemical Sciences, Indian Institute of Chemical Technology, Tarnaka, Hyderabad 500 607, AP, India
| | - A. Subha Mahadevi
- Molecular Modeling Group, Organic Chemical Sciences, Indian Institute of Chemical Technology, Tarnaka, Hyderabad 500 607, AP, India
| | - G. Narahari Sastry
- Molecular Modeling Group, Organic Chemical Sciences, Indian Institute of Chemical Technology, Tarnaka, Hyderabad 500 607, AP, India
| |
Collapse
|
17
|
Wang B, Truhlar DG. Including Charge Penetration Effects in Molecular Modeling. J Chem Theory Comput 2010; 6:3330-42. [DOI: 10.1021/ct1003862] [Citation(s) in RCA: 73] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Bo Wang
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431, United States
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
18
|
|
19
|
Bryantsev VS, Diallo MS, van Duin ACT, Goddard WA. Evaluation of B3LYP, X3LYP, and M06-Class Density Functionals for Predicting the Binding Energies of Neutral, Protonated, and Deprotonated Water Clusters. J Chem Theory Comput 2009; 5:1016-26. [PMID: 26609610 DOI: 10.1021/ct800549f] [Citation(s) in RCA: 272] [Impact Index Per Article: 18.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Vyacheslav S. Bryantsev
- Materials and Process Simulation Center, Beckman Institute, MC 139-74, California Institute of Technology, Pasadena, California 91125, and Department of Mechanical and Nuclear Engineering, The Pennsylvania State University, University Park, Pennsylvania 16801
| | - Mamadou S. Diallo
- Materials and Process Simulation Center, Beckman Institute, MC 139-74, California Institute of Technology, Pasadena, California 91125, and Department of Mechanical and Nuclear Engineering, The Pennsylvania State University, University Park, Pennsylvania 16801
| | - Adri C. T. van Duin
- Materials and Process Simulation Center, Beckman Institute, MC 139-74, California Institute of Technology, Pasadena, California 91125, and Department of Mechanical and Nuclear Engineering, The Pennsylvania State University, University Park, Pennsylvania 16801
| | - William A. Goddard
- Materials and Process Simulation Center, Beckman Institute, MC 139-74, California Institute of Technology, Pasadena, California 91125, and Department of Mechanical and Nuclear Engineering, The Pennsylvania State University, University Park, Pennsylvania 16801
| |
Collapse
|
20
|
Dinescu A, Benson MT. Electronic Structure Studies on Deprotonation of Dithiophosphinic Acids in Water Clusters. J Phys Chem A 2008; 112:12270-80. [DOI: 10.1021/jp806568k] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Adriana Dinescu
- Idaho National Laboratory, P.O. Box 1625, Idaho Falls, Idaho 83415-2208
| | - Michael T. Benson
- Idaho National Laboratory, P.O. Box 1625, Idaho Falls, Idaho 83415-2208
| |
Collapse
|
21
|
|