A practical method for the use of curvilinear coordinates in calculations of normal-mode-projected displacements and Duschinsky rotation matrices for large molecules

Excimer-induced high-efficiency fluorescence due to pairwise anthracene stacking in a crystal with long lifetime

Herein, we report an anthracene-based material, 2-(anthracen-9-yl)thianthrene (), whose crystal exhibits excimer fluorescence with an unexpected high luminous efficiency (up to 80%) and long lifetime (163.75 ns), due to pairwise anthracene stacking. These results will update the traditional view that excimers are poorly efficient in photoluminescence.

The effect of the embedded o-carborane ligand on the photophysical properties of a cyclometalated Pt(ii) complex: a theoretical investigation

The embedded o-carborane ligands will have an impact on the photophysical properties and quantum efficiency of these highly efficient blue phosphorescent materials.

New insight into the effects of N^N ligand isomerization and methyl modification on the phosphorescence properties of Cu(i) complexes with (1-(2-pyridyl)pyrazole/imidazole) ligands

The influence of ligand isomerization and methylation modification on the phosphorescence properties of metal complexes has been attracting great interest owing to the significant improvement of phosphorescence efficiency.

Theoretical study of radiative and nonradiative decay rates for Cu(i) complexes with double heteroleptic ligands

Phosphorescence quantum yield increases significantly because k_r increases and k_nr decreases upon introducing electron-donating groups to N^N ligands.

Enhancing charge mobilities in organic semiconductors by selective fluorination: a design approach based on a quantum mechanical perspective

Selective fluorination of organic semiconducting molecules is proposed as a means to achieving enhanced hole mobility. Naphthalene is examined here as a root molecular system with fluorination performed at various sites. Our quantum chemical calculations show that selective fluorination can enhance attractive intermolecular interactions while reducing charge trapping. Those observations suggest a design principle whereby fluorination is utilized for achieving high charge mobilities in the crystalline form. The utility of this design principle is demonstrated through an application to perylene, which is an important building block of organic semiconducting materials. We also show that a quantum mechanical perspective of nuclear degrees of freedom is crucial for a reliable description of charge transport.

Theoretical study on electron structure and charge transport properties of tetraazapentacene derivatives

By Means of Marcus electron transfer theory, the charge transport properties of tetraazapentacene (4N-PEN) derivatives were systematically explored. The reorganization energies were studied by both adiabatic potential-energy surfaces and normal mode analysis. The charge diffusion constants were evaluated from the random walk simulation. From the perspective of homology modeling, a selected 4N-PEN derivative without experimental crystal structure was built into three kinds of possible packing modes with reference to its relative analogues and then fully optimized. The calculated results show that the charge transport property for the same kind of systems strongly depends on the packing mode, and the π···stacking is more beneficial for electron transport of 4N-PEN derivatives. Meanwhile, the 4N-PEN derivatives have larger electron transfer integrals and lower energy levels of the lowest unoccupied molecular orbitals as well as smaller electron reorganization energies, which provides a three-in-one advantage for electron transport. Fascinatingly, the data obtained from the hopping and band models both suggest that the 4N-PEN derivatives have the intrinsic property of electron transport. Thus, the 4N-PEN derivatives have the potential for preparing n-type organic semiconductors.

Charge Transport in Molecular Materials: An Assessment of Computational Methods

The booming field of molecular electronics has fostered a surge of computational research on electronic properties of organic molecular solids. In particular, with respect to a microscopic understanding of transport and loss mechanisms, theoretical studies assume an ever-increasing role. Owing to the tremendous diversity of organic molecular materials, a great number of computational methods have been put forward to suit every possible charge transport regime, material, and need for accuracy. With this review article we aim at providing a compendium of the available methods, their theoretical foundations, and their ranges of validity. We illustrate these through applications found in the literature. The focus is on methods available for organic molecular crystals, but mention is made wherever techniques are suitable for use in other related materials such as disordered or polymeric systems.

Simulation of Vibronic Spectra of Flexible Systems: Hybrid DVR-Harmonic Approaches

Our general framework for the simulation of vibrational signatures in electronic spectra has been extended to treat one large-amplitude motion (LAM) at the anharmonic level, coupled to the other small-amplitude motions (SAM) treated as harmonic. The coupling between LAM and SAM is minimized thanks to the use of delocalized internal coordinates, which are built automatically from the molecular topology. General LAMs can be employed, ranging from intrinsic reaction coordinates to rigid or flexible paths based on the distinguished coordinate approach. The anharmonic model is based on a fully numerical method based on the discrete variable representation (DVR) theory, supporting different types of boundary conditions. The inclusion of this model in a general-purpose electronic structure code makes available to the user a large panel of quantum chemistry models, for both isolated and condensed phases. The flexibility and reliability of the new framework are illustrated by some case studies, covering various types of LAMs, ranging from a small test case, the photoelectron spectrum of ammonia, to larger systems, such as phenylanthracene and cyclobutanone.

A DFT Study on the Electronic Structures and Conducting Properties of Rubrene and its Derivatives in Organic Field-Effect Transistors

We systematically studied the electronic structures and conducting properties of rubrene and its derivatives reported recently, and disscussed the influences of electron-withdrawing groups and chemical oxidation on the reorganization energies, crystal packing, electronic couplings, and charge injection barrier of rubrene. Hirshfeld surface analysis and quantum-chemical calculations revealed that the introduction of CF₃ groups into rubrene decreases the H···H repulsive interaction and increases intermolecular F···H/H···F attractive interactions, which resulted in the tight packing arrangement and the increase of the electronic couplings, and finally cause the higer intrinsic hole-mobility in bis(trifluoromethyl)-dimethyl-rubrene crystal (μ_h = 19.2 cm² V⁻¹ s⁻¹) than in rubrene crystal (μ_h = 15.8 cm² V⁻¹ s⁻¹). In comparison, chemical oxidation reduces charge-carrier mobility of rubrene crystal by 2~4 orders of magnitude and increased the hole and electron injection barrier, which partly explains the rubrene-based field-effect transistor performance degrades upon exposure to air. Furthermore, we also discussed the influence of structural parameters of carbon nanotube (CNT) electrode on charge injection process, which suggests that the regulation of CNT diameters and increasing in thickness is an effective strategy to optimize CNT work functions and improve n-type OFET performances based on these organic materials.

Molecular Crystal Engineering: Tuning Organic Semiconductor from p-type to n-type by Adjusting Their Substitutional Symmetry

Focusing on the bottleneck of molecularly engineered organic semiconductors, a breakthrough is made to tune the electronic properties of organic semiconductors from p-type to n-type by using fluorinated metal phthalocyanines as examples. The experimentally observed p-type to n-type transition characteristics of single-crystal field-effect devices result from a combination of extrinsic and intrinsic properties of materials with different fluoridation substitution.

Applying strong external electric field to thiophene-based oligomers: A promising approach to upgrade semiconducting performance

A key parameter dictating the rate of charge transfer (CT) is reorganization energy (λ), an energy associated with geometry changes during hole/electron transfer. We show that "ironing" the inter-ring dihedral angles of oligothiophenes via proper substitutions or insertions (e.g., -OR, -F or -C≡C-), decreases the λ and thus promotes CT according to Marcus equation. Our results demonstrate, to attain a smaller λ, extending oligomer length is only significant if the flattened backbone structure is realized. Of great interest is that external electric fields, which are ubiquitous in electronic devices yet commonly overlooked in the computation of λ, can have a significantly greater impact than conventional substitutions. It is important to emphasize, the responses of λ to external fields is system-dependent. Compared to fused-ring conjugated systems, single-bond connected thiophenes are more sensitive to external fields. F_x lowers the λ (552 meV) of quaterthiophene by almost 80% at the intensity of 1 V/Å, down to a value (125 meV) which is even lower than that of pentacene (154 meV) and rubrene (219 meV) at the same level of theory. © 2016 Wiley Periodicals, Inc.

Theoretical perspective of the excited state intramolecular proton transfer for a compound with aggregation induced emission in the solid phase

The ESIPT process and AIE mechanism are theoretically investigated by a QM/MM method.

Understanding the interplay between the solvent and nuclear rearrangements in the negative solvatochromism of a push–pull flexible quinolinium cation

The main effects (solvation, vibronic progression) affecting the band position and shape of a push–pull flexible quinolinium cation OPA are highlighted.

Excited state dynamics for hybridized local and charge transfer state fluorescent emitters with aggregation-induced emission in the solid phase: a QM/MM study

Investigation on the excited state dynamics to reveal the AIE and HLCT mechanisms by a QM/MM method.

A theoretical study on the electronic properties of two ring-fused derivatives of 9,10-diphenylanthracene

Two new contorted polycyclic aromatic hydrocarbons (PAHs) 1 and 2 have been synthesized by Perepichka and coworkers (Org. Lett., 2015, 17, 4224).

Charge transport in organic donor–acceptor mixed-stack crystals: the role of nonlocal electron–phonon couplings

The influence of nonlocal electron–phonon couplings on charge transport is found to be very small in organic donor–acceptor mixed-stack crystals.

Coherent intermolecular proton transfer in the acid–base reaction of excited state pyranine

The acidic proton in pyranine is transferred coherently to acetate through the stretching motion of the whole molecule.

Revisiting Vertical Models To Simulate the Line Shape of Electronic Spectra Adopting Cartesian and Internal Coordinates

Vertical models for the simulation of spectroscopic line shapes expand the potential energy surface (PES) of the final state around the equilibrium geometry of the initial state. These models provide, in principle, a better approximation of the region of the band maximum. At variance, adiabatic models expand each PES around its own minimum. In the harmonic approximation, when the minimum energy structures of the two electronic states are connected by large structural displacements, adiabatic models can breakdown and are outperformed by vertical models. However, the practical application of vertical models faces the issues related to the necessity to perform a frequency analysis at a nonstationary point. In this contribution we revisit vertical models in harmonic approximation adopting both Cartesian (x) and valence internal curvilinear coordinates (s). We show that when x coordinates are used, the vibrational analysis at nonstationary points leads to a deficient description of low-frequency modes, for which spurious imaginary frequencies may even appear. This issue is solved when s coordinates are adopted. It is however necessary to account for the second derivative of s with respect to x, which here we compute analytically. We compare the performance of the vertical model in the s-frame with respect to adiabatic models and previously proposed vertical models in x- or Q₁-frame, where Q₁ are the normal coordinates of the initial state computed as combination of Cartesian coordinates. We show that for rigid molecules the vertical approach in the s-frame provides a description of the final state very close to the adiabatic picture. For sizable displacements it is a solid alternative to adiabatic models, and it is not affected by the issues of vertical models in x- and Q₁-frames, which mainly arise when temperature effects are included. In principle the G matrix depends on s, and this creates nonorthogonality problems of the Duschinsky matrix connecting the normal modes of initial and final states in adiabatic approaches. We highlight that such a dependence of G on s is also an issue in vertical models, due to the necessity to approximate the kinetic term in the Hamiltonian when setting up the so-called GF problem. When large structural differences exist between the initial and the final-state minima, the changes in the G matrix can become too large to be disregarded.

Effect of five-membered ring and heteroatom substitution on charge transport properties of perylene discotic derivatives: A theoretical approach

Density functional theory calculations were carried out to investigate the evolvement of charge transport properties of a set of new discotic systems as a function of ring and heteroatom (B, Si, S, and Se) substitution on the basic structure of perylene. The replacement of six-membered rings by five-membered rings in the reference compound has shown a prominent effect on the electron reorganization energy that decreases ∼0.2 eV from perylene to the new carbon five-membered ring derivative. Heteroatom substitution with boron also revealed to lower the LUMO energy level and increase the electron affinity, therefore lowering the electron injection barrier compared to perylene. Since the rate of the charge transfer between two molecules in columnar discotic systems is strongly dependent on the orientation of the stacked cores, the total energy and transfer integral of a dimer as a disc is rotated with respect to the other along the stacking axis have been predicted. Aimed at obtaining a more realistic approach to the bulk structure, the molecular geometry of clusters made up of five discs was fully optimized, and charge transfer rate and mobilities were estimated for charge transport along a one dimensional pathway. Heteroatom substitution with selenium yields electron transfer integral values ∼0.3 eV with a relative disc orientation of 25°, which is the preferred angle according to the dimer energy profile. All the results indicate that the tetraselenium-substituted derivative, not synthetized so far, could be a promising candidate among those studied in this work for the fabrication of n-type semiconductors based on columnar discotic liquid crystals materials.

General formulation of vibronic spectroscopy in internal coordinates

Our general platform integrating time-independent and time-dependent evaluations of vibronic effects at the harmonic level for different kinds of absorption and emission one-photon, conventional and chiral spectroscopies has been extended to support various sets of internal coordinates. Thanks to the implementation of analytical first and second derivatives of different internal coordinates with respect to cartesian ones, both vertical and adiabatic models are available, with the inclusion of mode mixing and, possibly, Herzberg-Teller contributions. Furthermore, all supported non-redundant sets of coordinates are built from a fully automatized algorithm using only a primitive redundant set derived from a bond order-based molecular topology. Together with conventional stretching, bending, and torsion coordinates, the availability of additional coordinates (including linear and out-of-plane bendings) allows a proper treatment of specific systems, including, for instance, inter-molecular hydrogen bridges. A number of case studies are analysed, showing that cartesian and internal coordinates are nearly equivalent for semi-rigid systems not experiencing significant geometry distortions between initial and final electronic states. At variance, delocalized (possibly weighted) internal coordinates become much more effective than their cartesian counterparts for flexible systems and/or in the presence of significant geometry distortions accompanying electronic transitions.

Using the isotope effect to probe an aggregation induced emission mechanism: theoretical prediction and experimental validation

Aggregation-induced emission (AIE) has become a hot topic for a variety of potential applications, but the understanding of its working mechanism is still under scrutiny. Herein, we proposed the use of the isotope effect (IE) to identify the AIE mechanism: under the restriction of an internal motion mechanism, the IE is pronouncedly different in excited-state decay rates when contrasting AIE luminogens (AIEgens) and non-AIEgens in theoretical calculations. For the complete deuteration of AIEgens, the IE of nonradiative decay rate in solution (<-10%) is much weaker than that (-65% to -95%) in aggregate, because the former stems from the overall results of competitive vibronic coupling and the severe mixing of low-frequency modes while the latter mainly comes from the vibronic coupling only. The experimental results confirm the isotopic "jump" behaviors in AIEgens well. However, non-AIEgens exhibit equivalent IEs (-40% to -90%) in both solution and solid phases. Further partial deuteration schemes for the 6-ring AIE analogues show positional dependence.

Simulation of Vacuum UV Absorption and Electronic Circular Dichroism Spectra of Methyl Oxirane: The Role of Vibrational Effects

Vibrationally resolved one-photon absorption and electronic circular dichroism spectra of (R)-methyl oxirane were calculated with different electronic and vibronic models selecting, through an analysis of the convergence of the results, the best compromise between reliability and computational cost. Linear-response TD-DFT/CAM-B3LYP/SNST electronic computations in conjunction with the simple vertical gradient vibronic model were chosen and employed for systematic comparison with the available experimental data. Remarkable agreement between simulated and experimental spectra was achieved for both one-photon absorption and circular dichroism concerning peak positions, relative intensities, and general spectral shapes considering the computational efficiency of the chosen theoretical approach. The significant improvement of the results with respect to smearing of vertical electronic transitions by phenomenological Gaussian functions and the possible inclusion of solvent effects by polarizable continuum models at a negligible additional cost paves the route toward the simulation and analysis of spectral shapes of complex molecular systems in their natural environment.

"Doping" pentacene with sp(2)-phosphorus atoms: towards high performance ambipolar semiconductors

Recent research progress in black phosphorus sheets strongly encourages us to employ pentacene as a parent system to systematically investigate how the "doping" of sp(2)-phosphorus atoms onto the backbone of pentacene influences its optical and charge transport properties. Our theoretical investigations proved that increasing the contribution of the pz atomic orbital of the sp(2)-phosphorus to the frontier molecular orbital of phosphapentacenes could significantly decrease both hole and electron reorganization energies and dramatically red-shift the absorption of pentacene. The record smallest hole and electron reorganization energies of 69.80 and 95.74 meV for heteropentacene derivatives were obtained. These results suggest that phosphapentacenes (or phosphaacenes) could be potential promising candidates to achieve both higher and balanced mobilities in organic field effect transistors and realize a better power conversion efficiency in organic photovoltaics.

Nuclear quantum tunnelling and carrier delocalization effects to bridge the gap between hopping and bandlike behaviors in organic semiconductors

The experimental carrier mobility value of organic semiconductors has been increasing rapidly in recent years to well exceed the theoretical limit based on the hopping model calculated using the semi-classical Marcus theory, calling for better understanding and evaluation of carrier mobility. On the other hand, bandlike transport behavior has been observed for some ultra-pure and closely-packed organic single crystals. In this work, we identify the roles of quantum nuclear tunnelling and the charge delocalization effects, leading to a comprehensive computational approach to assess the carrier mobility of organic semiconductors. We present the first-principles evaluated mobility results for some representative organic transport materials at four levels ranging from semiclassical hopping to quantum nuclear enabled hopping and to quantum wavepacket diffusion, and eventually to complete bandlike descriptions. We provide a comprehensive tool to assess the carrier mobility in organic semiconductors based on such improved understanding.

Theoretical investigation of the broad one-photon absorption line-shape of a flexible symmetric carbazole derivative

The one-photon absorption spectrum of a carbazole derivative has been studied by employing density functional response theory combined with a mixed quantum/classical approach to simulate the spectral shape.

Controlling electronic effects and intermolecular packing in contorted polyaromatic hydrocarbons (c-PAHs): towards high mobility field effect transistors

Classical MD simulations followed by DFT calculations for computationally screened contorted polyaromatic hydrocarbons reveal encouraging OFET potential.

Hydrostatic pressure effect on charge transport properties of phenacene organic semiconductors

A detailed DFT study on the effect of applied pressure on the hole and electron mobility of phenacene organic semiconductors using Marcus classical charge transfer theory.

Magnesium cofactor produces unpaired electrons confined by triplet nucleotide loops in a full-turn DNA fragment

The B-DNA curvature arising from pairing between nucleotides of the two curved complementary DNA strands affects the oxidation number of magnesium cofactor.