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.

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.

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.

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.

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.

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.

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 ClNO2 product yield in competition with hydrolysis following N2O5 uptake to aqueous NaCl solutions. For NaCl-D2O solutions spanning 0.0054-0.21 M, the ClNO2 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 N2O5 in D2O 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 D2O attack on N2O5 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+400, we conclude that the strong preference of N2O5 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 N2O5. 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.

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-/Å2·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.

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.

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.

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-photochemical 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.