Screening of novel serum biomarkers for gastric cancer in coastal populations using a protein microarray.

Gastric cancer (GC) has high rates of morbidity and mortality, and this phenomenon is particularly evident in coastal regions where local dietary habits favor the consumption of pickled foods such as salted fish and vegetables. In addition, the diagnosis rate of GC remains low due to the lack of diagnostic serum biomarkers. Therefore, in this study, we aimed to identify potential serum GC biomarkers for use in clinical practice. To identify candidate biomarkers of GC, 88 serum samples were first screened using a high-throughput protein microarray to measure the levels of 640 proteins. Then, 333 samples were used to validate the potential biomarkers using a custom antibody chip. ELISA, western blot, and immunohistochemistry were then used to verify the expression of the target proteins. Finally, logistic regression was performed to select serum proteins for the diagnostic model. As a result, five specific differentially expressed proteins, TGFß RIII, LAG-3, carboxypeptidase A2, Decorin and ANGPTL3, were found to have the ability to distinguish GC. Logistic regression analysis showed that the combination of carboxypeptidase A2 and TGFß RIII had superior potential for diagnosing GC (area under the ROC curve [AUC] = 0.801). The results suggested that these five proteins alone and the combination of carboxypeptidase A2 and TGFß RIII may be used as serum markers for the diagnosis of GC.

Computational study on the drug resistance mechanism of HCV NS3 protease to BMS-605339.

NS3 protease plays a vital role in the replication of the hepatitis C virus (HCV). BMS-605339 is a novel linear tetra-peptide α-ketoamide inhibitor of NS3 protease and shows specificity for HCV NS3 protease genotype 1a and genotype 1b. Mutation at the key site 168 of the HCV NS3 protease can induce resistance to BMS-605339, which greatly affects the antiviral therapy efficacy to hepatitis C. In the present study, we employed molecular dynamics simulations, free energy calculations, and free energy decomposition to explore the drug resistance mechanism of BMS-605339 due to the three representative mutations D168C/Y/V. The free energy decomposition analysis indicates that the decrease in the binding affinity is mainly attributed to the decrease in both van der Waals and electrostatic interactions. After detailed analysis of our calculated results, we observed that the break of the salt bridge between residues 155 and 168 caused by the mutations D168C/Y/V is the original reason for the decrease in the binding ability between BMS-605339 and the mutant NS3 proteases. The obtained results will reveal the drug resistance mechanism between BMS-605339 and the mutant NS3 proteases, and provide valuable clue for designing novel and more potent drugs to HCV NS3 protease.

Molecular modeling study on the drug resistance mechanism of NS5B polymerase to TMC647055.

NS5B polymerase plays an important role in viral replication machinery. TMC647055 (TMC) is a novel and potent non-nucleoside inhibitor of the HCV NS5B polymerase. However, mutations that result in drug resistance to TMC have been reported. In this study, we used molecular dynamics (MD) simulations, binding free energy calculations, and free energy decomposition to investigate the drug resistance mechanism of HCV to TMC resulting from L392I, P495T, P495S, and P495L mutations in NS5B polymerase. From the calculated results we determined that the decrease in the binding affinity between TMC and NS5B(L392I) polymerase is mainly caused by the extra methyl group at the CB atom of Ile. The polarity of the side-chain of residue 495 has no distinct influence on residue 495 binding with TMC, whereas the smaller size of the side-chain of residue 495 causes a substantial decrease in the van der Walls interaction between TMC and residue 495. Moreover, the longer length of the side-chain of residue 495 has a significant effect on the electrostatic interaction between TMC and Arg-503. Finally, we performed the same calculations and detailed analysis on other 3 mutations (L392V, P495V, and P495I). The results further confirmed our conclusions. The computational results not only reveal the drug resistance mechanism between TMC647055 and NS5B polymerase, but also provide valuable information for the rational design of more potent non-nucleoside inhibitors targeting HCV NS5B polymerase.

Computational study on the molecular mechanisms of drug resistance of Narlaprevir due to V36M, R155K, V36M+R155K, T54A, and A156T mutations of HCV NS3/4A protease.

Narlaprevir is a novel NS3/4A protease inhibitor of hepatitis C virus (HCV), and it has been tested in a phase II clinical trial recently. However, distinct drug-resistance of Narlaprevir has been discovered. In our study, the molecular mechanisms of drug-resistance of Narlaprevir due to the mutations V36M, R155K, V36M+R155K, T54A, and A156T of NS3/4A protease have been investigated by molecular dynamics (MD) simulations, free energy calculations, and free energy decomposition analysis. The predicted binding free energies of Narlaprevir towards the wild-type and five mutants show that the mutations V36M, R155K, and T54A lead to low-level drug resistance and the mutations V36M+R155K and A156T lead to high-level drug resistance, which is consistent with the experimental data. The analysis of the individual energy terms indicates that the van der Waals contribution is important for distinguishing the binding affinities of these six complexes. These findings again show that the combination of different molecular modeling techniques is an efficient way to interpret the molecular mechanism of drug-resistance. Our work mainly elaborates the molecular mechanism of drug-resistance of Narlaprevir and further provides valuable information for developing novel, safer, and more potent HCV antiviral drugs in the near future.

The C-terminal selenenylsulfide of extracellular/non-reduced thioredoxin reductase endows this protein with selectivity to small-molecule electrophilic reagents under oxidative conditions.

Mammalian cytosolic thioredoxin reductase (TrxR1) serves as an antioxidant protein by transferring electrons from NADPH to various substrates. The action of TrxR1 is achieved via reversible changes between NADPH-reduced and non-reduced forms, which involves C-terminal selenolthiol/selenenylsulfide exchanges. TrxR1 may be released into extracellular environment, where TrxR1 is present mainly in the non-reduced form with active-site disulfide and selenenylsulfide bonds. The relationships between extracellular TrxR1 and tumor metastasis or cellular signaling have been discovered, but there are few reports on small-molecule compounds in targeted the non-reduced form of TrxR1. Using eight types of small-molecule thiol-reactive reagents as electrophilic models, we report that the selenenylsulfide bond in the non-reduced form of TrxR1 functions as a selector for the thiol-reactive reagents at pH 7.5. The non-reduced form of TrxR1 is resistant to hydrogen peroxide/oxidized glutathione, but is sensitive to certain electrophilic reagents in different ways. With 5,5'-dithiobis-(2-nitrobenzoic acid) (DTNB) and S-nitrosoglutathione (GSNO), the polarized selenenylsulfide bond breaks, and selenolate anion donates electron to the dynamic covalent bond in DTNB or GSNO, forming TNB-S-Se-TrxR1 complex or ON-Se-TrxR1 complex. The both complexes lose the ability to transfer electrons from NADPH to substrate. For diamide, the non-reduced TrxR1 actually prevents irreversible damage by this oxidant. This is consistent with the regained activity of TrxR1 through removal of diamide via dialysis. Diamide shows effective in the presence of human cytosolic thioredoxin (hTrx1), Cys residue(s) of which is/are preferentially affected by diamide to yield disulfide, hTrx1 dimer and the mixed disulfide between TrxR1-Cys497/Sec498 and hTrx1-Cys73. In human serum samples, the non-reduced form of TrxR1 exists as dithiothreitol-reducible polymer/complexes, which might protect the non-reduced TrxR1 from inactivation by certain electrophilic reagents under oxidative conditions, because cleavage of these disulfides can lead to regain the activity of TrxR1. The details of the selective response of the selenenylsulfide bond to electrophilic reagents may provide new information for designing novel small-molecule inhibitors (drugs) in targeted extracellular/non-reduced TrxR1.

Zwitterionic base-stabilized digermadistannacyclobutadiene and tetragermacyclobutadiene.

The syntheses of a zwitterionic base-stabilized digermadistannacyclobutadiene and tetragermacyclobutadiene supported by amidinates and low-valent germanium amidinate substituents are described. The reaction of the amidinate Ge(I) dimer, [LGe:]2 (1, L=PhC(NtBu)2 ), with two equivalents of the amidinate tin(II) chloride, [LSnCl] (2), and KC8 in tetrahydrofuran (THF) at room temperature afforded a mixture of the zwitterionic base-stabilized digermadistannacyclobutadiene, [L2 Ge2 Sn2 L'2 ] (3; L'=LGe:), and the bis(amidinate) tin(II) compound, [L2 Sn:] (4). Compound 3 can also be prepared by the reaction of 1 with [L(Ar) SnCl] (5, L(Ar) =tBuC(NAr)2 , Ar=2,6-iPr2 C6 H3 ) in THF at room temperature. Moreover, the reaction of 1 with the "onio-substituent transfer" reagent [4-NMe2 -C5 H4 NSiMe3 ]OTf (8) in THF and 4-(N,N-dimethylamino)pyridine (DMAP) at room temperature afforded a mixture of the zwitterionic base-stabilized tetragermacyclobutadiene, [L4 Ge6 ] (9), the amidinium triflate, [PhC(NHtBu)2 ]OTf (10), and Me3 SiSiMe3 (11). X-ray structural data and theoretical studies show conclusively that compounds 3 and 9 have a planar and rhombic charge-separated structure. They are also nonaromatic.

Photodissociation mechanisms of the CO2(2+) dication studied using multi-state multiconfiguration second-order perturbation theory.

Employing the multi-state multiconfiguration second-order perturbation theory (MS-CASPT2) and complete active space self-consistent field (CASSCF) methods, the geometries, relative energies (T(v)') to the ground state (X(3)Σg(-)), adiabatic excited energies, and photodissociation mechanisms and corresponding kinetic energy releases for the lower-lying 14 electronic states of the CO2 (2+) ion are studied. The T(v)' values are calculated at the experimental geometry of the ground state CO2 molecule using MS-CASPT2 method and highly close to the latest threshold photoelectrons coincidence and time-of-flight photoelectron photoelectron coincidence spectrum observations. The O-loss dissociation potential energy curves (PECs) for these 14 states are drawn using MS-CASPT2 partial optimization method at C(∞v) symmetry with one C-O bond length ranging from 1.05 to 8.0 Å. Those 14 states are confirmed to be correlated to the lowest four dissociation limits [CO(+)(X(2)Σ(+)) + O(+)((4)S(u)), CO(+)(A(2)Π) + O(+)((4)S(u)), CO(+)(X(2)Σ(+)) + O(+)((2)D(u)), and CO(+)(X(2)Σ(+)) + O(+)((2)P(u))] by analyzing Coulomb interaction energies, charges, spin densities, and bond lengths for the geometries at the C-O bond length of 8.0 Å. On the basis of these 14 MS-CASPT2 PECs, several state/state pairs are selected to optimize the minimum energy crossing points (MECPs) at the CASSCF level. And then the CASSCF spin-orbit couplings and CASPT2 state/state energies are calculated at these located MECPs. Based on all of the computational results, the photodissociation mechanisms of CO2(2+) are proposed. The relationships between the present theoretical studies and the previous experiments are discussed.

The Bergman cyclizations of the enediyne and its N-substituted analogs using multiconfigurational second-order perturbation theory.

The Bergman cyclizations of the enediyne and its four N-substituted analogs [(Z)-pent-2-en-4-ynenitrile, 3-azahex-3-en-1,5-diyne, malenotrile, and 3,4-azahex-3-en-1,5-diyne] have been studied using the complete active space self-consistent field and multiconfigurational second-order perturbation theory methods in conjunction with the atomic natural orbital basis sets. The geometries and energies of the reactants, transition states, and products along both the S(0) (the ground state) and T(1) (the lowest-lying triplet state) potential energy surfaces (PESs) were calculated. The calculated geometries are in good agreement with the available experimental data. The distance between two terminal carbons in enediyne, which was considered as an important parameter governing the Bergman cyclization, was predicted to be 4.319 Å, in agreement with the experimental value of 4.321 Å. Our calculations indicate that the replacements of the terminal C atom(s) or the middle C atom(s) in the C=C bond by the N atom(s) increase or decrease the energy barrier values, respectively. There exist stable ring biradical products on the T(1) PESs for the five reactions. However, on the S(0) PESs the ring biradical products exist only for the reactions of enediyne, (Z)-pent-2-en-4-ynenitrile, and 3-azahex-3-en-1,5-diyne.

Computational insights into the selectivity mechanism of APP-IP over matrix metalloproteinases.

In this work, selectivity mechanism of APP-IP inhibitor (ß-amyloid precursor protein-derived inhibitory peptide) over matrix metalloproteinases (MMPs including MMP-2, MMP-7, MMP-9 and MMP-14) was investigated by molecular modeling methods. Among MMPs, MMP-2 is the most favorable one for APP-IP interacting based on our calculations. The predicted binding affinities can give a good explanation of the activity difference of inhibitor APP-IP. In Comparison with MMP-2/APP-IP complex, the side chain of Tyr214(MMP-7) makes the binding pocket so shallow that the whole side chain of Tyr3(APP-IP) can not be fully embraced, thus unfavorable for the N-terminal of APP-IP binding to MMP-7. The poor selectivity of APP-IP toward MMP-9 is mainly related with the decrease of interaction between the APP-IP C-terminal and MMP-9 due to the bulky side chains of Pro193 and Gln199, which is in agreement with experiment. The mutations at residues P193A and Q199G of MMP-9 alternate the binding pattern of the C-terminal of APP-IP by forming two new hydrogen bonds and hydrophobic interactions with MMP-9. The mutants favor the binding affinity of MMP-9 largely. For MMP-14/APP-IP, the large steric effect of Phe204(MMP-14) and the weak contributions of the polar residues Asn231(MMP-14) and Thr190(MMP-14) could explain why MMP-14 is non-selective for APP-IP interacting. Here, the molecular modeling methods were successfully employed to explore the selective inhibitor of MMPs, and our work gives valuable information for future rational design of selective peptide inhibitors toward individual MMP.

Promising and efficient lignin degradation versatile strategy based on DFT calculations.

The extraction of higher-value products from lignin degradations under mild conditions is a challenge. Previous research reported efficient two-step oxidation and reduction strategies for lignin degradation, which has great significance to lignin degradation. In this paper, the mechanism about the C-O bond cleavage of lignin with and without Cα oxidations has been studied systematically. Our calculation results show that the degradation of anionized lignin with Cα oxidations is kinetically and thermodynamically feasible. In addition, the calculations predict that the anionized lignin compounds without Cα oxidation also could be degraded under mild conditions. Moreover, we propose special lignin catalytic degradation systems containing the characteristic structure of "double hydrogen bonds." The double hydrogen bonds structure could further decrease the energy barriers of the C-O bond cleavage reaction. This provides a versatile strategy to design novel lignin degradation.

Coulomb explosion of nitrogen and oxygen molecules through non-Coulombic states.

We have systematically studied Coulomb explosion of nitrogen and oxygen molecules in intense 8 and 24 fs laser pulses. In the experiment, we explicitly separated all explosion pathways through coincident measurements. The high resolution kinetic energy releases (KERs) and the exotic angular distributions of atomic ions provide direct evidence that Coulomb explosion occurs through non-Coulombic states. In the theory, we calculated dissociation potential energy curves (PECs) of nitrogen and oxygen molecules and their multicharged molecular ions using multiconfiguration second-order perturbation theory. The results indicate that Coulomb potentials are close to the accurate PECs of multicharged molecular ions only when the internuclear distance is larger than 3 Å. In comparison with the experimental observations and the theoretical calculations, we determined the internuclear distance when Coulombic explosion occurs. It is near the equilibrium distance of the neutral molecules in the case of 8 fs laser pulses and expands gradually with the increase of the charge state of the molecular ions in the case of 24 fs laser pulses.

Insight into the drug resistance mechanisms of GS-9669 caused by mutations of HCV NS5B polymerase via molecular simulation.

GS-9669 is a non-nucleos(t)ide inhibitor (NNI) binding to the thumb site II of the Hepatitis C virus (HCV) NS5B polymerase and has advanced into phase II trials. To clarify the drug resistance mechanisms of GS-9669 caused by M423T/I/V, L419M, R422K, and I482L mutations of NS5B polymerase (GT1b) and the receptor-ligand interactions during the binding process, a series of molecular simulation methods including molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed for the wild-type (WT) and six mutant NS5B/GS-9669 complexes. The calculated results indicate that the binding free energies of the mutant systems are less negative than that of the WT system, indicating that these mutations will indeed cause NS5B to produce different degrees of resistance to GS-9669. The mutation-induced drug resistances are mainly caused by the loss of binding affinities of Leu419 and Trp528 with GS-9669 or the formation of multiple solvent bridges. Moreover, the ASMD calculations show that GS-9669 binds to the thumb II sites of the seven NS5B polymerases in distinct pathways without any obvious energy barriers. Although the recognition methods and binding pathways are distinct, the binding processes of GS-9669 with the WT and mutant NS5B polymerases are basically controlled thermodynamically. This study clearly reveals the resistance mechanisms of GS-9669 caused by M423T/I/V, L419M, R422K, and I482L mutations of HCV NS5B polymerase and provides some valuable clues for further optimization and design of novel NS5B inhibitors.

Computational insight into the mechanisms of action and selectivity of Afraxis PAK inhibitors.

Aim: The p21-activated kinases (PAKs) are involved in many important biological activity regulations. FRAX019, FRAX414, FRAX597, FRAX1036 and G-5555 were identified as PAKs inhibitors. Their detailed inhibitory mechanisms deserve further investigation. Results: Molecular dynamics simulations and further calculations for the PAK1/inhibitor and PAK4/inhibitor complexes indicate that their binding free energies are basically consistent with the trend of experimental activity data. Conclusion: The anchoring of residues Leu347PAK1 and Leu398PAK4 is the structural basis for designing Afraxis PAK inhibitors. This study discloses the inhibitory mechanisms of FRAX019, FRAX414, FRAX597, FRAX1036 and G-5555 toward PAK1 and PAK4 and some clues to enhance kinase activities and selectivities, which will provide valuable information to the development of more potent and selective PAK inhibitors.

A Study on Modified Bitumen with Metal Doped Nano-TiO2 Pillared Montmorillonite.

TiO2 pillared montmorillonite (T/M) modifiers have been studied to alleviate the aging of asphalt pavement and degrade automobile exhaust, but the photocatalytic activity of ordinary TiO2 is not good enough. In this study, in order to improve the photocatalytic performance of T/M, different metal (Ce, Cu, Fe) doped modifiers were prepared based on T/M. Metal doped TiO2 pillared montmorillonite was prepared by the sol-gel method. The modifier was characterized by X-ray diffraction (XRD) and an Ultraviolet visible (UV-Vis) spectrophotometer. The results show that TiO2 with different metal ions successfully entered into the layer of organic montmorillonite (OMMT) to form a pillared structure. Compared with the undoped TiO2 pillared montmorillonite (T/M), the optical absorption edge of the metal doped TiO2 pillared montmorillonite has an obvious red shift. In addition, the influences of the different content of modifiers on the properties of the original bitumen and catalytic capacities for automobile exhaust were also investigated. The results show that Ce doped TiO2 (Ce-T/M) pillared montmorillonite has the best improvement in high temperature performance and ultraviolet (UV) resistance of bitumen. In the experiment of automobile exhaust degradation, the degradation law of NO and HC showed Cu-T/M > Ce-T/M > Fe-T/M > T/M. These three kinds of metal ions can effectively improve the photocatalytic degradation efficiency of T/M.

Low-lying electronic states and Cl-loss photodissociation of the C2H3Cl+ ion studied using multiconfiguration second-order perturbation theory.

We studied the 1(2)A' '(X2A' '), 1(2)A' (A2A'), 2(2)A' ' (B2A' '), and 2(2)A' (C2A') states of the C2H3Cl+ ion using the complete active space self-consistent field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) methods. For the four ionic states, we calculated the equilibrium geometries, adiabatic (T0) and vertical (Tv) excitation energies, and relative energies (Tv') at the geometry of the molecule at the CASPT2 level and the Cl-loss dissociation potential energy curves (PECs) at the CASPT2//CASSCF level. The computed oscillator strength f value for the X2A' ' <-- A2A' transition is very small, which is in line with the experimental fact that the A state has a long lifetime. The CASPT2 geometry and T0 value for the A2A' state are in good agreement with experiment. The CASPT2 Tv' values for the A2A', B2A' ', and C2A' states are in good agreement with experiment. The Cl-loss PEC calculations predict that the X2A' ', A2A', and C2A' states correlate to C2H3+ (XA1) and the BA' ' state to C2H3+ (1A' ') (the B2A' ' and C2A' PECs cross at R(C-Cl) approximately 2.24 A). Our calculations indicate that at 357 nm the X2A' ' state can undergo a transition to B2A' ' followed by a predissociation of B2A' ' by the repulsive C2A' state (via the B/C crossing), leading to C2H3+ (X1A1), and therefore confirm the experimentally proposed pathway for the photodissociation of X2A' ' at 357 nm. Our CASPT2 D0 calculations support the experimental fact that the X state does not undergo dissociation in the visible spectral region and imply that a direct dissociation of the A state to C2H3+ (X1A1) is energetically feasible.

Computational study on the drug resistance mechanism of HCV NS5B RNA-dependent RNA polymerase mutants V494I, V494A, M426A, and M423T to Filibuvir.

Filibuvir, a potent non-nucleoside inhibitor of the hepatitis C virus (HCV) NS5B RNA-dependent RNA polymerase (RdRp), has shown great promise in phase IIb clinical trial. However, drug resistant mutations towards Filibuvir have been identified. In the present study, the drug resistance mechanism of wild-type (WT) and mutant NS5B polymerases (including V494I, V494A, M426A, and M423T) toward Filibuvir was investigated by molecular modeling methods. The predicted binding free energy of these five complexes is highly consistent with the experimental EC50 values of Filibuvir to the wild-type and mutant NS5B RdRps, V494I

Divalent metals can reside on bonds in fullerenes.

DFT calculations consistently suggest that a lanthanide will sit on either the 6/6 bond inside C60 having a divalent state or the hexagonal center having a trivalent state. Some lanthanides can stay only above the 6/6 bond inside C60 to form stabilized structures, despite the greatly reduced metal-cage coordination numbers. The preference for C-C bonds by a divalent metal has been confirmed by revisiting the structures of Yb@C2v(3)-C80, Yb@Cs(6)-C82 and Yb@C2v(9)-C82, for which the calculations suggest that the Yb atoms are indeed situated above the C-C bonds, close to the reported structures obtained by single crystal XRD experiments. The result will guide the characterization of structures and electronic structures of endohedral metallofullerenes (EMFs), especially C60 mono-EMFs, in the future.

Dissociation of the OCS+ ion in low-lying electronic states studied using multiconfiguration second-order perturbation theory.

Complete active space self-consistent-field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) calculations with atomic natural orbital basis sets were performed to investigate the S-loss direct dissociation of the 1 2Pi(X 2Pi), 2 2Pi(A 2Pi), 1 2Sigma+(B 2Sigma+), 1 4Sigma-, 1 2Sigma-, and 1 2Delta states of the OCS+ ion and the predissociations of the 1 2Pi, 2 2Pi, and 1 2Sigma+ states. Our calculations indicate that the S-loss dissociation products of the OCS(+) ion in the six states are the ground-state CO molecule plus the S+ ion in different electronic states. The CASPT2//CASSCF potential energy curves were calculated for the S-loss dissociation from the six states. The calculations indicate that the dissociation of the 1 4Sigma- state leads to the CO + S+ (4Su) products representing the first dissociation limit; the dissociations of the 1 2Pi, 1 2Sigma-, and 1 2Delta states lead to the CO + S+(2Du) products representing the second dissociation limit; and the dissociations of the 2 2Pi and 1 2Sigma+ states lead to the CO + S+(2Pu) products representing the third dissociation limit. Seams of the 1 2Pi-1 4Sigma-, 2 2Pi-1 4Sigma-, 2 2Pi-1 2Sigma-, 2 2Pi-1 2Delta, and 1 2Sigma(+)-1 4Sigma- potential energy surface intersections were calculated at the CASPT2 level, and the minima along the seams were located. The calculations indicate that within the experimental energy range (15.07-16.0 eV) the 2 2Pi(A 2Pi) state can be predissociated by 1 4Sigma- forming the S+(4Su) ion and can undergo internal conversion to 1 2Pi followed by the direct dissociation of 1 2Pi forming S+(2Du) and that within the experimental energy range (16.04-16.54 eV) the 1 2Sigma+(B 2Sigma+) state can be predissociated by 1 4Sigma- forming the S+(4Su) ion and can undergo internal conversion to 2 2Pi followed by the predissociation of 2 2Pi by 1 2Sigma- and 1 2Delta forming the S+(2Du) ion. These indications are in line with the experimental fact that both the 4Su and 2Du states of the S+ ion can be formed from the 2 2Pi and 1 2Sigma+ states of the OCS+ ion.

Cl-loss and H-loss dissociations in low-lying electronic states of the CH3Cl+ ion studied using multiconfiguration second-order perturbation theory.

To examine the experimentally suggested scheme of the pathways for Cl- and H-loss dissociations of the CH(3)Cl(+) ion in the X(2)E (1(2)A', 1(2)A' '), A(2)A(1) (2(2)A'), and B(2)E (3(2)A', 2(2)A") states, the complete active space-self-consistent field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) calculations with an atomic natural orbital (ANO) basis were performed for the 1(2)A' (X(2)A'), 1(2)A", 2(2)A', and 2(2)A'" states. The potential energy curves describing dissociation from the four C(s) states were obtained on the basis of the CASSCF partial geometry optimization calculations at fixed C-Cl or C-H distance values, followed by the CASPT2 energy calculations. The electronic states of the CH3(+) and CH(2)Cl(+) ions produced by Cl-loss and H-loss dissociation, respectively, were carefully determined. Our calculations confirm the following experimental facts: Cl-loss dissociation occurs from the 1(2)A' (X(2)A'), 1(2)A", and 2(2)A' states (all leading to CH3(+) (X(1)A(1)') + Cl), and H-loss dissociation does not occur from 2(2)A'. The calculations indicate that H-loss dissociation occurs from the 1(2)A' and 1(2)A' ' states (leading to CH(2)Cl(+) (X(1)A(1)) + H and CH(2)Cl(+) (1(3)A") + H, respectively). The calculations also indicate that H-loss dissociation occurs (with a barrier) from the 2(2)A" state (leading to CH(2)Cl(+) (1(1)A") + H), supporting the observation of direct dissociation from the B state to CH(2)Cl(+) and that Cl-loss dissociation occurs from the 2(2)A" state (leading to CH3(+) (1(3)A") + Cl), not supporting the previously proposed Cl-loss dissociation of the B state via internal conversion of B to A. The predicted appearance potential values for CH3(+) (X(1)A(1)') and CH(2)Cl(+) (X(1)A(1)) are in good agreement with the experimental values.

F-loss and H-loss dissociations in low-lying electronic states of the CH3F+ ion studied using multiconfiguration second-order perturbation theory.

Complete active space self-consistent field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) calculations with an atomic natural orbital basis were performed for the 1(2)A'', 1(2)A', 2(2)A', 2(2)A'', and 3(2)A' (X2E, A2A1, and B2E) states of the CH3F+ ion. The 1(2)A'' state is predicted to be the ground state, and the C(s)-state energy levels are different from those of the CH3Cl+ ion. The 2(2)A' (A2A1) state is predicted to be repulsive, and the calculated adiabatic excitation energies for 2(2)A'' and 3(2)A' are very close to the experimental value for the B state. The CASPT2//CASSCF potential energy curves (PECs) were calculated for F-loss dissociation from the five C(s) states and H-loss dissociation from the 1(2)A'', 1(2)A', and 2(2)A'' states. The electronic states of the CH3+ and CH2F+ ions as the dissociation products were carefully determined by checking the energies and geometries of the asymptote products, and appearance potentials for the two ions in different states are predicted. The F-loss PEC calculations for CH3F+ indicate that F-loss dissociation occurs from the 1(2)A'', 1(2)A', and 2(2)A' states [all correlating with CH3+(X1A1')], which supports the experimental observations of direct dissociation from the X and A states, and that direct F-loss dissociation can occur from the two Jahn-Teller component states of B2E, 2(2)A'' and 3(2)A' [correlating with CH3+(1(3)A'') and CH3+(1(3)A'), respectively]. Some aspects of the 3(2)A' Cl-loss PEC of the CH3Cl+ ion are inferred on the basis of the calculation results for CH3F+. The H-loss PEC calculations for CH3F+ indicate that H-loss dissociation occurs from the 1(2)A'', 1(2)A', and 2(2)A'' states [correlating with CH2F+(1(3)A''), CH2F+(X1A1), and CH2F+(1(1)A''), respectively], which supports the observations of direct dissociation from the X and B states. As the 2(2)A' H-loss PEC of CH3Cl+, the 2(2)A' H-loss PEC of CH3F+ does not lead to H + CH2X+, but the PECs of the two ions represent different types of reactions.