Theoretical study for the HNSi–HSiN isomerization: ab initio and DFT calculations

Ab initio and density functional theory (DFT) calculations using the GAUSSIAN 94 program have been performed to investigate the molecular structures of HNSi and HSiN in the ground state as well as the transition state for the HNSi–HSiN isomerization reaction at the 6-311G(d,p), 6-311+G(2d,p) and 6-311+G(2df,p) basis sets. The results show that DFT calculations at higher levels of theory reproduce experimental vibrational frequencies of both HNSi and HSiN better than ab initio methods including electron correlation effects. Those calculated geometries are accurate enough to predict the rotational constant of HNSi. The barrier height for the isomerization reaction is found to be about 10 kcal/mol.     

Theoretical study on second‐order nonlinear optical properties of spin crossover Fe(III) phenolate‐pyridyl Schiff base complexes

The density functional and ab initio theory were used to investigate the second‐order nonlinear optical (NLO) properties of Schiff base ligands, open‐shell Fe(III), and closed‐shell Ni(II) complexes. The effect of the metal center in complexes is thus manifold: it templates the formation of acentric structures, imparts high thermal stability to the chelate ring, and display higher second‐order NLO response than their ligands. The second‐order NLO response of metal complexes are intensively sensitive to the exchange donor/acceptor because the differences of the extent of charge separation and the intraligand charge transfer processes. Thus, substituted metal complexes could realize “switches on” the second‐order NLO response by exchange donor/acceptor. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Ring opening of boriranes vis‐à‐vis aziridines: An ab initio and DFT probe on the mechanisms

Borirane undergoes ring opening reaction with NOCl and HNF₂ yielding the corresponding alkenes. Ab initio and density functional investigations of this reaction with cis‐ and trans‐2,3‐dimethylboriranes reveal that these reactions take place in a single step through the formation of a prereactive complex and a transition state giving the alkene with the same stereochemistry. Calculations clearly show that the concerted cleavage of C? B bonds leads to retention of stereochemistry. Further, it shows that HNF₂ cleaves boriranes more efficiently than does NOCl. Intrinsic reaction coordinate analyses and bond order analysis describe the nature of the transition state very well and fix the reaction mechanism. Solvent effect calculations through PCM model, with acetonitrile and CCl₄ as solvents, do not alter the gas phase results significantly. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Density functional theory (DFT) and ab initio investigations of the electronic and molecular structures of the monomer and the dimer of trimethylaluminium

The electronic and molecular structures of the monomer and dimer of trimethylalu-minium have been studied using density functional theory and ab initio MP2 method. The optimized geometry of the monomer Al(CH3)3 is of C3h symmetry, whereas that of the dimer [A1(CH3)3]2 contains a carbon-bridged four-membered ring structure with C2h symmetry. The hydrogen-bridged six-membered ring structure is found to be unstable. The calculated dimerization energy for the four-membered ring structure is 78 kJ/mol, in close proximity to the experimental value of 85.27 kJ/mol.     

SIESTA‐SIPs: Massively parallel spectrum‐slicing eigensolver for an ab initio molecular dynamics package

Integration of Shift‐and‐Invert Parallel Spectral Transformation (SIPs) eigensolver (as implemented in the SLEPc library) into an ab initio molecular dynamics package, SIESTA, is described. The effectiveness of the code is demonstrated on applications to polyethylene chains, boron nitride sheets, and bulk water clusters. For problems with the same number of orbitals, the performance of the SLEPc eigensolver depends on the sparsity of the matrices involved, favoring reduced dimensional systems such as polyethylene or boron nitride sheets in comparison to bulk systems like water clusters. For all problems investigated, performance of SIESTA‐SIPs exceeds the performance of SIESTA with default solver (ScaLAPACK) at the larger number of cores and the larger number of orbitals. A method that improves the load‐balance with each iteration in the self‐consistency cycle by exploiting the emerging knowledge of the eigenvalue spectrum is demonstrated. © 2018 Wiley Periodicals, Inc.     

Theoretical study of the adsorption of isobutane over H-mordenite zeolite by ab initio and DFT methods

Values of relative interaction energy between an isobutane molecule and the surface of the H-mordenite zeolite were obtained. Sixteen tetrahedral sites were taken into account for the surface. Results showed that the adsorption energies were 4 and 8 kJ/mol, for ab initio restricted Hartree–Fock (RHF) method and density functional Becke–Lee–Yang Parr (B3-LYP) approaches, respectively, when C₂---H₂ bond is located in front of H⁺ atom of H-mordenite surface. The minimum value for the interaction in the adsorption was around 2.25 and 2.50 Å distance from the DFT and ab initio method, respectively. The interaction of isobutane by methyl groups side over the catalytic surface caused a repulsive interaction. Therefore, adsorption of isobutane is a selective process by the H₂ side of the molecule. Calculations were done by employing the 6-31G^** basis set.     

Electronic properties of the low‐lying spin states of dimethylnitrosamine coordinated to Fe(III) heme models: An ab initio study

The unusual O‐coordination mode of nitrosamines to Fe(III) heme models has been observed in the bis(dimethylnitrosamine)(meso‐tetraphenylporphyrinate)iron(III) cation. For the first time, this latter as well as the simpler bis(dimethylnitrosamine)(porphinate)iron(III) heme model cations have been studied through ab initio methods. The sextet, quartet, and doublet spin states of both cations have been studied through single‐point calculations based on the experimental (X‐ray) geometry. Their energies, charges, and spin densities have been analyzed. The obtained results (at the UHF/cc‐pVDZ and ROHF/cc‐pVDZ levels) indicate that the peripheral benzene rings are of secondary importance for the coordination of dimethylnitrosamine to the Fe(III) porphyrin core. The obtained energy ordering is sextet < quartet < doublet, at all computational levels. The UHF, ROHF, and UMP2 results indicate an excess of alpha spin density around the Fe atom, a low covalency for the Fe? O bond and a substantial charge transfer to the Fe atom. Our best estimates [obtained at ROMP2 level with the mixed cc‐pVDZ/cc‐pVTZ‐DK(Fe) basis set] for the energy differences (in eV) between the three spin states considered are 0.929 for the sextet‐quartet gap and 0.812 for the quartet‐doublet gap, which indicate that the spin crossover (at room temperature) is very unlikely. These results represent the substantial decrease in the uncorrelated values. The implications of spin contaminations at the UHF and UMP2 levels for subsequent geometry optimizations to be performed in the smaller cation have also been discussed. © 2013 Wiley Periodicals, Inc.     

High‐level ab initio computational study of doublet and quartet nitrogen reaction with methane

High levels of ab initio calculations were performed with the target of exploring the potential‐energy surface for the doublet and the quartet nitrogen with methane. There is a considerable difference between these two reaction paths in light of the formation of the reactant complex. Doublet nitrogen is an excited state of nitrogen, but it forms a stronger complex with methane, making the activation barrier with the doublet and quartet nitrogen reaction with methane. Activation barriers, heats of the reaction, and bond dissociation energies for minima located in the potential surface of nitrogen and methane transformations into the hydrogen radical and H₂CNH were evaluated, and the most probable reaction pathways were suggested. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 71: 481–490, 1999     

FTIR, FT-Raman spectra and ab initio DFT vibrational analysis of 2,4-dichloro-6-nitrophenol

The FTIR and FT-Raman spectra of 2,4-dichloro-6-nitrophenol (2,4-DC6NP) has been recorded in the region 4000-400 cm(-1) and 3500-100 cm(-1), respectively. The optimized geometry, frequency and intensity of the vibrational bands of (2,4-DC6NP) were obtained by the ab initio and DFT levels of theory with complete relaxation in the potential energy surface using 6-31G(d,p) and 6-311+G(d,p) basis sets. The harmonic vibrational frequencies were calculated and the scaled values have been compared with experimental FTIR and FT-Raman spectra. The observed and the calculated frequencies are found to be in good agreement. The experimental spectra also coincide satisfactorily with those of theoretically constructed bar type spectrograms.     

Reaction of OH radical with mono‐, di‐, and trichloroacetaldehyde: an ab initio study

The hydrogen abstraction reactions between chlorine‐substituted acetaldehydes and OH radicals have been investigated by using ab initio molecular orbital theory. Equilibrium geometries and transition‐state structures have been optimized at the (U)MP2/6‐311G(d,p) level. Activation barriers and heats of reaction for different reaction channels have been estimated from the single‐point calculations at the (U)MP2/6‐311G(2df,2p) level. Three, two, and one hydrogen abstraction channel have been found for the mono‐, di‐, and trichloroacetaldehyde, respectively. At a higher temperature region, hydrogen abstraction from the formyl group is found to be the major reaction channel for all the three chloroacetaldehydes. The effect of halogen substitution on reactivity toward hydrogen abstraction has been discussed. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1509–1521, 2001     

An ab initio and DFT study on the hydrolysis of carbonyl dichloride

Hydrolysis of carbonyl dichloride or phosgene (Cl₂CO) in gas phase has been investigated at Hartree–Fock, density functional and ab initio levels of theory. The effects of basis sets on the energetics of the reaction have also been explored. Calculations reveal that initially carbonyl dichloride and water form a weak complex and this complex can react further in two ways. In Path 1, water adds on to carbonyl dichloride across carbonyl bond in a concerted fashion to give dichloromethane diol, and this diol decomposes to form chloro formic acid by syn-1,2-elimination of HCl and forms CO₂ and HCl as final products. Path 2 is the concerted addition of water across carbon chlorine bond and elimination of HCl in a single step leading to the formation of chloro formic acid directly. This second path that skips the formation of dichloromethane diol is observed to be very low lying and hence is kinetically favored. Addition of second water molecule to the reacting system is found to catalyze the reaction by stabilizing the complex, intermediate and transition states and reduces the activation energy to 24.6 kcal mol⁻¹ compared to 29.9 kcal mol⁻¹ for a single water molecule.     

An ab initio and DFT study of structure and vibrational spectra of disiloxane H3SiOSiH3 conformers. Comparison to experimental data

The structural and vibrational properties of siloxane monomers may account in the physical and chemical properties of silicone polymers. Because disiloxane (H(3)SiOSiH(3)) is the smallest molecule in the set which runs through small siloxanes like hexamethyldisiloxane (CH(3))3SiOSi(CH(3))3 to silicone polymers, its energetic, structural and vibrational features have been investigated in detail using density functional theory (B3LYP), post Hartree-Fock methods (MP2 and CCSD(T)) and basis sets up to spdfg quality. Five conformations were considered: three bent structures with C2v (double staggered, SS, and double eclipsed, EE, conformations) and Cs symmetries, and two linear forms with D3d and D3h symmetries. At all levels of theory, the relative stability was C2v(SS) approximately C2v(EE)>Cs>D3h>D3d. The difference of energy between the two C2v conformers is lower than 0.04 kcal/mol. At the highest level of theory (CCSD(T)/cc-pVQZ), the barrier to linearisation from C(2v) to D(3h) conformers was calculated at 0.43 kcal/mol, which is extremely low. Most of the structural and vibrational features of the disiloxane do not depend on the conformation of the molecule but are strongly influenced by the SiOSi angle. Anharmonic calculations allowed, without any scaling factor, an exhaustive reinvestigation of the assignments of observed wavenumbers in the infrared and Raman spectra of gaseous disiloxane. Particularly, in the gas phase spectrum, the SiOSi symmetric and antisymmetric stretches have been assigned at 599 and 1105, 596 and 1060, 527 and 1093 cm(-1) for H(3)SiOSiH(3), H(3)Si(18)OSiH(3) and D(3)SiOSiD(3), respectively. The experimental wavenumber splitting of SiOSi symmetric and antisymmetric stretches of H(3)SiOSiH(3) gave an estimation of the SiOSi angle at around 145 degrees . Ab initio methods were revealed more accurate for structural parameters, when DFT/B3LYP was enough for spectral assignments, even at the harmonic level using a single scaling factor.     

苯并氧化呋咱稳定性和异构化的DFT和ab initio研究

运用B3LYP/6-31G(d)密度泛函理论（DFT）方法对苯并氧化呋咱、邻二亚硝基苯及其间的异构化反应进行了计算研究。结果表明,苯并氧化呋咱的分子总能量比邻二亚硝基苯的低;由苯并氧化呋咱异构为邻二亚硝基苯的正向反应活化能（Ea+=51.0kJ/mol）,与文献实测值（58.6kJ/mol）较接近,而其逆向反应活化能（Ea-=4.6kJ/mol）很小,从而揭示了苯并氧化呋咱比邻二亚硝基苯更稳定·此外,进行了HF/3-21G、HF/6-31G(d)和MP2/6-31G(d)//6-31G(d)水平下相应的计算,发现B3LYP-DFT的结果较abinitio为优。谐振动频率的B3LYP/6-31G(d)计算还支持了邻二亚硝基苯为苯并氧化呋咱“自-自”互变重排反应的中间体。     

Hydrogen bonding interaction in sarcosine–water complex using ab initio and DFT method

The hydrogen bonding interaction in the Sarcosine (N‐methylglycine)–water complex is studied using ab initio, MP2, and density functional theory (DFT/B3LYP). For this complex, binding energies, dipole–dipole interactions, chemical hardness, and chemical potential have been calculated. Three different basis sets, viz. 6‐311+G, 6‐311++G, and 6‐311++G*, have been used to optimize the geometries by all three methods. The basis set superposition errors are also calculated, and the corrected binding energies are reported for this complex. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

DFT plus PCM calculation for pairing specificity of Watson–Crick‐type bases in aqueous solutions

This study is aimed at explaining the preference for AT and CG pairings and the possible insertion of other tautomeric DNA base pairs such as G_enolT, that respect energetic and steric requirements including at least two hydrogen bonds and 11 ± 0.5Å distance between the 9‐CH₃ of purine and 5‐CH₃ of pyrimidine. The calculated free energy of formation ΔΔG at the DFT B3LYP/6‐31G*‐PCM/BEM level pointed out the CG and AT pairs as the most favored, followed closely by G_enolT, in good agreement with Michaelis–Menten first order kinetics (CG ≈ AT > G_enolT). Unusual DNA base pairs complexes such as AG (BEM) and CT (PCM) resulted to be stable, but it is very difficult to assume that they are likely to be included in the double strand DNA. The calculated enthalpy and dipole moments of isolated DNA bases agree well with experiment. The free energy of hydration, ΔG_hyd, was found to depend on the electrostatic term, while cavitation‐dispersion components are almost constant. The stability of DNA complexes in water resulted from PCM calculations is markedly influenced by the free energy of hydration. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

A comparative study of DFT and traditional ab initio methodologies on the OsO4 molecule

The performance of different conventional ab initio methodologies and density functional procedures is compared through its application to the theoretical calculation of the bond distance and harmonic vibrational frequencies of the OsO₄ molecule. The problem of the basis set is first considered, with up to nine different basis sets being tested in calculations using the hybrid Becke3LYP density functional, and the most appropriate basis set is used in the comparison of Hartree–Fock, post‐Hartree–Fock, and density functional methods. The post‐Hartree–Fock methods analyzed are MP2, CISD, and CCSD(T), and the density functionals tested are SVWN, BLYP, BPW91, and Becke3LYP. The results show that for this particular system density functional methods perform better than do HF‐based methods with the exception of CCSD(T), which gives the best overall results. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 544–551, 2000     

Electronic structure and molecular orbital study of the first excited state of the high‐efficiency blue OLED material bis(2‐methyl‐8‐quinolinolato)aluminum(III) hydroxide complex from ab initio and TD‐B3LYP

Bis(2‐methyl‐8‐quinolinolato)aluminum(III) hydroxide complex (AlMq₂OH) is used in organic light‐emitting diodes (OLEDs) as an electron transport material and emitting layer. By means of ab initio Hartree–Fock (HF) and density functional theory (DFT) B3LYP methods, the structure of AlMq₂OH was optimized. The frontier molecular orbital characteristics and energy levels of AlMq₂OH have been analyzed systematically to study the electronic transition mechanism in AlMq₂OH. For comparison and calibration, bis(8‐quinolinolato)aluminum(III) hydroxide complex (Alq₂OH) has also been examined with these methods using the same basis sets. The lowest singlet excited state (S₁) of AlMq₂OH has been studied by the singles configuration interaction (CIS) method and time‐dependent DFT (TD‐DFT) using a hybrid functional, B3‐LYP, and the 6‐31G* basis set. The lowest singlet electronic transition (S₀ → S₁) of AlMq₂OH is π → π* electronic transitions and primarily localized on the different quinolate ligands. The emission of AlMq₂OH is due to the electron transitions from a phenoxide donor to a pyridyl acceptor from another quinolate ligand including C → C and O → N transference. Two possible electron transfer pathways are presented, one by carbon, oxygen, and nitrogen atoms and the other via metal cation Al³⁺. The comparison between the CIS‐optimized excited‐state structure with the HF ground‐state structure indicates that the geometric shift is mainly confined to the one quinolate and these changes can be easily understood in terms of the nodal patterns of the highest occupied and lowest unoccupied molecular orbitals. On the basis of the CIS‐optimized structure of the excited state, TD‐B3‐LYP calculations predict an emission wavelength of 499.78 nm. An absorption wavelength at 380.79 nm on the optimized structure of B3LYP/6‐31G* was predicted. They are comparable to AlMq2OH 485 and 390 nm observed experimentally for photoluminescence and UV‐vis absorption spectra of AlMq₂OH solid thin film on quartz, respectively. Lending theoretical corroboration to recent experimental observations and supposition, the reasons for the blue‐shift of AlMq2OH were revealed. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Auxiliary basis sets for density‐fitting second‐order Møller–Plesset perturbation theory: Weighted core‐valence correlation consistent basis sets for the 4d elements Y–Pd

Auxiliary basis sets (ABS) specifically matched to the cc‐pwCVnZ‐PP and aug‐cc‐pwCVnZ‐PP orbital basis sets (OBS) have been developed and optimized for the 4d elements Y‐Pd at the second‐order Møller‐Plesset perturbation theory level. Calculation of the core‐valence electron correlation energies for small to medium sized transition metal complexes demonstrates that the error due to the use of these new sets in density fitting is three to four orders of magnitude smaller than that due to the OBS incompleteness, and hence is considered negligible. Utilizing the ABSs in the resolution‐of‐the‐identity component of explicitly correlated calculations is also investigated, where it is shown that i‐type functions are important to produce well‐controlled errors in both integrals and correlation energy. Benchmarking at the explicitly correlated coupled cluster with single, double, and perturbative triple excitations level indicates impressive convergence with respect to basis set size for the spectroscopic constants of 4d monofluorides; explicitly correlated double‐ζ calculations produce results close to conventional quadruple‐ζ, and triple‐ζ is within chemical accuracy of the complete basis set limit. © 2013 Wiley Periodicals, Inc.     

Molecular structure, IR spectra of 2-mercaptobenzothiazole and 2-mercaptobenzoxazole by density functional theory and ab initio Hartree–Fock calculations

Quantum chemistry calculations have been performed using Gaussian03 program to compute optimized geometry, harmonic vibrational frequency along with intensities in IR and Raman spectra and atomic charges at RHF/6-31+G*, B3LYP/6-31+G* and B3LYP/6-31++G* levels for 2-mercaptobenzothiazole (MBT, C₇H₅NS₂) and 2-mercaptobenzoxazole (MBO, C₇H₅NOS) in the ground state. The scaled harmonic vibrational frequencies have been compared with experimental FT-IR and FT-Raman spectra. The results show that the scaled theoretical vibrational frequencies is very good agreement with the experimental values. A detailed interpretation of the infrared and Raman spectra of 2-mercaptobenzothiazole and 2-mercaptobenzoxazole was reported. Comparison of calculated spectra with the experimental spectra provides important information about the ability of the computational method to describe the vibrational modes.     

DFT and ab initio direct dynamics studies on the hydrogen abstraction reactions of chlorine atoms with CH(4-n)F(n) (n = 1-3)

A direct dynamics study is carried out for the hydrogen abstraction reactions Cl + CH(4-n)F(n) (n = 1-3) in the temperature range of 200-1,000 K. The minimum energy paths (MEPs) of these reactions are calculated at the BH&H-LYP/6-311G(d,p) level, and the energies along the MEPs are further refined at the QCISD(T)/6-311+G(2df,2p) and QCISD(T)/6-311+G(d,p) (single-point) level. The rate constants obtained by using the improved canonical variational transition state theory incorporating small-curvature tunneling correction (ICVT/SCT) are in good agreement with the available experimental results. It is shown that the vibrational adiabatic potential energy curves for these reactions have two barriers, a situation similar to the analogous reactions CH(3)X+Cl (X=Cl, Br). The theoretical results show that for the title reactions the variational effect should not be neglected over the whole considered temperature range, while the small-curvature tunneling effect is only important in the lower temperature range. The effects of fluorine substitution on the rate of this kind of reactions are also examined.