Novel WRN Helicase Inhibitors Selectively Target Microsatellite Unstable Cancer Cells

Microsatellite-unstable (MSI) cancers require WRN helicase to resolve replication stress due to expanded DNA (TA)n-dinucleotide repeats. WRN is a promising synthetic lethal target for MSI tumours, and WRN inhibitors are in development. Here, we used CRISPR-Cas9 base editing to map WRN residues critical for MSI cells, validating the helicase domain as the primary drug target. Fragment-based screening led to the development of potent and highly selective WRN helicase covalent inhibitors. These compounds selectively suppressed MSI model growth In vitro and In vivo by mimicking WRN loss, inducing DNA double-strand breaks at expanded TA-repeats and DNA damage. Assessment of biomarkers in preclinical models linked TA-repeat expansions and mismatch repair (MMR) alterations to compound activity. Efficacy was confirmed in immunotherapy-resistant organoids and patient-derived xenograft (PDX) models. The discovery of potent, selective covalent WRN inhibitors provides proof of concept for synthetic-lethal targeting of WRN in MSI cancer and tools to dissect WRN biology.

Partial Molar Solvation Volume of the Hydrated Electron Simulated Via DFT

Different simulation models of the hydrated electron produce different solvation structures, but it has been challenging to determine which simulated solvation structure, if any, is the most comparable to experiment. In a recent work, Neupane et al. [J. Phys. Chem. B 2023, 127, 5941-5947] showed using Kirkwood-Buff theory that the partial molar volume of the hydrated electron, which is known experimentally, can be readily computed from an integral over the simulated electron-water radial distribution function. This provides a sensitive way to directly compare the hydration structure of different simulation models of the hydrated electron with experiment. Here, we compute the partial molar volume of an ab-initio-simulated hydrated electron model based on density-functional theory (DFT) with a hybrid functional at different simulated system sizes. We find that the partial molar volume of the DFT-simulated hydrated electron is not converged with respect to the system size for simulations with up to 128 waters. We show that even at the largest simulation sizes, the partial molar volume of DFT-simulated hydrated electrons is underestimated by a factor of 2 with respect to experiment, and at the standard 64-water size commonly used in the literature, DFT-based simulations underestimate the experimental solvation volume by a factor of ∼3.5. An extrapolation to larger box sizes does predict the experimental partial molar volume correctly; however, larger system sizes than those explored here are currently intractable without the use of machine-learned potentials. These results bring into question what aspects of the predicted hydrated electron radial distribution function, as calculated by DFT-based simulations with the PBEh-D3 functional, deviate from the true solvation structure.

Using Machine Learning to Understand the Causes of Quantum Decoherence in Solution-Phase Bond-Breaking Reactions

Decoherence is a fundamental phenomenon that occurs when an entangled quantum state interacts with its environment, leading to collapse of the wave function. The inevitability of decoherence provides one of the most intrinsic limits of quantum computing. However, there has been little study of the precise chemical motions from the environment that cause decoherence. Here, we use quantum molecular dynamics simulations to explore the photodissociation of Na2+ in liquid Ar, in which solvent fluctuations induce decoherence and thus determine the products of chemical bond breaking. We use machine learning to characterize the solute-solvent environment as a high-dimensional feature space that allows us to predict when and onto which photofragment the bonding electron will localize. We find that reaching a requisite photofragment separation and experiencing out-of-phase solvent collisions underlie decoherence during chemical bond breaking. Our work highlights the utility of machine learning for interpreting complex solution-phase chemical processes as well as identifies the molecular underpinnings of decoherence.

Does the Traditional Band Picture Describe the Electronic Structure of Doped Conjugated Polymers? TD-DFT and Natural Transition Orbital Study of Doped P3HT

Polarons and bipolarons are created when one or two electrons are removed from the π-system of a p-type conjugated polymer, respectively. In the traditional band picture, the creation of a polaron causes two electronic energy levels to move into the band gap. The removal of a second electron to form a bipolaron causes the two intragap states to move further into the gap. Several groups, however, who looked at the energies of the Kohn-Sham orbitals from DFT calculations, have recently argued that the traditional band picture is incorrect for explaining the spectroscopy of doped conjugated polymers. Instead, the DFT calculations suggest that polaron creation causes only one unoccupied state to move into the band gap near the valence band edge while half-filled state in the valence band and the conduction band bend downward in energy. To understand the discrepancy, we performed TD-DFT calculations of polarons and bipolarons on poly(3-hexylthiophene) (P3HT). Not only do the TD-DFT-calculated absorption spectra match the experimental absorption spectra, but an analysis using natural transitional orbitals (NTOs), which provides an approximate one-electron picture from the many-electron TD-DFT results, supports the traditional band picture. Our TD-DFT/NTO analysis indicates that the traditional band picture also works for bipolarons, a system for which DFT calculations were unable to determine the electronic structure.

Facile synthesis of 2-aza-9,10-diphenylanthracene and the effect of precise nitrogen atom incorporation on OLED emitters performance

Polycyclic aromatic hydrocarbons (PAHs) are important compounds in materials chemistry, particularly for optoelectronic applications. One strategy for tuning PAH properties involves the net exchange of carbon atoms for heteroatoms, such as nitrogen. We report a comparative study of the well-known fluorophore 9,10-diphenylanthracene with an aza analog. The latter compound is accessed using a short sequence involving the use of two strained cyclic alkynes, benzyne and a 3,4-piperidyne, in Diels-Alder cycloaddition sequences. Comparative studies of 9,10-diphenylanthracene and the aza-analog show how the addition of a single nitrogen atom impacts electrochemical and optical properties. Organic light-emitting diode (OLED) devices were prepared using both compounds, which showed that nitrogen substitution leads to an unexpected red shift in electroluminescence, likely due to exciplex formation between the active layer and the 4,4'-N,N'-bis[N-(1-naphthyl)-N-phenylamino]biphenyl (NPB) hole-transport layer. These studies highlight a unique approach to accessing heteroatom-containing PAHs, while underscoring the impact of heteroatoms on OLED device performance.

How Ions Break Local Symmetry: Simulations of Polarized Transient Hole Burning for Different Models of the Hydrated Electron in Contact Pairs with Na. J Phys Chem Lett 2023;14:3014-3022. [PMID: 36943261 DOI: 10.1021/acs.jpclett.3c00220]

The hydrated electron (e_aq^-) is known via polarized transient hole-burning (pTHB) experiments to have a homogeneously broadened absorption spectrum. Here, we explore via quantum simulation how the pTHB spectroscopy of different e_aq^- models changes in the presence of electrolytes. The idea is that cation-e_aq^- pairing can break the local symmetry and, thus, induce persistent inhomogeneity. We find that a "hard" cavity model shows a modest increase in the pTHB recovery time in the presence of salt, while a "soft" cavity model remains homogeneously broadened independent of the salt concentration. We also explore the orientational anisotropy of a fully ab initio density functional theory-based model of the e_aq^-, which is strongly inhomogeneously broadened without salt and which becomes significantly more inhomogeneously broadened in the presence of salt. The results provide a direct prediction for experiments that can distinguish between different models and, thus, help pin down the hydration structure and dynamics of the e_aq^-.

Ab Initio Studies of Hydrated Electron/Cation Contact Pairs: Hydrated Electrons Simulated with Density Functional Theory Are Too Kosmotropic

We have performed the first DFT-based ab initio MD simulations of a hydrated electron (e_aq^-) in the presence of Na⁺, a system chosen because ion-pairing behavior in water depends sensitively on the local hydration structure. Experiments show that e_aq^-'s interact weakly with Na⁺; the e_aq^-'s spectrum blue shifts by only a few tens of meV upon ion pairing without changing shape. We find that the spectrum of the DFT-simulated e_aq^- red shifts and changes shape upon interaction with Na⁺, in contrast with experiment. We show that this is because the hydration structure of the DFT-simulated e_aq^- is too ordered or kosmotropic. Conversely, simulations that produce e_aq^-'s with a less ordered or chaotropic hydration structure form weaker ion pairs with Na⁺, yielding predicted spectral blue shifts in better agreement with experiment. Thus, ab initio simulations based on hybrid GGA DFT functionals fail to produce the correct solvation structure for the hydrated electron.

Competitive Ion Pairing and the Role of Anions in the Behavior of Hydrated Electrons in Electrolytes

Experiments have shown that in the presence of electrolytes, the hydrated electron's absorption spectrum experiences a blue shift whose magnitude depends on both the salt concentration and chemical identity. Previous computer simulations have suggested that the spectral blue shift results from the formation of (cation, electron) contact pairs and that the concentration dependence arises because the number of cations simultaneously paired with the electron increases with increasing concentration. In this work, we perform new simulations to build an atomistic picture that explains the effect of salt identity on the observed hydrated electron spectral shifts. We simulate hydrated electrons in the presence of both monovalent (Na⁺) and divalent (Ca²⁺) cations paired with both Cl^- and a spherical species representing ClO₄^- anions. Our simulations reproduce the experimental observations that divalent ions produce larger blue shifts of the hydrated electron's spectrum than monovalent ions with the same anion and that perchlorate salts show enhanced blue shifts compared to chloride salts with the same cation. We find that these observations can be explained by competitive ion pairing. With small kosmotropic cations such as Na⁺ and Ca²⁺, aqueous chloride salts tend to form (cation, anion) contact pairs, whereas there is little ion pairing between these cations and chaotropic perchlorate anions. Hydrated electrons also strongly interact with these cations, but if the cations are also paired with anions, this affects the free energy of the electron-cation interaction. With chloride salts, hydrated electrons end up in complexes containing multiple cations plus a few anions as well as the electron. Repulsive interactions between the electron and the nearby Cl^- anions reduce the cation-induced spectral blue shift of the hydrated electron. With perchlorate salts, hydrated electrons pair with multiple cations without any associated anions, leading to the largest possible cation-induced spectral blue shift. We also see that the reason multivalent cations produce larger spectral blue shifts than monovalent cations is because hydrated electrons are able to simultaneously pair with a larger number of multivalent cations due to a larger free energy of interaction. Overall, the interaction of hydrated electrons with electrolytes fits well with the Hofmeister series, where the electron behaves as an anion that is slightly more able to break water's H-bond structure than chloride.

Trap-Seeking or Trap-Digging? Photoinjection of Hydrated Electrons into Aqueous NaCl Solutions

It is well-known that when excess electrons are injected into an aqueous solution, they localize and solvate in ∼1 ps. Still debated is whether localization occurs via "trap-digging", in which the electron carves out a suitable localization site, or by "trap-seeking", where the electron prefers to localize at pre-existing low-energy trap sites in solution. To distinguish between these two possible mechanisms, we study the localization dynamics of excess electrons in aqueous NaCl solutions using both ultrafast spectroscopy and mixed quantum-classical molecular dynamics simulations. By introducing pre-existing traps in the form of Na⁺ ions, we can use the cation-induced blue-shift of the hydrated electron's absorption spectrum to directly monitor the site of electron localization. Our experimental and computational results show that the electron prefers to localize directly at the sites of Na⁺ traps; the presence of concentrated electrolytes otherwise has little impact on the way trap-seeking hydrated electrons relax following injection.

Solvent Control of Chemical Identity Can Change Photodissociation into Photoisomerization

In solution-phase chemistry, the solvent is often considered to be merely a medium that allows reacting solutes to encounter each other. In this work, however, we show that moderate locally specific solute-solvent interactions can affect not only the nature of the solute but also the types of reactive chemistry. We use quantum simulation methods to explore how solvent participation in solute chemical identity alters reactions involving the breaking of chemical bonds. In particular, we explore the photoexcitation dynamics of Na₂⁺ dissolved in liquid tetrahydrofuran. In the gas phase, excitation of Na₂⁺ directly leads to dissociation, but in solution, photoexcitation leads to an isomerization reaction involving rearrangement of the first-shell solvent molecules; this isomerization must go to completion before the solute can dissociate. Despite the complexity, the solution-phase reaction dynamics can be captured by a two-dimensional energy surface where one dimension involves only the isomerization of the first-shell solvent molecules.

Bond-Breaking Reactions Encounter Distinct Solvent Environments Causing Breakdown of Linear Response

Solvent effects are important for understanding solution-phase chemical reactions. Surprisingly, very few studies have explored how solvent dynamics change during the course of a reaction with solutes that encounter a wide range of configurations. Here, we use quantum simulation methods to explore the solvent dynamics during a solution-phase bond-breaking reaction: the photodissociation of Na₂⁺ in liquid Ar. We find that the solute experiences a small number of distinct solvent environments that change in a discrete fashion as the bond lengthens. In characterizing the solvent environments, we show also that linear response fails by all measures, even when nonstationarity of solvent dynamics is considered. This observation of distinct solvent response environments with a solvent that can undergo only translational motions highlights the complexity of solute-solvent interactions, but that there are only a few environments gives hope to the idea that solvation dynamics can be understood for solution-phase reactions that explore a wide configuration space.

Understanding the Temperature Dependence and Finite Size Effects in Ab Initio MD Simulations of the Hydrated Electron

The hydrated electron is of interest to both theorists and experimentalists as a paradigm solution-phase quantum system. Although the bulk of the theoretical work studying the hydrated electron is based on mixed quantum/classical (MQC) methods, recent advances in computer power have allowed several attempts to study this object using ab initio methods. The difficulty with employing ab initio methods for this system is that even with relatively inexpensive quantum chemistry methods such as density functional theory (DFT), such calculations are still limited to at most a few tens of water molecules and only a few picoseconds duration, leaving open the question as to whether the calculations are converged with respect to either system size or dynamical fluctuations. Moreover, the ab initio simulations of the hydrated electron that have been published to date have provided only limited analysis. Most works calculate the electron's vertical detachment energy, which can be compared to experiment, and occasionally the electronic absorption spectrum is also computed. Structural features, such as pair distribution functions, are rare in the literature, with the majority of the structural analysis being simple statements that the electron resides in a cavity, which are often based only on a small number of simulation snapshots. Importantly, there has been no ab initio work examining the temperature-dependent behavior of the hydrated electron, which has not been satisfactorily explained by MQC simulations. In this work, we attempt to remedy this situation by running DFT-based ab initio simulations of the hydrated electron as a function of both box size and temperature. We show that the calculated properties of the hydrated electron are not converged even with simulation sizes up to 128 water molecules and durations of several tens of picoseconds. The simulations show significant changes in the water coordination and solvation structure with box size. Our temperature-dependent simulations predict a red-shift of the absorption spectrum (computed using TD-DFT with an optimally tuned range-separated hybrid functional) with increasing temperature, but the magnitude of the predicted red-shift is larger than that observed experimentally, and the absolute position of the calculated spectra are off by over half an eV. The spectral red-shift at high temperatures is accompanied by both a partial loss of structure of the electron's central cavity and an increased radius of gyration that pushes electron density onto and beyond the first solvation shell. Overall, although ab initio simulations can provide some insights into the temperature-dependent behavior of the hydrated electron, the simulation sizes and level of quantum chemistry theory that are currently accessible are inadequate for correctly describing the experimental properties of this fascinating object.

Molecular Dynamics Study of the Thermodynamics of Integer Charge Transfer vs Charge-Transfer Complex Formation in Doped Conjugated Polymers

Molecular dopants such as 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ) can interact with conjugated polymers such as poly(3-hexylthiophene-2,5-diyl) (P3HT) in two different ways: they can undergo integer charge transfer (ICT) or they can form a partial-charge-transfer complex (CTC). Both are seen experimentally, but the CTC has been challenging to characterize, making it difficult to answer questions such as the following. Which polymorph is more stable? Do they have similar barriers for formation? Is there a thermodynamic route to convert one to the other? Here, we study the structure and the thermodynamics of bulk F4TCNQ-doped P3HT with all-atom molecular dynamics simulations, using thermodynamic integration to calculate the relative free energies. We find that the ICT and CTC polymorphs have similar thermodynamic stabilities. The barrier to create the ICT polymorph, however, is lower than that to make the CTC polymorph, because the ICT polymorph has a small critical nucleus, but the critical nucleus for the CTC polymorph is larger than what we can simulate. Moreover, simulated thermal annealing shows that the activation barrier for converting the CTC polymorph to the ICT polymorph is relatively modest. Overall, the simulations explain both the observed structures and the thermodynamics of F4TCNQ-doped P3HT and offer guidelines for targeting the production of a desired polymorph for different applications.

Hydrated Electrons in High-Concentration Electrolytes Interact with Multiple Cations: A Simulation Study

Experimental studies have demonstrated that the hydrated electron's absorption spectrum undergoes a concentration-dependent blue-shift in the presence of electrolytes such as NaCl. The blue-shift increases roughly linearly at low salt concentration but saturates as the solubility limit of the salt is approached. Previous attempts to understand the origin of the concentration-dependent spectral shift using molecular simulation have only examined the interaction between the hydrated electron and a single sodium cation, and these simulations predicted a spectral blue-shift that was an order of magnitude larger than that seen experimentally. Thus, in this paper, we first explore the reasons for the exaggerated spectral blue-shift when a simulated hydrated electron interacts with a single Na⁺. We find that the issue arises from nonpairwise additivity of the Na⁺-e^- and H₂O-e^- pseudopotentials used in the simulation. This effect arises because the solvating water molecules donate charge into the empty orbitals of Na⁺, lowering the effective charge of the cation and thus reducing the excess electron-cation interaction. Careful analysis shows, however, that although this nonpairwise additivity changes the energetics of the electron-Na⁺ interaction, the forces between the electron, Na⁺, and water are unaffected, so that mixed quantum/classical (MQC) simulations produce the correct structure and dynamics. With this in hand, we then use MQC simulations to explore the behavior of the hydrated electron as an explicit function of NaCl salt concentration. We find that the simulations correctly reproduce the observed experimental spectral shifting behavior. The reason for the spectral shift is that as the electrolyte concentration increases, the average number of cations simultaneously interacting in contact pairs with the hydrated electron increases from 1.0 at low concentrations to ∼2.5 near the saturation limit. As the number of cations that interact with the electron increases, the cation/electron interactions becomes slightly weaker, so that the corresponding Na⁺-e^- distance increases with increasing salt concentration. We also find that the dielectric constant of the solution plays little role in the observed spectroscopy, so that the electrolyte-dependent spectral shifts of the hydrated electron are directly related to the concentration-dependent number of closely interacting cations.

How Water-Ion Interactions Control the Formation of Hydrated Electron:Sodium Cation Contact Pairs

Although solvated electrons are a perennial subject of interest, relatively little attention has been paid to the way they behave in aqueous electrolytes. Experimentally, it is known that the hydrated electron's (e_aq^-) absorption spectrum shifts to the blue in the presence of salts, and the magnitude of the shift depends on the ion concentration and the identities of both the cation and anion. Does the blue-shift result from some type of dielectric effect from the bulk electrolyte, or are there specific interactions between the hydrated electron and ions in solution? Previous work has suggested that e_aq^- forms contact pairs with aqueous ions such as Na⁺, leading to the question of what controls the stability of such contact pairs and their possible connection to the observed spectroscopy. In this work, we use mixed quantum/classical simulations to examine the nature of Na⁺:e^- contact pairs in water, using a novel method for quantum umbrella sampling to construct e_aq^--ion potentials of mean force (PMF). We find that the nature of the contact pair PMF depends sensitively on the choice of the classical interactions used to describe the Na⁺-water interactions. When the ion-water interactions are slightly stronger, the corresponding cation:e^- contact pairs form at longer distances and become free energetically less stable. We show that this is because there is a delicate balance between solvation of the cation, solvation of e_aq^- and the direct electronic interaction between the cation and the electron, so that small changes in this balance lead to large changes in the formation and stability of e^--ion contact pairs. In particular, strengthening the ion-water interactions helps to maintain a favorable local solvation environment around Na⁺, which in turn forces water molecules in the first solvation shell of the cation to be unfavorably oriented toward the electron in a contact pair; stronger solvation of the cation also reduces the electronic overlap of e_aq^- with Na⁺. We also find that the calculated spectra of different models of Na⁺:e^- contact pairs do not shift monotonically with cation-electron distance, and that the calculated spectral shifts are about an order of magnitude larger than experiment, suggesting that isolated contact pairs are not the sole explanation for the blue-shift of the hydrated electron's spectrum in the presence of electrolytes.

Ab Initio Simulations of Poorly and Well Equilibrated (CH3CN)n- Cluster Anions: Assigning Experimental Photoelectron Peaks to Surface-Bound Electrons and Solvated Monomer and Dimer Anions

Excess electrons in liquid acetonitrile are of particular interest because they exist in two different forms in equilibrium: they can be present as traditional solvated electrons in a cavity, and they can form some type of solvated molecular anion. Studies of small acetonitrile cluster anions in the gas phase show two isomers with distinct vertical detachment energies, and it is tempting to presume that the two gas-phase cluster anion isomers are precursors of the two excess electron species present in bulk solution. In this paper, we perform DFT-based ab initio molecular dynamics simulations of acetonitrile cluster anions to understand the electronic species that are present and why they have different binding energies. Using a long-range-corrected density functional that was optimally tuned to describe acetonitrile cluster anion structures, we have theoretically explored the chemistry of (CH₃CN)_n^- cluster anions with sizes n = 5, 7, and 10. Because the temperature of the experimental cluster anions is not known, we performed two sets of simulations that investigated how the way in which the cluster anions are prepared affects the excess electron binding motif: one set of simulations simply attached excess electrons to neutral (CH₃CN)_n clusters, providing little opportunity for the clusters to relax in the presence of the excess electron, while the other set allowed the cluster anions to thermally equilibrate near room temperature. We find that both sets of simulations show three distinct electron binding motifs: electrons can attach to the surface of the cluster (dipole-bound) or be present either as solvated monomer anions, CH₃CN^-, or as solvated molecular dimer anions, (CH₃CN)₂^-. All three species have higher binding energies at larger cluster sizes. Thermal equilibration strongly favors the formation of the valence-bound molecular anions relative to surface-bound excess electrons, and the dimer anion becomes more stable than the monomer anion and surface-bound species as the cluster size increases. The calculated photoelectron spectra from our simulations in which there was poor thermal equilibration are in good agreement with experiment, suggesting assignment of the two experimental cluster anion isomers as the surface-bound electron and the solvated molecular dimer anion. The simulations also suggest that the shoulder seen experimentally on the low-energy isomer's detachment peak is not part of a vibronic progression but instead results from molecular monomer anions. Nowhere in the size range that we explore do we see evidence for a nonvalence, cavity-bound interior-solvated electron, indicating that this species is likely only accessible at larger sizes with good thermal equilibration.

Driving Force and Optical Signatures of Bipolaron Formation in Chemically Doped Conjugated Polymers

Molecular dopants are often added to semiconducting polymers to improve electrical conductivity. However, the use of such dopants does not always produce mobile charge carriers. In this work, ultrafast spectroscopy is used to explore the nature of the carriers created following doping of conjugated push-pull polymers with both F₄ TCNQ (2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane) and FeCl₃ . It is shown that for one particular push-pull material, the charge carriers created by doping are entirely non-conductive bipolarons and not single polarons, and that transient absorption spectroscopy following excitation in the infrared can readily distinguish the two types of charge carriers. Based on density functional theory calculations and experiments on multiple push-pull conjugated polymers, it is argued that the size of the donor push units determines the relative stabilities of polarons and bipolarons, with larger donor units stabilizing the bipolarons by providing more area for two charges to co-reside.

Nonequilibrium Solvent Effects during Photodissociation in Liquids: Dynamical Energy Surfaces, Caging, and Chemical Identity

In the gas phase, potential energy surfaces can be used to provide insight into the details of photochemical reaction dynamics. In solution, however, it is unclear what potential energy surfaces, if any, can be used to describe even simple chemical reactions such as the photodissociation of a diatomic solute. In this paper, we use mixed quantum/classical (MQC) molecular dynamics (MD) to study the photodissociation of Na₂⁺ in both liquid Ar and liquid tetrahydrofuran (THF). We examine both the gas-phase potential surfaces and potentials of mean force (PMF), which assume that the solvent remains at equilibrium with the solute throughout the photodissociation process and show that neither resemble a nonequilibrium dynamical energy surface that is generated by taking the time integral of work. For the photodissociation of Na₂⁺ in liquid Ar, the dynamical energy surface shows clear signatures of solvent caging, and the degree of caging is directly related to the mass of the solvent atoms. For Na₂⁺ in liquid THF, local specific interactions between the solute and solvent lead to changes in chemical identity that create a kinetic trap that effectively prevents the molecule from dissociating. The results show that nonequilibrium effects play an important role even in simple solution-phase reactions, requiring the use of dynamical energy surface to understand such chemical events.

Evaluating Simple Ab Initio Models of the Hydrated Electron: The Role of Dynamical Fluctuations

Despite its importance in electron transfer reactions and radiation chemistry, there has been disagreement over the fundamental nature of the hydrated electron, such as whether or not it resides in a cavity. Mixed quantum/classical simulations of the hydrated electron give different structures depending on the pseudopotential employed, and ab initio models of computational necessity use small numbers of water molecules and/or provide insufficient statistics to compare to experimental observables. A few years ago, Kumar et al. (J. Phys. Chem. A 2015, 119, 9148) proposed a minimalist ab initio model of the hydrated electron with only a small number of explicitly treated water molecules plus a polarizable continuum model (PCM). They found that the optimized geometry had four waters arranged tetrahedrally around a central cavity, and that the calculated vertical detachment energy and radius of gyration agreed well with experiment, results that were largely independent of the level of theory employed. The model, however, is based on a fixed structure at 0 K and does not explicitly incorporate entropic contributions or the thermal fluctuations that should be associated with the room-temperature hydrated electron. Thus, in this paper, we extend the model of Kumar et al. by running Born-Oppenheimer molecular dynamics (BOMD) of a small number of water molecules with an excess electron plus PCM at room temperature. We find that when thermal fluctuations are introduced, the level of theory chosen becomes critical enough when only four waters are used that one of the waters dissociates from the cluster with certain density functionals. Moreover, even with an optimally tuned range-separated hybrid functional, at room temperature the tetrahedral orientation of the 0 K first-shell waters is entirely lost and the central cavity collapses, a process driven by the fact that the explicit water molecules prefer to make H-bonds with each other more than with the excess electron. The resulting average structure is quite similar to that produced by a noncavity mixed quantum/classical model, so that the minimalist 4-water BOMD models suffer from problems similar to those of noncavity models, such as predicting the wrong sign of the hydrated electron's molar solvation volume. We also performed BOMD with 16 explicit water molecules plus an extra electron and PCM. We find that the inclusion of an entire second solvation shell of explicit water leads to little change in the outcome from when only four waters were used. In fact, the 16-water simulations behave much like those of water cluster anions, in which the electron localizes at the cluster surface, showing that PCM is not acceptable for use in minimalist models to describe the behavior of the bulk hydrated electron. For both the 4- and 16-water models, we investigate how the introduction of thermal motions alters the predicted absorption spectrum, vertical detachment energy, and resonance Raman spectrum of the simulated hydrated electron. We also present a set of structural criteria that can be used to numerically determine how cavity-like (or not) a particular hydrated electron model is. All of the results emphasize that the hydrated electron is a statistical object whose properties are inadequately captured using only a small number of explicit waters, and that a proper treatment of thermal fluctuations is critical to understanding the hydrated electron's chemical and physical behavior.

The Role of the Solvent in the Condensed-Phase Dynamics and Identity of Chemical Bonds: The Case of the Sodium Dimer Cation in THF

When a solute molecule is placed in solution, is it acceptable to presume that its electronic structure is essentially the same as that in the gas phase? In this paper, we address this question from a simulation perspective for the case of the sodium dimer cation (Na₂⁺) molecule in both liquid Ar and liquid tetrahydrofuran (THF). In previous work, we showed that, when local specific interactions between a solute and solvent are energetically on the order of a hydrogen bond, the solvent can become part of the chemical identity of the solute. Here, using mixed quantum/classical molecular dynamics simulations, we see that, for the Na₂⁺ molecule, solute-solvent interactions lead to two stable, chemically distinct coordination states (Na(THF)₄-Na(THF)₅⁺ and Na(THF)₅-Na(THF)₅⁺) that are not only stable themselves as gas-phase molecules but that also have a completely new electronic structure with important implications for the excited-state photodissociation of this molecule in the condensed phase. Furthermore, we show through a set of comparative classical simulations that treating the solute's bonding electron explicitly quantum mechanically is necessary to understand both the ground-state dynamics and chemical identity of this simple diatomic molecule; even use of the quantum-derived potential of mean force is insufficient to describe the behavior of the molecule classically. Finally, we calculate the results of a proposed transient hole-burning experiment that could be used to spectroscopically disentangle the presence of the different coordination states.

Tunable Dopants with Intrinsic Counterion Separation Reveal the Effects of Electron Affinity on Dopant Intercalation and Free Carrier Production in Sequentially Doped Conjugated Polymer Films

Carrier mobility in doped conjugated polymers is limited by Coulomb interactions with dopant counterions. This complicates studying the effect of the dopant's oxidation potential on carrier generation because different dopants have different Coulomb interactions with polarons on the polymer backbone. Here, dodecaborane (DDB)-based dopants are used, which electrostatically shield counterions from carriers and have tunable redox potentials at constant size and shape. DDB dopants produce mobile carriers due to spatial separation of the counterion, and those with greater energetic offsets produce more carriers. Neutron reflectometry indicates that dopant infiltration into conjugated polymer films is redox-potential-driven. Remarkably, X-ray scattering shows that despite their large 2-nm size, DDBs intercalate into the crystalline polymer lamellae like small molecules, indicating that this is the preferred location for dopants of any size. These findings elucidate why doping conjugated polymers usually produces integer, rather than partial charge transfer: dopant counterions effectively intercalate into the lamellae, far from the polarons on the polymer backbone. Finally, it is shown that the IR spectrum provides a simple way to determine polaron mobility. Overall, higher oxidation potentials lead to higher doping efficiencies, with values reaching 100% for driving forces sufficient to dope poorly crystalline regions of the film.

The Fluxional Nature of the Hydrated Electron: Energy and Entropy Contributions to Aqueous Electron Free Energies

There has been a great deal of recent controversy over the structure of the hydrated electron and whether it occupies a cavity or contains a significant number of interior waters (noncavity). The questions we address in this work are, from a free energy perspective, how different are these proposed structures? Do the different structures all lie along a single continuum, or are there significant differences (i.e., free energy barriers) between them? To address these questions, we have performed a series of one-electron calculations using umbrella sampling with quantum biased molecular dynamics along a coordinate that directly reflects the number of water molecules in the hydrated electron's interior. We verify that a standard cavity model of the hydrated electron behaves essentially as a hard sphere: the model is dominated by repulsion at short range such that water is expelled from a local volume around the electron, leading to a water solvation shell like that of a pseudohalide ion. The repulsion is much larger than thermal energies near room temperature, explaining why such models exhibit properties with little temperature dependence. On the other hand, our calculations reveal that a noncavity model is highly fluxional, meaning that thermal motions cause the number of interior waters to fluctuate from effectively zero (i.e., a cavity-type electron) to potentially above the bulk water density. The energetic contributions in the noncavity model are still repulsive in the sense that they favor cavity formation, so the fluctuations in structure are driven largely by entropy: the entropic cost for expelling water from a region of space is large enough that some water is still driven into the electron's interior. As the temperature is lowered and entropy becomes less important, the noncavity electron's structure is predicted to become more cavity-like, consistent with the observed temperature dependence of the hydrated electron's properties. Thus, we argue that although the specific noncavity model we study overestimates the preponderance of fluctuations involving interior water molecules, with appropriate refinements to correctly capture the true average number of interior waters and molar solvation volume, a fluxional model likely makes the most sense for understanding the various experimental properties of the hydrated electron.

Dopant-Induced Ordering of Amorphous Regions in Regiorandom P3HT

Despite the fact that molecular doping of semiconducting polymers has emerged as a valuable strategy for improving the performance of organic electronic devices, the fundamental dopant-polymer interactions are not fully understood. Here we use 2-D grazing incidence wide-angle X-ray scattering (GIWAXS) to demonstrate that adding oxidizing small-molecule dopants, such as 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F₄TCNQ) and FeCl₃, into the amorphous conjugated polymer, regiorandom poly(3-hexylthiophene-2,5-diyl) (RRa-P3HT), improves polymer ordering and induces a change in domain orientation from isotropic to mostly edge-on. Doping thus causes RRa-P3HT to behave similarly to the more ordered regioregular P3HT. By comparing the optical, electrical, and structural properties of RRa-P3HT films doped with F₄TNCQ and FeCl₃ and those infiltrated with 7,7,8,8-tetracyanoquinodimethane (TCNQ), which occupies a similar volume as F₄TCNQ but does not dope RRa-P3HT, we show that the increased ordering results not from the ability of the dopant to fill space but instead from the need to delocalize charge on the polymer in more than one dimension.

Dodecaborane-Based Dopants Designed to Shield Anion Electrostatics Lead to Increased Carrier Mobility in a Doped Conjugated Polymer

One of the most effective ways to tune the electronic properties of conjugated polymers is to dope them with small-molecule oxidizing agents, creating holes on the polymer and molecular anions. Undesirably, strong electrostatic attraction from the anions of most dopants localizes the holes created on the polymer, reducing their mobility. Here, a new strategy utilizing a substituted boron cluster as a molecular dopant for conjugated polymers is employed. By designing the cluster to have a high redox potential and steric protection of the core-localized electron density, highly delocalized polarons with mobilities equivalent to films doped with no anions present are obtained. AC Hall effect measurements show that P3HT films doped with these boron clusters have conductivities and polaron mobilities roughly an order of magnitude higher than films doped with F₄ TCNQ, even though the boron-cluster-doped films have poor crystallinity. Moreover, the number of free carriers approximately matches the number of boron clusters, yielding a doping efficiency of ≈100%. These results suggest that shielding the polaron from the anion is a critically important aspect for producing high carrier mobility, and that the high polymer crystallinity required with dopants such as F₄ TCNQ is primarily to keep the counterions far from the polymer backbone.

Ultrafast transient absorption spectroscopy of doped P3HT films: distinguishing free and trapped polarons

It is generally presumed that the vast majority of carriers created by chemical doping of semiconducting polymer films are coulombically trapped by the counteranion, with only a small fraction that are free and responsible for the increased conductivity essential for organic electronic applications. At higher doping levels, it is also possible for bipolarons to form, which are expected to be less conductive than single polarons. Unfortunately, there is no simple way to distinguish free polarons, trapped polarons and bipolarons using steady-state spectroscopy. Thus, in this work, we use ultrafast transient absorption spectroscopy to study the dynamics of polarons in 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TNCQ)-doped films of poly(3-hexylthiophene-2,5-diyl) (P3HT) as a function of dopant concentration and excitation wavelength. When exciting on the red side of the polaron P1 transition, our transient absorption spectra and kinetics match well with what is expected for free 2-D-delocalized polarons; the measurements are not consistent with a recent theory of doped conjugated polymer electronic structure that suggests that the half-filled state lies deeper in the conduction band rather than in the bandgap. As we tune the excitation wavelength to the blue, our measurements reveal an increasing amount of slower transient kinetics that are consistent with the presence of coulombically-trapped polarons rather than bipolarons. Taking advantage of their distinct ultrafast relaxation kinetics as a type of action spectroscopy, we are able to extract the steady-state absorption spectra of free and trapped polarons as a function of dopant concentration. By comparing the results to theoretical models, we determine that in F4TCNQ-doped P3HT films, trapped polarons sit ∼0.4 nm away from the anion while free polarons reside between 0.7 and 0.9 nm from the counteranion. Perhaps counterintuitively, the ratio of trapped to free polarons increases at higher doping levels, an observation that is consistent with a plateau in the concentration-dependent conductivity of F4TCNQ-doped P3HT films.

Thermal Equilibration Controls H-Bonding and the Vertical Detachment Energy of Water Cluster Anions

One of the outstanding puzzles in the photoelectron spectroscopy of water anion clusters, which serve as precursors to the hydrated electron, is that the excess electron has multiple vertical detachment energies (VDEs), with different groups seeing different distributions of VDEs. We have studied the photoelectron spectroscopy of water cluster anions using simulation techniques designed to mimic the different ways that water cluster anions are produced experimentally. Our simulations take advantage of density functional theory-based Born-Oppenheimer molecular dynamics with an optimally tuned range-separated hybrid functional that is shown to give outstanding accuracy for calculating electron binding energies for this system. We find that our simulations are able to accurately reproduce the experimentally observed VDEs for cluster anions of different sizes, with different VDE distributions observed depending on how the water cluster anions are prepared. For cluster anion sizes up to 20 water molecules, we see that the excess electron always resides on the surface of the cluster and that the different discrete VDEs result from the discrete number of hydrogen bonds made to the electron by water molecules on the surface. Clusters that are less thermally equilibrated have surface waters that tend to make single H-bonds to the electron, resulting in lower VDEs, while clusters that are more thermally equilibrated have surface waters that prefer to make two H-bonds to the electron, resulting in higher VDEs.

Short-Range Electron Correlation Stabilizes Noncavity Solvation of the Hydrated Electron

The hydrated electron, e^-_(aq), has often served as a model system to understand the influence of condensed-phase environments on electronic structure and dynamics. Despite over 50 years of study, however, the basic structure of e^-_(aq) is still the subject of controversy. In particular, the structure of e^-_(aq) was long assumed to be an electron localized within a solvent cavity, in a manner similar to halide solvation. Recently, however, we suggested that e^-_(aq) occupies a region of enhanced water density with little or no discernible cavity. The potential we developed was only subtly different from those that give rise to a cavity solvation motif, which suggests that the driving forces for noncavity solvation involve subtle electron-water attractive interactions at close distances. This leads to the question of how dispersion interactions are treated in simulations of the hydrated electron. Most dispersion potentials are ad hoc or are not designed to account for the type of close-contact electron-water overlap that might occur in the condensed phase, and where short-range dynamic electron correlation is important. To address this, in this paper we develop a procedure to calculate the potential energy surface between a single water molecule and an excess electron with high-level CCSD(T) electronic structure theory. By decomposing the electron-water potential into its constituent energetic contributions, we find that short-range electron correlation provides an attraction of comparable magnitude to the mean-field interactions between the electron and water. Furthermore, we find that by reoptimizing a popular cavity-forming one-electron model potential to better capture these attractive short-range interactions, the enhanced description of correlation predicts a noncavity e^-_(aq) with calculated properties in better agreement with experiment. Although much attention has been placed on the importance of long-range dispersion interactions in water cluster anions, our study reveals that largely unexplored short-range correlation effects are crucial in dictating the solvation structure of the condensed-phase hydrated electron.

Free Energies of Cavity and Noncavity Hydrated Electrons Near the Instantaneous Air/Water Interface

The properties of the hydrated electron at the air/water interface are computed for both a cavity and a noncavity model using mixed quantum/classical molecular dynamics simulation. We take advantage of our recently developed formalism for umbrella sampling with a restrained quantum expectation value to calculate free-energy profiles of the hydrated electron's position relative to the water surface. We show that it is critical to use an instantaneous description of the air/water interface rather than the Gibbs' dividing surface to obtain accurate potentials of mean force. We find that noncavity electrons, which prefer to encompass several water molecules, avoid the interface where water molecules are scarce. In contrast, cavity models of the hydrated electron, which prefer to expel water, have a local free-energy minimum near the interface. When the cavity electron occupies this minimum, its absorption spectrum is quite red-shifted, its binding energy is significantly lowered, and its dynamics speed up quite a bit compared with the bulk, features that have not been found by experiment. The surface activity of the electron therefore serves as a useful test of cavity versus noncavity electron solvation.

Discovery of a Novel 2,6-Disubstituted Glucosamine Series of Potent and Selective Hexokinase 2 Inhibitors

A novel series of potent and selective hexokinase 2 (HK2) inhibitors, 2,6-disubstituted glucosamines, has been identified based on HTS hits, exemplified by compound 1. Inhibitor-bound crystal structures revealed that the HK2 enzyme could adopt an "induced-fit" conformation. The SAR study led to the identification of potent HK2 inhibitors, such as compound 34 with greater than 100-fold selectivity over HK1. Compound 25 inhibits in situ glycolysis in a UM-UC-3 bladder tumor cell line via (13)CNMR measurement of [3-(13)C]lactate produced from [1,6-(13)C2]glucose added to the cell culture.

Overcoming Film Quality Issues for Conjugated Polymers Doped with F4TCNQ by Solution Sequential Processing: Hall Effect, Structural, and Optical Measurements

We demonstrate that solution-sequential processing (SqP) can yield heavily doped pristine-quality films when used to infiltrate the molecular dopant 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ) into pure poly(3-hexylthiophene) (P3HT) polymer layers. Profilometry measurements show that the SqP method produces doped films with essentially the same surface roughness as pristine films, and 2-D grazing-incidence wide-angle X-ray scattering (GIWAXS) confirms that SqP preserves both the size and orientation of the pristine polymer's crystallites. Unlike traditional blend-cast F4TCNQ/P3HT doped films, our sequentially processed layers have tunable and reproducible conductivities reaching as high as 5.5 S/cm even when measured over macroscopic (>1 cm) distances. The high conductivity and superb film quality allow for meaningful Hall effect measurements, which reveal p-type conduction and carrier concentrations tunable from 10(16) to 10(20) cm(-3) and hole mobilities ranging from ∼0.003 to 0.02 cm(2) V(-1) s(-1) at room temperature over the doping levels examined.

Extensive Penetration of Evaporated Electrode Metals into Fullerene Films: Intercalated Metal Nanostructures and Influence on Device Architecture

Although it is known that evaporated metals can penetrate into films of various organic molecules that are a few nanometers thick, there has been little work aimed at exploring the interaction of the common electrode metals used in devices with fullerene derivatives, such as organic photovoltaics (OPVs) or perovskite solar cells that use fullerenes as electron transport layers. In this paper, we show that when commonly used electrode metals (e.g., Au, Ag, Al, Ca, etc.) are evaporated onto films of fullerene derivatives (such as [6,6]-phenyl-C61-butyric acid methyl ester (PCBM)), the metal penetrates many tens of nanometers into the fullerene layer. This penetration decreases the effective electrical thickness of fullerene-based sandwich structure devices, as measured by the device's geometric capacitance, and thus significantly alters the device physics. For the case of Au/PCBM, the metal penetrates a remarkable 70 nm into the fullerene, and we see penetration of similar magnitude in a wide variety of fullerene derivative/evaporated metal combinations. Moreover, using transmission electron microscopy to observed cross-sections of the films, we show that when gold is evaporated onto poly(3-hexylthiophene) (P3HT)/PCBM sequentially processed OPV quasi-bilayers, Au nanoparticles with diameters of ∼3-20 nm are formed and are dispersed entirely throughout the fullerene-rich overlayer. The plasmonic absorption and scattering from these nanoparticles are readily evident in the optical transmission spectrum, demonstrating that the interpenetrated metal significantly alters the optical properties of fullerene-rich active layers. This opens a number of possibilities in terms of contact engineering and light management so that metal penetration in devices that use fullerene derivatives could be used to advantage, making it critical that researchers are aware of the electronic and optical consequences of exposing fullerene-derivative films to evaporated electrode metals.

POLARON DYNAMICS. Long-lived photoinduced polaron formation in conjugated polyelectrolyte-fullerene assemblies

The efficiency of biological photosynthesis results from the exquisite organization of photoactive elements that promote rapid movement of charge carriers out of a critical recombination range. If synthetic organic photovoltaic materials could mimic this assembly, charge separation and collection could be markedly enhanced. We show that micelle-forming cationic semiconducting polymers can coassemble in water with cationic fullerene derivatives to create photoinduced electron-transfer cascades that lead to exceptionally long-lived polarons. The stability of the polarons depends on the organization of the polymer-fullerene assembly. Properly designed assemblies can produce separated polaronic charges that are stable for days or weeks in aqueous solution.

Photo-redox activated drug delivery systems operating under two photon excitation in the near-IR

We report the design and synthesis of a nano-container consisting of mesoporous silica nanoparticles with the pore openings covered by "snap-top" caps that are opened by near-IR light. A photo transducer molecule that is a reducing agent in an excited electronic state is covalently attached to the system. Near IR two-photon excitation causes inter-molecular electron transfer that reduces a disulfide bond holding the cap in place, thus allowing the cargo molecules to escape. We describe the operation of the "snap-top" release mechanism by both one- and two-photon activation. This system presents a proof of concept of a near-IR photoredox-induced nanoparticle delivery system that may lead to a new type of photodynamic drug release therapy.

Long residence time inhibition of EZH2 in activated polycomb repressive complex 2

EZH2/PRC2 catalyzes transcriptionally repressive methylation at lysine 27 of histone H3 and has been associated with numerous cancer types. Point mutations in EZH2 at Tyr641 and Ala677 identified in non-Hodgkin lymphomas alter substrate specificity and result in increased trimethylation at histone H3K27. Interestingly, EZH2/PRC2 is activated by binding H3K27me3 marks on histones, and this activation is proposed as a mechanism for self-propagation of gene silencing. Recent work has identified GSK126 as a potent, selective, SAM-competitive inhibitor of EZH2 capable of globally decreasing H3K27 trimethylation in cells. Here we show that activation of PRC2 by an H3 peptide trimethylated at K27 is primarily an effect on the rate-limiting step (kcat) with no effect on substrate binding (Km). Additionally, GSK126 is shown to have a significantly longer residence time of inhibition on the activated form of EZH2/PRC2 as compared to unactivated EZH2/PRC2. Overall inhibition constant (Ki*) values for GSK126 were determined to be as low as 93 pM and appear to be driven by slow dissociation of inhibitor from the activated enzyme. The data suggest that activation of EZH2 allows the enzyme to adopt a conformation that possesses greater affinity for GSK126. The long residence time of GSK126 may be beneficial in vivo and may result in durable target inhibition after drug systemic clearance.

To be or not to be in a cavity: the hydrated electron dilemma

The hydrated electron-the species that results from the addition of a single excess electron to liquid water-has been the focus of much interest both because of its role in radiation chemistry and other chemical reactions, and because it provides for a deceptively simple system that can serve as a means to confront the predictions of quantum molecular dynamics simulations with experiment. Despite all this interest, there is still considerable debate over the molecular structure of the hydrated electron: does it occupy a cavity, have a significant number of interior water molecules, or have a structure somewhere in between? The reason for all this debate is that different computer simulations have produced each of these different structures, yet the predicted properties for these different structures are still in reasonable agreement with experiment. In this Feature Article, we explore the reasons underlying why different structures are produced when different pseudopotentials are used in quantum simulations of the hydrated electron. We also show that essentially all the different models for the hydrated electron, including those from fully ab initio calculations, have relatively little direct overlap of the electron's wave function with the nearby water molecules. Thus, a non-cavity hydrated electron is better thought of as an "inverse plum pudding" model, with interior waters that locally expel the surrounding electron's charge density. Finally, we also explore the agreement between different hydrated electron models and certain key experiments, such as resonance Raman spectroscopy and the temperature dependence and degree of homogeneous broadening of the optical absorption spectrum, in order to distinguish between the different simulated structures. Taken together, we conclude that the hydrated electron likely has a significant number of interior water molecules.

Response to "Comment on 'Going beyond the frozen core approximation: development of coordinate-dependent pseudopotentials and application to Na2(+)'" [J. Chem. Phys. 139, 147101 (2013)]

Stoll, Fuentealba, and Szentpály (SFS) argue that the coordinate-dependent pseudopotential we developed for the sodium dimer cation molecule is inferior to other potentials that have been presented in the literature for this molecule. The goal of our work, however, was to present a novel method for the development of rigorous coordinate-dependent pseudopotentials. Our method is designed to reproduce all-electron Hartree-Fock calculations without the inclusion of adjustable parameters. Moreover, our method starts from the superposition of unoptimized, non-norm-conserved atomic potentials, so that when complete, the resulting norm-conserving potential can reproduce an all-electron Hartree-Fock calculation without the inclusion of adjustable parameters. We chose the sodium dimer cation system as a proof of principle for our method, and showed that our method does indeed allow a one-electron calculation to correctly reproduce the all-electron Hartree-Fock calculation from bonding to the dissociation limit. Our purpose in developing this method is to use such potentials in condensed-phase mixed quantum/classical molecular dynamics simulations, where inclusion of valence polarization effects is unimportant or can be added on after the fact. Thus we do not claim that our method provides a potential that is superior to potentials that have been specifically constructed to go beyond the static exchange approximation and/or include valence polarization effects-such potentials are beyond the scope of our work. We also note that although we made a numerical error in the application of our method to Na2(+) in our original work [A. Kahros and B. J. Schwartz, J. Chem. Phys. 138, 054110 (2013)] that led to an overestimation of the magnitude of core polarization effects for this particular molecule, out method does work as derived for this molecule and the error does not affect the significance of our method or its general applicability.

Activation of snap-top capped mesoporous silica nanocontainers using two near-infrared photons

Photoactivation of "snap-top" stoppers over the pore openings of mesoporous silica nanoparticles releases intact cargo molecules from the pores. The on-command release can be stimulated by either one UV photon or two coherent near-IR photons. Two-photon activation is particularly desirable for use in biological systems because it enables good tissue penetration and precise spatial control. Stoppers were assembled by first binding photolabile coumarin-based molecules to the nanoparticle surface. Then, after the particles were loaded with cargo, bulky β-cyclodextrin (CD) molecules were noncovalently associated with the substituted coumarin molecule, blocking the pores and preventing the cargo from escaping. One-photon excitation at 376 nm or two-photon excitation at 800 nm cleaves the bond holding the coumarin to the nanopore, releasing both the CD cap and the cargo. The dynamics of both the cleavage of the cap and the cargo release was monitored using fluorescence spectroscopy. This system traps intact cargo molecules without the necessity of chemical modification, releases them with tissue-penetrating near-IR light, and has possible applications in photostimulated drug delivery.

Self-assembling semiconducting polymers--rods and gels from electronic materials

In an effort to favor the formation of straight polymer chains without crystalline grain boundaries, we have synthesized an amphiphilic conjugated polyelectrolyte, poly(fluorene-alt-thiophene) (PFT), which self-assembles in aqueous solutions to form cylindrical micelles. In contrast to many diblock copolymer assemblies, the semiconducting backbone runs parallel, not perpendicular, to the long axis of the cylindrical micelle. Solution-phase micelle formation is observed by X-ray and visible light scattering. The micelles can be cast as thin films, and the cylindrical morphology is preserved in the solid state. The effects of self-assembly are also observed through spectral shifts in optical absorption and photoluminescence. Solutions of higher-molecular-weight PFT micelles form gel networks at sufficiently high aqueous concentrations. Rheological characterization of the PFT gels reveals solid-like behavior and strain hardening below the yield point, properties similar to those found in entangled gels formed from surfactant-based micelles. Finally, electrical measurements on diode test structures indicate that, despite a complete lack of crystallinity in these self-assembled polymers, they effectively conduct electricity.

Ultrafast Studies of Exciton Migration and Polaron Formation in Sequentially Solution-Processed Conjugated Polymer/Fullerene Quasi-Bilayer Photovoltaics

We examine the ultrafast dynamics of exciton migration and polaron production in sequentially processed 'quasi-bilayer' and preblended 'bulk heterojunction' (BHJ) solar cells based on conjugated polymer films that contain the same total amount of fullerene. We find that even though the polaron yields are similar, the dynamics of polaron production are significantly slower in quasi-bilayers than BHJs. We argue that the different polaron production dynamics result from the fact that (1) there is significantly less fullerene inside the polymer in quasi-bilayers than in BHJs and (2) sequential processing yields polymer layers that are significantly more ordered than BHJs. We also argue that thermal annealing improves the performance of quasi-bilayer solar cells not because annealing drives additional fullerene into the polymer but because annealing improves the fullerene crystallinity. All of the results suggest that sequential processing remains a viable alternative for producing polymer/fullerene solar cells with a nanometer-scale architecture that differs from BHJs.