Unraveling the mechanisms of triplet state formation in a heavy-atom free photosensitizer

Triplet excited state generation plays a pivotal role in photosensitizers, however the reliance on transition metals and heavy atoms can limit the utility of these systems. In this study, we demonstrate that an interplay of competing quantum effects controls the high triplet quantum yield in a prototypical boron dipyrromethene-anthracene (BD-An) donor-acceptor dyad photosensitizer, which is only captured by an accurate treatment of both inner and outer sphere reorganization energies. Our ab initio-derived model provides excellent agreement with experimentally measured spectra, triplet yields and excited state kinetic data, including the triplet lifetime. We find that rapid triplet state formation occurs primarily via high-energy triplet states through both spin-orbit coupled charge transfer and El-Sayed's rule breaking intersystem crossing, rather than direct spin-orbit coupled charge transfer to the lowest lying triplet state. Our calculations also reveal that competing effects of nuclear tunneling, electronic state recrossing, and electronic polarizability dictate the rate of non-productive ground state recombination. This study sheds light on the quantum effects driving efficient triplet formation in the BD-An system, and offers a promising simulation methodology for diverse photochemical systems.

On the Statistical Mechanics of Mass Accommodation at Liquid-Vapor Interfaces

We propose a framework for describing the dynamics associated with the adsorption of small molecules to liquid-vapor interfaces using an intermediate resolution between traditional continuum theories that are bereft of molecular detail and molecular dynamics simulations that are replete with them. In particular, we develop an effective single particle equation of motion capable of describing the physical processes that determine thermal and mass accommodation probabilities. The effective equation is parametrized with quantities that vary through space away from the liquid-vapor interface. Of particular importance in describing the early time dynamics is the spatially dependent friction, for which we propose a numerical scheme to evaluate from molecular simulation. Taken together with potentials of mean force computable with importance sampling methods, we illustrate how to compute the mass accommodation coefficient and residence time distribution. Throughout, we highlight the case of ozone adsorption in aqueous solutions and its dependence on electrolyte composition.

Molecular Insights into Chemical Reactions at Aqueous Aerosol Interfaces

Atmospheric aerosols facilitate reactions between ambient gases and dissolved species. Here, we review our efforts to interrogate the uptake of these gases and the mechanisms of their reactions both theoretically and experimentally. We highlight the fascinating behavior of N2O5 in solutions ranging from pure water to complex mixtures, chosen because its aerosol-mediated reactions significantly impact global ozone, hydroxyl, and methane concentrations. As a hydrophobic, weakly soluble, and highly reactive species, N2O5 is a sensitive probe of the chemical and physical properties of aerosol interfaces. We employ contemporary theory to disentangle the fate of N2O5 as it approaches pure and salty water, starting with adsorption and ending with hydrolysis to HNO3, chlorination to ClNO2, or evaporation. Flow reactor and gas-liquid scattering experiments probe even greater complexity as added ions, organic molecules, and surfactants alter the interfacial composition and reaction rates. Together, we reveal a new perspective on multiphase chemistry in the atmosphere. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 75 is April 2024. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

In Situ X-ray Scattering Reveals Coarsening Rates of Superlattices Self-Assembled from Electrostatically Stabilized Metal Nanocrystals Depend Nonmonotonically on Driving Force

Self-assembly of colloidal nanocrystals (NCs) into superlattices (SLs) is an appealing strategy to design hierarchically organized materials with promising functionalities. Mechanistic studies are still needed to uncover the design principles for SL self-assembly, but such studies have been difficult to perform due to the fast time and short length scales of NC systems. To address this challenge, we developed an apparatus to directly measure the evolving phases in situ and in real time of an electrostatically stabilized Au NC solution before, during, and after it is quenched to form SLs using small-angle X-ray scattering. By developing a quantitative model, we fit the time-dependent scattering patterns to obtain the phase diagram of the system and the kinetics of the colloidal and SL phases as a function of varying quench conditions. The extracted phase diagram is consistent with particles whose interactions are short in range relative to their diameter. We find the degree of SL order is primarily determined by fast (subsecond) initial nucleation and growth kinetics, while coarsening at later times depends nonmonotonically on the driving force for self-assembly. We validate these results by direct comparison with simulations and use them to suggest dynamic design principles to optimize the crystallinity within a finite time window. The combination of this measurement methodology, quantitative analysis, and simulation should be generalizable to elucidate and better control the microscopic self-assembly pathways of a wide range of bottom-up assembled systems and architectures.

Iodide oxidation by ozone at the surface of aqueous microdroplets

The oxidation of iodide by ozone occurs at the sea-surface and within sea spray aerosol, influencing the overall ozone budget in the marine boundary layer and leading to the emission of reactive halogen gases. A detailed account of the surface mechanism has proven elusive, however, due to the difficulty in quantifying multiphase kinetics. To obtain a clearer understanding of this reaction mechanism at the air-water interface, we report pH-dependent oxidation kinetics of I- in single levitated microdroplets as a function of [O3] using a quadrupole electrodynamic trap and an open port sampling interface for mass spectrometry. A kinetic model, constrained by molecular simulations of O3 dynamics at the air-water interface, is used to understand the coupled diffusive, reactive, and evaporative pathways at the microdroplet surface, which exhibit a strong dependence on bulk solution pH. Under acidic conditions, the surface reaction is limited by O3 diffusion in the gas phase, whereas under basic conditions the reaction becomes rate limited on the surface. The pH dependence also suggests the existence of a reactive intermediate IOOO- as has previously been observed in the Br- + O3 reaction. Expressions for steady-state surface concentrations of reactants are derived and utilized to directly compute uptake coefficients for this system, allowing for an exploration of uptake dependence on reactant concentration. In the present experiments, reactive uptake coefficients of O3 scale weakly with bulk solution pH, increasing from 4 × 10-4 to 2 × 10-3 with decreasing solution pH from pH 13 to pH 3.

Kinetics of the xanthophyll cycle and its role in photoprotective memory and response

Efficiently balancing photochemistry and photoprotection is crucial for survival and productivity of photosynthetic organisms in the rapidly fluctuating light levels found in natural environments. The ability to respond quickly to sudden changes in light level is clearly advantageous. In the alga Nannochloropsis oceanica we observed an ability to respond rapidly to sudden increases in light level which occur soon after a previous high-light exposure. This ability implies a kind of memory. In this work, we explore the xanthophyll cycle in N. oceanica as a short-term photoprotective memory system. By combining snapshot fluorescence lifetime measurements with a biochemistry-based quantitative model, we show that short-term memory arises from the xanthophyll cycle. In addition, the model enables us to characterize the relative quenching abilities of the three xanthophyll cycle components. Given the ubiquity of the xanthophyll cycle in photosynthetic organisms the model described here will be of utility in improving our understanding of vascular plant and algal photoprotection with important implications for crop productivity.

Nonequilibrium design strategies for functional colloidal assemblies

We use a nonequilibrium variational principle to optimize the steady-state, shear-induced interconversion of self-assembled nanoclusters of DNA-coated colloids. Employing this principle within a stochastic optimization algorithm allows us to identify design strategies for functional materials. We find that far-from-equilibrium shear flow can significantly enhance the flux between specific colloidal states by decoupling trade-offs between stability and reactivity required by systems in equilibrium. For isolated nanoclusters, we find nonequilibrium strategies for amplifying transition rates by coupling a given reaction coordinate to the background shear flow. We also find that shear flow can be made to selectively break detailed balance and maximize probability currents by coupling orientational degrees of freedom to conformational transitions. For a microphase consisting of many nanoclusters, we study the flux of colloids hopping between clusters. We find that a shear flow can amplify the flux without a proportional compromise on the microphase structure. This approach provides a general means of uncovering design principles for nanoscale, autonomous, functional materials driven far from equilibrium.

Structural evolution of water-in-propylene carbonate mixtures revealed by polarized Raman spectroscopy and molecular dynamics

The liquid structure of systems wherein water is limited in concentration or through geometry is of great interest in various fields such as biology, materials science, and electrochemistry. Here, we present a combined polarized Raman and molecular dynamics investigation of the structural changes that occur as water is added incrementally to propylene carbonate (PC), a polar, aprotic solvent that is important in lithium-ion batteries. Polarized Raman spectra of PC solutions were collected for water mole fractions 0.003 ≤ χwater ≤ 0.296, which encompasses the solubility range of water in PC. The novel approach taken herein provides additional hydrogen bond and solvation characterization of this system that has not been achievable in previous studies. Analysis of the polarized carbonyl Raman band in conjunction with simulations demonstrated that the bulk structure of the solvent remained unperturbed upon the addition of water. Experimental spectra in the O-H stretching region were decomposed through Gaussian fitting into sub-bands and comparison to studies of dilute HOD in D2O. With the aid of simulations, we identified these different bands as water arrangements having different degrees of hydrogen bonding. The observed water structure within PC indicates that water tends to self-aggregate, forming a hydrogen bond network that is distinctly different from the bulk and dependent on concentration. For example, at moderate concentrations, the most likely aggregate structures are chains of water molecules, each with two hydrogen bonds.

On the Mechanism of Polaritonic Rate Suppression from Quantum Transition Paths

Polariton chemistry holds promise for facilitating mode-selective chemical reactions, but the underlying mechanism behind the rate modifications observed under strong vibrational coupling is not well understood. Using the recently developed quantum transition path theory, we have uncovered a mechanism of resonant suppression of a thermal reaction rate in a simple model polaritonic system consisting of a reactive mode in a bath confined to a lossless microcavity with a single photon mode. We observed the formation of a polariton during rate-limiting transitions on reactive pathways and identified the concomitant rate suppression as being due to hybridization between the reactive mode and the cavity mode, which inhibits bath-mediated tunneling. The transition probabilities that define the quantum master equation can be directly translated into a visualization of the corresponding polariton energy landscape. This landscape exhibits a double funnel structure with a large barrier between the initial and final states.

Variational deep learning of equilibrium transition path ensembles

We present a time-dependent variational method to learn the mechanisms of equilibrium reactive processes and efficiently evaluate their rates within a transition path ensemble. This approach builds off of the variational path sampling methodology by approximating the time-dependent commitment probability within a neural network ansatz. The reaction mechanisms inferred through this approach are elucidated by a novel decomposition of the rate in terms of the components of a stochastic path action conditioned on a transition. This decomposition affords an ability to resolve the typical contribution of each reactive mode and their couplings to the rare event. The associated rate evaluation is variational and systematically improvable through the development of a cumulant expansion. We demonstrate this method in both over- and under-damped stochastic equations of motion, in low-dimensional model systems, and in the isomerization of a solvated alanine dipeptide. In all examples, we find that we can obtain quantitatively accurate estimates of the rates of the reactive events with minimal trajectory statistics and gain unique insights into transitions through the analysis of their commitment probability.

Size-Dependent Lattice Symmetry Breaking Determines the Exciton Fine Structure of Perovskite Nanocrystals

The order of bright and dark excitonic states in lead-halide perovskite nanocrystals is debated. It has been proposed that the Rashba effect, driven by lattice-induced symmetry breaking, causes a bright excitonic ground state. Direct measurements of excitonic spectra, however, show the signatures of a dark ground state, bringing the role of the Rashba effect into question. We use an atomistic theory to model the exciton fine structure of perovskite nanocrystals, accounting for realistic lattice distortions. We calculate optical gaps and excitonic features that compare favorably with experimental works. The exciton fine structure splittings show a nonmonotonic size dependence due to a structural transition between cubic and orthorhombic phases. Additionally, the excitonic ground state is found to be dark with spin triplet character, exhibiting a small Rashba coupling. We additionally explore the effects of nanocrystal shape on the fine structure, clarifying observations on polydisperse nanocrystals.

Spin selective charge recombination in chiral donor-bridge-acceptor triads

In this paper, we outline a physically motivated framework for describing spin-selective recombination processes in chiral systems, from which we derive spin-selective reaction operators for recombination reactions of donor-bridge-acceptor molecules, where the electron transfer is mediated by chirality and spin-orbit coupling. In general, the recombination process is selective only for spin-coherence between singlet and triplet states, and it is not, in general, selective for spin polarization. We find that spin polarization selectivity only arises in hopping-mediated electron transfer. We describe how this effective spin-polarization selectivity is a consequence of spin-polarization generated transiently in the intermediate state. The recombination process also augments the coherent spin dynamics of the charge separated state, which is found to have a significant effect on the recombination dynamics and to destroy any long-lived spin polarization. Although we only consider a simple donor-bridge-acceptor system, the framework we present here can be straightforwardly extended to describe spin-selective recombination processes in more complex systems.

Intrinsic Interface Adsorption Drives Selectivity in Atomically Smooth Nanofluidic Channels

Specific molecular interactions underlie unexpected and useful phenomena in nanofluidic systems, but these require descriptions that go beyond traditional macroscopic hydrodynamics. In this letter, we demonstrate how equilibrium molecular dynamics simulations and linear response theory can be synthesized with hydrodynamics to provide a comprehensive characterization of nanofluidic transport. Specifically, we study the pressure driven flows of ionic solutions in nanochannels comprised of two-dimensional crystalline substrates made from graphite and hexagonal boron nitride. While simple hydrodynamic descriptions do not predict a streaming electrical current or salt selectivity in such simple systems, we observe that both arise due to the intrinsic molecular interactions that act to selectively adsorb ions to the interface in the absence of a net surface charge. Notably, this emergent selectivity indicates that these nanochannels can serve as desalination membranes.

Inferring equilibrium transition rates from nonequilibrium protocols

We develop a theory for inferring equilibrium transition rates from trajectories driven by a time-dependent force using results from stochastic thermodynamics. Applying the Kawasaki relation to approximate the nonequilibrium distribution function in terms of the equilibrium distribution function and the excess dissipation, we formulate a nonequilibrium transition state theory to estimate the rate enhancement over the equilibrium rate due to the nonequilibrium protocol. We demonstrate the utility of our theory in examples of pulling of harmonically trapped particles in one and two dimensions, as well as a semiflexible polymer with a reactive linker in three dimensions. We expect our purely thermodynamic approach will find use in both molecular simulation and force spectroscopy experiments.

Odd Mobility of a Passive Tracer in a Chiral Active Fluid

Chiral active fluids break both time-reversal and parity symmetry, leading to exotic transport phenomena unobservable in ordinary passive fluids. We develop a generalized Green-Kubo relation for the anomalous lift experienced by a passive tracer suspended in a two-dimensional chiral active fluid subjected to an applied force. This anomalous lift is characterized by a transport coefficient termed the odd mobility. We validate our generalized response theory using molecular dynamics simulations, and we show that the asymmetric tracer mobility may be understood mechanically in terms of asymmetric deformations of the tracer-fluid density distribution function. We show that the even and odd components of the mobility decay at different rates with tracer size, suggesting the possibility of size-based particle separation using a chiral active working fluid.

Direct Observation of Transient Structural Dynamics of Atomically Thin Halide Perovskite Nanowires

Halide perovskite is a unique dynamical system, whose structural and chemical processes happening across different timescales have significant impact on its physical properties and device-level performance. However, due to its intrinsic instability, real-time investigation of the structure dynamics of halide perovskite is challenging, which hinders the systematic understanding of the chemical processes in the synthesis, phase transition, and degradation of halide perovskite. Here, we show that atomically thin carbon materials can stabilize ultrathin halide perovskite nanostructures against otherwise detrimental conditions. Moreover, the protective carbon shells enable atomic-level visualization of the vibrational, rotational, and translational movement of halide perovskite unit cells. Albeit atomically thin, protected halide perovskite nanostructures can maintain their structural integrity up to an electron dose rate of 10,000 e^-/Ų·s while exhibiting unusual dynamical behaviors pertaining to the lattice anharmonicity and nanoscale confinement. Our work demonstrates an effective method to protect beam-sensitive materials during in situ observation, unlocking new solutions to study new modes of structure dynamics of nanomaterials.

Weak Temperature Dependence of the Relative Rates of Chlorination and Hydrolysis of N2O5 in NaCl-Water Solutions

We have measured the temperature dependence of the ClNO₂ product yield in competition with hydrolysis following N₂O₅ uptake to aqueous NaCl solutions. For NaCl-D₂O solutions spanning 0.0054-0.21 M, the ClNO₂ product yield decreases on average by only 4 ± 3% from 5 to 25 °C. Less reproducible measurements at 0.54-2.4 M NaCl also fall within this range. The ratio of the rate constants for chlorination and hydrolysis of N₂O₅ in D₂O is determined on average to be 1150 ± 90 at 25 °C up to 0.21 M NaCl, favoring chlorination. This ratio is observed to decrease significantly at the two highest concentrations. An Arrhenius analysis reveals that the activation energy for hydrolysis is just 3.0 ± 1.5 kJ/mol larger than for chlorination up to 0.21 M, indicating that Cl^- and D₂O attack on N₂O₅ has similar energetic barriers despite the differences in charge and complexity of these reactants. In combination with the measured preexponential ratio favoring chlorination of 300_-200⁺⁴⁰⁰, we conclude that the strong preference of N₂O₅ to undergo chlorination over hydrolysis is driven by dynamic and entropic, rather than enthalpic, factors. Molecular dynamics simulations elucidate the distinct solvation between strongly hydrated Cl^- and the hydrophobically solvated N₂O₅. Combining this molecular picture with the Arrhenius analysis implicates the role of water in mediating interactions between such distinctly solvated species and suggests a role for diffusion limitations on the chlorination reaction.

Kinetic frustration by limited bond availability controls the LAT protein condensation phase transition on membranes

LAT is a membrane-linked scaffold protein that undergoes a phase transition to form a two-dimensional protein condensate on the membrane during T cell activation. Governed by tyrosine phosphorylation, LAT recruits various proteins that ultimately enable condensation through a percolation network of discrete and selective protein-protein interactions. Here, we describe detailed kinetic measurements of the phase transition, along with coarse-grained model simulations, that reveal that LAT condensation is kinetically frustrated by the availability of bonds to form the network. Unlike typical miscibility transitions in which compact domains may coexist at equilibrium, the LAT condensates are dynamically arrested in extended states, kinetically trapped out of equilibrium. Modeling identifies the structural basis for this kinetic arrest as the formation of spindle arrangements, favored by limited multivalent binding interactions along the flexible, intrinsically disordered LAT protein. These results reveal how local factors controlling the kinetics of LAT condensation enable formation of different, stable condensates, which may ultimately coexist within the cell.

Nonadiabatic transition paths from quantum jump trajectories

We present a means of studying rare reactive pathways in open quantum systems using transition path theory and ensembles of quantum jump trajectories. This approach allows for the elucidation of reactive paths for dissipative, nonadiabatic dynamics when the system is embedded in a Markovian environment. We detail the dominant pathways and rates of thermally activated processes and the relaxation pathways and photoyields following vertical excitation in a minimal model of a conical intersection. We find that the geometry of the conical intersection affects the electronic character of the transition state as defined through a generalization of a committor function for a thermal barrier crossing event. Similarly, the geometry changes the mechanism of relaxation following a vertical excitation. Relaxation in models resulting from small diabatic coupling proceeds through pathways dominated by pure dephasing, while those with large diabatic coupling proceed through pathways limited by dissipation. The perspective introduced here for the nonadiabatic dynamics of open quantum systems generalizes classical notions of reactive paths to fundamentally quantum mechanical processes.

Coupled charge and energy transfer dynamics in light harvesting complexes from a hybrid hierarchical equations of motion approach

We describe a method for simulating exciton dynamics in protein-pigment complexes, including effects from charge transfer as well as fluorescence. The method combines the hierarchical equations of motion, which are used to describe quantum dynamics of excitons, and the Nakajima-Zwanzig quantum master equation, which is used to describe slower charge transfer processes. We study the charge transfer quenching in light harvesting complex II, a protein postulated to control non-photochemcial quenching in many plant species. Using our hybrid approach, we find good agreement between our calculation and experimental measurements of the excitation lifetime. Furthermore our calculations reveal that the exciton energy funnel plays an important role in determining quenching efficiency, a conclusion we expect to extend to other proteins that perform protective excitation quenching. This also highlights the need for simulation methods that properly account for the interplay of exciton dynamics and charge transfer processes.

Renormalization of excitonic properties by polar phonons

We employ quasiparticle path integral molecular dynamics to study how theexcitonic properties of model semiconductors are altered by electron-phononcoupling. We describe ways within a path integral representation of the systemto evaluate the renormalized mass, binding energy, and radiative recombinationrate of excitons in the presence of a fluctuating lattice. To illustrate thisapproach, we consider Fr\"ohlich-type electron-phonon interactions and employan imaginary time influence functional to incorporate phonon-induced effectswithout approximation. The effective mass and binding energies are comparedwith perturbative and variational approaches, which provide qualitativelyconsistent trends. We evaluate electron-hole recombination rates as mediatedthrough both trap-assisted and bimolecular processes, developing a consistentstatistical mechanical approach valid in the reaction limited regime. Thesecalculations demonstrate how phonons screen electron-hole interactions,generically reducing exciton binding energies and increasing their radiativelifetimes.

Xanthophyll-cycle based model of the rapid photoprotection of Nannochloropsis in response to regular and irregular light/dark sequences

Abstract <p>We explore the photoprotection dynamics of Nannochloropsis oceanica using time-correlated single photon counting under regular and irregular actinic light sequences. The varying light sequences mimic natural conditions, allowing us to probe the real-time response of non-photochemical quenching (NPQ) pathways. Durations of fluctuating light exposure during a fixed total experimental time and prior light exposure of the algae are both found to have a profound effect on NPQ. These observations are rationalized with a quantitative model based on the xanthophyll cycle and the protonation of LHCX1. The model is able to accurately describe the dynamics of non-photochemical quenching across a variety of light sequences. The combined model and observations suggest that the accumulation of a quenching complex, likely zeaxanthin bound to a protonated LHCX1, is responsible for the gradual rise in NPQ. Additionally, the model makes specific predictions for the light sequence dependence of xanthophyll concentrations that are in reasonable agreement with independent chromatography measurements taken during a specific light/dark sequence.

Defect-mediated ripening of core-shell nanostructures

Understanding nanostructure ripening mechanisms is desirable for gaining insight on the growth and potential applications of nanoscale materials. However, the atomic pathways of nanostructure ripening in solution have rarely been observed directly. Here, we report defect-mediated ripening of Cd-CdCl₂ core-shell nanoparticles (CSN) revealed by in-situ atomic resolution imaging with liquid cell transmission electron microscopy. We find that ripening is initiated by dissolution of the nanoparticle with an incomplete CdCl₂ shell, and that the areas of the Cd core that are exposed to the solution are etched first. The growth of the other nanoparticles is achieved by generating crack defects in the shell, followed by ion diffusion through the cracks. Subsequent healing of crack defects leads to a highly crystalline CSN. The formation and annihilation of crack defects in the CdCl₂ shell, accompanied by disordering and crystallization of the shell structure, mediate the ripening of Cd-CdCl₂ CSN in the solution.

Understanding the ripening of core-shell nanostructures is challenging. Here, the authors use liquid cell transmission electron microscopy to show that the atomic ripening pathway for Cd-CdCl₂ core-shell nanoparticles is mediated by crack defects.

Self-assembly of nanocrystals into strongly electronically coupled all-inorganic supercrystals

Colloidal nanocrystals of metals, semiconductors, and other functional materials can self-assemble into long-range ordered crystalline and quasicrystalline phases, but insulating organic surface ligands prevent the development of collective electronic states in ordered nanocrystal assemblies. We reversibly self-assembled colloidal nanocrystals of gold, platinum, nickel, lead sulfide, and lead selenide with conductive inorganic ligands into supercrystals exhibiting optical and electronic properties consistent with strong electronic coupling between the constituent nanocrystals. The phase behavior of charge-stabilized nanocrystals can be rationalized and navigated with phase diagrams computed for particles interacting through short-range attractive potentials. By finely tuning interparticle interactions, the assembly was directed either through one-step nucleation or nonclassical two-step nucleation pathways. In the latter case, the nucleation was preceded by the formation of two metastable colloidal fluids.

Nonlocal Screening Dictates the Radiative Lifetimes of Excitations in Lead Halide Perovskites

We use path integral molecular dynamics simulations and theory to elucidate the interactions between charge carriers, as mediated by a lead halide perovskite lattice. We find that the charge-lattice coupling of MAPbI₃ results in a repulsive interaction between electrons and holes at intermediate distances. The effective interaction is understood using a Gaussian field theory, whereby the underlying soft, polar lattice contributes a nonlocal screening between quasiparticles. Path integral calculations of this nonlocal screening model are used to rationalize the small exciton binding energy and low radiative recombination rate observed experimentally and are compared to traditional Wannier-Mott and Fröhlich models, which fail to do so. These results clarify the origin of the high power conversion efficiencies in lead halide perovskites. Emergent repulsive electron-hole interactions provide a design principle for optimizing soft, polar semiconductors.

Peptide Isomerization is Suppressed at the Air-Water Interface

We use molecular dynamics simulations to study the thermodynamics and kinetics of alanine dipeptide isomerization at the air-water interface. Thermodynamically, we find an affinity of the dipeptide to the interface. This affinity arises from stabilizing intramolecular interactions that become unshielded as the dipeptide is desolvated. Kinetically, we consider the rate of transitions between the α_L and β conformations of alanine dipeptide and evaluate it as a continuous function of the distance from the interface using a recent extension of transition path sampling, TPS+U. The rate of isomerization at the Gibbs dividing surface is suppressed relative to the bulk by a factor of 3. Examination of the ensemble of transition states elucidates the role of solvent degrees of freedom in mediating favorable intramolecular interactions along the reaction pathway of isomerization. Near the air-water interface, water is less effective at mediating these intramolecular interactions.

Direct Evaluation of Rare Events in Active Matter from Variational Path Sampling

Active matter represents a broad class of systems that evolve far from equilibrium due to the local injection of energy. Like their passive analogs, transformations between distinct metastable states in active matter proceed through rare fluctuations; however, their detailed balance violating dynamics renders these events difficult to study. Here, we present a simulation method for evaluating the rate and mechanism of rare events in generic nonequilibrium systems and apply it to study the conformational changes of a passive solute in an active fluid. The method employs a variational optimization of a control force that renders the rare event a typical one, supplying an exact estimate of its rate as a ratio of path partition functions. Using this method we find that increasing activity in the active bath can enhance the rate of conformational switching of the passive solute in a manner consistent with recent bounds from stochastic thermodynamics.

Learning intermolecular forces at liquid-vapor interfaces

By adopting a perspective informed by contemporary liquid-state theory, we consider how to train an artificial neural network potential to describe inhomogeneous, disordered systems. We find that neural network potentials based on local representations of atomic environments are capable of describing some properties of liquid-vapor interfaces but typically fail for properties that depend on unbalanced long-ranged interactions that build up in the presence of broken translation symmetry. These same interactions cancel in the translationally invariant bulk, allowing local neural network potentials to describe bulk properties correctly. By incorporating explicit models of the slowly varying long-ranged interactions and training neural networks only on the short-ranged components, we can arrive at potentials that robustly recover interfacial properties. We find that local neural network models can sometimes approximate a local molecular field potential to correct for the truncated interactions, but this behavior is variable and hard to learn. Generally, we find that models with explicit electrostatics are easier to train and have higher accuracy. We demonstrate this perspective in a simple model of an asymmetric dipolar fluid, where the exact long-ranged interaction is known, and in an ab initio water model, where it is approximated.

Nonuniqueness of generalized quantum master equations for a single observable

When deriving exact generalized master equations for the evolution of a reduced set of degrees of freedom, one is free to choose what quantities are relevant by specifying projection operators. However, obtaining a reduced description does not always need to be achieved through projections-one can also use conservation laws for this purpose. Such an operation should be considered as distinct from any kind of projection; that is, projection onto a single observable yields a different form of master equation compared to that resulting from a projection followed by the application of a constraint. We give a simple example to show this point and give relationships that the different memory kernels must satisfy to yield the same dynamics.

Reinforcement learning of rare diffusive dynamics

We present a method to probe rare molecular dynamics trajectories directly using reinforcement learning. We consider trajectories that are conditioned to transition between regions of configuration space in finite time, such as those relevant in the study of reactive events, and trajectories exhibiting rare fluctuations of time-integrated quantities in the long time limit, such as those relevant in the calculation of large deviation functions. In both cases, reinforcement learning techniques are used to optimize an added force that minimizes the Kullback-Leibler divergence between the conditioned trajectory ensemble and a driven one. Under the optimized added force, the system evolves the rare fluctuation as a typical one, affording a variational estimate of its likelihood in the original trajectory ensemble. Low variance gradients employing value functions are proposed to increase the convergence of the optimal force. The method we develop employing these gradients leads to efficient and accurate estimates of both the optimal force and the likelihood of the rare event for a variety of model systems.

Distinct Chemistries Explain Decoupling of Slip and Wettability in Atomically Smooth Aqueous Interfaces

Despite essentially identical crystallography and equilibrium structuring of water, nanoscopic channels composed of hexagonal boron nitride and graphite exhibit an order-of-magnitude difference in fluid slip. We investigate this difference using molecular dynamics simulations, demonstrating that its origin is in the distinct chemistries of the two materials. In particular, the presence of polar bonds in hexagonal boron nitride, absent in graphite, leads to Coulombic interactions between the polar water molecules and the wall. We demonstrate that this interaction is manifested in a large typical lateral force experienced by a layer of oriented hydrogen atoms in the vicinity of the wall, leading to the enhanced friction in hexagonal boron nitride. The fluid adhesion to the wall is dominated by dispersive forces in both materials, leading to similar wettabilities. Our results rationalize recent observations that the difference in frictional characteristics of graphite and hexagonal boron nitride cannot be explained on the basis of the minor differences in their wettabilities.

Origin of Chirality Induced Spin Selectivity in Photoinduced Electron Transfer

Here we propose a mechanism by which spin-polarization can be generated dynamically in chiral molecular systems undergoing photoinduced electron transfer. The proposed mechanism explains how spin-polarization emerges in systems where charge transport is dominated by incoherent hopping, mediated by spin-orbit and electronic exchange couplings through an intermediate charge transfer state. We derive a simple expression for the spin-polarization that predicts a nonmonotonic temperature dependence, consistent with recent experiments, and a maximum spin-polarization that is independent of the magnitude of the spin-orbit coupling. We validate this theory using approximate quantum master equations and the numerically exact hierarchical equations of motion. The proposed mechanism of chirality induced spin selectivity should apply to many chiral systems, and the ideas presented here have implications for the study of spin transport at temperatures relevant to biology and provide simple principles for the molecular control of spins in fluctuating environments.

Molecular Properties and Chemical Transformations Near Interfaces

The properties of bulk water and aqueous solutions are known to change in the vicinity of an interface and/or in a confined environment, including the thermodynamics of ion selectivity at interfaces, transition states and pathways of chemical reactions, and nucleation events and phase growth. Here we describe joint progress in identifying unifying concepts about how air, liquid, and solid interfaces can alter molecular properties and chemical reactivity compared to bulk water and multicomponent solutions. We also discuss progress made in interfacial chemistry through advancements in new theory, molecular simulation, and experiments.

Vibronic coupling in energy transfer dynamics and two-dimensional electronic-vibrational spectra

We introduce a heterodimer model in which multiple mechanisms of vibronic coupling and their impact on energy transfer can be explicitly studied. We consider vibronic coupling that arises through either Franck-Condon activity in which each site in the heterodimer has a local electron-phonon coupling or Herzberg-Teller activity in which the transition dipole moment coupling the sites has an explicit vibrational mode-dependence. We have computed two-dimensional electronic-vibrational (2DEV) spectra for this model while varying the magnitude of these two effects and find that 2DEV spectra contain static and dynamic signatures of both types of vibronic coupling. Franck-Condon activity emerges through a change in the observed excitonic structure, while Herzberg-Teller activity is evident in the appearance of significant side-band transitions that mimic the lower-energy excitonic structure. A comparison of quantum beating patterns obtained from analysis of the simulated 2DEV spectra shows that this technique can report on the mechanism of energy transfer, elucidating a means of experimentally determining the role of specific vibronic coupling mechanisms in such processes.

On the molecular correlations that result in field-dependent conductivities in electrolyte solutions

Employing recent advances in response theory and nonequilibrium ensemble reweighting, we study the dynamic and static correlations that give rise to an electric field-dependent ionic conductivity in electrolyte solutions. We consider solutions modeled with both implicit and explicit solvents, with different dielectric properties, and at multiple concentrations. Implicit solvent models at low concentrations and small dielectric constants exhibit strongly field-dependent conductivities. We compare these results to Onsager-Wilson theory of the Wien effect, which provides a qualitatively consistent prediction at low concentrations and high static dielectric constants but is inconsistent away from these regimes. The origin of the discrepancy is found to be increased ion correlations under these conditions. Explicit solvent effects act to suppress nonlinear responses, yielding a weakly field-dependent conductivity over the range of physically realizable field strengths. By decomposing the relevant time correlation functions, we find that the insensitivity of the conductivity to the field results from the persistent frictional forces on the ions from the solvent. Our findings illustrate the utility of nonequilibrium response theory in rationalizing nonlinear transport behavior.

Entropy production fluctuations encode collective behavior in active matter

We derive a general lower bound on distributions of entropy production in interacting active matter systems. The bound is tight in the limit that interparticle correlations are small and short-ranged, which we explore in four canonical active matter models. In all models studied, the bound is weak where collective fluctuations result in long-ranged correlations, which subsequently links the locations of phase transitions to enhanced entropy production fluctuations. We develop a theory for the onset of enhanced fluctuations and relate it to specific phase transitions in active Brownian particles. We also derive optimal control forces that realize the dynamics necessary to tune dissipation and manipulate the system between phases. In so doing, we uncover a general relationship between entropy production and pattern formation in active matter, as well as ways of controlling it.

Solid-State Ionic Rectification in Perovskite Nanowire Heterostructures

Halide perovskites have attracted increasing research attention with regard to their potential for optoelectronic applications. Because of its low activation energy, ion migration is implicated in the long-term stability and many unusual transport behaviors of halide perovskite devices. However, direct observation and precise control of the ionic transport in halide perovskite crystals remain challenging. Here, we have designed an axial CsPbBr₃-CsPbCl₃ nanowire heterostructure, in which electric-field-induced halide ion migration was clearly visualized and quantified. We demonstrated that halide ion migration is dependent on the applied electric field and exhibits ionic rectification in this solid-state system, which is due to the nonuniform distribution of the ionic vacancies in the nanowire that results from a competition between electrical screening and their creation/destruction at the electrodes' interfaces. The asymmetric heterostructure characteristics add an additional knob to control the ion movement in the design of advanced ionic circuits with halide perovskites as building blocks.

Photoinduced phase separation in the lead halides is a polaronic effect

We present a perspective on recent observations of the photoinduced phase separation of halides in multi-component lead-halide perovskites. The spontaneous phase separation of an initial homogeneous solid solution under steady-state illumination conditions is found experimentally to be reversible, stochastic, weakly dependent on morphology, yet strongly dependent on composition and thermodynamic state. Regions enriched in a specific halide species that form upon phase separation are self-limiting in size, pinned to specific compositions, and grow in number in proportion to the steady-state carrier concentration until saturation. These empirical observations of robustness rule out explanations based on specific defect structures and point to the local modulation of an existing miscibility phase transition in the presence of excess charge carriers. A model for rationalizing existing observations based on the coupling between composition, strain, and charge density fluctuations through the formation of polarons is reviewed.

Field-Dependent Ionic Conductivities from Generalized Fluctuation-Dissipation Relations

We derive a relationship for the electric field dependent ionic conductivity in terms of fluctuations of time integrated microscopic variables. We demonstrate this formalism with molecular dynamics simulations of solutions of differing ionic strength with implicit solvent conditions and molten salts. These calculations are aided by a novel nonequilibrium statistical reweighting scheme that allows for the conductivity to be computed as a continuous function of the applied field. In strong electrolytes, we find the fluctuations of the ionic current are Gaussian, and subsequently, the conductivity is constant with applied field. In weaker electrolytes and molten salts, we find the fluctuations of the ionic current are strongly non-Gaussian, and the conductivity increases with applied field. This nonlinear behavior, known phenomenologically for dilute electrolytes as the Onsager-Wien effect, is general and results from the suppression of ionic correlations at large applied fields, as we elucidate through both dynamic and static correlations within nonequilibrium steady states.

Charge fluctuations from molecular simulations in the constant-potential ensemble

Statistical mechanics of constant-potential molecular simulations yields a new fluctuation–dissipation relation for the differential capacitance.

Nonlinear transport coefficients from large deviation functions

Nonlinear response occurs naturally when a strong perturbation takes a system far from equilibrium. Despite its omnipresence in nanoscale systems, it is difficult to predict in a general and efficient way. Here, we introduce a way to compute arbitrarily high order transport coefficients of stochastic systems, using the framework of large deviation theory. Leveraging time reversibility in the microscopic dynamics, we relate nonlinear response to equilibrium multitime correlation functions among both time reversal symmetric and asymmetric observables, which can be evaluated from derivatives of large deviation functions. This connection establishes a thermodynamiclike relation for nonequilibrium response and provides a practical route to its evaluation, as large deviation functions are amenable to importance sampling. We demonstrate the generality and efficiency of this method in predicting transport coefficients in single particle systems and an interacting system exhibiting thermal rectification.