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.

Photophysical Outcomes of Water-Solvated Heterocycles: Single-Conformation Ultraviolet and Infrared Spectroscopy of Microsolvated 2-Phenylpyrrole.

The molecular chromophores within brown carbon (BrC) aerosols absorb solar radiation at visible and near-ultraviolet wavelengths. This contributes to the overall warming of the troposphere and the photochemical aging of aerosols. In this investigation, we combine a suite of experimental and theoretical methods to reveal the conformation-specific ultraviolet and infrared spectroscopy of 2-phenylpyrrole (2PhPy)─an extended π-conjugated pyrrole derivative and a model BrC chromophore─along with its water microsolvated molecular complexes (2PhPy:nH2O, n = 1-3). Using resonant two-photon ionization and double-resonance holeburning techniques alongside MP3 (ground state) and ADC(3) (excited state) torsional potential energy surfaces and discrete variable representation simulations, we characterized the ultraviolet spectra of 2PhPy and 2PhPy:1H2O. This analysis revealed evidence for Herzberg-Teller vibronic coupling along the CH wagging and NH stretching coordinates of the aromatic rings. Conformation-specific infrared spectroscopy revealed extended hydrogen-bonding networks of the 2PhPy:nH2O complexes. Upon stepwise addition of H2O solvation, the nearest H2O acceptor forms a strong, noncovalent interaction with the pyrrole NH donor, while the second and third H2O partners interface with the phenyl and pyrrole aromatic rings through growing van der Waals π/H atom stabilization. A local-mode Hamiltonian approach was employed for comparison with the experimental spectra, thus identifying the vibrational spectral signatures to specific 2PhPy:nH2O oscillators.

Generative organic electronic molecular design informed by quantum chemistry.

Generative molecular design strategies have emerged as promising alternatives to trial-and-error approaches for exploring and optimizing within large chemical spaces. To date, generative models with reinforcement learning approaches have frequently used low-cost methods to evaluate the quality of the generated molecules, enabling many loops through the generative model. However, for functional molecular materials tasks, such low-cost methods are either not available or would require the generation of large amounts of training data to train surrogate machine learning models. In this work, we develop a framework that connects the REINVENT reinforcement learning framework with excited state quantum chemistry calculations to discover molecules with specified molecular excited state energy levels, specifically molecules with excited state landscapes that would serve as promising singlet fission or triplet-triplet annihilation materials. We employ a two-step curriculum strategy to first find a set of diverse promising molecules, then demonstrate the framework's ability to exploit a more focused chemical space with anthracene derivatives. Under this protocol, we show that the framework can find desired molecules and improve Pareto fronts for targeted properties versus synthesizability. Moreover, we are able to find several different design principles used by chemists for the design of singlet fission and triplet-triplet annihilation molecules.

Reorganization Energy Predictions with Graph Neural Networks Informed by Low-Cost Conformers.

A critical bottleneck for the design of high-conductivity organic materials is finding molecules with low reorganization energy. To enable high-throughput virtual screening campaigns for many types of organic electronic materials, a fast reorganization energy prediction method compared to density functional theory is needed. However, the development of low-cost machine-learning-based models for calculating the reorganization energy has proven to be challenging. In this paper, we combine a 3D graph-based neural network (GNN) recently benchmarked for drug design applications, ChIRo, with low-cost conformational features for reorganization energy predictions. By comparing the performance of ChIRo to another 3D GNN, SchNet, we find evidence that the bond-invariant property of ChIRo enables the model to learn from low-cost conformational features more efficiently. Through an ablation study with a 2D GNN, we find that using low-cost conformational features on top of 2D features informs the model for making more accurate predictions. Our results demonstrate the feasibility of reorganization energy predictions on the benchmark QM9 data set without needing DFT-optimized geometries and demonstrate the types of features needed for robust models that work on diverse chemical spaces. Furthermore, we show that ChIRo informed with low-cost conformational features achieves comparable performance with the previously reported structure-based model on π-conjugated hydrocarbon molecules. We expect this class of methods can be applied to the high-throughput screening of high-conductivity organic electronics candidates.

The role of the electrolyte in non-conjugated radical polymers for metal-free aqueous energy storage electrodes.

Metal-free aqueous batteries can potentially address the projected shortages of strategic metals and safety issues found in lithium-ion batteries. More specifically, redox-active non-conjugated radical polymers are promising candidates for metal-free aqueous batteries because of the polymers' high discharge voltage and fast redox kinetics. However, little is known regarding the energy storage mechanism of these polymers in an aqueous environment. The reaction itself is complex and difficult to resolve because of the simultaneous transfer of electrons, ions and water molecules. Here we demonstrate the nature of the redox reaction for poly(2,2,6,6-tetramethylpiperidinyloxy-4-yl acrylamide) by examining aqueous electrolytes of varying chao-/kosmotropic character using electrochemical quartz crystal microbalance with dissipation monitoring at a range of timescales. Surprisingly, the capacity can vary by as much as 1,000% depending on the electrolyte, in which certain ions enable better kinetics, higher capacity and higher cycling stability.

Data-Driven Models for Predicting Intrinsically Disordered Protein Polymer Physics Directly from Composition or Sequence.

The molecular-level understanding of intrinsically disordered proteins is challenging due to experimental characterization difficulties. Computational understanding of IDPs also requires fundamental advances, as the leading tools for predicting protein folding (e.g., AlphaFold), typically fail to describe the structural ensembles of IDPs. The focus of this paper is to 1) develop new representations for intrinsically disordered proteins and 2) pair these representations with classical machine learning and deep learning models to predict the radius of gyration and derived scaling exponent of IDPs. Here, we build a new physically-motivated feature called the bag of amino acid interactions representation, which encodes pairwise interactions explicitly into the representation. This feature essentially counts and weights all possible non-bonded interactions in a sequence and thus is, in principle, compatible with arbitrary sequence lengths. To see how well this new feature performs, both categorical and physically-motivated featurization techniques are tested on a computational dataset containing 10,000 sequences simulated at the coarse-grained level. The results indicate that this new feature outperforms the other purely categorical and physically-motivated features and possesses solid extrapolation capabilities. For future use, this feature can potentially provide physical insights into amino acid interactions, including their temperature dependence, and be applied to other protein spaces.

Charge Transfer in Spatially Defined Organic Radical Polymers.

Charge transfer in nonconjugated redox-active polymers is influenced by redox site proximity and polymer flexibility, but it is challenging to observe these effects independently. In this work, spatially defined radical-containing polymers are synthesized by using acyclic diene metathesis (ADMET) polymerization of α,ω-dienes bearing a central activated ester. Postpolymerization functionalization with 4-amino-2,2,6,6-tetramethylpiperidine-1-oxyl (4-amino-TEMPO) introduces TEMPO radical groups onto the polymer backbone through amide linkages to yield spatially defined polymers with radical units every 9, 11, 15, and 21 carbons. Increased radical spacing leads to reduced spin-spin coupling and increased chain flexibility. The glass transition temperatures (Tg) range from 47.6 to -13.8 °C, depending on the radical spacing. The spatially defined TEMPO-substituted polymer with a spacing length of 15 carbons displays the lowest Tg and the shortest hopping distance, as shown through molecular dynamics simulations. Also, this polymer displays kinetics 1000 times faster than the commonly studied TEMPO-containing polymer poly(2,2,6,6-tetramethylpiperidinyloxy-4-ylacrylamide) (PTAm). Remarkably, comparison of the diffusion and kinetics attributed to the redox reaction reveals that both the apparent diffusion coefficient and the self-exchange reaction rate constant are correlated to the polymer's Tg as log[Dapp] and log[kex,app] ∼ Tg, respectively. Critically, these data demonstrate that controlling the spacing of redox-active groups along a polymer backbone strongly influences backbone flexibility and radical packing, which leads to synergetic improvements in the charge transfer kinetics of nonconjugated redox-active polymers.

Solvent-Mediated Charge Transfer Dynamics of a Model Brown Carbon Aerosol Chromophore: Photophysics of 1-Phenylpyrrole Induced by Water Solvation.

Nitrogen heterocycles are known to be important light-absorbing chromophores in a newly discovered class of aerosols, commonly referred to as "brown carbon" (BrC) aerosols. Due to their significant absorption and spectral overlap with the solar actinic flux, these BrC chromophores steer the physical and optical properties of aerosols. To model the local aqueous solvation environment surrounding BrC chromophores, we generated cold molecular complexes with water and a prototypical BrC chromophore, 1-phenylpyrrole (1PhPy), using supersonic jet-cooling and explored their intermolecular interactions using single-conformation spectroscopy. Herein, we utilized resonant two-photon ionization (R2PI) and UV holeburning (UV HB) double-resonance spectroscopies to obtain a molecular-level understanding of the role of water microsolvation in charge transfer upon photoexcitation of 1PhPy. Quantum chemical calculations and one-dimensional discrete variable representation simulations revealed insights into the charge transfer efficacy of 1PhPy with and without addition of a single water molecule. Taken together, our results indicate that the intermolecular interactions with water guide the geometry of 1PhPy to adopt a more twisted intramolecular charge transfer (TICT) configuration, thus facilitating charge transfer from the pyrrole donor to the phenyl ring acceptor. Furthermore, the water network surrounding 1PhPy reports on the charge transfer such that the H2O solvent primarily interacts with the pyrrole ring donor in the ground state, whereas it preferentially interacts with the phenyl ring acceptor in the excited state. Large Franck-Condon activity is evident in the 1PhPy + 1H2O excitation spectrum for the water-migration vibronic bands, supporting H2O solvent reorganization upon excitation of the 1PhPy chromophore. Fluorescence measurements with increasing H2O % volume corroborated our gas-phase studies by indicating that a polar water solvation environment stabilizes the TICT configuration of 1PhPy in the excited electronic state, from which emission is observed at a lower energy compared to the locally excited configuration.

Interface chemistry of an amide electrolyte for highly reversible lithium metal batteries.

Metallic lithium is a promising anode to increase the energy density of rechargeable lithium batteries. Despite extensive efforts, detrimental reactivity of lithium metal with electrolytes and uncontrolled dendrite growth remain challenging interconnected issues hindering highly reversible Li-metal batteries. Herein, we report a rationally designed amide-based electrolyte based on the desired interface products. This amide electrolyte achieves a high average Coulombic efficiency during cycling, resulting in an outstanding capacity retention with a 3.5 mAh cm-2 high-mass-loaded LiNi0.8Co0.1Mn0.1O2 cathode. The interface reactions with the amide electrolyte lead to the predicted solid electrolyte interface species, having favorable properties such as high ionic conductivity and high stability. Operando monitoring the lithium spatial distribution reveals that the highly reversible behavior is related to denser deposition as well as top-down stripping, which decreases the formation of porous deposits and inactive lithium, providing new insights for the development of interface chemistries for metal batteries.

Design rules for high mobility xanthene-based hole transport materials.

Tunable and highly conductive hole transport materials are crucial for the performance of organic electronics applications such as organic light emitting diodes and perovskite solar cells. For commercial applications, these materials' requirements include easy synthesis, high hole mobility, and highly tuned and compatible electronic energy levels. Here, we present a systematic study of a recently discovered, easy-to-synthesize class of spiro[fluorene-9,9'-xanthene]-based organic hole transport materials. Systematic side group functionalization allows us to control the HOMO energy and charge carrier mobility. Analysis of the bulk simulations enables us to derive design rules for mobility enhancement. We show that larger functional groups (e.g. methyl) decrease the conformational disorder due to steric effects and thus increase the hole mobility. Highly asymmetric or polar side groups (e.g. fluorine), however, increase the electrostatic disorder and thus reduce the hole mobility. These generally applicable design rules will help in the future to further optimize organic hole transport materials.

Discovery of blue singlet exciton fission molecules via a high-throughput virtual screening and experimental approach.

Singlet exciton fission is a mechanism that could potentially enable solar cells to surpass the Shockley-Queisser efficiency limit by converting single high-energy photons into two lower-energy triplet excitons with minimal thermalization loss. The ability to make use of singlet exciton fission to enhance solar cell efficiencies has been limited, however, by the sparsity of singlet fission materials with triplet energies above the bandgaps of common semiconductors such as Si and GaAs. Here, we employ a high-throughput virtual screening procedure to discover new organic singlet exciton fission candidate materials with high-energy (>1.4 eV) triplet excitons. After exploring a search space of 4482 molecules and screening them using time-dependent density functional theory, we identify 88 novel singlet exciton fission candidate materials based on anthracene derivatives. Subsequent purification and characterization of several of these candidates yield two new singlet exciton fission materials: 9,10-dicyanoanthracene (DCA) and 9,10-dichlorooctafluoroanthracene (DCOFA), with triplet energies of 1.54 eV and 1.51 eV, respectively. These materials are readily available and low-cost, making them interesting candidates for exothermic singlet exciton fission sensitization of solar cells. However, formation of triplet excitons in DCA and DCOFA is found to occur via hot singlet exciton fission with excitation energies above ∼3.64 eV, and prominent excimer formation in the solid state will need to be overcome in order to make DCA and DCOFA viable candidates for use in a practical device.

Extending the Lifetime of Organic Flow Batteries via Redox State Management.

Redox flow batteries based on quinone-bearing aqueous electrolytes have emerged as promising systems for energy storage from intermittent renewable sources. The lifetime of these batteries is limited by quinone stability. Here, we confirm that 2,6-dihydroxyanthrahydroquinone tends to form an anthrone intermediate that is vulnerable to subsequent irreversible dimerization. We demonstrate quantitatively that this decomposition pathway is responsible for the loss of battery capacity. Computational studies indicate that the driving force for anthrone formation is greater for anthraquinones with lower reduction potentials. We show that the decomposition can be substantially mitigated. We demonstrate that conditions minimizing anthrone formation and avoiding anthrone dimerization slow the capacity loss rate by over an order of magnitude. We anticipate that this mitigation strategy readily extends to other anthraquinone-based flow batteries and is thus an important step toward realizing renewable electricity storage through long-lived organic flow batteries.

Precise Control of Thermal and Redox Properties of Organic Hole-Transport Materials.

We report design principles of the thermal and redox properties of synthetically accessible spiro-based hole transport materials (HTMs) and show the relevance of these findings to high-performance perovskite solar cells (PSCs). The chemical modification of an asymmetric spiro[fluorene-9,9'-xanthene] core is amenable to selective placement of redox active triphenylamine (TPA) units. We therefore leveraged computational techniques to investigate five HTMs bearing TPA groups judiciously positioned about this asymmetric spiro core. It was determined that TPA groups positioned about the conjugated fluorene moiety increase the free energy change for hole-extraction from the perovskite layer, while TPAs about the xanthene unit govern the Tg values. The synergistic effects of these characteristics resulted in an HTM characterized by both a low reduction potential (≈0.7 V vs. NHE) and a high Tg value (>125 °C) to yield a device power conversion efficiency (PCE) of 20.8 % in a PSC.

Excitonics: A Set of Gates for Molecular Exciton Processing and Signaling.

Regulating energy transfer pathways through materials is a central goal of nanotechnology, as a greater degree of control is crucial for developing sensing, spectroscopy, microscopy, and computing applications. Such control necessitates a toolbox of actuation methods that can direct energy transfer based on user input. Here we introduce a proposal for a molecular exciton gate, analogous to a traditional transistor, for regulating exciton flow in chromophoric systems. The gate may be activated with an input of light or an input flow of excitons. Our proposal relies on excitation migration via the second excited singlet (S2) state of the gate molecule. It exhibits the following features, only a subset of which are present in previous exciton switching schemes: picosecond time scale actuation, amplification/gain behavior, and a lack of molecular rearrangement. We demonstrate that the device can be used to produce universal binary logic or amplification of an exciton current, providing an excitonic platform with several potential uses, including signal processing for microscopy and spectroscopy methods that implement tunable exciton flux.

Infrared-Enhanced Fluorescence-Gain Spectroscopy: Conformation-Specific Excited-State Infrared Spectra of Alkylbenzenes.

An ultraviolet-infrared (UV-IR) double-resonance method for recording conformation-specific excited-state infrared spectra is described. The method takes advantage of an increase in fluorescence signal in phenylalkanes produced by infrared excitation of the S1 origin levels of different conformational isomers. The shorter lifetimes of these IR-excited molecules, combined with their red-shifted emission, provides a way to discriminate the fluorescence due to the infrared-excited molecules from the S1 origin fluorescence, resulting in spectra with high signal-to-noise ratios. Spectra for a series of phenylalkanes and a capped phenylalanine derivative (Ac-Phe-NHMe) demonstrate the potential of the method. The excited-state spectrum in the alkyl CH stretch region of ethylbenzene is well-fit by an anharmonic model developed for the ground electronic state, which explicitly takes into account stretch-bend Fermi resonance.

Fingerprints of inter- and intramolecular hydrogen bonding in saligenin-water clusters revealed by mid- and far-infrared spectroscopy.

Saligenin (2-(hydroxymethyl)phenol) exhibits both strong and weak intramolecular electrostatic interactions. The bonds that result from these interactions compete with intermolecular hydrogen bonds once saligenin binds to one or more water molecules. Infrared (IR) ultraviolet (UV) ion-dip spectroscopy was used to study isolated saligenin-(H2O)n clusters (n = 1-3) in the far- and mid-IR regions of the spectrum. Both harmonic and anharmonic (coupled local modes and Born-Oppenheimer molecular dynamics) quantum chemical calculations were applied to assign cluster geometries to the measured spectra, and to assign vibrational modes to all spectral features measured for each cluster. The hydrated clusters with n = 1 and 2 have geometries that are quite similar to benzyl alcohol-water clusters, whereas the larger clusters with n = 3 show structures equivalent to the isolated water pentamer. Systematic shifts in the frequencies of three hydrogen bond (H-bond) deforming modes, namely OH stretching, OH torsion and H-bond stretching, were studied as a function of the hydrogen bond strength represented by either the OH bond length or the H-bond length. The shifts of the frequencies of these three modes correlate linearly to the OH length, despite both intra- and intermolecular H-bonds being included in this analysis. The OH torsion vibration displays the largest frequency shift when H-bonded, followed by the OH stretching vibrations and finally the H-bond stretching frequency. The frequency shifts of these H-bond deforming modes behave non-linearly as a function of the H-bond length, asymptotically approaching the frequency expected for the non H-bonded modes. The nonlinear behavior was quantified using exponential functions.

Identifying the first folded alkylbenzene via ultraviolet, infrared, and Raman spectroscopy of pentylbenzene through decylbenzene.

The conformational preferences of pentyl- through decylbenzene are studied under jet-cooled conditions in the gas phase. Laser-induced fluorescence excitation spectra, fluorescence-dip infrared spectra in the alkyl CH stretch region, and Raman spectra are combined to provide assignments for the observed conformers. Density functional theory calculations at the B3LYP-D3BJ/def2TZVP level of theory provide relative energies and normal mode vibrations that serve as inputs for an anharmonic local mode theory introduced in earlier work on alkylbenzenes with n = 2-4. This model explicitly includes anharmonic mixing of the CH stretch modes with the overtones of scissors/bend modes of the CH2 and CH3 groups in the alkyl chain, and is used to assign and interpret the single-conformation IR spectra. In octylbenzene, a pair of LIF transitions shifted -92 and -78 cm-1 from the all-trans electronic origin have unique alkyl CH stretch transitions that are fit by the local model to a g1g3g4 conformation in which the alkyl chain folds back over the aromatic ring π cloud. Its calculated energy is only 1.0 kJ mol-1 above the all-trans global minimum. This fold is at an alkyl chain length less than half that of the pure alkanes (n = 18), consistent with a smaller energy cost for the g1 dihedral and the increased dispersive interaction of the chain with the π cloud. Local site frequencies for the entire set of conformers from the local mode model show 'edge effects' that raise the site frequencies of CH2(1) and CH2(2) due to the phenyl ring and CH2(n - 1) due to the methyl group. The g1g3g4 conformer also shows local sites shifted up in frequency at CH2(3) and CH2(6) due to interaction with the π cloud.

Infrared laser spectroscopy of the n-propyl and i-propyl radicals: Stretch-bend Fermi coupling in the alkyl CH stretch region.

The n-propyl and i-propyl radicals were generated in the gas phase via pyrolysis of n-butyl nitrite [CH3(CH2)3ONO] and i-butyl nitrite [(CH3)2CHCH2ONO], respectively. Nascent radicals were promptly solvated by a beam of He nanodroplets, and the infrared spectra of the radicals were recorded in the CH stretching region. Several previously unreported bands are observed between 2800 and 3150 cm-1. The CH stretching modes observed above 3000 cm-1 are in excellent agreement with CCSD(T) anharmonic frequencies computed using second-order vibrational perturbation theory. However, between 2800 and 3000 cm-1, the spectra of n- and i-propyl radicals become congested and difficult to assign due to the presence of multiple anharmonic resonance polyads. To model the spectrally congested region, Fermi and Darling-Dennison resonances are treated explicitly using "dressed" Hamiltonians and CCSD(T) quartic force fields in the normal mode representation, and the agreement with experiment is less than satisfactory. Computations employing local mode effective Hamiltonians reveal the origin of the spectral congestion to be strong coupling between the high frequency CH stretching modes and the lower frequency CHn bending/scissoring motions. The most significant coupling is between stretches and bends localized on the same CH2/CH3 group. Spectral simulations using the local mode approach are in excellent agreement with experiment.

Conformation-specific spectroscopy of alkyl benzyl radicals: Effects of a radical center on the CH stretch infrared spectrum of an alkyl chain.

An important initial step in the combustion of gasoline and diesel fuels is the abstraction of hydrogen from alkylbenzenes to form resonance-stabilized alkyl benzyl radicals. This work uses, for the first time, double resonance spectroscopy methods to explore the conformation-specific vibronic and infrared spectroscopy of the α-ethylbenzyl (αEtBz) and α-propylbenzyl (αPrBz) radicals. Local mode Hamiltonian modeling enables assignment of the alkyl CH stretch IR spectra, accounting for Fermi resonance that complicates aliphatic alkyl CH stretch IR spectroscopy. The ground state conformational preferences of the ethyl and propyl chains are changed from those in the alkylbenzenes themselves, with global minima occurring for an in-plane orientation of the alkyl chain (trans) about its first dihedral angle (ϕf123, numbers are alkyl C atoms. C1 is CH radical site). This in-plane structure is the only observed conformer for the α-EtBz radical, while two conformers, tt and tg' share this orientation at the first dihedral, but differ in the second (ϕ1234) for the αPrBz radical. The in-plane orientation lowers the local site frequencies of the CH2 group stretches immediately adjacent to the benzylic radical site by about 50 cm-1 relative to those in pure alkyl chains or alkylbenzenes. This effect of the radical site is localized on the first CH2 group, with little effect on subsequent members of the alkyl chain. In the D1 excited electronic state, an out-of-plane orientation is preferred for the alkyl chains, leading to torsional mode Franck-Condon activity in the D0-D1 spectra that is both conformer-specific and diagnostic of the conformational change.