1
|
Purtscher FRS, Hofer TS. Probing the range of applicability of structure- and energy-adjusted QM/MM link bonds III: QM/MM MD simulations of solid-state systems at the example of layered carbon structures. J Comput Chem 2024; 45:2186-2197. [PMID: 38795379 DOI: 10.1002/jcc.27428] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/29/2024] [Revised: 04/11/2024] [Accepted: 04/29/2024] [Indexed: 05/27/2024]
Abstract
The previously introduced workflow to achieve an energetically and structurally optimized description of frontier bonds in quantum mechanical/molecular mechanics (QM/MM)-type applications was extended into the regime of computational material sciences at the example of a layered carbon model systems. Optimized QM/MM link bond parameters at HSEsol/6-311G(d,p) and self-consistent density functional tight binding (SCC-DFTB) were derived for graphitic systems, enabling detailed investigation of specific structure motifs occurring in graphene-derived structures v i a quantum-chemical calculations. Exemplary molecular dynamics (MD) simulations in the isochoric-isothermic (NVT) ensemble were carried out to study the intercalation of lithium and the properties of the Stone-Thrower-Wales defect. The diffusivity of lithium as well as hydrogen and proton adsorption on a defective graphene surface served as additional example. The results of the QM/MM MD simulations provide detailed insight into the applicability of the employed link-bond strategy when studying intercalation and adsorption properties of graphitic materials.
Collapse
Affiliation(s)
- Felix R S Purtscher
- Institute of General, Inorganic and Theoretical Chemistry Center for Chemistry and Biomedicine, University of Innsbruck, Innsbruck, Austria
| | - Thomas S Hofer
- Institute of General, Inorganic and Theoretical Chemistry Center for Chemistry and Biomedicine, University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
2
|
Purtscher FRS, Hofer TS. Towards hybrid quantum mechanical/molecular mechanical simulations of Li and Na intercalation in graphite - force field development and DFTB parametrisation. Phys Chem Chem Phys 2024; 26:1729-1740. [PMID: 38165417 DOI: 10.1039/d3cp02765a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2024]
Abstract
In this work a previously established QM/MM simulation protocol for the treatment of solid-state interfaces was extended towards the treatment of layered bulk materials enabling for instance investigation of metal intercalation in graphitic carbon materials. In order to study the intercalation of Li in graphite, new density functional tight binding (DFTB) parameters for Li have been created. Molecular dynamics (MD) simulations at constant temperatures (273.15, 298.15 and 323.15 K) have been carried out to assess the performance of the presented DFTB MD simulation approach. The intercalation of variable lithium and sodium content was investigated via z-distribution functions and analysis of the diffusivity in the direction parallel to the graphene plane. Both the calculated diffusion coefficients and the activation energy in case of lithium are in good agreement with experimental data. The comparison of the QM/MM MD simulation results provide detailed insights into the structural and dynamical properties of intercalated metal ions.
Collapse
Affiliation(s)
- Felix R S Purtscher
- Institute of General, Inorganic and Theoretical Chemistry Center for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria.
| | - Thomas S Hofer
- Institute of General, Inorganic and Theoretical Chemistry Center for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria.
| |
Collapse
|
3
|
Florian J, Agarwal H, Singh N, Goldsmith BR. Why halides enhance heterogeneous metal ion charge transfer reactions. Chem Sci 2021; 12:12704-12710. [PMID: 34703556 PMCID: PMC8494035 DOI: 10.1039/d1sc03642d] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2021] [Accepted: 08/26/2021] [Indexed: 11/21/2022] Open
Abstract
The reaction kinetics of many metal redox couples on electrode surfaces are enhanced in the presence of halides (i.e., Cl-, Br-, I-). Using first-principles metadynamics simulations, we show a correlation between calculated desorption barriers of V3+-anion complexes bound to graphite via an inner-sphere anion bridge and experimental V2+/V3+ kinetic measurements on edge plane pyrolytic graphite in H2SO4, HCl, and HI. We extend this analysis to V2+/V3+, Cr2+/Cr3+, and Cd0/Cd2+ reactions on a mercury electrode and demonstrate that reported kinetics in acidic electrolytes for these redox couples also correlate with the predicted desorption barriers of metal-anion complexes. Therefore, the desorption barrier of the metal-anion surface intermediate is a descriptor of kinetics for many metal redox couple/electrode combinations in the presence of halides. Knowledge of the metal-anion surface intermediates can guide the design of electrolytes and electrocatalysts with faster kinetics for redox reactions of relevance to energy and environmental applications.
Collapse
Affiliation(s)
- Jacob Florian
- Department of Chemical Engineering, University of Michigan Ann Arbor Michigan 48109-2136 USA
- Catalysis Science and Technology Institute, University of Michigan Ann Arbor Michigan 48109-2136 USA
| | - Harsh Agarwal
- Department of Chemical Engineering, University of Michigan Ann Arbor Michigan 48109-2136 USA
- Catalysis Science and Technology Institute, University of Michigan Ann Arbor Michigan 48109-2136 USA
| | - Nirala Singh
- Department of Chemical Engineering, University of Michigan Ann Arbor Michigan 48109-2136 USA
- Catalysis Science and Technology Institute, University of Michigan Ann Arbor Michigan 48109-2136 USA
| | - Bryan R Goldsmith
- Department of Chemical Engineering, University of Michigan Ann Arbor Michigan 48109-2136 USA
- Catalysis Science and Technology Institute, University of Michigan Ann Arbor Michigan 48109-2136 USA
| |
Collapse
|
4
|
Does uranyl-TBP complex formation happen at the aqueous-organic interface? Revelation by molecular dynamics simulations. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.115621] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
5
|
Chopra M, Choudhury N. Adsorption of uranyl ions from its aqueous solution by functionalized carbon nanotubes: A molecular dynamics simulation study. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.111569] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
6
|
Das A, Ali SM. Molecular Dynamics Simulation Studies on Structure, Dynamics, and Thermodynamics of Uranyl Nitrate Solution at Various Acid Concentrations. J Phys Chem B 2019; 123:4571-4586. [PMID: 31070371 DOI: 10.1021/acs.jpcb.9b01498] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The structural and dynamical characteristics of uranyl ions in an aqueous acidic environment are of immense importance in the field of nuclear fuel reprocessing. In view of that, the structural and dynamical behavior of the uranyl ion in water has been investigated by performing molecular dynamics (MD) simulations using different force fields. All the force fields have depicted similar structural and dynamical properties except the free energy of hydration where the Guilbaud-Wipff (GW) model performs well over the others. The calculated density using MD simulations is found to be in excellent agreement with the measured experimental density, which ensures the accuracy of the adopted GW force field. The calculated surface tension and shear viscosity are seen to be increased with uranyl nitrate concentrations. At a higher concentration of about 4.0 mol/L, the supersaturation effect has been captured by an inflection in the plot of surface tension and shear viscosity against concentration because of the solution heterogeneity, which was correlated by an inflection in the scattering intensity observed by performing the dynamic light scattering experiment. The binding mode of nitrate ions with the uranyl ion is found to be concentration-dependent, and at higher concentration, it is predominantly monodentate.
Collapse
Affiliation(s)
- Arya Das
- Nuclear Recycle Board , Bhabha Atomic Research Centre , Mumbai 400094 , India.,Homi Bhabha National Institute , Mumbai 400094 , India
| | - Sk Musharaf Ali
- Chemical Engineering Division , Bhabha Atomic Research Centre , Mumbai 400085 , India.,Homi Bhabha National Institute , Mumbai 400094 , India
| |
Collapse
|
7
|
Arumugam K, Burton NA. Density functional theory (DFT) calculations of VI/V reduction potentials of uranyl coordination complexes in non-aqueous solutions. Phys Chem Chem Phys 2019; 21:3227-3241. [PMID: 30681090 DOI: 10.1039/c8cp05412f] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Of particular interest within the +6 uranium complexes is the linear uranyl(vi) cation and it forms numerous coordination complexes in solution and exhibits incongruent redox behavior depending on coordinating ligands. In this study, to determine the reduction potentials of uranyl complexes in non-aqueous solutions, a hybrid density functional theory (DFT) approach was used in which two different DFT functionals, B3LYP and M06, were applied. Bulk solvent effects were invoked through the conductor-like polarizable continuum model. The solute cavities were described with the united-atom Kohn-Sham (UAKS) cavity definition. Inside the cavity the dielectric constant matches the value of a vacuum and outside the cavity the dielectric constant value is the same as that of the solvent of interest, for example, dimethyl sulfoxide (DMSO), dimethylformamide (DMF), dichloromethane (DCM), acetonitrile and pyridine. With the help of the Nernst equation the calculated reduction potentials with respect to the ferrocene (Fc) reference electrode are converted into reduction free energies (RFEs). Uranyl complexes of organic ligands which range from mono- to hexa-dentate coordination modes were investigated in non-aqueous solutions of DMSO, DMF, DCM, acetonitrile and pyridine solutions. The effect of the spin-orbit correction and the reference electrode correction on the RFEs and various methods such as the direct method and the isodesmic reaction model were explored. Overall, our computational determination of RFEs of uranyl complexes in various non-aqueous solutions demonstrates that the RFEs can be obtained within ∼0.2 eV of experimental values.
Collapse
Affiliation(s)
- Krishnamoorthy Arumugam
- School of Chemistry, The University of Manchester, Brunswick Street, Manchester M13 9PL, UK.
| | | |
Collapse
|
8
|
Liu X, Cheng J, He M, Lu X, Wang R. Acidity constants and redox potentials of uranyl ions in hydrothermal solutions. Phys Chem Chem Phys 2016; 18:26040-26048. [DOI: 10.1039/c6cp03469a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report a first principles molecular dynamics (FPMD) study of the structures, acidity constants (pKa) and redox potentials (E0) of uranyl (UO22+) from ambient conditions to 573 K.
Collapse
Affiliation(s)
- Xiandong Liu
- State Key Laboratory for Mineral Deposits Research
- School of Earth Sciences and Engineering
- Nanjing University
- Nanjing 210093
- P. R. China
| | - Jun Cheng
- College of Chemistry and Chemical Engineering
- Xiamen University
- Xiamen 361005
- P. R. China
- Department of Chemistry
| | - Mengjia He
- State Key Laboratory for Mineral Deposits Research
- School of Earth Sciences and Engineering
- Nanjing University
- Nanjing 210093
- P. R. China
| | - Xiancai Lu
- State Key Laboratory for Mineral Deposits Research
- School of Earth Sciences and Engineering
- Nanjing University
- Nanjing 210093
- P. R. China
| | - Rucheng Wang
- State Key Laboratory for Mineral Deposits Research
- School of Earth Sciences and Engineering
- Nanjing University
- Nanjing 210093
- P. R. China
| |
Collapse
|
9
|
Chopra M, Choudhury N. Molecular dynamics simulation study of distribution and dynamics of aqueous solutions of uranyl ions: the effect of varying temperature and concentration. Phys Chem Chem Phys 2015; 17:27840-50. [PMID: 26439497 DOI: 10.1039/c5cp03769g] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Investigating the characteristics of actinyl ions has been of great interest due to their direct relevance in the nuclear fuel cycle. All-atom molecular dynamics simulations have been employed to study the orientational structure and dynamics of aqueous solutions of uranyl ions of various concentrations. The orientational structure of water around a uranyl ion has been thoroughly investigated by calculating different orientational probability distributions corresponding to different molecular axes of water. The orientational distribution of water molecules in the first coordination shell of a uranyl ion is found to be markedly different from that in bulk water. Analysis of counterion distribution around the uranyl ion reveals the presence of nitrate ions along with water molecules in the first solvation shell. From the comparison of the number of coordinated water and nitrate ions at various uranyl nitrate concentrations, it is evident that these two species compete for occupying the first solvation shell of the uranyl ion. Orientational dynamics of water molecules about different molecular axes of water in the vicinity of uranyl ions have also been investigated and decreasing orientational mobility of water with increasing uranyl concentration has been found. However, it is observed that the orientational dynamics remains more or less the same whether we consider all the water molecules in the aqueous solution or only the solvation shell water molecules. The effect of temperature on the translational and orientational characteristics of the aqueous uranyl solutions has also been studied in detail.
Collapse
Affiliation(s)
- Manish Chopra
- Radiation Safety Systems Division, Bhabha Atomic Research Centre, Mumbai - 400 085, India
| | | |
Collapse
|
10
|
Chopra M, Choudhury N. Effect of Uranyl Ion Concentration on Structure and Dynamics of Aqueous Uranyl Solution: A Molecular Dynamics Simulation Study. J Phys Chem B 2014; 118:14373-81. [DOI: 10.1021/jp506477s] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Manish Chopra
- Radiation
Safety Systems Division and ‡Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai, 400 085, India
| | - Niharendu Choudhury
- Radiation
Safety Systems Division and ‡Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai, 400 085, India
| |
Collapse
|
11
|
Density functional theory and molecular dynamics study of the uranyl ion (UO2)2+. J Mol Model 2014; 20:2150. [DOI: 10.1007/s00894-014-2150-7] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2013] [Accepted: 01/21/2014] [Indexed: 10/25/2022]
|
12
|
Hofer TS. Perspectives for hybrid ab initio/molecular mechanical simulations of solutions: from complex chemistry to proton-transfer reactions and interfaces. PURE APPL CHEM 2014. [DOI: 10.1515/pac-2014-5019] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Abstract
As a consequence of the ongoing development of enhanced computational resources, theoretical chemistry has become an increasingly valuable field for the investigation of a variety of chemical systems. Simulations employing a hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) technique have been shown to be a particularly promising approach, whenever ultrafast (i.e., picosecond) dynamical properties are to be studied, which are in many cases difficult to access via experimental techniques. Details of the quantum mechanical charge field (QMCF) ansatz, an advanced QM/MM protocol, are discussed and simulation results for various systems ranging from simple ionic hydrates to solvated organic molecules and coordination complexes in solution are presented. A particularly challenging application is the description of proton-transfer reactions in chemical simulations, which is a prerequisite to study acidified and basic systems. The methodical requirements for a combination of the QMCF methodology with a dissociative potential model for the description of the solvent are discussed. Furthermore, the possible extension of QM/MM approaches to solid/liquid interfaces is outlined.
Collapse
|
13
|
Réal F, Trumm M, Schimmelpfennig B, Masella M, Vallet V. Further insights in the ability of classical nonadditive potentials to model actinide ion-water interactions. J Comput Chem 2012; 34:707-19. [DOI: 10.1002/jcc.23184] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2012] [Revised: 09/12/2012] [Accepted: 10/19/2012] [Indexed: 11/06/2022]
|
14
|
Atta-Fynn R, Johnson DF, Bylaska EJ, Ilton ES, Schenter GK, de Jong WA. Structure and Hydrolysis of the U(IV), U(V), and U(VI) Aqua Ions from Ab Initio Molecular Simulations. Inorg Chem 2012; 51:3016-24. [DOI: 10.1021/ic202338z] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Raymond Atta-Fynn
- Environmental
Molecular Sciences Laboratory, and ‡Chemical and Material Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Donald F. Johnson
- Environmental
Molecular Sciences Laboratory, and ‡Chemical and Material Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Eric J. Bylaska
- Environmental
Molecular Sciences Laboratory, and ‡Chemical and Material Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Eugene S. Ilton
- Environmental
Molecular Sciences Laboratory, and ‡Chemical and Material Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Gregory K. Schenter
- Environmental
Molecular Sciences Laboratory, and ‡Chemical and Material Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Wibe A. de Jong
- Environmental
Molecular Sciences Laboratory, and ‡Chemical and Material Sciences Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| |
Collapse
|
15
|
Wang D, van Gunsteren WF, Chai Z. Recent advances in computational actinoid chemistry. Chem Soc Rev 2012; 41:5836-65. [DOI: 10.1039/c2cs15354h] [Citation(s) in RCA: 108] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
|
16
|
Bühl M, Wipff G. Insights into Uranyl Chemistry from Molecular Dynamics Simulations. Chemphyschem 2011; 12:3095-105. [DOI: 10.1002/cphc.201100458] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2011] [Revised: 07/28/2011] [Indexed: 11/10/2022]
|