Quantum-Chemical Investigation of Structure and Vibration Spectrum of [C(NPCl<Subscript>3</Subscript>)<Subscript>3</Subscript>]<Superscript>+</Superscript> Cation

Structural parameters of the cation [C(NPCl₃)₃]⁺ in vacuum and in acetonitrile are calculated by the methods RHF/6-311G(3d₆), RHF/6-31(3d₅) and B3LYP/6-311(3d₅f₇). Formation energy of the free adduct (MeCN) 2[C(NPCl₃)₃]⁺ is calculated and nonspecific character of interaction of the cation with liquid acetonitrile is established. Vibration spectrum of the cation is calculated and theoretical interpretation of IR and Raman spectra of the salt [C(NPCl₃)₃]⁺[SbCl₆]^? is refined.     

Theoretical investigation on multinuclear NMR chemical shifts of some tris(trifluoromethyl)boron complexes

Tris(trifluoromethyl)boron complexes have unusual properties and may find applications in many fields of chemistry, biology, and physics. To gain insight into their NMR properties, the isotropic ¹¹B, ¹³C, and ¹⁹F NMR chemical shifts of a series of tris(trifluoromethyl)boron complexes were systematically studied using the gauge‐included atomic orbitals (GIAO) method at the levels of B3LYP/6‐31 + G(d,p)//B3LYP/6‐31G* and B3LYP/6‐311 + G(d,p)//B3LYP/6‐311 + G(d,p). Solvent effects were taken into account by polarizable continuum models (PCM). The calculated results were compared with the experimental values. The reason that the structurally inequivalent fluorine atoms in a specific species give a same chemical shift in experimental measurements is attributed to the fast rotation of CF₃ group around the B? C(F₃) bond because of the low energy barrier. The calculated ¹¹B, ¹³C(F₃), and ¹⁹F chemical shifts are in good agreement with the experimental measurements, while the deviations of calculated ¹³C(X, X = O, N) chemical shifts are slightly large. For the latter, the average absolute deviations of the results from B3LYP/6‐311 + G(d,p)//B3LYP/6‐311 + G(d,p) are smaller than those from B3LYP/6‐31 + G(d,p)//B3LYP/6‐31G*, and the inclusion of PCM reduces the deviation values. The calculated ¹⁹F and ¹¹B chemical shieldings of (CF₃)₃BCO are greatly dependent on the optimized structures, while the influence of structural parameters on the calculated ¹³C chemical shieldings is minor. Copyright © 2009 John Wiley & Sons, Ltd.     

Theoretical study on the insertion reactions of the germylenoid H2GeClMgCl with RH (R = F, OH, NH2)

The insertion reactions of the germylenoid H₂GeClMgCl with RH (R = F, OH, NH₂) have been studied by using the DFT B3LYP and QCISD methods. The geometries of the stationary points on the potential energy surfaces of the reactions were optimized at the B3LYP/6-311+G(d, p) level of theory. The calculated results indicate that all the mechanisms of the three insertion reactions are identical to each other. The QCISD/6-311++G(d, p)//B3LYP/6-311+G(d, p) calculated potential energy barriers for the three insertion reactions of R = F, OH, and NH₂ are 164.62, 193.30, and 200.73 kJ mol^?1, and the reaction energies for the three reactions are ?57.46, ?35.65, and ?22.22 kJ mol^?1, respectively. Under the same situation, the insertion reactions should occur easily in the following order H-F > H-OH > H-NH₂. In THF solvent the insertion reactions get more difficult than in gas phase.     

Theoretical studies on a possible synthesis reaction pathway on N8 (CS) clusters

The potential energy surfaces of N₈ clusters were investigated by density functional theory (DFT) and a possible synthesis reaction pathway for N₈ (C_S) was suggested. The species involved were fully optimized up to the B3LYP/6‐311+G* level of theory. Relative energies were further calculated at the QCISD/6‐311+G*//B3LYP/6‐311+G* level. The reaction rate constants of these steps from the 1 (N₅⁺?N₃^?, complex, C_S) to 2 (N₈, C_S), 2 (N₈, C_S) to 3 (N₈, C_S), 3 (N₈, C_S) to 4 (N₈, D_2d), and 4 (N₈, D_2d) to 5 (N₈, C_S) reactions were predicted by the VTST theory. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1334–1339, 2001     

Theoretical investigations on nonlinear optical and spectroscopic properties of 6-(3,3,4,4,4-pentafluoro-2-hydroxy-1-butenyl)-2,4-pyrimidinedione: An efficient NLO material

In this study, quantum chemical calculations of geometric parameters, conformational, natural bond orbital (NBO) and nonlinear optical (NLO) properties, vibrational frequencies, ¹H and ¹³C NMR chemical shifts of the title molecule [C₉H₇F₅N₂O₃] in the ground state have been calculated with the help of Density Functional Theory (DFT-B3LYP/6-311++G(d,p)) and Hartree-Fock (HF/6-311++G(d,p)) methods. The optimized geometric parameters, vibrational frequencies, ¹H and ¹³C NMR chemical shifts values are compared with experimental values of the investigated molecules. Comparison between experimental and theoretical results showed that B3LYP/6-311++G(d,p) method is able to provide more satisfactory results. In order to understand this phenomenon in the context of molecular orbital picture, we examined the molecular frontier orbital energies (HOMO, HOMO-1, LUMO, and LUMO + 1), the energy difference (ΔE) between E _HOMO and E _LUMO, electronegativity (χ), hardness (η), softness (S) calculated by HF/6-311++G(d,p) and B3LYP/6-311++G(d,p) levels. The molecular surfaces, Mulliken, NBO, and Atomic polar tensor (APT) charges of the investigated molecule have also been calculated by using the same methods.     

HCN exchange on [Cu(HCN)4]+: a quantum chemical investigation

Density functional (B3LYP, B3PW91, X3LYP, BP86, PBEPBE, PW91PW91, and M06) and ab initio (MP2, MP4sdq, CCSD, and CCSD(T)) calculations with extended basis sets (6-311+G**, TZVP, LANL2DZ+p, and SDD+p, the latter including extra polarization and diffuse functions) indicate that HCN exchange on [Cu(HCN)₄]⁺ proceeds via an associative interchange (I_a) mechanism and a D_3h transition structure {[Cu(HCN)₅]⁺}^?. The activation barrier, relative to the model complex [Cu(HCN)₄]⁺·HCN, varies modestly, depending on the computational level. Typical values are 8.0?kcal?M^?1 (B3LYP/6-311+G**), 6.0?kcal?M^?1 (M06/6-311+G**), and 4.8?kcal?M^?1 (CCSD(T)/6-311+G**//MP2(full)/6-311+G**). Inclusion of an implicit solvent model (B3LYP(CPCM)/6-311+G**) leads to an activation barrier of 5.8?kcal?mol^?1. Comparison of the HCN exchange mechanisms on [Li(HCN)₄]⁺ (limiting associative, A) and [Cu(HCN)₄]⁺ (associative interchange, I_a) reveals that π back donation in the equatorial Cu–N bonds in the transition state determines the mechanism.     

Steric and Electronic Structure of C6H5XCF3 Molecules (X = O or S): A Quantum-Chemical Study

The potential functions of internal rotation around the Csp ²-X bond in C₆H₅XCF₃ molecules (X = O or S) were obtained by quantum-chemical calculations in the HF/6-31G(d), MP2(f)/6-31G(d), and B3LYP/6-31G(d) approximations. The calculations were performed in the range of torsion angles (angle between the planes of the benzene ring and Csp ²-X-Csp ³ bonds) from 0° to 90° with a 15° step. The barriers to rotation around the Csp ²-X bonds (kJ mol^{- 1}) were evaluated: for C₆H₅XCF₃, 7.60 (HF), 3.04 (MP2), and 1.04 (B3LYP); for C₆H₅XCF₃, 16.57 (HF), 14.67 (MP2), and 8.73 (B3LYP). The geometries (bond angles and bond lengths), Koopmans ionization potentials, and dipole moments of the molecules were calculated. The hybridization, energy, and population of the lone electron pairs of the heteroatoms, and also the energy of their resonance interaction with antibonding orbitals and the natural atomic charges were evaluated using the NBO approach.     

Theoretical studies on the reaction of pentafulvenone with water in gas phase and aqueous solvent

In gas phase, the hydrations of pentafulvenone to generate three types of cyclopentadienyl carboxylic acids are studied theoretically at the MP2/6-311+G**//B3LYP/6-311+G** level. A water molecule attacking the C=O double bond of pentafulvenone can yield cyclopentadienyl carboxylic acids via the formation of fulvenediols, and attacking the C=C double bond of pentafulvenone can directly yield cyclopentadienyl carboxylic acid. The barriers of rate-determining transition states are 42.2 and 30.4 kcal mol⁻¹, respectively. The barriers of rate-determining transition states for two water molecules system are 20.2 and 19.6 kcal mol⁻¹, respectively. The products can isomerize to each other. In aqueous solvent, the hydrations of pentafulvenone are investigated using PCM-UAHF model at the MP2 (PCM)/6-311+G**// B3LYP (PCM)/6-311+G** and MP2 (PCM)/6-311+G**// B3LYP/6-311+G** levels. The barriers of all rate-determining transition states are decreased. The added water molecule acts as catalyst in both gas phase and aqueous solvent. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Structural and electronic properties of heptachlor

In this work, the molecular geometry of heptachlor is investigated using ab initio HF, DFT, LDA, and GGA methods. The natural bond orbital (NBO) analysis is performed at the B3LYP/6-311++G(d,p) level of theory. The first order hyperpolarizability β_total, the mean polarizability Δα, the anisotropy of the polarizability Δα, and the dipole moment μ, are calculated by B3LYP/6-311++G(d,p) and HF/6- 311++G(d,p) methods. The first order hyperpolarizability (β_total) is calculated based on the finite field approach. UV spectral parameters along with HOMO, LUMO energies for heptachlor are determined in vacuum and the solvent phase using HF, DFT, and TD-DFT/B3LYP methods implemented with the 6-311++G(d,p) basis set. Atomic charges and electron density of heptachlor in vacuum and ethanol are calculated using DFT/B3LYP and TD-DFT/B3LYP methods and the 6-311++G(d,p) basis set. In addition, after the frontier molecular orbitals (FMOs), the molecular electrostatic potential (MEP), the electrostatic potential (ESP), the electron density (ED), and the solvent accessible surface of heptachlor are visualized as a results of the B3LYP/6-311++G(d,p) calculation. Densities of states (DOS), the external electric field (EF) effect on the HOMO-LUMO gap, and the dipole moment are investigated by LDA and GGA methods.     

Theoretical and experimental study of the absolute configuration of helical structure of (2R,3S)-Rubiginone A2 analog

A novel analogue of (2R,3S)-Rubiginone A₂ was synthesized as a chiral helical model compound via an eight-step procedure (2.7% overall yield). Quantum methods, such as density functional theory (DFT) at different basis sets of 6-311+(d), 6-311++G(2d,p), were used to compute its optical rotation and electronic circular dichroism at the B3LYP/6-311++G(2d,p) level in the gas phase and in solution using PCM model, respectively. UV corrections were performed in electronic circular dichroism (ECD) simulations to match the experimental ECD well. The suitable computational methods, e.g., B3LYP/6-311++G(2d,p)//B3LYP/6-311++G(2d,p) in the gas phase using zero-point energy in Boltzmann statistics, were found and suggested for optical rotation and circular dichroism computations that can be used for absolute configuration determination of chiral helical compounds.     

Ab initio study of dehalohydrogenation reaction of 2-halo-2,3-dihydrophosphinine

Decomposition of 2-fluoro-2,3-dihydrophosphinine (1), 2-chloro-2,3-dihydrophosphinine (3), 2-bromo-2,3-dihydrophosphinine (5) to phosphinine was investigated using Molecular orbital and density functional theory. Study on the B3LYP/6-311+G** level of theory revealed that the required energy for the decomposition of compounds 1, 3, and 5 to phosphinine is 30.56 kcal·mol^?1, 28.23 kcal·mol^?1, and 24.03 kcal·mol^?1, respectively. HF/6-311+G**//B3LYP/6-311+G** calculated barrier height for the decomposition of compound 1, 3, and 5 to phosphinine is 57.56 kcal·mol^?1, 37.26 kcal·mol^?1, and 30.77 kcal·mol^?1, respectively. Also, MP2/6-311+G**//B3LYP/6-311+G** results indicated that the barrier height for the decomposition of compound 1, 3, and 5 to phosphinine is 46.59 kcal·mol^?1, 47.28 kcal·mol^?1, and 42.57 kcal·mol^?1, respectively. Natural bond orbital (NBO) population analysis and nuclear independent chemical shift (NICS) results showed that, reactants are non-aromatic but products of elimination reaction are aromatic, C-H and C-X bonds are broken and H-X bond is appear.     

Octacoordinate carbon atom in tetra(metalloamino)methanes CN<Subscript>4</Subscript>M<Subscript>4</Subscript> (M = Be,Mg, Ca): Quantum-chemical investigation

Quantum-chemical methods, ab initio [MP2(full)/6-311+G**] and density functional theory (B3LYP/6-311+G**), were applied to investigation of the electronic and spatial structure of tetraaminomethane derivatives CN₄M₄ (M = Be, Mg, Ca) both in classic and nonclassic isomeric forms. For all compounds structurally stable nonclassic isomeric forms were found with an octacoordinate carbon atom. The CN₄Be₄ is energetically most stable in a classic structure with a tetracoordinate carbon atom, whereas CN₄Mg₄ and CN₄Ca₄ are more stable in nonclassic structures with an octacoordinate carbon center.     

Highly symmetrical phthalocyanines and perfluorophthalocyanines: The quantum-chemical study

By the quantum-chemical method (U)B3LYP/6-31G(d ₅,p) are determined point symmetry group (D _4h ) and equilibrium structure of phthalocyanine (PcH₂), phthalocyaninates PcBe, PcMg, PcCo, PcNi, PcCu, PcZn, perfluorophthalocyaninates ^FPcNi, ^FPcCu, ^FPcZn, cations Pc⁺Mt, ^FPc⁺Mt and anion PcCo⁻. In the approximation (U)B3LYP/6-311++G(3d ₅ f ₇,p)//6-31G(d ₅,p) is achieved the satisfactory accuracy of the calculation of ionization potentials of the studied molecules. Effect of nuclear relaxation at the ionization is 0.07±0.06 eV; correction for “zero” vibrations does not exceed 0.01 eV. Perfluorination increases ionization potentials by 0.7–0.8 eV.     

CH3SH和CH3SLi分子间的锂键和氢键共存型相互作用

在CH₃SLi＋CH₃SH势能面上求得锂键和氢键共存型复合物的两种稳定构型. 频率分析表明, 与单体相比复合物中S(5)—Li(6)键伸缩振动频率发生红移, 而C(8)—H(10)键伸缩振动频率发生蓝移. 经B3LYP/6-311＋＋G**, MP2/6-311＋＋G**及MP2/AUG-CC-PVDZ水平计算的含基组重叠误差(BSSE)校正的复合物?中相互作用能分别为－58.99, －57.87和－62.89 kJ•mol^－1. 采用自然键轨道(NBO)理论, 分析了复合物中单体轨道间的电荷转移, 电子密度重排及其与相关键键长变化的本质等. 采用分子中的原子(AIM)理论分析了复合物中氢键和锂键的电子密度拓扑性质．在极化连续模型(PCM)下, 考察了溶剂化效应. 结果表明, 所考察的水、二甲亚砜、乙醇和乙醚等四种溶剂均使单体间的相互作用能增大, 且溶剂对复合物中的锂键结构及其振动频率具有显著的影响, 而对复合物中的氢键的振动频率影响不大.     

Density functional theory insight into Eu(III) and Am(III) complexes with two 2,6-dicarboxypyridine diamide-type ligands

Extraction complexes of Eu(III) and Am(III) with two 2,6-dicarboxypyridine diamide-type ligands L–A and L–B (Fig. 1) are studied by density functional theory (DFT). At both B3LYP/6-31G(d)/RECP and MP2/6-31G(d)/RECP levels of theory, the geometrical optimizations of the structures of the complexes can achieve the same accuracy and obtain the same geometrical configuration. At the B3LYP/6-311G(d,p)/RECP level of theory Eu³⁺ and Am³⁺ prefer to form [ML]³⁺ complexes under the solvation conditions, and the Am(III) complexes with L–A are more stable than the corresponding Eu(III) complexes. In the system with the ligand L–B, both [ML]³⁺ and [ML(NO₃)₃] species are very unstable.     

Quantum-Chemical Investigation of Pentafluorophenylxenonium and (Heptafluorocyclohexa-1,4-dien-1-yl)xenonium Cations and Related Molecules

Quantum-chemical DFT B3LYP/(6-311G** &; Xe 3-111G*) package was used for calculation of structural parameters and NMR spectra of cations [MeCN ? XeC₆F₅]⁺, [XeC₆F₅]⁺ and [XeC₆F₇]⁺ free and in acetonitrile as solvent (solv., C-PCM). Structures and relative energies of hypothetical intermediates (perfluorocyclohexa-3-en-1-yliden)xenon(II) XeC₆F₈(solv) and (perfluorocyclohexa-3-en-1,1-diyl)xenon(0) Xe: C₆F₈(solv) in the reaction of cation [XeC₆F₇]⁺(solv) with anion F^?(solv), and those of perfluoro(3-methylenecyclopentene), the final compound formed at decomposition of [XeC₆F₇]⁺ with fluorides in acetonitrile, were determined.     

Quantumchemical investigation of borabenzene adduct with pyridine

Atomic charges and structural parameters of borabenzene, pyridine and their adduct in free state, in liquid argon and in tetrahydrofuran are calculated by the quantum-chemical method B3LYP/6-311G(3d5f7,p) &; PCM. Mutual polarization of the adduct and medium results in small increase in boron-nitrogen interatomic distance and dihedral angle between the aromatic heterocycles. Calculated dipole moment of the adduct (7.17 D) is by 4.55 D over than the sum of dipole moments of the free components. Internal rotation barriers are not high: 1 kcal mol^?1 (0° and 180°) and 4 kcal mol^?1 (90°). The borabenzene-pyridine bonding energy (46 kcal mol^?1) is higher than that with dinitrogen (19 kcal mol^?1) and xenon atom (6 kcal mol^?1). The B-N bond length in the little stable adduct with dinitrogen is by 0.08 Å shorter than in the stable adduct with pyridine. The Lewis acid properties inherent in borabenzene are transferred on the π-electron system of pyridine fragment in the adduct. The electron transfer wavelength from borabenzene to pyridine fragment in argon matrix is by 109 nm higher than in tetrahydrofuran, as calculated by CIS CNDO/S method.     

E/Z(C=C)-Isomerization of enamines of 3-formyl-4-hydroxycoumarin induced by organic solvents

According to ¹? and ¹³? NMR data, enamines of 3-formyl-4-hydroxycoumarin exist in the keto enamine tautomeric form and undergo Z/E-isomerization around the C=C bond in CDCl₃, DMSO-d₆, and CD₃OD at room temperature. The activation energies of ?/Z-isomerization were measured experimentally and calculated by the B3LYP/6-311++G(d,p) method. An X-ray diffraction study showed that 3-(benzyliminomethyl)chromane-2,4-dione in the crystalline state exists as a mixture of two keto enamine isomers.     

Geminal interactions in ClZ(CH3)2X molecules (Z = C, Si, Ge) according to ab initio calculations

The structure and electronic parameters of ClZ(CH₃)₂X molecules (Z = C, Si, Ge, X = CH₃, OCH₃) were calculated by the RHF/6–31G(d) and RHF/6–311G(d,p) methods with full geometry optimization; calculations of ClZ(CH₃)₂OCH₃ molecules were also performed by the RHF/6–31G(d) method with partial geometry optimization. The ³⁵Cl NQR frequencies calculated from the populations of less diffuse 3p constituents of valence p orbitals of chlorine [RHF/6–31G(d)] were in agreement with the experimental values. The ³⁵Cl NQR frequencies for molecules with X = OCH₃ are lower than those for molecules with X = CH₃ (the Z atom being the same), due mainly to direct through-field polarization of the Z-Cl bond, induced by the effect of unshared electron pair of the oxygen atom in the trans position with respect to that bond. The difference in the ³⁵Cl NQR frequencies decreases in going from Z = C to Z = Si, Ge, in parallel with variation of the Z-Cl bond polarization as the size of Z increases.     

Analysis of methods for calculating the transition frequencies of the torsional vibration of acrolein isomers in the ground (S 0) electronic state

B3LYP, MP2, CCSD(T), and MP4/MP2 in the 6-311G(d, p), 6-311++G(d, p), cc-pVTZ, aug-cc-pVTZ bases used to calculate the transition frequencies of torsional vibration of trans- and cis-isomers of acrolein in the ground electronic state (S ₀) are analyzed. It is found that for trans-isomers, all methods of calculation except for B3LYP in the cc-pVTZ basis yield good agreement between the calculated and experimental values. It is noted that for the cis-isomer of acrolein, no method of calculation confirms the experimental value of the frequency of torsional vibration (138 cm^?1). It is shown that the calculated and experimental values for obertones at 273.0 cm^?1 and other transitions of torsional vibration are different for this isomer in particular. However, it is established that in some calculation methods (B3LYP, MP2), the frequency of the torsional vibration of the cis-isomer coincides with another experimental value of this frequency (166.5 cm^?1). It is concluded that in analyzing the vibrational structure of the UV spectrum, the calculated and experimental values of its obertone (331.3 cm^?1) coincide, along with its frequency. It is also noted that the frequency of torsional vibration for the cis-isomer (166.5 cm^?1) can also be found in other experimental works if we change the allocation of torsional transition 18 ₁ ¹ .