Calculated Barrier Height Research Articles

Unlocking a new hydrogen-bonding marker: C-O bond shortening in vicinal diols revealed by rotational spectroscopy.

The conformational space of cis-1,2-cyclohexanediol, a model molecule for cyclic vicinal diols, was investigated using rotational spectroscopy and density functional theory calculations. Four low energy conformers within an energy window of 5 kJ mol-1 were identified computationally. A rotational spectrum of jet-cooled cis-1,2-cyclohexanediol was recorded with a chirped pulse Fourier transform microwave spectrometer. Two sets of rotational transitions were observed and could be assigned to conformers of cis-1,2-cyclohexanediol. The non-observation of other low energy conformers was explained by conformational conversion barrier height calculations and results from experimental spectra recorded with different carrier gases. Eight isotopologues, including those with 13C and 18O, of the lowest energy conformer were observed, allowing the determination of the semi-experimental equilibrium structure, reSE. Interestingly, the structural analysis revealed that the C-O bond length of the intramolecular hydrogen-bond donor is shorter than that of the acceptor. This appears to be a general characteristic of vicinal diols and can be used as a novel hydrogen-bond marker in such compounds.

DLPNO-CCSD(T) and DFT study of the acetate-assisted C-H activation of benzaldimine at [RuCl2(p-cymene)]2: the relevance of ligand exchange processes at ruthenium(II) complexes in polar protic media.

To gain mechanistic insights into the acetate-assisted cyclometallations of arylimines promoted by [RuCl2(p-cymene)]2 in polar protic media, DFT geometry optimizations (with M06 and ωB97X-D3 functionals and the cc-pVDZ-PP[Ru] basis set) followed by DLPNO-CCSD(T)/CBS energy evaluations were performed using benzaldimine as a model substrate and methanol as the solvent (with CPCM or SMD models). The calculation results show that coordination of the imine to an acetate ruthenium precursor is followed by anion (chloride or acetate) dissociation as the rate-determining step of the process. H-Bonding of two explicit MeOH to the anion reduces the calculated activation energy to ca. 23 kcal mol-1, in good agreement with the experimental half-life at room temperature. Subsequent AMLA/CMD C-H activation of the intermediate cationic complexes is a faster, reversible process. Alternative reaction pathways involving neutral diacetate ruthenium complexes offer AMLA/CMD transition state structures of lower energy but are precluded due to higher energy barriers for the initial ligand exchange processes at ruthenium. Solvent assistance accelerates the final chloride/acetate exchange processes on the cycloruthenate intermediates, particularly when compression in the condensed phase is taken into consideration. The performance of six DFT functionals (with the aug-pVTZ-PP[Ru] basis set) was assessed using the DLPNO-CCSD(T)/CBS reference energies. Neutral diacetate ruthenium complexes were incorrectly predicted as being kinetically relevant when using hybrid DFT methods (PBE0-D3(BJ), M06-2X or ωB97M-V). Good agreement between the calculated barrier heights and our benchmark energy results was obtained by using double-hybrid DFT methods. PWPB95 with D3(BJ) or D4 dispersion energy corrections was found to be the most accurate (ΔG≠ MUE of ca. 1 kcal mol-1). This study may aid our understanding of and help with further experimental investigations of synthetically useful carboxylate-assisted C-H bond functionalizations involving (N,C)-cyclometallated (p-cymene)Ru(II) intermediate complexes in sustainable polar protic solvents.

Optical and structural characterizations of NiO-doped PMMA films and fabrication of Ag/NiO-doped PMMA/n-Si/AuSb for analogue-digital switching applications

Herein, the authors report the NiO-doped PMMA nanocomposites fabrication process with different Ni concentrations of 0.001–10 wt% on pre-cleaned glass, ITO, and n-type Si via spin-coating technique. X-ray diffraction (XRD) was used to analyze the crystal structure of NiO in powder form and the NiO-doped PMMA nanocomposite thin films. The Rietveld refinement method was used to analyze the manufactured NiO's powder-form structural characteristics. The analysis data revealed a single phase of NiO nanoparticles with face-centred cubic of Fm-3m space group. The average crystallites' mean size and the microstrain value were determined based on the Williamson-Hall model and were equal to 32.47 ± 1.210 nm and 5 × 10−5, respectively. The FESEM micrographs illustrate dispersed cubic-shaped NiO grains in the PMMA matrix. The estimated energy gaps and the width of the localized energy states (Urbach tails) are plotted against the photon energy. The estimated energy gap, EgOp is increasing from 3.649 ± 0.003 eV to 3.782 ± 0.005 eV with increasing the NiO nanoparticles' doping ratio in the PMMA as a host insulating polymer. The fabricated Ag/NiO-PMMA (S4)/n-Si/AuSb heterojunction device revealed two different (HRS and LRS) resistance-switching behaviours. For second and third electrical (I–V) paths, the extracted magnitudes of the ideality factor are 2.397 ± 0.006 and 2.691 ± 0.003 and the calculated barrier height are 0.727 ± 0.004 and 0.886 ± 0.005, respectively. The interfacial electronics states NSS increasing from 6.58 × 1012 to 3.7 × 1014 cm−2.eV−1 and 3.48 × 1013 to 3.65 × 1014 cm−2 .eV−1 for the 2nd HRS path and 3rd LRS, respectively. The investigated hysteresis loop and the variation in interfacial electronics states reveal that the NiO-PMMA nanocomposite may be promising in non-volatile memory and electrical switching electronic devices.

Bending the rules of PAH hydrogenation: the case of corannulene

ABSTRACT The interaction of a curved polycyclic aromatic hydrocarbon (PAH), corannulene (C20H10), with H-atoms leading to the formation of highly superhydrogenated species was studied. In particular, we determined the H-atom addition sequence to a monolayer of corannulene deposited on a graphite surface in order to understand the effect of curvature on the interaction of PAHs with H-atoms. We also investigated why superhydrogenated corannulene species with a certain number of additional H-atoms are more stable than the others. A combination of thermal desorption mass spectrometry measurements and density functional theory calculations was used to reveal and understand the presence of the most stable configurations in the superhydrogenation sequence of corannulene. Specifically, the experiments demonstrate the presence of stable configurations of superhydrogenated corannulene at specific hydrogenation levels of 2, 6, 8, 10, 12, 14, 16, 17, and 18 extra H-atoms. Density functional theory calculations of binding energies and barrier heights explain why some of these configurations are particularly stable and provide new insights into the superhydrogenation of curved PAH molecules under interstellar conditions.

SBH17: Benchmark Database of Barrier Heights for Dissociative Chemisorption on Transition Metal Surfaces.

Accurate barriers for rate controlling elementary reactions on metal surfaces are key to understanding, controlling, and predicting the rate of heterogeneously catalyzed processes. While barrier heights for gas phase reactions have been extensively benchmarked, dissociative chemisorption barriers for the reactions of molecules on metal surfaces have received much less attention. The first database called SBH10 and containing 10 entries was recently constructed based on the specific reaction parameter approach to density functional theory (SRP-DFT) and experimental results. We have now constructed a new and improved database (SBH17) containing 17 entries based on SRP-DFT and experiments. For this new SBH17 benchmark study, we have tested three algorithms (high, medium, and light) for calculating barrier heights for dissociative chemisorption on metals, which we have named for the amount of computational effort involved in their use. We test the performance of 14 density functionals at the GGA, GGA+vdW-DF, and meta-GGA rungs. Our results show that, in contrast with the previous SBH10 study where the BEEF-vdW-DF2 functional seemed to be most accurate, the workhorse functional PBE and the MS2 density functional are the most accurate of the GGA and meta-GGA functionals tested. Of the GGA+vdW functionals tested, the SRP32-vdW-DF1 functional is the most accurate. Additionally, we found that the medium algorithm is accurate enough for assessing the performance of the density functionals tested, while it avoids geometry optimizations of minimum barrier geometries for each density functional tested. The medium algorithm does require metal lattice constants and interlayer distances that are optimized separately for each functional. While these are avoided in the light algorithm, this algorithm is found not to give a reliable description of functional performance. The combination of relative ease of use and demonstrated reliability of the medium algorithm will likely pave the way for incorporation of the SBH17 database in larger databases used for testing new density functionals and electronic structure methods.

Description of Si and Al Release from Aluminosilicate in the Acidic Condition Using Density Functional Theory: Protonated Terminal Oxygen

The molecular clusters ((HO)3Si-O-Si(OH)3 and (HO)3Al-O-Si(OH)3) representative of aluminosilicate mineral surface were employed to study the dissolution of aluminosilicate in acidic condition via density functional theory (DFT) with the M06-2X+G(d,p) methodology. The surface termination sites (Si and Al) were both tetra-coordinated and the terminal oxygen was protonated in an acidic condition. In the dissolution reaction, the calculated barrier height of the six-membered ring transition state complex containing two water molecules was predicted to be 76.13 kJ/mol, lower than that of the four-membered ring transition state complex containing one water molecule. The barrier height of the reaction decreased to 6.17 kJ/mol and was 91.90% lower than that for the Siter-O-Si without protonation. In addition, the calculated barrier heights for Al-terminated sites were predicted to be 22.23 kJ/mol, lower than those for the Si-terminated sites, suggesting that breaking the Al-O bond is easier than the Si-O bond in the aluminosilicate mineral surface. With the fracture of Si-O and Al-O bonds, the Si and Al release from the aluminosilicate. The results indicate that the acidic condition facilitates the release of Si and Al from the aluminosilicate, and the concentration of Al leaching from the aluminosilicate is higher than the Si.

Best-Practice DFT Protocols for Basic Molecular Computational Chemistry.

Nowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum-chemical calculations applies density functional theory (DFT) evaluated in atomic-orbital basis sets. This work provides best-practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step-by-step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi-level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.

Best‐Practice DFT Protocols for Basic Molecular Computational Chemistry**

AbstractNowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum‐chemical calculations applies density functional theory (DFT) evaluated in atomic‐orbital basis sets. This work provides best‐practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step‐by‐step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi‐level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.

Simple Approximation for the Ideal Reference State of Gases Adsorbed on Solid-State Surfaces

Reference states are useful as models for facilitating calculations of equilibrium constants, and they may also serve as standard states that are convenient for organizing and tabulating thermodynamic data; however, standard state conventions and appropriate reference states for adsorbed species have received less attention than those for pure substances and solutes. Here, we compare seven choices of reference states for calculations of equilibrium constants and transition state theory rate constants for flat surfaces, in particular (1) an ideal 2D harmonic oscillator, (2) an ideal rigid-molecule harmonic oscillator, (3) an ideal 2D harmonic oscillator with separable surface modes, (4) a 2D ideal gas, (5) an ideal 2D hindered translator, (6) an ideal 2D hindered translator with lowest-order barriers, and (7) a simple ideal 2D hindered translator proposed in this work. The advantage of models 5-7 is that they can treat both mobile and localized adsorbates in a consistent way, whereas models 1-3 are only appropriate for localized adsorbates, and model 4 is only appropriate for a freely translating adsorbate. Furthermore, models 6 and 7 reduce the computational cost without the user having to calculate barrier heights for diffusion. An advantage of the simple ideal 2D hindered translator is that it has a physical high-temperature limit. We also propose a reference state for nonflat surfaces. The user is encouraged to choose a reference state based on the appropriateness of the model and the practicality of the calculations.

Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials.

Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.

Solvent Effects on the Menshutkin Reaction.

The Menshutkin reaction is a methyl transfer reaction relevant in fields ranging from biochemistry to chemical synthesis. In the present work, the energetics and solvent distributions for NH3+MeCl and Pyr+MeBr reactions were investigated in explicit solvent (water, methanol, acetonitrile, benzene, cyclohexane) by means of reactive molecular dynamics simulations. For polar solvents (water, methanol, and acetonitrile) and benzene, strong to moderate catalytic effects for both reactions were found, whereas apolar and bulky cyclohexane interacts weakly with the solute and does not show pronounced barrier reduction. The calculated barrier heights for the Pyr+MeBr reaction in acetonitrile and cyclohexane are 23.2 and 28.1 kcal/mol compared with experimentally measured barriers of 22.5 and 27.6 kcal/mol, respectively. The solvent distributions change considerably between reactant and TS but comparatively little between TS and product conformations of the solute. As the system approaches the transition state, correlated solvent motions occur which destabilize the solvent-solvent interactions. This is required for the system to surmount the barrier. Finally, it is found that the average solvent-solvent interaction energies in the reactant, TS, and product state geometries are correlated with changes in the solvent structure around the solute.

DFT Study on Side Chain Detachment of Perfluorosulfonic Acid Ionomers by Radical-Assisted Nucleophilic Attack of Water

Chemical degradation of perfluorosulfonic acid (PFSA) ionomer membranes is still an unresolved critical issue in the development of polymer electrolyte fuel cells (PEFCs). Its reaction mechanism has not been settled yet even after experimental and computational studies for decades. In this study, the role of hydroxyl radical to side chain degradation is revisited and a new reaction mechanism is proposed. Its nucleophilic attack on the carbon atom in the ether group induces detachment of the side chain. The reactions are supposed to be feasible as their reaction barrier heights are only slightly higher than those of tertiary fluorine atom abstraction reactions by hydrogen radical in absence of the solvation effect. When a single water molecule is present with hydroxyl or hydroperoxyl radical it is the oxygen in the water that undergoes the nucleophilic attack on the carbon in the ether groups and the calculated barrier heights are lowered to 70∼100 kJ mol−1 as the radical assists the reaction as a hydrogen acceptor. Possible reaction pathways of the products are also calculated to explain the occurrence of the experimentally observed degradation products of Nafion after chemical durability tests.

Nuclear surface parameters revisited

The atomic nucleus is a leptodermous system in which the related densities fade gradually through the surface region. The surface is usually described by two parameters, one to describe the radius and the second to describe the surface width. In this study, we present a clear definition of the nucleons distribution surface based on the investigation of the surface region independent of the nuclear core. We also present an improved method for calculating the surface width relying only on the surface nucleons, avoiding the effect of the core nucleons. The proposed surface parameters are tested for centrally depressed, and centrally raised, distributions and the results confirmed the improvement in the surface description and the elimination of the exaggerated effect of core density on the values of surface parameters. The surface radii and surface widths, based on the proposed method, are presented for experimental proton distributions of 64 nuclei. Including the surface width by the proposed method in the nuclear proximity model improved the calculation of fusion barrier height and position compared to the experimental results.

BH9, a New Comprehensive Benchmark Data Set for Barrier Heights and Reaction Energies: Assessment of Density Functional Approximations and Basis Set Incompleteness Potentials.

The calculation of accurate reaction energies and barrier heights is essential in computational studies of reaction mechanisms and thermochemistry. To assess methods regarding their ability to predict these two properties, high-quality benchmark sets are required that comprise a reasonably large and diverse set of organic reactions. Due to the time-consuming nature of both locating transition states and computing accurate reference energies for reactions involving large molecules, previous benchmark sets have been limited in scope, the number of reactions considered, and the size of the reactant and product molecules. Recent advances in coupled-cluster theory, in particular local correlation methods like DLPNO-CCSD(T), now allow the calculation of reaction energies and barrier heights for relatively large systems. In this work, we present a comprehensive and diverse benchmark set of barrier heights and reaction energies based on DLPNO-CCSD(T)/CBS called BH9. BH9 comprises 449 chemical reactions belonging to nine types common in organic chemistry and biochemistry. We examine the accuracy of DLPNO-CCSD(T) vis-a-vis canonical CCSD(T) for a subset of BH9 and conclude that, although there is a penalty in using the DLPNO approximation, the reference data are accurate enough to serve as a benchmark for density functional theory (DFT) methods. We then present two applications of the BH9 set. First, we examine the performance of several density functional approximations commonly used in thermochemical and mechanistic studies. Second, we assess our basis set incompleteness potentials regarding their ability to mitigate basis set incompleteness errors. The number of data points, the diversity of the reactions considered, and the relatively large size of the reactant molecules make BH9 the most comprehensive thermochemical benchmark set to date and a useful tool for the development and assessment of computational methods.

Ab initio -based approach for the oxidation mechanisms at SiO2 /4H-SiC interface: Interplay of dry and wet oxidants during interfacial reaction

Effects of the coexistence of dry and wet oxidants on the reaction process at ${\mathrm{SiO}}_{2}$/4H-SiC(0001) and ($000\overline{1}$) interfaces during SiC oxidation are systematically investigated by performing ab initio calculations. We find characteristic features of the interfacial reaction mechanisms, which are dependent on the plane orientation and wet oxidation condition. The incorporation of wet oxidants leads to lower barrier heights for the reaction of ${\mathrm{O}}_{2}$ molecules, resulting in the promotion of interfacial reaction at the ${\mathrm{SiO}}_{2}$/4H-SiC(0001) interface. In contrast, the coexisting ${\mathrm{O}}_{2}$ molecule assists the reactions by ${\mathrm{H}}_{2}\mathrm{O}$ molecule and OH groups at the ${\mathrm{SiO}}_{2}$/4H-SiC($000\overline{1}$) interface. Furthermore, we estimate the linear rate constants in the Deal-Grove model using the calculated barrier heights and reveal that the preferable wet ambient condition at the ${\mathrm{SiO}}_{2}$/4H-SiC(0001) interface is different from that at the ${\mathrm{SiO}}_{2}$/4H-SiC($000\overline{1}$) interface. The difference in the rate constants is caused by the difference in the rate-limiting reaction pathway between ${\mathrm{SiO}}_{2}$/4H-SiC(0001) and ($000\overline{1}$) interfaces. These calculated results offer better understanding of the atom-scale mechanisms for the significant enhancement of SiC oxidation by the interplay of dry and wet oxidants.

Structural Dynamics and Catalytic Mechanism of ATP13A2 (PARK9) from Simulations.

ATP13A2 is a gene encoding a protein of the P5B subfamily of ATPases and is a PARK gene. Molecular defects of the gene are mainly associated with variations of Parkinson's disease (PD). Despite the established importance of the protein in regulating neuronal integrity, the three-dimensional structure of the protein currently remains unresolved crystallographically. We have modeled the structure and reactivity of the full-length protein in its E1-ATP state. Using molecular dynamics (MD), quantum cluster, and quantum mechanical/molecular mechanical (QM/MM) methods, we aimed at describing the main catalytic reaction, leading to the phosphorylation of Asp513. Our MD simulations suggest that two positively charged Mg2+ cations are present at the active site during the catalytic reaction, stabilizing a specific triphosphate binding mode. Using QM/MM calculations, we subsequently calculated the reaction profiles for the phosphoryl transfer step in the presence of one and two Mg2+ cations. The calculated barrier heights in both cases are found to be ∼12.5 and 7.5 kcal mol-1, respectively. We elucidated details of the catalytically competent ATP conformation and the binding mode of the second Mg2+ cofactor. We also examined the role of the conserved Arg686 and Lys859 catalytic residues. We observed that by significantly lowering the barrier height of the ATP cleavage reaction, Arg686 had major effect on the reaction. The removal of Arg686 increased the barrier height for the ATP cleavage by more than 5.0 kcal mol-1 while the removal of key electrostatic interactions created by Lys859 to the γ-phosphate and Asp513 destabilizes the reactant state. When missense mutations occur in close proximity to an active site residue, they can interfere with the barrier height of the reaction, which can halt the normal enzymatic rate of the protein. We also found large binding pockets in the full-length structure, including a transmembrane domain pocket, which is likely where the ATP13A2 cargo binds.

Theoretical investigation on the addition reaction of the germylenoid H2GeLiCl with acetone.

In this work, theoretical calculations were performed on the addition reaction of the germylenoid H2GeLiCl with acetone. The DFT M06-2X method was used to optimize the geometries of the whole stationary points on the potential energy surfaces and the QCISD method to calculate the single-point energy. The results reveal that the addition reaction of H2GeLiCl with acetone firstly generates an oxagermacyclopropane c-H2GeOC(CH3)2 and then c-H2GeOC(CH3)2 further reacts with acetone along two possible pathways, pathway I and pathway II, in which the 2,4-dioxagermolane is formed at the end of pathway I and 2,5-dioxagermolane is formed at the end of pathway II, respectively. According to the calculated barrier heights, we can deduce that the pathway I is more favorable than pathway II. The computational results suggest that this reaction model can provide new inspiration for the synthesis of heterocyclic germanium compounds.

Synthesis and thermal stability, electrical conductivity and dielectric spectroscopic studies of poly (ethylene-co-vinyl alcohol)/graphene oxide nanocomposite

A new nanocomposite was successfully synthesized by incorporating the prepared graphene oxide into poly(ethylene-co-vinyl alcohol) (EVOH) copolymer via simple solution casting method. XRD, FT-IR, SEM and TGA analysis were used to characterize the current nanocomposite. AC electrical transport properties of the nanocomposite were studied in the frequency range (10 Hz-1MHz) over the temperature range (298–423 K). The AC conductivity results were interpreted using jump relaxation model and were found to obey the universal double power law. The conduction mechanism inside the current samples has changed from non-overlapping small polaron tunneling to correlated barrier hopping. The calculated barrier height has indicated that the predominant charge carriers in the EVOH/GO nanocomposite are the single polarons and not bipolarons. The enhancement in σac of the EVOH-based nanocomposite, making these samples highly recommended for use in energy storage applications such as; electronic capacitors where materials with high dielectric constant and low dielectric loss are required.

Memristiv effect on W/Ti/p-Si structure: Aging phenomena and one of the origin of barrier inhomogeneity

Memristor as claimed to be fourth two terminal passive devices in electronics. Experimental realization memristor by HP research group is getting much more attention due to its potential applications. However, in reality this device with its unique properties always exist in MOS, MIS or even MS junctions. But its visibility has a dependence to oxide or interface thickness with mobile ions. As reported in this study, current–voltage (I-V) relation of W/Ti/p-Si continuously modified by the forward bias current. Departure form linear relation between ln(I)-V commonly analyzed with barrier inhomogeneity assumption. The appearance of barrier inhomogeneity for MS junction reported to be one of the reasons for non-ideal ln(I)-V relation. Investigation in this study shows that total charge Q(t) controlled by voltage sweep rate effecting forward ln(I)-V relation. If this situation exists, calculated barrier heights and ideality factors are not true and a pure barrier inhomogeneity approximation may not be applicable. However, under the reverse bias condition, the total charge during entire measurement is so small which it’s effect negligible over I-V relation as observed. The one of the solution is the use Missous-Rhoderick method for finding Schottky diode parameters or creating current stress over MS junction which saturates mobile ion drift. Voltage sweep rate during I-V measurement and hence total charge, affecting current transport and in parallel appearance of hysteresis loop indicate interface layer which behave as a memristor.

Thermal desorption and pyrolysis direct analysis in real time mass spectrometry of Nafion membrane

AbstractPerfluorosulfonic acid ionomer membranes have been widely used as proton conducting membranes in various electrochemical processes such as polymer electrolyte fuel cells and water electrolysis. While their thermal stability has been studied by thermogravimetry and analysis of low molecular weight products, their decomposition mechanism is little understood. In this study a newly developed methodology of thermal desorption and pyrolysis in combination with direct analysis in real time mass spectrometry is applied for Nafion membrane. An ambient ionization source and a high‐resolution time‐of‐flight mass spectrometer enabled unambiguous assignment of gaseous products. Thermal decomposition is initiated by side chain detachment above 350°C, which leaves carbonyls on the main chain at the locations of the side chains. Perfluoroalkanes are released above 400°C by main chain scission and their further decomposition products dominate above 500 °C. DFT calculation of reaction energies and barrier heights of model compounds support proposed decomposition reactions.