Predicting the acidity constant of a goethite hydroxyl group from first principles

Area of Science:

• Environmental Chemistry
• Surface Science
• Computational Chemistry

Background:

• Accurate prediction of hydroxyl group acid-base behavior on mineral surfaces is essential for understanding contaminant ion sorption in soils.
• Goethite, a common iron oxide mineral, plays a significant role in environmental processes.

Purpose of the Study:

• To calculate the pK(a) of a specific hydroxyl group on the goethite (101) surface using advanced computational methods.
• To compare ab initio molecular dynamics (AIMD) with empirical multi-site complexation (MUSIC) models for surface complexation studies.

Main Methods:

• Ab initio molecular dynamics (AIMD) simulations were employed to model the goethite surface.
• Potential-of-mean-force (PMF) techniques were utilized to determine the pK(a) values.
• Formic acid served as a reference system for calibration.

Main Results:

• A pK(a) of 7.0 was predicted for the doubly protonated oxygen atom bonded to a single Fe atom (Fe(I)OH(2)) on the goethite (101) surface.
• AIMD simulations provided insights into hydrogen bonding and proton sharing, which are critical for accurate surface complexation modeling.
• The study highlighted the importance of electronic structure methods capable of accurately predicting transition metal ion properties.

Conclusions:

• The calculated pK(a) suggests that positively charged hydroxyl groups on goethite are marginally stable at neutral pH.
• AIMD offers a more detailed treatment of surface interactions compared to empirical models.
• Accurate computational methods are vital for predicting mineral surface reactivity and its environmental implications.