Photodissociation and formation of an ion-pair in CH2 FCl (HCFC-31)

Although CH2 FCl (HCFC-31) recently became of great atmospheric importance, studies concerning its excited states are almost nonexistent. Several excited singlet states were studied (valence nσ* and Rydberg n3s, n3p, σ3s, and σ3p) through highly correlated multireference configuration interaction with singles and doubles, including extensivity correction. Comparison with the states of CH3 Cl indicates a strong influence of the F atom. Potential energy curves suggest formation of an electrostatically bound complex that relaxes to a hydrogen-bonded contact ion-pair (HBCIP) which can decay yielding CH2 F + Cl or to the ground state minimum of CH2 FCl. The HBCIP has a dipole moment of 9.57 D, a CI wavefunction described as 0.65ionic + 0.20biradical and it is strongly bonded by 4.72 eV. Its H bond has characteristics of moderate and strong H bonds. The simulated absorption spectrum confirms the nσ* assignment for the first and suggests the n3s + n3pσ assignment for the second band.

Surface hopping modeling of charge and energy transfer in active environments

An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.

Photochemistry of Monohydrated Chloromethane: Formation of Free and Hydrated Cl- and CH3+ Ions from a Solvent-Shared Semi-Ion-Pair

The effect of water molecule on the excited states of CH₃Cl(H₂O), as compared to those of the isolated chloromethane, has been studied at the multireference configuration interaction with singles and doubles (MR-CISD), including extensivity corrections. Eight new Rydberg states are due to the water molecule but the common states of both systems are not severely altered. Potential energy curves of 23 singlet states along the C-Cl coordinate have also been computed at the MR-CISD level. The dissociation energy of the C-Cl bond decreases from ∼0.4 to 0.5 eV due to the water molecule. As for CH₃Cl (de Medeiros, V. C., J. Am. Chem. Soc. 2016, 138, 272-280), a stable ion-pair has also been characterized. However, for CH₃Cl(H₂O), this ion-pair is better described as a solvent-shared semi-ion-pair, CH₃^+δ(H₂O)Cl^-δ. This species is connected with three ionic dissociation channels, with two being due to the water molecule. The presence of these new ionic channels, particularly the lowest energy one, [H₃C-O]⁺ + Cl^-, raises a very important question of atmospheric relevance: can the interaction of chloroalkanes with water decrease its deleterious effect on the ozone layer? Several potentially new competing dissociation channels are also studied. The latter results can help to set up the most important states to be included in nonadiabatic dynamic calculations to study how the yields of the ionic channels change due to the water molecule.

Can a gas phase contact ion pair containing a hydrocarbon carbocation be formed in the ground state?

So far, no conclusive evidence of a ground-state contact ion-pair containing a hydrocarbon carbocation has been given in the gas phase.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

Spin-Forbidden Branching in the Mechanism of the Intrinsic Haber-Weiss Reaction

The mechanism of the O₂^⋅− and H₂O₂ reaction (Haber–Weiss) under solvent‐free conditions has been characterized at the DFT and CCSD(T) level of theory to account for the ease of this reaction in the gas phase and the formation of two different set of products (Blanksby et al., Angew. Chem. Int. Ed. 2007, 46, 4948). The reaction is shown to proceed through an electron‐transfer process from the superoxide anion to hydrogen peroxide, along two pathways. While the O₃^⋅− + H₂O products are formed from a spin‐allowed reaction (on the doublet surface), the preferred products, O^⋅−(H₂O)+³O₂, are formed through a spin‐forbidden reaction as a result of a favorable crossing point between the doublet and quartet surface. Plausible reasons for the preference toward the latter set are given in terms of the characteristics of the minimum energy crossing point (MECP) and the stability of an intermediate formed (after the MECP) in the quartet surface. These unique results show that these two pathways are associated with a bifurcation, yielding spin‐dependent products.

Challenges encountered during development of Mn porphyrin-based, potent redox-active drug and superoxide dismutase mimic, MnTnBuOE-2-PyP5+, and its alkoxyalkyl analogues

We disclose here the studies that preceded and guided the preparation of the metal-based, redox-active therapeutic Mn(III) meso-tetrakis(N-n-butoxyethylpyridyl)porphyrin, MnTnBuOE-2-PyP⁵⁺ (BMX-001), which is currently in Phase I/II Clinical Trials at Duke University (USA) as a radioprotector of normal tissues in cancer patients. N-substituted pyridylporphyrins are ligands for Mn(III) complexes that are among the most potent superoxide dismutase mimics thus far synthesized. To advance their design, thereby improving their physical and chemical properties and bioavailability/toxicity profiles, we undertook a systematic study on placing oxygen atoms into N-alkylpyridyl chains via alkoxyalkylation reaction. For the first time we show here the unforeseen structural rearrangement that happens during the alkoxyalkylation reaction by the corresponding tosylates. Comprehensive experimental and computational approaches were employed to solve the rearrangement mechanism involved in quaternization of pyridyl nitrogens, which, instead of a single product, led to a variety of mixed N-alkoxyalkylated and N-alkylated pyridylporphyrins. The rearrangement mechanism involves the formation of an intermediate alkyl oxonium cation in a chain-length-dependent manner, which subsequently drives differential kinetics and thermodynamics of competing N-alkoxyalkylation versus in situ N-alkylation. The use of alkoxyalkyl tosylates, of different length of alkyl fragments adjacent to oxygen atom, allowed us to identify the set of alkyl fragments that would result in the synthesis of a single compound of high purity and excellent therapeutic potential.

Photochemistry of CH3Cl: Dissociation and CH···Cl Hydrogen Bond Formation

State-of-the-art electronic structure calculations (MR-CISD) are used to map five different dissociation channels of CH3Cl along the C-Cl coordinate: (i) CH3(X̃(2)A2″) + Cl((2)P), (ii) CH3(3s(2)A1') + Cl((2)P), (iii) CH3(+)((1)A1') + Cl(-)((1)S), (iv) CH3(3p(2)E') + Cl((2)P), and (v) CH3(3p(2)A2″) + Cl((2)P). By the first time these latter four dissociation channels, accessible upon VUV absorption, are described. The corresponding dissociation limits, obtained at the MR-CISD+Q level, are 3.70, 9.50, 10.08, 10.76, and 11.01 eV. The first channel can be accessed through nσ* and n3s states, while the second channel can be accessed through n(e)3s, n(e)3p(σ), and σ3s states. The third channel, corresponding to the CH3(+) + Cl(-) ion-pair, is accessed through n(e)3p(e) states. The fourth is accessed through n(e)3p(e), n(e)3p(σ), and σ3p(σ), while the fifth through σ3p(e) and σ(CH)σ* states. The population of the diverse channels is controlled by two geometrical spots, where intersections between multiple states allow a cascade of nonadiabatic events. The ion-pair dissociation occurs through formation of CH3(+)···Cl(-)and H2CH(+)···Cl(-) intermediate complexes bound by 3.69 and 4.65 eV. The enhanced stability of the H2CH(+)···Cl(-) complex is due to a CH···Cl hydrogen bond. A time-resolved spectroscopic setup is proposed to detect those complexes.

Assessment of density-functionals for describing the X(-) + CH3ONO2 gas-phase reactions with X = F, OH, CH2CN

The energetics of the ECO2, SN2@C and SN2@N channels of X(-) + CH3ONO2 (X = F, OH, CH2CN) gas-phase reactions were computed using the CCSD(T)/CBS method. This benchmark extends a previous study with X = OH [M. A. F. de Souza et al., J. Am. Chem. Soc., 2012, 134, 19004] and was used to ascertain the accuracy and robustness of nineteen density-functionals for describing these potential energy profiles (PEP) as well as the kinetic product distributions obtained from RRKM calculations. Assessments were based on the mean unsigned error (MUE), the mean signed error (MSE), the #best : #worst (BW) criterion and the statistical confidence interval (CI) for the MSE. In general, double-hybrid (DH) functionals perform better than the range-separated ones, and both are better than the global-hybrid functionals. Based on the MUE and CI criteria the B2GPPLYP, B2PLYP, M08-SO, BMK, ωB97X-D, CAM-B3LYP, M06, M08-HX, ωB97X and B97-K functionals show the best performance in the description of these PEPs. Within this set, the B2GPPLYP functional is the most accurate and robust. The RRKM results indicate that the DHs are the best for describing the selectivities of these reactions. Compared to CCSD(T), the B2PLYP method has a relative error of only ca. 1% for the selectivity and the accuracy to provide the correct conclusion concerning the nonstatistical behavior of these reactions.

Dynamic effects dictate the mechanism and selectivity of dehydration-rearrangement reactions of protonated alcohols [Me2 (R)CCH(OH2 )Me](+) (R=Me, Et, iPr) in the gas phase

The gas-phase dehydration-rearrangement (DR) reactions of protonated alcohols [Me2 (R)CCH(OH2 )Me](+) [R=Me (ME), Et (ET), and iPr (I-PR)] were studied by using static approaches (intrinsic reaction coordinate (IRC), Rice-Ramsperger-Kassel-Marcus theory) and dynamics (quasiclassical trajectory) simulations at the B3LYP/6-31G(d) level of theory. The concerted mechanism involves simultaneous water dissociation and alkyl migration, whereas in the stepwise reaction pathway the dehydration step leads to a secondary carbocation intermediate followed by alkyl migration. Internal rotation (IR) can change the relative position of the migrating alkyl group and the leaving group (water), so distinct products may be obtained: [Me(R)CCH(Me)Me⋅⋅⋅OH2 ](+) and [Me(Me)CCH(R)Me⋅⋅⋅OH2 ](+) . The static approach predicts that these reactions are concerted, with the selectivity towards these different products determined by the proportion of the conformers of the initial protonated alcohols. These selectivities are explained by the DR processes being much faster than IR. These results are in direct contradiction with the dynamics simulations, which indicate a predominantly stepwise mechanism and selectivities that depend on the alkyl groups and dynamics effects. Indeed, despite the lifetimes of the secondary carbocations being short (<0.5 ps), IR can take place and thus provide a rich selectivity. These different selectivities, particularly for ET and I-PR, are amenable to experimental observation and provide evidence for the minor role played by potential-energy surface and the relevance of the dynamics effects (non-IRC pathways, IR) in determining the reaction mechanisms and product distribution (selectivity).

Valence and Rydberg states of CH3Cl: a MR-CISD study

In this work ten singlet and nine triplet states are studied through multi-reference configuration interactions with singles and doubles (MR-CISD), including Davidson extensivity correction (MR-CISD+Q).

A theoretical study of hydrogen complexes of the XH-pi type between propyne and HF, HCL or HCN

The present manuscript reports a systematic investigation of the basis set dependence of some properties of hydrogen-bonded (pi type) complexes formed by propyne and a HX molecule, where X=F, Cl and CN. The calculations have been performed at Hartree-Fock, MP2 and B3LYP levels. Geometries, H-bond energies and vibrational have been considered. The more pronounced effects on the structural parameters of the isolated molecules, as a result of complexation, are verified on RCtriple bondC and HX bond lengths. As compared to double-zeta (6-31G**), triple-zeta (6-311G**) basis set leads to an increase of RCtriple bondC bond distance, at all three computational levels. In the case where diffuse functions are added to both hydrogen and 'heavy' atoms, the effect is more pronounced. The propyne-HX structural parameters are quite similar to the corresponding parameters of acetylene-HX complexes, at all levels. The largest difference is obtained for hydrogen bond distance, RH, with a smaller value for propyne-HX complex, indicating a stronger bond. Concerning the electronic properties, the results yield the following ordering for H-bond energies, DeltaE: propynecdots, three dots, centeredHF>propynecdots, three dots, centeredHCl>propynecdots, three dots, centeredHCN. It is also important to point out that the inclusion of BSSE and zero-point energies (ZPE) corrections cause significant changes on DeltaE. The smaller effect of ZPE is obtained for propynecdots, three dots, centeredHCN at HF/6-311++G** level, while the greatest difference is obtained at MP2/6-31G** level for propynecdots, three dots, centeredHF system. Concerning the IR vibrational it was obtained that larger shift can be associated with stronger hydrogen bonds. The more pronounced effect on the normal modes of the isolated molecule after the complexation is obtained for HX stretching frequency, which is shifted downward.

An ab initio study of the C2H2-HF, C2H(CH3)-HF and C2(CH3)2-HF hydrogen-bonded complexes

MP2/6-31++G** and B3LYP/6-31++G** ab initio molecular orbital calculations have been performed in order to obtain molecular geometries, binding energies and vibrational properties of the C2H2-HF, C2H(CH3)-HF and C2(CH3)2-HF H-bonded complexes. As expected, the more pronounced effects on the structural properties of the isolated molecules due to complexation was verified for the C[triple bond]C and H-F bond lengths, which are directly involved in the H-bond formation. These bond distances increased after complexation. BSSE uncorrected B3LYP binding energies are always lower than the corresponding MP2 values. However, the opposite trend has been verified after BSSE correction by the counterpoise method since it is much lower at B3LYP than at MP2 level. The binding energies for these complexes as well as for the HF acid submolecule modes (the HF stretching and vibrational frequency modes) showed an increasing hydrogen-bonding strength with increasing methyl substitution. The splitting in the HF in-plane and out-of-plane bending modes reflects the anisotropy in the hydrogen-bonding interaction with the pi system of the C[triple bond]C bond. The H-F stretching frequency is shifted downward after complexation and it increases with the methyl substitution. The IR intensities of the HF acid submolecule fundamentals are adequately interpreted through the atomic polar tensor of the hydrogen atom using the charge-charge flux-overlap model. The skeletal stretching modes of the Alkyne submolecule are decreased in the complex. The new vibrational modes arising from complexation show several interesting features.