Quantum-Quantum and Quantum-Quantum-Classical Schemes for Near-Gap Excitations with Projection-Based-Embedded -Bethe-Salpeter Equation.

We present quantum-quantum and quantum-quantum-classical schemes based on many-body Green's functions theory in the approximation with the Bethe-Salpeter equation (-BSE) employing projection-based-embedding (PbE). Such approaches allow defining active and inactive subsystems of larger, complex molecular systems, with only the smaller active subsystem being explicitly treated by -BSE offering significant computational advantages. However, as PbE can modify the single-particle states in the Kohn-Sham (KS) ground state calculation and screening effects from the inactive region are not automatically included in -BSE, results from such PbE--BSE calculations can deviate from a full-system reference.

Electronic Couplings and Conversion Dynamics between Localized and Charge Transfer Excitations from Many-Body Green's Functions Theory.

We investigate the determination of electronic coupling between localized excitations (LEs) and charge-transfer (CT) excitations based on many-body Green's functions theory in the approximation with the Bethe-Salpeter equation (-BSE). Using a small molecule dimer system, we first study the influence of different diabatization methods, as well as different model choices within -BSE, such as the self-energy models or different levels of self-consistency, and find that these choices affect the LE-CT couplings only minimally. We then consider a large-scale low-donor morphology formed from rubrene and fullerene and evaluate the LE-CT couplings based on coupled -BSE-molecular mechanics calculations.

Effect of Solvent Removal Rate and Annealing on the Interface Properties in a Blend of a Diketopyrrolopyrrole-Based Polymer with Fullerene.

We study the effect of solvent-free annealing and explicit solvent evaporation protocols in classical molecular dynamics simulations on the interface properties of a blend of a diketopyrrolopyrrole (DPP) polymer with conjugated substituents (DPP2Py2T) and PCBM[60]. We specifically analyze the intramolecular segmental mobility of the different polymer building blocks as well as intermolecular radial and angular distribution functions between donor and acceptor. The annealing simulations reveal an increase of the glass-transition temperature of 45 K in the polymer-fullerene blend compared to that of pure DPP2Py2T.

Machine Learning of Quasiparticle Energies in Molecules and Clusters.

We present a Δ-machine learning approach for the prediction of GW quasiparticle energies (ΔMLQP) and photoelectron spectra of molecules and clusters, using orbital-sensitive representations (OSRs) based on molecular Cartesian coordinates in kernel ridge regression-based supervised learning. Coulomb matrix, bag-of-bond, and bond-angle-torsion representations are made orbital-sensitive by augmenting them with atom-centered orbital charges and Kohn-Sham orbital energies, both of which are readily available from baseline calculations at the level of density functional theory (DFT). We first illustrate the effects of different constructions of the OSRs on the prediction of frontier orbital energies of 22k molecules of the QM8 data set and show that it is possible to predict the full photoelectron spectrum of molecules within the data set using a single model with a mean absolute error below 0.

Excited-State Geometry Optimization of Small Molecules with Many-Body Green's Functions Theory.

We present a benchmark study of gas phase geometry optimizations in the excited states of carbon monoxide, acetone, acrolein, and methylenecyclopropene using many-body Green's functions theory within the approximation and the Bethe-Salpeter equation (BSE) employing numerical gradients. We scrutinize the influence of several typical approximations in the -BSE framework; we used one-shot or eigenvalue self-consistent ev, employing a fully analytic approach or plasmon-pole model for the frequency dependence of the electron self-energy, or performing the BSE step within the Tamm-Dancoff approximation. The obtained geometries are compared to reference results from multireference perturbation theory (CASPT2), variational Monte Carlo (VMC) method, second-order approximate coupled cluster (CC2) method, and time-dependent density-functional theory (TDDFT).

Development and Testing of an All-Atom Force Field for Diketopyrrolopyrrole Polymers with Conjugated Substituents.

We develop an all-atom force field for a series of diketopyrrolopyrrole polymers with two aromatic pyridine substituents and a variable number of π-conjugated thiophene units in the backbone (DPP2PyT), used as donor materials in organic photovoltaic devices. Available intrafragment parameterizations of the individual fragment building blocks are combined with interfragment bonded and nonbonded parameters explicitly derived from density functional theory calculations. To validate the force field, we perform classical molecular dynamics simulations of single polymer chains with = 1, 2, 3 in good and bad solvents and of melts.

Excited-state electronic structure of molecules using many-body Green's functions: Quasiparticles and electron-hole excitations with VOTCA-XTP.

We present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green's function theory in the GW approximation with the Bethe-Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA-XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations.

Ultrafast Formation of the Charge Transfer State of Prodan Reveals Unique Aspects of the Chromophore Environment.

Lipophilic dyes such as laurdan and prodan are widely used in membrane biology due to a strong bathochromic shift in emission that reports the structural parameters of the membrane such as area per molecule. Disentangling of the factors which control the spectral shift is complicated by the stabilization of a charge-transfer-like excitation of the dye in polar environments. Predicting the emission therefore requires modeling both the relaxation of the environment and the corresponding evolution of the excited state.

Backbone Chemical Composition and Monomer Sequence Effects on Phenylene Polymer Persistence Lengths.

Despite a vast body of the literature devoted to the use of phenylene polymers in the fabrication of graphene nanoribbons, the study of the physical properties of these precursors still poses open questions whose answers will certainly contribute to the design of more efficient/precise synthesis protocols. Particularly, persistence length measurements combined with size exclusion chromatography techniques assign both semiflexible to semirigid structures depending on the molecular weight of the precursor (NaritaNat. Chem.

Insights into the Kinetics of Supramolecular Comonomer Incorporation in Water.

Multicomponent supramolecular polymers are a versatile platform to prepare functional architectures, but a few studies have been devoted to investigate their noncovalent synthesis. Here, we study supramolecular copolymerizations by examining the mechanism and time scales associated with the incorporation of new monomers in benzene-1,3,5-tricarboxamide (BTA)-based supramolecular polymers. The BTA molecules in this study all contain three tetra(ethylene glycol) chains at the periphery for water solubility but differ in their alkyl chains that feature either 10, 12 or 13 methylene units.

Evolutionary Approach to Constructing a Deep Feedforward Neural Network for Prediction of Electronic Coupling Elements in Molecular Materials.

We present a general framework for the construction of a deep feedforward neural network (FFNN) to predict distance and orientation dependent electronic coupling elements in disordered molecular materials. An evolutionary algorithm automatizes the selection of an optimal architecture of the artificial neural network within a predefined search space. Systematic guidance, beyond minimizing the model error with stochastic gradient descent based backpropagation, is provided by simultaneous maximization of a model fitness that takes into account additional physical properties, such as the field-dependent carrier mobility.

Improved general-purpose five-point model for water: TIP5P/2018.

A new five point potential for liquid water, TIP5P/2018, is presented along with the techniques used to derive its charges from per-molecule electrostatic potentials in the liquid phase using the split charge equilibration of Nistor [J. Chem. Phys.

Electronic Excitations in Complex Molecular Environments: Many-Body Green's Functions Theory in VOTCA-XTP.

Many-body Green's functions theory within the GW approximation and the Bethe-Salpeter Equation (BSE) is implemented in the open-source VOTCA-XTP software, aiming at the calculation of electronically excited states in complex molecular environments. Based on Gaussian-type atomic orbitals and making use of resolution of identity techniques, the code is designed specifically for nonperiodic systems. Application to a small molecule reference set successfully validates the methodology and its implementation for a variety of excitation types covering an energy range from 2 to 8 eV in single molecules.

Intermolecular Singlet and Triplet Exciton Transfer Integrals from Many-Body Green's Functions Theory.

A general approach to determine orientation and distance-dependent effective intermolecular exciton transfer integrals from many-body Green's functions theory is presented. On the basis of the GW approximation and the Bethe-Salpeter equation (BSE), a projection technique is employed to obtain the excitonic coupling by forming the expectation value of a supramolecular BSE Hamiltonian with electron-hole wave functions for excitations localized on two separated chromophores. Within this approach, accounting for the effects of coupling mediated by intermolecular charge transfer (CT) excitations is possible via perturbation theory or a reduction technique.

Getting excited: challenges in quantum-classical studies of excitons in polymeric systems.

A combination of classical molecular dynamics (MM/MD) and quantum chemical calculations based on the density functional theory (DFT) was performed to describe the conformational properties of diphenylethyne (DPE), methylated-DPE and poly para phenylene ethynylene (PPE). DFT calculations were employed to improve and develop force field parameters for MM/MD simulations. Many-body Green's function theory within the GW approximation and the Bethe-Salpeter (GW-BSE) equation were utilized to describe the excited states of the systems.

Impact of mesoscale order on open-circuit voltage in organic solar cells.

Structural order in organic solar cells is paramount: it reduces energetic disorder, boosts charge and exciton mobilities, and assists exciton splitting. Owing to spatial localization of electronic states, microscopic descriptions of photovoltaic processes tend to overlook the influence of structural features at the mesoscale. Long-range electrostatic interactions nevertheless probe this ordering, making local properties depend on the mesoscopic order.

Electronic Excitations in Push-Pull Oligomers and Their Complexes with Fullerene from Many-Body Green's Functions Theory with Polarizable Embedding.

We present a comparative study of excited states in push-pull oligomers of PCPDTBT and PSBTBT and prototypical complexes with a C60 acceptor using many-body Green's functions theory within the GW approximation and the Bethe-Salpeter equation. We analyze excitations in oligomers up to a length of 5 nm and find that for both materials the absorption energy practically saturates for structures larger than two repeat units due to the localized nature of the excitation. In the bimolecular complexes with C60, the transition from Frenkel to charge transfer excitons is generally exothermic and strongly influenced by the acceptor's position and orientation.

Parametrization of Extended Gaussian Disorder Models from Microscopic Charge Transport Simulations.

Simulations of organic semiconducting devices using drift-diffusion equations are vital for the understanding of their functionality as well as for the optimization of their performance. Input parameters for these equations are usually determined from experiments and do not provide a direct link to the chemical structures and material morphology. Here we demonstrate how such a parametrization can be performed by using atomic-scale (microscopic) simulations.

Can lattice models predict the density of states of amorphous organic semiconductors?

We extend existing lattice models of small-molecule amorphous semiconductors by accounting for changes in molecular polarizability upon charging or excitation. A compact expression of this contribution to the density of states is provided. Although the lattice model and the description based on a microscopic morphology both qualitatively predict an additional broadening, shift, and an exponential tail (traps) of the density of states, a quantitative agreement between the two cannot be achieved.

Frenkel and Charge-Transfer Excitations in Donor-acceptor Complexes from Many-Body Green's Functions Theory.

Excited states of donor-acceptor dimers are studied using many-body Green's functions theory within the GW approximation and the Bethe-Salpeter equation. For a series of prototypical small-molecule based pairs, this method predicts energies of local Frenkel and intermolecular charge-transfer excitations with the accuracy of tens of meV. Application to larger systems is possible and allowed us to analyze energy levels and binding energies of excitons in representative dimers of dicyanovinyl-substituted quarterthiophene and fullerene, a donor-acceptor pair used in state of the art organic solar cells.

Design rules for charge-transport efficient host materials for phosphorescent organic light-emitting diodes.

The use of blue phosphorescent emitters in organic light-emitting diodes (OLEDs) imposes demanding requirements on a host material. Among these are large triplet energies, the alignment of levels with respect to the emitter, the ability to form and sustain amorphous order, material processability, and an adequate charge carrier mobility. A possible design strategy is to choose a π-conjugated core with a high triplet level and to fulfill the other requirements by using suitable substituents.

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.

Excited States of Dicyanovinyl-Substituted Oligothiophenes from Many-Body Green's Functions Theory.

Excited states of dicyanovinyl-substituted oligothiophenes are studied using many-body Green's functions theory within the GW approximation and the Bethe-Salpeter equation. By varying the number of oligomer repeat units, we investigate the effects of resonant-antiresonant transition coupling, dynamical screening, and molecular conformations on calculated excitations. We find that the full dynamically screened Bethe-Salpeter equation yields absorption and emission energies in good agreement with experimental data.

Microscopic Simulations of Charge Transport in Disordered Organic Semiconductors.

Charge carrier dynamics in an organic semiconductor can often be described in terms of charge hopping between localized states. The hopping rates depend on electronic coupling elements, reorganization energies, and driving forces, which vary as a function of position and orientation of the molecules. The exact evaluation of these contributions in a molecular assembly is computationally prohibitive.