CO2 Hydrogenation on Ru Single-Atom Catalyst Encapsulated in Silicalite: a DFT and Microkinetic Modeling Study.

The critical levels of CO2 emissions reached in the past decade have encouraged researchers into finding techniques to reduce the amount of anthropogenic CO2 expelled to the atmosphere. One possibility is to capture the produced CO2 from the source of emission or even from air (i.e., direct air capture) by porous materials (e.g., zeolites and MOFs). Among the different usages of captured CO2, its conversion into light fuels such as methane, methanol, and formic acid is essential for ensuring the long-awaited circular economy. In the last years, single-atom catalysts encapsulated in zeolites have been considered to this purpose since they exhibit a high selectivity and activity with the minimum expression of catalytic species. In this study, a detailed mechanism composed by 47 elementary reactions, 42 of them in both forward and reverse directions and 5 of them that correspond to the desorption of gas products just forwardly studied), has been proposed for catalytic CO2 hydrogenation over Ru SAC encapsulated in silicate (Ru1@S-1). Periodic density functional theory (DFT) calculations along with microkinetic modeling simulations at different temperatures and pressures were performed to evaluate the evolution of species over time. The analysis of the results shows that carbon monoxide is the main gas produced, followed by formic acid and formaldehyde. The rate analysis shows that CO(g) is formed mainly through direct dissociation of CO2 (i.e., redox mechanism), whereas COOH formation is assisted by OH. Moreover, the Campbell's degree of rate control analysis suggests that the determining steps for the formation of CO(g) and CH2O(g) gas species are their own desorption processes. The results obtained are in line with recent experimental and theoretical results showing that Ru1 SACs are highly selective to CO(g), whereas few atom clusters as Ru4 increase selectivity toward methane formation.

On the CO2 Harvesting from N2 Using Grazyne Membranes.

The separation of carbon dioxide (CO2) from nitrogen (N2) is at the core of any global warming remediation technology aimed at reducing the CO2 content in the atmosphere. Chemical membranes designed to differentially permeate both molecules have become quite appealing due to their simple use, although many membrane-based separations stand out as a promising solution for CO2 separation. These are environmentally friendly, with high active surface areas, compact design, easy to maintain and cost-effective, although the field is still growing due to the difficulties in the CO2/N2 separation. The present study poses grazynes, two-dimensional C-based materials with sp and sp2 C atoms, aligned along stripes, as suited membranes for the CO2/N2 separation. The combination of density functional theory (DFT) and molecular dynamics (MD) simulations allow tackling the energetics, kinetics, and dynamics of the membrane effectiveness of grazynes with engineered pores for such a separation in a holistic fashion. The explored grazynes are capable of physisorbing CO2 and N2, thus avoiding material poisoning by molecular decoration, while the diffusion of CO2 through the pores is found to be rapid, yet easier than that of N2, in the rate order of the s-1 in the 100-500 K temperature range. In particular, low-temperature CO2 separation even for CO2 contents below 0.5 % are found for [1],[2]{2}-grazyne when controlling the membrane exposure contact to the gas mixture, paving the way for exploring and using grazynes for air CO2 remediation.

Comprehensive Density Functional and Kinetic Monte Carlo Study of CO2 Hydrogenation on a Well-Defined Ni/CeO2 Model Catalyst: Role of Eley-Rideal Reactions.

A detailed multiscale study of the mechanism of CO2 hydrogenation on a well-defined Ni/CeO2 model catalyst is reported that couples periodic density functional theory (DFT) calculations with kinetic Monte Carlo (kMC) simulations. The study includes an analysis of the role of Eley-Rideal elementary steps for the water formation step, which are usually neglected on the overall picture of the mechanism, catalytic activity, and selectivity. The DFT calculations for the chosen model consisting of a Ni4 cluster supported on CeO2 (111) show large enough adsorption energies along with low energy barriers that suggest this catalyst to be a good option for high selective CO2 methanation. The kMC simulations results show a synergic effect between the two 3-fold hollow sites of the supported Ni4 cluster with some elementary reactions dominant in one site, while other reactions prefer the another, nearly equivalent site. This effect is even more evident for the simulations explicitly including Eley-Rideal steps. The kMC simulations reveal that CO is formed via the dissociative pathway of the reverse water-gas shift reaction, while methane is formed via a CO2 → CO → HCO → CH → CH2 → CH3 → CH4 mechanism. Overall, our results show the importance of including the Eley-Rideal reactions and point to small Ni clusters supported on the CeO2 (111) surface as potential good catalysts for high selective CO2 methanation under mild conditions, while very active and selective toward CO formation at higher temperatures.

Novel inhibitors that target bacterial virulence identified via HTS against intra-macrophage survival of Shigella flexneri.

Shigella flexneri is a facultative intracellular pathogen that causes shigellosis, a human diarrheal disease characterized by the destruction of the colonic epithelium. Novel antimicrobial compounds to treat infections are urgently needed due to the proliferation of bacterial antibiotic resistance and lack of new effective antimicrobials in the market. Our approach to find compounds that block the Shigella virulence pathway has three potential advantages: (i) resistance development should be minimized due to the lack of growth selection pressure, (ii) no resistance due to environmental antibiotic exposure should be developed since the virulence pathways are not activated outside of host infection, and (iii) the normal intestinal microbiota, which do not have the targeted virulence pathways, should be unharmed. We chose to utilize two phenotypic assays, inhibition of Shigella survival in macrophages and Shigella growth inhibition (minimum inhibitory concentration), to interrogate the 1.7 M compound screening collection subset of the GlaxoSmithKline drug discovery chemical library. A number of secondary assays on the hit compounds resulting from the primary screens were conducted, which, in combination with chemical, structural, and physical property analyses, narrowed the final hit list to 44 promising compounds for further drug discovery efforts. The rapid development of antibiotic resistance is a critical problem that has the potential of returning the world to a "pre-antibiotic" type of environment, where millions of people will die from previously treatable infections. One relatively newer approach to minimize the selection pressures for the development of resistance is to target virulence pathways. This is anticipated to eliminate any resistance selection pressure in environmental exposure to virulence-targeted antibiotics and will have the added benefit of not affecting the non-virulent microbiome. This paper describes the development and application of a simple, reproducible, and sensitive assay to interrogate an extensive chemical library in high-throughput screening format for activity against the survival of Shigella flexneri 2457T-nl in THP-1 macrophages. The ability to screen very large numbers of compounds in a reasonable time frame (~1.7 M compounds in ~8 months) distinguishes this assay as a powerful tool in further exploring new compounds with intracellular effect on S. flexneri or other pathogens with similar pathways of pathogenesis. The assay utilizes a luciferase reporter which is extremely rapid, simple, relatively inexpensive, and sensitive and possesses a broad linear range. The assay also utilized THP-1 cells that resemble primary monocytes and macrophages in morphology and differentiation properties. THP-1 cells have advantages over human primary monocytes or macrophages because they are highly plastic and their homogeneous genetic background minimizes the degree of variability in the cell phenotype (1). The intracellular and virulence-targeted selectivity of our methodology, determined via secondary screening, is an enormous advantage. Our main interest focuses on hits that are targeting virulence, and the most promising compounds with adequate physicochemical and drug metabolism and pharmacokinetic (DMPK) properties will be progressed to a suitable in vivo shigellosis model to evaluate the therapeutic potential of this approach. Additionally, compounds that act via a host-directed mechanism could be a promising source for further research given that it would allow a whole new, specific, and controlled approach to the treatment of diseases caused by some pathogenic bacteria.

Safe drugs with high potential to block malaria transmission revealed by a spleen-mimetic screening.

Malaria parasites like Plasmodium falciparum multiply in red blood cells (RBC), which are cleared from the bloodstream by the spleen when their deformability is altered. Drug-induced stiffening of Plasmodium falciparum-infected RBC should therefore induce their elimination from the bloodstream. Here, based on this original mechanical approach, we identify safe drugs with strong potential to block the malaria transmission. By screening 13 555 compounds with spleen-mimetic microfilters, we identified 82 that target circulating transmissible form of P. falciparum. NITD609, an orally administered PfATPase inhibitor with known effects on P. falciparum, killed and stiffened transmission stages in vitro at nanomolar concentrations. Short exposures to TD-6450, an orally-administered NS5A hepatitis C virus inhibitor, stiffened transmission parasite stages and killed asexual stages in vitro at high nanomolar concentrations. A Phase 1 study in humans with a primary safety outcome and a secondary pharmacokinetics outcome ( https://clinicaltrials.gov , ID: NCT02022306) showed no severe adverse events either with single or multiple doses. Pharmacokinetic modelling showed that these concentrations can be reached in the plasma of subjects receiving short courses of TD-6450. This physiologically relevant screen identified multiple mechanisms of action, and safe drugs with strong potential as malaria transmission-blocking agents which could be rapidly tested in clinical trials.

Effect of Nanoparticles on the Thermal Stability and Reaction Kinetics in Ionic Nanofluids.

Nowadays, the incorporation of nanoparticles into thermal fluids has become one of the most suitable strategies for developing high-performance fluids. An unconventional improvement of thermo-physical properties was observed with the addition of 1% wt. of nanoparticles in different types of fluids, such as molten salts, allowing for the design of more thermally efficient systems using nanofluids. Despite this, there is a lack of knowledge about the effect that nanoparticles produce on the thermal stability and the decomposition kinetics of the base fluid. The present study performs IR- and UV-vis spectroscopy along with thermogravimetric analysis (TGA) of pure nitrate and nitrate based nanofluids with the presence of SiO2 and Al2O3 nanoparticles (1% wt.). The results obtained support that nanoparticles accelerate the nitrate to nitrite decomposition at temperatures below 500 °C (up to 4%), thus confirming the catalytic role of nanoparticles in nanofluids.

The small-molecule SMARt751 reverses Mycobacterium tuberculosis resistance to ethionamide in acute and chronic mouse models of tuberculosis.

The sensitivity of Mycobacterium tuberculosis, the pathogen that causes tuberculosis (TB), to antibiotic prodrugs is dependent on the efficacy of the activation process that transforms the prodrugs into their active antibacterial moieties. Various oxidases of M. tuberculosis have the potential to activate the prodrug ethionamide. Here, we used medicinal chemistry coupled with a phenotypic assay to select the N-acylated 4-phenylpiperidine compound series. The lead compound, SMARt751, interacted with the transcriptional regulator VirS of M. tuberculosis, which regulates the mymA operon encoding a monooxygenase that activates ethionamide. SMARt751 boosted the efficacy of ethionamide in vitro and in mouse models of acute and chronic TB. SMARt751 also restored full efficacy of ethionamide in mice infected with M. tuberculosis strains carrying mutations in the ethA gene, which cause ethionamide resistance in the clinic. SMARt751 was shown to be safe in tests conducted in vitro and in vivo. A model extrapolating animal pharmacokinetic and pharmacodynamic parameters to humans predicted that as little as 25 mg of SMARt751 daily would allow a fourfold reduction in the dose of ethionamide administered while retaining the same efficacy and reducing side effects.

Acetylene-Mediated Electron Transport in Nanostructured Graphene and Hexagonal Boron Nitride.

The discovery of graphene has catalyzed the search for other 2D carbon allotropes, such as graphynes, graphdiynes, and 2D π-conjugated polymers, which have been theoretically predicted or experimentally synthesized during the past decade. These materials exhibit a conductive nature bound to their π-conjugated sp2 electronic system. Some cases include sp-hybridized moieties in their nanostructure, such as acetylenes in graphynes; however, these act merely as electronic couplers between the conducting π-orbitals of sp2 centers. Herein, via first-principles calculations and quantum transport simulations, we demonstrate the existence of an acetylene-meditated transport mechanism entirely hosted by sp-hybridized orbitals. For that we propose a series of nanostructured 2D materials featuring linear arrangements of closely packed acetylene units which function as sp-nanowires. Because of the very distinct nature of this unique transport mechanism, it appears to be highly complementary with π-conjugation, thus potentially becoming a key tool for future carbon nanoelectronics.

Quantum Dynamics of Nonadiabatic Renner-Teller Effects in Atom + Diatom Collisions.

We review the quantum nonadiabatic dynamics of atom + diatom collisions due to the Renner-Teller (RT) effect, i.e., to the Hamiltonian operators that contain the total spinless electronic angular momentum L̂. As is well-known, this rovibronic effect is large near collinear geometries when at least one of the interacting states is doubly degenerate. In general, this occurs in insertion reactions and at short-range, where the potential wells exhibit deep minima and support metastable complexes. Initial-state-resolved reaction probabilities, integral cross sections, and thermal rate constants are calculated via the real wavepacket method, solving the equation of motion with an approximated or with an exact spinless RT Hamiltonian. We present the dynamics of 10 single-channel or multichannel reactions showing how RT effects depend on the product channels and comparing with the Born-Oppenheimer (BO) approximation or coexisting conical-intersection (CI) interactions. RT effects not only can significantly modify the adiabatic dynamics or correct purely CI results, but also they can be very important in opening collision channels which are closed at the BO or CI level, as in electronic-quenching reactions. In the OH(A2Σ+) + Kr electronic quenching, where both nonadiabatic effects (CI and RT) coexist, they are in competition because CI dominates the reactivity but RT couplings reduce the large CI cross section and open a CI-forbidden evolution toward products, so that CI + RT quantum results are in good agreement with experimental or semiclassical findings. The different roles of these couplings are due to the unlike nuclear geometries where they are large: rather far from or near to linearity for CI or RT, respectively. The OH(A2Σ+) + Kr electronic quenching was investigated with the exact RT Hamiltonian, validating the approximated one, which was employed for all other collisions.

Rate-of-Kill (RoK) assays to triage large compound sets for Chagas disease drug discovery: Application to GSK Chagas Box.

Chagas disease (CD) is a human disease caused by Trypanosoma cruzi. Whilst endemic in Latin America, the disease is spread around the world due to migration flows, being estimated that 8 million people are infected worldwide and over 10,000 people die yearly of complications linked to CD. Current chemotherapeutics is restricted to only two drugs, i.e. benznidazole (BNZ) and nifurtimox (NIF), both being nitroaromatic compounds sharing mechanism of action and exerting suboptimal efficacy and serious adverse effects. Recent clinical trials conducted to reposition antifungal azoles have turned out disappointing due to poor efficacy outcomes despite their promising preclinical profile. This apparent lack of translation from bench models to the clinic raises the question of whether we are using the right in vitro tools for compound selection. We propose that speed of action and cidality, rather than potency, are properties that can differentiate those compounds with better prospect of success to show efficacy in animal models of CD. Here we investigate the use of in vitro assays looking at the kinetics of parasite kill as a valuable surrogate to tell apart slow- (i.e. azoles targeting CYP51) and fast-acting (i.e. nitroaromatic) compounds. Data analysis and experimental design have been optimised to make it amenable for high-throughput compound profiling. Automated data reduction of experimental kinetic points to tabulated curve descriptors in conjunction with PCA, k-means and hierarchical clustering provide drug discoverers with a roadmap to guide navigation from hit qualification of a screening campaign to compound optimisation programs and assessment of combo therapy potential. As an example, we have studied compounds belonging to the GSK Chagas Box stemmed from the HTS campaign run against the full GSK 1.8 million compounds collection [1].

Understanding the abnormal thermal behavior of nanofluids through infrared thermography and thermo-physical characterization.

Nanofluids (NFs) are colloidal suspensions of nanoparticles (NPs) within a base fluid. Unlike conventional mixtures, NFs exhibit dramatically enhanced properties, such as an abnormal increase in heat capacity at low concentration of NPs (e.g., Cp values 30% higher than the base material value). Understanding the thermo-physical behavior of NFs is essential for their application as thermal energy storage systems. In this study, we analyze a sodium nitrate ionic system containing 1 wt%, 3 wt% and 7 wt% of SiO2 NPs with different techniques like infrared thermography, infrared spectroscopy and differential scanning calorimetry (DSC) in order to shed light on the mechanism behind the increase of Cp. The themographies reveal the presence of a colder layer on top of the NF with 1 wt% of NPs whereas this layer does not appear at higher concentrations of NPs. The IR spectrum of this foamy top layer evidences the high amount of SiO2 bonds suggesting the clustering of the NPs into this layer linked by the nitrate ions. The linking is enhanced by the presence of hydroxyls in the NPs' surface (i.e., hydroxilated NPs) that once mixed in the NF suffer ionic exchange between OH- and NO3- species, leading to O2-Si-O-NO2 species at the interface where a thermal boundary resistance or Kapitza resistance appears (RT = 2.2 m2 K kW-1). Moreover, the presence of an exothermic reactive processes in the calorimetry of the mixture with 1 wt% of NPs evidences a reactive process (ionic exchange). These factors contribute to the heat capacity increase and thus, they explain the anomalous behavior of the heat capacity in nanofluids.

Non-adiabatic quantum dynamics of the electronic quenching OH(A2Σ+) + Kr.

We present the dynamics of the electronic quenching OH(A2Σ+) + Kr(1S) → OH(X2Π) + Kr(1S), with OH(A2Σ+) in the ground ro-vibrational state. This study relies on a new non-adiabatic quantum theory that uses three diabatic electronic states Σ+, Π', and Π'', coupled by one conical-intersection and nine Renner-Teller matrix elements, all of which are explicitly considered in the equation of the motion. The time-dependent mechanism and initial-state-resolved quenching probabilities, integral cross sections, thermal rate constants, and thermally-averaged cross sections are calculated via the real wavepacket method. The results point out a competition among three non-adiabatic pathways: Σ+ ↔ Π', Σ+ ↔ Π'', and Π' ↔ Π''. In particular, the conical-intersection effects Σ+-Π' are more important than the Renner-Teller couplings Σ+-Π', Σ+-Π'', and Π'-Π''. Therefore, Π' is the preferred product channel. The quenching occurs via an indirect insertion mechanism, opening many collision complexes, and the probabilities thus present many oscillations. Some resonances are still observable in the cross sections, which are in good agreement with previous experimental and quasi-classical findings. We also discuss the validity of more approximate quantum methods.

Gametocytes from K13 Propeller Mutant Plasmodium falciparum Clinical Isolates Demonstrate Reduced Susceptibility to Dihydroartemisinin in the Male Gamete Exflagellation Inhibition Assay.

Mutations in the kelch propeller domain (K13 propeller) of Plasmodium falciparum parasites from Southeast Asia are associated with reduced susceptibility to artemisinin. We exposed in vitro-cultured stage V gametocytes from Cambodian K13 propeller mutant parasites to dihydroartemisinin and evaluated the inhibition of male gamete formation in an in vitro exflagellation inhibition assay (EIA). Gametocytes with the R539T and C580Y K13 propeller alleles were less susceptible to dihydroartemisinin and had significantly higher 50% inhibitory concentrations (IC50s) than did gametocytes with wild-type alleles.

A lexicon based method to search for extreme opinions.

Studies in sentiment analysis and opinion mining have been focused on many aspects related to opinions, namely polarity classification by making use of positive, negative or neutral values. However, most studies have overlooked the identification of extreme opinions (most negative and most positive opinions) in spite of their vast significance in many applications. We use an unsupervised approach to search for extreme opinions, which is based on the automatic construction of a new lexicon containing the most negative and most positive words.

Nonadiabatic Renner-Teller quantum dynamics of OH(X2Π) + H+ reactive collisions.

Following previous studies on the O(3P) + H2+(X2Σg+) collisions, we present the nonadiabatic quantum dynamics of the reactions OH(X2Π) + H'+ → OH'(X2Π) + H+, exchange (e), → OH+(X3Σ-) + H'(2S), quenching (q), and → OH'+ (X3Σ-) + H(2S), exchange-quenching (eq). The reactants and products correlate via the ground X[combining tilde]2A'' and first excited Ã2A' electronic states of OH2+, which are the degenerate components of linear 2Π species. Therefore, they are strongly perturbed by nonadiabatic Renner-Teller (RT) effects, opening the (q) and (eq) channels that are closed in the Born-Oppenheimer approximation. Using accurate potential energy surfaces (PESs) and RT matrix elements, initial-state-resolved reaction probabilities, real-time dynamics, cross sections, and rate constants of the product channels are obtained through the time-dependent real wavepacket (WP) method and full coupled-channel calculations. Owing to the nonadiabatic couplings, the WP jumps from the excited Ã2A' surface to the X[combining tilde]2A'' ground PES, avoiding any barrier, opening the quenching channels, and giving many collision complexes into the deep minima of both PESs, as it is clearly shown by the oscillations of the reaction probabilities and by the time-dependent WP dynamics. All the results show that the nonadiabatic-RT channels (q) and (eq) are highly reactive, much more than the adiabatic one (e), pointing out large RT effects. The reactivity of the quenching channels is similar, accounting for 97% of the overall reactivity. In fact, the maximum values of the (q) and (eq) cross sections σq and σeq are equal to 31.6 Å2, whereas the maximum σe value equals 1.34 Å2, and the maximum values of the rate constants kq, keq, and ke are 2.07 × 10-10, 2.45 × 10-10, and 0.23 × 10-10 cm3 s-1. Some calculations show that the centrifugal-sudden and the truncated coupled-channel approximations cannot be employed for the (q) channel. After a sharp increase at the threshold, σq and σeq decrease at larger collision energies while σe and the rate constants depend slightly on the collision energy and temperature, respectively. These findings are consistent with the barrierless nature of both PESs and the exoergicity of the quenching products, with the small role played by the centrifugal and RT barriers in the reactant channel, and with the large RT couplings in the OH2+ intermediates. Finally, we contrast the present results with those for the opposite reactions O + H2+ and for the nonadiabatic-induced quenchings NH + H' and OH + H'.

Dynamics of the O + H2+ → OH+ + H, OH + H+ proton and hydrogen atom transfer reactions on the two lowest potential energy surfaces.

The dynamics of the title reaction was studied using mainly the quasiclassical trajectory (QCT) method on the ground 12A'' (OH+ channel) and first excited 12A' (OH channel) potential energy surfaces (PESs) employing ab initio analytical representations of the PESs developed by us. Both PESs correspond to exoergic reactions, are barrierless and present a deep minimum along the minimum energy path (MEP). Some extra calculations (cross sections) were also performed with the time dependent quantum real wave packet method at the centrifugal sudden level (RWP-CS method). A broad set of properties as a function of collision energy (Ecol ≤ 0.5 eV) was considered using the QCT method: cross sections, average fractions of energy, product rovibrational distributions, two- and three-vector properties, and the microscopic mechanisms analyzing their influence on the dynamics. The proton transfer channel dominates the reactivity of the system and significant differences between the two reaction channels are found for the vibrational distributions and microscopic mechanisms. The results were interpreted according to the properties of the ground and excited PESs. Moreover, the QCT and RWP-CS cross sections are in rather good agreement for both reaction channels. We hope that this study will encourage the experimentalists to investigate the dynamics of this interesting but scarcely studied system, whose two lowest PESs include the ground and first excited electronic states of the H2O+ cation.

Extraction and analysis of signatures from the Gene Expression Omnibus by the crowd.

Gene expression data are accumulating exponentially in public repositories. Reanalysis and integration of themed collections from these studies may provide new insights, but requires further human curation. Here we report a crowdsourcing project to annotate and reanalyse a large number of gene expression profiles from Gene Expression Omnibus (GEO). Through a massive open online course on Coursera, over 70 participants from over 25 countries identify and annotate 2,460 single-gene perturbation signatures, 839 disease versus normal signatures, and 906 drug perturbation signatures. All these signatures are unique and are manually validated for quality. Global analysis of these signatures confirms known associations and identifies novel associations between genes, diseases and drugs. The manually curated signatures are used as a training set to develop classifiers for extracting similar signatures from the entire GEO repository. We develop a web portal to serve these signatures for query, download and visualization.

Unexpectedly large impact of van der Waals interactions on the description of heterogeneously catalyzed reactions: the water gas shift reaction on Cu(321) as a case example.

The molecular mechanisms of the water gas shift reaction on Cu(321) have been chosen to investigate the effect of dispersion terms on the description of the energy profile and reaction rates. The present study based on periodic DFT calculations shows that including dispersion terms does not change the qualitative picture of the overall reaction, maintaining the rate determining step and the predominant route. However, the effect of dispersion is different for different adsorbates - reactants, intermediates or products - with a clear net effect and with no compensation of errors. Thus, in the OH + OH → H2O + O process the dispersion effects imply up to three orders of magnitude in the calculated reaction rates; the formation of carboxyl is highly disfavoured when dispersion terms are explicitly included and finally, the reaction rate for CO2 production (at 463 K) through cis-COOH dissociation is enhanced by three orders of magnitude by including dispersion terms in the calculation of the energy barrier. Consequently, the inclusion of dispersion terms largely affects the overall potential energy profile and produces tremendous changes in the predicted reaction rates. Therefore, dispersion terms must be included when aiming at obtaining information from macroscopic simulations employing for instance microkinetic or kinetic Monte Carlo approaches, where these effects should be clearly shown.

Born-Oppenheimer and Renner-Teller Quantum Dynamics of CH(X(2)Π) + D((2)S) Reactions on Three CHD Potential Surfaces.

The quantum dynamics of three CH(X(2)Π) + D((2)S) reactions is studied by means of the coupled-channel time-dependent real-wavepacket (WP) and flux methods at collision energy Ecol ≤ 0.6 eV and on three potential energy surfaces (PESs): the Born-Oppenheimer (BO) ground PES X̃(3)A″ and the excited ones ã(1)A' and b̃(1)A″, coupled by nonadiabatic (NA) Renner-Teller (RT) effects. This three-state model is suitable for obtaining initial-state-resolved observables, is based on a complete analysis of the correlation diagram of the lowest electronic states of the CHD intermediate and of their NA interactions, and neglects the smaller coupling effects due to the asymptotic electronic angular momenta that become important in state-to-state dynamics. WPs are propagated on each PES at total angular momentum values J ≤ 70, with CH in the two lowest vibrational states v0 and in the ground rotational state j0 = 1. Reaction probabilities are obtained for three possible final products (f): (dP) CH decay and C((3)P) + HD(X(1)Σ(+)) formation that occurs on the uncoupled ground PES, (dD) CH decay and C((1)D) + HD(X(1)Σ(+)) formation that depends on the RT-coupled singlet species, and (ex) exchange to CD(X(2)Π) + H((2)S) available adiabatically from the X̃(3)A″ PES and nonadiabatically from ã(1)A' and b̃(1)A″. Observable cross sections σf,v0j0 and rate constants kf,v0j0 in the temperature range T = 100-2000 K are obtained for (dP), (dD), and (ex) channels. Comparing BO with RT probabilities, we show that NA effects are important at high J values for the (ex) channel at v0 = 1. Real time mechanisms on the three PESs show that RT couplings are opened after some time and clearly point out the formation of the product channels. Both cross sections and rate constants present the same sequence, for example σex,11 > σdP,01 ∼ σex,01 > σdP,11 ≫ σdD,11 ≫ σdD,01, and the CH vibrational excitation enhances the total removal CH+D reactivity by a factor of ∼1.7, mainly due to the increase of the (ex) channel contribution from ∼47% at v0 = 0 to ∼76% at v0 = 1. This fact implies a considerable vibrational enhancement of combustion processes at high temperature. In agreement with the probability results, the ã(1)A'/b̃(1)A″ RT coupling increases both σex,11 and kex,11 up to ∼30%. Moreover, including the three PESs in the dynamics simulation of CH+D increase by far the (ex)/(dP) branching ratio with respect to the CH + H' reaction. Thus, at room temperature, kdP,01 changes from 10.8 × 10(-11) to 3.4 × 10(-11) cm(3) s(-1) substituting H atom by D.

Born-Oppenheimer and Renner-Teller coupled-channel quantum reaction dynamics of O((3)P) + H2(+)(X(2)Σg(+)) collisions.

We present Born-Oppenheimer (BO) and Renner-Teller (RT) time dependent quantum dynamics studies of the reactions O((3)P) + H2(+)(X(2)Σg(+)) → OH(+)(X(3)Σ(-)) + H((2)S) and OH(X(2)Π) + H(+). We consider the OH2(+) X[combining tilde](2)A'' and Ã(2)A' electronic states that correlate with a linear (2)Π species. The electronic angular momenta operators L[combining circumflex] and L[combining circumflex](2) are considered in nonadiabatic coupled-channel calculations, where the associated RT effects are due to diagonal V(RT) potentials that add up to the PESs and to off-diagonal C(RT) couplings between the potential energy surfaces (PESs). Initial-state-resolved reaction probabilities PI, integral cross sections σI, and rate constants kI are obtained using recent ab initio PESs and couplings and the real wavepacket formalism. Because the PESs are strongly attractive, PI have no threshold energy and are large, σI decrease with collision energy, and kI depend little on the temperature. The X[combining tilde](2)A'' PES is up to three times more reactive than the Ã(2)A' PES and H2(+) rotational effects (j0 = 0, 1) are negligible. The diagonal V(RT) potentials are strongly repulsive at the collinearity and nearly halve all low-energy observables with respect to the BO ones. The off-diagonal C(RT) couplings are important at low partial waves, where they mix the X[combining tilde](2)A'' and Ã(2)A' states up to ∼20%. However, V(RT) effects predominate over the C(RT) ones that change at most by ∼19% the BO values of σI and kI. The reaction O((3)P) + H2(+)(X(2)Σg(+)) → OH(+)(X(3)Σ(-)) + H((2)S) is probably one of the most reactive atom + diatom collisions because its RT rate constant at room temperature is equal to 2.26 × 10(-10) cm(3) s(-1). Within the BO approximation, the present results agree rather well with recent quasiclassical and centrifugal-sudden data using the same PESs.