Benchmarking semiempirical and DFT methods for the interaction of thiophene and diethyl sulfide molecules with a
Alexander V Vorontsov1, Panagiotis G Smirniotis2
1N.N. Vorozhtsov Novosibirsk Institute of Organic Chemistry, Pr. Ak. Lavrentyeva 9, Novosibirsk, 630090, Russia. alex.voronts@gmail.com.
Abstract:
Semiempirical methods pm6 and pm7 as well as density functional theory functionals exchange LSDA, exchange-correlation PW91 and PBE, hybrid B3LYP1 and PBE0 were compared for energy and geometry of thiophene, diethyl sulfide (DES) molecules and their binding to a frozen Ti(OH)4(H2O) complex having one coordinatively unsaturated Ti5C site representing small fragment of TiO2 anatase (001) surface. PBE0/6-31G(d) with DFT-D3 dispersion correction was the best method for description of thiophene and DES molecules geometries as comparison with experimental data demonstrated. Semiempirical methods pm6 and pm7 resulted in only three of four possible binding configurations of thiophene with the Ti(OH)4(H2O) complex while pm7 described correctly the enthalpy and all configurations of DES binding with the Ti(OH)4(H2O) complex. SBKJC pseudopotential and LSDA with and without dispersion correction produced flawed results for many configurations. PBE0 and PBE with and without dispersion correction and PW91 with 6-31G(d) basis set systematically produced dependable results for thiophene and DES binding to the Ti(OH)4(H2O) complex. PBE0-D3/6-31G(d), B3LYP1-D3/6-31G(d), and PBE-D3/6-31G(d) gave best match of binding energy for thiophene while PBE0/6-31G(d) gave best match of DES binding energy as comparison with CCSD(T) energy demonstrated. On the basis of the superior results obtained with PBE0/6-31G(d), it is the recommended method for modeling of adsorption over TiO2 surfaces. Such a conclusion is in agreement with recent literature.
Related Concept Videos
Crystal Field Theory - Tetrahedral and Square Planar Complexes
Crystal field theory (CFT) is applicable to molecules in geometries other than octahedral. In octahedral complexes, the lobes of the dx2−y2 and dz2 orbitals point directly at the ligands. For tetrahedral complexes, the d orbitals remain in place, but with only four ligands located between the axes. None of the orbitals points directly at the tetrahedral ligands. However, the dx2−y2 and dz2 orbitals (along the Cartesian axes) overlap with the ligands less than the dxy,...
Crystal Field Theory - Octahedral Complexes
To explain the observed behavior of transition metal complexes (such as colors), a model involving electrostatic interactions between the electrons from the ligands and the electrons in the unhybridized d orbitals of the central metal atom has been developed. This electrostatic model is crystal field theory (CFT). It helps to understand, interpret, and predict the colors, magnetic behavior, and some structures of coordination compounds of transition metals.
CFT focuses on...
Hybridization of Atomic Orbitals II
Preparation and Reactions of Sulfides
Predicting Molecular Geometry
Preparation and Reactions of Thiols


