Molecular dynamics simulations of the interactions between water and inorganic solids

Kerisit, Sebastien N.; Cooke, David J.; Spagnoli, Dino; Parker, Stephen C.

doi:10.1039/b415633c

Cited by 99 publications

(118 citation statements)

References 42 publications

(68 reference statements)

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The green triangles show the growth of the number of acceptor bonds from a water molecule at x to a bulk water molecule. In the bulk region the total number of acceptor bonds is about 1.50, which is slightly less than the bulk water value, but is consistent with the acceptor bonding found for the bulk region in our previous Al 3+ + 64H 2 …”

Section: Hydrogen Bond Structure In the Hydration Shellssupporting

confidence: 74%

“…Highly charged ions in solution play key roles in mineral formation and extraction, 1 the stability of surfaces and nanoparticles, 2 as partially solvated ions ͑cofactors͒ activating biochemical reactions, 3,4 and in the transport of toxic materials. 5 Much of the chemistry of these ions in solution has been rationalized using phenomenological concepts such as Pearson hardness, 6 ion charge to size ratios, 7 electrophilicity, 8 nucleophilicity, 8 nucleofugality, 8 and other empirical structure-activity relationships.…”

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

Structure and dynamics of the hydration shells of the Zn2+ ion from ab initio molecular dynamics and combined ab initio and classical molecular dynamics simulations

Cauët

Bogatko

Weare

et al. 2010

The Journal of Chemical Physics

102

View full text Add to dashboard Cite

Results of ab initio molecular dynamics ͑AIMD͒ simulations ͑density functional theory+ PBE96͒ of the dynamics of waters in the hydration shells surrounding the Zn 2+ ion ͑T Ϸ 300 K, Ϸ 1 gm/ cm 3 ͒ are compared to simulations using a combined quantum and classical molecular dynamics ͓AIMD/molecular mechanical ͑MM͔͒ approach. Both classes of simulations were performed with 64 solvating water molecules ͑ϳ15 ps͒ and used the same methods in the electronic structure calculation ͑plane-wave basis set, time steps, effective mass, etc.͒. In the AIMD/MM calculation, only six waters of hydration were included in the quantum mechanical ͑QM͒ region. The remaining 58 waters were treated with a published flexible water-water interaction potential. No reparametrization of the water-water potential was attempted. Additional AIMD/MM simulations were performed with 256 water molecules. The hydration structures predicted from the AIMD and AIMD/MM simulations are found to agree in detail with each other and with the structural results from x-ray data despite the very limited QM region in the AIMD/MM simulation. To further evaluate the agreement of these parameter-free simulations, predicted extended x-ray absorption fine structure ͑EXAFS͒ spectra were compared directly to the recently obtained EXAFS data and they agree in remarkable detail with the experimental observations. The first hydration shell contains six water molecules in a highly symmetric octahedral structure is ͑maximally located at 2.13-2.15 Å versus 2.072 Å EXAFS experiment͒. The widths of the peak of the simulated EXAFS spectra agree well with the data ͑8.4 Å 2 versus 8.9 Å 2 in experiment͒. Analysis of the H-bond structure of the hydration region shows that the second hydration shell is trigonally bound to the first shell water with a high degree of agreement between the AIMD and AIMD/MM calculations. Beyond the second shell, the bonding pattern returns to the tetrahedral structure of bulk water. The AIMD/MM results emphasize the importance of a quantum description of the first hydration shell to correctly describe the hydration region. In these calculations the full d 10 electronic structure of the valence shell of the Zn 2+ ion is retained. The simulations show substantial and complex charge relocation on both the Zn 2+ ion and the first hydration shell. The dipole moment of the waters in the first hydration shell is 3.4 D ͑3.3 D AIMD/MM͒ versus 2.73 D bulk. Little polarization is found for the waters in the second hydration shell ͑2.8 D͒. No exchanges were seen between the first and the second hydrations shells; however, many water transfers between the second hydration shell and the bulk were observed. For 64 waters, the AIMD and AIMD/MM simulations give nearly identical results for exchange dynamics. However, in the larger particle simulations ͑256 waters͒ there is a significant reduction in the second shell to bulk exchanges.

show abstract

Section: Hydrogen Bond Structure In the Hydration Shellssupporting

confidence: 74%

Section: Introductionmentioning

confidence: 99%

Structure and dynamics of the hydration shells of the Zn2+ ion from ab initio molecular dynamics and combined ab initio and classical molecular dynamics simulations

Cauët

Bogatko

Weare

et al. 2010

The Journal of Chemical Physics

102

View full text Add to dashboard Cite

show abstract

“…For example, many calcite systems only exist at high temperatures and pressures in the earth's crust [10]. Properties of the calcite-water interface, such as the arrangement of water molecules on the surface, as well as surface and adsorption energies have been investigated using atomistic simulation techniques [6,9,33,34,7]. While the majority of the work done on the calcite-water interface has used classical methods, some work using electronic structure methods has been performed.…”

Section: The Calcite-water Interfacementioning

confidence: 99%

“…A large amount of theoretical work has been performed on the calcite-water interface [6,7,8]. This modelling has largely been through classical molecular dynamics (MD) simulations, modelling water in its associated form on the calcite {1014} surface, the most stable surface [8,9].…”

Section: Modelsmentioning

confidence: 99%

Investigation of the Interaction of Water with the Calcite (10.4) Surface Using Ab Initio Simulation

2009

View full text Add to dashboard Cite

Density functional theory calculations were employed to explore the interaction between water and the {1014} surface of calcite. In addition a defective {1014} surface and stepped surfaces in contact with water were investigated. A series of percentage water coverages and water configurations were explored, including dissociated water states. Static relaxations found associated water to be favourable on the {1014} surface, although a metastable dissociated state 1.77eV higher in energy was found.Molecular dynamics (MD) simulations of low water coverage reveal fluctuations in the H-O water bond when the H atom is directed towards a surface CO 3 ion. Desorption of an H 2 O molecule was observed in simulations above 900K. Water was found to be strongly bound to the perfect {1014} surface, with an adsorption energy of -0.91eV. MD simulations of a defective {1014} surface found water to favour dissociation at CO 3 vacancies. However, water at Ca vacancies diffused across the surface to form a bond with the nearest surface Ca ion. Water was also found to favour an associated state at both acute and obtuse steps. On all these imperfect surfaces water was found to adsorb strongly to the surface, with adsorption energies ranging from -0.99eV to -1.60eV.

show abstract

“…Key chemical processes such as adsorption, electron transfer, growth, and dissolution all depend principally on the atomic structure adopted at these interfaces. For example, dissolution and solute adsorption are regulated by the structure of interfacial water 7,8 . Surface acid/base chemistry 9,10 arises from the types and arrangement of terminal metal-coordinating aquo/hydroxyl groups 3,[11][12][13] .…”

mentioning

confidence: 99%

Dynamic Stabilization of Metal Oxide–Water Interfaces

et al. 2017

View full text Add to dashboard Cite

ABSTRACT:The interaction of water with metal oxide surfaces plays a crucial role in the catalytic and geochemical behavior of metal oxides. In a vast majority of studies, the interfacial structure is assumed to arise from a relatively static lowest energy configuration of atoms, even at room temperature. Using hematite (-Fe2O3) as a model oxide, we show through a direct comparison of in situ synchrotron X-ray scattering with density functional theorybased molecular dynamics (DFT-MD) simulations that the structure of the (11 ̅ 02) termination is dynamically stabilized by picosecond water exchange. Simulations show frequent exchanges between terminal aquo groups and adsorbed water in locations and with partial residence times consistent with experimentally determined atomic sites and fractional occupancies. Frequent water exchange occurs even for an ultrathin adsorbed water film persisting on the surface under a dry atmosphere. The resulting time-averaged interfacial structure consists of a ridged lateral arrangement of adsorbed water molecules hydrogen bonded to terminal aquo groups. Surface pKa prediction based on bond valence analysis suggests that water exchange will influence the proton transfer reactions underlying the acid/base reactivity at the interface. Our findings provide important new insights for understanding complex interfacial chemical processes at metal oxide-water interfaces.The interfaces between metal oxides and water are among the most important in nature and in emerging energy applications, with wide ranging impacts from photocatalytic water splitting [1][2][3][4] to the geochemical cycling of elements 5,6 . Key chemical processes such as adsorption, electron transfer, growth, and dissolution all depend principally on the atomic structure adopted at these interfaces. For example, dissolution and solute adsorption are regulated by the structure of interfacial water 7,8 . Surface acid/base chemistry 9,10 arises from the types and arrangement of terminal metal-coordinating aquo/hydroxyl groups 3,11-13 .At room temperature an interface is at dynamic equilibrium. In principle, the average interfacial structure depends on the interplay of relatively static atoms at the solid surface with relatively dynamic overlying water molecules. Simulations suggest that movement of overlying water molecules can play an essential role in stabilizing the interface and influencing its chemical behavior 14,15 . However, simulated [14][15][16][17] or spectroscopically probed 18 dynamics are seldom integrated with experimentally derived interface structure models to achieve comprehensive insight into interfacial structure 4 . To understand and predict chemical processes at dynamically active metal oxide-water interfaces, structure and dynamics must be considered as a unified whole.Accurate measurements of interface structure and water ordering rely on interface-sensitive synchrotron X-ray scattering methods 19 . The analysis of multiple crystal truncation rods (CTRs) provides a complete 3-dimensional interface mode...

show abstract

Molecular dynamics simulations of the interactions between water and inorganic solids

Cited by 99 publications

References 42 publications

Structure and dynamics of the hydration shells of the Zn2+ ion from ab initio molecular dynamics and combined ab initio and classical molecular dynamics simulations

Structure and dynamics of the hydration shells of the Zn2+ ion from ab initio molecular dynamics and combined ab initio and classical molecular dynamics simulations

Investigation of the Interaction of Water with the Calcite (10.4) Surface Using Ab Initio Simulation

Dynamic Stabilization of Metal Oxide–Water Interfaces

Contact Info

Product

Resources

About