Qadir K. Timerghazin

Key Role of End-Capping Groups in Optoelectronic Properties of Poly-p-phenylene Cation Radicals.

Poly-p-phenylenes (PPs) are prototype systems for understanding the charge transport in π-conjugated polymers. In a combined computational and experimental study, we demonstrate that the smooth evolution of redox and optoelectronic properties of PP cation radicals toward the polymeric limit can be significantly altered by electron-donating iso-alkyl and iso-alkoxy end-capping groups. A multiparabolic model (MPM) developed and validated here rationalizes this unexpected effect by interplay of the two modes of hole stabilization: due to the framework of equivalent p-phenylene units and due to the electron-donating end-capping groups. A symmetric, bell-shaped hole in unsubstituted PPs becomes either slightly skewed and shifted toward an end of the molecule in iso-alkyl-capped PPs or highly deformed and concentrated on a terminal unit in PPs with strongly electron-donating iso-alkoxy capping groups. The MPM shows that the observed linear 1/n evolution of the PP cation radical properties toward the polymer limit originates from the hole stabilization due to the growing chain of p-phenylene units, while shifting of the hole toward electron-donating end-capping groups leads to early breakdown of these 1/n dependencies. These insights, along with the readily applicable and flexible multistate parabolic model, can guide studies of complex donor–spacer–acceptor systems and doped molecular wires to aid the design of the next generation materials for long-range charge transport and photovoltaic applications.

Kinetics and mechanism of S-nitrosothiol acid-catalyzed hydrolysis: sulfur activation promotes facile NO+ release.

The denitrosation of three primary S-nitrosothiols (RSNO; S-nitrosocysteine, S-nitroso-N-acetylcysteine, and S-nitrosoglutathione) and two tertiary RSNOs (S-nitrosopenicillamine and S-nitroso-N-acetylpenicillamine) was investigated in 3.75 M H(2)SO(4) to probe the mechanism of acid-catalyzed RSNO hydrolysis and its dependence on RSNO structure. This reversible reaction was forced to proceed in the denitrosation direction by trapping the nitrosating agent with HN(3). The primary RSNOs exhibited hydrolysis k(obs) values of ∼2 × 10(-4) s(-1), and the tertiary RSNO k(obs) values were an order of magnitude higher. Product analysis by HPLC revealed that the parent thiols (RSHs) were formed in 90-100% yield on 79-99% RSNO denitrosation. Possible hydrolysis mechanisms were studied computationally at the CBS-QB3 level using S-nitrosomethanethiol (MeSNO) as a model RSNO. Consideration of RSNOs as a combination of conventional R-S-N═O, zwitterionic R-S(+)═N-O(-), and RS(-)/NO(+) ion-pair resonance structures was key in understanding the mechanistic details of acid-catalyzed hydrolysis. Protonation of the S-nitroso oxygen or nitrogen activates the sulfur and nucleophilic attack by H(2)O at this atom leads to the formation of the sulfoxide-protonated N-hydroxysulfinamide, MeS(+)(OH)NHOH, with barriers of 19 and 29 kcal/mol, respectively. Proton loss and reprotonation at the nitrogen lead to secondary hydrolysis that produces the sulfinic acid MeS(═O)OH and NH(2)OH. Notably, no low-energy RSNO hydrolysis pathway for HNO release was found in the computational analysis. Protonation of the S-nitroso sulfur gives rise to NO(+) release with a low activation barrier (ΔH(double dagger)(calc) ≈ 6 kcal/mol) and the formation of MeSH in agreement with experiment. The experimental k(obs) can be expressed as K(a)k(1), where K(a) is the acid dissociation constant for protonation of the S-nitroso sulfur and k(1) the pseudo-first-order hydrolysis rate constant. Given the low ΔH(double dagger)(calc) for denitrosation of the S-protonated isomer, the observed slow rates of acid-catalyzed RSNO hydrolysis must be controlled by the magnitude of K(a). The 10-fold higher K(a) calculated for Me(3)CS(H(+))NO (∼10(-15)) compared to MeS(H(+))NO (10(-16)) is consistent with the order of magnitude larger k(obs) reported here for the tertiary vs primary RSNOs.

Theoretical investigation of charge transfer to solvent in photoexcited iodide–acetonitrile clusters

Abstract Photoexcitation of the iodide–acetonitrile complex leading to charge-transfer-to-solvent precursor states has been investigated with electronic structure theory calculations. The resulting vertical excitation energies are in good agreement with experimental results. The picture that emerges from the calculations is that, upon photoexcitation, an electron is transferred from one of the iodide valence p orbitals to the acetonitrile molecule, where it is stabilized by interactions with the acetonitrile dipole and the electrophilic hydrogen atoms of the methyl group. A similar picture applies to the photoexcitation of iodide–water clusters, but charge transfer to solvent appears much more obvious for iodide–acetonitrile clusters.

Unprecedented External Electric Field Effects on S-Nitrosothiols: Possible Mechanism of Biological Regulation?

Reactions of S-nitrosothiols (RSNOs), ubiquitous carriers of nitric oxide NO and its physiological activity, are tightly regulated in biological systems, but the mechanisms of this regulation are not well understood. Here, we computationally demonstrate that RSNO properties can be dramatically altered by biologically accessible external electric fields (EEFs) by modulation of the two minor antagonistic resonance structures of RSNOs, which have opposite formal charge distributions and bonding patterns. As these resonance contributions relate to the two competing modes of RSNO reactivity with nucleophiles, via N- or S-atom directed nucleophilic attack, EEFs are predicted to be efficient in controlling biologically important RSNO reactions with thiols. For instance, EEF catalysis might be one of the mechanisms behind the high selectivity of protein trans-S-nitrosation reactions, or putative nitroxyl HNO formation via RSNO S-thiolation reactions.

Non-nuclear attractor of electron density as a manifestation of the solvated electron

Two or more polar molecules can trap an excess electron either in a dipole-bound fashion where it is located outside of the cluster (dipole-bound electron) or inside the cluster (solvated electron). The topology of the electron density in dipole-bound and solvated-electron clusters has been examined for the paradigm (HF)3- cluster. As spatial confinement of the excess electron increases, a non-nuclear maximum (or attractor) of the electron density eventually forms, which suggests that the solvated electron can be described as a topological atom with its own set of physicochemical properties.

Can a Dipole-Bound Electron Form a Pseudo-Atom? An Atoms-In-Molecules Study of the Hydrated Electron

Non-nuclear local maxima, or attractors, of electron density are a rare but very interesting feature of the electron density distribution in molecules and solids. Recently, non-nuclear attractors (NNAs) and the corresponding pseudoatoms of electron density have been identified with the quantum theory of atoms in molecules for some anionic clusters formed by several polar solvent molecules and an excess electron bound in either a solvated-electron or dipole-bound fashion. This contribution reports a detailed study of the topology of the electron density for a series of dipole-bound water cluster anions, as calculated with Hartree-Fock, Møller-Plesset perturbation theory, and coupled-cluster methods together with basis sets augmented with extra diffuse basis functions to accommodate the excess electron. For dipole-bound clusters, electron densities obtained with insufficient inclusion of electron correlation effects and tight basis sets feature a well-pronounced pseudoatom due to the excess electron, which ultimately disappears when a higher level of electronic structure theory and a more diffuse basis set are used. On the other hand, for solvated-electron clusters, where the excess electron is surrounded by solvent molecules, the existence of NNAs does not seem to be an artifact of the method employed, but rather a genuine feature of the electron density distribution. Pseudoatoms of electron density thus appear to be an exclusive feature of confined environments and are unlikely to be found on the tip of a cluster dipole or on solid surfaces.

Genetic Algorithm Optimization of Point Charges in Force Field Development: Challenges and Insights

Evolutionary methods, such as genetic algorithms (GAs), provide powerful tools for optimization of the force field parameters, especially in the case of simultaneous fitting of the force field terms against extensive reference data. However, GA fitting of the nonbonded interaction parameters that includes point charges has not been explored in the literature, likely due to numerous difficulties with even a simpler problem of the least-squares fitting of the atomic point charges against a reference molecular electrostatic potential (MEP), which often demonstrates an unusually high variation of the fitted charges on buried atoms. Here, we examine the performance of the GA approach for the least-squares MEP point charge fitting, and show that the GA optimizations suffer from a magnified version of the classical buried atom effect, producing highly scattered yet correlated solutions. This effect can be understood in terms of the linearly independent, natural coordinates of the MEP fitting problem defined by the eigenvectors of the least-squares sum Hessian matrix, which are also equivalent to the eigenvectors of the covariance matrix evaluated for the scattered GA solutions. GAs quickly converge with respect to the high-curvature coordinates defined by the eigenvectors related to the leading terms of the multipole expansion, but have difficulty converging with respect to the low-curvature coordinates that mostly depend on the buried atom charges. The performance of the evolutionary techniques dramatically improves when the point charge optimization is performed using the Hessian or covariance matrix eigenvectors, an approach with a significant potential for the evolutionary optimization of the fixed-charge biomolecular force fields.

Photoexcitation and charge-transfer-to-solvent relaxation dynamics of the I(-)(CH3CN) complex.

Photoexcitation of iodide-acetonitrile clusters, I(-)(CH3CN)n, to the charge-transfer-to-solvent (CTTS) state and subsequent cluster relaxation could result in the possible formation of cluster analogues of the bulk solvated electron. In this work, the relaxation process of the CTTS excited iodide-acetonitrile binary complex, [I(-)(CH3CN)]*, is investigated using rigorous ab initio quantum chemistry calculations and direct-dynamics simulations to gain insight into the role and motion of iodine and acetonitrile in the relaxation of CTTS excited I(-)(CH3CN)n. Computed potential energy curves and profiles of the excited electron vertical detachment energy for [I(-)(CH3CN)]* along the iodine-acetonitrile distance coordinate reveal for the first time significant dispersion effects between iodine and the excited electron, which can have a significant stabilizing effect on the latter. Results of direct-dynamics simulations demonstrate that [I(-)(CH3CN)]* undergoes dissociation to iodine and acetonitrile fragments, resulting in decreased stability of the excited electron. The present work provides strong evidence of solvent translational motion and iodine ejection as key aspects of the early time relaxation of CTTS excited I(-)(CH3CN)n that can also have a substantial impact on the subsequent electron solvation processes and further demonstrates that intricate details of the relaxation process of CTTS excited iodide-polar solvent molecule clusters make it heavily solvent-dependent.

Accurate ab initio potential for the Na+⋯I• complex

High-level ab initio calculations employing the multireference configuration interaction and coupled clusters methods with a correlation-consistent sequence of basis sets have been used to obtain accurate potential energy curves for the complex of the sodium cation with the iodine atom. Potential curves for the first two electronic Lambda-S states have very different characters: the potential for the 2pi state has a well depth of approximately 10 kcal/mol, while the 2sigma state is essentially unbound. This difference is rationalized in terms of the anisotropic interaction of the quadrupole moment of the iodine atom with the sodium cation, which is stabilizing in the case of the 2pi state and destabilizing in the case of the 2sigma state. The effects of spin-orbit coupling have been accounted for with both ab initio and semiempirical approaches, which have been found to give practically the same results. Inclusion of spin-orbit interactions does not affect the X(omega = 32) ground state, which retains its 2pi character, but it results in two omega = 12 spin-orbit states, with mixed 2sigma and 2pi characters and binding energies roughly half of that of the ground spin-orbit state. Complete basis set (CBS) extrapolations of potential curves, binding energies, and equilibrium geometries were also performed, and used to calculate a number of rovibronic parameters for the Na+...I* complex and to parameterize model potentials. The final CBS-extrapolated and zero-point vibrational energy-corrected binding energy is 10.2 kcal/mol. Applications of the present results for simulations of NaI photodissociation femtosecond spectroscopy are discussed.

Asymmetric solvation revisited: The importance of hydrogen bonding in iodide–acetonitrile clusters

The concept of asymmetric solvation in I−(CH3CN)2 clusters [C. E. H. Dessent, C. G. Bailey, and M. A. Johnson, J. Chem. Phys. 103, 2006 (1995)] was revisited in light of quantum chemistry calculations of the cluster structures, energetics and charge-transfer-to-solvent transition energies. In addition to the two “classical” linear structures—with a symmetrically and an asymmetrically solvated iodide—that were postulated previously, a nonlinear hydrogen-bonded isomer was also found. The calculated binding energies and vertical excitation energies for the linear symmetric and nonlinear hydrogen-bonded I−(CH3CN)2 isomers are in very good agreement with experimental data, while the linear asymmetric structure appears unlikely to be involved in experiments. The importance of hydrogen bonding in halide–acetonitrile clusters manifests itself, not only in the fact that a stable nonlinear hydrogen-bonded isomer exists, but also in slight geometrical distortions of the “classical” linear cluster structures.