Search

Electronic polarization and charge delocalization are important aspects that affect the charge-transport levels in organic materials. Here, using a quantum mechanical/embedded-charge (QM/EC) approach based on a combination of the long-range corrected ωB97X-D exchange-correlation functional (QM) and charge model 5 (CM5) point-charge model (EC), we evaluate the vertical detachment energies and polarization energies of various sizes of crystalline and amorphous anionic oligoacene clusters. Our results indicate that QM/EC calculations yield vertical detachment energies and polarization energies that compare well with the experimental values obtained from ultraviolet photoemission spectroscopy measurements. In order to understand the effect of charge delocalization on the transport levels, we considered crystalline naphthalene systems with QM regions including one or five-molecules. The results for these systems show that the delocalization and polarization effects are additive; therefore, allowing for electron delocalization by increasing the size of the QM region leads to the additional stabilization of the transport levels.

Modern X-ray spectroscopy has proven itself as a robust tool for probing the electronic structure of atoms in complex environments. Despite working on energy scales that are much larger than those corresponding to nuclear motions, taking nuclear dynamics and the associated nuclear correlations into account may be of importance for X-ray spectroscopy. Recently, we have developed an efficient protocol to account for nuclear dynamics in X-ray absorption and resonant inelastic X-ray scattering spectra [Karsten et al., J. Phys. Chem. Lett. 8, 992 (2017)], based on ground state molecular dynamics accompanied with state-of-the-art calculations of electronic excitation energies and transition dipoles. Here, we present an alternative derivation of the formalism and elaborate on the developed simulation protocol using gas phase and bulk water as examples. The specific spectroscopic features stemming from the nuclear motions are analyzed and traced down to the dynamics of electronic energy gaps and transition dipole correlation functions. The observed tendencies are explained on the basis of a simple harmonic model, and the involved approximations are discussed. The method represents a step forward over the conventional approaches that treat the system in full complexity and provides a reasonable starting point for further improvements.

In the present study, we examined the UV excitations of a newly introduced molecular set, Halons-9, composed of nine gaseous halon molecules. The performance of the density functional-based multi-reference configuration interaction method (DFT/MRCI) and time-dependent density functional theory with CAM-B3LYP functional (TD-CAM-B3LYP) in the computation of singlet and triplet excited states of this set was evaluated against coupled-cluster with singles and doubles (CCSD). Excited states up to the corresponding ionization limits, including both localized and delocalized excitations, have been benchmarked. TD-CAM-B3LYP significantly underestimates excitation energies of the higher mixed valence-Rydberg and Rydberg states, with computed mean absolute deviations from the equation of motion (EOM)-CCSD results 1.06 and 0.76 eV, respectively. DFT/MRCI gives a significantly better description of higher excited states, albeit still poor, compared to the TD-CAM-B3LYP. The mean absolute deviations of mixed valence-Rydberg and Rydberg states from the reference EOM-CCSD values are 0.66 and 0.47 eV, respectively. The performance of DFT/MRCI for description of strongly correlated states with valence-Rydberg mixing is still not satisfactory enough. On the other hand, oscillator strengths of most of singlet states obtained with both methods are close to the EOM-CCSD values. The largest deviations, occurring in the case of several high-lying multiconfigurational states, are of an order of magnitude.

X-ray photoelectron spectra provide a wealth of information on the electronic structure. The extraction of molecular details requires adequate theoretical methods, which in case of transition metal complexes has to account for effects due to the multi-configurational and spin-mixed nature of the many-electron wave function. Here, the restricted active space self-consistent field method including spin-orbit coupling is used to cope with this challenge and to calculate valence- and core-level photoelectron spectra. The intensities are estimated within the frameworks of the Dyson orbital formalism and the sudden approximation. Thereby, we utilize an efficient computational algorithm that is based on a biorthonormal basis transformation. The approach is applied to the valence photoionization of the gas phase water molecule and to the core ionization spectrum of the [Fe(H2O)6](2+) complex. The results show good agreement with the experimental data obtained in this work, whereas the sudden approximation demonstrates distinct deviations from experiments.

Equilibrium thermochemical measurements using the ion mobility drift cell technique have been utilized to investigate the binding energies and entropy changes for the stepwise association of HCN molecules with the pyridine and pyrimidine radical cations forming the C5H5N(+·)(HCN)n and C4H4N2 (+·)(HCN)n clusters, respectively, with n = 1-4. For comparison, the binding of 1-4 HCN molecules to the protonated pyridine C5H5NH(+)(HCN)n has also been investigated. The binding energies of HCN to the pyridine and pyrimidine radical cations are nearly equal (11.4 and 12.0 kcal/mol, respectively) but weaker than the HCN binding to the protonated pyridine (14.0 kcal/mol). The pyridine and pyrimidine radical cations form unconventional carbon-based ionic hydrogen bonds with HCN (CH(δ+)⋯NCH). Protonated pyridine forms a stronger ionic hydrogen bond with HCN (NH(+)⋯NCH) which can be extended to a linear chain with the clustering of additional HCN molecules (NH(+)⋯NCH··NCH⋯NCH) leading to a rapid decrease in the bond strength as the length of the chain increases. The lowest energy structures of the pyridine and pyrimidine radical cation clusters containing 3-4 HCN molecules show a strong tendency for the internal solvation of the radical cation by the HCN molecules where bifurcated structures involving multiple hydrogen bonding sites with the ring hydrogen atoms are formed. The unconventional H-bonds (CH(δ+)⋯NCH) formed between the pyridine or the pyrimidine radical cations and HCN molecules (11-12 kcal/mol) are stronger than the similar (CH(δ+)⋯NCH) bonds formed between the benzene radical cation and HCN molecules (9 kcal/mol) indicating that the CH(δ+) centers in the pyridine and pyrimidine radical cations have more effective charges than in the benzene radical cation.

The bonding energies of proton-bound homodimers BH(+)B were measured by ion mobility equilibrium studies and calculated at the DFT B3LYP/6-311++G** level, for a series of nitrogen heterocyclic molecules (B) with electron-withdrawing in-ring N and on-ring F substituents. The binding energies (ΔH°(dissoc)) of the proton-bound dimers (BH(+)B) vary significantly, from 29.7 to 18.1 kcal/mol, decreasing linearly with decreasing the proton affinity of the monomer (B). This trend differs significantly from the constant binding energies of most homodimers of other organic nitrogen and oxygen bases. The experimentally measured ΔH°(dissoc) for (1,3-diazine)2H(+), i.e., (pyrimidine)2H(+) and (3-F-pyridine)2H(+) are 22.7 and 23.0 kcal/mol, respectively. The measured ΔH°(dissoc) for the pyrimidine(·+)(3-F-pyridine) radical cation dimer (19.2 kcal/mol) is signifcantly lower than that of the proton-bound homodimers of pyrimidine and 3-F-pyridine, reflecting the stronger interaction in the ionic H-bond of the protonated dimers. The calculated binding energies for (1,2-diazine)2H(+), (pyridine)2H(+), (2-F-pyridine)2H(+), (3-F-pyridine)2H(+), (2,6-di-F-pyridine)2H(+), (4-F-pyridine)2H(+), (1,3-diazine)2H(+), (1,4-diazine)2H(+), (1,3,5-triazine)2H(+), and (pentafluoropyridine)2H(+) are 29.7, 24.9, 24.8, 23.3, 23.2, 23.0, 22.4, 21.9, 19.3, and 18.1 kcal/mol, respectively. The electron-withdrawing substituents form internal dipoles whose electrostatic interactions contribute to both the decreased proton affinities of (B) and the decreased binding energies of the protonated dimers BH(+)B. The bonding energies also vary with rotation about the hydrogen bond, and they decrease in rotamers where the internal dipoles of the components are aligned efficiently for inter-ring repulsion. For compounds substituted at the 3 or 4 (meta or para) positions, the lowest energy rotamers are T-shaped with the planes of the two rings rotated by 90° about the hydrogen bond, while the planar rotamers are weakened by repulsion between the ortho hydrogen atoms of the two rings. Conversely, in ortho-substituted (1,2-diazine)2H(+) and (2-F-pyridine)2H(+), attractive interactions between the ortho (C-H) hydrogen atoms of one ring and the electronegative ortho atoms (N or F) of the other ring are stabilizing, and increase the protonated dimer binding energies by up to 4 kcal/mol. In all of the dimers, rotation about the hydrogen bond can involve a 2-4 kcal/mol barrier due to the relative energies of the rotamers.

Equilibrium thermochemical measurements using an ion mobility drift cell technique have been utilized to investigate the binding energies and entropy changes associated with the stepwise hydration of the biologically significant ions pyrimidine radical cation and protonated pyrimidine. The binding energy of the hydrated pyrimidine radical cation is weaker than that of the proton-bound dimer pyrimidineH(+)(H2O) consistent with the formation of a weak carbon-based CH(δ+)··OH2 hydrogen bond (11.9 kcal/mol) and a stronger NH(+)··OH2 hydrogen bond (15.6 kcal/mol), respectively. Other proton-bound dimers such as pyrimidineH(+)(CH3OH) and pyrimidineH(+)(CH3CN) exhibit higher binding energies (18.2 kcal/mol and 22.8 kcal/mol, respectively) due to the higher proton affinities and dipole moments of acetonitrile and methanol as compared to water. The measured collisional cross sections of the proton-bound dimers provide experimental-based support for the DFT calculated structures at the M06-2x/6-311++G (d,p) level. The calculations show that the hydrated pyrimidine radical cation clusters form internally solvated structures in which the water molecules are bonded to the C4N2H4(●+) ion by weak CH(δ+)··OH2 hydrogen bonds. The hydrated protonated pyrimidine clusters form externally solvated structures where the water molecules are bonded to each other and the ion is external to the water cluster. Dissociative proton transfer reactions C4N2H4(●+)(H2O)(n-1) + H2O → C4N2H3(●) + (H2O)(n)H(+) and C4N2H5(+)(H2O)(n-1) + H2O → C4N2H4 + (H2O)(n)H(+) are observed for n ≥ 4 where the reactions become thermoneutral or exothermic. The absence of the dissociative proton transfer reaction within the C4N2H5(+)(CH3CN)n clusters results from the inability of acetonitrile molecules to form extended hydrogen bonding structures such as those formed by water and methanol due to the presence of the methyl groups which block the extension of hydrogen bonding networks.