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 ab initio per-molecule electrostatic potentials in the liquid phase using the split charge equilibration of Nistor et al. [J. Chem. Phys. 125, 094108 (2006)]. By taking the density and diffusion dependence on temperature as target properties, significant improvements to the behavior of isothermal compressibility were achieved along with improvements to other thermodynamic and rotational properties. While exhibiting a dipole moment close to ab initio values, TIP5P/2018 suffers from a too small quadrupole moment due to the charge assignment procedure and results in an overestimation of the dielectric constant.

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. Using a technique developed to address spatially aperiodic excitations in thin films and in bulk, we show how inclusion of mesoscale order resolves the controversy between experimental and theoretical results for the energy-level profile and alignment in a variety of photovoltaic systems, with direct experimental validation. Optimal use of long-range ordering also rationalizes the acceptor-donor-acceptor paradigm for molecular design of donor dyes. We predict open-circuit voltages of planar heterojunction solar cells in excellent agreement with experimental data, based only on crystal structures and interfacial orientation.

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. The reliability of the excitation energies based on the MM/MD conformations was examined and compared to the excitation energies from DFT conformations. The results show an overall agreement between the optical excitations based on MM/MD conformations and DFT conformations. This allows for the calculation of excitation energies based on MM/MD conformations.

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 GW approximation with the Bethe-Salpeter equation (GW-BSE). Using a small molecule dimer system, we first study the influence of different diabatization methods, as well as different model choices within GW-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 GW-BSE-molecular mechanics calculations. For these disordered systems of bulky molecules, we observe differences in the couplings based on the Edmiston-Ruedenberg diabatization compared to the more approximate Generalize Mulliken-Hush and fragment charge difference diabatization formalisms. In a kinetic model for the conversion between LE and CT states, these differences affect the details of state populations in an intermediate time scale but not the final populations.

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

We present quantum-quantum and quantum-quantum-classical schemes based on many-body Green's functions theory in the GW approximation with the Bethe-Salpeter equation (GW-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 GW-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 GW-BSE, results from such PbE-GW-BSE calculations can deviate from a full-system reference. Here, we scrutinize in detail, e.g., the individual and combined effects of different choices of active regions, the influence of omitting the screening from the inactive region, and strategies for basis set truncation on frontier orbital and near-gap electron-hole excitation energies. As prototypical systems, we consider a diketopyrrolopyrrole bicyclic ring including side-chains, a polarity-sensitive dye (prodan) in aqueous environment, and a π-stacked dimer of benzene and tetracyanoethylene in water, respectively, covering a variety of excitation characters in molecular systems with complex chemical environments and photoinduced processes. Our results suggest that to obtain agreement of approximately 0.1 eV between near-gap excitation energies from embedded and full calculations, the active region should be chosen based on the Mulliken population of the full highest-occupied molecular orbital and that careful benchmarking should be done on the KS level before the actual GW-BSE steps when basis set truncation is used. We find that PbE-GW-BSE offers significant reductions in computation times and, more importantly, memory requirements, making calculations for considerably larger systems tractable.

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. Bulky substituents, however, induce large spatial separations between conjugated cores, can substantially reduce intermolecular electronic couplings, and decrease the charge mobility of the host. In this work we analyze charge transport in amorphous 2,8-bis(triphenylsilyl)dibenzofuran, an electron-transporting material synthesized to serve as a host in deep-blue OLEDs. We show that mesomeric effects delocalize the frontier orbitals over the substituents recovering strong electronic couplings and lowering reorganization energies, especially for electrons, while keeping energetic disorder small. Admittance spectroscopy measurements reveal that the material has indeed a high electron mobility and a small Poole-Frenkel slope, supporting our conclusions. By linking electronic structure, molecular packing, and mobility, we provide a pathway to the rational design of hosts with high charge mobilities.

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.

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.

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. Our results show that the effective solvent evaporation rates at room temperature only have a minor influence on the segmental mobility and intermolecular orientation, characterized in all cases by a preferential arrangement of PCBM[60] close to the electron-donating substituents in DPP2Py2T. In contrast, solvent-free annealing from a liquid yields clustering of the fullerene close to the electron-withdrawing DPP, generally considered to be detrimental for application in organic solar cells. We find that the difference can be attributed to differences in the behavior of 2-hexyldecyl side-chains, which collapse toward DPP when solvent is explicitly removed, thereby blocking access of PCBM[60].

The Goos-Hänchen effect for surface plasmon polaritons.

By means of an impedance boundary condition and numerical solution of integral equations for the scattering amplitudes to which its use gives rise, we study as a function of its angle of incidence the reflection of a surface plasmon polariton beam propagating on a metal surface whose dielectric function is ɛ1(ω) when it is incident on a planar interface with a coplanar metal surface whose dielectric function is ɛ2(ω). When the surface of incidence is optically more dense than the surface of scattering, i.e. when |ɛ2(ω)|≫|ɛ1(ω)|, the reflected beam undergoes a lateral displacement whose magnitude is several times the wavelength of the incident beam. This displacement is the surface plasmon polariton analogue of the Goos-Hänchen effect. Since this displacement is sensitive to the dielectric properties of the surface, this effect can be exploited to sense modifications of the dielectric environment of a metal surface, e.g. due to adsorption of atomic or molecular layers on it.

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.1 eV. We further demonstrate that the OSR-based ΔMLQP captures the effects of intra- and intermolecular conformations in application to water monomers and dimers. Finally, we show that the approach can be embedded in multiscale simulation workflows, by studying the solvatochromic shifts of quasiparticle and electron-hole excitation energies of solvated acetone in a setup combining molecular dynamics, DFT, the GW approximation, and the Bethe-Salpeter equation. Our findings suggest that the ΔMLQP model allows us to predict quasiparticle energies and photoelectron spectra of molecules and clusters with GW accuracy at DFT cost.

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 GW approximation and the Bethe-Salpeter equation (BSE) employing numerical gradients. We scrutinize the influence of several typical approximations in the GW-BSE framework; we used one-shot G0W0 or eigenvalue self-consistent evGW, 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). We find overall a good agreement of the structural parameters optimized with the GW-BSE calculations with CASPT2, with an average relative error of around 1% for the G0W0 and 1.5% for the evGW variants based on a PBE0 ground state, respectively, while the other approximations have negligible influence. The relative errors are also smaller than those for CC2 and TDDFT with different functionals and only larger than VMC, indicating that the GW-BSE method does not only yield excitation energies but also geometries in good agreement with established higher-order wave function methods.

Charge transport in organic crystals: role of disorder and topological connectivity.

We analyze the relationship among the molecular structure, morphology, percolation network, and charge carrier mobility in four organic crystals: rubrene, indolo[2,3-b]carbazole with CH(3) side chains, and benzo[1,2-b:4,5-b']bis[b]benzothiophene derivatives with and without C(4)H(9) side chains. Morphologies are generated using an all-atom force field, while charge dynamics is simulated within the framework of high-temperature nonadiabatic Marcus theory or using semiclassical dynamics. We conclude that, on the length scales reachable by molecular dynamics simulations, the charge transport in bulk molecular crystals is mostly limited by the dynamic disorder, while in self-assembled monolayers the static disorder, which is due to the slow motion of the side chains, enhances charge localization and influences the transport dynamics. We find that the presence of disorder can either reduce or increase charge carrier mobility, depending on the dimensionality of the charge percolation network. The advantages of charge transporting materials with two- or three-dimensional networks are clearly shown.

Density-functional based determination of intermolecular charge transfer properties for large-scale morphologies.

Theoretical studies of charge transport in organic conducting systems pose a unique challenge since they must describe both extremely short-ranged and fast processes (charge tunneling) and extremely long-ranged and slow ones (molecular ordering). The description of the mobility of electrons and holes in the hopping regime relies on the determination of intermolecular hopping rates in large-scale morphologies. Using Marcus theory these rates can be calculated from intermolecular transfer integrals and on-site energies. Here we present a detailed computational study on the accuracy and efficiency of density-functional theory based approaches to the determination of intermolecular transfer integrals. First, it is demonstrated how these can be obtained from quantum-chemistry calculations by forming the expectation value of a dimer Fock operator with frontier orbitals of two neighboring monomers based on a projective approach. We then consider the prototypical example of one pair out of a larger morphology of tris(8-hydroxyquinolinato)aluminium (Alq(3)) and study the influence of computational parameters, e.g. the choice of basis sets, exchange-correlation functional, and convergence criteria, on the calculated transfer integrals. The respective results are compared in order to derive an optimal strategy for future simulations based on the full morphology.

Charge transport in columnar mesophases of carbazole macrocycles.

Charge transport properties of a columnar mesophase of carbazole macrocycles are analyzed. Realistic morphologies are sampled using all-atom molecular dynamics simulations while charge transport is simulated using the kinetic Monte Carlo method with transfer rates obtained from the high temperature nonadiabatic limit of Marcus theory. It is shown that the molecular design with side chains pointing inside the macrocycle allows close approach between molecules of neighboring columns, thus enabling three-dimensional transport and helping to circumvent charge trapping on structural defects.

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 (DPP2PymT), 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 m = 1, 2, 3 in good and bad solvents and of melts. We observe the expected dependence of the chain conformation on the solvent quality, with the chain collapsing in water, and swelling in chloroform. The glass-transition temperature for the polymer melts is found to be in the range of 340-370 K. Analysis of the mobility of the conjugated segments in the polymer backbone reveals two relaxation processes: a fast one with a characteristic time at room temperature on the order of 10 ps associated with nearly harmonic vibrations and a slow one on the order of 100 ns associated with temperature-activated cis-trans transitions.

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. Here, an approach is presented in which (i) the local environment is sampled by a classical molecular dynamics (MD) simulation of the dye and solvent, (ii) the electronically excited state of prodan upon light absorption is predicted by numerical quantum mechanics (QM), (iii) the iterative relaxation of the environment around the excited dye by MD coupled with the evolution of the excited state is performed, and (iv) the emission properties are predicted by QM. The QM steps are computed using the many-body Green's function in the GW approximation and the Bethe-Salpeter equation with the environment modeled as fixed point charges, sampled in the MD simulation steps. The comparison to ultrafast time-resolved transient absorption measurements demonstrates that the iterative molecular mechanics (MM)/QM approach agrees quantitatively with both the polarity-dependent shift in emission and the time scale over which the charge transfer state is stabilized. Together the simulations and experimental measurements suggest that the evolution into the charge transfer state is slower in amphiphilic solvents.

Cloaking from surface plasmon polaritons by a circular array of point scatterers.

In recent years it has been demonstrated both theoretically and experimentally that it is possible to cloak a predefined region of space from interaction with external volume electromagnetic waves, rendering an arbitrary object inside this region invisible to an outside observer. The several strategies that have been developed for achieving such cloaking cannot be applied directly to the cloaking of a surface feature from surface plasmon polaritons propagating on that surface. Here we demonstrate that it is possible to generate an arrangement of two concentric rings of point scatterers on a metal surface that significantly reduces the scattering of surface plasmon polaritons from an object enclosed within this circular structure.

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. As a prototypical system, we consider hole transport in amorphous tris(8-hydroxyquinolinato)aluminum. Reference data for training and validation is obtained from multiscale ab initio simulations, in which coupling elements are evaluated using density-functional theory, for a system containing 4096 molecules. The Coulomb matrix representation is chosen to encode the explicit molecular pair coordinates into a rotation and translation invariant feature set for the FFNN. The final optimized deep feedforward neural network is tested for transport models without and with energetic disorder. It predicts electronic coupling elements and mobilities in excellent agreement with the reference data. Such a FFNN is readily applicable to much larger systems at negligible computational cost, providing a powerful surrogate model to overcome the size limitations of the ab initio approach.

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.2014, 6, 126-132). Peculiarly, these results suggest an apparent structural change upon increasing the length of the polymers. To address this puzzle, we use single-chain models to study the stiffness of polyphenylene precursors in a theta-like solvent as a function of chain composition and monomer sequence. Steric effects are isolated by considering random walk chains with segment length distributions and the position of monomers determined by the nature of the arene substitution along the backbone. Moreover, two homopolymer limiting cases are defined, that is, meta and para sequences, by associating two types of monomers to each possible substitution pattern. We consider, within these two limiting cases, chains with different compositions and monomer sequences. We compute persistence lengths, mean square end-to-end distances, and gyration and hydrodynamic radii. We find that distinct values of the persistence length for apparently the same chain chemistry are the result of different mixing ratios and the arrangement along the chain of the two positional isomers of the same monomer. Finally, we discuss the relation between two-dimensional density of the number of crossings and the length of polyphenylene segments as they would occur upon strong chain adsorption onto a substrate.