Modeling energy landscapes of proton motion in nonaqueous, tethered proton wires.

We have modeled structures and energetics of anhydrous proton-conducting wires: tethered hydrogen-bonded chains of the form ···HX···HX···HX···, with functional groups HX = imidazole, triazole, and formamidine; formic, sulfonic, and phosphonic acids. We have applied density functional theory (DFT) to model proton wires up to 19 units long, where each proton carrier is linked to an effective backbone to mimic polymer tethering. This approach allows the direct calculation of hydrogen bond strengths. The proton wires were found to be stabilized by strong hydrogen bonds (up to 50 kJ/mol) whose strength correlates with the proton affinity of HX [related to pK(b)(HX)] and not to pK(a)(HX) as is often assumed. Geometry optimizations and ab initio molecular dynamics near 400 K on imidazole-based proton wires both predict that adding a proton to the end of such wires causes the excess charge to embed into the interior segments of these wires. Proton translocation energy landscapes for imidazole-based wires are sensitive to the imidazole attachment point (head or feet) and to wire architecture (linear or interdigitated). Linear imidazole wires with head-attachment exhibit low barriers for intrawire proton motion, rivaling proton diffusion in liquid imidazole. Excess charge relaxation from the edge of wires is found to be dominated by long-range Grotthuss shuttling for distances as long as 42 Å, especially for interdigitated wires. For imidazole, we predict that proton translocation is controlled by the energetics of desorption from the proton wire, even for relatively long wires (600 imidazole units). Proton desorption energies show no correlation with functional group properties, suggesting that proton desorption is a collective process in proton wires.

PSI3: an open-source Ab Initio electronic structure package.

PSI3 is a program system and development platform for ab initio molecular electronic structure computations. The package includes mature programming interfaces for parsing user input, accessing commonly used data such as basis-set information or molecular orbital coefficients, and retrieving and storing binary data (with no software limitations on file sizes or file-system-sizes), especially multi-index quantities such as electron repulsion integrals. This platform is useful for the rapid implementation of both standard quantum chemical methods, as well as the development of new models. Features that have already been implemented include Hartree-Fock, multiconfigurational self-consistent-field, second-order Møller-Plesset perturbation theory, coupled cluster, and configuration interaction wave functions. Distinctive capabilities include the ability to employ Gaussian basis functions with arbitrary angular momentum levels; linear R12 second-order perturbation theory; coupled cluster frequency-dependent response properties, including dipole polarizabilities and optical rotation; and diagonal Born-Oppenheimer corrections with correlated wave functions. This article describes the programming infrastructure and main features of the package. PSI3 is available free of charge through the open-source, GNU General Public License.

Modeling Proton Transfer in Zeolites: Convergence Behavior of Embedded and Constrained Cluster Calculations.

We have studied the convergence properties of embedded and constrained cluster models of proton transfer in zeolites. We applied density functional theory to describe clusters and ONIOM to perform the embedding. We focused on converging the reaction energy and barrier of the O(1) to O(4) jump in H-Y zeolite as well as vibrational and structural aspects of this jump. We found that using successively larger clusters in vacuo gives convergence of this reaction energy to 14 ± 2 kJ mol(-)(1) and the barrier to 135 ± 5 kJ mol(-)(1) at a cluster size of 5 Å, which contains 11 tetrahedral (Si or Al) atoms. We embedded quantum clusters of various sizes in larger clusters with total radii in the range 7-20 Å, using the universal force field as the lower level of theory in ONIOM. We found convergence to the same values as the constrained clusters, without the use of reactive force fields or periodic boundary conditions in the embedding procedure. For the reaction energy, embedded cluster calculations required smaller clusters than in vacuo calculations, reaching converged reaction energies for quantum systems containing at least 8 tetrahedral atoms. In addition, optimizations on embedded clusters required many fewer cycles, and hence much less CPU time, than did optimizations on comparable constrained clusters.

Are gas-phase reactions of five-coordinate divalent metal ion complexes affected by coordination geometry?

Five-coordinate metal complex ions of the type [ML](2+) [where M = Mn(II), Fe(II), Co(II), Ni(II), Cu(II), Zn(II) and L= 1,9-bis(2-pyridyl)-2,5,8-triazanonane (DIEN-(pyr)(2)) and 1,9-bis(2-imidazolyl)-2,5,8-triazanonane (DIEN-(imi)(2)] have been reacted with acetonitrile in the gas phase using a modified quadrupole ion trap mass spectrometer. The kinetics and thermodynamics of these reactions show that the reactivity of these complexes is affected by metal electronic structure and falls into three groups: Mn(II) and Ni(II) complexes are the most reactive, Fe(II) and Co(II) complexes exhibit intermediate reactivity, and Cu(II) and Zn(II) complexes are the least reactive. To help explain the experimental trends in reactivity, theoretical calculations have been used. Due to the relatively large size of the metal complexes involved, we have utilized a two-layered ONIOM method to perform geometry optimizations and single point energy calculations for the [ML](2+) and [ML + CH(3)CN](2+) systems. The calculations show that the reactant five-coordinate complexes ([ML](2+)) exhibit structures that are slightly distorted trigonal bipyramidal geometries, while the six-coordinate complexes ([ML + CH(3)CN](2+)) have geometries that are close to octahedral. The Delta G values obtained from the ONIOM calculations roughly agree with the experimental data, but the calculations fail to completely explain the trends for the different metal complexes. The failure to consider all possible isomers as well as adequately represent pi-d interactions for the metal complexes is the likely cause of this discrepancy. Using the angular overlap model (AOM) to obtain molecular orbital stabilization energies (MOSE) also fails to reproduce the experimental trends when only sigma interactions are considered but succeeds in explaining the trends when pi interactions are taken into account. These results indicate that the pi-donor character of the CH(3)CN plays a subtle, yet important, role in controlling the reactivity of these five-coordinate complexes. Also, the AOM calculations are consistent with the experimental data when the [ML](2+) complexes have high-spin trigonal bipyramidal configurations. Generally, these results suggest that ion-molecule reactions can be very sensitive to metal complex coordination geometry and thus may have some promise for providing gas-phase coordination structure.