Operando x-ray studies of an earth-abundant nanostructured catalyst for solar-driven water oxidation have shown that the excellent properties of this catalyst are due to the specific arrangement of Ni and Fe in the material. This catalyst offers an improved potential for scalable applications than noble metal catalysts.

The function of nanoporous electrodes involves chemical reactions coupled to transport of reactants and products into and out of the solid matrix. To understand how this coupling can depend on electrode structure and chemistry, a general multiscale, molecular-level computational framework is introduced and applied to a model system containing a nanoporous photoanode, the dye-sensitized solar cell. The pore structure and electrolyte composition are systematically varied to probe this coupling, revealing that in most cases the photocurrent density depends on dye coverage within the anode matrix, the specific architecture of an individual pore, and electrolyte gradients in between the anode and cathode. The large electrolyte gradients that develop in the anode are not kinetically controlling. The computational framework demonstrated is general and can be applied to other electrochemical, photochemical, and thermally driven nanoporous systems to evaluate the fundamental characteristics of their function.

Photovoltage and photocurrents below theoretical limits in dye-sensitized photoelectrochemical solar energy conversion systems are usually attributed to electron loss processes such as dye-electron and electrolyte-electron recombination reactions within the porous photoanode. Whether recombination is a major loss mechanism is examined here, using a multiscale reaction-diffusion computational model to evaluate system characteristics. The dye-sensitized solar cell with an I^-/I₃ ^- redox couple is chosen as a simple, representative model system because of the extensive information available for it. Two photoanode architectures with dye excitation frequencies spanning 1-25 s^-1 are examined, assuming two distinct recombination mechanisms. The simulation results show that although electrolyte-electron reactions are very efficient, they do not significantly impact photoanode performance within the system as defined. This is because the solution-phase electrolyte chemistry plays a key role in mitigating electron losses through coupled reactions that produce I^- within the photoanode pores, thereby cycling the electrolyte species without requiring that all electrolyte reduction reactions take place at the more distantly located cathode. This is a functionally adaptive response of the chemistry that may be partly responsible for the great success of this redox couple for dye-sensitized solar cells. The simulation results provide predictions that can be tested experimentally.

Diffusion within polymer electrolyte membranes is often coincident with time-dependent processes such as swelling and polymer relaxation, which are factors that limit their ability to block molecular crossover during use. The solution-diffusion model of membrane permeation, which is the accepted theory for dense polymers, applies only to steady-state processes and does not address dynamic internal structural changes that can accompany permeation. To begin discovery of how such changes can be coupled to the permeation process, we have constructed a stochastic multiscale reaction-diffusion model that examines time-dependent methanol uptake into and swelling of hydrated Nafion. Several potential mechanisms of diffusion and polymer response are tested. The simulation predictions are compared to real-time Fourier transform infrared attenuated total reflectance spectroscopy (FTIR-ATR) absorbance reported in the literature [ Hallinan , D. T. , Jr. ; Elabd , Y. A. J. Phys. Chem. B 2007 , 111 , 13221 - 13230 ]. Of the proposed polymer response mechanisms, only one, a reaction-limited, local response to increasing methanol concentration that takes the entire experimental time frame of 600 s, produces simulated FTIR-ATR data consistent with experiment. The simulations show that water diffusion out of the membrane is minimal during methanol sorption and that changes in the measured infrared absorbances are due primarily to the increase in methanol concentration accompanied by dilution of water during swelling. Swelling involves densification of the polymer structure even as there is an overall volume expansion of the film. Potential connections between the polymer densification and molecular-level structural changes of Nafion in methanol are discussed. These results indicate that the interaction between methanol and Nafion serves to increase Nafion's capacity to accommodate large volumes of methanol-water solutions, facilitating increased permeation across the membrane relative to pure water.

Stable hydration in perfluorinated polyelectrolyte membranes such as Nafion is essential to maintain good ion conductivity and manage permeation, especially in vapor-fed devices where water content depends on relative humidity in a gas stream. Extensive studies in the literature have shown that Nafion hydration in water vapor is controlled by its interfacial transport resistance. Nafion forms a fluorine-rich layer at the polymer−gas interface, and it has been proposed that this layer blocks water transport due to its hydrophobicity. To develop a molecular-level description of the physics underlying transport resistance in this system, we have performed a computational reaction−diffusion kinetics study of water evaporation from Nafion. Two distinct models are examined, one mimicking the blocking function proposed in the literature and the other assuming that there is no blocking, treating instead water evaporation as a dynamic balance between uptake from the gas and desorption from the polymer surface. Simulation results are compared to time-dependent infrared data over a range of 100−0% relative humidity from the literature. Only the dynamic model successfully reproduces experimental observations. This indicates that the physical nature of interfacial transport resistance is not slow diffusion across an interfacial layer; rather, it is due to the competition between dehydration and rehydration. The simulation data provide details on the accompanying water distributions throughout the membrane and on interfacial kinetics, showing that they are characterized by strong fluctuations.

Solar harvesting devices using dyes convert the sun's energy to usable forms. The photophysics involved are generally investigated using time-resolved spectroscopic experiments with femtosecond to nanosecond resolution. We show that a kinetic framework constructed from transient and linear absorption measurements of metal-ligand charge transfer states for a set of ruthenium complexes in solution can be used to simulate the steady-state dynamics of dyes adsorbed on a substrate under diffuse solar radiation. Even though the intensity of sunlight is relatively low, double excitations to higher excited states can occur. The steady-state populations show that the dyes' triplet state is the main species present besides the ground state. While small, these persistent excited populations can influence reactivity over the extended periods of time that the systems operate. The results show that non-radiative and optical events (dye^-1 s^-1) within the singlet manifold and from the triplet state exhibit a dependence on ligand substituents.

An error in computing the optical rate coefficients for dye RuP wad discovered in Table S4 of the Supporting Information Section S8 odel S for this article, published in J. Phys. Chem. A, 2021 125 4365-372, 10.1021/acs.jpca.1c02386. The optical rate coefficients include absorption, stimulated emission, and ground-state bleach associated with the ground state and the singlet and triplet excited states as well as for the excited-state absorption between the excited states and implicit higher-energy excited states. The optical rate coefficients for continuous solar irradiance were computed using eq 3 of the main text in which the product of the individual signal components and the solar spectrum are integrated and multiplied by the fundamental optical rate coefficient. We determined the incorrect signal components were used for RuP in the published Supporting Information Table S4. The corrected Table S4 column for RuP is presented with the reported values from the original paper in the corrected Supporting Information. New simulations have been performed using the corrected values, and the results were used to update Figures 5 and 6 and also Figure S9 in the corrected Supporting Information. The published conclusions are unchanged. The populations of RuP excited states, Figure 5, are only marginally smaller, which should be expected with minor changes to the optical rate coefficients. There is not a significant change in the results for the number of radiative and nonradiative interactions per dye per second (the optical transitions and nonradiative transitions increased, respectively, by approximately 2% and 10%, Figure 6A). The occurrence of excited-state absorption (ESA) and excited-state emission (ESE) events, Figure 6C, both decreased from the original simulation results by 20% and 33%. The most notable difference is exhibited in Figure 6B; the fraction of population moving from the lowest-energy singlet state (i.e., intersystem crossing or ultrafast nonradiative relaxation) and nonradiative relaxation from the triplet state have decreased from ∼0.45 to ∼0.4 and increased from 0.35 to 0.4, respectivel.

We examine in a simple organic aerosol the transition between heterogeneous chemistry under well-mixed conditions to chemistry under interfacial confinement. A single reaction mechanism, shown to reproduce observed OH oxidation chemistry for liquid and semisolid C₃₀H₆₂, is used in reaction-diffusion simulations to explore reactivity over a broad viscosity range. The results show that when internal mixing of the aerosol is fast and the particle interface is enriched in C-H groups, ketone and alcohol products, formed via peroxy radical disproportionation, predominate. As viscosity increases the reactions become confined to a shell at the gas-aerosol interface. The confinement is accompanied by emergence of acyloxy reaction pathways that are particularly active when the shell is 1 nm or less. We quantify this trend using a reaction-diffusion index, allowing the parts of the mechanism that control reactivity as viscosity increases to be identified.

Criegee intermediates (CI), formed in alkene ozonolysis, are central for controlling the multiphase chemistry of organic molecules in both indoor and outdoor environments. Here, we examine the heterogeneous ozonolysis of squalene, a key species in indoor air chemistry. Aerosol mass spectrometry is used to investigate how the ozone (O₃) concentration, relative humidity (RH), and particle size control reaction rates and mechanisms. Although the reaction rate is found to be independent of RH, the reaction products and particle size depend upon H₂O. Under dry conditions (RH = 3%) the reaction produces high-molecular-weight secondary ozonides (SOZ), which are known skin irritants, and a modest change in particle size. Increasing the RH reduces the aerosol size by 30%, while producing mainly volatile aldehyde products, increases potential respiratory exposure. Chemical kinetics simulations link the elementary reactions steps of CI to the observed kinetics, product distributions, and changes in particle size. The simulations reveal that ozonolysis occurs near the surface and is O₃-transport limited. The observed secondary ozonides are consistent with the formation of mainly secondary CI, in contrast to gas-phase ozonolysis mechanisms.

To fully utilize the results of laboratory-based studies of the chemistry of model atmospheric aerosol reactions, it is important to understand how to relate them to the conditions found in nature. In this study, we have taken a validated reaction-diffusion mechanism for oxidation of C₃₀H₆₂ aerosol by OH under flow tube conditions and examined its predictions for another experimental regime (continuous flow stirred tank reactor) and for the atmosphere, spanning alkane aerosol viscosities from liquid to semisolid. The results show that under OH-concentration-limited and aerosol-mixing-limited conditions, it should be possible to select laboratory experimental conditions where many aspects of the particle phase and volatile product chemistry under atmospheric conditions can be revealed. If the OH collision and organic diffusion rates are comparable, however, reactivity is highly sensitive to the details of both OH concentration and internal mixing. The characteristics of the transition between limiting conditions provide key insights into which parts of the reaction mechanism dominate in the various kinetic regimes.