Accessing the electronic structure of liquid crystalline semiconductors with bottom-up electronic coarse-graining.

Understanding the relationship between multiscale morphology and electronic structure is a grand challenge for semiconducting soft materials. Computational studies aimed at characterizing these relationships require the complex integration of quantum-chemical (QC) calculations, all-atom and coarse-grained (CG) molecular dynamics simulations, and back-mapping approaches. However, these methods pose substantial computational challenges that limit their application to the requisite length scales of soft material morphologies. Here, we demonstrate the bottom-up electronic coarse-graining (ECG) of morphology-dependent electronic structure in the liquid-crystal-forming semiconductor, 2-(4-methoxyphenyl)-7-octyl-benzothienobenzothiophene (BTBT). ECG is applied to construct density functional theory (DFT)-accurate valence band Hamiltonians of the isotropic and smectic liquid crystal (LC) phases using only the CG representation of BTBT. By bypassing the atomistic resolution and its prohibitive computational costs, ECG enables the first calculations of the morphology dependence of the electronic structure of charge carriers across LC phases at the ∼20 nm length scale, with robust statistical sampling. Kinetic Monte Carlo (kMC) simulations reveal a strong morphology dependence on zero-field charge mobility among different LC phases as well as the presence of two-molecule charge carriers that act as traps and hinder charge transport. We leverage these results to further evaluate the feasibility of developing mesoscopic, field-based ECG models in future works. The fully CG approach to electronic property predictions in LC semiconductors opens a new computational direction for designing electronic processes in soft materials at their characteristic length scales.

Assessing molecular doping efficiency in organic semiconductors with reactive Monte Carlo.

The addition of molecular dopants into organic semiconductors (OSCs) is a ubiquitous augmentation strategy to enhance the electrical conductivity of OSCs. Although the importance of optimizing OSC-dopant interactions is well-recognized, chemically generalizable structure-function relationships are difficult to extract due to the sensitivity and dependence of doping efficiency on chemistry, processing conditions, and morphology. Computational modeling for an integrated OSC-dopant design is an attractive approach to systematically isolate fundamental relationships, but requires the challenging simultaneous treatment of molecular reactivity and morphology evolution. We present the first computational study to couple molecular reactivity with morphology evolution in a molecularly doped OSC. Reactive Monte Carlo is employed to examine the evolution of OSC-dopant morphologies and doping efficiency with respect to dielectric, the thermodynamic driving for the doping reaction, and dopant aggregation. We observe that for well-mixed systems with experimentally relevant dielectric constants, doping efficiency is near unity with a very weak dependence on the ionization potential and electron affinity of OSC and dopant, respectively. At experimental dielectric constants, reaction-induced aggregation is observed, corresponding to the well-known insolubility of solution-doped materials. Simulations are qualitatively consistent with a number of experimental studies showing a decrease of doping efficiency with increasing dopant concentration. Finally, we observe that the aggregation of dopants lowers doping efficiency and thus presents a rational design strategy for maximizing doping efficiency in molecularly doped OSCs. This work represents an important first step toward the systematic integration of molecular reactivity and morphology evolution into the characterization of multi-scale structure-function relationships in molecularly doped OSCs.

Distilling coarse-grained representations of molecular electronic structure with continuously gated message passing.

Bottom-up methods for coarse-grained (CG) molecular modeling are critically needed to establish rigorous links between atomistic reference data and reduced molecular representations. For a target molecule, the ideal reduced CG representation is a function of both the conformational ensemble of the system and the target physical observable(s) to be reproduced at the CG resolution. However, there is an absence of algorithms for selecting CG representations of molecules from which complex properties, including molecular electronic structure, can be accurately modeled. We introduce continuously gated message passing (CGMP), a graph neural network (GNN) method for atomically decomposing molecular electronic structure sampled over conformational ensembles. CGMP integrates 3D-invariant GNNs and a novel gated message passing system to continuously reduce the atomic degrees of freedom accessible for electronic predictions, resulting in a one-shot importance ranking of atoms contributing to a target molecular property. Moreover, CGMP provides the first approach by which to quantify the degeneracy of "good" CG representations conditioned on specific prediction targets, facilitating the development of more transferable CG representations. We further show how CGMP can be used to highlight multiatom correlations, illuminating a path to developing CG electronic Hamiltonians in terms of interpretable collective variables for arbitrarily complex molecules.

Excited State Dynamics of a Conformationally Fluxional Copper Coordination Complex.

The conversion of solar energy into chemical fuel represents a capstone goal of the 21st century and has the potential to supply terawatts of power in a globally distributed manner. However, the disparate time scales of photodriven charge separation (∼fs) and steps in chemical reactions (∼µs) represent an inherent bottleneck in solar-to-fuels technology. To address this discrepancy, we are developing earth-abundant coordination complexes that undergo light-induced conformational rearrangements such that charge separation (CS) is hastened, while charge recombination (CR) is slowed. To these ends, we report the preparation and characterization of a new series of conformationally fluxional copper coordination complexes that contain a twisted intramolecular charge transfer (TICT) fluorophore as part of their ligand scaffold. Structural and spectroscopic characterization of the Cu(I) and Cu(II) complexes formed with these ligands in their ground states establish oxidation state-dependent conformational dynamicity, while time-resolved emission and transient absorption spectroscopies define the photophysical parameters of photo-induced excited states. Building on initial reports with a related set of molecules, the improved ligand design presented here greatly simplifies the observed photophysics, effectively shutting down unwanted ligand-centered excited states previously observed. Time-dependent density functional theory (TDDFT) analyses reveal an unusual metal-to-TICT electronic transition only reported once before, and though the formation of a CS state is not observed directly through experiments, TDDFT geometry optimizations in the excited states support the formation of transient Cu(II) CS species, lending credence to the potential success of our approach. These studies establish a clear model for the excited state dynamics at play in proof-of-concept systems and clarify key design parameters for future optimizations toward achieving long-lived CS via photoinduced conformational gating.

Open Macromolecular Genome: Generative Design of Synthetically Accessible Polymers.

A grand challenge in polymer science lies in the predictive design of new polymeric materials with targeted functionality. However, de novo design of functional polymers is challenging due to the vast chemical space and an incomplete understanding of structure-property relations. Recent advances in deep generative modeling have facilitated the efficient exploration of molecular design space, but data sparsity in polymer science is a major obstacle hindering progress. In this work, we introduce a vast polymer database known as the Open Macromolecular Genome (OMG), which contains synthesizable polymer chemistries compatible with known polymerization reactions and commercially available reactants selected for synthetic feasibility. The OMG is used in concert with a synthetically aware generative model known as Molecule Chef to identify property-optimized constitutional repeating units, constituent reactants, and reaction pathways of polymers, thereby advancing polymer design into the realm of synthetic relevance. As a proof-of-principle demonstration, we show that polymers with targeted octanol-water solubilities are readily generated together with monomer reactant building blocks and associated polymerization reactions. Suggested reactants are further integrated with Reaxys polymerization data to provide hypothetical reaction conditions (e.g., temperature, catalysts, and solvents). Broadly, the OMG is a polymer design approach capable of enabling data-intensive generative models for synthetic polymer design. Overall, this work represents a significant advance, enabling the property targeted design of synthetic polymers subject to practical synthetic constraints.

Identifying Coarse-Grained Representations for Electronic Predictions.

Coarse-grained (CG) simulations are an important computational tool in chemistry and materials science. Recently, systematic "bottom-up" CG models have been introduced to capture electronic structure variations of molecules and polymers at the CG resolution. However, the performance of these models is limited by the ability to select reduced representations that preserve electronic structure information, which remains a challenge. We propose two methods for (i) identifying important electronically coupled atomic degrees of freedom and (ii) scoring the efficacy of CG representations used in conjunction with CG electronic predictions. The first method is a physically motivated approach that incorporates nuclear vibrations and electronic structure derived from simple quantum chemical calculations. We complement this physically motivated approach with a machine learning technique based on the marginal contribution of nuclear degrees of freedom to electronic prediction accuracy using an equivariant graph neural network. By integrating these two approaches, we can both identify critical electronically coupled atomic coordinates and score the efficacy of arbitrary CG representations for making electronic predictions. We leverage this capability to make a connection between optimized CG representations and the future potential for "bottom-up" development of simplified model Hamiltonians incorporating nonlinear vibrational modes.

AL4GAP: Active learning workflow for generating DFT-SCAN accurate machine-learning potentials for combinatorial molten salt mixtures.

Machine learning interatomic potentials have emerged as a powerful tool for bypassing the spatiotemporal limitations of ab initio simulations, but major challenges remain in their efficient parameterization. We present AL4GAP, an ensemble active learning software workflow for generating multicomposition Gaussian approximation potentials (GAP) for arbitrary molten salt mixtures. The workflow capabilities include: (1) setting up user-defined combinatorial chemical spaces of charge neutral mixtures of arbitrary molten mixtures spanning 11 cations (Li, Na, K, Rb, Cs, Mg, Ca, Sr, Ba and two heavy species, Nd, and Th) and 4 anions (F, Cl, Br, and I), (2) configurational sampling using low-cost empirical parameterizations, (3) active learning for down-selecting configurational samples for single point density functional theory calculations at the level of Strongly Constrained and Appropriately Normed (SCAN) exchange-correlation functional, and (4) Bayesian optimization for hyperparameter tuning of two-body and many-body GAP models. We apply the AL4GAP workflow to showcase high throughput generation of five independent GAP models for multicomposition binary-mixture melts, each of increasing complexity with respect to charge valency and electronic structure, namely: LiCl-KCl, NaCl-CaCl2, KCl-NdCl3, CaCl2-NdCl3, and KCl-ThCl4. Our results indicate that GAP models can accurately predict structure for diverse molten salt mixture with density functional theory (DFT)-SCAN accuracy, capturing the intermediate range ordering characteristic of the multivalent cationic melts.

Bypassing backmapping: Coarse-grained electronic property distributions using heteroscedastic Gaussian processes.

We employ deep kernel learning electronic coarse-graining (DKL-ECG) with approximate Gaussian processes as a flexible and scalable framework for learning heteroscedastic electronic property distributions as a smooth function of coarse-grained (CG) configuration. The appropriateness of the Gaussian prior on predictive CG property distributions is justified as a function of CG model resolution by examining the statistics of target distributions. The certainties of predictive CG distributions are shown to be limited by CG model resolution with DKL-ECG predictive noise converging to the intrinsic physical noise induced by the CG mapping operator for multiple chemistries. Further analysis of the resolution dependence of learned CG property distributions allows for the identification of CG mapping operators that capture CG degrees of freedom with strong electron-phonon coupling. We further demonstrate the ability to construct the exact quantum chemical valence electronic density of states (EDOS), including behavior in the tails of the EDOS, from an entirely CG model by combining iterative Boltzmann inversion and DKL-ECG. DKL-ECG provides a means of learning CG distributions of all-atom properties that are traditionally "lost" in CG model development, introducing a promising methodological alternative to backmapping algorithms commonly employed to recover all-atom property distributions from CG simulations.

Efficient Intermolecular Charge Transport in π-Stacked Pyridinium Dimers Using Cucurbit[8]uril Supramolecular Complexes.

Intermolecular charge transport through π-conjugated molecules plays an essential role in biochemical redox processes and energy storage applications. In this work, we observe highly efficient intermolecular charge transport upon dimerization of pyridinium molecules in the cavity of a synthetic host (cucurbit[8]uril, CB[8]). Stable, homoternary complexes are formed between pyridinium molecules and CB[8] with high binding affinity, resulting in an offset stacked geometry of two pyridiniums inside the host cavity. The charge transport properties of free and dimerized pyridiniums are characterized using a scanning tunneling microscope-break junction (STM-BJ) technique. Our results show that π-stacked pyridinium dimers exhibit comparable molecular conductance to isolated, single pyridinium molecules, despite a longer transport pathway and a switch from intra- to intermolecular charge transport. Control experiments using a CB[8] homologue (cucurbit[7]uril, CB[7]) show that the synthetic host primarily serves to facilitate dimer formation and plays a minimal role on molecular conductance. Molecular modeling using density functional theory (DFT) reveals that pyridinium molecules are planarized upon dimerization inside the host cavity, which facilitates charge transport. In addition, the π-stacked pyridinium dimers possess large intermolecular LUMO-LUMO couplings, leading to enhanced intermolecular charge transport. Overall, this work demonstrates that supramolecular assembly can be used to control intermolecular charge transport in π-stacked molecules.

Coarse-Grained Density Functional Theory Predictions via Deep Kernel Learning.

Scalable electronic predictions are critical for soft materials design. Recently, the Electronic Coarse-Graining (ECG) method was introduced to renormalize all-atom quantum chemical (QC) predictions to coarse-grained (CG) resolutions using deep neural networks (DNNs). While DNNs can learn complex representations that prove challenging for kernel-based methods, they are susceptible to overfitting and the overconfidence of uncertainty estimations. Here, we develop ECG within a GPU-accelerated Deep Kernel Learning (DKL) framework to enable CG QC predictions using range-separated hybrid density functional theory (DFT), obtaining a 107 speedup relative to naive all-atom QC. By treating the predicted electronic properties as random Gaussian Processes, DKL incorporates CG mapping degeneracy by learning the distribution of electronic energies as a function of CG configuration. DKL-ECG accurately reproduces molecular orbital energies from range-separated DFT while facilitating efficient training via active learning using the uncertainties provided by DKL. We show that while active learning algorithms enable efficient sampling of a more diverse configurational space relative to random sampling, all explored query methods exhibit comparable performance for the examined system. We attribute this result to the significant overlap of the feature space and output property distributions across multiple temperatures.

Coarse-Graining Organic Semiconductors: The Path to Multiscale Design.

Four decades of molecular theory and computation have helped form the modern understanding of the physical chemistry of organic semiconductors. Whereas these efforts have historically centered around characterizations of electronic structure at the single-molecule or dimer scale, emerging trends in noncrystalline molecular and polymeric semiconductors are motivating the need for modeling techniques capable of morphological and electronic structure predictions at the mesoscale. Provided the challenges associated with these prediction tasks, the community has begun to evolve a computational toolkit for organic semiconductors incorporating techniques from the fields of soft matter, coarse-graining, and machine learning. Here, we highlight recent advances in coarse-grained methodologies aimed at the multiscale characterization of noncrystalline organic semiconductors. As organic semiconductor performance is dependent on the interplay of mesoscale morphology and molecular electronic structure, specific emphasis is placed on coarse-grained modeling approaches capable of both structural and electronic predictions without recourse to all-atom representations.

Targeted sequence design within the coarse-grained polymer genome.

The chemical design of polymers with target structural and/or functional properties represents a grand challenge in materials science. While data-driven design approaches are promising, success with polymers has been limited, largely due to limitations in data availability. Here, we demonstrate the targeted sequence design of single-chain structure in polymers by combining coarse-grained modeling, machine learning, and model optimization. Nearly 2000 unique coarse-grained polymers are simulated to construct and analyze machine learning models. We find that deep neural networks inexpensively and reliably predict structural properties with limited sequence information as input. By coupling trained ML models with sequential model-based optimization, polymer sequences are proposed to exhibit globular, swollen, or rod-like behaviors, which are verified by explicit simulations. This work highlights the promising integration of coarse-grained modeling with data-driven design and represents a necessary and crucial step toward more complex polymer design efforts.

Over What Length Scale Does an Inorganic Substrate Perturb the Structure of a Glassy Organic Semiconductor?

While the bulk structure of vapor-deposited glasses has been extensively studied, structure at buried interfaces has received little attention, despite being important for organic electronic applications. To learn about glass structure at buried interfaces, we study the structure of vapor-deposited glasses of the organic semiconductor DSA-Ph (1,4-di-[4-(N,N-diphenyl)amino]styrylbenzene) as a function of film thickness; the structure is probed with grazing incidence X-ray scattering. We deposit on silicon and gold substrates and span a film thickness range of 10-600 nm. Our experiments demonstrate that interfacial molecular packing in vapor-deposited glasses of DSA-Ph is more disordered compared to the bulk. At a deposition temperature near room temperature, we estimate ∼8 nm near the substrate can have modified molecular packing. Molecular dynamics simulations of a coarse-grained representation of DSA-Ph reveal a similar length scale. In both the simulations and the experiments, deposition temperature controls glass structure beyond this interfacial layer of a few nanometers.

Spatiotemporal Formation and Growth Kinetics of Polyelectrolyte Complex Micelles with Millisecond Resolution.

We have directly observed the in situ self-assembly kinetics of polyelectrolyte complex (PEC) micelles by synchrotron time-resolved small-angle X-ray scattering, equipped with a stopped-flow device that provides millisecond temporal resolution. A synthesized neutral-charged diblock polycation and homopolyanion that we have previously investigated as a model charge-matched, core-shell micelle system were selected for this work. The initial micellization of the oppositely charged polyelectrolytes was completed within the dead time of mixing of 100 ms, followed by micelle growth and equilibration up to several seconds. By combining the structural evolution of the radius of gyration (Rg) with complementary molecular dynamics simulations, we show how the self-assemblies evolve incrementally in size over time through a two-step kinetic process: first, oppositely charged polyelectrolyte chains pair to form nascent aggregates that immediately assemble into spherical micelles, and second, these PEC micelles grow into larger micellar entities. This work has determined one possible kinetic pathway for the initial formation of PEC micelles, which provides useful physical insights for increasing fundamental understanding self-assembly dynamics, driven by polyelectrolyte complexation that occurs on ultrafast time scales.

Crossover from Rouse to Reptation Dynamics in Salt-Free Polyelectrolyte Complex Coacervates.

Considerable interest in the dynamics and rheology of polyelectrolyte complex coacervates has been motivated by their industrial application as viscosity modifiers. A central question is the extent to which classical Rouse and reptation models can be applied to systems where electrostatic interactions play a critical role on the thermodynamics. By relying on molecular simulations, we present a direct analysis of the crossover from Rouse to reptation dynamics in salt-free complex coacervates as a function of chain length. This crossover shifts to shorter chain lengths as electrostatic interactions become stronger, which corresponds to the formation of denser coacervates. To distinguish the roles of Coulomb interactions and density, we compare the dynamics of coacervates to those of neutral, semidilute solutions at the same density. Both systems exhibit a universal dynamical behavior in the connectivity-dominated (subdiffusion and normal diffusion) regimes, but the monomer relaxation time in coacervates is much longer and increases with increasing Bjerrum length. This is similar to the cage effect observed in glass-forming polymers, but the local dynamical slowdown is caused here by strong Coulomb attractions (ion pairing) between oppositely charged monomers. Our findings provide a microscopic framework for the quantitative understanding of coacervate dynamics and rheology.

Electronic structure at coarse-grained resolutions from supervised machine learning.

Computational studies aimed at understanding conformationally dependent electronic structure in soft materials require a combination of classical and quantum-mechanical simulations, for which the sampling of conformational space can be particularly demanding. Coarse-grained (CG) models provide a means of accessing relevant time scales, but CG configurations must be back-mapped into atomistic representations to perform quantum-chemical calculations, which is computationally intensive and inconsistent with the spatial resolution of the CG models. A machine learning approach, denoted as artificial neural network electronic coarse graining (ANN-ECG), is presented here in which the conformationally dependent electronic structure of a molecule is mapped directly to CG pseudo-atom configurations. By averaging over decimated degrees of freedom, ANN-ECG accelerates simulations by eliminating backmapping and repeated quantum-chemical calculations. The approach is accurate, consistent with the CG spatial resolution, and can be used to identify computationally optimal CG resolutions.

Controlling Complex Coacervation via Random Polyelectrolyte Sequences.

The utilization of chemical sequence control in polymeric materials is key to enabling material design on par with biomacromolecular systems. One important avenue for scalable sequence-controlled polymers leverages the random copolymerization of distinct monomers, with the statistical distribution of the monomeric sequence arising from reaction kinetics following a first-order Markov process. Here we utilize the framework of the random phase approximation (RPA) to develop a theory for the phase behavior of symmetric polyelectrolyte coacervates whose chemical sequences are dictated by simple statistical distributions. We find that a high charge "blockiness" within the random sequences favors the formation of denser and more salt-resistant coacervates while simultaneously increasing the width of the two-phase region. We trace these physical effects to the increased cooperativity of Coulomb interactions that results from increased charge blockiness in oppositely charged polyelectrolytes.

Origin of Anisotropic Molecular Packing in Vapor-Deposited Alq3 Glasses.

Anisotropic molecular packing is a key feature that makes glasses prepared by physical vapor deposition (PVD) unique materials, warranting a mechanistic understanding of how a PVD glass attains its structure. To this end, we use X-ray scattering and ellipsometry to characterize the structure of PVD glasses of tris(8-hydroxyquinoline) aluminum (Alq3), a molecule used in organic electronics, and compare our results to simulations of its supercooled liquid. X-ray scattering reveals a tendency for molecular layering in Alq3 glasses that depends upon the substrate temperature during deposition and the deposition rate. Simulations reveal that the Alq3 supercooled liquid, like liquid metals, exhibits surface layering. We propose that the layering in Alq3 glasses observed here as well as the previously reported bulk dipole orientation are inherited from the surface structure of the supercooled liquid. This work significantly advances our understanding of the mechanism governing the formation of anisotropic structure in PVD glasses.