Unlocking Phase Diagrams for Molybdenum and Tungsten Nanoclusters and Prediction of their Formation Constants

POMSimulator: An open-source tool for predicting the aqueous speciation and self-assembly mechanisms of polyoxometalates

Elucidating the speciation (in terms of concentration versus pH) and understanding the formation mechanisms of polyoxometalates remains a significant challenge, both in experimental and computational domains. POMSimulator is a new methodology that tackles this problem from a purely computational perspective. The methodology uses results from quantum mechanics based methods to automatically set up the chemical reaction network, and to build speciation models. As a result, it becomes possible to predict speciation and phase diagrams, as well as to derive new insights into the formation mechanisms of large molecular clusters. In this work we present the main features of the first open-source version of the software. Since the first report [Chem. Sci. 2020, 11, 8448-8456], POMSimulator has undergone several improvements to keep up with the growing challenges that were tackled. After four years of research, we recognize that the source code is sufficiently stable to share a polished and user-friendly version. The Python code, manual, examples, and install instructions can be found at https://github.com/petrusen/pomsimulator.

Computational Study into the Effects of Countercations on the [P8W48O184]40- Polyoxometalate Wheel

Porous metal oxide materials have been obtained from a ring-shaped macrocyclic polyoxometalate (POM) structural building unit, [P8W48O184]40-. This is a tungsten oxide building block with an integrated "pore" of 1 nm in diameter, which, when connected with transition metal linkers, can assemble frameworks across a range of dimensions and which are generally referred to as POMzites. Our investigation proposes to gain a better understanding into the basic chemistry of this POM, specifically local electron densities and locations of countercations within and without the aforementioned pore. Through a rigorous benchmarking process, we discovered that 8 potassium cations, located within the pore, provided us with the most accurate model in terms of mimicking empirical properties to a sufficient degree of accuracy while also requiring a relatively small number of computer cores and hours to successfully complete a calculation. Additionally, we analyzed two other similar POMs from the literature, [As8W48O184]40- and [Se8W48O176]32-, in the hopes of determining whether they could be similarly incorporated into a POMzite network; given their close semblance in terms of local electron densities and interaction with potassium cations, we judge these POMs to be theoretically suitable as POMzite building blocks. Finally, we experimented with substituting different cations into the [P8W48O184]40- pore to observe the effect on pore dimensions and overall reactivity; we observed that the monocationic structures, particularly the Li8[P8W48O184]32- framework, yielded the least polarized structures. This correlates with the literature, validating our methodology for determining general POM characteristics and properties moving forward.

Multi-Time-Scale Simulation of Complex Reactive Mixtures: How Do Polyoxometalates Form?

Understanding the dynamics of reactive mixtures still challenges both experiments and theory. A relevant example can be found in the chemistry of molecular metal-oxide nanoclusters, also known as polyoxometalates. The high number of species potentially involved, the interconnectivity of the reaction network, and the precise control of the pH and concentrations needed in the synthesis of such species make the theoretical/computational treatment of such processes cumbersome. This work addresses this issue relying on a unique combination of recently developed computational methods that tackle the construction, kinetic simulation, and analysis of complex chemical reaction networks. By using the Bell-Evans-Polanyi approximation for estimating activation energies, and an accurate and robust linear scaling for correcting the computed pKa values, we report herein multi-time-scale kinetic simulations for the self-assembly processes of polyoxotungstates that comprise 22 orders of magnitude, from tens of femtoseconds to months of reaction time. This very large time span was required to reproduce very fast processes such as the acid/base equilibria (at 10-12 s), relatively slow reactions such as the formation of key clusters such as the metatungstate (at 103 s), and the very slow assembly of the decatungstate (at 106 s). Analysis of the kinetic data and of the reaction network topology shed light onto the details of the main reaction mechanisms, which explains the origin of kinetic and thermodynamic control followed by the reaction. Simulations at alkaline pH fully reproduce experimental evidence since clusters do not form under those conditions.

Speciation atlas of polyoxometalates in aqueous solutions

Speciation is the key parameter in solution chemistry that describes the composition, concentration, and oxidation state of each chemical form of an element present in a sample. The speciation study of complex polyatomic ions has remained challenging because of the large number of factors affecting stability and the limited number of direct methods. To address these challenges, we developed the speciation atlas of 10 polyoxometalates commonly used in catalytic and biological applications in aqueous solutions, where the speciation atlas provides both a species distribution database and a predictive model for other polyoxometalates to be used. Compiled for six different polyoxometalate archetypes with three types of addenda ions based on 1309 nuclear magnetic resonance spectra under 54 different conditions, the atlas has revealed a previously unknown behavior of polyoxometalates that may account for their potency as biological agents and catalysts. The atlas is intended to promote the interdisciplinary use of metal oxides in various scientific fields.

Computational Prediction of Speciation Diagrams and Nucleation Mechanisms: Molecular Vanadium, Niobium, and Tantalum Oxide Nanoclusters in Solution

Understanding the aqueous speciation of molecular metal-oxo-clusters plays a key role in different fields such as catalysis, electrochemistry, nuclear waste recycling, and biochemistry. To describe the speciation accurately, it is essential to elucidate the underlying self-assembly processes. Herein, we apply a computational method to predict the speciation and formation mechanisms of polyoxovanadates, -niobates, and -tantalates. While polyoxovanadates have been widely studied, polyoxoniobates and -tantalates lack the same level of understanding. First, we propose a pentavanadate cluster ([V₅O₁₄]^3-) as a key intermediate for the formation of the decavanadate. Our computed phase speciation diagram is in particularly good agreement with the experiments. Second, we report the formation constants of the heptaniobate, [Nb₇O₂₂]^9-, decaniobate, [Nb₁₀O₂₈]^6-, and tetracosaniobate [H₉Nb₂₄O₇₂]^15-. Additionally, we compute the speciation and phase diagram of niobium, which so far was restricted to Lindqvist derivates. Finally, we predict the formation constant of the decatantalate ([Ta₁₀O₂₆]^6-) in water, even though it had only been synthesized in toluene. Furthermore, we also calculate the corresponding speciation and phase diagrams for polyoxotantalates. Overall, we show that our method can be successfully applied to different families of molecular metal oxides without any need for readjustments; therefore, it can be regarded as a trustworthy tool for exploring polyoxometalates' chemistry.

Composition-driven archetype dynamics in polyoxovanadates

Molecular metal oxides often adopt common structural frameworks (i.e. archetypes), many of them boasting impressive structural robustness and stability. However, the ability to adapt and to undergo transformations between different structural archetypes is a desirable material design feature offering applicability in different environments. Using systems thinking approach that integrates synthetic, analytical and computational techniques, we explore the transformations governing the chemistry of polyoxovanadates (POVs) constructed of arsenate and vanadate building units. The water-soluble salt of the low nuclearity polyanion [V₆As₈O₂₆]⁴⁻ can be effectively used for the synthesis of the larger spherical (i.e. kegginoidal) mixed-valent [V₁₂As₈O₄₀]⁴⁻ precipitate, while the novel [V₁₀As₁₂O₄₀]⁸⁻ POVs having tubular cyclic structures are another, well soluble product. Surprisingly, in contrast to the common observation that high-nuclearity polyoxometalate (POM) clusters are fragmented to form smaller moieties in solution, the low nuclearity [V₆As₈O₂₆]⁴⁻ anion is in situ transformed into the higher nuclearity cluster anions. The obtained products support a conceptually new model that is outlined in this article and that describes a continuous evolution between spherical and cyclic POV assemblies. This new model represents a milestone on the way to rational and designable POV self-assemblies.

Systems-based elucidation of the polyoxovanadate speciation reveals that heterogroup substitution can transform spherical kegginoids into tubular architectures in a programmable manner.