Potential energy profiles of the geometric isomerization and the thermal decomposition of diphosphene HP=PH in the ground and excited electronic states

Summary The geometric isomerization and the dehydrogenation of HP=PH in the ground and some low-lying excited states are investigated by theoretical calculations. The reaction paths are traced by either the CASSCF or UHF-SCF calculations using the 6-31G(d,p) basis functions, and the accompanying energy changes are calculated by the MRD-CI method employing the [5s3p1d]/[2s1p] basis functions. The barrier heights for the trans-to-cis isomerization, by the planar inversion and the nonplanar twisting, in the ground state are calculated to be 265 and 144 kJ/mol (with the vibrational zero-point energy corrections), respectively. The latter barrier is noticeably lower than the H-P and the P-P bond dissociation energies oftrans-HP=PH (¹A_g), which are 304 and 271 kJ/mol, respectively. The ground-state HP₂ radical (²A'), which is to be formed by the dehydrogenation of HP=PH, should suffer further decomposition into P₂ (¹ _g ⁺ ) and H with an activation energy of 139 kJ/mol. The lowest excited state of HP₂ is found to be a hydrogen-bridged 3-electron system (²A₂) having an isosceles triangle structure. It has proved to be formed by the dehydrogenation of the lowest excited singlet state (¹B) of HP=PH via a transition state which lies 194 kJ/mol above the¹B state. The excited HP₂ (²A₂) is state-correlated with P₂ (³_u)+H.     

Theoretical study of hydrogen bonded complexes of dimethyl disulfide or dimethyl peroxide with nitric acid

Ab initio and DFT calculations performed on the title systems revealed two types of structures for both DMDS-HNO₃ and DMDO-HNO₃ complexes. In both structures two hydrogen bonds are formed between the OH group interacting with one of sulfur (or oxygen) atoms and methyl CH group being a proton donor to one of the oxygen atoms of the NO₂ group of nitric acid. Depending on the location of the interacting methyl group with respect to the S or O acceptor of the main O-H?S(O) bond, the seven or eight-membered ring structures are formed. For all the structures, the most pronounced changes in geometric parameters upon interaction are observed for the proton donor molecule. The calculated binding energies are between −20.86 and −29.95 kJ/mol at MP2 and between −17.52 and −27.47 kJ/mol at B3LYP using the 6-311++G(2d,2p) basis set. The complexes involving disulfide are slightly weaker by ca. 6.7-8.6 kJ/mol than the corresponding peroxide complexes. The performed NBO analysis reveals that the charge transferred to σ*(OH) orbital of the nitric acid molecule comes mainly from the high p-character lone pair orbital of sulfur or oxygen atom being the hydrogen bond acceptor site in the disulfide or peroxide molecule.     

Topological analysis of the charge density response of a Ni4 cluster to a probe H2 molecule

Local density self-consistent field (SCF) discrete variationalX calculations are performed on a Ni₄ tetrahedron interacting with a probe H₂ molecule in special geometries. Optimized basis functions generated from the spherically averaged SCF potential are used. Topological charge-density analyses and binding energy calculations are used to study a portion of the energy surface for the approach of the H₂ molecule toward the Ni₄ tetrahedron. The effect of the H₂ molecule on Ni-Ni, Ni-H bonds and changes in the H-H covalent bond are investigated with the help of the field and various data at its critical points. The qualitative relationship between these data and the calculated binding energies is exploited.     

Molecular Structure of Silatrane Determined by Gas Electron Diffraction and Quantum-Mechanical Calculations

The geometry of silatrane HSi(OCH₂CH₂)₃N has been determined by gas electron diffraction, ab initio calculations, and vibrational spectroscopy of crystal. Using the scaled force field from DFT calculations the amplitudes and perpendicular corrections were calculated. It was assumed that the silatrane molecule has C ₃ symmetry. The following values (r _g bond lengths in Å and _a bond angles in deg. with three standard deviations from the least-squared refinements using a diagonal weight matrix) are: SiN 2.406(27); NC 1.443(7); OC 1.399(11); SiO 1.648(3); CC 1.504(15); NSiO 78.8(21); SiOC 128.1(11); SiNC 105.4(14); CCO 117.0(26); CCN 108.2(30); CNC 113.2(17); OSiO 116.3(13). The 5-membered rings are flattened. The sum of its bond angles is equal to 537.5(42). It is shown that a very large difference is found for Si—N distance from ab initio and DFT calculating.     

Formation of crystalline TiO2−xNx and its photocatalytic activity

Amorphous precursors to nitrogen-doped TiO₂ (NTP) and pure TiO₂ (ATP) powders were synthesized by hydrolytic synthesis and sol-gel method (SGM), respectively. Corresponding crystalline phases were obtained by thermally induced transformation of these amorphous powders. From FT-IR and XPS data, it was concluded that a complex containing titanium and ammonia was formed in the precipitate stage while calcination drove weakly adsorbed ammonium species off the surface, decomposed ammonia bound on surface of precipitated powder and led to substitution of nitrogen atom into the lattice of TiO₂ during the crystallization. The activation energies required for grain growth in amorphous TiO_2−xN_x and TiO₂ samples were determined to be 1.6 and 1.7 kJ/mol, respectively. Those required for the phase transformation from amorphous to crystalline TiO_2−xN_x and TiO₂ were determined to be 129 and 142 kJ/mol, respectively. A relatively low temperature was required for the phase transformation in NTP sample than in ATP sample. The fabricated N-doped TiO₂ photocatalyst absorbed the visible light showing two absorption edges; one in UV range due to titanium oxide as the main edge and the other due to nitrogen doping as a small shoulder. TiO_2−xN_x photocatalyst demonstrated its photoactivity for photocurrent generation and decomposition of 2-propanol (IPA) under visible light irradiation ().     

Electronic structure and the hydrogen-shift isomerization of hydrogen nitryl HNO2

Summary Electronic structure of hydrogen nitryl HNO₂, a yet not identified entity, and the path of its possible isomerization totrans-HONO have been investigated byab initio SCF and MRD-CI computations using the 6-31G** basis set. HNO₂ isC _2v -symmetric and its ground state (¹ A ₁) is less stable thantrans-HONO by 66 kJ/mol (with the SCF vibrational zero-point energy correction). The lowest two excited singlet states (¹ A ₂ and¹ B ₁) are nearly degenerate, their vertical excitation energies being predicted to be 4.8 eV. The isomerization path is traced by the CASSCF procedure and the activation barrier height is evaluated by the CI treatment. HNO₂ in its ground state isomerizes totrans-HONO by maintaining the planar (C _s-symmetric) structure. The activation energy is calculated to be 171 kJ/mol, which is clearly lower than the calculated H-N bond energy (253 kJ/mol). The transition state seems to be more adequately described as an interacting system of proton and the nitrite anion rather than as a pair of two fragment radicals.     

Molecular structure,conformational mobility,vibrational spectra,and thermochemistry of peroxyacetic acid: an ab initio and density functional study

The structure of the peroxyacetic acid (PAA) molecule and its conformational mobility under rotation about the peroxide bond was studied by ab initio and density functional methods. The free rotation is hindered by the trans-barrier of height 22.3 kJ mol^–1. The equilibrium molecular structure of AcOOH (C _s symmetry) is a result of intramolecular hydrogen bond. The high energy of hydrogen bonding (46 kJ mol^–1 according to natural bonding orbital analysis) hampers formation of intermolecular associates of AcOOH in the gas and liquid phases. The standard enthalpies of formation for AcOOH (–353.2 kJ mol^–1) and products of radical decomposition of the peroxide — AcO^· (–190.2 kJ mol^–1) and AcOO^· (–153.4 kJ mol^–1) — were determined by the G2 and G2(MP2) composite methods. The O—H and O—O bonds in the PAA molecule (bond energies are 417.8 and 202.3 kJ mol^–1, respectively) are much stronger than in alkyl hydroperoxide molecules. This provides an explanation for substantial contribution of non-radical channels of the decomposition of peroxyacetic acid. The electron density distribution and gas-phase acidity of PAA were determined. The transition states of the ethylene and cyclohexene epoxidation reactions were located (E _a = 71.7 and 50.9 kJ mol^–1 respectively).     

A comparison of quantitative theoretical results of the bonding in Ni(CO)4 and Ni(N2)4

Using non-empirical calculations the details of bonding in Ni(CO)₄ and in the analogous Ni(N₂)₄ are investigated.For Ni(CO)₄ some previous results are confirmed. In the calculation on Ni(N₂)₄ the close resemblance with Ni(CO)₄ is quite remarkable. The main difference is contained in the fact that carbon has a lower -electron density than nitrogen and that therefore the *-orbital in CO is lower in energy and geometrically more favourable for back donation.From the calculations we find a difference in metal-ligand bond energy between the carbonyl complex and the dinitrogen complex of approximately 18 kcal/mol.     

Aromaticity or non aromaticity in germanazenes? Theoretical studies on germanazenes, the N-cyano analogs and their carbodiimide isomers

Theoretical studies of germanazene rings [(Ge_II-NR)_2,3; R = H, Me, CN, Ph] have been performed at the DFT/B3LYP level. The fully optimized geometrical structures display four or six-membered planar rings of alternating germanium and nitrogen, in good agreement with the available X-ray experimental data. The hypothetical molecule (GeN-H)₂ presents only a small distortion from planarity. Although the planar conformation could indicate some degree of delocalization, the stabilization energy - estimated using the concept of homodesmotic reactions - indicates very little or no aromatic character in these molecules. The easy experimental formation of these germanazenes can be explained by di- (or tri-)merisation of the transient monomeric germylene-imine GeNR in its triplet state. When R = CN, in conformity with the experimental results, the most stable species is the isomeric carbodiimide form (GeNCN)_n, a result which is easily explained by the maximum spin density on the terminal nitrogen in the calculated monomer.     

Molecule dynamics and combined QM/MM study on one-carbon unit transfer reaction catalyzed by GAR transformylase

Both a molecule dynamic study and a combined quantum mechanics and molecule mechanics (QM/MM) study on Glycinamide ribonucleotide transformylase (GAR Tfase) catalytic mechanism are presented. The results indicate a direct one-carbon unit transfer process but not a stepwise mechanism in this reaction. The residues near the active center can fix the cofactor (N10-formyltetrahydrofolate) and GAR in proper relative positions by a H-bond network. The transition state and the minimum energy pathway are located on the potential energy surface. After all the residues (including H₂O molecules) are removed from the system the activation energy has increased from 145.1 kJ/mol to 243.3 kJ/mol, and the formly transfer reaction is very hard to achieve. The interactions between coenzyme, GAR and residues near the reactive center are discussed as well.     

Conformation analysis of diphosphine

Extended basis set ab initio computations are performed on HF, PNO-CI and CEPA level to determine the structure of P₂H₄ and the potential curve E() for rotation around the P-P axis. The structure parameters are optimized for dihedral angles of 0 ° (cis), 50 °, 80 ° (gauche or semi-eclipsed), 130 °, and 180 ° (trans). It turns out that P₂H₄ has a gauche equilibrium structure, a local minimum for trans which is 2.5 kJ/mol above gauche, a rather large cis barrier of 20 kJ/mol and a gauche trans barrier of 3.5 kJ/mol. The potential E() is extremely flat in the region 50 ° < < 310 °, where E() varies by less than 5 kJ/mol. Electron correlation tends to reduce the barriers but has no drastic effect on E().     

DFT Study of Methylene Blue Adsorption on ZnTiO3 and TiO2 Surfaces (101)

The search for alternative materials with high dye adsorption capacity, such as methylene blue (MB), remains the focus of current studies. This computational study focuses on oxides ZnTiO₃ and TiO₂ (anatase phase) and on their adsorptive properties. Computational calculations based on DFT methods were performed using the Viena Ab initio Simulation Package (VASP) code to study the electronic properties of these oxides. The bandgap energy values calculated by the Hubbard U (GGA + U) method for ZnTiO₃ and TiO₂ were 3.17 and 3.21 eV, respectively, which are consistent with the experimental data. The most favorable orientation of the MB adsorbed on the surface (101) of both oxides is semi-perpendicular. Stronger adsorption was observed on the ZnTiO₃ surface (−282.05 kJ/mol) than on TiO₂ (–10.95 kJ/mol). Anchoring of the MB molecule on both surfaces was carried out by means of two protons in a bidentate chelating (BC) adsorption model. The high adsorption energy of the MB dye on the ZnTiO₃ surface shows the potential value of using this mixed oxide as a dye adsorbent for several technological and environmental applications.     

Stability of functionalized C60 paramagnetic dimers and monomers

Density functional theory is used to calculate the bond dissociation energy to cleave the C₆₀C₆₀ bond of the paramagnetic X-C₆₀C₆₀-X and X-C₆₀C₆₀ dimers where X is F, OH, O and H. The results show that these dimers would not be stable much above room temperature and therefore cannot constitute the paramagnetic phase needed to form the observed ferromagnetism which has been shown to be stable up to 800 K. The calculated bond dissociation energies to remove an F, OH or H from a single C₆₀ are large suggesting that they could be the source of the unpaired spin needed for the high temperature ferromagnetism.     

Homodesmotic method of determining the O–H bond dissociation energies in phenols

The dissociation energy of the O–H bond has been calculated by the homodesmotic reaction method for phenolic compounds, which are well-known antioxidants, including for natural phenols. Use of moderately complex computational levels, such as B3LYP/6-31G(d), is sufficient for reliably estimating the D(O–H) value for phenols within the homodesmotic approach. The O–H bond dissociation energy for monosubstituted phenols has been calculated, and the additive character of the effect of methyl groups on D(O–H) in methylphenols has been demonstrated: the introduction of a CH₃ group into the aromatic ring decreases the D value by 7.8 kJ/mol (ortho position), 1.8 kJ/mol (meta position), and 7.6 kJ/mol (para position). The O–H bond strength has been calculated for a number of ubiquinols, selenophens, flavonoids, and chromanols. The D(O–H) value recommended for α-tocopherol is 328.0 ± 1.3 kJ/mol.     

Synthesis and characterization of three new fluorovanadate complexes: N(C2H5)4[VOF4], N(CH3)4[VOF3Cl], N(C4H9)4[VOF3Br] and theoretical calculations of VOF4, VOF3Cl and VOF3Br ions

The reaction of VOF₃ with (C₂H₅)₄NF, (CH₃)₄NCl and (C₄H₉)₄NBr salts in anhydrous CH₃CN produced new complexes with the anion general formula [VOF₃X]⁻ in that (X = F⁻, Cl⁻, Br⁻). These were characterized by elemental analysis, IR, UV/Visible and ¹⁹F NMR spectroscopy. The optimized geometries and frequencies of the stationary point are calculated at the B3LYP/6-311G level of theory. Theoretical results showed that the VX (X = F, Cl, Br) bond length values for the [VOF₃X]⁻ in compounds 1-3 are 1.8247, 2.4031 and 2.5595 Å, respectively. Also, the VF₅ bond length values in [VOF₃X]⁻ are 1.824, 1.812 and 1.802 Å, respectively. These results reveal that the bond order for VX bonds decrease from compounds 1 to 3, while for VF₅ bonds, the bond orders increase. It can be concluded that the decrease of VX bonds lengths and the increase of VF₅ bond lengths in compounds 1-3 result from the increase of the hyperconjugation from compounds 1 to 3. Harmonic vibrational frequencies and infrared intensities for VOF₄⁻, VOF₃Cl⁻ and VOF₃Br⁻ are studied by means of theoretical and experimental methods. The calculated frequencies are in reasonable agreement with the experiment values. These data can be used in models of phosphoryl transfer enzymes because vanadate can often bind to phosphoryl transfer enzymes to form a trigonal-bipyramidal structure at the active site.     

DFT quantum-chemical calculations of nitrogen oxide chemisorption and reactivity on the Cu(100) surface

A DFT quantum-chemical study of NO adsorption and reactivity on the Cu₂₀ and Cu₁₆ metal clusters showed that only the molecular form of NO is stabilized on the copper surface. The heat of monomolecular adsorption was calculated to be ΔH _m = ?49.9 kJ/mol, while dissociative adsorption of NO is energetically unfavorable, ΔH _d = + 15.7 kJ/mol, and dissociation demands a very high activation energy, E _a = + 125.4 kJ/mol. Because of the absence of NO dissociation on the copper surface, the formation mechanism of the reduction products, N₂ and N₂O, is debatable since the surface reaction ultimately leads to N-O bond cleavage. As the reaction occurs with a very low activation energy, E _a = 7.3 kJ/mol, interpretation of the NO direct reduction mechanism is both an important and intriguing problem because the binding energy in the NO molecule is high (630 kJ/mol) and the experimental studies revealed only physically adsorbed forms on the copper surface. It was found that the formation mechanism of the N₂ and N₂O reduction products involves formation (on the copper surface) of the (O_adN-NO_ad) dimer intermediate that is chemisorbed via the oxygen atoms and characterized by a stable N-N bond (r _N-N ～1.3 Å). The N-N binding between the adsorbed NO molecules occurs through electron-accepting interaction between the oxygen atoms in NO and the metal atoms on the “defective” copper surface. The electronic structure of the (O_adN-NO_ad) surface dimer is characterized by excess electron density (ON-NO)^δ? and high reactivity in N-O_ad bond dissociation. The calculated activation energy of the destruction of the chemisorbed intermediate (O_adN-NO_ad) is very low (E _a = 5–10 kJ/mol), which shows that it is kinetically unstable against the instantaneous release of the N₂ and N₂O reduction products into the gas phase and cannot be identified by modern experimental methods of metal surface studies. At the same time, on the MgO surface and in the individual (Ph₃P)₂Pt(O₂N₂) complex, a stable (O_adN-NO_ad) dimer was revealed experimentally.     

Characteristics of the hydrogen bond and the structure of Н-complexes of <Emphasis Type="Italic">p</Emphasis>-n-propyloxybenzoic acid and <Emphasis Type="Italic">p</Emphasis>-n-propyloxy-<Emphasis Type="Italic">p</Emphasis>′-cyanobiphenyl

The conformational properties of p-n-propyloxybenzoic acid and p-n-propyloxy-p′-cyanobiphenyl molecules, which can exhibit liquid crystalline properties in the formation of Н-complexes, are studied (DFT/B3LYP)/cc-pVTZ method). It is found that a molecule of p-n-propyloxybenzoic acid has 16 conformers that can be divided into four groups with respect to relative energies (0 kcal/mol, 1.6 kcal/mol, 6.5 kcal/mol, and 8.1 kcal/mol), and a molecule of p-n-propyloxy-p′-cyanobiphenyl has six conformers with relative energies of 0 kcal/mol (two conformers, φ(С_Ph–O–C–C)=180°) and 1.6 kcal/mol (four conformers, φ(С_Ph–O–C–C)=64.4°). In all conformers of the 3-AOCB molecule, phenyl rings are turned at 35° relative to each other. A conformation with the planar arrangement of two rings has a higher energy by 1.5 kcal/mol. Barriers to the internal rotation of different groups are determined and it is established that the structural nonrigidity of the molecules is mainly due to the possible rotation of the–C₂Н₅ moiety about the C–C bond. It is shown that with increasing temperature the vibrational amplitudes of the OC₃H₇ substituent, which enhance the probabilities of transitions between the conformers, become appreciably larger. It is found that p-n-propyloxybenzoic acid and p-n-propyloxy-p′-cyanobiphenyl can form Н-complexes with the medium hydrogen bond. Two types of the structural organization of Н-complexes are considered: linear and angular. The similarity of energies of Н-complexes with different structures (NBO analysis) can be the reason for the occurrence of two liquid crystalline subphases of p-n-propyloxybenzoic acid and p-n-propyloxy-p′-cyanobiphenyl system.     

Theoretical studies on the structures and isomerization of the LiSiF3 system

Various possible isomers of LiSiF₃ system and isomerization between them have been studied at G2(MP2) level usingab initio calculations. The relative energies of four minimum points on the potential energy surface are-128.6,-194.3,-12.7 and-122.8 kJ/mol (taking the sum of the energies of LiF and SiF₂ as zero). The structural energy of the four-membered ring that contains three F-Si-F-Li four-membered rings with C_3v symmetry is the lowest. The highest potential barrier for the isomerization of the remaining three- or four-membered structure is 12.5 kJ/ mol. Project supported by the National Natural Science Foundation of China (Grant No. 29673026).     

Matrix-isolation FT-IR and theoretical investigation of the competitive intramolecular hydrogen bonding in 5-methyl-3-nitro-2-hydroxyacetophenone

The paper presents the conformational, vibrational and hydrogen bond characteristics of 5-methyl-3-nitro-2-hydroxyacetophenone studied with the combined matrix-isolation FT-IR spectroscopic and theoretical (DFT/B3LYP/6-31++G^**) technique. Theoretical calculations predict three stable conformations of the studied compound. Only two of these conformations could be identified experimentally using the matrix-isolation FT-IR technique. The conformation with the intramolecular hydrogen bond OH^…ON has been found to be more stable than the conformation with the OH^…OC type of hydrogen bond by 7.28 kJ/mol. The complete assignment of the experimental spectra could be performed based on the theoretical calculations including the normal coordinate analysis and isotopic substitution.