Time‐dependent density functional theory study on direction‐dependent electron and hole transfer processes in molecular systems

We report on real‐time time‐dependent density functional theory calculations on direction‐dependent electron and hole transfer processes in molecular systems. As a model system, we focus on α‐sulfur. It is shown that time scale of the electron transfer process from a negatively charged S₈ molecule to a neighboring neutral monomer is comparable to that of a strong infrared‐active molecular vibrations of the dimer with one negatively charged monomer. This results in a strong coupling between the electrons and the nuclei motion which eventually leads to S₈ ring opening before the electron transfer process is completed. The open‐ring structure is found to be stable. The similar infrared‐active peak in the case of hole transfer, however, is shown to be very weak and hence no significant scattering by the nuclei is possible. The presented approach to study the charge transfer processes in sulfur has direct applications in the increasingly growing research field of charge transport in molecular systems. © 2017 Wiley Periodicals, Inc.     

Theoretical study on the influence of ancillary ligand on the spectroscopic properties and electronic structures of phosphorescent Pt(II) complexes

The geometries, energies, and electronic properties of a series of phosphorescent Pt(II) complexes including FPt, CFPt, COFPt, and NFPt have been characterized within density functional theory DFT calculations which can reproduce and rationalize experimental results. The properties of excited‐states of the Pt(II) complexes were characterized by configuration interaction with singles (CIS) method. The ground‐ and excited‐state geometries were optimized at the B3LYP/LANL2DZ and CIS/LANL2DZ levels, respectively. In addition, we also have performed a triplet UB3LYP optimization for complex FPt and compared it with CIS method in the emission properties. The datum (562.52 nm) of emission wavelength for complex FPt, which were computed based on the triplet UB3LYP optimization excited‐state geometry, is not agreement with the experiment value (500 nm). The absorption and phosphorescence wavelengths were computed based on the optimized ground‐ and excited‐state geometries, respectively, by the time‐dependent density functional theory (TD‐DFT) methods. The results revealed that the nature of the substituent at the phenylpyridine ligand can influence the distributions of HOMO and LUMO and their energies. Moreover, the auxiliary ligand pyridyltetrazole can make the molecular structure present a solid geometry. In addition, the charge transport quality has been estimated approximately by the predicted reorganization energy (λ). Our result also indicates that the substitute groups and different auxiliary ligand not only change the nature of transition but also affect the rate and balance of charge transfer. By summarizing the results, we can conclude that the NFPt is good OLED materials with a solid geometry and a balanced charge transfer rate. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Tuning the Protein‐Induced Absorption Shifts of Retinal in Engineered Rhodopsin Mimics

Rational design of light‐capturing properties requires understanding the molecular and electronic structure of chromophores in their native chemical or biological environment. We employ here large‐scale quantum chemical calculations to study the light‐capturing properties of retinal in recently designed human cellular retinol binding protein II (hCRBPII) variants (Wang et al. Science, 2012 , 338, 1340–1343). Our calculations show that these proteins absorb across a large part of the visible spectrum by combined polarization and electrostatic effects. These effects stabilize the ground or excited state energy levels of the retinal by perturbing the Schiff‐base or β‐ionone moieties of the chromophore, which in turn modulates the amount of charge transfer within the molecule. Based on the predicted tuning principles, we design putative in silico mutations that further shift the absorption properties of retinal in hCRBPII towards the ultraviolet and infrared regions of the spectrum.     

Influence of the Temperature on the Crystal Structure of Poly(p‐phenylene sulfide): A Study by X‐Ray Measurements and Molecular Mechanics

The crystal structure of poly(p‐phenylene sulfide) (PPS) has been investigated by X‐ray analysis on fiber samples and by molecular mechanics calculations over a wide range of temperatures, from 0 K to 548 K, showing that all the structural parameters remain substantially constant. The thermal expansion coefficients of the a and b axes have been evaluated. Structural parameters experimentally obtained at the various temperatures have been used in calculations by molecular mechanics. The crystal structures calculated by various methods and using several potential functions are in very good agreement. New parameters are proposed for the nonbonded terms of the potential functions.     

First‐Principles Investigation of Anisotropic Electron and Hole Mobility in Heterocyclic Oligomer Crystals

Based on quantum chemistry calculations combined with the Marcus–Hush electron transfer theory, we investigated the charge‐transport properties of oligothiophenes (nTs) and oligopyrroles (nPs) (n=6, 7, 8) as potential p‐ or n‐type organic semiconductor materials. The results of our calculations indicate that 1) the nPs show intrinsic hole mobilities as high as or even higher than those of nTs, and 2) the vertical ionization potentials (VIPs) of the nPs are about 0.6–0.7 eV smaller than the corresponding VIPs of the nTs. Based on their charge‐transport ability and hole‐injection efficiency, the nPs have potential as p‐type organic semiconducting materials. Furthermore, it was also found that the maximum values of the electron‐transfer mobility for the nTs are larger by one‐to‐two orders of magnitude than the corresponding maximum values of hole‐transfer mobility, which suggests that the nTs have the potential to be developed as promising n‐type organic semiconducting materials owing to their electron mobility.     

Theoretical study on the spectroscopic properties and electronic structures of heteroleptic phosphorescent Ir(III) complexes

The geometries, spectroscopic and electronic structures properties of a series of heteroleptic phosphorescent Ir(III) complexes including N981, N982, N983, N984 have been characterized by density functional theory calculations. The excited‐state properties of the Ir(III) complexes have been characterized by CIS method. The ground‐ and excited‐state geometries were optimized at the B3LYP/LANL2DZ and CIS/LANL2DZ levels, respectively. By using the time‐dependent density functional theory method, the absorption and phosphorescence spectra were calculated based on the optimized ground‐ and excited‐state geometries, respectively. The results show that the absorption and emission data agree well with the corresponding experimental results. The calculated results also revealed that the nature of the substituent at the 4‐position of the pyridyl moiety can influence the distributions of HOMO and LUMO and their energies. In addition, the charge transport quality has been estimated approximately by the calculated reorganization energy (λ). Our result also indicates that the positions of the substitute groups not only change the transition characters but also affect the charge transfer rate and balance, and complex N982 is a very good charge transfer material for green OLEDs. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Conformational studies of dammarane‐type triterpenoids using computational and NMR spectroscopic methods

Natural triterpenoids are of great interest to researchers of various fields as they possess diverse physicochemical and biological properties. In medicinal chemistry, detailed information about the chemical structures of bioactive triterpenoids often helps find new lead compounds. Herein, the low‐energy structures of (20S)‐protopanaxadiol and (20S)‐protopanaxatriol, the aglycones of various triterpenoid saponins found in Panax ginseng, and their (20R)‐epimers have been predicted by the geometry optimization of the conformers extracted from molecular dynamics simulations with the self‐consistent‐charge density functional tight‐binding method. By performing quantum mechanical calculations on the low‐energy conformers, we have estimated the NMR chemical shifts of the compounds, which display good agreement with the most recently reported experimental values within an expected range of errors. Our results indicate that theoretical estimation of the NMR parameters of a relatively large molecule with a molecular mass of 500 is feasible. Copyright © 2015 John Wiley & Sons, Ltd.     

Intramolecular charge transfer in π-conjugated oligomers: a theoretical study on the effect of temperature and oxidation state

Intramolecular charge transfer (ICT) of gaseous π-conjugated oligo-phenyleneethynylenes (OPE) induced by a homogeneous applied electric field has been theoretically investigated using a combined approach integrating molecular dynamics (MD) simulations and Perturbed Matrix Method calculations. In line with recent investigations, our results indicate the peculiar role of conformational transitions on OPE electronic properties which reflects on a strong temperature effect on ICT. Electron transfer reactions inducing chemical alteration on OPE, also taken into account in this study, revealed extremely important for explaining non-linear ICT effects and probably plays a central role in the mechanisms underlying molecular transport junctions. Our study further points out the necessity of using MD-based approach for modelling molecular electronics, even when relatively rigid molecular systems are concerned.     

Insights into the dynamics of evaporation and proton migration in protonated water clusters from Large‐scale Born–Oppenheimer direct dynamics

Large‐scale on‐the‐fly Born–Oppenheimer molecular dynamics simulations using recent advances in linear scaling electronic structure theory and trajectory integration techniques have been performed for protonated water clusters around the magic number (H₂O)_nH⁺, for n = 20 and 21. Besides demonstrating the feasibility and efficiency of the computational approach, the calculations reveal interesting dynamical details. Elimination of water molecules is found to be fast for both cluster sizes but rather insensitive to the initial geometry. The water molecules released acquire velocities compatible with thermal energies. The proton solvation shell changes between the well‐known Eigen and Zundel motifs and is characterized by specific low‐frequency vibrational modes, which have been quantified. The proton transfer mechanism largely resembles that of bulk water but one interesting variation was observed. © 2012 Wiley Periodicals, Inc.     

Monte Carlo morphological modelling of a P3HT:PCBM bulk heterojunction organic solar cell

To deepen the understanding of morphology evolution in bulk heterojunction P3HT:PCBM organic photovoltaics system by thermal treatment, domain‐size‐dependent interfacial energies were first determined by coarse‐grained molecular dynamics modelling and then used in Monte Carlo simulations of the morphology evolution. Thereby initial conditions associated with optimal interfacial surface area, continuous volume, as well as domain sizes, and spatial distributions of the phase separated domains were identified. In line with earlier studies, a 1:1 P3HT:PCBM blend ratio is found to exhibit the most efficient morphology for exciton dissociation and charge transport. Our simulations reveal that preseeding of P3HT crystal at the anode side prior to the annealing process will be instrumental to pin the formation of P3HT at the favorable electrode especially when seeding exceeds a threshold of 10% surface coverage, whereas denser seeding patterns beyond the threshold did not improve the active layer morphology further. The observed trilayer depth profile (in the absence of preseeded P3HT crystals) implies that the commonly used thickness 100 nm of the active layer is not ideal for ensuring that donor and acceptor phases dominate at opposite ends of the active layer. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53, 270–279     

Influence of thermal fluctuations on interfacial electron transfer in functionalized TiO2 semiconductors

The influence of thermal fluctuations on the dynamics of interfacial electron transfer in sensitized TiO2-anatase semiconductors is investigated by combining ab initio DFT molecular dynamics simulations and quantum dynamics propagation of transient electronic excitations. It is shown that thermal nuclear fluctuations speed up the underlying interfacial electron transfer dynamics by introducing nonadiabatic transitions between electron acceptor states, localized in the vicinity of the photoexcited adsorbate, and delocalized states extended throughout the semiconductor material, creating additional relaxation pathways for carrier diffusion. Furthermore, it is shown that room-temperature thermal fluctuations reduce the anisotropic character of charge diffusion along different directions in the anatase crystal and make similar the rates for electron injection from adsorbate states of different character. The reported results are particularly relevant to the understanding of temperature effects on surface charge separation mechanisms in molecular-based photo-optic devices.     

Computational study of structure,electronic, and microscopic charge transport properties of small conjugated diketopyrrolopyrrole‐thiophene molecules

π‐Conjugated small molecules containing diketopyrrolopyrrole (DPP) and thiophene moieties represent a modern class of functional materials that exhibit promising charge transport properties and therefore have great potential as building blocks of active elements of electronic devices. As a starting point of this computational study, the molecular structure, electronic characteristics, and reorganization energies associated with electron or hole transfer are considered. Prediction of molecular crystal packing is followed by the calculation of couplings between adjacent molecules and detection of the effective charge transfer pathways. Finally, the rates of charge transfer process are evaluated. The obtained results shed light not only on the properties of materials containing low‐molecular species but also serve as a benchmark for further classical force‐field simulations of DPP‐based polymers.     

Theoretical comparative studies on transport properties of pentacene,pentathienoacene, and 6,13‐dichloropentacene

Pentacene derivative 6,13‐dichloropentacene (DCP) is one of the latest additions to the family of organic semiconductors with a great potential for use in transistors. We carry out a detailed theoretical calculation for DCP, with systematical comparison to pentacene, pentathienoacene (PTA, the thiophene equivalent of pentacene), to gain insights in the theoretical design of organic transport materials. The charge transport parameters and carrier mobilities are investigated from the first‐principles calculations, based on the widely used Marcus electron transfer theory and quantum nuclear tunneling model, coupled with random walk simulation. Molecular structure and the crystal packing type are essential to understand the differences in their transport behaviors. With the effect of molecule modification, significant one‐dimensional π‐stacks are found within the molecular layer in PTA and DCP crystals. The charge transport along the a‐axis plays a dominant role for the carrier mobilities in the DCP crystal due to the strong transfer integrals within the a‐axis. Pentacene shows a relatively large 3D mobility. This is attributed to the relatively uniform electronic couplings, which thus provides more transport pathways. PTA has a much smaller 3D mobility than pentacene and DCP for the obvious increase of the reorganization energy with the introduction of thiophene. It is found that PTA and DCP exhibit lower HOMO (highest occupied molecular orbital) levels and better environmental stability, indicating the potential applications in organic electronics. © 2015 Wiley Periodicals, Inc.     

Correlation between thermal conductivity and bond length alternation in carbon nanotubes: A combined reverse nonequilibrium molecular dynamics—Crystal orbital analysis

The thermal conductivity (λ) of carbon nanotubes (CNTs) with chirality indices (5,0), (10,0), (5,5), and (10,10) has been studied by reverse nonequilibrium molecular dynamics (RNEMD) simulations as a function of different bond length alternation patterns (Δr_i). The Δr_i dependence of the bond force constant (k_rx) in the molecular dynamics force field has been modeled with the help of an electronic band structure approach. These calculations show that the Δr_i dependence of k_rx in tubes with not too small a diameter can be mapped by a simple linear bond length–bond order correlation. A bond length alternation with an overall reduction in the length of the nanotube causes an enhancement of λ, whereas an alternation scheme leading to an elongation of the tube is coupled to a decrease of the thermal conductivity. This effect is more pronounced in carbon nanotubes with larger diameters. The formation of a polyene‐like structure in the direction of the longitudinal axis has a negligible influence on λ. A comparative analysis of the RNEMD and crystal orbital results indicates that Δr_i‐dependent modifications of λ and the electrical conductivity are uncorrelated. This behavior is in‐line with a heat transfer that is not carried by electrons. Modifications of λ as a function of the bond alternation in the (10,10) nanotube are explained with the help of power spectra, which provide access to the density of vibrational states. We have suggested longitudinal low‐energy modes in the spectra that might be responsible for the Δr_i dependence of λ. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Comparative study of microscopic charge dynamics in crystalline acceptor-substituted oligothiophenes

By performing microscopic charge transport simulations for a set of crystalline dicyanovinyl-substituted oligothiophenes, we find that the internal acceptor-donor-acceptor molecular architecture combined with thermal fluctuations of dihedral angles results in large variations of local electric fields, substantial energetic disorder, and pronounced Poole-Frenkel behavior, which is unexpected for crystalline compounds. We show that the presence of static molecular dipoles causes large energetic disorder, which is mostly reduced not by compensation of dipole moments in a unit cell but by molecular polarizabilities. In addition, the presence of a well-defined π-stacking direction with strong electronic couplings and short intermolecular distances turns out to be disadvantageous for efficient charge transport since it inhibits other transport directions and is prone to charge trapping.     

Binding of Cationic and Neutral Phenanthridine Intercalators to a DNA Oligomer Is Controlled by Dispersion Energy: Quantum Chemical Calculations and Molecular Mechanics Simulations

Correlated ab initio as well as semiempirical quantum chemical calculations and molecular dynamics simulations were used to study the intercalation of cationic ethidium, cationic 5‐ethyl‐6‐phenylphenanthridinium and uncharged 3,8‐diamino‐6‐phenylphenanthridine to DNA. The stabilization energy of the cationic intercalators is considerably larger than that of the uncharged one. The dominant energy contribution with all intercalators is represented by dispersion energy. In the case of the cationic intercalators, the electrostatic and charge‐transfer terms are also important. The ΔG of ethidium intercalation to DNA was estimated at ?4.5 kcal mol^?1 and this value agrees well with the experimental result. Of six contributions to the final free energy, the interaction energy value is crucial. The intercalation process is governed by the non‐covalent stacking (including charge‐transfer) interaction while the hydrogen bonding between the ethidium amino groups and the DNA backbone is less important. This is confirmed by the evaluation of the interaction energy as well as by the calculation of the free energy change. The intercalation affects the macroscopic properties of DNA in terms of its flexibility. This explains the easier entry of another intercalator molecule in the vicinity of an existing intercalation site.     

Crystallographic and Dynamic Aspects of Solid‐State NMR Calibration Compounds: Towards ab Initio NMR Crystallography

The excellent results of dispersion‐corrected density functional theory (DFT‐D) calculations for static systems have been well established over the past decade. The introduction of dynamics into DFT‐D calculations is a target, especially for the field of molecular NMR crystallography. Four ¹³C ss‐NMR calibration compounds are investigated by single‐crystal X‐ray diffraction, molecular dynamics and DFT‐D calculations. The crystal structure of 3‐methylglutaric acid is reported. The rotator phases of adamantane and hexamethylbenzene at room temperature are successfully reproduced in the molecular dynamics simulations. The calculated ¹³C chemical shifts of these compounds are in excellent agreement with experiment, with a root‐mean‐square deviation of 2.0 ppm. It is confirmed that a combination of classical molecular dynamics and DFT‐D chemical shift calculation improves the accuracy of calculated chemical shifts.