Dioxirane oxidation of nitrosoamines. An ab initio study

1999 ◽  
Vol 77 (1) ◽  
pp. 74-85 ◽  
Author(s):  
Gennady V Shustov ◽  
Arvi Rauk

Pathways for oxidation of the parent nitrosoamine 2, nitrosodimethylamine 3, and anti- nitrosoethylmethylamine 4 by the parent dioxirane 1 have been explored computationally using the B3LYP hybrid density functional theoretical method in conjunction with the 6-31G* and 6-311+G** basis sets. Oxygen transfer from 1 to the nitrogen of the NO group (NO oxidation), yielding nitroamines 5,7, and 10, has the lowest activation barrier (15, 12.8, and 12 kcal mol-1 for 2, 3, and anti-4, respectively). Oxidation of the amine nitrogen (N oxidation) in 2, 3 leads to nitric oxide and nitroxyl radicals 6, 8 and is characterized by the highest activation energy: 28.5 kcal mol-1 for 2, 22.3 kcal mol-1 for 3. The potential barrier to hydroxylation of the methyl groups in 3 (CH oxidation) is intermediate - ca. 19 kcal mol-1. Introduction of a methyl group to the carbon reaction centre decreases the activation barrier of the anti-CH oxidation by 2-3 kcal mol-1. In 3, a gas-phase small preference for anti-CH oxidation over syn-CH oxidation, 0.1 kcal mol-1, is predicted to increase dramatically to 4.7 kcal mol-1 when the reaction is carried out in a dielectric medium (acetone, IPCM model). In general, polar solvents (CH2Cl2, acetone) lower the activation barriers for the NO, N, and CH oxidations, more so for the species of higher polarity than those of lower polarity (reactive complexes (1+3), (1+anti-4), transition state anti,syn-TSCH(1+4)]. However, the chemoselectivity (EaNO <Eaanti-CH <Easyn-CH <EaN) in the polar solvents is predicted to be the same as in the gas phase.Key words: dioxirane, nitrosoamines, oxidation, ab initio, omputation.

2018 ◽  
Vol 71 (3) ◽  
pp. 102
Author(s):  
Emma Persoon ◽  
Yuekui Wang ◽  
Gerhard Raabe

Quantum-chemical ab initio, time-independent, as well as time-dependent density functional theory (TD-DFT) calculations were performed on the so far elusive heterocycles inda- and thallabenzene (C5H5In and C5H5Tl), employing several different methods (MP2, CISD, CCSD, CCSD(T), BD, BD(T), QCISD, QCISD(T), CASSCF, DFT/B3LYP), effective core potentials, and different basis sets. While calculations on the MP2 level predict the ground states of the title compounds to be singlets with the first triplet states between 13 and 15 kcal mol−1 higher in energy, single point calculations with the QCISD(T), CCSD(T), and BD(T) methods at CCSD-optimized structures result in energy differences between the singlet and the triplet states in the range between 0.3 and 2.1 kcal mol−1 in favour of the triplet states. According to a CASSCF(8,8) calculation the triplets are also more stable by about 2.5–2.9 kcal mol−1. Calculations were also performed for the C5v-symmetric η5 structural isomers (cyclopentadienylindium, CpIn, and cyclopentadienylthallium, CpTl, Cp = C5H5) of the title compounds. At the highest level of theory employed in this study, C5H5In is between 79 and 88 kcal mol−1 higher in energy than CpIn, while this energy difference is even larger for thallabenzene where C5H5Tl is energetically between 94 and 102 kcal mol−1 above CpTl. In addition we report on the UV/vis spectra calculated with a TD-DFT method as well as on the spectra of the normal modes of C5H5In and C5H5Tl. Both types of spectra might facilitate identification of the title compounds eventually formed in photolysis or pyrolysis experiments.


2013 ◽  
Vol 3 ◽  
pp. 69-73
Author(s):  
Kapil Adhikari ◽  
Asok K. Ray

Ab initio calculations of the electronic structures of silicon carbide (SiC) nanotubes represented by clusters are presented. The nanotube clusters of chiralities (3,3) and (5,5) are studied using the hybrid density functional B3LYP (Becke’s 3-parameter and the Lee-Yang-Parr exchange-correlation) and LANL2DZ (Los Alamos National Laboratory double ?) and 3-21G* basis sets. Evolution of electronic properties of silicon carbide (SiC) nanotubes (3, 3) and (5,5) with their length is discussed. The results suggest that the electronic properties of nanotubes change for short tubes of unit cells with lengths varying from 1 to 5. However, the properties do not seem to change significantly after this. Therefore, an infinite silicon carbide (SiC) nanotube can be approximated by a nanotube cluster of 5 unit cells.The Himalayan PhysicsVol. 3, No. 3, July 2012Page : 69-73


2020 ◽  
Vol 2020 ◽  
pp. 1-9 ◽  
Author(s):  
P. O. Jomo ◽  
C. O. Otieno ◽  
P. W. O. Nyawere

We report the results of pressure-induced semiconductor-metal phase transition of the semiconducting chalcogenide compound KPSe6 under high pressure using the ab initio methods. The ground-state energy calculations were performed within density functional theory and the generalized gradient approximation using the pseudopotential method with plane-wave basis sets. The projector augmented-wave (PAW) pseudopotentials were used in our calculation. The optimized lattice parameters were found from total energy calculations as 13 Bohr, 1.6 Bohr, and 1.8 Bohr for cell dimensions one, two, and three, respectively, which are in good agreement with experimental calculations. At zero pressure, the material portrayed a semiconducting property with a direct bandgap of ≈1.7 eV. As we subjected the material to pressure, the band gap was observed to reduce until it disappeared. The phase transition from the semiconductor to metal was found to occur at ∼45 GPa, implying that the material underwent metallization as pressure was increased further.


2014 ◽  
Vol 69 (7) ◽  
pp. 349-359 ◽  
Author(s):  
Stefanie Mersmann ◽  
Halima Mouhib ◽  
Matthias Baldofski ◽  
Gerhard Raabe

1Quantum-chemical ab initio and time-dependent density functional theory (TD-DFT) calculations employing various basis sets were used to elucidate the spatial as well as the electronic structure of C5H5Al () and C5H5Ga (2) (ala- and galabenzene). The lowest closed shell singlet states of both compounds were found to have a non-planar structure of CS symmetry with C-X-C bond angles of about 116° (MP2/6-311++G**) and 125° (CCSD/aug-cc-pVDZ). At approximately 103°, the corresponding angles of the lowest triplets are significantly smaller. The lowest triplet state of alabenzene is also non-planar (CS) at the MP2 level while optimization with the CCSD and the CASPT2 method resulted in planar structures with C2v symmetry. The corresponding state of galabenzene has C2v symmetry at all levels of optimization. The relative stability of the lowest closed shell singlet and the lowest triplet (ΔE(T1-S0)) state is small and its sign even strongly method-dependent. However, according to the highest levels of theory applied in this study the singlet states of both molecules are slightly lower in energy than the corresponding triplets with singlet/triplet gaps between about 0.5 and 5.8 kcal/mol in favour of the singlet states. Most of the applied methods give a slightly smaller splitting for ala- than for galabenzene. Independent of the applied method (TD-DFT/CAM-B3LYP/6-311++G(3df,3pd)//MP2/6- 311++G** or SAC-CI/6-31++G(3df,3pd)//MP2/6-311++G**), the general shape of the calculated UV/VIS spectral curves are quite similar for the lowest singlet states of ala- and galabenzene, and the same applies to the spectra of the normal modes. The calculated UV/VIS spectra of C5H5Al and C5H5Ga are featured by long wavelength bands of moderate intensity around 900 nm at the TD-DFT and between 1300 and 1500 nm at the SAC-CI level. According to both methods these bands are predominantly due to HOMO(π)→LUMO(σ*) transitions. The results of isodesmic bond separation reactions for the singlet states indicate some degree of stabilization due to delocalization in both of the title compounds. With our best values between 29 and 32 kcal/mol this stabilization appears to be only slightly less than the previously reported value for borabenzene (∼38 kcal/mol).


1999 ◽  
Vol 77 (5-6) ◽  
pp. 540-549 ◽  
Author(s):  
Gennady V Shustov ◽  
Michael TH Liu ◽  
K N Houk

The reactions of the singlet methylene (1a) and dimethylcarbene (1b), with their diazirine precursors, diazirine (2a), and dimethyldiazirine (2b), have been studied theoretically using ab initio and density functional theory. The reaction has no activation barriers for the parent system (1a + 2a) and proceeds via a reactive complex and a transition state with a small negative enthalpy of activation Δ Hnot =298 = -1.1 kcal mol-1, ΔSnot =298 = -34.4 cal mol-1 K-1, ΔG°298 = 9.2 kcal mol-1) for the dimethyl derivatives (1b + 2b). The formation of N-methylene diazirinium ylides (3a,b) is exothermic by 64-80 kcal mol-1. The isomer, 1,3-diazabicyclo[1.1.0]butane (4a), is more stable (5-12 kcal mol-1) than isomer 3a, but can neither be formed by direct thermal reaction of 1a with 2a nor undergo the direct rearrangement into formaldazine (5a). The rearrangement of ylides 3a,b into azines 5a,b proceeds by conrotatory C3-N1 ring opening. The predicted activation barrier of ca. 15 kcal mol-1 for the ring opening in ylide 3b is in excellent agreement with experimental data. The formation of pyridinium ylides from carbenes and pyridine is also studied.Key words: diazirinium ylide, ab initio MO (molecular orbital) theory, density functional theory, pyridinium ylide, CIS (singles configuration interaction) transition energies.


2021 ◽  
Author(s):  
Kazuumi Fujioka ◽  
Yuheng Luo ◽  
Rui Sun

Ab initio molecular dymamics (AIMD) simulation studies are a direct way to visualize chemical reactions and help elucidate non-statistical dynamics that does not follow the intrinsic reaction coordinate. However, due to the enormous amount of the ab initio energy gradient calculations needed for AIMD, it has been largely restrained to limited sampling and low level of theory (i.e., density functional theory with small basis sets). To overcome this issue, a number of machine learning (ML) methods have been employed to predict the energy gradient of the system of interest. In this manuscript, we outline the theoretical foundations of a novel ML method which trains from a varying set of atomic positions and their energy gradients, called interpolating moving ridge regression (IMRR), and directly predicts the energy gradient of a new set of atomic positions. Several key theoretical findings are presented regarding the inputs used to train IMRR and the predicted energy gradient. A hyperparameter used to guide IMRR is rigorously examined as well. The method is then applied to three bimolecular reactions studied with AIMD, including HBr+ + CO2, H2S + CH, and C4H2 + CH, to demonstrate IMRR’s performance on different chemical systems of different sizes. This manuscript also compares the computational cost of the energy gradient calculation with IMRR vs. ab initio, and the results highlight IMRR as a viable option to greatly increase the efficiency of AIMD.


2002 ◽  
Vol 80 (11) ◽  
pp. 1435-1443 ◽  
Author(s):  
J Komulainen ◽  
R S Laitinen ◽  
R J Suontamo

The structures and spectroscopic properties of SenS8–n ring molecules have been studied by the use of ab initio molecular orbital techniques and density functional techniques involving Stuttgart relativistic large core effective core potential approximation with double zeta basis sets for valence orbitals augmented by two polarization functions for both sulfur and selenium. Full geometry optimizations have been carried out for all 30 isomers at the Hartree-Fock level of theory. The optimized geometries and the calculated fundamental vibrations and Raman intensities of the SenS8–n molecules agree closely with experimental information where available. The nuclear magnetic shielding tensor calculations have been carried out by the Gauge-independent atomic orbital method at the DFT level using Becke's three-parameter hybrid functional with Perdew/Wang 91 correlation. The isotropic shielding tensors correlate well with the observed chemical shift data. The calculated chemical shifts provide a definite assignment of the observed 77Se NMR spectroscopic data and can be used in the prediction of the chemical shifts of unknown SenS8–n ring molecules.Key words: selenium sulfides, ab initio, DFT, effective core potentials, geometry optimization, energetics, fundamental vibrations, 77Se chemical shifts.


2021 ◽  
Author(s):  
Kazuumi Fujioka ◽  
Rui Sun

Ab initio molecular dymamics (AIMD) simulation studies are a direct way to visualize chemical reactions and help elucidate non-statistical dynamics that does not follow the intrinsic reaction coordinate. However, due to the enormous amount of the ab initio energy gradient calculations needed for AIMD, it has been largely restrained to limited sampling and low level of theory (i.e., density functional theory with small basis sets). To overcome this issue, a number of machine learning (ML) methods have been employed to predict the energy gradient of the system of interest. In this manuscript, we outline the theoretical foundations of a novel ML method which trains from a varying set of atomic positions and their energy gradients, called interpolating moving ridge regression (IMRR), and directly predicts the energy gradient of a new set of atomic positions. Several key theoretical findings are presented regarding the inputs used to train IMRR and the predicted energy gradient. A hyperparameter used to guide IMRR is rigorously examined as well. The method is then applied to three bimolecular reactions studied with AIMD, including HBr+ + CO2, H2S + CH, and C4H2 + CH, to demonstrate IMRR’s performance on different chemical systems of different sizes. This manuscript also compares the computational cost of the energy gradient calculation with IMRR vs. ab initio, and the results highlight IMRR as a viable option to greatly increase the efficiency of AIMD.


2021 ◽  
Author(s):  
◽  
Zhi Xiang Wong

<p>The reaction mechanism of the gold(III)-catalysed hydroalkoxylation of alkynes is studied to provide a deeper understanding of homogeneous gold catalysis. The study is conducted computationally using Density Functional Theory (DFT), with the PBE0 and BP86 functionals and basis sets of triple-ζ quality (aug-cc-pVTZ and aug-cc-pVTZ-PP for the gold atom). It emphasises the mechanisms undergone by various alkynes when they are activated by gold(III) catalysts towards nucleophilic attack to first form an enol ether and followed by a second nucleophilic attack to form a ketal as the final product. Hydrogen bonding networks formed by the solvent methanols are found to play a crucial role in the mechanism especially in the hydrogen migration steps that follow after the nucleophilic attacks. The first nucleophilic attacks are predicted to have rather low activation energies and hence they are expected to proceed fast while the second additions vary in activation barriers, depending on the steric effects in the substrates. The activation barrier for the last hydrogen migration is highest for all of the three reactions investigated and is expected to be the rate determining step. Investigations of internal alkyne reactions reveal that each elementary step requires a higher activation energy compared to terminal alkynes, which explains the low experimental rate of such reactions. Due to the regioselectivity problem in internal alkyne reactions, this results in a mixture of products which is difficult to isolate due to the similarities in their reaction energies. The study also highlights the calculated thermodynamics and kinetics of the reactions, which can be useful in predicting experimental outcomes. Arrhenius plots of concentration of each intermediate species against time were produced to further help the understanding of these mechanisms, whether or not the reactions go to full completion or stop at the formation of enol ether.</p>


1998 ◽  
Vol 53 (10) ◽  
pp. 1223-1235
Author(s):  
Inge Warttmann ◽  
Günter Häfelinger

AbstractAb initio Hartree-Fock (HF) and density functional (DFT) optimizations on the test m olecule osmiumtetracarbonyldihydride (13) with various basis sets show that the lanl2mb pseudopotential basis set for osmium leads in the HF approximation to more reliable molecular geometries than the DFT calculations. This HF procedure was used for the optimizations of molecular geometries of three isomeric 4,4,4,4,17,17,17,17-octacarbonyl-4,17-diosma[7.7]ortho-, meta- and paracyclophanes 1 to 3, of which 3 was found to be predestined for formation of various host-guest complexes with possible guests benzene (4), fluorobenzene (5), 1,3,5- trifluorobenzene (6), 1,2,4,5-tetrafluorobenzene (7), hexafluorobenzene (8), fluoroanil (9), tetrafluoroethene (10), tetracyanoethene (11) and aniline (12). Results of optimized hostguest geometries are presented graphically for inclusions and associations of guest 4 to 12 with 3. Calculated lanl2mb interaction energies, after correction for basis set superposition error (BSSE), remain favourable only for inclusion of 5 and associations of 5, 10, 11 and 12. Additionally lanl2dz single point calculations for inclusion, which may not need BSSE correction because of the improved basis set, are favourable for 6 and 12. According to lanl2mb HOMO and LUMO energies, 3 may as well easily donate or accept electrons. This may be an interpretation to the surprising effect, that Mulliken total charges are positive on the electron accepting guest molecules 4 to 11. There are geometrical peculiarities in the optimized host-guest complexes for inclusion and association. Fluorine atoms of 5 to 10 and nitrogen atoms of a cyano group of 11 and the amino group of 12 like to come close to one or two carbonyl groups. Similar distances of 2.70 Å to 3.57 Å between the O atom of the carbonyl group and the F atom or N atom appear in all optimizations of inclusion and association of 5 to 12 except in the case of association of tetrafluoroethene (10).


Sign in / Sign up

Export Citation Format

Share Document