This review covers a family of atomistic, mostly quantum chemistry (QC) based semiempirical methods for the fast and reasonably accurate description of large molecules in gas and condensed phase. The theory is derived from a density functional (DFT) perturbation expansion of the electron density in fluctuation terms to various orders similar to the original density functional tight binding model. The term “eXtended” in their name (xTB) emphasizes the parameter availability for almost the entire periodic table of elements (Z ≤ 86) and improvements of the underlying theory regarding, for example, the atomic orbital basis set, the level of multipole approximation and the treatment of the important electrostatic and dispersion interactions. A common feature of most members is their consistent parameterization on accurate gas phase theoretical reference data for geometries, vibrational frequencies and noncovalent interactions, which are the primary properties of interest in typical applications to systems composed of up to a few thousand atoms. Further specialized versions were developed for the description of electronic spectra and corresponding response properties. Besides a provided common theoretical background with some important implementation details in the efficient and free xtb program, various benchmarks for structural and thermochemical properties including (transition‐)metal systems are discussed. The review is completed by recent extensions of the model to the force‐field (FF) level as well as its application to solids under periodic boundary conditions. The general applicability together with the excellent cost‐accuracy ratio and the high robustness make the xTB family of methods very attractive for various fields of computer‐aided chemical research. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Electronic Structure Theory > Semiempirical Electronic Structure Methods Software > Quantum Chemistry

We present a composite procedure for the quantum‐chemical computation of spin–spin‐coupled 1H NMR spectra for general, flexible molecules in solution that is based on four main steps, namely conformer/rotamer ensemble (CRE) generation by the fast tight‐binding method GFN‐xTB and a newly developed search algorithm, computation of the relative free energies and NMR parameters, and solving the spin Hamiltonian. In this way the NMR‐specific nuclear permutation problem is solved, and the correct spin symmetries are obtained. Energies, shielding constants, and spin–spin couplings are computed at state‐of‐the‐art DFT levels with continuum solvation. A few (in)organic and transition‐metal complexes are presented, and very good, unprecedented agreement between the theoretical and experimental spectra was achieved. The approach is routinely applicable to systems with up to 100–150 atoms and may open new avenues for the detailed (conformational) structure elucidation of, for example, natural products or drug molecules.

Special-purpose classical force fields (FFs) provide good accuracy at very low computational cost, but their application is limited to systems for which potential energy functions are available. This excludes most metal-containing proteins or those containing cofactors. In contrast, the GFN2-xTB semiempirical quantum chemical method is parametrized for almost the entire periodic table. The accuracy of GFN2-xTB is assessed for protein structures with respect to experimental X-ray data. Furthermore, the results are compared with those of two special-purpose FFs, HF-3c, PM6-D3H4X, and PM7. The test sets include proteins without any prosthetic groups as well as metalloproteins. Crystal packing effects are examined for a set of smaller proteins to validate the molecular approach. For the proteins without prosthetic groups, the special purpose FF OPLS-2005 yields the smallest overall RMSD to the X-ray data but GFN2-xTB provides similarly good structures with even better bond-length distributions. For the metalloproteins with up to 5000 atoms, a good overall structural agreement is obtained with GFN2-xTB. The full geometry optimizations of protein structures with on average 1000 atoms in wall-times below 1 day establishes the GFN2-xTB method as a versatile tool for the computational treatment of various biomolecules with a good accuracy/computational cost ratio.

A fully quantum mechanical (QM) treatment to calculate electronic absorption (UV-vis) and circular dichroism (CD) spectra of typical biomolecules with thousands of atoms is presented. With our highly efficient sTDA-xTB method, spectra averaged along structures from molecular dynamics (MD) simulations can be computed in a reasonable time frame on standard desktop computers. This way, nonequilibrium structure and conformational, as well as purely quantum mechanical effects like charge-transfer or exciton-coupling, are included. Different from other contemporary approaches, the entire system is treated quantum mechanically and neither fragmentation nor system-specific adjustment is necessary. Among the systems considered are a large DNA fragment, oligopeptides, and even entire proteins in an implicit solvent. We propose the method in tandem with experimental spectroscopy or X-ray studies for the elucidation of complex (bio)molecular structures including metallo-proteins like myoglobin.

A series of formally triply bonded diplumbyne analogues of alkynes of the general formula ArPbPbAr (Ar = terphenyl ligand with different steric properties) was synthesized by two routes. All diplumbyne products were synthesized by a simple reduction of the corresponding Pb(II) halide precursor ArPb(Br) by DIBAL-H with yields in the range 8−48%. For one of the diplumbynes Ar Pri 4 PbPbAr Pri 4 (Ar Pri 4 = C 6 H 3 -2,6-(C 6 H 3 -2,6-Pr i 2 ) 2 ) it was shown that reduction of Ar Pri 4 Pb(Br) using a magnesium(I) beta-diketiminate afforded a much improved yield in comparison (29 vs 8%) to that obtained by reduction with DIBAL-H. The more sterically crowded diplumbyne Ar Pri 8 PbPbAr Pri 8 (Ar Pri 8 = C 6 H-3,5-Pr i 2 -2,6-(C 6 H 2 -2,4,6-Pr i 3 ) 2 ) displayed a shortened Pb−Pb bond with a length of 3.0382(5) Å and wide Pb−Pb−C angles of 114.73(7)°and 116.02( 6)°consistent with multiple-bond character with a bond order of up to 1.5. The others displayed longer metal−metal distances and narrower Pb−Pb−C angles that were consistent with a lower bond order that approached one. Computational studies of the diplumbynes yielded detailed insight of the unusual bonding and explained their similar electronic spectra arising from the flexibility of the C−Pb−Pb−C core in solution. Furthermore, the importance of London dispersion interactions for the stabilization of the diplumbynes was demonstrated.