In-Depth Ab Initio Study of Thermally Activated Delayed Fluorescence in 4,5-Di(9H-carbazol-9-yl)-phthalonitrile.

Molecules capable of thermally activated delayed fluorescence (TADF) are promising as emitters in organic light-emitting devices. Processes leading to and competing with TADF in 4,5-di(9H-carbazol-9-yl)-phthalonitrile are analyzed in detail. It is demonstrated that the key features of an efficient TADF emitter include the presence of two triplet states of different natures with potential energy surfaces crossing between the T1 and S1 minima and a noticeable dependence of the S1 → S0 oscillator strength on molecular deformations from low-frequency antisymmetric vibrational modes. These conclusions can be useful in the targeted design of efficient TADF emitters.

Thermodynamic stability of non-stoichiometric SrFeO3-δ: a hybrid DFT study.

SrFeO3-δ is a mixed ionic-electronic conductor with a complex magnetic structure that reveals a colossal magnetoresistance effect. This material and its solid solutions are attractive for various spintronic, catalytic and electrochemical applications, including cathodes for solid oxide fuel cells and permeation membranes. Its properties strongly depend on oxygen non-stoichiometry. Ab initio hybrid functional approach was applied herein to study the thermodynamic stability of a series of SrFeO3-δ compositions with several non-stoichiometries δ, ranging from 0 to 0.5 (SrFeO3-SrFeO2.875-SrFeO2.75-SrFeO2.5) as a function of temperature and oxygen pressure. The results obtained by two approaches, in which either (i) all electrons at Fe atoms explicitly described or (ii) inner core electrons at Fe atoms are replaced by effective core potential, are compared. Based on our calculations, phase diagrams were constructed, allowing the determination of environmental conditions for the existence of stable phases. It is shown that (within an employed model) only the SrFeO2.5 phase appears to be stable. The stability region for this phase was re-drawn on the contour map of oxygen chemical potential, presented as a function of temperature and oxygen partial pressure. A similar analysis was also performed using experimental Gibbs energies of perovskite formation from the elements. The present modelling strongly suggest a significant attraction between neutral oxygen vacancies. These vacancies are created during a series of the abovementioned SrFeO3-δ mutual transformations accompanied by oxygen release.

pyEFP: Automatic decomposition of the complex molecular systems into rigid polarizable fragments.

We present an open source tool able to describe intermolecular electrostatic interactions within the framework of the effective fragment potential (EFP) method. Complex molecular structure is subdivided into compact rigid fragments and parameters of their interactions are obtained from ab initio calculations. Automatic procedure allows for searching of these parameters into the existing database and merge new fragments into it. A set of standard fragments useful for the studies of organic semiconductors is also provided. Input files both for purely EFP and hybrid QM/MM calculations can be generated. The program is written in python and freely available on GitHub: https://github.com/ale-odinokov/pyEFP © 2017 Wiley Periodicals, Inc.

First principles crystal engineering of nonlinear optical materials. I. Prototypical case of urea.

The crystalline materials with nonlinear optical (NLO) properties are critically important for several technological applications, including nanophotonic and second harmonic generation devices. Urea is often considered to be a standard NLO material, due to the combination of non-centrosymmetric crystal packing and capacity for intramolecular charge transfer. Various approaches to crystal engineering of non-centrosymmetric molecular materials were reported in the literature. Here we propose using global lattice energy minimization to predict the crystal packing from the first principles. We developed a methodology that includes the following: (1) parameter derivation for polarizable force field AMOEBA; (2) local minimizations of crystal structures with these parameters, combined with the evolutionary algorithm for a global minimum search, implemented in program USPEX; (3) filtering out duplicate polymorphs produced; (4) reoptimization and final ranking based on density functional theory (DFT) with many-body dispersion (MBD) correction; and (5) prediction of the second-order susceptibility tensor by finite field approach. This methodology was applied to predict virtual urea polymorphs. After filtering based on packing similarity, only two distinct packing modes were predicted: one experimental and one hypothetical. DFT + MBD ranking established non-centrosymmetric crystal packing as the global minimum, in agreement with the experiment. Finite field approach was used to predict nonlinear susceptibility, and H-bonding was found to account for a 2.5-fold increase in molecular hyperpolarizability to the bulk value.

Thermodynamic stability of stoichiometric LaFeO3 and BiFeO3: a hybrid DFT study.

BiFeO3 perovskite attracts great attention due to its multiferroic properties and potential use as a parent material for Bi1-xSrxFeO3-δ and Bi1-xSrxFe1-yCoyO3-δ solid solutions in intermediate temperature cathodes of oxide fuel cells. Another iron-based LaFeO3 perovskite is the end member for well-known solid solutions (La1-xSrxFe1-yCoyO3-δ) used for oxide fuel cells and other electrochemical devices. In this study an ab initio hybrid functional approach was used for the study of the thermodynamic stability of both LaFeO3 and BiFeO3 with respect to decompositions to binary oxides and to elements, as a function of temperature and oxygen pressure. The localized (LCAO) basis sets describing the crystalline electron wave functions were carefully re-optimized within the CRYSTAL09 computer code. The results obtained by considering Fe as an all-electron atom and within the effective core potential technique are compared in detail. Based on our calculations, the phase diagrams were constructed allowing us to predict the stability region of stoichiometric materials in terms of atomic chemical potentials. This permits determining the environmental conditions for the existence of stable BiFeO3 and LaFeO3. These conditions were presented as contour maps of oxygen atoms' chemical potential as a function of temperature and partial pressure of oxygen gas. A similar analysis was also performed using the experimental Gibbs energies of formation. The obtained phase diagrams and contour maps are compared with the calculated ones.

Synthesis and photophysical properties of halogenated derivatives of (dibenzoylmethanato)boron difluoride.

A series of (dibenzoylmethanato)boron difluoride (BF2DBM) derivatives with a halogen atom in one of the phenyl rings at the para-position were synthesized and used to elucidate the effects of changing the attached halogen atom on the photophysical properties of BF2DBM. The room-temperature absorption and fluorescence maxima of fluoro-, chloro-, bromo- and iodo-substituted derivatives of BF2DBM in THF are red-shifted by about 2-10nm relative to the corresponding peaks of the parent BF2DBM. The fluorescence quantum yields of the halogenated BF2DBMs (except the iodinated derivative) are larger than that of the unsubstituted BF2DBM. All the synthesized compounds are able to form fluorescent exciplexes with benzene and toluene (emission maxima at λem=433 and 445nm, respectively). The conformational structure and electronic spectral properties of halogenated BF2DBMs have been modeled by DFT/TDDFT calculations at the PBE0/SVP level of theory. The structure and fluorescence spectra of exciplexes were calculated using the CIS method with empirical dispersion correction.

Electronic Structure and Energy Transfer in Europium(III)-Ciprofloxacin Complexes: A Theoretical Study.

The structure and ligand-localized excited states of [Eu(cfqH) (cfq)(H2O)4]Cl2 (cfqH is ciprofloxacin) are studied by XMCQDPT2/CASSCF with full geometry optimization. The complex includes one anionic and one zwitterionic ligand. Two low-lying triplet states, both localized on the anionic ligand, are found. One of them has sufficient energy to transfer to the (5)D1 sublevel of Eu(3+), because its T-S0 vertical transition energy is equal (or very close) to the (7)F0-(5)D1 Eu(3+) excitation energy. The other triplet state has a very small S0-T1 gap, which favors fast nonradiative relaxation. Two other triplet states are localized on the zwitterionic ligand. One low-lying excited singlet state (S1) is localized on the anionic ligand; the other excited singlet is localized on the zwitterionic one. Spin-orbit coupling constants were calculated for the relaxed geometry of each state (ground state, two low-lying triplets, and one low-lying excited singlet) by spin-orbit configuration interaction (CI) with Pauli-Breit Hamiltonian. Large spin-orbit coupling constants between S1 and both triplets together with small energy gaps are indicative of fast intersystem crossing (ISC) from the excited singlet state to the triplet manifold. This ISC process is followed by energy transfer from the ligand-localized triplet states to the (5)D1 sublevel of Eu(3+). However, relatively large spin-orbit coupling constants between S0 and one of the triplet states together with the small T-S0 energy gap shows that this state can decay without transferring its energy to Eu(3+). This mechanism is expected to be common for other Ln(3+)-fluoroquinolone complexes.

Vibronic bandshape of the absorption spectra of dibenzoylmethanatoboron difluoride derivatives: analysis based on ab initio calculations.

The nature of absorption bandshapes of dibenzoylmethanatoboron difluoride (DBMBF2) dye substituted in ortho-, meta-, and para-positions of the phenyl ring is investigated using DFT and TDDFT with the range-separated hybrid CAM-B3LYP functional and the 6-311G(d,p) basis set. The solvent effects are taken into account within the polarized continuum model. The vibronic bandshape is simulated using a time-dependent linear coupling model with a vertical gradient approach through an original code. For flexible chromophores, the spectra of individual conformers are summed up with Boltzmann factors. It is shown that the long-wavelength absorption bandshape of DBMBF2 derivatives is determined by three factors: the relative statistical weights of conformers with different electronic absorption patterns, the relative position and intensity of the second low-energy electronic transition, and the vibronic structure of individual electronic peaks. The latter is governed by the relationship between the hard vibrational modes, which contribute to vibronic progression, and soft modes, which provide broadening of the peaks. The simulated spectra of the dyes in the study are generally consistent with the available experimental data and explain the observed spectral features.

Specific Interactions of Neutral Side Chains of an Adsorbed Protein with the Surface of α-Quartz and Silica Gel.

Many key features of the protein adsorption on the silica surfaces still remain unraveled. One of the open questions is the interaction of nonpolar side chains with siloxane cavities. Here, we use nonequilibrium molecular dynamics simulations for the detailed investigation of the binding of several hydrophobic and amphiphilic protein side chains with silica surface. These interactions were found to be a possible driving force for protein adsorption. The free energy gain was larger for the disordered surface of amorphous silica gel as compared to α-quartz, but the impact depended on the type of amino acid. The dependence was analyzed from the structural point of view. For every amino acid an enthalpy-entropy compensation behavior was observed. These results confirm a hypothesis of an essential role of hydrophobic interactions in protein unfolding and irreversible adsorption on the silica surface.

Symmetry-Breaking in Cationic Polymethine Dyes: Part 2. Shape of Electronic Absorption Bands Explained by the Thermal Fluctuations of the Solvent Reaction Field.

The electronic absorption spectra of the symmetric cyanines exhibit dramatic dependence on the conjugated chain length: whereas short-chain homologues are characterized by the narrow and sharp absorption bands of high intensity, the long-chain homologues demonstrate very broad, structureless bands of low intensity. Spectra of the intermediate homologues combine both features. These broad bands are often explained using spontaneous symmetry-breaking and charge localization at one of the termini, and the combination of broad and sharp features was interpreted as coexistence of symmetric and asymmetric species in solution. These explanations were not supported by the first principle simulations until now. Here, we employ a combination of time-dependent density functional theory, a polarizable continuum model, and Franck-Condon (FC) approximation to predict the absorption line shapes for the series of 2-azaazulene and 1-methylpyridine-4-substituted polymethine dyes. To simulate inhomogeneous broadening by the solvent, the molecular structures are optimized in the presence of a finite electric field of various strengths. The calculated FC line shapes, averaged with the Boltzmann weights of different field strengths, reproduce the experimentally observed spectra closely. Although the polarizable continuum model accounts for the equilibrium solvent reaction field at absolute zero, the finite field accounts for the thermal fluctuations in the solvent, which break the symmetry of the solute molecule. This model of inhomogeneous broadening opens the possibility for computational studies of thermochromism. The choice of the global hybrid exchange-correlation functional SOGGA11-X, including 40% of the exact exchange, plays the critical role in the success of our model.

Ab initio study of energy transfer pathways in dinuclear lanthanide complex of europium(III) and terbium(III) ions.

An ab initio XMCQDPT2/CASSCF study of energy transfer processes in the dinuclear lanthanide complex [(Acac)3Eu(µ-Bpym)Tb(Acac)3] (Acac is acetylacetonate, and Bpym is 2,2'-bipyrimidine) and a corresponding computational procedure are presented. Because ligands in lanthanide complexes weakly interact with each other, the large dinuclear complex bearing seven organic ligands is divided into fragments that reproduce the electrostatic effects of the ions on the electronic and geometrical structure of the ligands. The multireference XMCQDPT2/CASSCF approach is directly applied to these relatively small fragments with reasonable computational cost. The calculated energies of the singlet and triplet excited states agree well with the experiment. Based on the calculated energies, the energy level diagrams of the complex are constructed and intramolecular energy transfer channels are determined.

Charge localization and charge transfer in the Bebq2 monomer and dimer.

The geometrical structure and electronic properties of bis(10-hydroxybenzo[h]quinolinato)beryllium (Bebq2) molecule and its dimer both in the neutral and in the positively and negatively charged states were studied using quantum-chemical calculations. It is found that the excess charge in the charged systems is localized on one of the hydroxybenzoquinoline ligands. Structural changes in charged Bebq2 are pronounced in the charged ligand and nearly negligible in the neutral ligand. Charge transfer from the charged ligand to a neutral one can proceed either within a single Bebq2 monomer molecule or between the different monomers in the Bebq2 dimer. The corresponding hopping integrals were estimated as half the excitation energy from the ground to the first excited state of either the monomer or the dimer calculated at the avoided crossing point.

Ab initio study of phosphorescent emitters based on rare-earth complexes with organic ligands for organic electroluminescent devices.

An ab initio approach is developed for calculation of low-lying excited states in Ln(3+) complexes with organic ligands. The energies of the ground and excited states are calculated using the XMCQDPT2/CASSCF approximation; the 4f electrons of the Ln(3+) ion are included in the core, and the effects of the core electrons are described by scalar quasirelativistic 4f-in-core pseudopotentials. The geometries of the complexes in the ground and triplet excited states are fully optimized at the CASSCF level, and the resulting excited states have been found to be localized on one of the ligands. The efficiency of ligand-to-lanthanide energy transfer is assessed based on the relative energies of the triplet excited states localized on the organic ligands with respect to the receiving and emitting levels of the Ln(3+) ion. It is shown that ligand relaxation in the excited state should be properly taken into account in order to adequately describe energy transfer in the complexes. It is demonstrated that the efficiency of antenna ligands for lanthanide complexes used as phosphorescent emitters in organic light-emitting devices can be reasonably predicted using the procedure suggested in this work. Hence, the best antenna ligands can be selected in silico based on theoretical calculations of ligand-localized excited energy levels.

Atomistic simulations of materials for optical chemical sensors: DFT-D calculations of molecular interactions between gas-phase analyte molecules and simple substrate models.

The structures of complexes of some small molecules (formaldehyde, acetaldehyde, ammonia, methylamine, methanol, ethanol, acetone, benzene, acetonitrile, ethyl acetate, chloroform, and tetrahydrofuran, considered as possible analytes) with ethylbenzene and silanol (C(6)H(5)C(2)H(5) and SiH(3)OH, considered as models of polystyrene and silica gel substrates) and with acridine (C(13)H(9)N, considered as a model of an indicator dye molecule of the acridine series) and the corresponding interaction energies have been calculated using the DFT-D approximation. The PBE exchange-correlation potential was used in the calculations. The structures of complexes between the analyte and the substrate were determined by optimizing their ground-state geometry using the SVP split-valence double-zeta plus polarization basis set. The complex formation energies were refined by single-point calculations at the calculated equilibrium geometries using the sufficiently large triple-zeta TZVPP basis set. The calculated interaction energies are used to assess the possibility of using dyes of the acridine series adsorbed on a polystyrene or silica substrate for detecting the small molecules listed above.