Martini 3: a general purpose force field for coarse-grained molecular dynamics.

The coarse-grained Martini force field is widely used in biomolecular simulations. Here we present the refined model, Martini 3 ( http://cgmartini.nl ), with an improved interaction balance, new bead types and expanded ability to include specific interactions representing, for example, hydrogen bonding and electronic polarizability. The updated model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.

High-efficiency stretchable light-emitting polymers from thermally activated delayed fluorescence.

Stretchable light-emitting materials are the key components for realizing skin-like displays and optical biostimulation. All the stretchable emitters reported to date, to the best of our knowledge, have been based on electroluminescent polymers that only harness singlet excitons, limiting their theoretical quantum yield to 25%. Here we present a design concept for imparting stretchability onto electroluminescent polymers that can harness all the excitons through thermally activated delayed fluorescence, thereby reaching a near-unity theoretical quantum yield. We show that our design strategy of inserting flexible, linear units into a polymer backbone can substantially increase the mechanical stretchability without affecting the underlying electroluminescent processes. As a result, our synthesized polymer achieves a stretchability of 125%, with an external quantum efficiency of 10%. Furthermore, we demonstrate a fully stretchable organic light-emitting diode, confirming that the proposed stretchable thermally activated delayed fluorescence polymers provide a path towards simultaneously achieving desirable electroluminescent and mechanical characteristics, including high efficiency, brightness, switching speed and stretchability as well as low driving voltage.

Comparing Dimerization Free Energies and Binding Modes of Small Aromatic Molecules with Different Force Fields.

Dimerization free energies are fundamental quantities that describe the strength of interaction of different molecules. Obtaining accurate experimental values for small molecules and disentangling the conformations that contribute most to the binding can be extremely difficult, due to the size of the systems and the small energy differences. In many cases, one has to resort to computational methods to calculate such properties. In this work, we used molecular dynamics simulations in conjunction with metadynamics to calculate the free energy of dimerization of small aromatic rings, and compared three models from popular online servers for atomistic force fields, namely G54a7, CHARMM36 and OPLS. We show that, regardless of the force field, the profiles for the dimerization free energy of these compounds are very similar. However, significant care needs to be taken when studying larger molecules, since the deviations from the trends increase with the size of the molecules, resulting in force field dependent preferred stacking modes; for example, in the cases of pyrene and tetracene. Our results provide a useful background study for using topology builders to model systems which rely on stacking of aromatic moieties, and are relevant in areas ranging from drug design to supramolecular assembly.

Molecular versus Excitonic Disorder in Individual Artificial Light-Harvesting Systems.

Natural light-harvesting antennae employ a dense array of chromophores to optimize energy transport via the formation of delocalized excited states (excitons), which are critically sensitive to spatio-energetic variations of the molecular structure. Identifying the origin and impact of such variations is highly desirable for understanding and predicting functional properties yet hard to achieve due to averaging of many overlapping responses from individual systems. Here, we overcome this problem by measuring the heterogeneity of synthetic analogues of natural antennae-self-assembled molecular nanotubes-by two complementary approaches: single-nanotube photoluminescence spectroscopy and ultrafast 2D correlation. We demonstrate remarkable homogeneity of the nanotube ensemble and reveal that ultrafast (∼50 fs) modulation of the exciton frequencies governs spectral broadening. Using multiscale exciton modeling, we show that the dominance of homogeneous broadening at the exciton level results from exchange narrowing of strong static disorder found for individual molecules within the nanotube. The detailed characterization of static and dynamic disorder at the exciton as well as the molecular level presented here opens new avenues in analyzing and predicting dynamic exciton properties, such as excitation energy transport.

Crystal Field in Rare-Earth Complexes: From Electrostatics to Bonding.

The flexibility of first-principles (ab initio) calculations with the SO-CASSCF (complete active space self-consistent field theory with a treatment of the spin-orbit (SO) coupling by state interaction) method is used to quantify the electrostatic and covalent contributions to crystal field parameters. Two types of systems are chosen for illustration: 1) The ionic and experimentally well-characterized PrCl3 crystal; this study permits a revisitation of the partition of contributions proposed in the early days of crystal field theory; and 2) a series of sandwich molecules [Ln(ηn -Cn Hn )2 ]q , with Ln=Dy, Ho, Er, and Tm and n=5, 6, and 8, in which the interaction between LnIII and the aromatic ligands is more difficult to describe within an electrostatic approach. It is shown that a model with three layers of charges reproduces the electrostatic field generated by the ligands and that the covalency plays a qualitative role. The one-electron character of crystal field theory is discussed and shown to be valuable, although it is not completely quantitative. This permits a reduction of the many-electron problem to a discussion of the energy of the seven 4f orbitals.

Can the Dielectric Constant of Fullerene Derivatives Be Enhanced by Side-Chain Manipulation? A Predictive First-Principles Computational Study.

The low efficiency of organic photovoltaic (OPV) devices has often been attributed to the strong Coulombic interactions between the electron and hole, impeding the charge separation process. Recently, it has been argued that by increasing the dielectric constant of materials used in OPVs, this strong interaction could be screened. In this work, we report the application of periodic density functional theory together with the coupled perturbed Kohn-Sham method to calculate the electronic contribution to the dielectric constant for fullerene C60 derivatives, a ubiquitous class of molecules in the field of OPVs. The results show good agreement with experimental data when available and also reveal an important undesirable outcome when manipulating the side chain to maximize the static dielectric constant: in all cases, the electronic contribution to the dielectric constant decreases as the side chain increases in size. This information should encourage both theoreticians and experimentalists to further investigate the relevance of contributions to the dielectric constant from slower processes like vibrations and dipolar reorientations for facilitating the charge separation, because electronically, enlarging the side chain of conventional fullerene derivatives only lowers the dielectric constant, and consequently, their electronic dielectric constant is upper bound by the one of C60.

Bulk Heterojunction Morphologies with Atomistic Resolution from Coarse-Grain Solvent Evaporation Simulations.

Control over the morphology of the active layer of bulk heterojunction (BHJ) organic solar cells is paramount to achieve high-efficiency devices. However, no method currently available can predict morphologies for a novel donor-acceptor blend. An approach which allows reaching relevant length scales, retaining chemical specificity, and mimicking experimental fabrication conditions, and which is suited for high-throughput schemes has been proven challenging to find. Here, we propose a method to generate atom-resolved morphologies of BHJs which conforms to these requirements. Coarse-grain (CG) molecular dynamics simulations are employed to simulate the large-scale morphological organization during solution-processing. The use of CG models which retain chemical specificity translates into a direct path to the rational design of donor and acceptor compounds which differ only slightly in chemical nature. Finally, the direct retrieval of fully atomistic detail is possible through backmapping, opening the way for improved quantum mechanical calculations addressing the charge separation mechanism. The method is illustrated for the poly(3-hexyl-thiophene) (P3HT)-phenyl-C61-butyric acid methyl ester (PCBM) mixture, and found to predict morphologies in agreement with experimental data. The effect of drying rate, P3HT molecular weight, and thermal annealing are investigated extensively, resulting in trends mimicking experimental findings. The proposed methodology can help reduce the parameter space which has to be explored before obtaining optimal morphologies not only for BHJ solar cells but also for any other solution-processed soft matter device.

Structural, Ionic, and Electronic Properties of Solid-State Phthalimide-Containing Polymers for All-Organic Batteries.

Redox-active polymers serving as the active materials in solid-state electrodes offer a promising path toward realizing all-organic batteries. While both cathodic and anodic redox-active polymers are needed, the diversity of the available anodic materials is limited. Here, we predict solid-state structural, ionic, and electronic properties of anodic, phthalimide-containing polymers using a multiscale approach that combines atomistic molecular dynamics, electronic structure calculations, and machine learning surrogate models. Importantly, by combining information from each of these scales, we are able to bridge the gap between bottom-up molecular characteristics and macroscopic properties such as apparent diffusion coefficients of electron transport (D app). We investigate the impact of different polymer backbones and of two critical factors during battery operation: state of charge and polymer swelling. Our findings reveal that the state of charge significantly influences solid-state packing and the thermophysical properties of the polymers, which, in turn, affect ionic and electronic transport. A combination of molecular-level properties (such as the reorganization energy) and condensed-phase properties (such as effective electron hopping distances) determine the predicted ranking of electron transport capabilities of the polymers. We predict D app for the phthalimide-based polymers and for a reference nitroxide radical-based polymer, finding a 3 orders of magnitude increase in D app (≈10-6 cm2 s-1) with respect to the reference. This study underscores the promise of phthalimide-containing polymers as highly capable redox-active polymers for anodic materials in all-organic batteries, due to their exceptional predicted electron transport capabilities.

Bartender: Martini 3 Bonded Terms via Quantum Mechanics-Based Molecular Dynamics.

Coarse-grained (CG) molecular dynamics (MD) simulations have grown in applicability over the years. The recently released version of the Martini CG force field (Martini 3) has been successfully applied to simulate many processes, including protein-ligand binding. However, the current ligand parametrization scheme is manual and requires an a priori reference all-atom (AA) simulation for benchmarking. For systems with suboptimal AA parameters, which are often unknown, this translates into a CG model that does not reproduce the true dynamical behavior of the underlying molecule. Here, we present Bartender, a quantum mechanics (QM)/MD-based parametrization tool written in Go. Bartender harnesses the power of QM simulations and produces reasonable bonded terms for Martini 3 CG models of small molecules in an efficient and user-friendly manner. For small, ring-like molecules, Bartender generates models whose properties are indistinguishable from the human-made models. For more complex, drug-like ligands, it is able to fit functional forms beyond simple harmonic dihedrals and thus better captures their dynamical behavior. Bartender has the power to both increase the efficiency and the accuracy of Martini 3-based high-throughput applications by producing numerically stable and physically realistic CG models.

Polyply; a python suite for facilitating simulations of macromolecules and nanomaterials.

Molecular dynamics simulations play an increasingly important role in the rational design of (nano)-materials and in the study of biomacromolecules. However, generating input files and realistic starting coordinates for these simulations is a major bottleneck, especially for high throughput protocols and for complex multi-component systems. To eliminate this bottleneck, we present the polyply software suite that provides 1) a multi-scale graph matching algorithm designed to generate parameters quickly and for arbitrarily complex polymeric topologies, and 2) a generic multi-scale random walk protocol capable of setting up complex systems efficiently and independent of the target force-field or model resolution. We benchmark quality and performance of the approach by creating realistic coordinates for polymer melt simulations, single-stranded as well as circular single-stranded DNA. We further demonstrate the power of our approach by setting up a microphase-separated block copolymer system, and by generating a liquid-liquid phase separated system inside a lipid vesicle.

Strategies for Enhancing the Dielectric Constant of Organic Materials.

High dielectric constant organic semiconductors, often obtained by the use of ethylene glycol (EG) side chains, have gained attention in recent years in the efforts of improving the device performance for various applications. Dielectric constant enhancements due to EGs have been demonstrated extensively, but various effects, such as the choice of the particular molecule and the frequency and temperature regime, that determine the extent of this enhancement require further understanding. In this work, we study these effects by means of polarizable molecular dynamics simulations on a carefully selected set of fullerene derivatives with EG side chains. The selection allows studying the dielectric response in terms of both the number and length of EG chains and also the choice of the group connecting the fullerene to the EG chain. The computed time- and frequency-dependent dielectric responses reveal that the experimentally observed rise of the dielectric constant within the kilo/megahertz regime for some molecules is likely due to the highly stretched dielectric response of the EGs: the initial sharp increase over the first few nanoseconds is followed by a smaller but persistent increase in the range of microseconds. Additionally, our computational protocol allows the separation of different factors that contribute to the overall dielectric constant, providing insights to make several molecular design guides for future organic materials in order to enhance their dielectric constant further.

The Martini Model in Materials Science.

The Martini model, a coarse-grained force field initially developed with biomolecular simulations in mind, has found an increasing number of applications in the field of soft materials science. The model's underlying building block principle does not pose restrictions on its application beyond biomolecular systems. Here, the main applications to date of the Martini model in materials science are highlighted, and a perspective for the future developments in this field is given, particularly in light of recent developments such as the new version of the model, Martini 3.

Fullerene derivatives with oligoethylene-glycol side chains: an investigation on the origin of their outstanding transport properties.

For many years, fullerene derivatives have been the main n-type material of organic electronics and optoelectronics. Recently, fullerene derivatives functionalized with ethylene glycol (EG) side chains have been showing important properties such as enhanced dielectric constants, facile doping and enhanced self-assembly capabilities. Here, we have prepared field-effect transistors using a series of these fullerene derivatives equipped with EG side chains of different lengths. Transport data show the beneficial effect of increasing the EG side chain. In order to understand the material properties, full structural determination of these fullerene derivatives has been achieved by coupling the X-ray data with molecular dynamics (MD) simulations. The increase in transport properties is paired with the formation of extended layered structures, efficient molecular packing and an increase in the crystallite alignment. The layer-like structure is composed of conducting layers, containing of closely packed C60 balls approaching the inter-distance of 1 nm, that are separated by well-defined EG layers, where the EG chains are rather splayed with the chain direction almost perpendicular to the layer normal. Such a layered structure appears highly ordered and highly aligned with the C60 planes oriented parallel to the substrate in the thin film configuration. The order inside the thin film increases with the EG chain length, allowing the systems to achieve mobilities as high as 0.053 cm2 V-1 s-1. Our work elucidates the structure of these interesting semiconducting organic molecules and shows that the synergistic use of X-ray structural analysis and MD simulations is a powerful tool to identify the structure of thin organic films for optoelectronic applications.

How Ethylene Glycol Chains Enhance the Dielectric Constant of Organic Semiconductors: Molecular Origin and Frequency Dependence.

Incorporating ethylene glycols (EGs) into organic semiconductors has become the prominent strategy to increase their dielectric constant. However, EG's contribution to the dielectric constant is due to nuclear relaxations, and therefore, its relevance for various organic electronic applications depends on the time scale of these relaxations, which remains unknown. In this work, by means of a new computational protocol based on polarizable molecular dynamics simulations, the time- and frequency-dependent dielectric constant of a representative fullerene derivative with EG side chains is predicted, the origin of its unusually high dielectric constant is explained, and design suggestions are made to further increase it. Finally, a dielectric relaxation time of ∼1 ns is extracted which suggests that EGs may be too slow to reduce the Coulombic screening in organic photovoltaics but are definitely fast enough for organic thermoelectrics with much lower charge carrier velocities.

Protein-ligand binding with the coarse-grained Martini model.

The detailed understanding of the binding of small molecules to proteins is the key for the development of novel drugs or to increase the acceptance of substrates by enzymes. Nowadays, computer-aided design of protein-ligand binding is an important tool to accomplish this task. Current approaches typically rely on high-throughput docking essays or computationally expensive atomistic molecular dynamics simulations. Here, we present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein-ligand interactions of small drug-like molecules. Remarkably, we achieve high accuracy without the need of any a priori knowledge of binding pockets or pathways. Our approach is applied to a range of systems from the well-characterized T4 lysozyme over members of the GPCR family and nuclear receptors to a variety of enzymes. The presented results open the way to high-throughput screening of ligand libraries or protein mutations using the coarse-grained Martini model.

Multiscale modeling of molecular structure and optical properties of complex supramolecular aggregates.

Supramolecular aggregates of synthetic dye molecules offer great perspectives to prepare biomimetic functional materials for light-harvesting and energy transport. The design is complicated by the fact that structure-property relationships are hard to establish, because the molecular packing results from a delicate balance of interactions and the excitonic properties that dictate the optics and excited state dynamics, in turn sensitively depend on this packing. Here we show how an iterative multiscale approach combining molecular dynamics and quantum mechanical exciton modeling can be used to obtain accurate insight into the packing of thousands of cyanine dye molecules in a complex double-walled tubular aggregate in close interaction with its solvent environment. Our approach allows us to answer open questions not only on the structure of these prototypical aggregates, but also about their molecular-scale structural and energetic heterogeneity, as well as on the microscopic origin of their photophysical properties. This opens the route to accurate predictions of energy transport and other functional properties.

N-type organic thermoelectrics: demonstration of ZT > 0.3.

The 'phonon-glass electron-crystal' concept has triggered most of the progress that has been achieved in inorganic thermoelectrics in the past two decades. Organic thermoelectric materials, unlike their inorganic counterparts, exhibit molecular diversity, flexible mechanical properties and easy fabrication, and are mostly 'phonon glasses'. However, the thermoelectric performances of these organic materials are largely limited by low molecular order and they are therefore far from being 'electron crystals'. Here, we report a molecularly n-doped fullerene derivative with meticulous design of the side chain that approaches an organic 'PGEC' thermoelectric material. This thermoelectric material exhibits an excellent electrical conductivity of >10 S cm-1 and an ultralow thermal conductivity of <0.1 Wm-1K-1, leading to the best figure of merit ZT = 0.34 (at 120 °C) among all reported single-host n-type organic thermoelectric materials. The key factor to achieving the record performance is to use 'arm-shaped' double-triethylene-glycol-type side chains, which not only offer excellent doping efficiency (~60%) but also induce a disorder-to-order transition upon thermal annealing. This study illustrates the vast potential of organic semiconductors as thermoelectric materials.

Pitfalls of the Martini Model.

The computational and conceptual simplifications realized by coarse-grain (CG) models make them a ubiquitous tool in the current computational modeling landscape. Building block based CG models, such as the Martini model, possess the key advantage of allowing for a broad range of applications without the need to reparametrize the force field each time. However, there are certain inherent limitations to this approach, which we investigate in detail in this work. We first study the consequences of the absence of specific cross Lennard-Jones parameters between different particle sizes. We show that this lack may lead to artificially high free energy barriers in dimerization profiles. We then look at the effect of deviating too far from the standard bonded parameters, both in terms of solute partitioning behavior and solvent properties. Moreover, we show that too weak bonded force constants entail the risk of artificially inducing clustering, which has to be taken into account when designing elastic network models for proteins. These results have implications for the current use of the Martini CG model and provide clear directions for the reparametrization of the Martini model. Moreover, our findings are generally relevant for the parametrization of any other building block based force field.

Enhancing Molecular n-Type Doping of Donor-Acceptor Copolymers by Tailoring Side Chains.

In this contribution, for the first time, the molecular n-doping of a donor-acceptor (D-A) copolymer achieving 200-fold enhancement of electrical conductivity by rationally tailoring the side chains without changing its D-A backbone is successfully improved. Instead of the traditional alkyl side chains for poly{[N,N'-bis(2-octyldodecyl)-naphthalene-1,4,5,8-bis(dicarboximide)-2,6-diyl](NDI)-alt-5,5'-(2,2'-bithiophene)} (N2200), polar triethylene glycol type side chains is utilized and a high electrical conductivity of 0.17 S cm-1 after doping with (4-(1,3-dimethyl-2,3-dihydro-1H-benzoimidazol-2-yl)phenyl)dimethylamine is achieved, which is the highest reported value for n-type D-A copolymers. Coarse-grained molecular dynamics simulations indicate that the polar side chains can significantly reduce the clustering of dopant molecules and favor the dispersion of the dopant in the host matrix as compared to the traditional alkyl side chains. Accordingly, intimate contact between the host and dopant molecules in the NDI-based copolymer with polar side chains facilitates molecular doping with increased doping efficiency and electrical conductivity. For the first time, a heterogeneous thermoelectric transport model for such a material is proposed, that is the percolation of charge carriers from conducting ordered regions through poorly conductive disordered regions, which provides pointers for further increase in the themoelectric properties of n-type D-A copolymers.