RESUMEN
The definition of a hydrogen bond (H-bond) is intimately related to the topological and dynamic properties of the hydrogen bond network within liquid water. The development of a universal H-bond definition for water is an active area of research as it would remove many ambiguities in the network properties that derive from the fixed definition employed to assign whether a water dimer is hydrogen bonded. This work investigates the impact that an electronic-structure based definition, an energetic, and a geometric definition of the H-bond has upon both topological and dynamic network behavior of simulated water. In each definition, the use of a cutoff (either geometric or energetic) to assign the presence of a H-bond leads to the formation of transiently bonded or broken dimers, which have been quantified within the simulation data. The relative concentration of transient species, and their duration, results in two of the three definitions sharing similarities in either topological or dynamic features (H-bond distribution, H-bond lifetime, etc.), however no two definitions exhibit similar behavior for both classes of network properties. In fact, two networks with similar local network topology (as indicated by similar average H-bonds) can have dramatically different global network topology (as indicated by the defect state distributions) and altered H-bond lifetimes. A dynamics based correction scheme is then used to remove artificially transient H-bonds and to repair artificially broken bonds within the network such that the corrected network exhibits the same structural and dynamic properties for two H-bond definitions (the properties of the third definition being significantly improved). The algorithm described represents a significant step forward in the development of a unified hydrogen bond network whose properties are independent of the original hydrogen bond definition that is employed.
RESUMEN
Many intermolecular chemical interactions persist across length and timescales and can be considered to form a "network" or "graph." Obvious examples include the hydrogen bond networks formed by polar solvents such as water or alcohols. In fact, there are many similarities between intermolecular chemical networks like those formed by hydrogen bonding and the complex and distributed networks found in computer science. Contemporary network analyses are able to dissect the complex local and global changes that occur within the network over multiple time and length scales. This work discusses the ChemNetworks software, whose purpose is to process Cartesian coordinates of chemical systems into a network/graph formalism and apply topological network analyses that include network neighborhood, the determination of geodesic paths, the degree census, direct structural searches, and the distribution of defect states of network. These properties can help to understand the network patterns and organization that may influence physical properties and chemical reactivity. The focus of ChemNetworks is to quantitatively describe intermolecular chemical networks of entire systems at both the local and global levels and as a function of time. The code is highly general, capable of converting a wide variety of systems into a chemical network formalism, including complex solutions, liquid interfaces, or even self-assemblies.
RESUMEN
The ruthenium "blue dimer" [(bpy)(2)Ru(III)OH(2)](2)O(4+) is best known as the first well-defined molecular catalyst for water oxidation. It has been subject to numerous computational studies primarily employing density functional theory. However, those studies have been limited in the functionals, basis sets, and continuum models employed. The controversy in the calculated electronic structure and the reaction energetics of this catalyst highlights the necessity of benchmark calculations that explore the role of density functionals, basis sets, and continuum models upon the essential features of blue-dimer reactivity. In this paper, we report Kohn-Sham complete basis set (KS-CBS) limit extrapolations of the electronic structure of "blue dimer" using GGA (BPW91 and BP86), hybrid-GGA (B3LYP), and meta-GGA (M06-L) density functionals. The dependence of solvation free energy corrections on the different cavity types (UFF, UA0, UAHF, UAKS, Bondi, and Pauling) within polarizable and conductor-like polarizable continuum model has also been investigated. The most common basis sets of double-zeta quality are shown to yield results close to the KS-CBS limit; however, large variations are observed in the reaction energetics as a function of density functional and continuum cavity model employed.
RESUMEN
Density functional theory (DFT) has been used to investigate the plausibility of water addition to the simple mononuclear ruthenium complexes, [(NH(3))(3)(bpy)RuâO](2+/3+) and [(NH(3))(3)(bpy)RuOH](3+), in which the OH fragment adds to the 2,2'-bipyridine (bpy) ligand. Activation of bpy toward water addition has frequently been postulated within the literature, although there exists little definitive experimental evidence for this type of "covalent hydration". In this study, we examine the energetic dependence of the reaction upon metal oxidation state, overall spin state of the complex, as well as selectivity for various positions on the bipyridine ring. The thermodynamic favorability is found to be highly dependent upon all three parameters, with free energies of reaction that span favorable and unfavorable regimes. Aqueous addition to [(NH(3))(3)(bpy)RuâO](3+) was found to be highly favorable for the S = 1/2 state, while reduction of the formal oxidation state on the metal center makes the reaction highly unfavorable. Examination of both facial and meridional isomers reveals that when bipyridine occupies the position trans to the ruthenyl oxo atom, reactivity toward OH addition decreases and the site preferences are altered. The electronic structure and spectroscopic signatures (EPR parameters and simulated spectra) have been determined to aid in recognition of "covalent hydration" in experimental systems. EPR parameters are found to uniquely characterize the position of the OH addition to the bpy as well as the overall spin state of the system.