Mechanistic Theoretical Investigation of Self-Discharge Reactions in a Vanadium Redox Flow Battery

Controlling Charge Percolation in Solutions of Metal Redox Active Polymers: Implications of Microscopic Polyelectrolyte Dynamics on Macroscopic Energy Storage

Soluble redox-active polymers (RAPs) enable size-exclusion nonaqueous redox flow batteries (NaRFBs) which promise high energy density. Pendants along the RAPs not only store charge but also engage in electron transfer to varying extents based on their designs. Here, we explore these phenomena in Metal-containing Redox Active Polymers (M-RAPs, M = Ru, Fe, Co). We assess by using cyclic voltammetry and chronoamperometry with ultramicroelectrodes the current response to electrolyte concentration spanning 3 orders of magnitude. Currents scaled as Ru-RAP > Fe-RAP ≫ Co-RAP, consistent with electron self-exchange trends in the small molecule analogues of the MII/III redox pair. Varying the ionic strength of the electrolyte also revealed nonmonotonic behavior, evidencing the impact of polyelectrolytic dynamics on M-RAP redox response. We developed a model to account for the behavior by combining kinetic Monte Carlo and Brownian dynamics near a boundary representing an electrode. While 1D pendant-to-pendant charge transfer along the chain is not a strong function of electrolyte concentration, the microstructure of the RAP at different electrolyte concentrations is decisively impacted, yielding qualitative trends to those observed experimentally. M-RAP size-exclusion NaRFBs using a poly viologen as negolyte varied in average potential with ∼1.54 V for Ru-RAP, ∼1.37 V for Fe-RAP, and ∼0.52 V for Co-RAP. Comparison of batteries at their optimal and suboptimal solution conditions as gauged from analytical experiments showed clear correlations in performance. This work provides a blueprint for understanding the factors underpinning charge transfer in solutions of RAPs for batteries and beyond.

Combined effects of the lewis acidity and electric field of proximal redox innocent metal ions on the redox potential of vanadyl Schiff base complexes: an experimental and theoretical study

The reactivity of biological or synthetic metalloenzymes is modulated in the presence of redox innocent Lewis acidic metal ions as they change the redox potential of the redox active metal ions present in the active site of metalloenzymes. To study this effect, we synthesised a mono-nuclear V(IV) complex (VOL, 1) with an N₂O₄ donor bicompartmental ligand, characterized it by single-crystal X-ray crystallography and recorded its cyclic voltammogram in acetonitrile. The CV revealed a reversible redox process for the V(IV)/V(V) couple. The potential of the V(IV)/V(V) couple shifted to a more positive value when equivalent amounts of Li⁺, Na⁺, K⁺, Mg²⁺, Ca²⁺ and Ba²⁺ ions were added separately to its acetonitrile solution, but the extent of shift for Li⁺ and Mg²⁺ was much less than that of the other metal ions. The guest metal ions except Li⁺ and Mg²⁺ were accommodated in the outer compartment of VOL as confirmed by IR and UV-Vis spectral analysis. Single-crystal structural analysis of [(VOL)KPF₆]₂, (1·K) and [(VOL)Ba(ClO₄)₂(H₂O)]_n, (1·Ba) also confirmed the hetero-metallic adduct formation. The correlation of the shift of the V(IV/V) redox potential with the Lewis acidity of respective metal ions deviated appreciably from linearity. DFT calculations suggest that the shift in potential is probably controlled by local electric fields induced by those ions, as indicated by 2D vector electric field maps.

Effect of nitrogen molecules on the growth kinetics at the interface between a (111) plane of cubic ice and water

The molecular-scale growth kinetics of ice from water in the presence of air molecules are still poorly understood, despite their importance for understanding ice particle formation in nature. In this study, a molecular dynamics simulation is conducted to elucidate the molecular-scale growth kinetics at the interface between a (111) plane of cubic ice and water in the presence of N2 molecules. Two potential models of N2 molecules with and without atomic charges are examined. For both models, N2 molecules bind stably to the interface for a period of 1 ns or longer, and the stability of the binding is higher for the charged model than for the noncharged model. Free-energy surfaces of an N2 molecule along the interface and along an ideal (111) plane surface of cubic ice suggest that for both models, the position where an N2 molecule binds stably is different at the interface and on the ideal plane surface, and the stability of the binding is much higher for the interface than for the ideal plane surface. For both models, stacking-disordered ice grows at the interface, and the formation probability of a hexagonal ice layer in the stacking-disordered ice is higher for the charged model than for the uncharged model. The formation probability for the hexagonal ice layer in the stacking-disordered ice depends not only on the stability of binding but also on the positions where N2 molecules bind on the underlying ice, and the number of N2 molecules that bind stably to the underlying ice.

Stable Binding Conformations of Polymaleic and Polyacrylic Acids at a Calcite Surface in the Presence of Countercations: A Metadynamics Study

Elucidating the stable binding conformations of additives at the surface of CaCO3 crystals is essential to biomineralization, scale inhibition, and materials technology. However, accomplishing this by experimental means is rather difficult. In this study, molecular dynamics simulations based on a metadynamics approach were conducted to elucidate the stable binding conformations of a deprotonated polymaleic acid (PMA) additive and two deprotonated poly(acrylic acid) (PAA) additives with different polymerization degrees in the presence of various countercations at a hydrated calcite (104) surface. The simulated free-energy surfaces suggested the existence of several slightly different stable binding conformations for each additive. The appearance of these distinct binding conformations is speculated to originate from different balances of interactions between the additive, the calcite surface, and the countercations. The binding conformations and binding stabilities at the calcite surface were affected by the countercations, with Ca2+ ions producing a more pronounced effect than Na+ ions. Furthermore, the simulation results suggested that the binding stability at the calcite surface was higher for the PMA additive than for the PAA additives, and the PAA additive with a polymerization degree of 10 displayed a binding stability that was similar to or lower than that of the PAA additive with a polymerization degree of 5. The present simulation method provides a new strategy for analyzing the binding conformations of complex additives at material surfaces, developing additives that stably bind to these surfaces, and designing additives to control crystal growth.

Pathways for the formation of ice polymorphs from water predicted by a metadynamics method

The mechanism of how ice crystal form has been extensively studied by many researchers but remains an open question. Molecular dynamics (MD) simulations are a useful tool for investigating the molecular-scale mechanism of crystal formation. However, the timescale of phenomena that can be analyzed by MD simulations is typically restricted to microseconds or less, which is far too short to explore ice crystal formation that occurs in real systems. In this study, a metadynamics (MTD) method was adopted to overcome this timescale limitation of MD simulations. An MD simulation combined with the MTD method, in which two discrete oxygen–oxygen radial distribution functions represented by Gaussian window functions were used as collective variables, successfully reproduced the formation of several different ice crystals when the Gaussian window functions were set at appropriate oxygen–oxygen distances: cubic ice, stacking disordered ice consisting of cubic ice and hexagonal ice, high-pressure ice VII, layered ice with an ice VII structure, and layered ice with an unknown structure. The free-energy landscape generated by the MTD method suggests that the formation of each ice crystal occurred via high-density water with a similar structure to the formed ice crystal. The present method can be used not only to study the mechanism of crystal formation but also to search for new crystals in real systems.