def2-TZVPP Level Research Articles

Elucidation of the influence of trinitro-diazinotriazine isomerism on the energetic properties and stability: Insights from DFT approach

The exploration of isomeric energetic molecules with high nitrogen content and the examination of their structure–property relationships remains compelling areas of interest within the field of energetic materials. Energetic isomers comprising the same molecular composition but distinct molecular structures potentially lead to notable variations in performance and properties. This study focused on designing and evaluating different forms of diazinotriazine using two molecular compositions: C5H3N5 and C5N8O6, aiming to generate various isomeric structures and to guide the rational design of new molecules. The electronic configurations, energetic properties, and sensitivity prediction of the optimized structures were investigated using the M06-2X/def2-TZVPP level of theory. The trinitro-diazinotriazine isomers exhibit high heat of formation (>486 kJ/mol), high detonation properties (D > 8.7 km/s, P > 34 GPa), high density (>1.81 g/cm3), and better thermal stability (BDE > 230 kJ/mol) making them promising explosive contenders. Among the designed diazinotriazine backbones, the pyridazine-1,2,3-triazine framework is favourable for the heat of formation and detonation properties. The heats of formation, detonation velocity, pressure, and heats of detonation of A4 and B4 outperform the other trinitro-diazinotriazine isomers.

Nature of Charge Transfer Effects in Complexes of Dopamine Derivatives Adsorbed on Graphene-Type Nanostructures

Continuing the investigation started for dopamine (DA) and dopamine-o-quinone (DoQ) (see, the light absorption and charge transfer properties of the dopamine zwitterion (called dopamine-semiquinone or DsQ) adsorbed on the graphene nanoparticle surface is investigated using the ground state and linear-response time-dependent density functional theories, considering the ωB97X-D3BJ/def2-TZVPP level of theory. In terms of the strength of molecular adsorption on the surface, the DsQ form has 50% higher binding energy than that found in our previous work for the DA or DoQ cases (−20.24 kcal/mol vs. −30.41 kcal/mol). The results obtained for electronically excited states and UV-Vis absorption spectra show that the photochemical behavior of DsQ is more similar to DA than that observed for DoQ. Of the three systems analyzed, the DsQ-based complex shows the most active charge transfer (CT) phenomenon, both in terms of the number of CT-like states and the amount of charge transferred. Of the first thirty electronically excited states computed for the DsQ case, eleven are purely of the CT type, and nine are mixed CT and localized (or Frenkel) excitations. By varying the adsorption distance between the molecule and the surface vertically, the amount of charge transfer obtained for DA decreases significantly as the distance increases: for DoQ it remains stable, for DsQ there are states for which little change is observed, and for others, there is a significant change. Furthermore, the mechanistic compilation of the electron orbital diagrams of the individual components cannot describe in detail the nature of the excitations inside the complex.

Vibrational spectroscopy and structural characterization of neutral and cationic hydroxyacetone

In this study, we utilized the infrared (IR)–vacuum-ultraviolet (VUV) non-resonant ionization and fragmentation detected IR spectroscopy (NRIFD-IR) technique to explore the IR spectra (2700–7200 cm-1) and the structures of both neutral and cationic hydroxyacetone. By comparing the experimentally obtained IR spectra with the anharmonic frequencies calculated at the B2PLYP-D3/cc-pVTZ and B3LYP-D3(BJ)/def2-TZVPP level of theory, we assigned the structures of hydroxyacetone and its cation. Notably, we observed a discrepancy in the neutral spectrum within the CH stretch region, depending on whether the parent or fragment channel was monitored. This observation suggests that the near-threshold dissociative ionization process exhibits mode-selectivity, implying that certain vibrational modes play a crucial role in either facilitating or inhibiting the dissociation event.

Characterization and conformational analysis of a novel hydrazonoyltetrazole: Crystallographic and theoretical investigations of an unexpected reaction product

A new hydrazonoyltetrazole was isolated as a byproduct from a previously reported reaction, and characterized by 1H and 13C NMR spectroscopy, infrared spectroscopy, single-crystal X-ray diffraction, and high-resolution mass spectroscopy. Intermolecular interactions in the crystal state were examined by Hirshfeld surface analysis (HSA), 2D fingerprint plots, and noncovalent interaction analysis (NCI). Geometry optimization in the gaseous phase was executed at the ωB97X-D4/def2-TZVPP level of theory, and subsequent maps of molecular electrostatic potential (MEP) were generated and analyzed. The calculated MEP was in good agreement with the HSA and NCI results. Conformational analysis was performed using a CREST/CENSO protocol. Two low-energy solution conformers were obtained. To confirm the presence of the second conformer, 1H NMR spectra were calculated, the transition state connecting the two conformers was found, and their energies were obtained with the DLPNOCCSD(T) scheme. A barrier of only 4 kcal mol-1 was found between them, and their MEP maps were also created. The dispersion energies between fragments of the three structures were calculated using the local energy decomposition (LED) scheme.

Using Hybrid PDI-Fe3O4 Nanoparticles for Capturing Aliphatic Alcohols: Halogen Bonding vs. Lone Pair-π Interactions.

In this study, Fe3O4 nanoparticles (FeNPs) decorated with halogenated perylene diimides (PDIs) have been used for capturing VOCs (volatile organic compounds) through noncovalent binding. Concretely, we have used tetrachlorinated/brominated PDIs as well as a nonhalogenated PDI as a reference system. On the other hand, methanol, ethanol, propanol, and butanol were used as VOCs. Experimental studies along with theoretical calculations (the BP86-D3/def2-TZVPP level of theory) pointed to two possible and likely competitive binding modes (lone pair-π through the π-acidic surface of the PDI and a halogen bond via the σ-holes at the Cl/Br atoms). More in detail, thermal desorption (TD) experiments showed an increase in the VOC retention capacity upon increasing the length of the alkyl chain, suggesting a preference for the interaction with the PDI aromatic surface. In addition, the tetrachlorinated derivative showed larger VOC retention times compared to the tetrabrominated analog. These results were complemented by several state-of-the-art computational tools, such as the electrostatic surface potential analysis, the Quantum Theory of Atoms in Molecules (QTAIM), as well as the noncovalent interaction plot (NCIplot) visual index, which were helpful to rationalize the role of each interaction in the VOC···PDI recognition phenomena.

Bistriazolotriazole-tetramine: commendable energetic moiety and cation.

Various 7H,7'H-[6,6'-bi[1,2,4]triazolo[4,3-b][1,2,4]triazole]-3,3',7,7'-tetramine (A) based nitrogen-rich energetic salts were designed and their properties explored. All energetic salts possess relatively high nitrogen content (> 48%), positive heats of formation (> 429kJ/mol) and stability owing to a significant contribution from fused backbone. The cationic component shows a very high heat of formation (2516kJ/mol); therefore, it is highly suitable for enthalpy enhancement in new energetic salts. The cation was paired with the energetic anions nitrate (NO3-), perchlorate (ClO4-), dinitromethanide (CH(NO2)2-), trinitromethanide (C(NO2)3-), nitroamide (NHNO2-), and dinitroamide (N(NO2)2-) to improve oxygen balance and detonation performance. Designed salts show moderate detonation velocities (7.9-8.7 km/s) and pressures (23.8 - 33.1 GPa). The distribution of frontier molecular orbitals, molecular electrostatic surface potentials, QTAIM topological properties, and noncovalent interactions of designed salts were simulated to understand the electronic structures, charge distribution in molecules, hydrogen bonding, and other nonbond interactions. The predicted safety factor (SF) and impact sensitivity (H50) of designed salts suggest their insensitivity to mechanical stimuli. This work explored the 7H,7'H-[6,6'-bi[1,2,4]triazolo[4,3-b][1,2,4]triazole]-3,3',7,7'-tetramine as a suitable cationic component which could be promising and serve exemplarily in energetic materials. The optimization and energy calculations of all the designed compounds were carried out at the B3LYP/6-311 + + G(d,p) and M06-2X/def2-TZVPP levels, utilizing the Gaussian software package. The molecular surface electrostatic potential, quantum theory of atoms in molecules (QTAIM), reduced density gradient (RDG), and noncovalent interaction (NCI) analysis were performed by employing Multiwfn software. The EXPLO5 (v 7.01) thermochemical code and PILEM web application were used to predict the detonation properties.

Hydrogen bonds of a water molecule in the second coordination sphere of amino acid metal complexes: influence of amino acid types on different complex geometries

Quantum chemical calculations at the M06L-GD3/def2-TZVPP level were done to investigate hydrogen bonds between amino acid metal complexes and a free water molecule. Octahedral nickel(II) and square planar palladium(II) complexes of glycine, cysteine, phenylalanine, and serine with different charges of metal complexes (+1, 0, and −1) were investigated. The following hydrogen bonds were considered: NH/O (amino acid is a H-donor), O1/HO (coordinated O1 oxygen from amino acid is a H-acceptor), and O2/HO (non-coordinated O2 oxygen from amino acid is a H-acceptor). Amino acid type has a small influence on interaction energies, both for octahedral nickel(II) and square planar palladium(II) complexes. The influence is the largest for NH/O and the smallest for O2/HO hydrogen bonds. For NH/O interaction, palladium(II) complexes showed stronger hydrogen bonds than nickel(II), up to −11.8 kcal mol−1 for singly positively charged complexes. Nickel(II) complexes demonstrated higher O1/HO hydrogen bond strength than palladium(II) with interaction energies up to −8.9 kcal mol−1 for singly negative complexes. With up to −9.0 kcal mol−1 interaction energy for singly negative complexes, O2/HO interactions were also stronger for nickel(II) complexes than palladium(II).

Analysis of the Unusual Chemical Bonds and Dipole Moments of AeF- (Ae=Be-Ba): A Lesson in Covalent Bonding.

Quantum chemical calculations of the anions AeF- (Ae=Be-Ba) have been carried out using ab initio methods at the CCSD(T)/def2-TZVPP level and density functional theory employing BP86 with various basis sets. The detailed bonding analyses using different charge- and energy partitioning methods show that the molecules possess three distinctively different dative bonds in the lighter species with Ae=Be, Mg and four dative bonds when Ae=Ca, Sr, Ba. The occupied 2p atomic orbitals (AOs) and to a lesser degree the occupied 2s AO of F- donate electronic charge into the vacant spx(σ) and p(π) orbitals of Be and Mg which leads to a triple bond Ae F-. The heavier Ae atoms Ca, Sr, Ba use their vacant (n-1)d AOs as acceptor orbitals which enables them to form a second σ donor bond with F- that leads to quadruply bonded Ae F- (Ae=Ca-Ba). The presentation of molecular orbitals or charge distribution using only one isodensity value may give misleading information about the overall nature of the orbital or charge distribution. Better insights are given by contour line diagrams. The ELF calculations provide monosynaptic and disynaptic basins of AeF- which nicely agree with the analysis of the occupied molecular orbitals and with the charge density difference maps. A particular feature of the covalent bonds in AeF- concerns the inductive interaction of F- with the soft valence electrons in the (n)s valence orbitals of Ae. The polarization of the (n)s2 electrons induces a (n)spx hybridized lone-pair orbital at atom Ae, which yields a large dipole moment with the negative end at Ae. The concomitant formation of a vacant (n)spx AO of atom Ae, which overlaps with the occupied 2p(σ) AO of F-, leads to a strong covalent σ bond.

The strongest dative bond in main-group compounds. Theoretical study of OAeF- (Ae = Be-Ba).

Quantum chemical calculations of the anions OAeF- (Ae = Be-Ba) have been carried out using ab initio methods at the CCSD(T)/def2-TZVPP level and density functional theory employing BP86 with various basis sets. The equilibrium structures have linear geometries for Ae = Be and Mg but they are strongly bent for Ae = Sr and Ba while the calcium species has a quasi-linear structure with a very low bending potential. The calculated bond dissociation energies suggest a record-high BDE of De = 144.08 kcal mol-1 for OBeF- at the CCSD(T)/def2-TZVPP level, which is the strongest BDE for a dative bond that has been found so far. The BDE of the heavier homologues have a continuously decreasing order for Ae with Be > Mg (113.01 kcal mol-1) > Ca (84.06 kcal mol-1) > Sr (72.06 kcal mol-1) > Ba (60.00 kcal mol-1). The calculation of the charge distribution reveals a significant charge donation OAe ← F- with a declining sequence for the heavier atoms Ae. The oxygen atom in OAeF- carries always a higher partial charge than the fluorine atom, which contradicts the standard electronegativities of the atoms. The surprising partial charges are explained with the bonding situation of the atoms in the actual electronic structure. The bonding analysis of the OAe-F- bonds using the EDA-NOCV method shows that the bonds have much more electrostatic character than the Ae-F- bonds in the diatomic anions. This finding is supported by the results of the LED partitioning approach. The dative interactions have three major and one minor component. The assignment of a quadruple bond for the heavier species with Ae = Ca, Sr, Ba is not reasonable. The driving force for the bent geometries is the accumulation of electronic charge in the lone-pair region at the Ae atoms, which enhances the electrostatic attraction with the other atoms. An adequate description of the bonding situation is given by the formula O--Ae+ ← F-.

Predicting valence tautomerism in diverse cobalt-dioxolene complexes: elucidation of the role of ligands and solvent.

The ability of molecular switches to reversibly interconvert between different forms promises potential applications at the scale of single molecules up to bulk materials. One type of molecular switch comprises cobalt-dioxolene compounds that exhibit thermally-induced valence tautomerism (VT) interconversions between low spin Co(iii)-catecholate (LS-CoIII-cat) and high spin Co(ii)-semiquinonate (HS-CoII-sq) forms. Two families of these compounds have been investigated for decades but have generally been considered separately: neutral [Co(diox)(sq)(N2L)] and cationic [Co(diox)(N4L)]+ complexes (diox = generic dioxolene, N2L/N4L = bidentate/tetradentate N-donor ancillary ligand). Computational identification of promising new candidate compounds prior to experimental exploration is beneficial for environmental and cost considerations but requires a thorough understanding of the underlying thermochemical parameters that influence the switching. Herein, we report a robust approach for the analysis of both cobalt-dioxolene families, which involved a quantitative density functional theory-based study benchmarked with reliable quasi-experimental references. The best-performing M06L-D4/def2-TZVPP level of theory has subsequently been verified by the synthesis and experimental investigation of three new complexes, two of which exhibit thermally-induced VT, while the third remains in the LS-CoIII-cat form across all temperatures, in agreement with prediction. Valence tautomerism in solution is markedly solvent-dependent, but the origin of this has not been definitively established. We have extended our computational approach to elucidate the correlation of VT transition temperature with solvent stabilisation energy and change in dipole moment. This new understanding may inform the development of VT compounds for applications in soft materials including films, gels, and polymers.

Computational and crystallographic study of hydrogen bonds in the second coordination sphere of chelated amino acids with a free water molecule: Influence of complex charge and metal ion

Hydrogen bonds of glycine complexes were calculated using quantum chemistry calculations at M06L-GD3/def2-TZVPP level and by analyzing the crystal structures from the Cambridge Structural Database (CSD). One hydrogen bond where amino acid plays the role of the H-donor (NH/O), and two where it plays the role of the H-acceptor (O1/HO, O1 is a coordinated oxygen atom, and, O2/HO, O2 is a non-coordinated oxygen atom) were investigated. The calculations were done on octahedral nickel(II), square pyramidal copper(II), square planar copper(II), palladium(II), and platinum(II) glycine complexes with different charges adjusted using water(s) and/or chlorine ion(s) as the remaining ligands. For NH/O hydrogen bond, interaction energies of neutral complexes are the weakest, from −5.2 to −7.2 kcal/mol for neutral, stronger for singly positive, from −8.3 to −12.1 kcal/mol, and the strongest for doubly positive complex, −16.9 kcal/mol. For O1/HO and O2/HO interactions, neutral complexes have weaker interaction energies (from −2.2 to −5.1 kcal/mol for O1/HO, and from −3.7 to −5.0 kcal/mol for O2/HO), than for singly negative complexes (from −6.9 to −8.2 kcal/mol for O1/HO, and from −8.0 to −9.0 kcal/mol for O2/HO). Additionally to the complex charge, metal oxidation number, coordination number, and metal atomic number also influence the hydrogen bond strength, however, the influence is smaller.

Vibrational spectroscopy of neutral and cationic 2-Methoxyethanol

Hydrogen bonding interactions, such as the O-H···O interaction, play a crucial role in stabilizing conformations of both organic and biological molecules. In this study, we employed the infrared (IR)–vacuum-ultraviolet (VUV) non-resonant ionization detected IR spectroscopy (NRID-IR) method to investigate the molecular structure of neutral and cationic 2-methoxyethanol (CH3OCH2CH2OH, 2-ME). The stable structures and anharmonic IR spectra of neutral and cationic 2-ME were calculated using density functional theory (DFT) at the B3LYP-D3(BJ)/def2-TZVPP level. Our results revealed that the two most stable conformers of neutral 2-ME exhibit a weak O-H···O intramolecular hydrogen bond, while the cationic 2-ME lacks the O-H···O interaction but contains a C-H···O intramolecular hydrogen bond. A comparison of the experimental and theoretical IR absorption spectra reveals that the most stable conformer (gauche-(anti-gauche)-trans) is the primary contributor to the observed IR spectra of neutral 2-ME in the 2700–7250 cm−1 range. Similarly, the second stable conformer is found to be the dominant contributor to the observed IR spectra of cationic 2-ME in the range of 2700–7100 cm−1. Additionally, the C-H fundamental stretching mode of cationic 2-ME is found to be blue-shifted by approximately 100 cm−1 compared to that of neutral 2-ME. Furthermore, the natural bond orbital (NBO) analysis suggests that in neutral 2-ME, the n(O)→σ*(CαH) interactions result in significant negative hyperconjugation, leading to the weakening of CαH bonds and a subsequent reduction in vibrational energy.

Investigation on the correlation of stability, reactivity, and structural properties of C/C2-doped neutral and charged Aln (n = 2–7) clusters

The equilibrium geometries, electronic structures, reactivity, and stability of C/C2-doped Aln0,± (n = 2–7) clusters were investigated computationally at B3LYP-D4/def2-TZVPP and CCSD(T)/def2-TZVPP levels of theory. The lowest energy structures of AlnC1,20,± clusters were the carbon atoms inserted inside the aluminium clusters. The stability of the neutral and charged clusters were investigated under various reaction conditions, such as dissociation of the clusters via Al0,+,C0,-,C20,--elimination and electron elimination from neutral (ionization) and anionic species (electron detachment). The dominant dissociation reaction channels are the Al-elimination from the neutral and anionic clusters, Al+-elimination from cationic clusters, rather than the C/C2-elimintion reaction. On the other hand, C-elimination are were more favorable compared to C2-elimination due to the stronger CC and Al-C bonds. In addition, binding energy, second-order difference of energy, chemical hardness, and electron affinity parameters were used to determine the stability and reactivity. The results show the singlet cluster (Al2,4,6C1,2 and Al3,5,7C0-20,±) have the greater stability than the doublet clusters (Al3,5,7C0-2 and Al2,4,6C1,20,±) and triplet clusters (Al2,4,6). It has been demonstrated that anionic clusters are more stable than the neutral species, and cationic species are less stable than that of neutral species. The current study could be useful in understanding the stability and reactivity of metal clusters after carbon doping, which is important in combustion, material science, and the astrochemical community.

Iridium/Zinc-Co-Catalyzed Ring-Opening Reactions of Oxabicyclic Alkenes with Indole Nucleophiles: A Combined Experimental and Theoretical Study

An experimental and theoretical investigation on the iridium/zinc-co-catalyzed ring-opening reactions of oxabicyclic alkenes with indole nucleophiles is reported. The reaction affords trans-3-indolyl-1,2-dihydronaphthalen-1-ol products in good yield with no N-alkylated products observed. The C–C bond-forming reaction does not require prior functionalization and is entirely atom-economic. The mechanism and origins of selectivity in the iridium-catalyzed ring-opening reaction have been examined at the M06-D3/Def2TZVPP level of theory. Orbital analysis and natural population analysis charges demonstrate that the Cα site of the π-allyliridium intermediate is the most electrophilic site of attack. Distortion/interaction analysis reveals that the chemoselectivity likely originates from an earlier-stage transition state between Cα and C3, which requires less distortion energy compared to Cα–N1 during the rate-determining intermolecular nucleophilic attack. Moreover, conceptual density functional theory was used to conceptualize the preferential reactive sites of the nucleophiles probed. The C–C bond-forming step is speculated to proceed through a Friedel–Crafts-type reaction.

Theoretical study of multi-coordinated Xe(AuF)n (n = 2–4): Intriguing bond-bending isomerism

This work theoretically investigates the novel family of multi-coordinated noble gas (Ng) complexes Xe(AuF)n (n = 2–4). Geometrical optimization at the LC-ωPBE-D3/def2TZVPP level locates two types of structures based on the bond-bending isomerism arising from an intramolecular competition between Xe-Au covalent bonding and Au···Au aurophilic interactions. The isomers all comprise a center Xe atom that serves as an electron donor to bind to AuF moieties, but differ pronouncedly in the Au-Au spacing and the Au-Xe-Au bending angle. Electron and bonding analysis demonstrate that aurophilic interactions exist and affect the strength of Xe-Au bond whose covalent nature is evident. By potential energy surface scanning of Xe(AuF)2 with arbitrary-one atom being substituted by its analogue, we found only highly electronegative F/Cl atom can induce strong Xe-Au covalent bond to compete against aurophilic interactions. This work could be an interesting addition to the Ng chemistry and call out more experimental attentions.

Kinetics of IO radicals with C1, C2 aliphatic alcohols in tropospherically relevant conditions.

Kinetics of the reaction of IO radicals with methanol (MeOH) and ethanol (EtOH) were experimentally studied in the gas phase using pulsed laser photolysis-cavity ring-down spectroscopy (PLP-CRDS). IO radicals were produced in situ at the reaction zone by photolysing a mixture of precursors (CH3I + O3 + N2) at 248nm and thereby electronically excited at 445.04nm. The rate coefficients for the reactions of (IO + MeOH) and (IO + EtOH) were measured at a total pressure of 60Torr/N2 in the range of 258-360K. At room temperature, the experimental rate coefficients of the title reactions were measured to be [Formula: see text] and [Formula: see text]. Dependencies of the kinetics with photolysis laser fluence and experimental pressures were verified. Effects of pressure over the kinetic behaviour of the studied systems were observed to be insignificant within the statistical uncertainties when studied in the range of ~ 30-150Torr/N2, whereas a minor and linear fluence dependency was observed within the studied limit. From the measured kinetic parameters, the atmospheric lifetimes of MeOH and EtOH were calculated in the tropospherically relevant conditions regarding their reactions with important atmospheric oxidants like Cl atom, OH and IO radicals. To complement experimental results, kinetics and thermochemistry for the title reactions were investigated theoretically via canonical variational transition state (CVT) theory in combination with small curvature tunnelling (SCT) corrections with a dual-level Interpolated Single Point Energy (ISPE) approach at the CCSD(T)/def2-QZVPP//M06-2X/def2-TZVPP level of theory/basis set in the temperatures between 200 and 400K. Good degree of agreement was encountered between experimentally measured and theoretically calculated rate coefficients. This article also discusses the thermochemical parameters and kinetic branching ratios (BRs) of all the pathways involved in the title reactions.

Investigations into the flotation of molybdenite in the presence of chalcopyrite using (3S,4S,5S,6R)-3,4,5,6-tetrahydroxyoxane-2-carboxylate acid as a novel selective depressant: An experimental and theoretical perspective

(3S,4S,5S,6R)-3,4,5,6-tetrahydroxyoxane-2-carboxylate (TTC) acid was experimentally and theoretically investigated as a chalcopyrite (CPY) depressant during molybdenite-chalcopyrite (MLT-CPY) beneficiation. It was observed that the adsorption of TTC onto CPY was more than that onto MLT. The zero-point corrected complexation energies for complexes on mineral surfaces was calculated by using M06-2X/6–311++g(d,p)/def2-TZVPP level of theory in gas phase. It was established that Mo(OH)4-∙∙∙FeOH+∙∙∙TTC (S = 1) and FeOH+∙∙∙TTC (S = 0) complexes on MLT and CPY, respectively, of complexation energies of −189.0 kJ/mol and −831.0 kJ/mol, respectively, were the most stable complexes. The analyzed nature of bonding in complexes using Non-Covalent Interaction and Quantum Theory of Atoms in Molecules (QTAIM) suggested that the complexes were governed by traditional hydrogen bonding, as well as covalent, and other non-covalent interactions. Flotation tests carried out using a synthetic mixture of MLT and CPY (1:1) at pH 5.4 ± 0.1 1 in the presence an economical TTC depressant yielded a grade, recovery and separation efficiency of 58.7 %, > 97.6 % and 88.8 %, respectively. Therefore, TTC could be used as a potential depressant of CPY during the flotation of MLT.

Iron-Catalyzed Cross-Electrophile Coupling of Inert C–O Bonds with Alkyl Bromides

Iron-Catalyzed Cross-Electrophile Coupling of Inert C–O Bonds with Alkyl Bromides

The unimolecular decomposition of dimethoxymethane: channel switching as a function of temperature and pressure.

Branching ratios of competing unimolecular reactions often exhibit a complicated temperature and pressure dependence that makes modelling of complex reaction systems in the gas phase difficult. In particular, the competition between steps proceeding via tight and loose transition states is known to present a problem. A recent example from the field of combustion chemistry is the unimolecular decomposition of CH3OCH2OCH3 (DMM), which is discussed as an alternative fuel accessible from sustainable sources. It is shown by a detailed master equation analysis with energy- and angular-momentum-resolved specific rate coefficients from RRKM theory and from the simplified statistical adiabatic channel model, how channel switching of DMM depends on temperature and pressure, and under which experimental conditions which channels prevail. The necessary molecular and energy data were obtained from quantum-chemical calculations at the CCSD(F12*)(T*)/cc-pVQZ-F12//B2PLYP-D3/def2-TZVPP level of theory. A parameterization describing the channel branching over extended ranges of temperature and pressure is derived, and the model is used to simulate shock tube experiments with detection by atomic resonance absorption spectroscopy and time-of-flight mass spectrometry. The agreement between the simulated and experimental concentration-time profiles is very good. The temperature and pressure dependence of the channel branching is rationalized, and the data are presented in a form that can be readily implemented into DMM combustion models.

Examination of How Well Long-Range-Corrected Density Functionals Satisfy the Ionization Energy Theorem.

We calculated the vertical ionization energies (VIE) of 99 species in two ways to examine the accuracy of several long-range-corrected (LC) hybrid meta functionals in comparison with a gradient approximation (GA), global hybrids, and doubly hybrids. In the category of LC functionals, we examined both those with meta ingredients (i.e., that depend on the kinetic energy density) and those without them. The LC-hybrid meta functionals examined are M11, revM11, M11plus, and ωB97M-V. The reference data used to assess accuracy consist of 95 molecules and 4 atoms in the GW100 set. The two methods studied are the ΔSCF method (involving the difference of neutral and cation self-consistent field (SCF) energies) and the ionization energy theorem (involving the orbital energy of the highest occupied molecular orbital, HOMO). We calculated linear correlation coefficients (r2) and mean absolute deviations (MADs) between each approach and the reference VIE value from the CCSD(T)/def2-TZVPP level of theory. We compared the new LC-hybrid meta calculations to calculations with the 10 functionals in a previous VIE study by Brémond et al. and to the calculations with LC-BLYP (LC-Becke, Lee-Yang-Parr), CAM-B3LYP (Coulomb-attenuating-method Becke-3-parameter Lee-Yang-Parr), LC-ωHPBE, and ωB97X-D. The results show that Minnesota LC-hybrid meta functionals have the smallest mean absolute deviation of ionization energy theorem VIEs with the reference data; the LC-ωHPBE functional also does quite well in this test. This is very encouraging and indicates that LC-hybrid meta functionals would be the best starting points for the tuning strategy that has been shown to be a very good procedure for improving time-dependent density functional calculations, and it also helps explain the good success of LC-hybrid meta functionals for molecular excitation energies.