TY - JOUR A1 - Boese, Adrian Daniel T1 - Assessment of coupled cluster theory and more approximate methods for hydrogen bonded systems JF - Journal of chemical theory and computation N2 - To assess the accuracy of post-Hartree-Fock methods like CCSD(T), MP3, MP2.5, MP2, SCS-MP2, SOS-MP2, and DFT-SAPT, we evaluated several effects going beyond valence-correlated CCSD(T). For 16 small hydrogen bonded systems, CCSD(T) achieves an RMS error of 0.17 kJ/mol in the dissociation energy compared to our best estimate, which is a composite method akin to W4 theory. The error of CCSD(T) is thus much lower than for atomization energies. MP2 is surprisingly accurate for these systems with an RMS error of 1.3 kJ/mol. MP2.5 yields a clear improvement over MP2 (RMS of 0.5 kJ/mol) but still has an error about 3 times as large as CCSD(T) for the absolute RMS and almost 10 times as large for the relative RMS. error. Neither SCS-MP2, SOS-MP2, nor DFT-SAPT yield lower errors than MP2. With a Delta CCSD(T) correction to MP2, the basis set limit is readily achieved when employing diffuse functions-without these, the convergence is rather slow. Y1 - 2013 U6 - https://doi.org/10.1021/ct400558w SN - 1549-9618 SN - 1549-9626 VL - 9 IS - 10 SP - 4403 EP - 4413 PB - American Chemical Society CY - Washington ER - TY - JOUR A1 - Codorniu-Hernandez, Edelsys A1 - Boese, Adrian Daniel A1 - Kusalik, Peter G. T1 - The hemibond as an alternative condensed phase structure for the hydroxyl radical JF - Canadian journal of chemistry = Revue canadienne de chimie N2 - Despite the critical importance of the hydroxyl radical in major scientific fields, there are still open questions on the behavior of this species in the aqueous phase. In particular, there has been much debate on the existence of a hemibonded interaction between the hydroxyl radical and water molecules. While some reports indicate that the hemibonded radical might explain some experimental data, others have claimed that this interaction is simply a density functional theory (DFT) artifact. Here, we provide results from high level (basis set limit of coupled-cluster levels up to single, double, triple excitations (CCSD(T)) and beyond) ab initio calculations of different OH center dot(H2O)(n) clusters in the gas phase to accurately explore the existence of the hemibonded interaction and its energy difference with respect to other well-defined hydrogen bond interactions. Additional comparisons with second order perturbation theory (MP2) and DFT are also presented. Constrained molecular dynamics was applied to determine the free energy for the formation/disruption and ice systems. Overall, our findings confirm that the hemibond can be an alternative structure for the hydroxyl radical in the condensed phase when the formation of hydrogen bonds is impeded. These results will aid the understanding of theoretical and experimental data and help future experimental designs for the detection of this important species. KW - hydroxyl radical KW - hemibond KW - Car-Parrinello KW - ab initio KW - density functional theory (DFT) Y1 - 2013 U6 - https://doi.org/10.1139/cjc-2012-0520 SN - 0008-4042 VL - 91 IS - 7 SP - 544 EP - 551 PB - NRC Research Press CY - Ottawa ER - TY - JOUR A1 - Boese, Adrian Daniel T1 - Density Functional Theory and Hydrogen Bonds: Are We There Yet? JF - ChemPhysChem : a European journal of chemical physics and physical chemistry N2 - Density functional theory (DFT) has become more successful at introducing dispersion interactions, and can be thus applied to a wide range of systems. Amongst these are systems that contain hydrogen bonds, which are extremely important for the biological regime. Here, the description of hydrogen-bonded interactions by DFT with and without dispersion corrections is investigated. For small complexes, for which electrostatics are the determining factor in the intermolecular interactions, the inclusion of dispersion with most functionals yields large errors. Only for larger systems, in which van der Waals interactions are more important, do dispersion corrections improve the performance of DFT for hydrogen-bonded systems. None of the studied functionals, including double hybrid functionals (with the exception of DSD-PBEP86 without dispersion corrections), are more accurate than MP2 for the investigated species. KW - ab initio calculations KW - basis sets KW - density functional calculations KW - hydrogen bonds KW - intermolecular interactions Y1 - 2015 U6 - https://doi.org/10.1002/cphc.201402786 SN - 1439-4235 SN - 1439-7641 VL - 16 IS - 5 SP - 978 EP - 985 PB - Wiley-VCH CY - Weinheim ER - TY - JOUR A1 - Tsendra, Oksana A1 - Scott, Andrea Michalkova A1 - Gorb, Leonid A1 - Boese, Adrian Daniel A1 - Hill, Frances C. A1 - Ilchenko, Mykola M. A1 - Leszczynska, Danuta A1 - Leszczynski, Jerzy T1 - Adsorption of Nitrogen-Containing Compounds on the (100) alpha-Quartz Surface: Ab Initio Cluster Approach JF - The journal of physical chemistry : C, Nanomaterials and interfaces N2 - A cluster approach extended to the ONIOM methodology has been applied using several density functionals and Moller-Plesset perturbation theory (MP2) to simulate the adsorption of selected nitrogen-containing compounds [NCCs, 2,4,6-trinitrotoluene (TNT), 2,4-dinitrotoluene (DNT), 2,4-dinitroanisole (DNAN), and 3-nitro-1,2,4-triazole-5-one (NTO)] on the hydroxyated (100) surface of a-quartz. The structural properties were calculated using the M06-2X functional and 6-31G(d,p) basis set. The M06-2X-D3, PBE-D3, and MP2 methods were used to calculate the adsorption energies. Results have been compared with the data from other studies of adsorption of compounds of similar nature on silica. Effect of deformation of the silica surface and adsorbates on the binding energy values was also studied. The atoms in molecules (AIM) analysis was employed to characterize the adsorbate-adsorbent binding and to calculate the bond energies. The silica surface shows different sorption affinity toward the chemicals considered depending on their electronic structure. All target NCCs are physisorbed on the modeled silica surface. Adsorption occurs due to the formation of multiple hydrogen bonds between the functional groups of NCCs and surface silanol groups. Parallel orientation of NCCs interacting with the silica surface was found to be favorable when compared with perpendicularly oriented NCCs. NTO was found to be the most strongly adsorbed on the silica surface among all of the considered compounds. Dispersion correction was shown to play an important role in the DFT calculations of the adsorption energies of silica-NCC systems. Y1 - 2014 U6 - https://doi.org/10.1021/jp406827h SN - 1932-7447 VL - 118 IS - 6 SP - 3023 EP - 3034 PB - American Chemical Society CY - Washington ER - TY - GEN A1 - Boese, Adrian Daniel T1 - Assessment of coupled cluster theory and more approximate methods for Hydrogen Bonded Systems (vol 9, pg 4403, 2013) T2 - Journal of chemical theory and computation Y1 - 2014 U6 - https://doi.org/10.1021/ct500041j SN - 1549-9618 SN - 1549-9626 VL - 10 IS - 2 SP - 893 EP - 893 PB - American Chemical Society CY - Washington ER - TY - JOUR A1 - Codorniu-Hernandez, Edelsys A1 - Hall, Kyle Wm. A1 - Boese, Adrian Daniel A1 - Ziemianowicz, Daniel A1 - Carpendale, Sheelagh A1 - Kusalik, Peter G. T1 - Mechanism of O(P-3) Formation from a Hydroxyl Radical Pair in Aqueous Solution JF - Journal of chemical theory and computation N2 - The reaction mechanism for the rapid formation of a triplet oxygen atom, O(P-3), from a pair of triplet-state hydroxyl radicals in liquid water is explored utilizing extensive Car-Parrinello MD simulations and advanced visualization techniques. The local solvation structures, the evolution of atomic charges, atomic separations, spin densities, electron localization functions, and frontier molecular orbitals, as well as free energy profiles, evidence that the reaction proceeds through a hybrid (hydrogen atom transfer and electron proton transfer) and hemibond-assisted reaction mechanism. A benchmarking study utilizing high-level ab initio calculations to examine the interactions of a hydroxyl radical pair in the gas phase and the influence of a hemibonded water is also provided. The results presented here should serve as a foundation for further experimental and theoretical studies aimed at better understanding the role and potential applications of the triplet oxygen atom as a potent reactive oxygen species. Y1 - 2015 U6 - https://doi.org/10.1021/acs.jctc.5b00783 SN - 1549-9618 SN - 1549-9626 VL - 11 IS - 10 SP - 4740 EP - 4748 PB - American Chemical Society CY - Washington ER - TY - JOUR A1 - Boese, Adrian Daniel T1 - Basis set limit coupled-cluster studies of hydrogen-bonded systems JF - Molecular physics N2 - As hydrogen-bonded systems are of utmost importance in especially biological and chemical systems, a new set of highly accurate reference dissociation energies, denoted HB49, is devised. For the molecules in this set, the basis set convergence of post-Hartree-Fock methods, including F12 methods, is investigated. Using combined Moller-Plesset perturbation theory (MP2) and CCSD(T) approaches for energies and MP2 and QCISD(T) for gradients, we achieve CCSD(T) accuracy, which has been determined before to yield an accuracy of 0.2 kJ/mol for a subset of HB49. Both conventional extrapolation techniques and F12 techniques are competitive with each other. By using MP2+Delta CCSD(T), a rather fast basis set convergence is obtained when both basis sets are carefully chosen. KW - coupled-cluster KW - hydrogen bonds KW - ab initio KW - F12 methods Y1 - 2015 U6 - https://doi.org/10.1080/00268976.2014.1001806 SN - 0026-8976 SN - 1362-3028 VL - 113 IS - 13-14 SP - 1618 EP - 1629 PB - Routledge, Taylor & Francis Group CY - Abingdon ER - TY - JOUR A1 - Boese, Adrian Daniel A1 - Boese, Roland T1 - Tetrahydrothiophene and Tetrahydrofuran, Computational and X-ray Studies in the Crystalline Phase JF - Crystal growth & design : integrating the fields of crystal engineering and crystal growth for the synthesis and applications of new materials N2 - Calculations at various levels of theory with different methods and respective evaluations confirm that the twist conformation (C-2) is preferred for tetrahydrothiophene (THT) in the gas phase. In the crystalline phase, achieved by a laser assisted crystallization device, THT has C-1 symmetry (slightly distorted C-2 symmetry) in the chiral space group P2(1)2(1)2(1). This is obviously a packing effect caused by the nonsymmetrical arrangement of neighboring molecules. The distortion from C-2 symmetry costs very little energy as confirmed by computational methods in the gas phase. Only one enantiomer of the chiral THT is found in the cell which requires spontaneous crystallization, which results in a racemic mixture of crystals, or a racemization occurs prior to/during nucleation or in the embryonic state. The racemization happens by a mechanism that can be described as a partial pseudo rotation within a five-membered mono-heterocycle with a C-2-C-S-C-2' transition (C-2 and C-2' are enantiomers) maintaining the heteroatom residing within the symmetry elements. While THT has the molecular symmetry of the gas phase almost also in the crystalline phase, THF has an envelope conformation (CS). This was also established by calculations at various levels of theory which agrees well with the previously experimentally found conformation by electron diffraction. However, in the X-ray crystal structure, previously determined by Luger & Buschmann, THF has C-2 symmetry in the centrosymmetric space group C2/c with the oxygen atom situated on the crystallographic C-2 polar axis, requesting a racemic crystal for the twisted conformers of the enantiomers. No solid-state phase transitions were detected within the experimental ranges for THT and THF. Following the stabilization by molecular clustering, and ending at the crystal lattice, we stepwise increased the number of molecules by calculation of the respective monomers, dimers, trimers, and tetramers for THF and THT. The starting point was taken from the arrangements as found in the respective crystal structures. Both conformational enantiomers are equal in energy. In such cases, a crystal may contain either a racemate of conformers or one of the conformational enantiomers only. The first case is observed in THF, the latter one in THT. It is quite likely that the selection of one enantiomeric conformer of THT from an equilibrium of conformers at the early stage of nucleation (embryonic stage) is responsible for the spontaneous crystallization. In order to check if THF could form a polymorph with the molecular packing of THT and vice versa, we first calculated THF and THT in their respective crystal lattices as determined by X-ray diffraction. Exchanging the compounds in the THT and THF crystal lattices (i.e., replacing O against S and vice versa) results in significantly worse lattice energies indicating that such a polymorph is not a probable option. Y1 - 2015 U6 - https://doi.org/10.1021/cg501228w SN - 1528-7483 SN - 1528-7505 VL - 15 IS - 3 SP - 1073 EP - 1081 PB - American Chemical Society CY - Washington ER - TY - JOUR A1 - Boese, Adrian Daniel A1 - Saalfrank, Peter T1 - CO Molecules on a NaCl(100) Surface: Structures, Energetics, and Vibrational Davydov Splittings at Various Coverages JF - The journal of physical chemistry : C, Nanomaterials and interfaces N2 - In this work, we study the adsorption of CO from low to high coverage at a defect-free NaCl(100) surface by means of duster and periodic models, using highly accurate wave function-based QM:QM embedding as well as density functional theory. At low coverages, the most accurate methods predict a zero-point-corrected adsorption energy of around 13 kJ/mol, and the CO molecules are found to be oriented perpendicular to the surface. At higher coverages, lower-energy phases with nonparallel/upright, tilted orientations emerge. Besides the well-known p(2 x 1)/antiparallel phase (T/A), we find other tilted phases (tilted/irregular, T/I; tilted/spiral, T/S) as local minima. Vibrational frequencies for CO adsorbed on NaCl(100) and Davydov splittings of the C-O stretch vibration are also determined. The IR spectra are characteristic fingerprints for the relative orientation of CO molecules and may therefore be used as sensitive probes to distinguish parallel/upright from various tilted adsorption phases. Y1 - 2016 U6 - https://doi.org/10.1021/acs.jpcc.6b03726 SN - 1932-7447 VL - 120 SP - 12637 EP - 12653 PB - American Chemical Society CY - Washington ER - TY - THES A1 - Boese, Adrian Daniel T1 - Theorie und Berechnung intermolekularer Wechselwirkungen T1 - Theory and computation of intermolecular interactions N2 - Die klassische Physik/Chemie unterscheidet zwischen drei Bindungstypen: Der kovalenten Bindung, der ionischen Bindung und der metallischen Bindung. Moleküle untereinander werden hingegen durch schwache Wechselwirkungen zusammen gehalten, sie sind trotz ihrer schwachen Kräfte weniger verstanden, aber dabei nicht weniger wichtig. In zukunftsweisenden Gebieten wie der Nanotechnologie, der Supramolekularen Chemie und Biochemie sind sie von elementarer Bedeutung. Um schwache, intermolekulare Wechselwirkungen zu beschreiben, vorauszusagen und zu verstehen, sind sie zunächst theoretisch zu erfassen. Hierzu gehören verschiedene quantenchemische Methoden, die in dieser Arbeit vorgestellt, verglichen, weiterentwickelt und schließlich auch exemplarisch auf Problemstellungen in der Chemie angewendet werden. Aufbauend auf einer Hierarchie von Methoden unterschiedlicher Genauigkeit werden sie für diese Ziele eingesetzt, ausgearbeitet und kombiniert. Berechnet wird die Elektronenstruktur, also die Verteilung und Energie von Elektronen, die im Wesentlichen die Atome zusammen halten. Da Ungenauigkeiten von der Beschreibung der Elektronenstruktur von den verwendeten Methoden abhängen, kann man die Effekte detailliert untersuchen, sie beschreiben und darauf aufbauend weiter entwickeln, um sie anschließend an verschiedenen Modellen zu testen. Die Geschwindigkeit der Berechnungen mit modernen Computern ist eine wesentliche, zu berücksichtigende Komponente, da im Allgemeinen die Genauigkeit mit der Rechenzeit exponentiell steigt, und die damit an die Grenzen der Möglichkeiten stoßen muss. Die genaueste der verwendeten Methoden basiert auf der Coupled-Cluster-Theorie, die sehr gute Voraussagen ermöglicht. Für diese wird eine sogenannte spektroskopische Genauigkeit mit Abweichungen von wenigen Wellenzahlen erzielt, was Vergleiche mit experimentellen Daten zeigen. Eine Möglichkeit zur Näherung von hochgenauen Methoden basiert auf der Dichtefunktionaltheorie: Hier wurde das „Boese-Martin for Kinetics“ (BMK)-Funktional entwickelt, dessen Funktionalform sich in vielen nach 2010 veröffentlichten Dichtefunktionalen wiederfindet. Mit Hilfe der genaueren Methoden lassen sich schließlich semiempirische Kraftfelder zur Beschreibung intermolekularer Wechselwirkungen für individuelle Systeme parametrisieren, diese benötigen weit weniger Rechenzeit als die Methoden, die auf der genauen Berechnung der Elektronenstruktur von Molekülen beruhen. Für größere Systeme lassen sich auch verschiedene Methoden kombinieren. Dabei wurden Einbettungsverfahren verfeinert und mit neuen methodischen Ansätzen vorgeschlagen. Sie verwenden sowohl die symmetrieadaptierte Störungstheorie als auch die quantenchemische Einbettung von Fragmenten in größere, quantenchemisch berechnete Systeme. Die Entwicklungen neuer Methoden beziehen ihren Wert im Wesentlichen durch deren Anwendung: In dieser Arbeit standen zunächst die Wasserstoffbrücken im Vordergrund. Sie zählen zu den stärkeren intermolekularen Wechselwirkungen und sind nach wie vor eine Herausforderung. Im Gegensatz dazu sind van-der-Waals Wechselwirkungen relativ einfach durch Kraftfelder zu beschreiben. Deshalb sind viele der heute verwendeten Methoden für Systeme, in denen Wasserstoffbrücken dominieren, vergleichsweise schlecht. Eine Untersuchung molekularer Aggregate mit Auswirkungen intermolekularer Wechselwirkungen auf die Schwingungsfrequenzen von Molekülen schließt sich an. Dabei wird auch über die sogenannte starrer-Rotor-harmonischer-Oszillator-Näherung hinausgegangen. Eine weitreichende Anwendung behandelt Adsorbate, hier die von Molekülen auf ionischen/metallischen Oberflächen. Sie können mit ähnlichen Methoden behandelt werden wie die intermolekularen Wechselwirkungen, und sind mit speziellen Einbettungsverfahren sehr genau zu beschreiben. Die Resultate dieser theoretischen Berechnungen stimulierten eine Neubewertung der bislang bekannten experimentellen Ergebnisse. Molekulare Kristalle sind ein äußerst wichtiges Forschungsgebiet. Sie werden durch schwache Wechselwirkungen zusammengehalten, die von van-der-Waals Kräften bis zu Wasserstoffbrücken reichen. Auch hier wurden neuentwickelte Methoden eingesetzt, die eine interessante, mindestens ebenso genaue Alternative zu den derzeit gängigen Methoden darstellen. Von daher sind die entwickelten Methoden, als auch deren Anwendung äußerst vielfältig. Die behandelten Berechnungen der Elektronenstruktur erstrecken sich von den sogenannten post-Hartree-Fock-Methoden über den Einsatz der Dichtefunktionaltheorie bis zu semiempirischen Kraftfeldern und deren Kombinationen. Die Anwendung reicht von einzelnen Molekülen in der Gasphase über die Adsorption auf Oberflächen bis zum molekularen Festkörper. N2 - In classical chemistry and physics, three basic categories of chemical bonds between atoms are usually distinguished, the covalent, ionic and metallic bond. However, molecules are held together by so-called weak intermolecular interactions. Despite the rather small forces working together in supramolecular chemistry, which is the chemistry beyond the molecules, they are not of minor importance. The solid state of organic material is based on such forces and in future subjects like nanotechnology or biochemistry, they are of utmost interest. Unfortunately, they are little understand For a better knowledge it is desired to theoretically characterise, to compute the interactions and finally to compare them with experimental data. Quantum mechanical methods are the choice of methods, which are introduced, developed, and compared among each other in this work. Based on a hierarchy of methods of different accuracy, they are applied, elaborated and combined. Finally, they are employed to chemical problems. Such calculations are performed on the electronic structure, i.e. the distribution and energy of the electrons, which basically keep the atoms together. Since the accuracy of the methods depends on the accurate description of the electronic structure, we can scrutinize the effects, describe and further develop it in order to test it on chemical model systems. The speed of such calculations with modern computers is very important, as the accuracy of electronic structure methods increases exponentially with computing time that is of course limited. One of the most accurate electronic structure methods is based on the coupled-cluster theory, this is considered as the “gold standard” in quantum chemistry. It yields ‘spectroscopic accuracy’ for intermolecular interactions with deviations of just a couple of wavenumbers, as comparisons with experiments demonstrate. Another possibility to approximate such high-accuracy methods is density functional theory: the functional form of the “Boese-Martin for Kinetics” functional, introduced in this work, is nowadays utilized in several newly developed functionals. With such methods, semiempirical force fields can be parametrized for the description of intermolecular interactions for certain systems. As semiempirical force fields use a fraction of the computer time needed compared to density functional theory or coupled-cluster methods, this provides a big advantage. For even larger systems, we can combine such electronic structure methods. Here, embedding methods have been refined and developed. These include symmetry adapted perturbation theory as well as the quantum mechanical embedding of fragments into larger systems which are also computed with more approximate quantum mechanical methods. The development of such novel methods is more valuable when they are applied to chemical systems and problems: In this work, we first focused on hydrogen bonds. These are relatively strong intermolecular interactions and are likewise theoretically rather challenging, whereas van-der-Waals interactions are described rather well by simple force fields. Thus, the former hydrogen bonds are still inadequately described by many modern methods. Furthermore, molecular aggregates have been investigated with special attention to their vibrational frequencies, for which we went beyond the rigid-rotor-harmonic-oscillator approximation. Another application is on adsorbates of molecules at ionic and metallic surfaces. Here, similar methods can be applied as for the intermolecular interactions, and they can be calculated to high accuracy by using special embedding techniques. The outcome of these computations stimulated a re-assessment of experimental benchmark values. Finally, the characterization of molecular crystals is an increasingly growing field of research. These crystals are held together by weak interactions, ranging from van-der-Waals interactions to hydrogen bonds. The newly developed methods were used to compute such crystals and are an alternative to the currently employed methods. Thus, the developed methods as well as their application span a wide range: from post-Hartree-Fock methods to density functional theory and semiempirical force fields and their combinations. Their application reaches from the intermolecular interactions from single gas phase molecules to surfaces up to the solid state in molecular crystals. KW - intermolekulare Wechselwirkungen KW - Dichtefunktionaltheorie KW - Post-Hartree-Fock Methoden KW - Molekülkristalle KW - Oberflächen KW - anharmonische Kraftfelder KW - intermolecular interactions KW - molecular crystals KW - surfaces KW - anharmonic force fields KW - Density Functional Theory KW - Post-Hartree Fock methods Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-412867 ER -