Theoretical calculation and prediction for experimental design to obtain spin crossover complexes

DFT methods were utilized to study SCO complexes. [Fe(2btz)₂(NCX)₂] (2btz = 2,2′‐bithiazoline, X = S ( 1 ) and Se ( 2 )), [Fe(phen)₂(NCX)₂] (phen = 1,10‐phenantroline, X = S ( 3 ) and Se ( 4 )), and [Fe(bpy)₂(NCS)₂] ( 5 ) (bpy = 2,2′‐bipyridine) compounds, which have experimentally shown SCO behavior, were calculated. B3LYP, B3LYP*, OPBE, and OLYP with 6‐31G* and 6‐311 + G** basis sets were employed to calculate the ΔE_HS/LS energy gap as a clue to find complexes with SCO behavior. It is found that calculated result by B3LYP* with c₃ = 0.14 and OPBE methods and 6‐31G* basis set are in agreement with experimentally observed SCO complexes. Then, newly designed Fe(N‐N)₂(X)₂ complexes, where N‐N are bidentate nitrogen donor chelating ligands and X= SCN^‐, SeCN^‐, Cl^‐, Br^‐, I^‐, were chosen to see their potential to be SCO compounds. ΔE_HS/LS for potential SCO complexes are estimated from 0.8 to 6.5 kcal/mol in B3LYP* and 0.6–5.7 kcal/mol in OPBE. These calculations suggest [Fe(bpy)₂(NCSe)₂], [Fe(5dmbpy)₂(NCS)₂], and [Fe(3‐BrPhen)₂(NCSe)₂] compounds have the ability to show SCO behavior. © 2016 Wiley Periodicals, Inc.     

FeIII Quinolylsalicylaldimine Complexes: A Rare Mixed‐Spin‐State Complex and Abrupt Spin Crossover

A new synthesis of (8‐quinolyl)‐5‐methoxysalicylaldimine (Hqsal‐5‐OMe) is reported and its crystal structure is presented. Two Fe^III complexes, [Fe(qsal‐5‐OMe)₂]Cl ? solvent (solvent=2 MeOH ? 0.5 H₂O ( 1 ) and MeCN ? H₂O ( 2 )) have been prepared and their structural, electronic and magnetic properties studied. [Fe(qsal‐5‐OMe)₂] Cl ? 2 MeOH ? 0.5 H₂O ( 1 ) exhibits rare crystallographically independent high‐spin and low‐spin Fe^III centres at 150 K, whereas [Fe(qsal‐5‐OMe)₂]Cl ? MeCN ? H₂O ( 2 ) is low spin at 100 K. In both structures there are extensive π–π and C? H???π interactions. SQUID magnetometry of 2 reveals an unusual abrupt stepped‐spin crossover with T_1/2=245 K and 275 K for steps 1 and 2, respectively, with a slight hysteresis of 5 K in the first step and a plateau of 15 K between the steps. In contrast, 1 is found to undergo an abrupt half‐spin crossover also with a hysteresis of 10 K. The two compounds are the first Fe^III complexes of a substituted qsal ligand to exhibit abrupt spin crossover. These conclusions are supported by ⁵⁷Fe Mössbauer spectroscopy. Both complexes exhibit reversible reduction to Fe^II at ?0.18 V and irreversible oxidation of the coordinated qsal‐5‐OMe ligand at +1.10 V.     

Improving on equilibrium isotropic nuclear shielding constants

Isotropic nuclear shielding constants at the equilibrium molecular structure σ_eq and zero‐point vibrational corrections (ZPVCs) to σ_eq are evaluated using the B3LYP/aug‐cc‐pVTZ level of theory, as well as the KT2/aug‐cc‐pVTZ level of theory. Various scaling factors and systematic corrections are obtained by linear regression to experimental shielding constants. Comparisons of the scaled and systematically corrected equilibrium and vibrationally averaged shielding constants reveal that, at the 99% confidence level, the ZPVCs via second‐order perturbation theory do not improve the agreement of B3LYP/aug‐cc‐pVTZ and KT2/aug‐cc‐pVTZ calculated shielding constants with experiment. This holds true when the same analysis is applied to CCSD(T)/aug‐cc‐pCV[TQ]Z calculated σ_eq of Teale et al. [Journal of Chemical Physics 2013, 138, 024111]. In addition, at the 99% confidence level, B3LYP/aug‐cc‐pVTZ and KT2/aug‐cc‐pVTZ scaled and systematically corrected shielding constants are found to be statistically no different from CCSD(T)/aug‐cc‐pCV[TQ]Z calculated σ_eq. The use of scaling factors and systematic corrections could thus provide a cheaper but yet reasonably accurate alternative for the study of nuclear shielding constants of larger systems.     

Synthesis,structures and magnetic properties of nickel bis (dithiolene) complexes with [Fe(qsal)2]+

Two new compounds containing the possible Fe(III) spin-crossover cation, [Fe(qsal)₂]⁺ (qsalH = N-(8-quinolyl)salicylaldimine), and nickel bis(dithiolene) anions have been synthesized. Both are 1 : 1 salts [Fe(qsal)₂][Ni(dddt)₂] · CH₃CN · CH₃OH (1) and [Fe(qsal)₂][Ni(pddt)₂] (2) (dddt = 5,6-dihydro-1,4-dithiin-2,3-dithiolate; pddt = 6,7-dihydro-5H-1,4-dithiepin-2,3-dithiolate). They have been characterized by X-ray crystal structure determination, elemental analysis, UV-Vis spectra and magnetic susceptibility measurements. The UV–Vis spectra are dominated by [Ni(L)₂]^? (1, L = dddt; 2, L = pddt). Magnetic studies show antiferromagnetic interaction in 1 from intermolecular S···S contacts and π–π stacking interactions, while the antiferromagnetic interaction in 2 is very weak.     

Pyridylmethylamines as Ligands in Iron Halide Complexes – Coordination Behaviour Depending on the Halide,the Denticity of the Amino Ligand and the Oxidation State of Iron

2‐Pyridylmethylamine (amp) and 8‐aminochinoline (ach) readily form the following complexes with iron halides in methanol: [(amp)₂FeCl₂] ( 1a ), [(amp)₂FeBr₂] ( 1b ), [(ach)₂Fe(MeOH)₂]Br₂ ( 1c ), and [(amp)FeCl₂(μ‐OMe)]₂ ( 2 ). Methanol was chosen as a solvent because these reactions are rather complex in ether. For example, FeCl₃ forms the ionic complex pair [(dme)₂FeCl₂] [FeCl₄] ( 3 ) with 1,2‐dimethoxyethane (dme). The reaction of FeBr₂ with tridentate di(2‐pyridylmethyl)amine (dpa) and tetradentate 1,2‐dipyridyl‐1,2‐diaminoethane (dpdae) yields the complexes [(dpa)₂Fe]Br₂·2 MeOH ( 4 ) and [(dpdae)₂Fe] [FeBr₄] ( 5 ), respectively. Crystallographic and magnetochemical investigations show the high‐spin configuration for the complexes 1 and 2 , whereas the short Fe‐N distances of 4 clearly indicate a low‐spin state. Compound 2 exhibits an antiferromagnetic exchange interaction with a coupling constant J = ?29.4 cm^?1 (H;af = ?J S;af_A·S;af_B).     

A mechanism of the 1,3‐dipolar cycloaddition between the hydrogen nitryl HNO2 and acetylene HCCH: The electron localization function study on evolution of the chemical bonds

The reaction between the simplest nitro compound HNO₂ (hydrogen nitryl) and acetylene HCCH ‐ formally proceeding via 1,3‐dipolar cycloaddition ‐ has been studied by means of the B3LYP, MPW1K and MP2 methods. The energy barrier of 20.74 ÷ 32.91 kcal/mol is similar to ΔE_a of the NNO + HCCH process but is essentially larger than computed for the reactions of HCCH with fulminic acid (HCNO) and NNCH₂. Whole process is exothermic with the reaction energy: ?10.87 ÷ ?17.94 kcal/mol. An evolution of the chemical bonding has been analyzed by means of the Bonding Evolution Theory (BET) at the B3LYP/6‐31+G(d) and B3LYP/cc‐pVTZ levels. Two approximations of the reaction path have been considered, namely: the IRC and pseudo‐reaction paths. The reaction requires five steps and seven catastrophes of the fold and cusp type. A different effect of first fold catastrophe has been noticed. At the B3LYP/6‐31+G(d) level one of two nonbonding V_i=1,2(N) attractors is annihilated (F), meanwhile at B3LYP/cc‐pVTZ new V(N) attractor is created (F^?). The chemical bonds are not formed/broken in TS. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Computational NMR coupling constants: Shifting and scaling factors for evaluating 1JCH

Optimized shifting and/or scaling factors for calculating one‐bond carbon–hydrogen spin–spin coupling constants have been determined for 35 combinations of representative functionals (PBE, B3LYP, B3P86, B97‐2 and M06‐L) and basis sets (TZVP, HIII‐su3, EPR‐III, aug‐cc‐pVTZ‐J, ccJ‐pVDZ, ccJ‐pVTZ, ccJ‐pVQZ, pcJ‐2 and pcJ‐3) using 68 organic molecular systems with 88 ¹J_CH couplings including different types of hybridized carbon atoms. Density functional theory assessment for the determination of ¹J_CH coupling constants is examined, comparing the computed and experimental values. The use of shifting constants for obtaining the calculated coupling improves substantially the results, and most models become qualitatively similar. Thus, for the whole set of couplings and for all approaches excluding those using the M06 functional, the root‐mean‐square deviations lie between 4.7 and 16.4 Hz and are reduced to 4–6.5 Hz when shifting constants are considered. Alternatively, when a specific rovibrational contribution of 5 Hz is subtracted from the experimental values, good results are obtained with PBE, B3P86 and B97‐2 functionals in combination with HIII‐su3, aug‐cc‐pVTZ‐J and pcJ‐2 basis sets. Copyright © 2013 John Wiley & Sons, Ltd.     

A theoretical analysis of substituted aromatic compounds

The substituents ? CH₃, ? F, ? NO₂, ? OCH₃, and ? CH₂?CH₂ were placed at the ortho, meta, and para positions on the aromatic molecules aniline, benzaldehdye, nitrobenzene, and phenol. MMFF94, AM1, B3LYP, M06, M06‐2X, ωB97X, ωB97X‐d, and RI‐MP2 using cc‐pVDZ and cc‐pVTZ and CCSD(T) with cc‐pVDZ basis sets were used to calculate the geometries and energies of all regiomers of the molecules. Relative energies of the ortho and meta regiomers relative to the para regiomers were calculated and compared to the CCSD(T) values. A good basis set correlation between cc‐pVDZ and cc‐pVTZ was observed in RI‐MP2. Overall, RI‐MP2 gave the best correlation with the CCSD(T) results. All of the hybrid functionals showed similar accuracy and could effectively describe the intramolecular hydrogen‐bonding interactions of these compounds. The methoxy group at the para position in methoxyaniline, methoxyphenol, methoxynitrobenzene, and methoxybenzaldehyde was rotated around the phenyl‐O bond. HF, along with the cc‐pVDZ basis with the other methods, generated inaccurate energy profiles for p‐methoxyphenol. For the density functional theory methods, it was necessary to use improved grids to get smooth curves. © 2013 Wiley Periodicals, Inc.     

ChemInform Abstract: Smallest Deltahedra Silicon Dicarbide: C2Si32‐.

A total of 20 singlet and 36 triplet C₂Si₃^2‐ isomers are obtained by quantum chemical calculations at the B3LYP/aug‐cc‐pVTZ level.     

DFT calculations on the spin-crossover complex Fe(salen)(NO): a quest for the best functional

DFT calculations on the spin-crossover complex Fe(salen)(NO) provide a striking illustration of the comparative performance of different exchange-correlation functionals vis-à-vis the issue of transition metal spin state energetics. Thus, although the "classic" pure functionals PW91 and BLYP favor the S = 1/2 state by about 10 kcal/mol, relative to the S = 3/2 state, the hybrid functional B3LYP favors the latter state by nearly the same margin. In contrast, the newer pure functionals OLYP and OPBE, based on the OPTX exchange functional, as well as the B3LYP* hybrid functional (which has 15% Hartree-Fock exchange, compared with 20% for B3LYP) predict nearly isoenergetic S = 1/2 and 3/2 states, as required for a spin-crossover complex. Intriguingly, the OLYP and B3LYP* spin density profiles for the S = 1/2 state of Fe(salen)(NO) are substantially dissimilar.     

A computational study on the ring stretching modes of halogen‐substituted pyridine involved in H‐bonding

Density functional method B3LYP plus the AUG‐cc‐pVDZ and AUG‐cc‐pVTZ basis sets is used to investigate ring normal modes of halogen‐substituted pyridines involved in the N ··· H? X H‐bonds with HX (X = F, Cl). The results demonstrated that the formation of hydrogen bond leads to an increase in the frequencies of the ring breathing mode v₁, the N‐para‐C stretching mode v_6a and the meta‐CC stretching mode v_8a, whereas there is no change in the triangle mode v₁₂ for free pyridine and a smaller blue shift for substituted pyridines. There is a strong coupling between the C? Y stretching vibration and the triangle mode (ortho‐ and para‐substituted) or the breathing mode (meta‐substituted) in substituted pyridines, which leads to the frequency decrease in the triangle or breathing modes. The natural bond orbital analysis suggests that electrostatic interaction and charge transfer caused by the intermolecular and intramolecular hyperconjugations are the origin of the frequency blue shift in the ring stretching modes. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Calibration of DFT Functionals for the Prediction of Fe Mössbauer Spectral Parameters in Iron-Nitrosyl and Iron-Sulfur Complexes: Accurate Geometries Prove Essential

Six popular density functionals in conjunction with the conductor-like screening (COSMO) solvation model have been used to obtain linear M?ssbauer isomer shift (IS) and quadrupole splitting (QS) parameters for a test set of 20 complexes (with 24 sites) comprised of nonheme nitrosyls (Fe-NO) and non-nitrosyl (Fe-S) complexes. For the first time in an IS analysis, the Fe electron density was calculated both directly at the nucleus, ρ(0)(N), which is the typical procedure, and on a small sphere surrounding the nucleus, ρ(0)(S), which is the new standard algorithm implemented in the ADF software package. We find that both methods yield (near) identical slopes from each linear regression analysis but are shifted with respect to ρ(0) along the x-axis. Therefore, the calculation of the Fe electron density with either method gives calibration fits with equal predictive value. Calibration parameters obtained from the complete test set for OLYP, OPBE, PW91, and BP86 yield correlation coefficients (r(2)) of approximately 0.90, indicating that the calibration fit is of good quality. However, fits obtained from B3LYP and B3LYP* with both Slater-type and Gaussian-type orbitals are generally found to be of poorer quality. For several of the complexes examined in this study, we find that B3LYP and B3LYP* give geometries that possess significantly larger deviations from the experimental structures than OLYP, OPBE, PW91 or BP86. This phenomenon is particularly true for the di- and tetranuclear Fe complexes examined in this study. Previous M?ssbauer calibration fit studies using these functionals have usually included mononuclear Fe complexes alone, where these discrepancies are less pronounced. An examination of spin expectation values reveals B3LYP and B3LYP* approach the weak-coupling limit more closely than the GGA exchange-correlation functionals. The high degree of variability in our calculated S(2) values for the Fe-NO complexes highlights their challenging electronic structure. Significant improvements to the isomer shift calibrations are obtained for B3LYP and B3LYP* when geometries obtained with the OLYP functional are used. In addition, greatly improved performance of these functionals is found if the complete test set is grouped separately into Fe-NO and Fe-S complexes. Calibration fits including only Fe-NO complexes are found to be excellent, while those containing the non-nitrosyl Fe-S complexes alone are found to demonstrate less accurate correlations. Similar trends are also found with OLYP, OPBE, PW91, and BP86. Correlations between experimental and calculated QSs were also investigated. Generally, universal and separate Fe-NO and Fe-S fit parameters obtained to determine QSs are found to be of good to excellent quality for every density functional examined, especially if [Fe(4)(NO)(4)(μ(3)-S)(4)](-) is removed from the test set.     

Synthesis and Characterization of the Trisulfate Pb[S3O10] and Theoretical Analysis of H2S3O10

The reaction of (NH₄)₂PbCl₆ and fuming sulfuric acid (65 % SO₃) in a sealed glass tube at 250 °C led to colorless single crystals of Pb[S₃O₁₀] (orthorhombic, Pbcn, Z = 4, a = 10.9908(4), b = 8.5549(3), c = 8.0130(3) Å, V = 753.42(5) ų). The compound shows a three‐dimensional linkage of the tenfold oxygen coordinated Pb²⁺ ions and exhibits the unusual trisulfate anion, [S₃O₁₀]^2–, that consists of three vertex connected [SO₄] tetrahedra. The distances S–O within the S–O–S bridges of the anion are remarkable asymmetric with distances of 155 and 169 pm, respectively. This structural feature is well reproduced by calculations on a PBE0/cc‐pVTZ and a MP2/cc‐pVTZ level of theory. Similar calculations allow also for an inspection of the yet unknown corresponding acid, H₂S₃O₁₀. Also for this acid non‐symmetric S–O–S bridges are predicted. The thermal behavior of Pb[S₃O₁₀] is characterized by the loss of two equivalents of SO₃ at low temperature and the decomposition of intermediate Pb[SO₄] at higher temperature.     

The gas phase hydrogen‐bonded dimers of HOCl: A high‐level quantum chemical study

The O···H? O and Cl···H? O hydrogen bonding interactions were analyzed for HOCl dimers by using B3LYP, MP2, CCSD, and MP4(SDTQ) methods in conjunction with the various basis sets. Five isomers were found for the HOCl dimer. The ZPE and BSSE corrected binding energies were computed at the different levels of theory. At the optimized geometries obtained at CCSD/AUG‐cc‐pVDZ level, energies were re‐evaluated at MP4(SDTQ)/AUG‐cc‐pVTZ and CCSD(T)/cc‐pVTZ levels of theory. We found an average of ?20.9 and ?9.6 kJ/mol for the strength of the O···H and Cl···H hydrogen bonding interactions, respectively. Excitation and vertical ionization energies as well as rotational constants were computed at different levels of theory. The quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analysis were used to elucidate the nature of the interactions of HOCl dimers. The interaction energies were decomposed by Morokuma methodology. We have computed Δ_fH°(HOCl) and Δ_fH°(HOCl⁺) using the atomization reactions. The Δ_fH°₂₉₈(HOCl) values are ?17.85 and ?18.05 kcal/mol by using CBS‐Q and CBS‐QB3 extrapolation models, respectively, in good agreement with the results given in JANAF tables. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical investigation on the IR spectra of four‐membered silicon oxide ring

To test the mass effect on the ring vibrational frequencies (SiO₂) shift of the four‐membered silicon oxide ring, the deuterium and the tritium substituted cyclodisiloxanes on the hydrogen positions are examined at the CCSD(T)/cc‐pVTZ level of theory. The SiO₂ ring vibrations for a silicon oxide surface model compound, substituted cyclodisiloxane (Si₃O₅? O₂? Si₃O₅) with two six‐membered ring, are also calculated at the B3LYP/cc‐pVTZ level of theory. Our results of 909 and 920 cm^?1 are in good agreement with the experimental result. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

IR Low‐Temperature Matrix and ab Initio Study on β‐Alanine Conformers

For the first time the argon‐matrix low‐temperature FTIR spectra of β‐alanine are recorded. They reveal a quite complicated spectral pattern which suggests the presence of several β‐alanine conformers in the matrix. To interpret the spectra, the eighteen β‐alanine conformers, stable in the gas phase, are estimated at the B3LYP and MP2 levels combined with the aug‐cc‐pVDZ. Ten low‐energy structures are reoptimized at the QCISD/aug‐cc‐pVDZ and B3LYP and MP2 levels by using the aug‐cc‐pVTZ basis sets. Assignment of the experimental spectra is undertaken on the basis of the calculated B3LYP/aug‐cc‐pVDZ anharmonic IR frequencies as well as careful estimation of the conformer population. The presence of at least three β‐alanine conformers is demonstrated. The detailed analysis of IR spectra points to the possible presence of five additional β‐alanine conformers.     

An electronic structure theory investigation of the physical chemistry of the intermolecular complexes of cyclopropenylidene with hydrogen halides

The proton accepting and donating abilities of cyclopropenylidene (c‐C₃H₂) on its complexation with hydrogen halides H? X (X = F, Cl, Br) are analyzed using density‐functional theory with three functionals (PBE0, B3LYP, and B3LYP‐D) and benchmarked against second‐order Møller–Plesset (MP2) theory. Standard signatures including, inter alia, dipole moment enhancement, charge transfer from the carbenic lone pair to the antibonding σ*(H? X) orbital, and H? X bond elongation are examined to ascertain the presence of hydrogen bonding in these complexes. The latter property is found to be accompanied with a pronounced red shift in the bond stretching frequency and with a substantial increase in the infrared intensity of the band on complex formation. The MP2/aug‐cc‐pVTZ c‐C₃H₂···H? F complex potential energy surface turns out to be an asymmetric deep single well, while asymmetric double wells are found for the c‐C₃H₂···H? Cl and c‐C₃H₂···H? Br complexes, with an energy barrier of 4.1 kcal mol^?1 for proton transfer along the hydrogen bond in the latter complex. Hydrogen‐bond energy decomposition, with the reduced variational space self‐consistent field approach, indicates that there are large polarization and charge‐transfer interactions between the interacting partners in c‐C₃H₂···H? Br compared to the other two complexes. The C···H bonds are found to be predominantly ionic with partial covalent character, unveiled by the quantum theory of atoms in molecules. The present results reveal that the c‐C₃H₂ carbene divalent carbon can act as a proton acceptor and is responsible for the formation of hydrogen bonds in the complexes investigated. © 2012 Wiley Periodicals, Inc.     

Three New Iron(II) Thiocyanato Coordination Polymers Based on 4,4′‐Bipyridine as Ligand and the Influence of Methanol on Their Structures

Reaction of iron(II) thiocyanate with 4,4‐bipyridine (bipy) in methanol leads to the formation of three new solvates of different composition depending on the reaction conditions: At room temperature two new ligand‐rich 1:2 (1:2 = ratio between metal and N‐donor ligand) polymorphic forms [Fe(NCS)₂(bipy)₂ · 2MeOH]_n ( 1I ) and [Fe(NCS)₂(bipy)(MeOH)₂ · (bipy)]_n ( 1II ) are obtained, whereas solvothermal conditions leads to the formation of the new ligand‐deficient 1:1 compound [{Fe(NCS)₂(bipy)(MeOH)}₂]_n ( 2 ). All crystal structures were determined by X‐ray single crystal structure analysis. In the crystal structure of modification 1I the metal atoms are coordinated by four bridging bipy ligands, which connect them into layers. The methanol molecules occupy voids in the structure. Compared to 1I in modification 1II the crystal structure contains of linear Fe–bipy–Fe chains, which are further connected by hydrogen bonds between coordinating MeOH and noncoordinated bipy ligands into layers. The ligand‐deficient 1:1 compound 2 shows a completely different coordination topology with linear Fe–bipy–Fe chains, which are connected by coordinating methanol molecules into double‐chains. In all compounds the thiocyanato anions are terminal N‐bonded to the metal atoms. Investigation of the thermal behavior of compound 1I shows a two‐step decomposition, in which ligand‐deficient intermediates are formed. Magnetic measurements on 1I reveal Curie–Weiss paramagnetism with increasing antiferromagnetic interactions on cooling.     

Theoretical study on HBC− and HCB− anions using multiconfiguration second‐order perturbation theory

The HBC^? and HCB^? anions have been studied using the complete active space self‐consistent field, CASPT2, and DFT/B3LYP methods with the contracted atomic natural orbital (ANO) and cc‐pVTZ basis sets. The geometries of all stationary points along the potential energy surfaces were optimized at the CASSCF/ANO and B3LYP/cc‐pVTZ levels. The ground state of HBC^? is predicted to be X²? state, which is different from the previously published results. The CASPT2/ANO potential energy curves (PECs) of isomerization reactions were calculated as a function of HBC bond angle and the PECs also show the ²? state is the ground state of HBC^? anion. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009