Hybrid Functional and Plane Waves based Ab Initio Molecular Dynamics Study of the Aqueous Fe2+ /Fe3+ Redox Reaction

Projected Wave Function Extrapolation Scheme to Accelerate Plane-Wave Hybrid Functional-Based Born-Oppenheimer Molecular Dynamics Simulations

Born-Oppenheimer molecular dynamics (BOMD) simulations are of great interest for the dynamic properties of molecular and solid systems. However, BOMD simulations necessitate not only an extensive period of dynamical evolution but also costly self-consistent-field (SCF) electronic structure calculations, especially for hybrid functional-based BOMD (H-BOMD) simulations within plane-wave basis sets. Here, we propose an improved always stable predictor-corrector (ASPC) method for the wave function extrapolation to accelerate the plane-wave H-BOMD simulations, named projected ASPC (PASPC), yielding a wave function closer to the actual solution space and efficiently reducing the number of SCF iterations at each MD step. We investigated the convergence properties of different extrapolation schemes for molecular and solid systems. Numerical results demonstrate that plane-wave H-BOMD simulations can be significantly faster than conventional cases by combining the accelerated algorithms with the PASPC method. The energy drift is also evaluated, showing that PASPC produces energy drift with smaller oscillations and can simulate a larger time step for systems containing heavy atoms, demonstrating the accuracy of the extrapolation schemes. Furthermore, H-BOMD simulations showcase more accurate power and infrared spectra of silicon dioxide and liquid water that are comparable to those of experimental measurements.

Evidence of Sharp Transitions between Octahedral and Capped Trigonal Prism States of the Solvation Shell of the Aqueous Fe3+ Ion

The structure of the solvation shell of the aqueous Fe3+ ion has been a subject of controversy due to discrepancies between experiments and different levels of theory. We address this issue by performing simulations for a wide range of ion concentrations, using various potential energy functions, supplemented by density functional theory calculations of selected configurations. The solvation shell undergoes abrupt transitions between two states: a hexacoordinated octahedral (OH) state and a capped trigonal prism (CTP) state with 7-fold coordination. The lifetime of these states is dependent on concentration. In dilute FeCl3 solutions, the lifetimes of both are similar (≈1 ns). However, the lifetime of the OH state increases with ion concentration, while that of the CTP state decreases slightly. When a uniform negative background charge is used instead of explicit counterions, the lifetime of the OH state is greatly overestimated. These findings underscore the need for further experimental measurements and higher-level simulations.

Characterization of iron(III) in aqueous and alkaline environments with ab initio and ReaxFF potentials

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.

Speeding-up Hybrid Functional-Based Ab Initio Molecular Dynamics Using Multiple Time-stepping and Resonance-Free Thermostat

Ab initio molecular dynamics (AIMD) based on density functional theory (DFT) has become a workhorse for studying the structure, dynamics, and reactions in condensed matter systems. Currently, AIMD simulations are primarily carried out at the level of generalized gradient approximation (GGA), which is at the second rung of DFT functionals in terms of accuracy. Hybrid DFT functionals, which form the fourth rung in the accuracy ladder, are not commonly used in AIMD simulations as the computational cost involved is 100 times or higher. To facilitate AIMD simulations with hybrid functionals, we propose here an approach using multiple time stepping with adaptively compressed exchange operator and resonance-free thermostat, that could speed up the calculations by ∼30 times or more for systems with a few hundred of atoms. We demonstrate that by achieving this significant speed up and making the compute time of hybrid functional-based AIMD simulations at par with that of GGA functionals, we are able to study several complex condensed matter systems and model chemical reactions in solution with hybrid functionals that were earlier unthinkable to be performed.