Relationship between electron-phonon interaction and low-frequency Raman anisotropy in high-mobility organic semiconductors.

Recent theoretical studies have shown that charge transport in high-mobility organic semiconductors is limited by low-frequency vibrations because of strong non-local electron-phonon interaction. Here we investigate two high-electron-mobility organic semiconductors with similar molecular structures but considerably different crystal packings, TCNQ and F2-TCNQ, and reveal the relationship between the experimental low-frequency Raman spectra and the calculated contributions of various vibrational modes to the electron-phonon interaction. We suggest that the combination of Raman spectroscopy with solid-state DFT is a powerful tool for probing electron-phonon interaction and focused search for high-mobility organic semiconductors.

The structure and IR signatures of the arginine-glutamate salt bridge. Insights from the classical MD simulations.

Salt bridges and ionic interactions play an important role in protein stability, protein-protein interactions, and protein folding. Here, we provide the classical MD simulations of the structure and IR signatures of the arginine (Arg)-glutamate (Glu) salt bridge. The Arg-Glu model is based on the infinite polyalanine antiparallel two-stranded ß-sheet structure. The 1 µs NPT simulations show that it preferably exists as a salt bridge (a contact ion pair). Bidentate (the end-on and side-on structures) and monodentate (the backside structure) configurations are localized [Donald et al., Proteins 79, 898-915 (2011)]. These structures are stabilized by the short (+)N-H⋯O(-) bonds. Their relative stability depends on a force field used in the MD simulations. The side-on structure is the most stable in terms of the OPLS-AA force field. If AMBER ff99SB-ILDN is used, the backside structure is the most stable. Compared with experimental data, simulations using the OPLS all-atom (OPLS-AA) force field describe the stability of the salt bridge structures quite realistically. It decreases in the following order: side-on > end-on > backside. The most stable side-on structure lives several nanoseconds. The less stable backside structure exists a few tenth of a nanosecond. Several short-living species (solvent shared, completely separately solvated ionic groups ion pairs, etc.) are also localized. Their lifetime is a few tens of picoseconds or less. Conformational flexibility of amino acids forming the salt bridge is investigated. The spectral signature of the Arg-Glu salt bridge is the IR-intensive band around 2200 cm(-1). It is caused by the asymmetric stretching vibrations of the (+)N-H⋯O(-) fragment. Result of the present paper suggests that infrared spectroscopy in the 2000-2800 frequency region may be a rapid and quantitative method for the study of salt bridges in peptides and ionic interactions between proteins. This region is usually not considered in spectroscopic studies of peptides and proteins.

Intermolecular hydrogen bond energies in crystals evaluated using electron density properties: DFT computations with periodic boundary conditions.

The hydrogen bond (H-bond) energies are evaluated for 18 molecular crystals with 28 moderate and strong O-H···O bonds using the approaches based on the electron density properties, which are derived from the B3LYP/6-311G** calculations with periodic boundary conditions. The approaches considered explore linear relationships between the local electronic kinetic G(b) and potential V(b) densities at the H···O bond critical point and the H-bond energy E(HB). Comparison of the computed E(HB) values with the experimental data and enthalpies evaluated using the empirical correlation of spectral and thermodynamic parameters (Iogansen, Spectrochim. Acta Part A 1999, 55, 1585) enables to estimate the accuracy and applicability limits of the approaches used. The V(b)-E(HB) approach overestimates the energy of moderate H-bonds (E(HB) < 60 kJ/mol) by ~20% and gives unreliably high energies for crystals with strong H-bonds. On the other hand, the G(b)-E(HB) approach affords reliable results for the crystals under consideration. The linear relationship between G(b) and E(HB) is basis set superposition error (BSSE) free and allows to estimate the H-bond energy without computing it by means of the supramolecular approach. Therefore, for the evaluation of H-bond energies in molecular crystals, the G(b) value can be recommended to be obtained from both density functional theory (DFT) computations with periodic boundary conditions and precise X-ray diffraction experiments.

QTAIM study of strong H-bonds with the O-H...a fragment (A=O, N) in three-dimensional periodical crystals.

The relationship between the d(H...A) distance (A=O, N) and the topological properties at the H...A bond critical point of 37 strong (short) hydrogen bonds occurring in 26 molecular crystals are analyzed using the quantum theory of atoms in molecules (QTAIM). Ground-state wave functions of the three-dimensional periodical structures representing the accurate experimental geometries calculated at the B3LYP/6-31G** level of approximation were used to obtain the QTAIM electron density characteristics. The use of an electron-correlated method allowed us to reach the quantitatively correct values of electron density rhob at the H...A bond critical point. However, quite significant differences can appear for small absolute values of the Laplacian (<0.5 au). The difference between the H...O and H...N interactions is described using the rhob versus d(H...A) dependence. It is demonstrated that the values of parameters in this dependence are defined by the nature of the heavy atom forming the H...A bond. An intermediate (or transit) region separating the shared and closed-shell interactions is observed for the H-bonded crystals in which the bridging proton can move from one heavy atom to another. The crystalline environment changes the location of the bridging proton in strong H-bonded systems; however, the d(O-H)/d(H...O) ratio is approximately the same for both the gas-phase complexes and molecular crystals with a linear or near-linear O-H...O bond.

Molecular simulations of outersphere reorganization energies in polar and quadrupolar solvents. The case of intramolecular electron and hole transfer.

Outersphere reorganization energies (lambda) for intramolecular electron and hole transfer are studied in anion- and cation-radical forms of complex organic substrates (p-phenylphenyl-spacer-naphthyl) in polar (water, 1,2-dichloroethane, tetrahydrofuran) and quadrupolar (supercritical CO2) solvents. Structure and charge distributions of solute molecules are obtained at the HF/6-31G(d,p) level. Standard Lennard-Jones parameters for solutes and the nonpolarizable simple site-based models of solvents are used in molecular dynamics (MD) simulations. Calculation of lambda is done by means of the original procedure, which treats electrostatic polarization of a solvent in terms of a usual nonpolarizable MD scheme supplemented by scaling of reorganization energies at the final stage. This approach provides a physically relevant background for separating inertial and inertialless polarization responses by means of a single parameter epsilon(infinity), optical dielectric permittivity of the solvent. Absolute lambda values for hole transfer in 1,2-dichloroethane agree with results of previous computations in terms of the different technique (MD/FRCM, Leontyev, I. V.; et al. Chem. Phys. 2005, 319, 4). Computed lambda values for electron transfer in tetrahydrofuran are larger than the experimental values by ca. 2.5 kcal/mol; for the case of hole transfer in 1,2-dichloroethane the discrepancy is of similar magnitude provided the experimental data are properly corrected. The MD approach gives nonzero lambda values for charge-transfer reaction in supercritical CO2, being able to provide a uniform treatment of nonequilibrium solvation phenomena in both quadrupolar and polar solvents.