Incorporation of Charge Transfer into Classical Molecular Dynamics Force Fields with Applications in Physical Chemistry.

The presence of charge transfer (CT) interactions is clear in a variety of systems. In CT, some electron density is shifted from one molecule to another (non-bonded) molecule. The importance of this CT interaction is unclear. Previous attempts to look at the conse- quences of CT required the use of ab initio molecular dynamics (AIMD), a computationally intensive method. Herein, a method for including CT in force field (FF) simulations is described. It is efficient, produces charges in agreement with AIMD, and prevents long- ranged CT. This CT MD method has been applied to monatomic ions in water. When solvated, ions do not have an integer charge. Anions give up some electron density to their ligands, and cations receive some electron density from their ligands. In bulk, the first solvation shell does not compensate for all CT, i.e. the charge is not smeared out over the first solvation shell. Rather, some charge is also found in the second solvation shell and further into the bulk. The charge of the first solvation shell depends on the balance between ion-water and water-water CT. When an interface is present, the charge outside of the second solvation shell will reside at the interface. This occurs even when the ion is over 15 Å away from the surface. The effect of long-ranged CT is mediated by changes in the hydrogen bonding patterns in water induced by the ions (not direct CT from the ions to distant waters). The model has also been applied to water’s ‘‘self-ions’’ hydronium and hydroxide. Trajectories from the multi-state empirical valence bond model (MS-EVB3) are analyzed. The differences between monatomic and molecular ions are explored. The direction of CT and the effect of hydrogen bonding with the ion are considered. The damping of CT as ligands are added is discussed and a method to improve the MD model, in order to account for damping, is proposed

Permeation of CO₂ and N₂ through glassy poly(dimethyl phenylene) oxide under steady- and presteady-state conditions

Glassy polymers are often used for gas separations because of their high selectivity. Although the dual‐mode permeation model correctly fits their sorption and permeation isotherms, its physical interpretation is disputed, and it does not describe permeation far from steady state, a condition expected when separations involve intermittent renewable energy sources. To develop a more comprehensive permeation model, we combine experiment, molecular dynamics, and multiscale reaction–diffusion modeling to characterize the time‐dependent permeation of N₂ and CO₂ through a glassy poly(dimethyl phenylene oxide) membrane, a model system. Simulations of experimental time‐dependent permeation data for both gases in the presteady‐state and steady‐state regimes show that both single‐ and dual‐mode reaction–diffusion models reproduce the experimental observations, and that sorbed gas concentrations lag the external pressure rise. The results point to environment‐sensitive diffusion coefficients as a vital characteristic of transport in glassy polymers

Predictive simulation of non-steady-state transport of gases through rubbery polymer membranes

A multiscale, physically-based, reaction-diffusion kinetics model is developed for non-steady-state transport of simple gases through a rubbery polymer. Experimental data from the literature, new measurements of non-steady-state permeation and a molecular dynamics simulation of a gas-polymer sticking probability for a typical system are used to construct and validate the model framework. Using no adjustable parameters, the model successfully reproduces time-dependent experimental data for two distinct systems: (1) O_2 quenching of a phosphorescent dye embedded in poly(n-butyl(amino) thionylphosphazene), and (2) O_2, N_2, CH_4 and CO_2 transport through poly(dimethyl siloxane). The calculations show that in the pre-steady-state regime, permeation is only correctly described if the sorbed gas concentration in the polymer is dynamically determined by the rise in pressure. The framework is used to predict selectivity targets for two applications involving rubbery membranes: CO_2 capture from air and blocking of methane cross-over in an aged solar fuels device

Predictive simulation of non-steady-state transport of gases through rubbery polymer membranes

A multiscale, physically-based, reaction-diffusion kinetics model is developed for non-steady-state transport of simple gases through a rubbery polymer. Experimental data from the literature, new measurements of non-steady-state permeation and a molecular dynamics simulation of a gas-polymer sticking probability for a typical system are used to construct and validate the model framework. Using no adjustable parameters, the model successfully reproduces time-dependent experimental data for two distinct systems: (1) O_2 quenching of a phosphorescent dye embedded in poly(n-butyl(amino) thionylphosphazene), and (2) O_2, N_2, CH_4 and CO_2 transport through poly(dimethyl siloxane). The calculations show that in the pre-steady-state regime, permeation is only correctly described if the sorbed gas concentration in the polymer is dynamically determined by the rise in pressure. The framework is used to predict selectivity targets for two applications involving rubbery membranes: CO_2 capture from air and blocking of methane cross-over in an aged solar fuels device

Incorporation of Charge Transfer into Classical Molecular Dynamics Force Fields with Applications in Physical Chemistry.

The presence of charge transfer (CT) interactions is clear in a variety of systems. In CT, some electron density is shifted from one molecule to another (non-bonded) molecule. The importance of this CT interaction is unclear. Previous attempts to look at the conse- quences of CT required the use of ab initio molecular dynamics (AIMD), a computationally intensive method. Herein, a method for including CT in force field (FF) simulations is described. It is efficient, produces charges in agreement with AIMD, and prevents long- ranged CT. This CT MD method has been applied to monatomic ions in water. When solvated, ions do not have an integer charge. Anions give up some electron density to their ligands, and cations receive some electron density from their ligands. In bulk, the first solvation shell does not compensate for all CT, i.e. the charge is not smeared out over the first solvation shell. Rather, some charge is also found in the second solvation shell and further into the bulk. The charge of the first solvation shell depends on the balance between ion-water and water-water CT. When an interface is present, the charge outside of the second solvation shell will reside at the interface. This occurs even when the ion is over 15 Å away from the surface. The effect of long-ranged CT is mediated by changes in the hydrogen bonding patterns in water induced by the ions (not direct CT from the ions to distant waters). The model has also been applied to water’s ‘‘self-ions’’ hydronium and hydroxide. Trajectories from the multi-state empirical valence bond model (MS-EVB3) are analyzed. The differences between monatomic and molecular ions are explored. The direction of CT and the effect of hydrogen bonding with the ion are considered. The damping of CT as ligands are added is discussed and a method to improve the MD model, in order to account for damping, is proposed

Swelling and Diffusion during Methanol Sorption into Hydrated Nafion.

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

Swelling and Diffusion during Methanol Sorption into Hydrated Nafion.

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