Elucidating the Effects of Chemisorbed Water Molecules on the Adhesive Interactions of Epoxy Resin to γ-Alumina Surfaces.

The adhesion mechanism of epoxy resin to the γ-alumina (110) surface was investigated using first-principles density functional theory (DFT). Aluminum materials are lightweight and are used in a wide range of industrial fields. Its surface is oxidized to alumina, and the stable surface is known as the γ-alumina (110) surface. The coverage of hydroxy groups by chemisorbed water molecules on this surface varied depending on the pretreatment temperature. In this study, we investigated the adhesive interactions of epoxy resin on four alumina surfaces with different densities of surface hydroxy groups (0, 3, 6, and 9 OH/nm2) and have discussed their effects. At each interface, the energy curves of the vertically displaced epoxy resin were calculated and the adhesive forces were estimated by differentiating these curves. As the coverage of the surface hydroxy groups increased from 0 to 6 OH/nm2, the adhesive strength gradually decreased. However, the adhesive strength at 9 OH/nm2 was relatively large and almost equal to that at 3 OH/nm2. This inverse volcano-type behavior was analyzed via the decomposition of adhesive forces and the crystal orbital Hamilton population (COHP). The decomposition of adhesive forces into DFT and dispersion components revealed that the inverse volcano-type behavior is derived from the DFT component, and the interfacial interactions owing to the DFT component are accompanied by charge transfer. These were investigated using a COHP analysis, which revealed that this behavior was caused by changes in the activity of the aluminum atoms on the surface and surface reconstruction by chemisorbed water molecules. It is noteworthy that the adhesive strength for 9 OH/nm2 was only 6.9% lower than that for 0 OH/nm2 wherein the chemisorbed water molecules were completely removed from the surface. These results are expected to provide a guideline for the adhesion of epoxy resin to aluminum materials.

Molecular Understanding of Adhesion of Epoxy Resin to Graphene and Graphene Oxide Surfaces in Terms of Orbital Interactions.

The adhesion mechanism of epoxy resin (ER) cured material consisting of diglycidyl ether of bisphenol A (DGEBA) and 4,4'-diaminodiphenyl sulfone (DDS) to pristine graphene and graphene oxide (GO) surfaces is investigated on the basis of first-principles density functional theory (DFT) with dispersion correction. Graphene is often used as a reinforcing filler incorporated into ER polymer matrices. The adhesion strength is significantly improved by using GO obtained by the oxidation of graphene. The interfacial interactions at the ER/graphene and ER/GO interfaces were analyzed to clarify the origin of this adhesion. The contribution of dispersion interaction to the adhesive stress at the two interfaces is almost identical. In contrast, the DFT energy contribution is found to be more significant at the ER/GO interface. Crystal orbital Hamiltonian population (COHP) analysis suggests the existence of hydrogen bonding (H-bonding) between the hydroxyl, epoxide, amine, and sulfonyl groups of the ER cured with DDS and the hydroxyl groups of the GO surface, in addition to the OH-π interaction between the benzene rings of ER and the hydroxyl groups of the GO surface. The H-bond has a large orbital interaction energy, which is found to contribute significantly to the adhesive strength at the ER/GO interface. The overall interaction at the ER/graphene is much weaker due to antibonding type interactions just below the Fermi level. This finding indicates that only dispersion interaction is significant when ER is adsorbed on the graphene surface.

Shear adhesive strength between epoxy resin and copper surfaces: a density functional theory study.

Adhesive strength varies greatly with direction; various adhesion tests have been conducted. In this study, the shear and tensile adhesive strength of epoxy resin for copper (Cu) and copper oxide (Cu2O) surfaces were estimated based on quantum chemical calculations. Here, density functional theory (DFT) calculations with dispersion correction were used. In the tensile process, the entire epoxy resin is peeled off vertically, whereas in the shear process, a force parallel to the adhesive surface is applied. Then, a bending moment acts on the adhesive layer, and a total force (stress) inclined at an angle θ with respect to the adherend surface is applied to the adhesive interface. We computed adhesive stress-displacement curves for each θ exhaustively and discussed the changes. When θ equals 90°, it corresponds to a tensile process. As θ decreases from 90°, the shear adhesive stress on both surfaces decreases slowly. When θ is less than 30°, the constraint to the surface causes periodic changes in the adhesive stress curves. The constraint to the Cu2O surface is especially strong, and this change is large. This periodicity is similar to the stick-slip phenomenon in tribology. To further understand the shear adhesive forces, force decomposition analysis was performed, revealing that the periodicity of the adhesive stress originates from the DFT contribution rather than the dispersion one. The procedure proposed in this study for estimating shear adhesive strength is expected to be useful in the evaluation and prediction of adhesive and adherend properties.

Peel Adhesion Strength between Epoxy Resin and Hydrated Silica Surfaces: A Density Functional Theory Study.

Adhesive strength is known to change significantly depending on the direction of the force applied. In this study, the peel and tensile adhesive forces between the hydroxylated silica (001) surface and epoxy resin are estimated based on quantum chemical calculations. Here, density functional theory (DFT) with dispersion correction is used. In the peel process, the epoxy resin is pulled off from the terminal part, while in the tensile process, the entire epoxy resin is pulled off vertically. As a result of these calculations, the maximum adhesive force in the peel process is decreased to be about 40% of that in the tensile process. The adhesion force-displacement curve for the peeling process shows two characteristic peaks corresponding to the process where the adhesive molecule horizontally oriented to the surface shifts to a vertical orientation to the surface and the process where the vertical adhesive molecule is dissociated from the surface. Force decomposition analysis is performed to further understand the peel adhesion force; the contribution of the dispersion force is found to be slightly larger than that of the DFT force. This feature is common to the tensile process as well. Each force in the peel process is about 40% smaller than the corresponding force in the tensile process.

Quantum Chemical Calculations to Trace Back Reaction Paths for the Prediction of Reactants.

The long-due development of a computational method for the ab initio prediction of chemical reactants that provide a target compound has been hampered by the combinatorial explosion that occurs when reactions consist of multiple elementary reaction processes. To address this challenge, we have developed a quantum chemical calculation method that can enumerate the reactant candidates from a given target compound by combining an exhaustive automated reaction path search method with a kinetics method for narrowing down the possibilities. Two conventional name reactions were then assessed by tracing back the reaction paths using this new method to determine whether the known reactants could be identified. Our method is expected to be a powerful tool for the prediction of reactants and the discovery of new reactions.

Palladium-Catalyzed C-H Iodination of Arenes by Means of Sulfinyl Directing Groups.

C-H iodination of aromatic compounds has been accomplished with the aid of sulfinyl directing groups under palladium catalysis. The reaction proceeds selectively at the peri-position of polycyclic aryl sulfoxides or at the ortho-position of phenyl sulfoxides. The iodination products can be further converted via iterative catalytic cross-coupling at the expense of the C-I and C-S bonds. Computational studies suggest that peri-C-H palladation would proceed via a non-directed pathway, wherein neither of the sulfur nor oxygen atom of the sulfinyl group coordinates to the palladium before and at the transition state.

Understanding CO oxidation on the Pt(111) surface based on a reaction route network.

Analysis of a reaction on a solid surface is an important task for understanding the catalytic reaction mechanism. In this study, we studied CO oxidation on the Pt(111) surface by using the artificial force induced reaction (AFIR) method. A systematic reaction path search was done, and the reaction route network was created. This network included not only bond rearrangement paths but also migration paths of adsorbed species. Then, the obtained network was analyzed using a kinetics method called rate constant matrix contraction (RCMC). It is found that the bottleneck of the overall reaction is the CO2 generation step from an adsorbed CO molecule and an O atom. This result is consistent with the Langmuir-Hinshelwood (LH) mechanism with O2 dissociation discussed in previous studies. The present procedure, i.e., construction of the reaction route network by the AFIR method followed by application of the RCMC kinetics method to the resultant reaction route network, was fully systematic and uncovered two aspects: the impact of the existence of multiple paths in each bond rearrangement step and an entropic contribution arising from short-range migration of adsorbed species.

Designing the Backbone of Hexasilabenzene Derivatives with a High Unimolecular Kinetic Stability.

It is an important subject to theoretically predict the kinetic stability of transient species. In this study, we have studied the kinetic stability of hexasilabenzene Si6 H6 and its derivatives, that is, decasilanaphthalene Si10 H8 and Li-substituted hexasilabenzene Si6 Li6 , theoretically by the artificial force induced reaction (AFIR) method combined with the rate constant matrix contraction (RCMC) method. Molecular design was further conducted to extend the unimolecular lifetime of hexasilabenzene derivatives. Although both Si10 H8 and Si6 Li6 were shown to possess shorter lifetimes than Si6 H6 , we found that the lifetimes of Si6 Li6 changed depending on arrangements of Li atoms around the monocyclic Si6 backbone. Based on this knowledge, we found that a compound of an atomic composition Si6 H4 Li2 with a planar, monocyclic Si6 backbone has a relatively long unimolecular lifetime. Moreover, substitution of the two Li atoms by Na atoms further increased the lifetime.

Implementation and performance of the artificial force induced reaction method in the GRRM17 program.

This article reports implementation and performance of the artificial force induced reaction (AFIR) method in the upcoming 2017 version of GRRM program (GRRM17). The AFIR method, which is one of automated reaction path search methods, induces geometrical deformations in a system by pushing or pulling fragments defined in the system by an artificial force. In GRRM17, three different algorithms, that is, multicomponent algorithm (MC-AFIR), single-component algorithm (SC-AFIR), and double-sphere algorithm (DS-AFIR), are available, where the MC-AFIR was the only algorithm which has been available in the previous 2014 version. The MC-AFIR does automated sampling of reaction pathways between two or more reactant molecules. The SC-AFIR performs automated generation of global or semiglobal reaction path network. The DS-AFIR finds a single path between given two structures. Exploration of minimum energy structures within the hypersurface in which two different electronic states degenerate, and an interface with the quantum mechanics/molecular mechanics method, are also described. A code termed SAFIRE will also be available, as a visualization software for complicated reaction path networks. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

An autocatalytic cycle in autoxidation of triethylborane.

An autocatalytic cycle was found in the mechanism of autoxidation of triethylborane using density functional theory calculations. The reaction starts with the generation of an ethyl radical via slow homolytic substitution. Fast radical propagation then takes place through a catalytic cycle in which the ethyl radical acts as a catalyst.

Full rate constant matrix contraction method for obtaining branching ratio of unimolecular decomposition.

The branching ratio of unimolecular decomposition can be evaluated by solving the rate equations. Recent advances in automated reaction path search methods have enabled efficient construction of the rate equations based on quantum chemical calculations. However, it is still difficult to solve the rate equations composed of hundreds or more elementary steps. This problem is especially serious when elementary steps that occur in highly different timescales coexist. In this article, we introduce an efficient approach to obtain the branching ratio from a given set of rate equations. It has been derived from a recently proposed rate constant matrix contraction (RCMC) method, and termed full-RCMC (f-RCMC). The f-RCMC gives the branching ratio without solving the rate equations. Its performance was tested numerically for unimolecular decomposition of C3 H5 and C4 H5 . Branching ratios obtained by the f-RCMC precisely reproduced the values obtained by numerically solving the rate equations. It took about 95 h to solve the rate equations of C4 H5 consisting of 234 elementary steps. In contrast, the f-RCMC gave the branching ratio in less than 1 s. The f-RCMC would thus be an efficient alternative of the conventional kinetic simulation approach. © 2016 Wiley Periodicals, Inc.

Kinetic Analysis for the Multistep Profiles of Organic Reactions: Significance of the Conformational Entropy on the Rate Constants of the Claisen Rearrangement.

The significance of kinetic analysis as a tool for understanding the reactivity and selectivity of organic reactions has recently been recognized. However, conventional simulation approaches that solve rate equations numerically are not amenable to multistep reaction profiles consisting of fast and slow elementary steps. Herein, we present an efficient and robust approach for evaluating the overall rate constants of multistep reactions via the recursive contraction of the rate equations to give the overall rate constants for the products and byproducts. This new method was applied to the Claisen rearrangement of allyl vinyl ether, as well as a substituted allyl vinyl ether. Notably, the profiles of these reactions contained 23 and 84 local minima, and 66 and 278 transition states, respectively. The overall rate constant for the Claisen rearrangement of allyl vinyl ether was consistent with the experimental value. The selectivity of the Claisen rearrangement reaction has also been assessed using a substituted allyl vinyl ether. The results of this study showed that the conformational entropy in these flexible chain molecules had a substantial impact on the overall rate constants. This new method could therefore be used to estimate the overall rate constants of various other organic reactions involving flexible molecules.