Nucleation mechanisms and speciation of metal oxide clusters

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.

Revealing the Electrocatalytic Self-Assembly Route from Building Blocks into Giant Mo-Blue Clusters

The assembly of single-core molybdate into hundreds of cores of giant molybdenum blue (Mo-blue) clusters has remained a long-standing unresolved scientific puzzle. To reveal this fascinating self-assembly behavior, we demonstrate an aqueous flowing in-operando Raman characterization system to capture the building blocks' evolution from the "black box" reaction process. We successfully visualized the sequential transformation of Na2MoO4 into Mo7O246- ({Mo7}), high nuclear Mo36O1128- ({Mo36}) cluster, and finally polymerization product of [H6K2Mo3O12(SO4)]n ({Mo3(SO4)}n) during the H2SO4 acidification. Notably, the facile conversion of {Mo3(SO4)}n back to the {Mo36} cluster by simple dilution is also discovered. Furthermore, we identified {Mo36} and {Mo3(SO4)}n as exclusive precursors responsible for driving the electrochemical self-assembly of {Mo154} and {Mo102}, respectively. The study also unravels a pivotal intermediate, the pentagonal reduced state fragment [H18MoVI4MoVO24]-, originating from {Mo36}, which catalyzes the autocatalytic self-assembly of {Mo154} with electron and proton injection during electrochemical processes. Concurrently, {Mo3(SO4)}n serves as the indispensable precursor for {Mo102} formation, generating sulfation pentagon building blocks of [H2Na2O2(H4MoVMoVI4O16SO4)4]4- that facilitate the consecutive assembly of giant {Mo102} sphere clusters. As a result, a complete elucidation of the assembly pathway of giant Mo-blue clusters derived from single-core molybdate was obtained, and H+/e- redox couple is revealed to play a critical role in catalyzing the deassembly of the precursor, leading to the formation of thermodynamically stable intermediates essential for further self-assembly of reduced state giant clusters.

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.

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.

Understanding Metal-Organic Framework Nucleation from a Solution with Evolving Graphs

A mechanistic understanding of metal-organic framework (MOF) synthesis and scale-up remains underexplored due to the complex nature of the interactions of their building blocks. In this work, we investigate the collective assembly of building units at the early stages of MOF nucleation, using MIL-101(Cr) as a prototypical example. Using large-scale molecular dynamics simulations, we observe that the choice of solvent (water and N,N-dimethylformamide), the introduction of ions (Na⁺ and F^-) and the relative populations of MIL-101(Cr) half-secondary building unit (half-SBU) isomers have a strong influence on the cluster formation process. Additionally, the shape, size, nucleation and growth rates, crystallinity, and short and long-range order largely vary depending on the synthesis conditions. We evaluate these properties as they naturally emerge when interpreting the self-assembly of MOF nuclei as the time evolution of an undirected graph. Solution-induced conformational complexity and ionic concentration have a dramatic effect on the morphology of clusters emerging during assembly. While pure solvents lead to the rapid formation of a small number of large clusters, the presence of ions in aqueous solutions results in smaller clusters and slower nucleation. This diversity is captured by the key features of the graph representation. Principle component analysis on graph properties reveals that only a small number of molecular descriptors is needed to deconvolute MOF self-assembly. Descriptors such as the average coordination number between half-SBUs and fractal dimension are of particulalr interest as they can be can be followed experimentally by techniques like by time-resolved spectroscopy. Ultimately, graph theory emerges as an approach that can be used to understand complex processes revealing molecular descriptors accessible by both simulation and experiment.

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.

A new route to polyoxometalates via mechanochemistry

Mechanochemistry offers a new route to polyoxometalates (POMs) under mild conditions. The molybdenum isoPOM heptamolybdate and the molybdenum heteroPOMs of the Strandberg- and Keggin-type could be achieved from grinding together molybdenum oxide, potassium or ammonium carbonate and phosphate. The reactions were controlled by the stoichiometric ratio of the starting materials and the liquid used, with reaction times between 30 min and 3 h. In situ investigations of the syntheses reveal the formation of intermediates during the reactions. Their identification helps explaining the mechanism of formation of the intermediates as well as the final POMs.

Under mild conditions, molybdenum POMs could be achieved mechanochemically from simple building blocks, within short reaction times. In situ investigations reveal the formation of intermediates and help explaining the mechanism behind the reaction.

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

Understanding and controlling aqueous speciation of metal oxides are key for the discovery and development of novel materials, and challenge both experimental and computational approaches. Here we present a computational method, called POMSimulator, which is able to predict speciation phase diagrams (Conc. vs pH) for multispecies chemical equilibria in solution, and which we apply to molybdenum and tungsten isopolyoxoanions (IPAs). Starting from the MO₄ monomers, and considering dimers, trimers, and larger species, the chemical reaction networks involved in the formation of [H₃₂Mo₃₆O₁₂₈]^8- and [W₁₂O₄₂]^12- are sampled in an automatic manner. This information is used for setting up ∼10⁵ speciation models, and from there, we generate the speciation phase diagrams, which show an insightful picture of the behavior of IPAs in aqueous solution. Furthermore, we predict the values of 107 formation constants for a diversity of molybdenum and tungsten molecular oxides. Among these species, we could include several pentagonal-shaped species and very reactive tungsten intermediates as well. Last but not least, the calibration employed for correcting the density functional theory (DFT) Gibbs energies is remarkably similar for both metals, which suggests that a general rule might exist for correcting computed free energies for other metals.

Atomically precise vanadium-oxide clusters

Polyoxovanadate (POV) clusters are an important subclass of polyoxometalates with a broad range of molecular compositions and physicochemical properties. One relatively underdeveloped application of these polynuclear assemblies involves their use as atomically precise, homogenous molecular models for bulk metal oxides. Given the structural and electronic similarities of POVs and extended vanadium oxide materials, as well as the relative ease of modifying the homogenous congeners, investigation of the chemical and physical properties of pristine and modified cluster complexes presents a method toward understanding the influence of structural modifications (e.g. crystal structure/phase, chemical makeup of surface ligands, elemental dopants) on the properties of extended solids. This review summarises recent advances in the use of POV clusters as atomically precise models for bulk metal oxides, with particular focus on the assembly of vanadium oxide clusters and the consequences of altering the molecular composition of the assembly via organofunctionalization and the incorporation of elemental "dopants".