From Gas‐Phase to Liquid‐Water Chemical Reactions: The Fluorine Atom Plus Water Trimer System

The potential energy profile for the F+(H₂O)₃→HF+(H₂O)₂OH reaction has been investigated using the “gold standard” CCSD(T) method with correlation‐consistent basis sets up to cc‐pVQZ. Four different reaction pathways have been found and these are related, both geometrically and energetically. The entrance complexes F???(H₂O)₃ for all four reaction pathways are found lying ca. 7 kcal mol^?1 below the separated reactants F+(H₂O)₃. The four reaction barriers on their respective reaction coordinates lie ca. 4 kcal mol^?1 below the reactants. There are also corresponding exit complexes HF???(H₂O)₂OH, lying about 13 kcal mol^?1 below the separated products HF+(H₂O)₂OH. Compared with analogous F+(H₂O)₂ and F+H₂O reactions, the F+(H₂O)₃ reaction is somewhat similar to the former but qualitatively different from the latter. It may be reasonable to predict that the reactions between atomic fluorine and water tetramer (or even larger water clusters) may be similar to the F+(H₂O)₃ reaction.     

Theoretical Studies on Mechanism and Rate Constant of Gas Phase Hydrolysis of Glyoxal Catalyzed by Sulfuric Acid

The gas phase hydration of glyoxal (HCOCHO) in the presence of sulfuric acid (H2SO4) were studied by the high-level quantum chemical calculations with M06-2X and CCSD(T) theoretical methods and the conventional transition state theory (CTST). The mechanism and rate constant of the ve di erent reaction paths are consid-ered corresponding to HCOCHO+H₂O, HCOCHO+H₂O H₂O, HCOCHO H₂O+H₂O, HCOCHO+H₂O H₂SO₄ and HCOCHO H₂O+H₂SO₄. Results show that H₂SO₄ has a strong catalytic ability, which can signi cantly reduce the energy barrier for the hydration reaction of glyoxal. The energy barrier of hydrolysis of glyoxal in gas phase is lowered to 7.08 kcal/mol from 37.15 kcal/mol relative to pre-reactive complexes at the CCSD(T)/6-311++G(3df, 3pd)//M06-2X/6-311++G(3df, 3pd) level of theory. The rate constant of the H₂SO₄ catalyzed hydrolysis of glyoxal is 1.34×10^-11cm³/(molecule s), about 10¹³ higher than that involving catalysis by an equal number of water molecules, and is greater than the reaction rate of glyoxal reaction with OH radicals of 1.10×10^-11cm³/(molecule s) at the room temperature, indicating that the gas phase hydrolysis of glyoxal of H₂SO₄ catalyst is feasible and could compete with the reaction glyoxal+OH under certain atmospheric condi-tions. This study may provide useful information on understanding the mechanistic features of inorganic acid-catalyzed hydration of glyoxal for the formation of oligomer     

Is There an Entrance Complex for the F+NH3 Reaction?

Challenges associated with the theoretical and experimental kinetics of the F+NH₃→HF+NH₂ reaction suggest the need for a more‐precise potential surface. We have investigated the reactants and the products of the reaction, as well as the transition state and two complexes, with rather rigorous ab initio methods. The F????NH₃ complex existing in the entrance valley is predicted to lie 13.7 kcal mol^?1 below the reactants. A small classical barrier of 2.0 kcal mol^?1 separates this entrance well from products HF+NH₂. These results explain the observation by Persky of unprecedented inverse temperature dependence for the F+NH₃ rate constants. The strong hydrogen‐bonded complex FH????NH₂ exists in the exit valley, and with a binding energy of 9.9 kcal mol^?1 relative to separated products. The vibrational frequencies of all stationary points are predicted with the CCSD(T)/aug‐cc‐pVQZ method.     

The proton‐coupled proton transfer mechanism,H2O catalysis,and hydrogen tunneling effects in the reaction of HNCH2 with HCOOH in the interstellar medium

The reaction HNCH₂ + HCOOH → H₂NCH₂COOH is supposed to be an important reaction related to the possible origin of amino acids on the early Earth. We find that it has an energy barrier of 87.37 kcal mol⁻¹ obtained with MP2/6‐311+G** in the gas phase, but it is likely enhanced to occur in the interstellar medium (ISM) through a proton‐coupled proton transfer reaction, initiated by HNCH₂ coupled with H₂⁺, H₃⁺, or H₃⁺O. H₂⁺, H₃⁺, and H₃⁺O serve as a donor of energy in the coupled reactions. H⁺, which is a key species to the coupled reactions, further, plays a catalytic role in reducing a barrier up to 14.14 kcal mol⁻¹. In the coupled reaction with H₃⁺O, H₂O, which can seize, transport, and deliver a proton from HCOOH to H₂NCH₂⁺, reduces a barrier up to 14.96 kcal mol⁻¹. A significant hydrogen‐tunneling pathway is predicted by the temperature dependences of k_H^CVT/SCT, calculated using the small curvature tunneling (SCT) approximation and canonical variational transition state theory (CVT). Hydrogen tunneling is another important mechanism to make the reaction happen in the ISM. The achieved results can be applied to discuss the origin of amino acids from the materials of the Earth itself. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Direct dynamics study on the mechanism and the kinetics of the reaction of CH3NH2 with OH

A theoretical study of the mechanism and the kinetics for the hydrogen abstraction reaction of methylamine by OH radical has been presented at the CCSD(T)/6‐311 ++G(2d,2p)//CCSD/6‐31G(d) level of theory. Our theoretical calculations suggest a stepwise mechanism involving the formation of a prereactant complex in the entrance channel and a preproduct complex in the exit channel, for the two hydrogen abstraction channels involving the methyl and amine groups. For clarity, the diagram of potential for the reaction is given. The calculated standard reaction enthalpies are ?98.48 and ?76.50 kJ mol^?1 and barrier heights are 0.36 and 25.25 kJ mol^?1, respectively. The rate constants are evaluated by means of the improved canonical variational transition state theory with small‐curvature tunneling correction (ICVT/SCT) in the temperature range of 299–3000 K. The calculated results show that the rate constants at experimentally measured temperatures are in good agreement with the experimental values. It is shown that the calculated rate constants exhibit a non‐Arrhenius behavior. Moreover, the variational effect is obvious in the calculated temperature range. The dominant product channel is to form CH₂NH₂ and H₂O via hydrogen abstraction from the CH₃ group of CH₃NH₂ by OH in the calculated temperature range. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

The fragmentations of [M-H]- anions derived from underivatised peptides. The side-chain loss of H2S from Cys. A joint experimental and theoretical study

Loss of H₂S is the characteristic Cys side‐chain fragmentation of the [M? H]^? anions of Cys‐containing peptides. A combination of experiment and theory suggests that this reaction is initiated from the Cys enolate anion as follows: RNH‐^?C(CH₂SH)CONHR′ Ø [RNHC(?CH₂)CONHR′ (HS^?)] Ø [RNHC(?CH₂)CO‐HNR′‐H]^?+H₂S. This process is facile. Calculations at the HF/6‐31G(d)//AM1 level of theory indicate that the initial anion needs only ≥20.1 kcal mol^?1 of excess energy to effect loss of H₂S. Loss of CH₂S is a minor process, RNHCH(CH₂SH)CON^?‐R′ Ø RNHCH(CH₂S^?)CONHR′ Ø RNH ^?CHCONHR+CH₂S, requiring an excess energy of ≥50.2 kcal mol^?1. When Cys occupies the C‐terminal end of a peptide, the major fragmentation from the [M–H]^? species involves loss of (H₂S+CO₂). A deuterium‐labelling study suggests that this could either be a charge‐remote reaction (a process which occurs remote from and uninfluenced by the charged centre in the molecule), or an anionic reaction initiated from the C‐terminal CO₂^? group. These processes have barriers requiring the starting material to have an excess energy of ≥79.6 (charge‐remote) or ≥67.1 (anion‐directed) kcal mol^?1, respectively, at the HF/6‐31G(d)//AM1 level of theory. The corresponding losses of CH₂O and H₂O from the [M? H]^? anions of Ser‐containing peptides require ≥35.6 and ≥44.4 kcal mol^?1 of excess energy (calculated at the AM1 level of theory), explaining why loss of CH₂O is the characteristic side‐chain loss of Ser in the negative ion mode. Copyright © 2003 John Wiley & Sons, Ltd.     

The Beryllium Pentamer: Trailing an Uneven Sequence of Dissociation Energies

Recent high‐resolution spectroscopic studies by Merritt, Bondybey, and Heaven (Science 2009 , 324, 1548) have heightened the anticipation that small beryllium clusters will soon be observed in the laboratory. Beryllium clusters are important discrete models for the theoretical study of metals. The trigonal bipyramidal Be₅ molecule is studied using high‐level coupled cluster methods. We obtain the optimized geometry, atomization and dissociation energies, and vibrational frequencies. The c～CCSDT(Q) method is employed to compute the atomization and dissociation energies. In this approach, complete basis set (CBS) extrapolations at the CCSD(T) level of theory are combined with an additive correction for the effect of iterative triple and perturbative quadruple excitations. Harmonic vibrational frequencies are obtained using analytic gradients computed at the CCSD(T) level of theory. We report an atomization energy of 129.6 kcal mol^?1 at the trigonal bipyramid global minimum geometry. The Be₅→Be₄+Be dissociation energy is predicted to be 39.5 kcal mol^?1. The analogous dissociation energies for the smaller beryllium clusters are 64.0 kcal mol^?1 (Be₄→Be₃+Be), 24.2 kcal mol^?1 (Be₃→Be₂+Be), and 2.7 kcal mol^?1 (Be₂→Be+Be). The trigonal bipyramidal Be₅ structure has an equatorial–equatorial bond length of 2.000 Å and an axial–equatorial distance of 2.060 Å. Harmonic frequencies of 730, 611, 456, 583, 488, and 338 cm^?1 are obtained at the CCSD(T)/cc‐pCVQZ level of theory. Quadruple excitations are found to make noticeable contributions to the energetics of the pentamer, which exhibits a significant level of static correlation.     

Tetragermacyclobutadiene: Energetically Disfavored with Respect to Its Structural Isomers

Germanium has been a central feature in the renaissance of main‐group inorganic chemistry. Herein, we present the stationary‐point geometries of tetragermacyclobutadiene and its related isomers on the singlet potential energy surface at the CCSD(T)/cc‐pVTZ level of theory. Three of these 12 structures are reported for the first time and one of them is predicted to lie only 0.4 kcal mol^?1 above the previously reported global minimum. Focal‐point analyses has provided electronic energies at the CCSD(T) level of theory, which are extrapolated to the complete basis‐set limit and demonstrate the convergence behavior of the electronic energies with improving levels of theory and increasing basis‐set size. The lowest‐energy structure is the bicyclic structure, which lies 35 kcal mol^?1 below the “all‐Ge” cyclobutadiene structure. The reaction energies for the association of known Ge hydrides (e.g., digermene) to form Ge₄H₄ indicate that Ge₄H₄ could be observed experimentally. We investigate the bonding patterns by examining the frontier molecular orbitals. Our results demonstrate that: 1) the cyclic isomers of (GeH)₄ distort to maximize the mixing of the p orbitals that are involved in the π system of tetragermacyclobutadiene and 2) the lowest‐energy isomers exhibit unusual bonding arrangements (e.g., bridging H bonds) that maximize the nonbonding electron density at the Ge centers.     

Theoretical Study of the Reaction Formalhydrazone with Singlet Oxygen. Fragmentation of the C=N Bond,Ene Reaction and Other Processes

Photobiologic and synthetic versatility of hydrazones has not yet been established with ¹O₂ as a route to commonly encountered nitrosamines. Thus, to determine whether the “parent” reaction of formalhydrazone and ¹O₂ leads to facile C=N bond cleavage and resulting nitrosamine formation, we have carried out CCSD(T)//DFT calculations and analyzed the energetics of the oxidation pathways. A [2 + 2] pathway occurs via diradicals and formation of 3‐amino‐1,2,3‐dioxazetidine in a 16 kcal/mol^?1 process. Reversible addition or physical quenching of ¹O₂ occurs either on the formalhydrazone carbon for triplet diradicals at 2–3 kcal mol^?1, or on the nitrogen (N(3)) atom forming zwitterions at ~15 kcal/mol^?1, although the quenching channel by charge‐transfer interaction was not computed. The computations also predict a facile conversion of formalhydrazone and ¹O₂ to hydroperoxymethyl diazene in a low‐barrier ‘ene’ process, but no 2‐amino‐oxaziridine‐O‐oxide (perepoxide‐like) intermediate was found. A Benson‐like analysis (group increment calculations) on the closed‐shell species are in accord with the quantum chemical results.     

Nearly Degenerate Isomers of C(BH)2: Cumulene,Carbene, or Carbone?

The ground electronic state of C(BH)₂ exhibits both a linear minimum and a peculiar angle‐deformation isomer with a central B‐C‐B angle near 90°. Definitive computations on these species and the intervening transition state have been executed by means of coupled‐cluster theory including single and double excitations (CCSD), perturbative triples (CCSD(T)), and full triples with perturbative quadruples (CCSDT(Q)), in concert with series of correlation‐consistent basis sets (cc‐pVXZ, X=D, T, Q, 5, 6; cc‐pCVXZ, X=T, Q). Final energies were pinpointed by focal‐point analyses (FPA) targeting the complete basis‐set limit of CCSDT(Q) theory with auxiliary core correlation, relativistic, and non‐Born–Oppenheimer corrections. Isomerization of the linear species to the bent form has a minuscule FPA reaction energy of 0.02 kcal mol^?1 and a corresponding barrier of only 1.89 kcal mol^?1. Quantum tunneling computations reveal interconversion of the two isomers on a timescale much less than 1 s even at 0 K. Highly accurate CCSD(T)/cc‐pVTZ and composite c～CCSDT(Q)/cc‐pCVQZ anharmonic vibrational frequencies confirm matrix‐isolation infrared bands previously assigned to linear C(BH)₂ and provide excellent predictions for the heretofore unobserved bent isomer. Chemical bonding in the C(BH)₂ species was exhaustively investigated by the atoms‐in‐molecules (AIM) approach, molecular orbital plots, various population analyses, local mode vibrations and force constants, unified reaction valley analysis (URVA), and other methods. Linear C(BH)₂ is a cumulene, whereas bent C(BH)₂ is best characterized as a carbene with little carbone character. Weak B–B attraction is clearly present in the unusual bent isomer, but its strength is insufficient to form a CB₂ ring with a genuine boron–boron bond and attendant AIM bond path.     

Anionic CH⋅⋅⋅X− Hydrogen Bonds: Origin of Their Strength,Geometry, and Other Properties

CF₃H as a proton donor was paired with a variety of anions, and its properties were assessed by MP2/aug‐cc‐pVDZ calculations. The binding energy of monoanions halide, NO₃^?, formate, acetate, HSO₄^?, and H₂PO₄^? lie in the 12–17 kcal mol^?1 range, although F^? is more strongly bound, by 26 kcal mol^?1. Dianions SO₄^2? and HPO₄^2? are bound by 27 kcal mol^?1, and trianion PO₄^3? by 45 kcal mol^?1. When two O atoms are available on the anion, the CH???O^? H‐bond (HB) is usually bifurcated, although asymmetrically. The CH bond is elongated and its stretching frequency redshifted in these ionic HBs, but the shift is reduced in the bifurcated structures. Slightly more than half of the binding energy is attributed to Coulombic attraction, with smaller contributions from induction and dispersion. The amount of charge transfer from the anions to the σ*(CH) orbital correlates with many of the other indicators of bond strength, such as binding energy, CH bond stretch, CH redshift, downfield NMR spectroscopic chemical shift of the bridging proton, and density at bond critical points.     

A theoretical study on the mechanism and thermodynamics of ozone-water gas phase reaction

Ozone water reaction including a complex was studied at the MP2/6-311++G(d,p) and CCSD/6-311++G(2df,2p)//MP2/6-311++G(d,p) levels of theory. The interaction between water oxygen and central oxygen of ozone produces stable H₂O-O₃ complex with no barrier. With decomposition of this complex through H-abstraction by O₃ and O-abstraction by H₂O, three possible product channels were found. Intrinsic reaction coordinate, topological analyses of atom in molecule, and vibrational frequency calculation have been used to confirm the preferred mechanism. Thermodynamic data at T = 298.15 K and atmospheric pressure have been calculated. The results show that the production of hydrogen peroxide is the main reaction channel with ΔG = ?21.112 kJ mol^-1.     

Nature of the interaction between ammonia derivatives and carbon disulfide. A theoretical investigation

The interactions between NH₃, its methylated and chlorinated derivatives and CS₂ are investigated by ab initio CCSD(T) and density functional BLYP‐D3 methods. The CCSD(T)/aug‐cc‐pVTZ calculated interaction energies of complexes characterized by the S···N chalcogen bonds range between ?1.71 and ?2.78 kcal mol^?1. The S···N bonds are studied by atoms in molecules, natural bond orbital, and noncovalent interaction methods. The lack of correlation between the interaction energies of methylated amines complexes and the electrostatic potential results from the lone pair effect in aliphatic amines. Different structures of CS₂ complexed with ammonia derivatives, stabilized by other than the S···N chalcogen bonds, are also predicted. These structures are characterized by interaction energies ranging between 1.15 and 3.46 kcal mol^?1. The results show that the complexing ability of CS₂ is not very high but this molecule is able to attack the electrophilic or nucleophilic sites of a guest molecule.     

A Theoretical Study of the Kinetics of the Hydrogen Atom Abstraction Reactions from Cyclopropane by H,O (3P), and Cl (2P3/2) Atoms and OH Radicals

The rate constants of the H‐abstraction reactions from cyclopropane by H, O (³P), Cl (²P_3/2), and OH radicals have been calculated over the temperature range of 250?2500 K using two different levels of theory. Calculations of optimized geometrical parameters and vibrational frequencies are performed using the MP2 method combined with the cc‐pVTZ basis set and the 6–311++G(d,p) basis set. Single‐point energy calculations have been carried out with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (perturbatively) electron excitations CCSD(T) using either the cc‐pVTZ, aug‐cc‐pVTZ, and aug‐cc‐pVQZ basis sets or the 6–311++G(3df,3pd) basis set. The CCSD(T) calculated potential energies have been extrapolated to the complete basis limit (CBS) limit. The Full Configuration Interaction (FCI) energies have been also estimated using the continued‐fraction approximation as proposed by Goodson (J. Chem. Phys., 2002, 116, 6948–6956). Canonical transition‐state theory combined with an Eckart tunneling correction has been used to predict the rate constants as a function of temperature using two kinetic models (direct abstraction or complex mechanism) at two levels of theory (CCSD(T)‐cf/CBS//MP2/cc‐pVTZ and CCSD(T)‐cf/6–311++G(3df,3pd)//MP2/6–311++G(d,p)). The calculated kinetic parameters are in reasonable agreement with their literature counterparts for all reactions. In the light of these trends, the use of the Pople‐style basis sets for studying the reactivity of other systems such as larger cycloalkanes or halogenated cycloalkanes is recommended because the 6–311++G(3df,3pd) basis set is less time consuming than the aug‐cc‐pVQZ basis set. Based on our calculations performed at the CCSD(T)‐cf/CBS//MP2/cc‐pVTZ level of theory, the standard enthalpy of formation at 298 K for the cyclopropyl radical has been reassessed and its value is (290.5 ± 1.6) kJ mol^?1.     

Characterization of the potential energy surfaces of two small but challenging noncovalent dimers: (P2)2 and (PCCP)2

This work characterizes eight stationary points of the P₂ dimer and six stationary points of the PCCP dimer, including a newly identified minimum on both potential energy surfaces. Full geometry optimizations and corresponding harmonic vibrational frequencies were computed with the second‐order Møller–Plesset (MP2) electronic structure method and six different basis sets: aug‐cc‐pVXZ, aug‐cc‐pV(X+d)Z, and aug‐cc‐pCVXZ where X = T, Q. A new L‐shaped structure with C₂ symmetry is the only minimum for the P₂ dimer at the MP2 level of theory with these basis sets. The previously reported parallel‐slipped structure with C_2h symmetry and a newly identified cross configuration with D₂ symmetry are the only minima for the PCCP dimer. Single point energies were also computed using the canonical MP2 and CCSD(T) methods as well as the explicitly correlated MP2‐F12 and CCSD(T)‐F12 methods and the aug‐cc‐pVXZ (X = D, T, Q, 5) basis sets. The energetics obtained with the explicitly correlated methods were very similar to the canonical results for the larger basis sets. Extrapolations were performed to estimate the complete basis set (CBS) limit MP2 and CCSD(T) binding energies. MP2 and MP2‐F12 significantly overbind the P₂ and PCCP dimers relative to the CCSD(T) and CCSD(T)‐F12 binding energies by as much as 1.5 kcal mol^?1 for the former and 5.0 kcal mol^?1 for the latter at the CBS limit. The dominant attractive component of the interaction energy for each dimer configuration was dispersion according to several symmetry‐adapted perturbation theory analyses. © 2014 Wiley Periodicals, Inc.     

Ab Initio Study on the Mechanism of the Atmospheric Reaction OH+O3→HO2+O2

The atmospheric reaction (1) OH + O₃→HO₂ + O₂ was investigated theoretically by using MP2, QCISD, QCISD(T), and CCSD(T) methods with various basis sets. At the highest level of theory, namely, QCISD, the reaction is direct, with only one transition state between reactants and products. However, at the MP2 level, the reaction proceeds through a two‐step mechanism and shows two transition states, TS1 and TS2 , separated by an intermediate, Int . The different methodologies employed in this paper consistently predict the barrier height of reaction (1) to be within the range 2.16–5.11 kcal mol^?1, somewhat higher than the experimental value of 2.0 kcal mol^?1.     

Atomistic understanding of the C·T mismatched DNA base pair tautomerization via the DPT: QM and QTAIM computational approaches

It was established that the cytosine·thymine (C·T) mismatched DNA base pair with cis‐oriented N1H glycosidic bonds has propeller‐like structure (|N3C4C4N3| = 38.4°), which is stabilized by three specific intermolecular interactions–two antiparallel N4H…O4 (5.19 kcal mol^?1) and N3H…N3 (6.33 kcal mol^?1) H‐bonds and a van der Waals (vdW) contact O2…O2 (0.32 kcal mol^?1). The C·T base mispair is thermodynamically stable structure (ΔG_int = ?1.54 kcal mol^?1) and even slightly more stable than the A·T Watson–Crick DNA base pair (ΔG_int = ?1.43 kcal mol^?1) at the room temperature. It was shown that the C·T ? C*·T* tautomerization via the double proton transfer (DPT) is assisted by the O2…O2 vdW contact along the entire range of the intrinsic reaction coordinate (IRC). The positive value of the Grunenberg's compliance constants (31.186, 30.265, and 22.166 Å/mdyn for the C·T, C*·T*, and TS_{C·T ? C*·T*}, respectively) proves that the O2…O2 vdW contact is a stabilizing interaction. Based on the sweeps of the H‐bond energies, it was found that the N4H…O4/O4H…N4, and N3H…N3 H‐bonds in the C·T and C*·T* base pairs are anticooperative and weaken each other, whereas the middle N3H…N3 H‐bond and the O2…O2 vdW contact are cooperative and mutually reinforce each other. It was found that the tautomerization of the C·T base mispair through the DPT is concerted and asynchronous reaction that proceeds via the TS_{C·T ? C*·T*} stabilized by the loosened N4? H? O4 covalent bridge, N3H…N3 H‐bond (9.67 kcal mol^?1) and O2…O2 vdW contact (0.41 kcal mol^?1). The nine key points, describing the evolution of the C·T ? C*·T* tautomerization via the DPT, were detected and completely investigated along the IRC. The C*·T* mispair was revealed to be the dynamically unstable structure with a lifetime 2.13·× 10^?13 s. In this case, as for the A·T Watson–Crick DNA base pair, activates the mechanism of the quantum protection of the C·T DNA base mispair from its spontaneous mutagenic tautomerization through the DPT. © 2013 Wiley Periodicals, Inc.     

Proton‐Conductive 3D LnIII Metal–Organic Frameworks for Formic Acid Impedance Sensing

Metal‐organic frameworks (MOFs) as new classes of proton‐conducting materials have been highlighted in recent years. Nevertheless, the exploration of proton‐conducting MOFs as formic acid sensors is extremely lacking. Herein, we prepared two highly stable 3D isostructural lanthanide(III) MOFs, {(M(μ₃‐HPhIDC)(μ₂‐C₂O₄)_0.5(H₂O))?2 H₂O}_n (M=Tb ( ZZU‐1 ); Eu ( ZZU‐2 )) (H₃PhIDC=2‐phenyl‐1H‐imidazole‐4,5‐dicarboxylic acid), in which the coordinated and uncoordinated water molecules and uncoordinated imidazole N atoms play decisive roles for the high‐performance proton conduction and recognition ability for formic acid. Both ZZU‐1 and ZZU‐2 show temperature‐ and humidity‐dependent proton‐conducting characteristics with high conductivities of 8.95×10^?4 and 4.63×10^?4 S cm^‐1 at 98 % RH and 100 °C, respectively. Importantly, the impedance values of the two MOF‐based sensors decrease upon exposure to formic acid vapor generated from formic aqueous solutions at 25 °C with good reproducibility. By comparing the changes of impedance values, we can indirectly determine the concentration of HCOOH in aqueous solution. The results showed that the lowest detectable concentrations of formic acid aqueous solutions are 1.2×10^?2 mol L^?1 by ZZU‐1 and 2.0×10^?2 mol L^?1 by ZZU‐2 . Furthermore, the two sensors can distinguish formic acid vapor from interfering vapors including MeOH, N‐hexane, benzene, toluene, EtOH, acetone, acetic acid and butane. Our research provides a new platform of proton‐conductive MOFs‐based sensors for detecting formic acid.     

Stacking of Benzene with Metal Chelates: Calculated CCSD(T)/CBS Interaction Energies and Potential‐Energy Curves

Accurate values for the energies of stacking interactions of nickel‐ and copper‐based six‐membered chelate rings with benzene are calculated at the CCSD(T)/CBS level. The results show that calculations made at the ωB97xD/def2‐TZVP level are in excellent agreement with CCSD(T)/CBS values. The energies of [Cu(C₃H₃O₂)(HCO₂)] and [Ni(C₃H₃O₂)(HCO₂)] chelates stacking with benzene are ?6.39 and ?4.77 kcal mol^?1, respectively. Understanding these interactions might be important for materials with properties that are dependent on stacking interactions.     

Mechanistic study on the atmospheric formation of acid rain base on the sulfur dioxide

The reaction pathways of acid rain formation from reaction of sulfur dioxide vapor and water vapor on the singlet potential energy surface have been investigated theoretically. The calculated results show that the reactants are initially associated with the adduct SO₂–H₂O through a barrier less process. Subsequently, via a variety of transformations of isomer SO₂–H₂O, three kinds of products H₂SO₃, SO₃ + H₂, and H₂O₂ + ³SO are obtained. The cleavage and formation of the chemical bonds in the reaction pathways have been discussed using the structural parameters. Also, by means of the transition states and their connected intermediates or products at the CCSD(T)//B3LYP level, mechanism of H₂O + SO₂ reaction on the singlet potential energy surface are plotted. The calculation results show that the most suitable reaction pathways are the formation of H₂SO₃. Finally, the rate constants have been calculated only for these suitable pathways by the RRKM and TST theories at temperature range of 250–2500 K.