Mechanism for plasmon-generated solvated electrons.

Solvated electrons are powerful reducing agents capable of driving some of the most energetically expensive reduction reactions. Their generation under mild and sustainable conditions remains challenging though. Using near-ultraviolet irradiation under low-intensity one-photon conditions coupled with electrochemical and optical detection, we show that the yield of solvated electrons in water is increased more than 10 times for nanoparticle-decorated electrodes compared to smooth silver electrodes. Based on the simulations of electric fields and hot carrier distributions, we determine that hot electrons generated by plasmons are injected into water to form solvated electrons. Both yield enhancement and hot carrier production spectrally follow the plasmonic near-field. The ability to enhance solvated electron yields in a controlled manner by tailoring nanoparticle plasmons opens up a promising strategy for exploiting solvated electrons in chemical reactions.

A Molecular Expression for "Line Tension".

"Line tension", a concept that features in an additional term to the Young's equation, was introduced to describe the size dependence of contact angles of nanodroplets on surfaces. Although this concept describes the observations in a succinct, elegant manner, theorists have long had misgivings about the physical interpretation of the phenomenon. Papers have been published that attempt to nail down its value, which is reportedly very small (∼10 pN) and evidently even the sign has been uncertain. Attempts to interpret it in a mechanical manner analogous to interfacial tension, i.e., due to the curvature of the three-phase contact line, have run into conceptual problems that require invocations of ever more complex models. In this work, we have used molecular simulations to systematically relate "line tension" to the additional free energy per unit length of the three-phase line and found no direct relation. However, when we rederived the Young's equation without ignoring the interfacial molecules, we found a physically satisfying explanation for the size dependence of the contact angle of nanodroplets without invoking the curvature of the three-phase contact line. The new model does not have the elegant form of the modified Young's equation, but each parameter in it has an unambiguous physical interpretation. An approximate form of this model, linearized in the inverse droplet radius, yields a quantity that is mathematically analogous to what is conventionally called "line tension", but unpacked at the molecular level, showing that it is unrelated to a restoring force associated with the curvature of the macroscopic three-phase contact line.

Imaging Heterogeneous 3D Dynamics of Individual Solutes in a Polyelectrolyte Brush.

Understanding molecular transport in polyelectrolyte brushes (PEBs) is crucial for applications such as separations, drug delivery, anti-fouling, and biosensors, where structural features of the polymer control intermolecular interactions. The complex structure and local heterogeneity of PEBs, while theoretically predicted, are not easily accessed with conventional experimental methods. In this work, we use 3D single-molecule tracking to understand transport behavior within a cationic poly(2-(N,N-dimethylamino)ethyl acrylate) (PDMAEA) brush using an anionic dye, Alexa Fluor 546, as the probe. The analysis is done by a parallelized, unbiased 3D tracking algorithm. Our results explicitly demonstrate that spatial heterogeneity within the brush manifests as heterogeneity of single-molecule displacements. Two distinct populations of probe motion are identified, with anticorrelated axial and lateral transport confinement, which we believe to correspond to intra- vs inter-chain probe motion.

Accumulation and ordering of P3HT oligomers at the liquid-vapor interface with implications for thin-film morphology.

The morphology of semiconducting polymer thin films is known to have a profound effect on their opto-electronic properties. Although considerable efforts have been made to control and understand the processes which influence the structures of these systems, it remains largely unclear what physical factors determine the arrangement of polymer chains in spin-cast films. Here, we investigate the role that the liquid-vapor interfaces in chlorobenzene solutions of poly(3-hexylthiophene) [P3HT] play in the conformational geometries adopted by the polymers. Using all-atom molecular dynamics (MD), and supported by toy-model simulations, we demonstrate that, with increasing concentration, P3HT oligomers in solution exhibit a strong propensity for the liquid-vapor interface. Due to the differential solubility of the backbone and side chains of the oligomers, in the vicinity of this interface, hexyl chains and the thiophene rings, have a clear orientational preference with respect to the liquid surface. At high concentrations, we additionally establish a substantial degree of inter-oligomer alignment and thiophene ring stacking near the interface. Our results broadly concur with the limited existing experimental evidence and we suggest that the interfacial structure can act as a template for film structure. We argue that the differences in solvent affinity of the side chain and backbone moieties are the driving force for the anisotropic orientations of the polymers near the interface. This finer grained description contrasts with the usual monolithic characterization of polymer units. Since this phenomenon can be controlled by concurrent chemical design and the choice of solvents, this work establishes a fabrication principle which may be useful to develop more highly functional polymer films.

Predicting optical spectra for optoelectronic polymers using coarse-grained models and recurrent neural networks.

Coarse-grained modeling of conjugated polymers has become an increasingly popular route to investigate the physics of organic optoelectronic materials. While ultraviolet (UV)-vis spectroscopy remains one of the key experimental methods for the interrogation of these materials, a rigorous bridge between simulated coarse-grained structures and spectroscopy has not been established. Here, we address this challenge by developing a method that can predict spectra of conjugated polymers directly from coarse-grained representations while avoiding repetitive procedures such as ad hoc back-mapping from coarse-grained to atomistic representations followed by spectral computation using quantum chemistry. Our approach is based on a generative deep-learning model: the long-short-term memory recurrent neural network (LSTM-RNN). The latter is suggested by the apparent similarity between natural languages and the mathematical structure of perturbative expansions of, in our case, excited-state energies perturbed by conformational fluctuations. We also use this model to explore the level of sensitivity of spectra to the coarse-grained representation back-mapping protocol. Our approach presents a tool uniquely suited for improving postsimulation analysis protocols, as well as, potentially, for including spectral data as input in the refinement of coarse-grained potentials.

Aggregation of Charge Acceptors on Nanocrystal Surfaces Alters Rates of Photoinduced Electron Transfer.

Semiconductor nanocrystals (NCs) interfaced with molecular ligands that function as charge and energy acceptors are an emerging platform for the design of light-harvesting, photon-upconverting, and photocatalytic materials. However, NC systems explored for these applications often feature high concentrations of bound acceptor ligands, which can lead to ligand-ligand interactions that may alter each system's ability to undergo charge and energy transfer. Here, we demonstrate that aggregation of acceptor ligands impacts the rate of photoinduced NC-to-ligand charge transfer between lead(II) sulfide (PbS) NCs and perylenediimide (PDI) electron acceptors. As the concentration of PDI acceptors is increased, we find the average electron transfer rate from PbS to PDI ligands decreases by nearly an order of magnitude. The electron transfer rate slowdown with increasing PDI concentration correlates strongly with the appearance of PDI aggregates in steady-state absorption spectra. Electronic structure calculations and molecular dynamics (MD) simulations suggest PDI aggregation slows the rate of electron transfer by reducing orbital overlap between PbS charge donors and PDI charge acceptors. While we find aggregation slows electron transfer in this system, the computational models we employ predict ligand aggregation could also be used to speed electron transfer by producing delocalized states that exhibit improved NC-molecule electronic coupling and energy alignment with NC conduction band states. Our results demonstrate that ligand aggregation can alter rates of photoinduced electron transfer between NCs and organic acceptor ligands and should be considered when designing hybrid NC:molecule systems for charge separation.

Single-Molecule Dynamics Reflect IgG Conformational Changes Associated with Ion-Exchange Chromatography.

Conformational changes of antibodies and other biologics can decrease the effectiveness of pharmaceutical separations. Hence, a detailed mechanistic picture of antibody-stationary phase interactions that occur during ion-exchange chromatography (IEX) can provide critical insights. This work examines antibody conformational changes and how they perturb antibody motion and affect ensemble elution profiles. We combine IEX, three-dimensional single-protein tracking, and circular dichroism spectroscopy to investigate conformational changes of a model antibody, immunoglobulin G (IgG), as it interacts with the stationary phase as a function of salt conditions. The results indicate that the absence of salt enhances electrostatic attraction between IgG and the stationary phase, promotes surface-induced unfolding, slows IgG motion, and decreases elution from the column. Our results reveal previously unreported details of antibody structural changes and their influence on macroscale elution profiles.

Direct observation of backbone planarization via side-chain alignment in single bulky-substituted polythiophenes.

The backbone conformation of conjugated polymers affects, to a large extent, their optical and electronic properties. The usually flexible substituents provide solubility and influence the packing behavior of conjugated polymers in films or in bad solvents. However, the role of the side chains in determining and potentially controlling the backbone conformation, and thus the optical and electronic properties on the single polymer level, is currently under debate. Here, we investigate directly the impact of the side chains by studying the bulky-substituted poly(3-(2,5-dioctylphenyl)thiophene) (PDOPT) and the common poly(3-hexylthiophene) (P3HT), both with a defined molecular weight and high regioregularity, using low-temperature single-chain photoluminescence (PL) spectroscopy and quantum-classical simulations. Surprisingly, the optical transition energy of PDOPT is significantly (∼2,000 cm-1 or 0.25 eV) red-shifted relative to P3HT despite a higher static and dynamic disorder in the former. We ascribe this red shift to a side-chain induced backbone planarization in PDOPT, supported by temperature-dependent ensemble PL spectroscopy. Our atomistic simulations reveal that the bulkier 2,5-dioctylphenyl side chains of PDOPT adopt a clear secondary helical structural motif and thus protect conjugation, i.e., enforce backbone planarity, whereas, for P3HT, this is not the case. These different degrees of planarity in both thiophenes do not result in different conjugation lengths, which we found to be similar. It is rather the stronger electronic coupling between the repeating units in the more planar PDOPT which gives rise to the observed spectral red shift as well as to a reduced calculated electron-hole polarization.

Impact of backbone fluorination on nanoscale morphology and excitonic coupling in polythiophenes.

Fluorination represents an important strategy in developing high-performance conjugated polymers for photovoltaic applications. Here, we use regioregular poly(3-ethylhexylthiophene) (P3EHT) and poly(3-ethylhexyl-4-fluorothiophene) (F-P3EHT) as simplified model materials, using single-molecule/aggregate spectroscopy and molecular dynamic simulations, to elucidate the impacts of backbone fluorination on morphology and excitonic coupling on the molecular scale. Despite its high regioregularity, regioregular P3EHT exhibits a rather broad distribution in polymer chain conformation due to the strong steric hindrance of bulky ethylhexyl side chains. This conformational variability results in disordered interchain morphology even between a few chains, prohibiting long-range effective interchain coupling. In stark contrast, the experimental and molecular dynamic calculations reveal that backbone fluorination of F-P3EHT leads to an extended rod-like single-chain conformation and hence highly ordered interchain packing in aggregates. Surprisingly, the ordered and close interchain packing in F-P3EHT does not lead to strong excitonic coupling between the chains but rather to dominant intrachain excitonic coupling that greatly reduces the molecular energetic heterogeneity.

Fluorescent Proteins Detect Host Structural Rearrangements via Electrostatic Mechanism.

The rational design of genetically encoded fluorescent biosensors, which can detect rearrangements of target proteins via interdomain allostery, is hindered by the absence of mechanistic understanding of the underlying photophysics. Here, we focus on the modulation of fluorescence by mechanical perturbation in a popular biological probe: enhanced Green Fluorescent Protein (eGFP). Using a combination of molecular dynamics (MD) simulations and quantum chemistry, and a set of physically motivated assumptions, we construct a map of fluorescence quantum yield as a function of a 2D electric field imposed by the protein environment on the fluorophore. This map is transferable between Tsien's Class 2 GFP's, and it allows one to estimate the shifts in fluorescence intensity due to mechanical perturbations directly from MD simulations. We use it in combination with steered MD simulations to put forward a hypothesis for the mechanism of a genetically encoded voltage probe (ArcLight) whose mechanism is currently under debate.

Isotope effects in aqueous solvation of simple halides.

We present a path-integral-molecular-dynamics study of the thermodynamic stabilities of DOH⋯ X- and HOD⋯ X- (X = F, Cl, Br, I) coordination in aqueous solutions at ambient conditions. In agreement with experimental evidence, our results for the F- case reveal a clear stabilization of the latter motif, whereas, in the rest of the halogen series, the former articulation prevails. The DOH⋯ X- preference becomes more marked the larger the size of the ionic solute. A physical interpretation of these tendencies is provided in terms of an analysis of the global quantum kinetic energies of the light atoms and their geometrical decomposition. The stabilization of the alternative ionic coordination geometries is the result of a delicate balance arising from quantum spatial dispersions along parallel and perpendicular directions with respect to the relevant O-H⋯X- axis, as the strength of the water-halide H-bond varies. This interpretation is corroborated by a complementary analysis performed on the different spectroscopic signals of the corresponding IR spectra.

Isotopic equilibria in aqueous clusters at low temperatures: Insights from the MB-pol many-body potential.

By combining path-integrals molecular dynamics simulations with the accurate MB-pol potential energy surface, we investigate the role of alternative potential models on isotopic fractionation ratios between H and D atoms at dangling positions in water clusters at low temperatures. Our results show clear stabilizations of the lighter isotope at dangling sites, characterized by free energy differences ΔG that become comparable to or larger than kBT for temperatures below ∼75 K. The comparison between these results to those previously reported using the empirical q-TIP4P/F water model [P. E. Videla et al., J. Phys. Chem. Lett. 5, 2375 (2014)] reveals that the latter Hamiltonian overestimates the H stabilization by ∼25%. Moreover, predictions from the MB-pol model are in much better agreement with measured results reported for similar isotope equilibria at ice surfaces. The dissection of the quantum kinetic energies into orthogonal directions shows that the dominant differences between the two models are to be found in the anharmonic characteristics of the potential energy surfaces along OH bond directions involved in hydrogen bonds.

Communication: Isotopic effects on tunneling motions in the water trimer.

We present results of ring polymer molecular dynamics simulations that shed light on the effects of nuclear quantum fluctuations on tunneling motions in cyclic [H2O]3 and [D2O]3, at the representative temperature of T = 75 K. In particular, we focus attention on free energies associated with two key isomerization processes: The first one corresponds to flipping transitions of dangling OH bonds, between up and down positions with respect to the O-O-O plane of the cluster; the second involves the interchange between connecting and dangling hydrogen bond character of the H-atoms in a tagged water molecule. Zero point energy and tunneling effects lead to sensible reductions of the free energy barriers. Due to the lighter nature of the H nuclei, these modifications are more marked in [H2O]3 than in [D2O]3. Estimates of the characteristic time scales describing the flipping transitions are consistent with those predicted based on standard transition-state-approximation arguments.

Dramatic vapor-phase modulation of the characteristics of graphene field-effect transistors.

Here we report on dramatic and favorable changes to the operating characteristics in monolayer graphene field-effect transistors (FETs) exposed to vapor-phase, polar organic molecules in ambient. These changes include significant reduction of the Dirac voltage, accompanied by both an increase in electron and hole mobility, µ, and a decrease in residual carrier density, N0, to < 3 × 10(11) cm(-2). In contrast to graphene FET modulation with various liquid- and solid-phase dielectric media present in the literature, we attribute these changes to screening by polar vapor-phase molecules of fields induced by charged impurities and defects, n(imp), in or near the active layer. The magnitude of the changes produced in the graphene FET parameters scales remarkably well with the dipole moment of the delivered molecules. These effects are reversible, a unique advantage of working in the vapor phase. The changes observed upon polar molecule delivery are analogous to those produced by depositing and annealing fluoropolymer coatings on graphene that have been reported previously, and we attribute these changes to similar charge screening or neutralization phenomena.

A new class of ensemble conserving algorithms for approximate quantum dynamics: Theoretical formulation and model problems.

We develop two classes of quasi-classical dynamics that are shown to conserve the initial quantum ensemble when used in combination with the Feynman-Kleinert approximation of the density operator. These dynamics are used to improve the Feynman-Kleinert implementation of the classical Wigner approximation for the evaluation of quantum time correlation functions known as Feynman-Kleinert linearized path-integral. As shown, both classes of dynamics are able to recover the exact classical and high temperature limits of the quantum time correlation function, while a subset is able to recover the exact harmonic limit. A comparison of the approximate quantum time correlation functions obtained from both classes of dynamics is made with the exact results for the challenging model problems of the quartic and double-well potentials. It is found that these dynamics provide a great improvement over the classical Wigner approximation, in which purely classical dynamics are used. In a special case, our first method becomes identical to centroid molecular dynamics.

Application of a new ensemble conserving quantum dynamics simulation algorithm to liquid para-hydrogen and ortho-deuterium.

We apply the Feynman-Kleinert Quasi-Classical Wigner (FK-QCW) method developed in our previous work [Smith et al., J. Chem. Phys. 142, 244112 (2015)] for the determination of the dynamic structure factor of liquid para-hydrogen and ortho-deuterium at state points of (T = 20.0 K, n = 21.24 nm(-3)) and (T = 23.0 K, n = 24.61 nm(-3)), respectively. When applied to this challenging system, it is shown that this new FK-QCW method consistently reproduces the experimental dynamic structure factor reported by Smith et al. [J. Chem. Phys. 140, 034501 (2014)] for all momentum transfers considered. This shows that FK-QCW provides a substantial improvement over the Feynman-Kleinert linearized path-integral method, in which purely classical dynamics are used. Furthermore, for small momentum transfers, it is shown that FK-QCW provides nearly the same results as ring-polymer molecular dynamics (RPMD), thus suggesting that FK-QCW provides a potentially more appealing algorithm than RPMD since it is not formally limited to correlation functions involving linear operators.

Excitonic energy migration in conjugated polymers: the critical role of interchain morphology.

Excitonic energy migration was studied using single molecule spectroscopy of individual conjugated polymer (CP) chains and aggregates. To probe the effect of interchain morphology on energy migration in CP, tailored interchain morphologies were achieved using solvent vapor annealing to construct polymer aggregates, which were then studied with single aggregate spectroscopy. We report that highly ordered interchain packing in regioregular poly(3-hexylthiophene) (rr-P3HT) enables long-range interchain energy migration, while disordered packing in regiorandom poly(3-hexylthiophene) (rra-P3HT), even in aggregates of just a few chains, can dramatically impede the interchain mechanism. In contrast to rr-P3HT, interchain energy migration in poly(3-(2'-methoxy-5'-octylphenyl)thiophene) (POMeOPT), a polythiophene derivative with bulky side chains, can be completely inhibited. We use simulated structures to show that the reduction in interchain coupling is not due simply to increased packing distance between backbones of different chains, but reflects inhibition of stacking due to side-chain-induced twisting of the contours of individual chains. A competition from intrachain coupling has also been demonstrated by comparing POMeOPT aggregates with different polymer chain sizes.

Hot charge-transfer excitons set the time limit for charge separation at donor/acceptor interfaces in organic photovoltaics.

Photocurrent generation in organic photovoltaics (OPVs) relies on the dissociation of excitons into free electrons and holes at donor/acceptor heterointerfaces. The low dielectric constant of organic semiconductors leads to strong Coulomb interactions between electron-hole pairs that should in principle oppose the generation of free charges. The exact mechanism by which electrons and holes overcome this Coulomb trapping is still unsolved, but increasing evidence points to the critical role of hot charge-transfer (CT) excitons in assisting this process. Here we provide a real-time view of hot CT exciton formation and relaxation using femtosecond nonlinear optical spectroscopies and non-adiabatic mixed quantum mechanics/molecular mechanics simulations in the phthalocyanine-fullerene model OPV system. For initial excitation on phthalocyanine, hot CT excitons are formed in 10(-13) s, followed by relaxation to lower energies and shorter electron-hole distances on a 10(-12) s timescale. This hot CT exciton cooling process and collapse of charge separation sets the fundamental time limit for competitive charge separation channels that lead to efficient photocurrent generation.

Refinement of the experimental dynamic structure factor for liquid para-hydrogen and ortho-deuterium using semi-classical quantum simulation.

The dynamic structure factor of liquid para-hydrogen and ortho-deuterium in corresponding thermodynamic states (T = 20.0 K, n = 21.24 nm(-3)) and (T = 23.0 K, n = 24.61 nm(-3)), respectively, has been computed by both the Feynman-Kleinert linearized path-integral (FK-LPI) and Ring-Polymer Molecular Dynamics (RPMD) methods and compared with Inelastic X Ray Scattering spectra. The combined use of computational and experimental methods enabled us to reduce experimental uncertainties in the determination of the true sample spectrum. Furthermore, the refined experimental spectrum of para-hydrogen and ortho-deuterium is consistently reproduced by both FK-LPI and RPMD results at momentum transfers lower than 12.8 nm(-1). At larger momentum transfers the FK-LPI results agree with experiment much better for ortho-deuterium than for para-hydrogen. More specifically we found that for k ∼ 20.0 nm(-1) para-hydrogen provides a test case for improved approximations to quantum dynamics.