Thomas Weymuth

New Benchmark Set of Transition-Metal Coordination Reactions for the Assessment of Density Functionals.

We present the WCCR10 data set of 10 ligand dissociation energies of large cationic transition metal complexes for the assessment of approximate exchange-correlation functionals. We analyze nine popular functionals, namely BP86, BP86-D3, B3LYP, B3LYP-D3, B97-D-D2, PBE, TPSS, PBE0, and TPSSh by mutual comparison and by comparison to experimental gas-phase data measured with well-known precision. The comparison of all calculated data reveals a large, system-dependent scattering of results with nonnegligible consequences for computational chemistry studies on transition metal compounds. Considering further the comparison with experimental results, the nonempirical functionals PBE and TPSS turn out to be among the best functionals for our reference data set. The deviation can be lowered further by including Hartree-Fock exchange. Accordingly, PBE0 and TPSSh are the two most accurate functionals for our test set, but also these functionals exhibit deviations from experimental results by up to 50 kJ mol(-1) for individual reactions. As an important result, we found no functional to be reliable for all reactions. Furthermore, for some of the ligand dissociation energies studied in this work, invoking semiempirical dispersion corrections yields results which increase the deviation from experimental results. This deviation increases further if structure optimization including such dispersion corrections is performed, although the contrary should be the case, pointing to the need to develop the currently available dispersion corrections further. Finally, we compare our results to other benchmark studies and highlight that the performance assessed for different density functionals depends significantly on the reference molecule set chosen.

MOVIPAC: Vibrational spectroscopy with a robust meta‐program for massively parallel standard and inverse calculations

We present the software package MOVIPAC for calculations of vibrational spectra, namely infrared, Raman, and Raman Optical Activity (ROA) spectra, in a massively parallelized fashion. MOVIPAC unites the latest versions of the programs SNF and AKIRA alongside with a range of helpful add‐ons to analyze and interpret the data obtained in the calculations. With its efficient parallelization and meta‐program design, MOVIPAC focuses in particular on the calculation of vibrational spectra of very large molecules containing on the order of a hundred atoms. For this purpose, it also offers different subsystem approaches such as Mode‐ and Intensity‐Tracking to selectively calculate specific features of the full spectrum. Furthermore, an approximation to the entire spectrum can be obtained using the Cartesian Tensor Transfer Method. We illustrate these capabilities using the example of a large π‐helix consisting of 20 (S)‐alanine residues. In particular, we investigate the ROA spectrum of this structure and compare it to the spectra of α‐ and 310‐helical analogs.

A Local-Mode Model for Understanding the Dependence of the Extended Amide III Vibrations on Protein Secondary Structure

The extended amide III region in vibrational spectra of polypeptides and proteins is particularly sensitive to changes in secondary structure. To investigate this structural sensitivity, we have performed density-functional calculations on the small model compound N-acetyl-l-alanine-N-methylamide, which are analyzed using the recently developed analysis in terms of localized modes [J. Chem. Phys. 2009, 130, 084106]. We find that the local modes obtained for different backbone conformations are actually rather similar. To probe the secondary structure sensitivity, we investigate the dependence of the local-mode frequencies and coupling constants on the torsional angles phi and psi. This enables us to set up a local-mode model of the extended amide III region for better understanding its structural sensitivity.

Inverse quantum chemistry: Concepts and strategies for rational compound design

The rational design of molecules and materials is becoming more and more important. With the advent of powerful computer systems and sophisticated algorithms, quantum chemistry plays a decisive role in the design process. While traditional quantum chemical approaches predict the properties of a predefined molecular structure, the goal of inverse quantum chemistry is to find a structure featuring one or more desired properties. Herein, we review inverse quantum chemical approaches proposed so far and discuss their advantages as well as their weaknesses.

Gradient-driven molecule construction: An inverse approach applied to the design of small-molecule fixating catalysts

Rational design of molecules and materials usually requires extensive screening of molecular structures for the desired property. The inverse approach to deduce a structure for a predefined property would be highly desirable, but is, unfortunately, not well defined. However, feasible strategies for such an inverse design process may be successfully developed for specific purposes. We discuss options for calculating “jacket” potentials that fulfill a predefined target requirement—a concept that we recently introduced (Weymuth and Reiher, MRS Proceedings 2013, 1524, DOI:10.1557/opl.2012.1764). We consider the case of small-molecule activating transition metal catalysts. As a target requirement we choose the vanishing geometry gradients on all atoms of a subsystem consisting of a metal center binding the small molecule to be activated. The jacket potential can be represented within a full quantum model or by a sequence of approximations of which a field of electrostatic point charges is the simplest. In a second step, the jacket potential needs to be replaced by a chemically viable chelate-ligand structure for which the geometry gradients on all of its atoms are also required to vanish. To analyze the feasibility of this approach, we dissect a known dinitrogen-fixating catalyst to study possible design strategies that must eventually produce the known catalyst.

Characteristic Raman Optical Activity Signatures of Protein β-Sheets

In this study, we compute and analyze theoretical Raman optical activity spectra of large model β-sheets in order to identify reliable signatures for this important secondary structure element. We first review signatures that have already been proposed to be indicative of β-sheets. From these signatures, we find that only the couplet in the amide I region can be regarded as a truly reliable signature. In addition, we propose a strong negative peak at ∼1350 cm(-1) to be another good signature for parallel as well as antiparallel β-sheets. We study the robustness of these signatures with respect to perturbations induced by the amino acid side chains, the overall conformation of the sheet structure, and microsolvation. It is found that the latter effects can be very well understood and separated employing the concept of localized modes. Finally, we investigate whether Raman optical activity is capable of discriminating between parallel and antiparallel β-sheets. The amide III region turns out to be most promising for this purpose.

How Many Chiral Centers Can Raman Optical Activity Spectroscopy Distinguish in a Molecule

To study the capabilities and limitations of Raman optical activity, (-)-(M)σ-[10]helicene and (-)-(M)σ-[4]helicene serve as scaffold molecules on which new chiral centers are introduced by substitution of hydrogen atoms with other functional groups. These functional groups are deuterium atoms, fluorine atoms, and methyl groups. Multiply deuterated species are compared. Then, results of singly deuterated derivatives are compared against results obtained from singly fluorinated and methylated derivatives. The analysis required the calculation of a total of 2433 Raman optical activity spectra. The method we propose for the comparison of the various Raman optical activity spectra is based on the total intensity of squared difference spectra. This allows a qualitative comparison of pairs of Raman optical activity spectra and the extraction of the pair of most similar Raman optical activity spectra for each group of stereoisomers. Different factors were accounted for, such as the spectral resolution (modeled by line broadening) and the range of vibrational frequencies considered. In the case of σ-[4]helicene all generated stereoisomers in each group can be distinguished from one another by Raman optical activity spectroscopy. For σ-[10]helicene this holds except for the lower one of the two resolutions considered. Here, the group consisting of stereoisomers with five chiral centers contains at least one pair of derivatives whose Raman optical activity spectra cannot be distinguished from one another. This indicates that an increased molecular size has a negative effect on the number of chiral centers which can be distinguished by Raman optical activity spectroscopy. Regarding the different substituents, stereoisomers are the better distinguishable in Raman optical activity spectroscopy, the more distinct the signals of the substituent are from the rest of the spectrum.

Identifying protein β-turns with vibrational Raman optical activity.

β-turns belong to the most important secondary structure elements in proteins. On the basis of density functional calculations, vibrational Raman optical activity signatures of different types of β-turns are established and compared as well as related to other signatures proposed in the literature earlier. Our findings indicate that there are much more characteristic ROA signals of β-turns than have been hitherto suggested. These suggested signatures are, however, found to be valid for the most important type of β-turns. Moreover, we compare the influence of different amino acid side chains on these signatures and investigate the discrimination of β-turns from other secondary structure elements, namely α- and 3(10)-helices.

Redox-Active Chiroptical Switching in Mono- and Bis-Iron Ethynylcarbo[6]helicenes Studied by Electronic and Vibrational Circular Dichroism and Resonance Raman Optical Activity

Introducing one or two alkynyl-iron moieties onto a carbo[6]helicene results in organometallic helicenes (2 a,b) that display strong chiroptical activity combined with efficient redox-triggered switching. The neutral and oxidized forms have been studied in detail by electronic and vibrational circular dichroism, as well as by Raman optical activity (ROA) spectroscopy. The experimental results were analyzed and spectra were assigned with the help of first-principles calculations. In particular, a recently developed method for ROA calculations under resonance conditions has been used to study the intricate resonance effects on the ROA spectrum of mono-iron ethynylhelicene 2 a.