Transferable atomic multipole machine learning models for small organic molecules

Accurate representation of the molecular electrostatic potential, which is often expanded in distributed multipole moments, is crucial for an efficient evaluation of intermolecular interactions. Here we introduce a machine learning model for multipole coefficients of atom types H, C, O, N, S, F, and Cl in any molecular conformation. The model is trained on quantum chemical results for atoms in varying chemical environments drawn from thousands of organic molecules. Multipoles in systems with neutral, cationic, and anionic molecular charge states are treated with individual models. The models' predictive accuracy and applicability are illustrated by evaluating intermolecular interaction energies of nearly 1,000 dimers and the cohesive energy of the benzene crystal.Comment: 11 pages, 6 figure

Finite-Size Scaling of Charge Carrier Mobility in Disordered Organic Semiconductors

Simulations of charge transport in amorphous semiconductors are often performed in microscopically sized systems. As a result, charge carrier mobilities become system-size dependent. We propose a simple method for extrapolating a macroscopic, nondispersive mobility from the system-size dependence of a microscopic one. The method is validated against a temperature-based extrapolation [Phys. Rev. B 82, 193202 (2010)]. In addition, we provide an analytic estimate of system sizes required to perform nondispersive charge transport simulations in systems with finite charge carrier density, derived from a truncated Gaussian distribution. This estimate is not limited to lattice models or specific rate expressions

Atomistic simulation of structure and dynamics of columnar phases of hexabenzocoronene derivatives

Using atomistic molecular dynamics simulations we study solid and liquid crystalline columnar discotic phases formed by alkyl-substituted hexabenzocoronene mesogens. Correlations between the molecular structure, packing, and dynamical properties of these materials are established

Twist-bend instability for toroidal DNA condensates

We propose that semiflexible polymers in poor solvent collapse in two stages. The first stage is the well known formation of a dense toroidal aggregate. However, if the solvent is sufficiently poor, the condensate will undergo a second structural transition to a twisted entangled state, in which individual filaments lower their bending energy by additionally orbiting around the mean path along which they wind. This ``topological ripening'' is consistent with known simulations and experimental results. It connects and rationalizes various experimental observations ranging from strong DNA entanglement in viral capsids to the unusually short pitch of the cholesteric phase of DNA in sperm-heads. We propose that topological ripening of DNA toroids could improve the efficiency and stability of gene delivery.Comment: 4 pages, 3 figures, RevTeX4 styl

Coarse-grained interaction potentials for polyaromatic hydrocarbons

Using Kohn-Sham density functional theory (KS-DFT), we have studied the interaction between various polyaromatic hydrocarbon molecules. The systems range from mono-cyclic benzene up to hexabenzocoronene (hbc). For several conventional exchange-correlation functionals potential energy curves of interaction of the $\pi$-$\pi$ stacking hbc dimer are reported. It is found that all pure local density or generalized gradient approximated functionals yield qualitatively incorrect predictions regarding structure and interaction. Inclusion of a non-local, atom-centered correction to the KS-Hamiltonian enables quantitative predictions. The computed potential energy surfaces of interaction yield parameters for a coarse-grained potential, which can be employed to study discotic liquid-crystalline mesophases of derived polyaromatic macromolecules

Charge mobility of discotic mesophases: A multiscale quantum/classical study

A correlation is established between the molecular structure and charge mobility of discotic mesophases of hexabenzocoronene derivatives by combining electronic structure calculations, Molecular Dynamics, and kinetic Monte Carlo simulations. It is demonstrated that this multiscale approach can provide an accurate ab-initio description of charge transport in organic materials

Stark Effect of Hybrid Charge Transfer States at Planar ZnO/Organic Interfaces

We investigate the bias-dependence of the hybrid charge transfer state emission at planar heterojunctions between the metal oxide acceptor ZnO and three donor molecules. The electroluminescence peak energy linearly increases with the applied bias, saturating at high fields. Variation of the organic layer thickness and deliberate change of the ZnO conductivity through controlled photo-doping allow us to confirm that this bias-induced spectral shifts relate to the internal electric field in the organic layer rather than the filling of states at the hybrid interface. We show that existing continuum models overestimate the hole delocalization and propose a simple electrostatic model in which the linear and quadratic Stark effects are explained by the electrostatic interaction of a strongly polarizable molecular cation with its mirror image

Liquid crystal director fluctuations and surface anchoring by molecular simulation

We propose a simple and reliable method to measure the liquid crystal surface anchoring strength by molecular simulation. The method is based on the measurement of the long-range fluctuation modes of the director in confined geometry. As an example, molecular simulations of a liquid crystal in slab geometry between parallel walls with homeotropic anchoring have been carried out using the Monte Carlo technique. By studying different slab thicknesses, we are able to calculate separately the position of the elastic boundary condition, and the extrapolation length

A general framework for consistent estimation of charge transport properties via random walks in random environments

A general framework is proposed for the study of the charge transport properties of materials via random walks in random environments (RWRE). The material of interest is modeled by a random environment, and the charge carrier is modeled by a random walker. The framework combines a model for the fast generation of random environments that realistically mimic materials morphology with an algorithm for efficient estimation of key properties of the resulting random walk. The model of the environment makes use of tools from spatial statistics and the theory of random geometric graphs. More precisely, the disordered medium is represented by a random spatial graph with directed edge weights, where the edge weights represent the transition rates of a Markov jump process (MJP) modeling the motion of the random walker. This MJP is a multiscale stochastic process. In the long term, it explores all vertices of the random graph model. In the short term, however, it becomes trapped in small subsets of the state space and makes many transitions in these small regions. This behavior makes efficient estimation of velocity by Monte Carlo simulations a challenging task. Therefore, we use aggregate Monte Carlo (AMC), introduced in [T. Brereton et al., Methodol. Comput. Appl. Probab., 16 (2014), pp. 465-484], for estimating the velocity of a random walker as it passes through a realization of the random environment. In this paper, we prove the strong consistency of the AMC velocity estimator and use this result to conduct a detailed case study, in which we describe the motion of holes in an amorphous mesophase of an organic semiconductor, dicyanovinyl-substituted oligothiophene (DCV4T). In particular, we analyze the effect of system size (i.e., number of molecules) on the velocity of single charge carriers

Defect structures and torque on an elongated colloidal particle immersed in a liquid crystal host

Combining molecular dynamics and Monte Carlo simulation we study defect structures around an elongated colloidal particle embedded in a nematic liquid crystal host. By studying nematic ordering near the particle and the disclination core region we are able to examine the defect core structure and the difference between two simulation techniques. In addition, we also study the torque on a particle tilted with respect to the director, and modification of this torque when the particle is close to the cell wall