An explanation for the superior proton conductivity of low equivalent weight (EW) short-side-chain (SSC) perfluorosulfonic acid membranes is pursued through the determination of hydrated morphology and hydronium ion diffusion coefficients using classical molecular dynamics (MD) simulations. A unique force field set for the SSC ionomer was derived from torsion profiles determined from ab initio electronic structure calculations of an oligomeric fragment consisting of two side chains. MD simulations were performed on a system consisting of a single macromolecule of the polymer (EW of 580) with the general formula F3C-[CF(OCF2CF2SO3H)-(CF2)7]40-CF3 at hydration levels corresponding to 3, 6, and 13 water molecules per sulfonic acid group. Examination of the hydrated morphology indicates the formation of hydrogen bond "bridges" between distant sulfonate groups without significant bending of the polytetrafluoroethylene backbone. Pair correlation functions of the system identify the presence of ion cages consisting of hydronium ions hydrogen-bonded to three sulfonate groups at the lowest water content. Such structures exhibit very low S-OH3+ separations, well below 4 A and severely inhibit vehicular diffusion of the protons. The number of sulfonate groups in the first solvation shell of a given hydronium ion correlates well with the differences between Nafion and the SSC polymer (Hyflon). The calculated hydronium ion diffusion coefficients of 2.84 x 10-7, 1.36 x 10-6, and 3.47 x 10-6 cm2/s for water contents of 3, 6, and 13, respectively, show only good agreement to experimentally measured values at the lowest water content, underscoring the increasing contribution of proton shuttling or hopping at the higher hydration levels. At the highest water content, the vehicular diffusion accounts for only about 1/5 of the total proton transport similar to that observed in Nafion.
Density functional theory (DFT) together with Car-Parrinello ab initio molecular dynamics (CP-AIMD) simulation has been used to investigate the free energy profiles of two representative S N 2 reactions: (A) Cl -+ CH 3 Cl f ClCH 3 + Cl -; (B) NH 3 + H 3 BNH 3 f H 3 NBH 3 + NH 3 . The free energy profiles along the reaction coordinates at 300 K and 600 K were determined directly by a pointwise thermodynamic integration (PTI) technique. Comparison between the well-known double-well potential energy profile (PEP) and the free energy profiles (FEP) has been made. The results show that, for reaction A, the double-well profile is maintained for the FEP at 300 K due to the stronger ion-dipole interaction between chloromethane and the chloride anion. In comparison with the PEP, the FEP has a higher central barrier and a more shallow well depth. However, at 600 K the double wells almost disappear on the FEP, whereas the central barrier increases further. For reaction B, the 300 K FEP also presents a higher central barrier peak and a more shallow well depth compared to the PEP. However, when the temperature increases to 600 K, a saddle shape FEP is obtained, which indicates that the reaction has changed mechanism from an associative S N 2 reaction to a dissociative S N 1 reaction. This change is driven by entropy.
The Catalytica process converts methane to methyl bisulfate in good yield at relatively low temperature in fuming sulfuric acid with (bipyrimidine)PtCl 2 as a catalyst. Previously we examined the first step, methane C-H activation, and here we look at the oxidation by SO 3 and the reductive elimination steps. In the oxidation step a Pt(II)-CH 3 complex (a) reacts with protonated SO 3 , which splits to form two new ligands, SO 2 and OH -, thus oxidizing a to a Pt(IV)-CH 3 complex. The final step in the cycle is the reductive elimination of methyl bisulfate from this complex.
Presented here is the application of a scheme for optimizing the structures of minima and transition states on the free energy surface (FES) for a path along a fixed reaction coordinate with the aid of ab initio molecular dynamics (AIMD) simulation. In the direction of the reaction coordinate, the values corresponding to the stationary points were optimized using the quasi-Newton method, in which the gradient of the free energy along the reaction coordinate was obtained by a constraint AIMD method, and the Bofill Hessian update scheme was used. The equilibrium values for the other directions were taken as the corresponding averages in the dynamic simulation. This scheme was applied to several elementary bimolecular addition reactions: (A) BH(3) + H(2)O --> H(2)O.BH(3); (B) BF(3) + NH(3) --> FB(3).NH(3); (C) SO(3) + NH(3) --> O(3)S.NH(3); (D) C(2)H(4) + CCl(2) --> H(4)C(2).CCl(2); (E) Ni(NH(2))(2) + PH(3) --> (NH(2))(2)Ni.PH(3); (F) W(CO)(5) + CO --> W(CO)(6). For reactions A, B, C, and F, no transition state (TS) exists on the potential energy surface (PES). However there is a TS on the FES. This stems from the curvature difference of the PES and -TDeltaS as a function of the reaction coordinate. For all reactions, it is found that the TS shifts toward the complexation product with increasing temperature because of the curvature increase of -TDeltaS. The equilibrium bond distances for the inactive coordinates perpendicular to the reaction coordinate always increase with temperature, which is due to the thermal excitation and anharmonicity of the PES.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.