@prefix vivo: .
@prefix edm: .
@prefix ns0: .
@prefix dcterms: .
@prefix skos: .
vivo:departmentOrSchool "Science, Faculty of"@en, "Chemistry, Department of"@en ;
edm:dataProvider "DSpace"@en ;
ns0:degreeCampus "UBCV"@en ;
dcterms:creator "Carlson, Colton Dean"@en ;
dcterms:issued "2021-04-27T17:18:52Z"@en, "2021"@en ;
vivo:relatedDegree "Master of Science - MSc"@en ;
ns0:degreeGrantor "University of British Columbia"@en ;
dcterms:description "It is well known that amino acids are the building blocks for proteins, and as such, likely played a critical role in the origin of biological life. When in the solid phase, or in solution, amino acids take on a zwitterionic form, enabling a plethora of stabilizing intermolecular interactions. However, prior to understanding amino acids in the bulk phase, one must first characterize them in their neutral form. Due to the conformational flexibly and thermal instability of amino acids, this is best accomplished through the use of high-resolution matrix isolation Fourier transform spectroscopy in conjunction with high level ab initio calculations. The vibrational spectrum of alpha-serine is reported for the first time within a solid parahydrogen matrix, in which ten conformers were identified, thereby supporting the claim that solid parahydrogen matrices are superior to other noble gas matrices for the analysis of highly flexible molecules. Assignment of the vibrational transitions was accomplished through the use of quantum chemical calculations at both the DFT and MP2 levels of theory. Additionally, preliminary TD-DFT calculations indicate that the first singlet excited state of alpha-alanine is a dissociative state, in which, upon excitation alpha-alanine forms the hydrocarboxyl (HOCO) and ethylamine radicals through a Norrish type I pathway. Interestingly, it is predicted that both the first and second singlet excited states are dissociative states for alpha-alanine conformers displaying strong hydrogen bonds between the carboxyl hydrogen and amine lone pair, indicating that noncovalent interactions play a substantial role in both the ground and excited state potential energy surfaces of amino acids."@en ;
edm:aggregatedCHO "https://circle.library.ubc.ca/rest/handle/2429/77997?expand=metadata"@en ;
skos:note "Exploring the Conformational Landscape of Amino Acids in Solid Parahydrogen Matrices through Ab Initio Calculations by Colton Dean Carlson B.Sc., University of Alberta, 2018 A THESIS SUBMITTED IN PARTIAL FULFULLMENT OF THE REQUIREMENTS FOR THE DEGREE OF Master of Science in THE FACULTY OF GRADUATE AND POSTDOCTORAL STUDIES (Chemistry) The University of British Columbia (Vancouver) April 2021 ยฉ Colton Dean Carlson, 2021 ii The following individuals certify that they have read, and recommend to the Faculty of Graduate and Postdoctoral Studies for acceptance, the thesis entitled: Exploring the Conformational Landscape of Amino Acids in Solid Parahydrogen Matrices through Ab Initio Calculations submitted by Colton Dean Carlson in partial fulfillment of the requirements for the degree of Master of Science in Chemistry Examining Committee: Takamasa Momose, Department of Chemistry, UBC Supervisor Mark Thachuk, Department of Chemistry, UBC Supervisory Committee Member Andrew MacFarlane, Department of Chemistry, UBC Supervisory Committee Member iii Abstract It is well known that amino acids are the building blocks for proteins, and as such, likely played a critical role in the origin of biological life. When in the solid phase, or in solution, amino acids take on a zwitterionic form, enabling a plethora of stabilizing intermolecular interactions. However, prior to understanding amino acids in the bulk phase, one must first characterize them in their neutral form. Due to the conformational flexibly and thermal instability of amino acids, this is best accomplished through the use of high-resolution matrix isolation Fourier transform spectroscopy in conjunction with high level ab initio calculations. The vibrational spectrum of ๐ผ-serine is reported for the first time within a solid parahydrogen matrix, in which ten conformers were identified, thereby supporting the claim that solid parahydrogen matrices are superior to other noble gas matrices for the analysis of highly flexible molecules. Assignment of the vibrational transitions was accomplished through the use of quantum chemical calculations at both the DFT and MP2 levels of theory. Additionally, preliminary TD-DFT calculations indicate that the first singlet excited state of ๐ผ-alanine is a dissociative state, in which, upon excitation ๐ผ-alanine forms the hydrocarboxyl (HOCO) and ethylamine radicals through a Norrish type I pathway. Interestingly, it is predicted that both the first and second singlet excited states are dissociative states for ๐ผ-alanine conformers displaying strong hydrogen bonds between the carboxyl hydrogen and amine lone pair, indicating that noncovalent interactions play a substantial role in both the ground and excited state potential energy surfaces of amino acids. iv Lay Summary Prior to understanding the role of amino acids in complex biological systems, one must first understand how they behave in isolated environments. The vibrational frequencies of these molecules enable one to assess the strength of bonds within the molecule. Thereby allowing for one to experimentally identify the physicochemical properties of the molecule. In this thesis the conformational landscape of gaseous ๐ผ-serine trapped in solid parahydrogen matrices is investigated through the use of Fourier transform infrared spectroscopy (FTIR) and quantum chemical calculations. In addition, the excited states of neutral ๐ผ-alanine are briefly discussed through the context of quantum chemical calculations. v Preface All of the work presented in this thesis was completed under the supervision of Dr. Takamasa Momose in the Chemistry Department at the University of British Columbia, Vancouver campus. Experimental spectra of ๐ผ-serine presented within this thesis were collected by Shin Toh and Alethea Lee during the summer of 2016. I was responsible for the analysis of said experimental data, including all relevant quantum chemical calculations. Additionally, all of the calculations presented for ๐ผ-alanine were performed by me. All ab initio calculations presented herein were completed on Computed Canada clusters. vi Table of Contents ABSTRACT ................................................................................................................. iii LAY SUMMARY .......................................................................................................... iv PREFACE ..................................................................................................................... v TABLE OF CONTENTS ................................................................................................. vi LIST OF TABLES ........................................................................................................ viii LIST OF FIGURES ........................................................................................................ ix GLOSSARY ................................................................................................................. xi ACKNOWLEDGMENTS .............................................................................................. xiii 1. INTRODUCTION ...................................................................................................... 1 1.1. AMINO ACIDS AND MOLECULAR SPECTROSCOPY .................................................................. 1 1.1.1. Introduction to ๐ผ-serine.................................................................................................. 2 1.1.2. Introduction to ๐ผ-alanine ............................................................................................... 6 1.2. MATRIX ISOLATION SPECTROSCOPY ................................................................................... 7 1.2.1. Fourier Transform Infrared Spectroscopy ....................................................................... 7 2. EXPERIMENTAL AND COMPUTATIONAL METHODS .................................................. 9 2.1. MATRIX ISOLATION FTIR SPECTROSCOPY ........................................................................... 9 2.1.1. Matrix Isolation Chamber ............................................................................................... 9 2.1.2. Production of Enriched Parahydrogen Gas ..................................................................... 9 2.1.3. Deposition and Sublimation .......................................................................................... 10 2.1.4. FTIR Spectrometer and Measurements ........................................................................ 11 2.2. ELECTRONIC STRUCTURE THEORY .................................................................................... 11 2.2.1. Hartree-Fock Theory ..................................................................................................... 12 2.2.2. Mรธller-Plesset Perturbation Theory.............................................................................. 18 2.2.3. Coupled Cluster Theory ................................................................................................. 21 2.2.4. Truncating the Cluster Operator ................................................................................... 23 2.3. QUANTUM CHEMISTRY: EXCITED STATE CALCULATIONS ...................................................... 27 2.3.1. Configuration Interaction-Singles ................................................................................. 28 vii 2.3.2. Equation of Motion Coupled Cluster Theory ................................................................. 29 2.4. THERMODYNAMICS IN MOLECULAR SPECTROSCOPY ............................................................ 32 2.4.1. Relative Electronic Energy ............................................................................................ 33 2.4.2. Molecular Partition Functions ...................................................................................... 33 2.4.3. Boltzmann Distribution Law ......................................................................................... 38 3. RESULTS AND DISCUSSION .................................................................................... 39 3.1. CONFORMATIONAL DYNAMICS OF ๐ผ-SERINE IN SOLID PARAHYDROGEN MATRICES ................... 39 3.2. EXPERIMENTAL AND COMPUTATIONAL PARAMETERS FOR ๐ผ-SERINE ....................................... 40 3.2.1. Noncovalent Interactions in ๐ผ-serine ........................................................................... 44 3.2.2. ๐ผ-serine in Solid Parahydrogen Matrices ..................................................................... 51 3.2.3. ๐ผ-serine Conformational Dependence on Sublimation Temperature .......................... 62 3.3. EXCITED STATE DYNAMICS OF ๐ผ-ALANINE......................................................................... 63 3.3.1. Computational Procedures for ๐ผ-alanine ..................................................................... 63 3.3.2. Lowest Lying Electronic States for ๐ผ-alanine................................................................ 67 4. CONCLUSION ........................................................................................................ 74 BIBLIOGRAPHY ......................................................................................................... 75 APPENDICES ............................................................................................................. 86 A. SUPPLEMENTAL MATERIAL FOR THE EXPERIMENTAL AND COMPUTATIONAL METHODS ................. 86 B. SUPPLEMENTAL MATERIAL FOR THE CONFORMATIONAL ANALYSIS OF SERINE ............................. 87 viii List of Tables Table 1.1 Summary of the torsion considered by Gronert and OโHair for ๐ผ-serine. .................................... 3 Table 1.2 Summary of the torsions considered by He and Allen for ๐ผ-serine. ............................................. 4 Table 3.1. Tabulation of the relative ZPE-corrected electronic energies, and relative Gibbs free energies (420 K) for the 15 lowest energy ๐ผ-serine conformers. .............................................................................. 45 Table 3.2 Summary of noncovalent interactions present in each of the 15 lowest energy conformers. ..... 51 Table 3.3 Experimental assignment of the mid-IR vibrational bands for ๐ผ-serine in solid parahydrogen. 58 Table 3.4 Tabulation of the relative Gibbs free energies at 400, 410, 420, and 430 K for the 15 lowest energy conformers of ๐ผ-serine..................................................................................................................... 63 Table 3.5 Vertical excitation energies for the first three singlet excited states for all 13 conformers of ๐ผ-alanine at the EOM-CCSD/aug-cc-pVTZ and CIS(D)/aug-cc-pVTZ levels of theory. .............................. 70 Table B.1 Nomenclature conversion table for the conformers of serine in select recent investigations. .... 90 Table B.2 Definition of internal coordinates, Si, for use in the normal mode analysis and total energy distribution reported for each normal mode. Atom numbering corresponds to that shown in Figure 3.1. . 91 Table B.3 Conversion of internal coordinates, Si, to those presented by He and Allen (Ref. 18). .............. 91 Table B.4 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 1 and 2, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. ........................ 92 Table B.5 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 5 and 4, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. ........................ 93 Table B.6 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 5 and 4, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. ........................ 94 Table B.7 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 7 and 8, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. ........................ 95 Table B.8 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 9 and 10, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. ...................... 96 Table B.9 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 11 and 12, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. ...................... 97 Table B.10 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 13 and 14, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. ................. 98 Table B.11 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformer 15, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. ............................. 99 ix List of Figures Figure 1.1 Skeletal formulas for the L- and D-enantiomers of ๐ผ-serine (top) and ๐ผ-alanine (bottom). L-enantiomer is on the left for both molecules, and ๐ผ denotes the chiral carbon. ............................................ 2 Figure 3.1 Three-dimensional molecular structures of the 15 lowest energy conformers of ๐ผ-serine. ....... 44 Figure 3.2 NCI analysis for the conformers of ๐ผ-serine. The colour bar on the right denotes whether the interaction is attractive, repulsive, or weak, denoted by red, blue, green respectively. .............................. 48 Figure 3.3 Projection of QTAIM critical points onto ๐ผ-serine. Nuclear critical points are purple, ring critical points are blue, and bond critical points are yellow. ................................................................................... 50 Figure 3.4 Zoomed in sections of the 700 โ 1000 cm-1 regions of the experimental ๐ผ-serine spectrum in solid parahydrogen. The TSub = 430 K spectrum has been artificially situated above the other spectra in C, D, and E such that the larger amounts of noise do not interfere with observing other spectral features. Additionally, the TSub = 430 K spectrum was not included in A and B due to an overwhelming amount of noise. ............................................................................................................................................................ 53 Figure 3.5 Zoomed in sections of the amine (N-H), A, and carbonyl (C=O), B, stretching regions for ๐ผ-serine in solid parahydrogen. The TSub = 430 K spectrum has been artificially situated above the other spectra in A such that the larger amounts of noise do not interfere with observing other spectral features. ..................................................................................................................................................................... 55 Figure 3.6 Zoomed in sections of the methyl stretching regions for ๐ผ-serine in solid parahydrogen. The TSub = 430 K spectrum has been artificially situated above the other spectra in A such that the larger amounts of noise do not interfere with observing other spectral features. ................................................................. 56 Figure 3.7 Zoomed in sections of the hydroxyl stretching bands for ๐ผ-serine in solid parahydrogen. The legend shown in Figure A (top right) is the same for all of the figures shown within. The TSub = 430 K spectrum has been artificially situated above the other spectra in B through E such that the larger amounts of noise do not interfere with observing other spectral features. ................................................................. 58 Figure 3.8 Three-dimensional molecular structures of the 13 lowest energy conformers of ๐ผ-alanine. ..... 66 Figure 3.9 Electron-hole analysis for the four lowest energy conformers of ๐ผ-alanine. Blue mesh corresponds to the โelectronโ, and red mesh corresponds to the โholeโ. .................................................... 69 Figure 3.10 NCI analysis for the conformers of ๐ผ-alanine. The colour bar on the right denotes whether the interaction is attractive, repulsive, or weak, denoted by red, blue, green respectively. .............................. 72 Figure 3.11 Projection of QTAIM critical points onto ๐ผ-serine. Nuclear critical points are purple, ring critical points are blue, and bond critical points are yellow. ....................................................................... 73 Figure B.1 FTIR spectrum of L-serine immediately after deposition in the solid parahydrogen matrix (flow rate = 5 ccm, TSub = 400 K, TDep. = 4 K). ................................................................................................... 87 x Figure B.2 FTIR spectrum of L-serine immediately after deposition in the solid parahydrogen matrix (flow rate = 5 ccm, TSub = 410 K, TDep. = 4 K). ................................................................................................... 88 Figure B.3 FTIR spectrum of L-serine immediately after deposition in the solid parahydrogen matrix (flow rate = 5 ccm, TSub = 420 K, TDep. = 4 K). ................................................................................................... 89 Figure B.4 FTIR spectrum of L-serine immediately after deposition in the solid parahydrogen matrix (flow rate = 5 ccm, TSub = 430 K, TDep. = 4 K). ................................................................................................... 90 xi Glossary AM1 Semiempirical Austin Model 1 Method aug-cc-pVTZ Dunningโs Correlation Consistent Augmented Triple Zeta Basis Set with Additional Polarization and Diffuse Functions. BJ Becke-Johnson Damping B3LYP Beckeโs 3-Parameter Lee-Yang-Parr Hybrid Functional B3LYP-D3BJ B3LYP Functional with D3 and BJ Empirical Corrections CBS Complete Basis Set CI Configuration Interaction CIS Configuration Interaction Singles CIS(D) Configuration Interaction Singles and Perturbative Doubles CC Coupled Cluster Theory ccm Cubic Centimetre per Minute cc-pVTZ Dunningโs Correlation Consistent Augmented Triple Zeta Basis Set with Additional Polarization Functions. cc-pVTZ-F12 Correlation Consistent Triple Zeta Auxiliary Basis Set with Additional Polarization Functions for Explicitly Correlated Calculations cc-pVQZ/C Correlation Consistent Quadrupole Zeta Auxiliary Basis Set for the Resolution of Identity Approximation CCSD Coupled Cluster Singles and Doubles CCSD(T) Coupled Cluster Singles, Doubles, and Perturbative Triples CCSD(T)-F12 Explicitly Correlated Coupled Cluster Singles, Doubles, and Perturbative Triples DC Direct Current def2-TZVPD Alhrichs Split-Valence Triple Zeta Property Optimized Basis Set with Additional Polarization and Diffuse Functions DFT Density Functional Theory DVPT2 Depurturbed Second Order Vibrational Perturbation Theory D3 Grimmeโs Empirical Atom-Pairwise Dispersion Correction EOM Equation of Motion FCI Full Configuration Interaction Theory FPA Focal Point Analysis xii FTIR Fourier Transform Infrared FWHM Full Width Half Max HB Hydrogen Bond HF Hartree-Fock IR Infrared MCT Mercury Cadmium Telluride MI Matrix Isolation MI-FTIR Matrix Isolation Fourier Transform Infrared MIR Mid-Infrared MPn nth Order Mรธller-Plesset Perturbation Theory MP2 Second Order Mรธller-Plesset Perturbation Theory MW Microwave NCI Noncovalent Interactions NIR Near-Infrared NTC Negative Temperature Coefficient PED Potential Energy Distribution PES Potential Energy Surface QTAIM Quantum Theory of Atoms in Molecules RI Resolution of Identity SPT Simple Perturbation Theory UV Ultra-Violet VPT2 Second Order Vibrational Perturbation Theory VEDA Vibrational Energy Distribution Analysis ZPE Zero-Point Vibrational Energy xiii Acknowledgments I would like to express my appreciation to my research supervisor, Dr. Takamasa Momose, for his guidance and support throughout my studies. In addition, I would like to express my appreciation to Dr. Manish Vashishta, who taught me a number of different laboratory skills, and provided a helpful hand whenever needed. I would also like to thank current and previous lab mates for their fruitful discussions and banter, including, in no particular order, Dr. Eric Miller, Jiahong Hu, Hatsuki Otani, and Brendan Moore. Finally, I offer my sincere gratitude to my family and girlfriend, Julia Lu, for their extended and unconditional support throughout my studies. 1 1. Introduction 1.1. Amino Acids and Molecular Spectroscopy Organic compounds containing an amino and carboxyl functional group are typically referred to as amino acids, and they are an essential component to biological life. In fact, it has been hypothesized that amino acids were deeply involved in the prebiotic processes responsible for triggering the origin of life.1 As a result of their biological importance, amino acids have been the subject of countless investigations spanning every scientific discipline. In total there are over 500 naturally occurring amino acids, however, only 20 of these appear in the genetic code of biological species.2 When amino acids are dissolved in solution or in the solid phase, they take on a zwitterionic form enabling stabilizing intermolecular interactions with the surrounding molecules. Unlike in the solid phase, amino acids in the gas phase are chemically neutral and stabilize themselves via the intramolecular interaction of partial charges. These stabilizing intramolecular interactions include hyperconjugation and hydrogen bonding, however, the presence of destabilizing steric crowding is also possible. Prior to understanding the complex behaviors of molecules in biological settings, one must first completely understand said molecules in an environment free of external perturbation. However, as a result of the thermal instability and low vapour pressure of amino acids,3 investigations into their neutral form are rather scarce, and as such, the conformational landscape of neutral amino acids remains largely unknown. In addition to providing crucial information regarding the physicochemical properties of amino acids, investigations into their neutral form enables the identification of amino acids in interstellar space.4 Identification of amino acids in interstellar space became a hot topic upon the identification of amino acids, diamino acids, nucleobases, and hydrocarbons on the Murchison meteorite, which is a carbonaceous chondrite meteorite, in 1969.1 The search for amino acids in interstellar space led to questions regarding the stability of amino acids in the presence of ultra-violet (UV) irradiation. Studies into the photostability of amino acids indicate that amino acids may undergo racemization, 5 conformational isomerism,6 or degradation upon UV irradiation.7 The latter is of particular interest as it indicates that amino acids may be indirectly detected in interstellar space through the detection of the photoproducts. However, just as was the case with the conformational landscape, the photoproducts from neutral amino acids are largely unknown as a result of experimental and computational complexities. 2 One popular technique for probing the conformational landscape of molecules is matrix isolation (MI) spectroscopy.8 Just as the name implies, MI spectroscopy involves the trapping of a target species in an inert matrix kept at cryogenic temperatures, in which the matrix is typically composed of an inert gas such as argon, neon, nitrogen, xenon, or krypton, but in the case of this thesis, solid parahydrogen is utilized. Unlike other gas phase experiments, the MI technique enables one to probe not only the ground state, but also the excited state dynamics of molecules via in-situ irradiation of said trapped molecules. However, coupling the use of MI spectroscopy with ab initio calculations enables one to accurately investigate both, the ground and excited states. The work presented in this thesis aims to extend the previous investigations into the physicochemical properties of ๐ผ-serine through the use of solid parahydrogen MI-FTIR spectroscopy and high-level ab initio calculations. In which the noncovalent intramolecular interactions present in the conformers of ๐ผ-serine will be used to aid in their experimental identification. Variation of the sublimation temperature was utilized to further aid in the experimental assignment. This work represents the first study to consider all known low energy ๐ผ-serine conformers during the assignment process. Additionally, preliminary results for the investigation into the excited states of ๐ผ-alanine are briefly introduced. Figure 1.1 Skeletal formulas for the L- and D-enantiomers of ๐ผ-serine (top) and ๐ผ-alanine (bottom). L-enantiomer is on the left for both molecules, and ๐ผ denotes the chiral carbon. 1.1.1. Introduction to ๐ถ-serine ๐ผ-serine, with a molecular weight of 105.09 g molโ1 is one of the smallest biologically relevant amino acids being only larger than glycine and alanine. As seen in Figure 1.1, five torsions in ๐ผ-serine exist where, upon the internal rotational about said torsion, a new conformer is formed. In 3 order for a given molecular configuration to be a conformer, the configuration must be located at a local minima on the ground state potential energy surface. The internal rotations of interest are summarized below in Table 1.1, along with an additional skeletal formula of ๐ผ-serine which serves as a pictorial representation of said torsions. Of the five torsions of interest, there are two 3-fold internal rotations, two 2-fold, and one 6-fold internal rotation, yielding 324 potential configurations of ๐ผ-serine. Unfortunately, as a result of the substantial number of intramolecular interactions present in ๐ผ-serine, those 324 configurations do not accurately represent the full ensemble of possible conformers. Due to the large number of possible conformers, many theoretical studies did not consider the full set of trial structures, and as result of this, often missed key low energy conformers throughout their analyses.9,10,11,12,13,14,15 Table 1.1 Summary of the torsion considered by Gronert and OโHair for ๐ผ-serine. ๐ Torsion Considered Angles Figure ๐1 C โ C(O) โ O โ H 0, 180 ๐2 N โ C โ C(O) โ O(H) -150, -90, -30, 30, 90, 150 ๐3 C โ C โ N โ H -60, 60, 180 ๐4 N โ C โ C โ O(H) -120, 0, 120 ๐5 C โ C(H2) โ O โ H -120, 0, 120 In 1995, Gronert and OโHair published the first large-scale study on ๐ผ-serine.16 They sampled the 324 trial structures mentioned above and shown in Table 1.1, with the semiempirical AM1 method identifying 73 unique configurations, which, upon refinement at the HF/6-31G* level of theory, was reduced to 51 unique conformations. All of the reported conformers were within 50 kJ mol-1. It is worth mentioning that HF/3-21G* was first used to refine the AM1 structures, reducing the number of stable configurations from 73 to 54. Additional refinement with HF/6-31G* further reduced the number of configurations from 54 to the final value of 51. As a result of the incomplete ensemble of conformers initially considered, a number of low energy conformers were missed. Unfortunately, a number of experimental studies used this publication as a theoretical basis, and as such, they did not consider all possible low energy conformations, which will be discussed below. In 2005, Miao et al. employed B3LYP/6-31G* to screen Gronert and OโHairโs 324 structures, yielding 83 unique configurations, which upon further refinement with B3LYP/6-311+G** reduced that number to 74.17 Interestingly, increasing the level of theory to B3LYP 4 resulted in the discovery of an additional 23 stable configurations, which is the first indication of the importance of electron correlation in ๐ผ-serine. Despite having explicitly mentioned the importance of including the full ensemble of possible conformers, this study suffered the same fate as that of Gronert and OโHair. That is, they missed a number of key low energy conformers as a result of using an incomplete set of trial structures. In 2016, He and Allen sampled 15,552 configurations of ๐ผ-serine at the HF/6-31G* level of theory, identifying 89 unique configurations.18 Hence, consideration of the full ensemble obtained an additional 38 unique configurations. Further refinement with MP2/cc-pVTZ reduced the number of stable configurations to 85. In short, the 15,552 starting configurations were constructed from the simultaneous variation of six large-amplitude modes, which are listed below in Table 1.2. Upon optimization with MP2, 12 configurations were identified to be within 12 kJ mol-1, and their relative energies were calculated with the highly accurate Focal Point Analysis (FPA). At the time this thesis is being written, He and Allen have presented the most comprehensive computational analysis on the ground state configurations of ๐ผ-serine. It is worth mentioning that a quick comparison of Table 1.1 and Table 1.2 show that Gronert and OโHair failed to account for the inversion of NH2, in addition to not considering many possible configurations resulting from rotation about the C โ N and C4 โ C6 torsions. Table 1.2 Summary of the torsions considered by He and Allen for ๐ผ-serine. ๐ Mode Rotation Figurea ๐1 NH2 inversion 2-fold ๐2 C4 โ C6 โ O10 โ H12 2-fold ๐3 C4 โ C7 โ O13 โ H14 2-fold ๐4 N1 โ C4 โ C7 โ O11 6-fold ๐5 C7 โ C4 โ C6 โ O10 12-fold ๐6 C6 โ C4 โN1 โ H2 12-fold a Red corresponds to oxygen, blue to nitrogen, grey to carbon, and white to hydrogen. As mentioned in the introduction, the low vapour pressure and thermal instability of neutral amino acids, the first experimental identification of neutral ๐ผ-serine was not achieved until 2004 when Lambie et al. successfully observed four conformers of ๐ผ-serine through the use of MI-FTIR spectroscopy with an argon host matrix.19 On the basis of the theoretical work performed by Gronert and OโHair, they reoptimized the 12 most stable structures at the B3LYP/6-31++G** level of theory and utilized scaled quantum mechanical force fields to predict the vibrational 5 frequencies. Identification of conformers 1, 2, 3, and 4 was achieved solely through comparison of the experimental spectrum with the predicted spectrum. As a result of this crude analysis, a number of conformers were likely missed during the assignment process. However, this article is significant in that it represents the first time neutral ๐ผ-serine was experimentally observed, which is an accomplishment in and of itself. Shortly afterwards, Jarmelo et al. reported two separate MI studies on ๐ผ-serine.20,21 In these reports they utilized isotopic substitution, and UV photoirradiation (200 nm) techniques to assess the conformational distribution of ๐ผ-serine in argon matrices, and its associated photoproducts. This enabled the assignment of 9 ๐ผ-serine conformers, namely conformers 1 through 9. In addition, they identified ethanolamine, CO2, CO, H2O, and acetamide as the primary photoproducts. Interestingly, two, conformer dependent, photodissociation pathways were identified. Photolysis of conformers exhibiting OHโฏN hydrogen bonding between the alcohol side chain and amine group, yielded ethanolamine and CO2 via decarboxylation. Whereas the photolysis of conformers exhibiting OHโฏN hydrogen bonding between the carboxyl hydrogen and amine group resulted in the formation of acetamide through the decarbonylation. In 2007 Blanco et al. reported the assignment of seven ๐ผ-serine conformers through the use of laser-ablation microwave spectroscopy.22 Namely, they identified conformers Ia, IIb, Iโb, IIc, III๏ขb, III๏ขc, and IIa. For the sake of comparison, a table of conversion for nomenclature used in relevant studies on ๐ผ-serine has been included in the Appendix. The absence of two conformers which were identified by Jarmelo et al. was attributed to collisional cooling during the supersonic expansion. This publication represents the first experimental identification of neutral ๐ผ-serine achieved with laser ablation. Lastly, in 2015 Najbauer et al. coupled MI-FTIR spectroscopy with selective near-infrared (NIR) irradiation to identify eight conformers of ๐ผ-serine.23 Namely, these conformers are 1 through 7 and 10. An argon matrix was utilized for the identification of conformers 1 through 6, with conformers 7 and 10 being confirmed through the use of a nitrogen matrix, which is known to stabilize unstable conformers. That being said, conformers 7 and 10 were assigned to the same experiment bands, hence, only conformer 7 or 10 was assigned as opposed to both. This assignment contradicts those made by Jarmelo et al., and as such, the goal of this thesis is to identify the configurations of ๐ผ-serine present in a solid parahydrogen matrix. 6 1.1.2. Introduction to ๐ถ-alanine ๐ผ-alanine, being only larger than glycine, is the second smallest biologically relevant amino acid and the smallest chiral amino acid. Since this thesis is primarily interested in the excited states of ๐ผ-alanine, the identification of conformers will only be briefly introduced. There have been many investigations into the excited states of ๐ผ-alanine, and as such, only the two most relevant publications will be introduced. Similar to ๐ผ-serine, ๐ผ-alanine has been the subject of many theoretical investigations, in fact, the first comprehensive investigation was published in 1995 by Gronert and OโHair alongside their investigation of ๐ผ-serine.16 In their study they identified ten stable configurations with HF/6-31G*. Shortly afterwards in the same year, Cao et al. identified thirteen conformers of ๐ผ-alanine with MP2/6-31G**.24 In addition to that, Schรคfer et al. utilized HF/6-311G** to optimize all of the trial structures presented by Gronert and OโHair, discovering an additional three stable configurations, for a total of thirteen. Hence, thirteen identical conformers were independently identified by Caro and Schรคfer via MP2 and HF respectively. Shortly thereafter, Csรกzรกr et al. identified 13 stable configurations of ๐ผ-alanine at the MP2(full)/6-311++G** level of theory and reported the CCSD(T) energies.25 It appears that ๐ผ-alanine was a hot topic in the theoretical chemistry community during the years of 1995 and 1996, as all of these publications came out within a few months of one another. After 1996, not many theoretical investigations into the conformational landscape of ๐ผ-alanine were reported. However, in 2011 Balabin et al. reported a FPA on the conformers of ๐ผ-alanine identifying twelve stable configurations at the B3LYP/aug-cc-pVTZ level of theory.26 Lastly, in 2010 Jaeger et al. optimized the two lowest energy conformers of ๐ผ-alanine at the CCSD(T)/cc-pVTZ level of theory, then subsequently utilized the FPA to determine their relative energies.27 However, unlike Balabin, this publication is primarily concerned with determining the semiexperimental structures through the use of rotational constants obtained from microwave spectroscopic experiments. In 2007, Maul et al. probed the electronic excitations for the 3 lowest energy conformers of ๐ผ-alanine through the use of gradient corrected (spin-) DFT with a projector-augmented wave scheme.28 Interestingly, they predict that, upon irradiation with UV light, ๐ผ-alanine undergoes photodissociation. Although it is not explicitly stated, the resulting product structures provided within their publication are that of the hydrocarboxyl (HOCO) and ethylamine radicals. 7 Additionally, the HOCO radical is trans for conformers 1 and 2, and cis for conformer 3. Unlike the experimental investigation on ๐ผ-serine, the same photoproducts are predicted for the three conformers of ๐ผ-alanine. It is worth mentioning that this publication does not include any of the necessary details regarding the photodissociation process. However, as far as the author is aware, this is the only publication showing the photodissociation of ๐ผ-alanine. In 2014 Tia et al. examined the photoionization of ๐ผ-alanine through the coupling of synchrotron radiation with a double imaging photoelectron coincidence spectrometer in the Vacuum UV (VUV) range.29 Through the use of symmetry adapted cluster configuration interaction (SAC-CI) calculations, they determined the vertical ionization potentials to be 9.581, 9.668, and 9.593 eV for the first 3 conformers of ๐ผ-alanine. Additionally, they reported the presence of the HOCO radical at 9.17 eV, and the ethylamine radical at 9.48 eV. Sample deposition was obtained through the use of resistive heating and aerosol thermodesportion. Interestingly, the ethylamine radical was only detected in experiments utilizing the resistive heating technique, whereas the HOCO radical was detected in both. 1.2. Matrix Isolation Spectroscopy 1.2.1. Fourier Transform Infrared Spectroscopy Through the use of the Fourier Transform (FT) any waveform can be rewritten as a sum of sinusoidal functions. Recall that a waveform with a period of ๐ก ๐โ is periodic at intervals of ๐ก. In short, the FT converts time-based waveforms into individual cycles, reporting the amplitude, offset, and rotation speed for each cycle. This will be touched upon after a brief introduction to the Michelson interferometer. The Michelson interferometer functions to modulate optical radiation. In short, the Michelson interferometer consists of a beam splitter, travelling mirror, and stationary mirror. An incoming beam of radiation is split via the beam splitter into two beams of nearly equal power and intensity. One beam is directed towards the stationary mirror, while the other to the travelling mirror. Upon reflection at the mirrors the twin beams recombine at the beam splitter, in such a way that changes in intensity can be measured as a function of the difference in the distance travelled by each beam. In other words, the travelling mirror serves to fluctuate the power incident on the detector in a reproducible manner. For example, when the travelling and stationary mirrors are 8 equal distance from the beam splitter (same path length), the split beams are said to be in phase and the resultant recombined beam is at the maximum signal power. Hence, differences in signal power at the detector are a function of the difference in path length, termed retardation. Plotting power against retardation yields an interferogram, which is a periodic function, and as such, can be analyzed with the FT. In brief, compared to conventional dispersive spectrometers, FT spectrometers yield higher signal power, signal-to-noise ratios, and higher resolution as a result of the multiplex advantage. Additional improvements can be achieved through the use of multiple interferometers, in which the additional interferometers serve to improve spectral quality through indirect measurements, as opposed to directly interacting with the source light. For example, most high-resolution FTIR spectrometers contain a laser-fringe reference system, in which a helium-neon (HeNe) laser is utilized to regulate the travelling mirror. Briefly, monochromatic sources, such as lasers, produce cosine interferograms as maximum power is always achieved when retardation is zero, which can be utilized to regulate the oscillation speed of the travelling mirror. Additionally, one could include a white-light system for the purpose of triggering data sampling. Unfortunately, the maximum resolution attainable by FTIR spectrometers is limited by the travelling mirror. That is, resolution is primarily determined by the step-size of the travelling mirror, putting a physical limit on the attainable resolution. One method to overcome this is through the use of a dual frequency comb spectrometer, in which, the resolution is determined by difference in wavelengths of the two frequency combs.30 Thereby, enabling one to attain resolutions equal to the line-width of the employed combs. 9 2. Experimental and Computational Methods 2.1. Matrix Isolation FTIR Spectroscopy 2.1.1. Matrix Isolation Chamber An in-house built doubly shielded cryogenic chamber was utilized for all the experiments presented in this thesis. The chamber is composed of a stainless-steel outer shell with inner and outer diameters of 20 and 22 cm respectively. Inside the chamber, a copper cold plate, which is 10 cm in diameter, has the sole responsibility of housing the BaF2 wedged window above the cold head. The cryostat utilized in this apparatus is capable of cooling the cold plate to 3.6 K, thereby cooling the window to ca. 4 K. A closed-cycle Gifford-McMahon refrigerator is utilized to cool the cold head. The BaF2 wedge window is 25 mm in diameter with a thickness of 2 mm (mean) and is attached to the cold plate with a 4 by 4 cm copper window holder. The temperature of the chamber was monitored through the use of two silicon diode temperature sensors, which were attached to the cold plate and wedged window. In all experiments the chamber was kept at roughly 3 ร 10โ8 torr through the combined effects of a turbo molecular pump, which is capable of attaining pressures of 3 ร 10โ6 torr, and the cold head which acts as a cryopump. Like cryopumps, the cold head is capable of โremovingโ water vapour at four times the rate turbomolecular pumps are. 2.1.2. Production of Enriched Parahydrogen Gas Enriched parahydrogen gas was produced from normal hydrogen gas (99.99%, Praxair Canada Inc.) through the use of an inhouse-built converter. The converter was constructed by Dr. Takamasa Momose, alongside associates from the University of Illinois and University of British Columbia.31 Unlike other parahydrogen converters, which reach cryogenic temperatures through the use of liquid N2 or H2 cryogens, this one utilizes a closed cycle He cryostat (liquid He cryogen).32 The cold head has been modified to attach a reactor coil onto it, which is composed of oxygen free high conductivity copper. Briefly, conversion from ortho to parahydrogen is catalyzed through the use of 3.3 grams of hydrous ferric oxide catalyst (30 x 50 mesh), which is positioned inside the copper tubing. 10 2.1.3. Deposition and Sublimation Deposition of parahydrogen onto the window was accomplished through the use of steel inlet tubing with a diameter of 1.59 mm, positioned 45ยฐ and 3.5 cm from the window. Flow rate of parahydrogen gas was controlled with a gas flow controller. Whereas deposition of solid samples is performed directly from the tip of the Knudsen cell, which is positioned 45ยฐ and 2.5 cm from the window. The fact that the matrix gas and sample inlets are both situated 45ยฐ from the window ensures that the two beams intersect at a 90ยฐ angle immediately in front of the window, ensuring ample mixing of the two gases. Prior to deposition in the experiments presented in this thesis, the window was cooled to 4 K. It is worth noting that the temperature of the window increases upon deposition of a gas simply due to the temperature of the gas being greater than that of the window. The Knudsen cell used in these experiments was built in-house and is similarly designed to standard resistance-heated Knudsen cells. The cell itself is composed of a 42 by 8 mm copper cell body, and an 18 by 1.6 mm copper nozzle. It is worth noting that the nozzle diameter is 1 mm. Heating is achieved through the use of a ceramic heater attached to the rear of the cell and was powered with a standard DC power supply. The temperature was accurately monitored through the use of a solid state NTC thermistor sensor capable of reporting temperatures with an accuracy of ยฑ1 K, placed directly behind the sample compartment. Due to the positioning of the sensor directly behind the sample compartment, one can only report the temperature directly beside the sensor. That is, the temperature of the entire sample compartment is unknown. However, due to the small size of the sample compartment and the fact that no extreme temperatures were utilized, it is unlikely that the temperature greatly deviates from one side of the compartment to the other. Additionally, the experiments investigating conformational dependence on sublimation temperature are more concerned with the relative temperature than the absolute temperature. That is, the point of modulating the sublimation temperature is to see a change in the relative intensity of peaks, which corresponds to a change in conformational distribution that can be predicted from theory and used to support the spectral assignment. Hence, this non-uniformity in temperature, and the inability of the experimental apparatus to measure the non-uniformity, will not have any effect on the results presented in this thesis. 11 2.1.4. FTIR Spectrometer and Measurements The FTIR employed in the experiments presented here is a Bruker IFS 125HR high-resolution spectrometer that is equipped with a KBr beam splitter, globar MIR source, and a reference He-Ne laser. The spectrometer itself is kept at an approximate pressure of 3 ร 10โ2 torr to maximize the intensity of IR light irradiated onto the sample and is accomplished through the use of a standard mechanical pump. A parabolic mirror was used to deflect the IR light into the matrix isolation chamber and thereby onto the sample. BaF2 windows with identical measurements to the sample window, were utilized to allow the IR beam to travel into and out of the matrix isolation chamber. Upon exiting the chamber, an additional parabolic mirror was used to deflect the beam into the liquid N2 cooled MCT detector situated in a separate optic chamber, which was kept at a pressure similar to that of the spectrometer. All spectra presented in this thesis were recorded with a resolution of 0.1 cmโ1. 2.2. Electronic Structure Theory Throughout the past century (ca. 1920 โ 2020) the application of quantum chemistry to research has exploded from a relatively small group of theorists to being utilized in most chemistry research laboratories. For instance, synthetic chemists commonly apply quantum chemistry for the rationalization and prediction of reaction pathways, spectroscopists use it to aid in spectral analysis, and pharmaceutical chemists apply it to drug discovery.33 The widespread use and advancement of quantum chemistry can largely be attributed to two factors: (1) rapid development of computational hardware (and the associated reduction in cost), and (2) the introduction of โblack boxโ methods for predicting molecular properties.33,34,35 A black box method is a one that can be easily performed with little to no knowledge of the programsโ inner-workings. However, the utilization of these black box methods does not come without cost. This is made abundantly clear by the frequent use of the B3LYP functional36,37 in the prediction of organic chemical reactions. Since 2007, multiple independent studies have shown that the B3LYP functional, without any empirical dispersion corrections, fails to qualitatively and quantitatively describe organic reaction pathways, and the associated thermochemistry for many organic reactions.38,39,40,41 Regardless of these findings, the B3LYP functional is still widely used as the sole theoretical basis supporting experimental results, often resulting in the publication of erroneous information.34 It is important to note that the B3LYP functional does not fail for every organic molecule. In fact, it describes the 12 ground state potential energy surface of some amino acids quite accurately.25 As such, when using any approximate quantum chemical method, it is of the utmost importance to employ at least two different techniques to confirm the validity of the theoretical results. When quantum chemistry is applied to molecular spectroscopy, this is typically accomplished through the use of both density functional theory, and Mรธller-Plesset perturbation theory. The widespread misuse of the B3LYP functional noted above, is largely a result of users not understanding the inner workings of the technique and as such, they do not realize its shortcomings. Due to the complexity of the analyses discussed in this thesis, it is necessary for one to fully understand the techniques applied, as opposed to treating them as a black box. As such, this section aims to briefly introduce the reader to the computational techniques employed for the ground state analysis of ๐ผ-serine. In particular, the Hartree-Fock (HF)42,43 method, density functional theory (DFT),44 Mรธller-Plesset (MPn) perturbation theory,45 and coupled cluster theory (CC)46,47 will be discussed. Additionally, it is worth mentioning that explicit electron correlation in coupled cluster theory,48 and empirical dispersion corrections in DFT, particularly with the B3LYP functional,49 were utilized within this thesis but will not be discussed here. Note that a description of the computational techniques used in the excited state analysis of ๐ผ-alanine can be found in the following section of this thesis (Section 2.3). The equations shown in the following sections were sourced from Frank Jensenโs Introduction to Computational Chemistry textbook, with additional references being provided when required.55 2.2.1. Hartree-Fock Theory Electronic structure theory aims to describe the motion of electrons within atoms and molecules most commonly via the time-independent Schrรถdinger equation, ?ฬ?|๐น๐โฉ = ๐ธ๐|๐น๐โฉ, (2.1) where ?ฬ?, ฮจ๐, and ๐ธ๐ are the Hamiltonian operator, wavefunction, and energy eigenstate respectively. The Hamiltonian operator for a given system corresponds to the total energy of that system. Hence, the molecular Hamiltonian, also called the Coulomb Hamiltonian, represents the kinetic and potential energy of all electrons, nuclei, and their interactions. In atomic units, the molecular Hamiltonian is given by the following expression, 13 ?ฬ? = โโ12โ๐2๐โ ?ฬ?๐โโ12๐๐ดโ๐ด2๐ดโ ?ฬ?๐+โ1|r๐ โ r๐|๐>๐โ ?ฬ?๐๐โโ๐๐ด|r๐ โ R๐ด|๐,๐ดโ ?ฬ?๐๐+โ๐๐ด๐๐ต|R๐ด โ R๐ต|๐ด>๐ตโ ?ฬ?๐๐. (2.2) Kinetic energy operators for electrons and nuclei are ?ฬ?๐ and ?ฬ?๐ respectively, whereas โ๐2 and โ๐ด2 denote the Laplacian with respect to electronic, ๐ซ๐, and nuclear, ๐๐ด, coordinates. ?ฬ?๐๐, ?ฬ?๐๐, and ?ฬ?๐๐, are the potential energy operators corresponding to the Coulombic interactions between electrons, nuclei, and electrons and nuclei respectively. Lastly, ๐๐ด and ๐๐ด denote the charge and mass of nuclei A in atomic units.50 It is worth mentioning that the above equation is non-relativistic, hence, it does not account for relativistic effects such as the relativistic mass, spin-orbit coupling, and the parity violation in chiral molecules. No relativistic calculations were performed within this thesis, and as such, no comments regarding the inclusion of relativistic effects will be made, with the exception of mentioning the article published by Wenjian Liu in 2020, which is a good reference for anyone interested in the subject.51 Solving the full molecular Hamiltonian, Equation (2.2), is currently infeasible for molecules of any size. In response to this a number of approximations will be made to simplify the Hamiltonian into one that we can solve. It is here that the Born-Oppenheimer (BO) approximation first comes into play.52 The BO approximation is a central assumption to molecular spectroscopy, and its origins can be found in the fact that electrons are substantially lighter than nuclei. On the basis of this mass difference, it is reasonable to assume that electron velocity is much greater than that of nuclei, and as such it can be said that, from the perspective of the nuclei, electrons adjust to any change in nuclear coordinates instantaneously.53 As a result of this it is possible to decouple the electron and nuclear terms in the molecular Hamiltonian, resulting in an electronic and nuclear Hamiltonian. Respectively, the electronic and nuclear Hamiltonians have the form, ?ฬ?๐ = โโ12โ๐2๐โ ?ฬ?๐+โ1|๐ซ๐ โ ๐ซ๐|๐>๐โ ?ฬ?๐๐โโ๐๐ด|๐ซ๐ โ ๐๐จ|๐,๐ดโ ?ฬ?๐๐, (2.3) and, ?ฬ?๐ = โโ12๐๐ดโ๐ด2๐ดโ ?ฬ?๐+โ๐๐ด๐๐ต|๐๐จ โ ๐๐ฉ|๐ด>๐ตโ ?ฬ?๐๐+๐ธ๐(R)โ From 2.6. (2.4) Under the BO approximation the total wavefunction, ฮจ, is given by the product of the nuclear, ๐๐, and electronic, ๐๐, wavefunctions: 14 ฮจ(๐, ๐ซ) = ๐๐(๐)๐๐(๐ซ,๐). (2.5) It is interesting to note that the nuclear wavefunction depends only on nuclear coordinates, but the electronic wavefunction depends on electronic and nuclear coordinates. The electronic wavefunctions dependence on nuclear coordinates is commonly interpreted as being parametrically dependent on electron โ nuclei separation distances arising from the ?ฬ?๐๐ term in Equation (2.3). Substituting the electronic Hamiltonian and wavefunction into the Schrรถdinger equation gives the following expression, ?ฬ?๐๐๐(๐ซ, ๐) = ๐ธ๐(๐)๐๐(๐ซ,๐). (2.6) Where the electronic energy is given by ๐ธ๐ and its dependence on nuclear position is the underlying basis of the adiabatic potential energy surface (PES), which is a central concept in molecular spectroscopy.54 The nuclear Schrรถdinger equation describes the rotational, vibrational, and translational motion of the nuclei, and can be obtained in a fashion identical to that done for the electronic Schrรถdinger equation, ?ฬ?๐๐๐(๐) = ๐ธ๐๐๐(๐). (2.7) Hence, the electronic energies, as a result of their parametric dependence on nuclear geometry, provide a potential on which the nuclei can move. As a result of the two-electron Coulombic potential term in the electronic Hamiltonian, ?ฬ?๐๐, Equation (2.6), obtaining an exact solution to the electronic Schrรถdinger equation for any molecule with more than one electron completely infeasible. That being the case we can approximate the two-electron potential as a sum of one-electron terms, making the Hamiltonian solvable for molecules with more than a single electron. A few different approaches to this approximation exist, however, the most popular method is termed the Hartree approximation.42 In this approximation the two-electron potentials are approximated through the use of a mean field, thereby allowing the total electronic wavefunction to be represented as the product of one-electron spatial orbitals.42,54,55 That is, it is assumed that the electrons in a given system are independent and interact only through a mean electric field, making each potential term dependent on the coordinates of a single electron as opposed to the coordinates of every electron. This method is also referred to as the mean field approximation since the individual fields are replaced by a mean field. It is worth noting that Hartree approximated the electronic wavefunction as simple one-electron wavefunctions. In order to appropriately represent the wavefunction as a product of one-electron spatial orbitals, one should recall that the total electronic wavefunction must adhere to the Pauli principle. 15 Additionally, since electrons are fermions, the wavefunction must be antisymmetric upon exchange.54 This is accomplished through the product of a spatial orbital ๐๐(๐ซ) and a spin function. As mentioned above, electrons are fermions, so spin can be either +12 or โ12. The resulting wavefunction from this product is termed the spin orbital, ๐๐(๐), which is a function of the electrons spatial and spin coordinates, ๐ซ and ๐, respectively: ๐๐(๐) = |๐โฉ = ๐๐(๐ซ๐ , ๐๐). (2.8) In other words, each one-electron spatial orbital generates two orthonormal spin orbitals that are a function of the electrons spatial and spin coordinates. Note that spin orbitals are sometimes referred to as molecular orbitals. Since each spatial orbital generates two spin orbitals, there will always be more spin orbitals than there are electrons to fill them. It is standard to refer to the unoccupied spin orbitals as virtual orbitals. Recall that the electronic wavefunction must also be antisymmetric upon exchange, which is defined as ๐(๐,โฆ , ๐, โฆ , ๐) = โ๐(๐, โฆ , ๐, โฆ , ๐). This condition is met by constructing the total electronic wavefunction from the antisymmetrized product of spin orbitals. That is, the electronic wavefunction can be constructed as a Slater determinant of ๐ orthonormal spin orbitals ๐๐, ๐๐(1, 2,โฏ , ๐) =1โ๐!det |๐1(1) โฏ ๐1(๐)โฎ โฑ โฎ๐๐(1) โฏ ๐๐(๐)|. (2.9) It is worth mentioning that through the procedure for creating spin orbitals and taking their antisymmetrized product, the Pauli Exclusion principle is automatically satisfied.56 The use of the above wavefunction was discovered by Fock whom improved upon Hartreeโs simple one-electron wavefunction by accounting for electron antisymmetry. The Hartree-Fock energy can be calculated by substituting the electronic Hamiltonian and wavefunction, Equations (2.3) and (2.9), into the Schrรถdinger equation, which yields the following expression for the HF energy, ๐ธHF, ๐ธHF = โจ๐๐|?ฬ?๐|๐๐โฉ =โโจ๐|โฬ|๐โฉ๐+12โ(โจ๐๐โ๐๐โฉ)๐๐, (2.10) where โฬ contains the kinetic energy, ?ฬ?๐, and electron-nuclei, ?ฬ?๐๐, interaction terms via integration over the electronic spatial and spin coordinates, 16 โจ๐|โฬ|๐โฉ = โซ๐๐โ(1) [12โ๐2 โโ๐๐ด|r๐ โ R๐ด|๐ด]โ โฬ(๐)๐๐(1)๐๐๐ = โ๐๐ , (2.11) where ๐๐ represents the spatial and spin coordinates of electron ๐. Whereas the antisymmetrized integrals, โจ๐๐|๐๐โฉ โ โจ๐๐|๐๐โฉ = โจ๐๐โ๐๐โฉ, refer to the Coulomb and exchange interaction, โจ๐๐|๐๐โฉ = โซ๐๐1๐๐โ(1)โซ๐๐2๐๐โ(2)1๐12๐๐(2)๐๐(1)โ ๐ฝ๐(1)๐๐(1), (2.12) and, โจ๐๐|๐๐โฉ = โซ๐๐1๐๐1๐๐โ(1)โซ๐๐2๐๐โ(2)1๐12๐๐(2)๐๐(1)โ ๐พ๐(1)๐๐(1). (2.13) In the above expressions the Coulomb and exchange operators are ๐ฝ๐ and ๐พ๐ respectively. Combining the above two equations: โจ๐๐โ๐๐โฉ = โซ๐๐1๐๐โ(1)๐ฝ๐(1)๐๐(1) โ โซ๐๐1๐๐โ(1)๐พ๐(1)๐๐(1). (2.14) The exchange contribution arises from the indistinguishability of electrons, whereas the Coulomb contribution arises from the mean-field approximation previously discussed in this section. Minimization of the HF energy is achieved through Self Consistent Field theory, see Algorithm 1 of the Appendix A, which is a variational method for solving the HF energy by means of optimizing the molecular orbitals. This requires the molecular orbitals to be orthonormal, โซ๐๐1๐๐โ(1)๐๐(1) = ๐ฟ๐๐ . (2.15) It can be shown that application of the variational procedure leads to the Hartree-Fock equation, ?ฬ?๐๐ = [โฬ +โ(๐ฝ๐ โ๐พ๐)๐] ๐๐ = ๐๐๐๐ , (2.16) where ?ฬ? is the Fock operator and ๐๐ are the orbital energies. It is worth mentioning that the Fock operator is a one-electron operator that depends on the molecular orbital ๐๐. It is useful to define the one-electron HF potential operator as ?ฬ?๐ป๐น, which from Equation (2.16) is: ?ฬ?๐ป๐น =โ(๐ฝ๐ โ ๐พ๐)๐๐, (2.17) where ๐ represents the orbital and ๐ is the total number of occupied orbitals. 17 The HF energy can, theoretically, be calculated directly from Equation (2.10). However, doing so requires numerical integration along a finely spaced three-dimensional grid, making it infeasible for a molecule with more than a few electrons. In the 1950s Roothaan57 and Hall58 discovered that the HF equation can be represented as a set of matrix equations by expanding the spatial orbitals into a basis of atomic orbitals (AO), thereby, reducing the HF equation to a single matrix eigenvalue problem that can routinely solved. Expansion of the spatial component, ๐๐, for each spin orbital in a Slater determinant of a finite basis set, {๐๐}, is given by: ๐๐ =โ๐ถ๐๐๐๐๐ , (2.18) where ๐ถ๐๐ is the expansion coefficient for spin orbital ๐. The expansion coefficients are determined through the HF-SCF procedure. It is now convenient to define the density matrix as: ๐๐๐ =โ2๐ถ๐๐โ ๐ถ๐๐๐. (2.19) From this, the HF energy can be written as: ๐ธHF =โ๐ท๐๐โ๐๐๐๐+12โ ๐ท๐๐๐ท๐๐ (โจ๐๐|๐๐โฉ โ12โจ๐๐|๐๐โฉ)๐๐๐๐. (2.20) Operating on the expanded basis set with the Fock operator, Equations (2.16) and (2.18), followed by projection onto the basis set yields the Roothaan-Hall equations: โ(โ๐๐ +โ๐ท๐๐(๐๐โจ๐๐|๐๐โฉ โ12โจ๐๐|๐๐โฉ)โ ๐น๐๐๐ถ๐๐๐=โโจ๐๐|๐๐โฉโ ๐๐๐๐ถ๐๐๐๐๐. (2.21) Where ๐น๐๐ and ๐๐๐ are the Fock and overlap matrix elements. The resulting eigenvectors and eigenvalues represent the expansion coefficients ๐ถ๐๐ and orbital energies ๐๐ respectively. It is worth noting that only the restricted HF procedure for closed shell molecules is shown.33 With this method, the HF-SCF procedure is variational in the density matrix, during which the HF energy is minimized through the optimization of the basis set expansion coefficients. The full algorithm is shown in the Appendix, however, in short, one inputs an initial guess ๐0 for the density matrix, substituting the guess into the Roothaan-Hall equations yields a new matrix ๐1 which is typically closer to the true solution. With each iteration one is able to improve the molecular orbital coefficients, thereby reducing the energy. The procedure is said to have 18 converged when the input coefficients are equal to the output coefficients. Since the procedure is variational in nature, the current energy is always greater than the exact ground state energy: ๐ธ๐๐ฅ๐๐๐ก โค ๐ธHF. (2.22) The two-electron interaction in Equation (2.14) treats the exchange interaction correctly. However, in HF theory the motion of an electron is described based on the mean field of all other electrons, as opposed to true ๐๐๐โ1 repulsion energy. As a result of this approximation, HF theory completely neglects the correlational component of the Coulomb interaction. Correlation energy, ๐ธ๐๐๐๐ , is defined as the difference of the exact non-relativistic energy and the HF energy, ๐ธ๐๐๐๐ = ๐ธ๐๐ฅ๐๐๐ก โ ๐ธ๐ป๐น. (2.23) Since dispersion forces arise from electron correlation, systems where van der Waals interactions play a significant role, such as non-covalently bound molecular complexes and complexes with substantial intramolecular hydrogen bonding, cannot be accurately depicted by HF theory.59 One can correct for electron correlation variationally or through the application of perturbation theory. The most popular method for this is Mรธller-Plesset perturbation theory, which is discussed in the sub-section below. 2.2.2. Mรธller-Plesset Perturbation Theory In the following section, little effort has been made to explain the general procedure for many-body perturbation theory, rather, the author directs the reader to a review article in the International Journal of Quantum Chemistry as an excellent reference on the topic.60 Mรธller-Plesset perturbation theory (MPn) is a special case of Rayleigh-Schrรถdinger perturbation theory, where the zeroth-order Hamiltonian is chosen to be the HF Hamiltonian, previously denoted as the Fock operator.45 The total Hamiltonian for any perturbative method is simply the sum of an unperturbed (zeroth order) Hamiltonian, and a perturbed Hamiltonian, ?ฬ? = ?ฬ?(0) + ๐?ฬ?โฒ, (2.24) where ?ฬ?(0), ?ฬ?โฒ, and ๐ is the unperturbed, and perturbed Hamiltonian, and perturbation parameter, respectively.45 Recall that the zeroth order Hamiltonian for MPn theory is the HF Hamiltonian, ?ฬ?(0) = ?ฬ?. (2.25) Whereas the perturbed Hamiltonian, ?ฬ?โฒ, is defined as the difference between the electron-electron Coulomb and HF potentials, 19 ?ฬ?โฒ = ?ฬ? โ ?ฬ?0 =12โ1|๐ซ๐ โ ๐ซ๐|๐๐>๐โโ?ฬ?๐ป๐น(๐ซ๐)๐๐=โ(1|๐ซ๐ โ ๐ซ๐|โ ?ฬ?๐ป๐น(๐ซ๐))๐>๐, (2.26) From this it is evident that the perturbative term describes the correlation energy, and as such, the more terms considered in the perturbation, the more accurate the correlation energy will be. Additionally, it is worth mentioning that the wavefunction, ฮฆ, and energy, ๐ธ, are expanded as a perturbative series in ๐, ฮฆ๐ = lim๐โโโ๐๐ฮฆ๐๐๐๐=0, (2.27) with energy being given by, ๐ธ๐ = lim๐โโโ๐๐๐ธ๐๐๐๐=0. (2.28) Recall that the unperturbed Hamiltonian is the Fock operator, it follows that the zeroth order wavefunction for the ground state is a single HF determinant, ฮฆ0, and the excited state wavefunctions are given by their respective excited determinants, ฮฆ๐. That is, ฮฆ๐๐ corresponds to the one electron, ๐, excitation from an occupied orbital, ๐, to some higher energy unoccupied orbital. Additionally, ฮฆ๐๐๐๐ and ฮฆ๐๐๐๐๐ ๐ก represent the double excitation of electrons ๐, ๐, and triple excitation of electrons ๐, ๐, and ๐ respectively. The ground state zeroth order energy, ๐ธ0(0), is simply the sum of all occupied orbitals energies, ๐ธ0(0) =โ๐๐๐๐, (2.29) where ๐ is the number of occupied orbitals. First order corrections, known as MP1, or more commonly, HF, to the energy are easily obtained as expectation values of the HF determinant over the perturbative operator, ๐ธ0(1) = โจฮฆ0|?ฬ?โฒ|ฮฆ0โฉ = โ12โโจ๐๐โ๐๐โฉ๐๐,๐. (2.30) From standard perturbation theory, the total first order energy is simply the sum of the zeroth energy and the first order correction, ๐ธ0(0) + ๐ธ0(1) =โ๐๐๐๐โ12โโจ๐๐โ๐๐โฉ๐๐,๐= ๐ธ0๐ป๐น. (2.31) 20 A quick comparison with Equation (2.10) verifies that the first order corrected energy is equal to the HF energy. Second order corrections, known as MP2, are obtained via the expression, ๐ธ0(2) = โโ|โจฮฆ0|?ฬ?โฒ|ฮฆ๐โฉ|2๐ธ๐(0) โ ๐ธ0(0)๐โ 0, (2.32) where ฮฆ๐ denotes single, double, triple, and higher excited determinants. However, ?ฬ?โฒ is a two-body operator, and as such, any excitation with respect to ฮฆ0 involving more than two-bodies yields vanishing matrix elements. That is, one need only consider single and double excited determinants, greatly simplifying Equation (2.32). Further simplification arises from the fact that single excitations also yield vanishing matrix elements: โจฮฆ0|?ฬ?โฒ|ฮฆ๐๐โฉ = โโจ๐๐โ๐๐โฉ๐๐โ โจ๐|?ฬ?๐ป๐น|๐โฉ =โโจ๐๐โ๐๐โฉ๐๐โโโจ๐๐โ๐๐โฉ๐๐= 0. (2.33) Thus, the second order contribution to energy only requires consideration of the doubly excited determinants, ฮฆ๐๐๐๐ , which simplifies Equation (2.32) to: ๐ธ0(2) = โโโ|โจ๐๐โ๐๐ โฉ|2๐๐ + ๐๐ โ ๐๐ โ ๐๐๐๐<๐ ๐๐<๐, (2.34) where ๐ and ๐ correspond to the occupied and unoccupied orbitals respectively. Furthermore, recall the antisymmetry requirement for the integralโs states that โจ๐๐โ๐๐ โฉ = โโจ๐๐โ๐ ๐โฉ =โโจ๐๐โ๐๐ โฉ, and โจ๐๐โ๐๐ โฉ = 0 if ๐ = ๐ or ๐ = ๐ . Applying these properties to Equation (2.34) yields an even simpler expression, ๐ธ0(2) = โ14โโ|โจ๐๐โ๐๐ โฉ|2๐๐ + ๐๐ โ ๐๐ โ ๐๐๐๐,๐ ๐๐,๐. (2.35) Thus, the ground state total MP2 energy is given by the sum of all contributions, ๐ธ0๐๐2 = ๐ธ0(0) + ๐ธ0(1) + ๐ธ0(2) = ๐ธ0๐ป๐น + ๐ธ0(2). (2.36) It is worth mentioning that the second order correction to the energy is always negative, and since this method is not variational, the MP2 total energy is not necessarily greater than the true energy. Additionally, MPn theory fails for systems with no (or an infinitely small) HOMO-LUMO gap, as ๐๐ + ๐๐ โ ๐๐ โ ๐๐ = 0 causing second order corrections to diverge to negative infinity. 21 In theory, one could calculate an infinite number of corrections to the energy, yielding the true energy. However, with each additional term, the computational cost increases exponentially. Additionally, with each additional term the correction to the energy decreases. That is, if one assumes that the HF energy contains 99% of the systems energy, then second order corrections will account for 0.6%, third order corrections accounting for 0.2%, fourth order corrections account for 0.05%, and so forth. As such, due to the steep increase in computational cost, truncating after the second correction is most commonly used.55 That is, MP2 is substantially more popular than all other forms of MPn theory, such as MP3, MP4, and so forth. In this thesis calculations were only performed up to MP2, and as such, the inclusion of the additional terms would be irrelevant. 2.2.3. Coupled Cluster Theory Unlike the methods discussed in the above subsections, the coupled cluster (CC) wavefunction can trace its origins to an exponential ansatz: |ฮจ๐ถ๐ถโฉ = ๐?ฬ?|ฮจ0โฉ, (2.37) where ฮจ0 is typically the HF wavefunction, and ?ฬ? is an excitation operator typically referred to as the cluster operator. The total cluster operator is obtained from the sum of cluster operators of different excitation levels, ?ฬ?๐ =โ?ฬ?๐๐๐= ?ฬ?1 + ?ฬ?2 +โฏ+ ?ฬ?๐, (2.38) where ๐ denotes the number of electrons excited from occupied to unoccupied orbitals. That is, ?ฬ?1 considers all single electron excitations, ?ฬ?2 considers all double electron excitations, and so forth.55 The single excitation operator, in the second quantization formalism, has the following form, ?ฬ?1 =โโ๐ก๐๐๐๐?ฬ?๐โ ?ฬ?๐๐๐, (2.39) where ?ฬ?๐โ and ?ฬ?๐ are the annihilation and creation operators for the spin orbitals ๐ and ๐, respectively. Additionally, the cluster amplitudes to be determined are given by ๐ก๐๐ . Additionally, the second excitation operator in the second quantization formalism is, ?ฬ?2 = โโ๐ก๐๐๐๐๐๐<๐?ฬ?๐โ ?ฬ?๐?ฬ?๐โ ?ฬ?๐๐๐<๐. (2.40) 22 Imposing the condition that the amplitudes, ๐ก๐๐๐๐, must be antisymmetric with respect to exchange, and recalling the anticommutation property for annihilation and creation operators, ?ฬ?๐?ฬ?๐ =โ?ฬ?๐?ฬ?๐, the second excitation operator has the form, ?ฬ?2 =14โโ๐ก๐๐๐๐๐๐,๐?ฬ?๐โ ?ฬ?๐?ฬ?๐โ ?ฬ?๐๐๐,๐. (2.41) It is useful to generalize the cluster operator for up to ๐-tuple excitations, ?ฬ?๐ =1(๐!)2โ โ ๐ก๐๐๐โฏ๐๐๐โฏ๐๐,๐,๐,โฏ?ฬ?๐โ ?ฬ?๐?ฬ?๐โ ?ฬ?๐?ฬ?๐โ ?ฬ?๐โฏ๐๐,๐,๐โฏ. (2.42) Substituting Equation (2.37) into the Schrรถdinger equation yields: ?ฬ?๐ ?ฬ?|ฮจ0โฉ = ๐ธ๐?ฬ?|ฮจ0โฉ. (2.43) An expression for the correlation energy can be obtained by subtracting the HF energy from both sides of the above equation: (?ฬ? โ โจฮจ0|?ฬ?|ฮจ0โฉ)โ ?ฬ?๐๐ ?ฬ?|ฮจ0โฉ = (๐ธ โ ๐ธ๐ป๐น)โ ๐ธ๐๐๐๐๐ ?ฬ?|ฮจ0โฉ, (2.44) where ๐ธ๐๐๐๐ represents the correlation energy, and ?ฬ?๐ is the normal-ordered Hamiltonian operator. By definition, the Hamiltonian is constructed as a normal-ordered product of annihilation and creation operators, in which all annihilation operators are to the right of all creation operators, thereby the product, when applied to the vacuum state, is zero.55 Contributions from one-electron terms in the normal-ordered Hamiltonian are given by the operator ?ฬ?๐, whereas the operator ?ฬ?๐ contains the two-electron terms: ?ฬ?๐ =โ๐
๐๐{?ฬ?๐โ ?ฬ?๐}๐๐, (2.45) ?ฬ?๐ =14โโจ๐๐โ๐๐โฉ{?ฬ?๐โ ?ฬ?๐โ ?ฬ?๐?ฬ?๐}๐๐๐๐. (2.46) The one electron operator, ?ฬ?๐, is simply the Fock operator, and as such, ๐
๐๐ are the Fock matrix elements, ๐
๐๐ = โจ๐|โฬ|๐โฉ +โโจ๐๐โ๐๐โฉ๐. (2.47) Solving these equations requires substantial computational effort, and as such, the CC equations are typically not solved variationally. However, the amplitudes for each excitation and their 23 contribution to the correlational energy can be obtained via multiplication from the left of Equation (2.44) with ๐โ?ฬ?, followed by projection onto the reference and excited determinants, โ๐ธ๐๐๐๐ = โจฮจ0|๐โ?ฬ??ฬ?๐๐?ฬ?|ฮจ0โฉ (2.48) โจฮจ๐๐|๐โ?ฬ??ฬ?๐๐?ฬ?|ฮจ0โฉ = 0 (2.49) โจฮจ๐๐๐๐|๐โ?ฬ??ฬ?๐๐?ฬ?|ฮจ0โฉ = 0 (2.50) โฎ โจฮจ๐๐๐โฏ๐๐๐โฏ|๐โ?ฬ??ฬ?๐๐?ฬ?|ฮจ0โฉ = 0. (2.51) The obtained energies are identical to those obtained at the full configuration interaction (FCI) level of theory, as all possible excited determinants are included. In other words, the CC wavefunction is a nonlinear reparameterization of the FCI wavefunction, and optimization of the cluster amplitudes yields the FCI wavefunction.55 However, the interest in CC theory does not stem from solving for all of the amplitudes, after all, FCI already exists. Rather, the CC approach allows for a more natural truncation of excitation amplitudes than that observed from the CI approach. Utilizing the Baker-Campbell-Hausdorff expansion, the similarity transformed Hamiltonian, Equation (2.48), is simplified to, ๐โ?ฬ??ฬ?๐๐?ฬ? = ?ฬ?๐ + [?ฬ?๐ , ?ฬ?] +12![?ฬ?๐ , ?ฬ?]2 +13![?ฬ?๐ , ?ฬ?]3 +14![?ฬ?๐ , ?ฬ?]4, (2.52) where, [?ฬ?๐ , ?ฬ?]๐ = [[?ฬ?๐ , ?ฬ?]๐โ1, ?ฬ?] , for ๐ > 1. (2.53) The expression terminates after the fifth term since only one- and two-electron terms are present in the normal-ordered Hamiltonian. To conclude, a set of nonlinear equations have been derived from the exponential ansatz enabling one to iteratively solve for the cluster amplitudes and their associated contributions to the correlational energy.55 2.2.4. Truncating the Cluster Operator As mentioned in the previous section, interest in coupled cluster theory arises from the fact that truncation of the cluster operator comes about more naturally than that of the CI excitation operator. This is achieved by taking the Taylor series expansion the cluster operator, which, after rearranging the operators by excitation level, yields the following expression, 24 ๐ ?ฬ? = 1 + ?ฬ?1 +?ฬ?22!+ โฏ+?ฬ?๐๐!= 1 + ?ฬ?1 + ?ฬ?2 +โฏ+ ?ฬ?๐, (2.54) where ?ฬ?๐ generates all the ๐-electron excitations. Using the fact that the cluster operators commute with one another, the following expressions for ?ฬ? in terms of ?ฬ? can be obtained for the first three excitation levels: ?ฬ?1 = ?ฬ?1, (2.55) ?ฬ?2 =12?ฬ?12 + ?ฬ?2, (2.56) ?ฬ?3 =16?ฬ?13 + ?ฬ?1?ฬ?2 + ?ฬ?3. (2.57) That is, the first excitation operator generates the single excitations only through the excitation of a single electron. However, for the case of the second excitation operator, ?ฬ?2, the excitation of two electrons can be achieved through the simultaneous excitation of two electrons, ?ฬ?2, or the independent excitation of two single electrons, ?ฬ?12. The third excitation operator, ?ฬ?3, describes the excitation of three electrons, which can be accomplished via the independent excitation of three electrons, ?ฬ?13, independent excitations of one and two electrons, ?ฬ?1?ฬ?2, or through the simultaneous excitation of three electrons, ?ฬ?3.55 Using this property, the CC wavefunction can be written as a sum of individual excitation operators: |ฮจ๐ถ๐ถโฉ = |ฮฆ0โฉ +โโ๐๐๐|ฮฆ๐๐โฉ๐ฃ๐๐๐ก๐๐๐๐๐+โโ๐๐๐๐๐ |ฮฆ๐๐๐๐ โฉ๐ฃ๐๐๐ก๐<๐ ๐๐๐๐<๐+โฏ, (2.58) where ๐๐๐ are coefficients which represent the cluster amplitudes for the expanded cluster operators. The coefficients for the first three excitations are related to the cluster amplitudes through, ๐๐๐ = ๐ก๐๐ , (2.59) ๐๐๐๐๐ = ๐ก๐๐๐๐ + ๐ก๐๐ ร ๐ก๐๐ , (2.60) ๐๐๐๐๐๐ ๐ก = ๐ก๐๐ ร ๐ก๐๐ ร ๐ก๐๐ก + ๐ก๐๐ ร ๐ก๐๐๐ ๐ก . (2.61) In the above expressions, ร denotes the antisymmetric product. The use of the antisymmetric product is required to ensure the antisymmetry principle holds for the exchange of any two occupied or unoccupied spin orbitals. It is of particular interest to notice the expression for the coefficient for triplet excitations, (2.61). That is, the coefficients are the product of single- and double- excitation amplitudes. Hence, contribution from the triples term can be calculated through 25 the use of only the singles and doubles amplitudes, as opposed to having to explicitly calculate the triples amplitudes. The most common use of CC theory involves truncation of the cluster operator after the doubles term. This method is referred to as Coupled Cluster Singles and Doubles (CCSD), in which the cluster operator is simply the sum of the two individual excitation operators, which were explicitly shown in Equations (2.39) and (2.41). Solving the CC equations for the singles amplitudes yields, ๐ก๐๐๐ท๐๐ = ๐๐๐ +โ๐ก๐๐๐น๐๐๐โโ๐ก๐๐ ๐น๐๐๐+โ๐ก๐๐๐๐๐น๐๐๐๐โโ๐ก๐๐โจ๐๐โ๐๐โฉ๐๐โ12โ ๐ก๐๐๐๐โจ๐๐โ๐๐โฉ๐๐๐โ12โ ๐ก๐๐๐๐ โจ๐๐โ๐๐โฉ๐๐๐. (2.62) Equivalently for the doubles amplitudes, ๐ก๐๐๐๐๐ท๐๐๐๐ = โจ๐๐โ๐๐โฉ + ๐(๐๐)โ[๐ก๐๐๐๐ (๐น๐๐ โ12โ๐ก๐๐ ๐น๐๐๐)]๐โ ๐(๐๐)โ[๐ก๐๐๐๐ (๐น๐๐ +12โ๐ก๐๐๐น๐๐๐)]๐+12โ๐๐๐๐๐๐๐๐๐๐๐๐+12โ๐๐๐๐๐๐๐๐๐๐๐๐+ ๐(๐๐)๐(๐๐)โ(๐ก๐๐๐๐๐๐๐๐๐ โ ๐ก๐๐๐ก๐๐ โจ๐๐โ๐๐โฉ)๐๐+ ๐(๐๐)โ(๐ก๐๐โจ๐๐โ๐๐โฉ)๐โ ๐(๐๐)โ(๐ก๐๐ โจ๐๐โ๐๐โฉ)๐. (2.63) Where ๐(๐๐) is a permutation operator of the form, ๐(๐๐)ฮฉ(๐ฅ, ๐ฆ) = ฮฉ(๐ฅ, ๐ฆ) โ ฮฉ(๐ฆ, ๐ฅ), (2.64) with ฮฉ(๐ฅ, ๐ฆ) being an arbitrary function. Additionally, the ๐น intermediates are defined as follows: ๐น๐๐ = (1 โ ๐ฟ๐๐)๐๐๐ โ12โ๐ก๐๐ ๐๐๐๐+โ๐ก๐๐ โจ๐๐โ๐๐โฉ๐๐โ12โ ?ฬ?๐๐๐๐ โจ๐๐โ๐๐โฉ๐๐๐ (2.65) ๐น๐๐ = (1 โ ๐ฟ๐๐)๐๐๐ +12โ๐ก๐๐๐๐๐๐+โ๐ก๐๐โจ๐๐โ๐๐โฉ๐๐+12โ?ฬ?๐๐๐๐โจ๐๐โ๐๐โฉ๐๐๐ (2.66) ๐น๐๐ = ๐๐๐ +โ๐ก๐๐โจ๐๐โ๐๐โฉ๐๐. (2.67) Whereas the ๐ intermediates are defined as follows: 26 ๐๐๐๐๐ = โจ๐๐โ๐๐โฉ + ๐(๐๐)โ๐ก๐๐โจ๐๐โ๐๐โฉ๐+14โ๐๐๐๐๐โจ๐๐โ๐๐โฉ๐๐ (2.68) ๐๐๐๐๐ = โจ๐๐โ๐๐โฉ โ ๐(๐๐)โ๐ก๐๐ โจ๐๐โ๐๐โฉ๐+14โ๐๐๐๐๐ โจ๐๐โ๐๐โฉ๐๐ (2.69) ๐๐๐๐๐ = โจ๐๐โ๐๐โฉ +โ๐ก๐๐โจ๐๐โ๐๐โฉ๐โโ๐ก๐๐โจ๐๐โ๐๐โฉ๐โโ(12๐ก๐๐๐๐+ ๐ก๐๐๐ก๐๐) โจ๐๐โ๐๐โฉ๐๐. (2.70) Note that ๐ and ?ฬ? were used to simplify the expressions, and correspond to the following combinations of excitation amplitudes: ๐๐๐๐๐ = ๐ก๐๐๐๐ + ๐ก๐๐๐ก๐๐ โ ๐ก๐๐๐ก๐๐ (2.71) ?ฬ?๐๐๐๐ = ๐ก๐๐๐๐ +12(๐ก๐๐๐ก๐๐ โ ๐ก๐๐๐ก๐๐). (2.72) Lastly, ๐ท is the diagonal elements from the Fock matrix (i.e. the inverse orbital energies), which are defined as follows: ๐ท๐๐ = ๐น๐๐ โ ๐น๐๐ , (2.73) ๐ท๐๐๐๐ = ๐น๐๐ + ๐น๐๐ โ ๐น๐๐ โ ๐น๐๐ . (2.74) Similarly, to MP2, the amplitudes are calculated iteratively. In general, the guess for the ๐ก๐๐ amplitudes are simply zero, whereas those for the ๐ก๐๐๐๐ amplitudes are the corresponding MP2 amplitudes. Hence, prior to performing a CCSD calculation one typically calculates the MP2 amplitudes to use as an initial guess. The guess and corrected amplitudes are inserted into the right-hand side of Equations (2.62) and (2.63) yielding the next set of amplitudes.55 In general, CCSD calculations can achieve an accuracy of 0.005 ร
for bond lengths, and 50 kJ mol-1 for electronic energies.55 This is reasonable, however, for applications such as high-resolution spectroscopy, higher accuracies are sometimes required to aid in the assignment of clustered spectra. One of the most well-known methods for achieving said accuracy increase is termed CCSD(T) and is identical to CCSD with the inclusion of triples amplitudes via non-iterative perturbation theory. The CCSD(T) Hamiltonian can be written in a way analogous to that for MP2, and the triples amplitudes are simply calculated using the exact same formula as that for MP4, however, the converged CCSD amplitudes are used instead of the perturbation coefficients. Corrections to the energy arising from the triples contributions are then the fourth and fifth order 27 terms from the perturbative corrections, describing the coupling between singles and triples. The triples corrections to the energy can simply be added to the CCSD energy, ๐ธCCSD(T) = ๐ธCCSD + ๐ธ(4) + ๐ธ(5). (2.75) Where the corrections are given by the following expressions: ๐ธ(4) =136โ ๐ก๐๐๐๐๐๐โ๐ท๐๐๐๐๐๐๐ก๐๐๐๐๐๐๐๐๐๐๐๐ (2.76) ๐ธ(5) =14โ โจ๐๐โ๐๐โฉ๐ก๐๐โ๐ก๐๐๐๐๐๐.๐๐๐๐๐๐ (2.77) Where ๐ท๐๐๐๐๐๐ are the diagonal elements of the Fock matrix. In a way analogous to Equations (2.62) and (2.63), the triples amplitudes can be obtained via, ๐ท๐๐๐๐๐๐๐ก๐๐๐๐๐๐ = ๐(๐|๐๐)๐(๐|๐๐) [โโจ๐๐โ๐๐โฉ๐ก๐๐๐๐๐โโโจ๐๐โ๐๐โฉ๐ก๐๐๐๐๐]. (2.78) Where the permutation operator is defined such that ๐(๐|๐๐)๐(๐๐๐) = ๐(๐๐๐) โ ๐(๐๐๐) โ ๐(๐๐๐). The widespread popularity of CCSD(T) can be attributed to the high-level of accuracy that it is capable of achieving, with a reasonable increase in computational cost.55 Unfortunately, as one includes higher excitation amplitudes, the methods dependence on basis set size increases exponentially. As such, one must employ at least a triple zeta basis set in order for CCSD(T) calculations to be reliable, with larger quadrupole zeta basis sets being preferred. This basis set requirement makes it such that only molecules of up to approximately 30 atoms can be calculated at the CCSD(T) level of theory.55 2.3. Quantum Chemistry: Excited State Calculations As discussed in the section above, quantum chemistry is most commonly applied to calculate the ground state properties of molecules, by solving the electronic Schrรถdinger equation for the lowest-energy solution of the electronic Hamiltonian operator. Solutions to the higher roots of the operator relate to the electronically excited states, and unfortunately, require the use of substantially more complex algorithms. The electronic wavefunction for molecules in the ground state is typically dominated by a single configuration, and as such, can be accurately represented by a single Slater determinant (Section 2.2.1). Excited states on the other hand, require one to examine the wavefunction to confirm whether it can be treated with a single-reference method, such as those mentioned in the paragraph above, or if they require the use of a multi- 28 configurational technique. Additionally, one must determine whether the excited state is a bound or dissociative state, which is an entirely different problem in and of itself. Unlike the previously discussed ground state calculations, calculations for electronically excited states are much more complex, and the results require a careful analysis. Vertical excitation energies, adiabatic excitation energies, geometry optimizations, and vibrational frequencies were performed for ๐ผ-alanine. A few different well-known approximations were used to ensure the reproducibility of results, those being the equation of motion coupled cluster (EOM-CC)61 and symmetry adapted cluster-configuration interaction (SAC-CI)62 theories. At the time this thesis is being written, only the EOM-CC calculations have finished, and as such, only EOM-CC will be discussed in detail. Additionally, a brief introduction to configuration interaction theory will be included.63 2.3.1. Configuration Interaction-Singles The basis of most single-configurational techniques for examining excited states is the configuration interaction-singles (CIS) approach. CIS, just like HF for the ground state, is the computationally and conceptually simplest method available for the calculation of excited state properties.55 One must first optimize the ground state Slater determinant, Equation (2.9), through the HF procedure. Here we will denote this determinant as, ฮฆ0(r) = |โ
1(๐)โ
2(๐)โฏโ
๐(๐)|, (2.79) where โ
๐(๐) refers to the doubly occupied spatial orbitals. Insertion of this wavefunction into the time-independent HF equation yields the optimized occupied and virtual molecular orbitals (Section 2.2.1). In this section, indices ๐, ๐, ๐ refer to occupied orbitals, ๐, ๐, ๐ refer to virtual orbitals, and ๐, ๐, ๐ refer to general orbitals. In the CIS method the excited Slater determinants are constructed via replacing occupied orbitals with virtual ones in the ground state Slater determinant. For example, replacing any occupied orbital ๐ with any virtual orbital ๐, then the single excited Slater determinant is denoted as ฮฆ๐๐(๐). Note that one must only include single excitations. The electronic wavefunction can now be constructed as a linear combination of these excited state determinants along with the ground state wavefunction through the following ansatz, 29 ฮจCIS =โ๐๐๐ฮฆ๐๐(๐)๐๐. (2.80) Substitution into the time-independent Schrรถdinger equation with the electronic Hamiltonian operator, followed by projection from the left onto the singly excited determinants, โจฮฆ๐๐|, gives the following expression: โโจฮฆ๐๐|?ฬ?|ฮฆ๐๐โฉ๐๐๐๐๐= ๐ธCISโ๐๐๐๐ฟ๐๐๐ฟ๐๐๐๐. (2.81) From the above equation, one can easily obtain an expression for the excitation energies, โ{(๐๐ โ ๐๐)๐ฟ๐๐๐ฟ๐๐ + โจ๐๐โ๐๐โฉ}๐๐๐๐๐= ๐CISโ๐๐๐๐ฟ๐๐๐ฟ๐๐๐๐, (2.82) which utilizes the relation: โจฮฆ๐๐|?ฬ?|ฮฆ๐๐โฉ = (๐ธ0 + ๐๐ โ ๐๐)๐ฟ๐๐๐ฟ๐๐ + โจ๐๐โ๐๐โฉ. (2.83) Recall that ๐๐ and ๐๐ are the orbital energies for the single-electron orbitals โ
๐ and โ
๐, and โจ๐๐โ๐๐โฉ represents the antisymmetrized two-electron integrals (defined in Section 2.2.1). Just as we did for the HF method, converting Equation (2.83) into a matrix eigenvalue equation: ๐๐ = ๐๐, (2.84) where ๐ is the Hamiltonian matrix projected onto the space of the singly excited determinants, ๐ is the diagonal matrix of excitation energies, and ๐ is the matrix of CIS expansion coefficients, ๐๐๐. By definition, the matrix elements of the Hamiltonian are given by, ๐๐๐,๐๐ = (๐๐ โ ๐๐)๐ฟ๐๐๐ฟ๐๐ + โจ๐๐โ๐๐โฉ. (2.85) Thus, diagonalization of ๐, (๐ โ ๐)๐ = 0, yields the excitation energies. In particular, just as was true for HF, the eigenvalues and eigenvectors from the diagonalization correspond to the excitation energies and expansion coefficients, respectively.55 2.3.2. Equation of Motion Coupled Cluster Theory EOM-CC is a well-known single reference technique that can be applied to the calculation of excited, electron attached, and ionization states. The excited state wave functions, ฮจ๐ธ๐, are obtained from the ground state wave function, ฮจ๐บ๐, through the use of a linear excitation operator, |ฮจ๐ธ๐โฉ = โฬ|ฮจ๐บ๐โฉ, (2.86) where โฬ denotes the total linear excitation operator.61 For excitation energies, โฬ is a CI-like particle conserving operator, however, for the ionization and electron attachment cases it is simply 30 a net annihilation or net creation operator, respectively.61 However, since we are only concerned with the excitation energies no details regarding the ionization and attachment cases will be discussed. For an excited state, ๐, โฬ๐ has the form, โฬ๐ =1๐!2โ๐๐๐๐โฏ๐๐๐๐โฏ๐โ ๐๐โ ๐๐โ ๐๐โ โฏ , (2.87) where ๐, ๐, ๐, ๐ โฏ and ๐, ๐, ๐ โฏ denote the unoccupied and occupied orbitals respectively. The total linear excitation operator is simply the sum of the individual operators, โฬ =โโฬ๐๐. (2.88) Recall that the ground state coupled cluster wavefunction was discussed in detail in Section 2.2.3 of this thesis and can be seen in Equation (2.58). Substituting the ground state coupled cluster wavefunction and Equation (2.86) into the Schrรถdinger equation yields, ?ฬ?โฬ๐ ?ฬ?|ฮฆ0โฉ = ๐ธโฬ๐?ฬ?|ฮฆ0โฉ = (๐โ?ฬ??ฬ?๐?ฬ? โ ๐ธ)โ ?ฬ?๐๐๐โฬ|ฮฆ0โฉ = 0. (2.89) The simplification of Equation (2.89) arises from the fact that โฬ and ?ฬ? must commute.61 If one denotes the effective Hamiltonian ?ฬ?๐๐๐ as ?ฬ?๐๐๐ = ๐โ?ฬ??ฬ?๐ ?ฬ? โ ๐ธ, then Equation (2.89) presents itself as a typical eigenvalue problem, ?ฬ?๐๐๐โฬ|ฮฆ0โฉ = 0. (2.90) EOM-CC is a formally exact method, meaning that the similarity transformation presented above does not change the eigenvalues of the full Hamiltonian.61 However, just as is true for standard coupled cluster theory, the treatment of all excitations is infeasible, and as such, the operator must be truncated. This is done in a fashion identical to that discussed in Section 2.2.4. Most commonly, single and double excitations are considered. This is known as the EOM-CCSD approximation, which is the equation of motion equivalent to CCSD. Just like CCSD, EOM-CCSD provides an exact treatment of two-electron systems, as no higher order excitations are possible (3 electrons cannot be excited if only 2 electrons exist). Expressing the effective Hamiltonian in the spin-orbital basis yields, ?ฬ?๐๐๐ =โ๐น๐๐๐โ ๐๐๐+14โ ๐๐๐๐๐ ๐โ ๐โ ๐ ๐๐๐๐๐ +โฏ, (2.91) where ๐น and ๐ are the one- and two-body amplitudes respectively, whose leading terms correspond to the Fock matrix elements and two-electron integrals, with remaining contributions 31 coming from the appropriate contractions between the Hamiltonian and the ?ฬ? amplitudes that parametrize the ground state wavefunction. Unlike the standard Hamiltonian, the effective Hamiltonian contains three- and higher-body operators. The effective Hamiltonian takes on the following matrix form, ?ฬ?๐๐๐ = (๐ธ๐ถ๐ถ โฏ โฏ0 โฑ โฏ0 โฎ โฑ), (2.92) where the zero terms arise from the condition that ?ฬ? amplitudes obey the coupled cluster equations, as presented in Section 2.2.3. Hence, the ground state coupled cluster solution can be obtained from the root of ?ฬ?๐๐๐ corresponding to โฬ = โฬ0 = 1. Additionally, since ?ฬ?๐๐๐ is not a Hermitian operator, each root is associated with two eigenvectors, โฬ and โฬ, which correspond to the distinct states, โจฮจฬ| = โจฮฆ0|โฬ๐?ฬ? , (2.93) and, |ฮจโฉ = ๐ ?ฬ?โฬ|ฮฆ0โฉ. (2.94) โฬ is the left eigenvector of ?ฬ?๐๐๐, and is a deexcitation operator, which has the same form as the excitation operator โฬ. That is, the individual and total operator is given by, โฬ๐ =1๐!2โ๐๐๐๐โฏ๐๐๐๐โฏ๐โ ๐๐โ ๐๐โ ๐๐โ โฏ, (2.95) and, โฬ =โโฬ๐๐. (2.96) Unfortunately, due to the non-Hermitian effective Hamiltonian and intermediate normalization, these states are never Hermitian conjugates. It is of interest to note that excitation energies are the eigenvalues for the excitation operator, and as such, one can neglect the deexcitation terms when calculation the excitation energies. However, for the calculation of all other properties, such as transition densities, the excitation and deexcitation operators must both be evaluated. As a result of the non-Hermitian nature of this eigenvalue problem, we can make the two sets of solutions orthogonal by multiplying them by a unity matrix, โจฮจฬ๐|ฮจ๐โฉ = ๐๐ฟ๐๐ , (2.97) 32 which is the EOM-CC normalization condition. It follows to write the energy in its functional form, ๐ธ = โจฮจฬ|?ฬ?|ฮจโฉ = โจฮฆ0|โฬ?ฬ?๐๐๐โฬ|ฮฆ0โฉ. (2.98) Hellmann-Feynman theorem states that the derivative of the stationary values of a functional with respect to a parameter on which it depends, can be calculated from the explicit dependence, disregarding any and all implicit dependence. Hence, due to the stationary nature of the solutions with respect to variation of the โฬ and โ vectors, one can write, ๐ = โจฮจฬ|ฮฬ|ฮจโฉ, (2.99) where ฮฬ is an arbitrary linear operator. Using Equation (2.99) as a generalized expectation value, we can calculate properties in the usual way, ๐ = Tr(๐ฮ), (2.100) where ๐ is the reduced n-particle density matrix, ๐๐๐โฏ๐๐ = โฉ๐โ ๐โ โฏ๐๐ โช = โจฮฆ0|โฬ(๐โ ๐โ โฏ๐๐ ๐ ?ฬ?)๐โฬ|ฮฆ0โฉ, (2.101) where the subscript c indicates that contractions between the quasiparticle creation and annihilation operators and the T amplitudes are restricted to connected diagrams. Elements of the reduced one-particle density matrix suitable for calculating one-electron properties are โฉ๐โ ๐โช (๐๐๐), โฉ๐โ ๐โช (๐๐๐), โฉ๐โ ๐โช (๐๐๐), and โฉ๐โ ๐โช (๐๐๐), where the shorthand notation is reported in the parentheses. It is worth noting that the density matrix is only Hermitian when the excitation and deexcitation states are adjoint. Interestingly, since this condition is always satisfied at the exact limit, the non-Hermiticity of the density matrix acts as a measure of the treatment of electron correlation. 2.4. Thermodynamics in Molecular Spectroscopy For flexible molecules with more than one stable configuration, it is useful to compare their relative energies. This is due to the fact that, in the absence of thermal contributions, the relative abundance of each conformer is directly related to their place in the relative energy order. That is, the lowest energy conformer will be the most abundant, with the second lowest being the second most abundant, and so forth.64 Hence, for experiments performed in the absence of significant thermal sources, such as spectroscopic experiments utilizing seeded supersonic expansions for sample injection, the relative electronic energies greatly simplify spectral assignment by enabling 33 researchers to better predict which conformers will be observable.64,65,66 However, for experiments involving temperature variation, such as those presented in this thesis, the absence of temperature dependence in the electronic energy makes it a poor predictor of relative abundance.67 For experiments utilizing a variation in temperature, the relative Gibbs free energy, โ๐บ๐, is much more insightful than the relative electronic energy. As a result of the temperature-dependent nature of the Gibbs free energy, the relative abundance of conformers can be calculated at multiple temperatures. This arises due to the fact that the Gibbs free energy takes the vibrational, rotational, translational, and electronic degrees of freedom into account, giving a more accurate prediction of the conformers relative energies.54 The following subsections highlight the steps required to calculate the relative electronic energies, and relative Gibbs free energies associated with the investigation of ๏ก-serine. 2.4.1. Relative Electronic Energy These days, the calculation of a molecules electronic energy is relatively straightforward and can be done with a number of different electronic structure theory techniques (see above Sections). The total electronic energy, ๐ธ0, is obtained from solving the Schrรถdinger equation through an approximate method, which gives the bottom of the well energy. One can account for the zero-point vibrational energy (ZPE) effects by simply adding the ZPE to the total electronic energy, ๐ธ๐๐๐ธ = ๐ธ0 + ๐๐๐ธ. (2.102) The ZPE corrected relative energy can then be calculated by simply subtracting the energy of the most stable conformer. That is, โ๐ธ๐๐๐ธ,๐ = ๐ธ๐๐๐ธ,๐ โ ๐ธ๐๐๐ธ,1, (2.103) where ๐ธ๐๐๐ธ,๐ is the zero-point corrected electronic energy of the ith conformer, and conformer 1 is the conformer with the lowest electronic energy. 2.4.2. Molecular Partition Functions The Gibbs free energy at temperature T can be calculated from the following expression, ๐บ๐ = ๐ธ0 + ๐บ๐๐๐๐ , (2.104) where ๐ธ0 is the total electronic energy and ๐บ๐๐๐๐ is the thermal corrections to the Gibbs free energy.68 Thermal corrections to the Gibbs free energy, ๐บ๐๐๐๐ , arise directly from the thermal corrections to enthalpy and entropy: 34 ๐บ๐๐๐๐ = ๐ป๐๐๐๐ โ ๐๐๐๐๐๐ (2.105) where ๐ป๐๐๐๐ and ๐๐๐๐๐ are the internal thermal energy contributions to enthalpy and entropy, respectively.54 Thermal corrections to enthalpy can be calculated with the expression, ๐ป๐๐๐๐ = ๐๐๐๐๐ + ๐๐ต๐, (2.106) where ๐๐๐๐๐ is the thermal correction to the internal thermal energy, and ๐๐ต is the Boltzmann constant. The thermal corrections to entropy and internal thermal energy can be directly calculated from the molecular partition functions using the following expressions, ๐๐๐๐๐ = ๐๐ก๐๐๐๐ + ๐๐๐๐ก + ๐๐ฃ๐๐ + ๐๐๐ , (2.107) and, ๐๐๐๐๐ = ๐๐ก๐๐๐๐ + ๐๐๐๐ก + ๐๐ฃ๐๐ + ๐๐๐ . (2.108) Each term in the above equations corresponds to a molecular partition function. These individual contributions can be directly calculated from the partition functions via: ๐๐๐๐๐ = ๐๐๐ต๐2 (๐ ln(๐)๐๐)๐, (2.109) and, ๐๐๐๐๐ = ๐
(ln(๐๐) + ๐ (๐ ln(๐)๐๐)๐), (2.110) where ๐
is the ideal gas constant (๐๐ด๐๐ต = ๐
), ๐ is the number of moles (๐ = ๐๐๐ด), ๐ is Eulerโs constant, ๐๐ด is Avogadroโs number, and ๐ is the product of the individual contributions from each partition function, ๐ = ๐๐ก๐๐๐๐ ๐๐๐๐ก๐๐ฃ๐๐๐๐๐ . (2.111) One must now calculate the individual contributions from each partition function. Firstly, the translational partition function, ๐๐ก๐๐๐๐ , arises from the movement (translation) of the molecule, and as such, is dependent on the temperature and mass of the molecule: ๐๐ก๐๐๐๐ = (2๐๐๐๐ต๐โ2)32โ๐ = (2๐๐๐๐ต๐โ2)32โ ๐๐ต๐๐โ ๐ผ๐๐๐๐ ๐ต๐โ๐๐ฃ๐๐๐ข๐. (2.112) Ideal behavior assumes that the system behaves ideally, that is, ๐๐ = ๐๐
๐ = (๐ ๐๐ดโ )๐๐ด๐๐ต๐, hence, ๐ = ๐๐ต๐ ๐โ . It should be noted that in Equation, ๐ is the total mass of the molecule, and โ is Planckโs constant. Taking the partial derivative of the translational partition function with respect to temperature at a constant volume, yields the following expression: 35 (๐ ln(๐๐ก๐๐๐๐ )๐๐)๐=32๐. (2.113) By substituting the above two equations into the conversion relations, Equation (2.109) and (2.110), one obtains expressions for the translational contribution to entropy and internal thermal energy, respectively, ๐๐ก๐๐๐๐ = ๐
(ln(๐๐ก๐๐๐๐ ๐) + ๐ (32๐)) = ๐
(ln(๐๐ก๐๐๐๐ ) +52), (2.114) and, ๐๐ก๐๐๐๐ = ๐๐ด๐๐ต๐2 (๐ ln(๐๐ก๐๐๐๐ )๐๐)๐=32๐
๐. (2.115) Hence, one can easily obtain the contributions from translation by substituting the known values into the above two expressions. Similarly, to the translational partition function, the rotational partition function describes the rotational motion of the molecule and is given by the following expression: ๐๐๐๐ก =๐1 2โ๐๐(๐3 2โ(ฮ๐,๐ฅฮ๐,๐ฆฮ๐,๐ง)1 2โ), (2.116) where ๐๐ is the rotational symmetry number, and ฮ๐,๐ is the rotational temperature. The rotational temperature can easily be determined from the point group of the molecule. Alternatively, it is equal to the number of times on can rotate the molecule about the molecular axis to obtain an equivalent orientation. Rotational temperature on the other hand, arises from the moment of inertia. That is, the rotational temperature is given by, ฮ๐,๐ =โ28๐2๐ผ๐๐๐ต, ๐ = ๐ฅ, ๐ฆ, ๐ง (2.117) where ๐ผ๐ is the moment of inertia with respect to the molecular axis. Notice that the rotational energy levels are given by Hence, for linear polyatomic molecules, the rotational partition function simplifies to ๐๐๐๐ก =1๐๐(๐ฮ๐), and for single atoms, it further simplifies to ๐๐๐๐ก = 1. Continuing on with the general case, the rotational contribution to entropy is obtained by substituting the partition function from (2.116) into (2.110), ๐๐๐๐ก = ๐
(ln(๐๐๐๐ก) + ๐ (๐ ln(๐๐๐๐ก)๐๐)๐) = ๐
(ln(๐๐๐๐ก) +32). (2.118) Substituting (2.116) into (2.109) yields rotational contribution to internal thermal energy, 36 ๐๐๐๐ก = ๐๐๐ต๐2 (๐ ln(๐๐๐๐ก)๐๐)๐=32๐
๐. (2.119) Since there are three rotational degrees of freedom, it can be said that the average contribution from each rotational degree of freedom is ๐
๐ 2โ . The concept of the normal mode can about to describe the vibrational motion of a molecule. Hence, it is only natural that the total vibrational partition function is composed of partition functions for each normal mode. Under the harmonic approximation, the total vibrational partition function is the sum of the contributions from each normal mode: ๐๐ฃ๐๐ =๐โโ๐๐ธ0 ๐๐ต๐โโ 1 โ exp(โโ๐๐๐ ๐๐ต๐โ )๐๐=1 , (2.120) where ๐๐ is the fundamental energy associated with mode i, and ๐ธ0 is the zero-point vibrational energy (ZPE). It is useful to express the partition function in terms of vibrational temperature, ฮ๐ฃ,๐ = โ๐ฃ๐ ๐๐ตโ , ๐๐ฃ๐๐ =๐โฮ๐ฃ,๐ 2๐โโ 1 โ exp(โฮ๐ฃ,๐ 2๐โ )๐๐=1, (2.121) where ฮ๐ฃ,๐ is the vibrational temperature of the ith mode (ฮ๐ฃ,๐ = โ๐ฃ๐ ๐๐ต)โ . In the above two expressions, the zero-point reference for energy is set as the bottom of the well. Interestingly, when calculating the vibrational contribution to entropy and internal thermal energy, one must take the derivative with respect to temperature, and as such, it does not matter if one uses the bottom of the well or zero-point energy as the zero-point. That is, the bottom of the well reference and zero-point reference yield the same contributions to entropy and enthalpy. However, for the sake of completion, when ZPE is taken as the reference zero-point the vibrational partition function has the form, ๐๐ฃ๐๐ =โ11 โ exp(โฮ๐ฃ,๐ 2๐โ )๐๐=1. (2.122) Naturally, it follows that the vibrational components for entropy and internal thermal energy are: ๐๐ฃ๐๐ = ๐
(ln(๐๐ฃ๐๐) + ๐ (๐ ln(๐๐ฃ๐๐)๐๐)๐)= ๐
โ(ฮ๐ฃ,๐ ๐โ๐ฮ๐ฃ,๐ ๐โ โ 1โ ln(1 โ ๐โ๐ฉ๐ฃ,๐ ๐โ ))๐, (2.123) and, 37 ๐๐ฃ๐๐ = ๐๐๐ต๐2 (๐ ln(๐๐ฃ๐๐)๐๐)๐= ๐
โ(ฮ๐ฃ,๐2+ฮ๐ฃ,๐๐โฮ๐ฃ,๐ ๐โ โ 1)๐. (2.124) It is worth noting that one must carefully inspect the low frequency normal modes, which correspond to internal rotations, and ensure that they are treated properly. This is especially important for symmetric top molecules, in which the hindered rotor approximation must be applied. Additionally, anharmonic corrections to the vibrational partition function can be accounted for through the use of simple perturbation theory (SPT).69 In this case, the harmonic vibrational partition function is still utilized, however, the ZPE and fundamental frequencies are calculated with vibrational perturbation theory (VPT), typically truncated at second order (VPT2). Within the harmonic approximation, vibrational contribution to entropy results in discrepancies around 1.5 and 0.5 J mol-1 K-1 for molecules with and without hindered internal rotations.70 Accounting for anharmonic corrections reduces these discrepancies to approximately 0.5 and 0.2 J mol-1 K-1, respectively.70 It is interesting to note that the improvement in accuracy largely results from the anharmonic ZPE, as opposed to the anharmonic corrections to each fundamental transition. Electronic motion within a molecule is characterized by the electronic partition function, ๐๐๐, which is given by, ๐๐๐ =โ๐๐๐(โ๐ธ๐๐๐ต๐โ )๐, ๐ = 0, 1, 2, โฆ (2.125) where ๐๐ and ๐ธ๐ are the degeneracy and energy of the ith electronic state. Making the assumption that all excited electronic states have excitation energies much greater than ๐๐ต๐, such that they are inaccessible at any temperature, reduces the number of excited states to zero. Additionally, if the ground state energy, ๐ธ0, is set to zero, then the electronic partition function further simplifies to, ๐๐๐ = ๐0, (2.126) which is equal to one for non-degenerate ground states. This assumption is reasonable for molecules at low temperatures, such as those presented in this thesis. The electronic contribution to entropy can now be determined by substituting (2.126) into (2.110), ๐๐๐ = ๐
(ln(๐๐๐) + ๐ (๐ ln(๐๐๐)๐๐)๐) = ๐
ln(๐๐๐) = ๐
ln(๐0). (2.127) 38 Notice that the above expression is equal to zero for molecules with non-degenerate ground states. A quick inspection of the electronic partition function indicates that there are no temperature dependent terms, and as a result of this, the contribution to internal thermal energy is zero, ๐๐๐ = ๐๐๐ต๐2 (๐ ln(๐๐๐)๐๐)๐= 0. (2.128) Finally, the thermal corrections to entropy and internal thermal energy are given by the following expressions: ๐๐๐๐๐ = ๐
[ln(๐๐ก๐๐๐๐ ร ๐๐๐๐ก ร ๐0) +โ(ฮ๐ฃ,๐ ๐โ๐ฮ๐ฃ,๐ ๐โ โ 1โ ln(1 โ ๐โ๐ฉ๐ฃ,๐ ๐โ ))๐+ 4], (2.129) and, ๐๐๐๐๐ = ๐
[3๐ +โ(ฮ๐ฃ,๐2+ฮ๐ฃ,๐๐โฮ๐ฃ,๐ ๐โ โ 1)๐], (2.130) assuming ideal behavior. The relative Gibbs free energy for conformer i at temperature T can now be calculated with the expression, โ๐บ๐,๐ = ๐บ๐,๐ โ ๐บ๐,1, (2.131) where ๐บ๐,๐ and ๐บ๐,1 are the Gibbs free energies for molecular conformers i and 1 at temperature T. 2.4.3. Boltzmann Distribution Law The Boltzmann distribution law states that for a system at temperature ๐, the probability of a state with energy ๐ธ๐ is, ๐(๐ธ๐) =๐โ๐ธ๐ ๐๐ต๐โโ ๐โ๐ธ๐ ๐๐ต๐โ๐, (2.132) where j is the total number of occupied states. As noted at the start of Section 2.4, we are primarily interested in the Gibbs free energy, and as such, we calculate the Boltzmann distribution using the Gibbs free energy. It is of interest to note that the Boltzmann distribution law says that the population of higher energy states increases with temperature. Therefore, by varying the sublimation temperature, one changes the conformational distribution of the molecule, further simplifying the spectral assignment procedure by enabling the cross-examination of spectra to identify the change in spectral intensity for a given transition with respect to sublimation temperature.67 39 3. Results and Discussion 3.1. Conformational Dynamics of ๐ถ-serine in Solid Parahydrogen Matrices The physicochemical properties of ๐ผ-serine were investigated through the use of solid parahydrogen MI-FTIR spectroscopy. This application of MI-IR spectroscopy for conformational studies is made possible by the fact that the conformational distribution at the sublimation temperature is preserved upon deposition into the matrix.71 However, in order to accurately predict the conformational distribution at some temperature, one must also know all the barriers to conversion for the molecular species. That is, the barrier for a given conformer to convert into a lower energy conformer through the rotation of a torsion. If the thermal energy of the environment is higher than this barrier, then the conformer will freely convert to the lower energy species, and not be observable. As a rule of thumb, it is said that if the barrier to conversion between two conformers is less than 5 kJ molโ1, then the higher energy conformer will tunnel to the lower energy species.71,72 For highly flexible molecules, such as ๐ผ-serine, calculating all possible barriers to conversion would require the consideration of hundreds of saddle-points, up to fifth order. As such, it is not plausible to calculate all of the barriers to conversion. Rather, one crudely estimates the conformational distribution based on the Boltzmann distribution, which is calculated from the Gibbโs free energy of the conformers with the understanding that all of the low energy conformers may not be stable enough to be observed. Hence, careful assignment of the experimental spectrum is required. For flexible molecules, it is expected that the experimental spectrum will contain many conformers, and as such, will be very crowded. In order to accurately assign vibrational bands, the use of highly accurate anharmonic vibrational calculations is required. These calculations enable one to accurately assign the low-intensity features for the most abundant conformers, further increasing confidence in the assignment of lower-abundant conformers, for which only the strongest fundamental transitions will be observed. Another common technique to improve confidence in the assignment is to experimentally shift the conformational ratio correlating observed changes in the relative intensity. This is possible since the intensity of transitions from a given conformer will always change uniformly. One method to shift the conformational ratio involves selective irradiation with a near-infrared laser to induce conformational changes and is generally referred to as a pump-probe experiment.73 Other methods include annealing the matrix, 40 in situ broad-band or selective UV irradiation, and variation of the sublimation temperature.74,75,76 The latter method is utilized here, and operates on the principle that the conformational distribution can be calculated from the relative Gibbโs free energy, which is temperature dependent. As such, changing the temperature will result in a change in the conformational distribution. In particular, variation of sublimation temperature was used to alter the conformational population of ๐ผ-serine conformers. This section aims to assign the experimental mid-IR spectrum of ๐ผ-serine in solid parahydrogen matrices at various sublimation temperatures. The first subsection will be dedicated to a detailed analysis of the spectrum taken with a sublimation temperature of 420 K. Whereas the following subsection will include a comparison of the assigned bands with those observed at 400, 410, and 430 K. 3.2. Experimental and Computational Parameters for ๐ถ-serine The basis of the experimental techniques and computational methods are described in detail in Section 2 of this thesis. Rather, this section details the exact experimental conditions and computational methods employed for the analysis of ๐ผ-serine. Analysis of the experimental data can be found in the following subsections. The following experiments were primarily completed by Alethea Lee and Shin Toh during the summer of 2016. Deposition of the ๐ผ-serine sample (L-serine, 99%, and D,L-serine, >98%, Millipore Sigma) was performed with sublimation temperatures of 400, 410, 420, and 430 K (Figure B.1 to Figure B.4). Spectra were registered at 4 K in the 700 โ 4800 wavenumber range with a resolution of 0.1 cm-1 over 1000 scans. Spectral band intensities, widths, and areas of overlapping peaks were obtained through the use of Gaussian and Lorentzian fitting functions. The most intense serine peaks were observed with a sublimation temperature of 420 K, and as such, this spectrum was used for the primary spectral assignment, with all other spectra being used for verification of said assignment. Interestingly, the 430 K spectrum displayed significant spectral features consistent with thermal decomposition. Additionally, the intensity of characteristic serine peaks in the 410 and 430 K deposition spectra are very similar. Unfortunately, due to the complexity of the experimental spectra, the exact concentration of ๐ผ-serine in the matrix is unknown for all experiments. 41 The 15 lowest energy conformers of ๐ผ-serine, as depicted below in Figure 3.1, were selected from both, the works of Gronert and OโHair, and He and Allen.16,18 Both works were considered when choosing the conformers due to the necessity of considering all low energy configurations. The nomenclature utilized in this thesis follows no logical order, rather, it is based upon the reference works from which the initial structures were obtained. Due to the large number of studies on ๐ผ-serine, a conversion table for the nomenclature presented in studies closely related to this work can be found in the appendix. The optimized geometry and anharmonic vibrational frequencies were calculated for all 15 conformers at the DFT and MP2 levels of theory, with the Gaussian 16 suite of programs.77 For the sake of consistency, all calculations were performed for the L-enantiomer of ๐ผ-serine, albeit, within the nonrelativistic limit, said calculations are identical to those completed for the D-enantiomer. Additionally, the frozen core approximation was utilized in all calculations. The Becke 3-parameter Lee-Yang-Parr hybrid functional (B3LYP)36,37 with the inclusion of Grimmeโs empirical atom-pairwise dispersion corrections, D3,78 and Becke-Johnson damping, BJ,79 was employed for all DFT calculations. Note that the above method is commonly denoted as B3LYP-D3BJ. All DFT calculations were performed with Alhrichs triple zeta property optimized split-valence basis set with additional polarization and diffuse functions, def2-TZVPD.80,81 Due to the floppy nature of ๏ก-serine, a pruned integration grid consisting of 175,974 for hydrogen and 250,974 points for carbon, nitrogen, and oxygen was used. MP2 geometry optimizations were performed with Dunningโs augmented correlation consistent triple zeta split-valence basis set with additional polarization and diffuse functions, aug-cc-pVTZ.82 Whereas the MP2 anharmonic frequencies were calculated with the slightly smaller correlation consistent triple zeta split-valence basis set with additional polarization functions, cc-pVTZ, which was also introduced by Dunning.82 It should be noted that the anharmonic frequencies were calculated using geometries optimized with the cc-pVTZ basis set. The smaller basis set was utilized due to the exponential scaling associated with the MP2 method. That is, calculation of the MP2 gradient scales as ๐2๐3, where ๐ and ๐ are the number of occupied and total basis functions respectively. However, the analysis of noncovalent interactions within ๐ผ-serine requires the use of diffuse basis functions, which is why the geometries were optimized with the larger aug-cc-pVTZ basis set. 42 The Berny algorithm developed by Bernhard Schlegel,83 was used for all DFT and MP2 geometry optimizations. Depurturbed second order vibrational perturbation theory (DVPT2) was utilized for the numerical differential of analytical second derivatives at 72 nuclear displacements (two displacements per normal mode).84 Additionally, anharmonic vibrational-rotational couplings, and anharmonic IR intensities were calculated within the DVPT2 framework.85,86,87 Convergence thresholds were set to 3.7074 ร 10โ8 for the gradient and Hessian, 3.7074 ร 10โ5 for the cubic and quartic force fields, 1.0 ร 10โ3 for Coriolis coupling, and 1.0 ร 10โ4ร
2 for the inertia. Note the higher convergence requirements for the gradient and Hessian compared to the cubic and quartic force fields. The cubic and quartic force fields are calculated via the numerical differentiation of the Hessian. As such, since the gradient and hessian are analytically calculated, it is optimal to set tighter convergence thresholds to minimize errors observed as a result of said numerical differentiation. Coriolis coupling and inertia refer to the calculation of vibrational-rotational coupling terms. The maximum frequency for 1-2 Fermi resonances, 1-1, and 2-2 Darling-Dennison resonances were set to 200 cmโ1 and 100 cmโ1, respectively. A maximum tolerance between VPT2 and variational results of 1 cmโ1 was allowed. A minimum value of 1 cmโ1 for off-diagonal 1-1 Darling-Dennison terms was utilized, whereas the minimum for 2-2 Darling-Dennison terms was 10 cmโ1. Accurate energies were calculated with the CCSD(T)-F1288 explicitly correlated coupled-cluster (CC) technique available within the ORCA 4 suite of programs.89,90 These calculations utilized the explicitly correlated cc-pVTZ-F12 basis set,91 along with the complementary auxiliary basis set (CABS), cc-pVTZ-F12/CABS,92 and auxiliary cc-pVQZ/C basis set.93 It should be noted that the resolution of identity (RI) approximation, also known as density fitting, was implemented for the integral transformation during the explicitly correlated portion of the calculation, which drastically reduces computational time.94 These energies were calculated using geometries optimized at the MP2/aug-cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory, which were discussed above. Relative stability of the ๏ก-serine conformers was assessed within the quantum theory of atoms-in-molecules (QTAIM)95 and non-covalent interaction (NCI)96 frameworks. The MultiWFN97 suite of programs was employed for the extraction of wavefunctions from the gaussian formatted checkpoint files, and the detection of QTAIM critical points. Reduced density 43 gradient (RDG) and sign(๐2)๐ cubic files required for the NCI analysis were generated with cubegen, a gaussian 16 utility program. Chimera98 and ChimeraX99,100 were utilized for molecular visualization, including the QTAIM critical points, NCI molecular overlays, and the standalone molecules. Scatter plots of the RDG and sign(๐2)๐ cubic functions were created with Igor Pro 8, a suite of programs for graphic and data analysis. A brief overview of NCI theory can be found in Section 3.2.1 of this thesis. The vibrational potential energy distribution (PED) was calculated for each normal mode through the use of the VEDA 4 program.101,102 In short, the PED represents the normal mode coordinate as a superposition of local mode coordinates, enabling one to determine the precise energetic contributions to the vibration. Thereby allowing for one to determine which vibrational frequencies will be affected by the presence of noncovalent interactions. This analysis is similar to simply looking at the internal coordinates for a given normal mode. However, using only the internal coordinates often leads to misinterpretation of the normal mode. That is, light atoms, such as hydrogen, may show significant movement for a given normal mode while contributing very little to the change in energy observed during the vibration. It is worth mentioning that, while not presented in this thesis, geometry optimizations and anharmonic frequency calculations for the conformers of ๐ผ-serine at the CCSD(T)/cc-pVQZ level of theory are currently in progress. Briefly, these calculations are being performed with the CFOUR 2.1 suite of programs,103,104 and utilize analytical gradients, and second derivatives as implemented within CFOUR.105,106,107,108 At the time of writing, these calculations have not yet finished, and only slight modifications to the experimental assignment are expected upon their completion. 44 Ser 1 Ser-2 Ser-3 Ser-4 Ser-5 Ser-6 Ser-7 Ser-8 Ser-9 Ser-10 Ser-11 Ser-12 Ser-13 Ser-14 Ser-15 Figure 3.1 Three-dimensional molecular structures of the 15 lowest energy conformers of ๐ผ-serine. 3.2.1. Noncovalent Interactions in ๐ถ-serine In this section the role of noncovalent interactions in the conformational distribution of ๐ผ-serine will be discussed in detail. The relative ZPE corrected, and Gibbโs free energies at 420 K were 45 calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. All calculations were carried out as discussed in Section 3.2. Anharmonic corrections to the vibrational partition function were calculated through the use of SPT, with the three lowest vibrational modes being treated under the hindered rotor approximation. Previous studies have shown that the relative energy of conformers 4, 11, 13, and 14 show significant reductions when calculated with MP2 compared to those calculated with HF.18 This indicates that accounting for an adequate amount of electron correlation is extremely important, and as such, single point energies were calculated with CCSD(T)-F12, which, as discussed in the previous section, produces electronic energies near the CBS limit. The results of these calculations are summarized below in Table 3.1, whereas the relative Gibbโs free energies at other sublimation temperatures can be found in Table 3.4. Table 3.1. Tabulation of the relative ZPE-corrected electronic energies, and relative Gibbs free energies (420 K) for the 15 lowest energy ๐ผ-serine conformers. a 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 โEDFT 0.00 -0.90 1.84 1.35 7.07 5.81 4.75 7.06 6.39 7.27 7.15 10.43 11.16 15.04 6.91 โEMP22 0.00 -0.96 1.81 1.48 7.01 6.01 4.78 6.96 6.45 7.32 7.31 10.38 11.37 15.15 7.01 BE,MP2 20.4 26.8 12.1 13.3 2.7 3.6 5.2 2.8 3.2 2.5 2.5 1.0 0.8 0.3 2.7 โGDFT 0.00 1.09 5.52 2.88 5.65 7.45 6.00 7.53 8.54 7.30 4.88 8.94 9.75 14.24 5.99 โGMP2 0.00 1.46 5.34 3.24 5.77 7.93 6.09 7.47 8.77 6.70 6.59 8.59 10.27 13.98 6.33 BG,MP2 28.1 18.5 6.1 11.1 5.4 2.9 4.9 3.3 2.3 4.1 4.3 2.4 1.5 0.5 4.6 a Anharmonic thermochemical corrections were calculated as described in the main text with MP2 and B3LYP-D3BJ as denoted by the subscript MP2 and DFT. Electronic energies were calculated at the CCSD(T)-F12/cc-pVTZ-F12 level of theory. All energies are reported in kJ mol-1. The Boltzmann percentages are reported as BE,MP2 and BG,MP2 for the electronic energies and Gibbโs free energies, respectively, all calculated with MP2 values. As seen in Table 3.1, there is a large difference between the relative electronic energies and Gibbโs free energies at 420 K, indicating that thermal effects play a substantial role in the conformational distribution. In order to fully understand the thermal effects, we must first understand the stabilizing and destabilizing forces present within the molecule. The geometrical structure of each conformer is largely determined by the presence and absence of NCIs. Namely, these interactions are hydrogen bonding (HB), steric repulsion, and hyperconjugation. The previously mentioned importance of electron correlation for the accurate calculation of electronic energies confirms that noncovalent interactions play a large role in ๐ผ-serine, since significant dependence on electron correlation for the electronic energy has been linked to the presence of NCIs.109 To gain an understanding into the nature of said noncovalent interactions, they will be 46 probed through the QTAIM95 and NCI96 frameworks. Additionally, it is important to mention that these interactions have an effect on vibrational frequencies as the presence of stabilizing or destabilizing interactions results in a blue or red shift of the observed vibrational frequency, respectively, where the magnitude of the shift is dependent on the strength of the interaction. Hence, understanding these interactions will greatly simplify the spectral assignment. A brief description is included below for the analysis of noncovalent interactions through the NCI framework, as presented by Johnson et al. in 2010.96 It is well known that all chemical properties can be calculated from the quantum-mechanical electron density, denoted as ๐. Additionally, the reduced density gradient (RDG) is determined from the gradient and first derivative of the electron density, RDG =|โ๐|2(3๐2)1 3โ ๐4 3โ. (3.1) RDG is a dimensionless quantity that describes deviations from a homogenous electron distribution. Hence, in the presence of a noncovalent interaction the RDG will approach zero, becoming equal to zero at the exact point of the interaction. In 2010, Johnson et al. discovered that the presence of noncovalent interactions is characterized by regions with low density and low reduced gradient.96 Unfortunately, these regions characterize all noncovalent interactions, regardless of whether they are stabilizing or destabilizing. However, the sign of the Laplacian of the density, โ2๐, is commonly used to distinguish between different types of weak interactions.96,110 The Laplacian is commonly expressed a sum of individual contributions along the three principal axes of maximal variation, โ2๐ = ๐1 + ๐2 + ๐3. (3.2) Local maxima regions of the electron density are characterized by negative values for all three eigenvalues. It has been shown that, for noncovalent interactions, the sign of the second eigenvalue, ๐2, will be negative. However, in the case of steric crowding, the sign of ๐2will be positive. Hence, the sign of the second eigenvalue of the Laplacian can be used to distinguish between regions of stabilizing and destabilizing intramolecular interactions. Additionally, the magnitude of the sign of ๐2 can be used to characterize the strength of the interaction.96,110 This procedure was carried out to analyze the noncovalent interactions for 15 low energy conformers of ๐ผ-serine, the results of which are presented in Figure 3.2. 47 As seen in Figure 3.2, each conformer of ๐ผ-serine has at least two intramolecular interactions, with most conformers having three or four. NCI indicates that most of these interactions are weak and stabilizing, which, due to the small size of ๐ผ-serine, is surprising. The most stable conformer at 420 K, conformer 1, appears to form the weakest intramolecular interactions. On the other hand, conformer 2, which is the second most stable at 420 K and the most stable when only considering the ZPE corrected electronic energies, appears to form the strongest intramolecular interactions. Hence, there is a clear temperature dependence for the steric hinderance observed in the conformers of ๐ผ-serine. In other words, in the absence of any thermal contributions, the additional stabilization provided by the noncovalent interactions in conformer 2 make it the lowest energy conformer. However, upon increasing the temperature, the steric effects quickly overcome the stabilization provided by the additional noncovalent interactions, resulting in conformer 1 becoming the lowest energy conformer. However, these results could simply be due to the approximations made within the NCI analysis. As such, it is only fitting that the QTAIM and NCI results be explored together. 48 Figure 3.2 NCI analysis for the conformers of ๐ผ-serine. The colour bar on the right denotes whether the interaction is attractive, repulsive, or weak, denoted by red, blue, green respectively. Briefly, within the QTAIM framework molecules are divided into atomic basins that are partitioned with respect to critical points in the electron density. A critical point is defined as the Ser-1 Ser-2 Ser-3 Ser-4 Ser-5 Ser-6 Ser-7 Ser-8 Ser-9 Ser-10 Ser-11 Ser-12 Ser-13 Ser-14 Ser-15 49 point in which the gradient norm of electron density is zero. There are four types of critical points, which are classified according to the sign of the eigenvalues in the electron density Hessian matrix. These are termed, Nuclear, Bond, Ring, and Cage critical points. As expected, nuclear critical points occur at the center of nuclei, and as described above, all three eigenvalues are negative (local maxima). Bond critical points occur in between attractive atom pairs and correspond to regions with two negative eigenvalues. Ring and cage critical points occur in regions with one and no negative eigenvalues respectively. Notice that bond critical points donโt distinguish between noncovalent and covalent interactions. Additionally, QTAIM typically only observes strong intramolecular interactions. Due to the shortcomings of QTAIM, it is advantageous to use it in conjunction with NCI. Of particular interest to this thesis are the bond critical points (BCPs) and ring critical points (RCPs). From the QTAIM analysis, Figure 3.3, it is clear that most of the interactions observed with NCI are weak as most conformers display only one or two critical points, not include the nuclear or bond critical points corresponding to covalent bonds. Conformers 2, 3, 6, 7, 8, 9, 10, and 12 all contain at least one ring critical point. Based on the NCI analysis, conformers 10 and 12 form the weakest intramolecular interactions of the conformers containing a ring critical point. This is important as the relative energies of conformers 2, 3, 6, 7, 8, and 9 increase when accounting for thermal contributions, whereas those from 10 and 12 slightly decrease. Further supporting the finding that the relative stability of ๐ผ-serine conformers is largely dependent on temperature. This is the basis for varying the sublimation temperature, as a slight variation in temperature yields significant changes to the conformation distribution, with a large reduction predicted to be observed for conformers exhibiting substantial intramolecular hydrogen bonding. In addition to the conformers mentioned above, conformer 4 also increases in relative energy upon an increase in temperature. Even though conformer 4 lacks the presence of a ring critical point, NCI shows the presence of three noncovalent interactions with relatively large isosurfaces. As such, the increase in relative energy upon increasing temperature can be attributed to steric hinderances results from the size of the noncovalent interactions present in conformer 4. 50 Figure 3.3 Projection of QTAIM critical points onto ๐ผ-serine. Nuclear critical points are purple, ring critical points are blue, and bond critical points are yellow. In total, there are sixteen distinct intramolecular interactions that can occur within ๐ผ-serine. As such any discussion about the noncovalent interactions will become complicated rather quickly. Ser-1 Ser-2 Ser-3 Ser-4 Ser-5 Ser-6 Ser-7 Ser-8 Ser-9 Ser-10 Ser-11 Ser-12 Ser-13 Ser-14 Ser-15 51 Prior to any additional discussion about the QTAIM and NCI results, it is useful to partition the conformers into groups based on the noncovalent interactions they display. The groupings are reported below in Table 3.2, in which they are partitioned based upon the โdonorโ and โacceptorโ functional groups for each intramolecular interaction. This will be revisited during the spectral assignment; however, it is worthwhile to introduce it now. The table simply serves to demonstrate the number of noncovalent interactions present in each conformer. Table 3.2 Summary of noncovalent interactions present in each of the 15 lowest energy conformers. Label Conformersa Donor Group Acceptor Group Functional Atomb Functional Atomb A 1, 4, 5 Hydroxyl H2C โ O๐ Amine XC โ ๐H2 B 2, 3, 7, 10 Hydroxyl H2C โ O๐ Carboxyl (HO)C = ๐ C 11, 12 Hydroxyl H2C โ O๐ Carboxyl (H๐)C = O D 6, 8, 9, 13, 14, 15 Hydroxyl H2C โ O๐ No Interaction E 2, 6, 7, 8, 9, 11, 13, 14, 15 Amine XC โ N๐2 Hydroxyl H2C โ ๐H F 5, 7, 10, 13 Amine XC โ N๐2 Carboxyl (H๐)C = O G 1, 4, 11, 12, 14, 15 Amine XC โ N๐2 Carboxyl (HO)C = ๐ - Amine XC โ N๐2 No Interaction - Carboxyl (๐O)C = O Hydroxyl H2C โ ๐H - Carboxyl (๐O)C = O Carboxyl (HO)C = ๐ H 2, 3, 6, 8, 9 Carboxyl (๐O)C = O Amine XC โ ๐H2 I 1 Carboxyl (๐O)C = O No Interaction J 11 Methyl XC โ ๐ Hydroxyl H2C โ ๐H K 1, 4, 14, 15 Methyl XC โ ๐ Carboxyl (H๐)C = O L 2, 8, 13 Methyl XC โ ๐ Carboxyl (HO)C = ๐ - Methyl XC โ ๐ Amine XC โ ๐H2 a Empty rows in the conformer column represent that none of the conformers of interest undergo this type of intramolecular interaction. However, the rows were left in for the sake of completion. b The bolded atom in the Atom columns represents the primary atom involved in the noncovalent interaction, may it be as the donor or the acceptor. For the carboxyl functional group, (HO) denotes the hydroxyl component, while C =O denotes the carbonyl component. 3.2.2. ๐ถ-serine in Solid Parahydrogen Matrices In order to reproducibly examine the vibrational frequencies in terms of the noncovalent interactions, the total energy distribution for each normal mode was calculated. The internal coordinates, ๐๐, are defined in Table B.2 of Appendix B. To make an accurate comparison of the results with those of previous studies, Table B.3 in Appendix B converts all of the internal 52 coordinates used for the analysis of this thesis to those presented by He and Allen.18 Additionally, the theoretical vibrational, intensities, and total energy distributions (TED) for all conformers can be found in Table B.4 to Table B.11 of Appendix B. Note that the TED for the vibrational modes calculated at the DFT and MP2 levels of theory are well within the margin of error of the analysis, and as such, only the TED calculated from the MP2 frequencies are presented. A comparison of the vibrational energy distribution for each normal mode across all conformers enables one to ascertain how the vibrational modes change character, thereby providing further insight into the relative stabilities. In addition to the total energy distribution, the internal modes were utilized to assign normal modes. Based on the relative Gibbโs free energies at 420 K, Table 3.1, it is expected that conformers 1, 2, and 4 will be the most abundant species in the spectrum. This is due to the fact they are 2 kJ molโ1 lower in energy than the next lowest energy conformer, which, according to the Boltzmann distribution, indicates that they should account for 28.1, 18.5, and 11.1 % of the distribution, respectively. Additionally, conformers 3, 5, 6, 7, and 10 have been experimentally observed in argon matrices, and as such, we expect to observe their stronger features. Notice that conformer 6 has a Boltzmann percentage of 2.9 %, whereas that for conformers 8 and 15 are 3.3 and 4.6 % respectively. This is the first experiment to consider conformer 15 and based on the Boltzmann percentages it is reasonable to say that it should be present in measurable quantities. Conformer 8 has been identified in one previous study,20,21 but more recent studies were unable to observe it.23 However, based on the Boltzmann distribution, we predict that it should be present in observable amounts. Experimental observation of conformers 9, 11, 12, 13, and 14 is not expected, but they are included in the analysis for the sake of completion. That being said, conformer 11 has the exact same Boltzmann percentage as conformer 10, however, As a result of the conformational flexibility of ๐ผ-serine, the experimental spectrum has many overlapping bands and is very crowded. However, the calculated frequencies indicates that vibrational modes in the 700 โ 1000 cmโ1 region should be sufficiently spaced out enabling an initial identification of the conformers present. There are a few different vibrational modes present in this region, namely, NH2 wagging, ๐(NH2), C โ O stretching, ๐(C โ O)OH and ๐(CO)COOH, C โ N stretching, ๐(CN), and hydroxyl torsions, ๐(OH)OH and ๐(OH)COOH. Based on the 53 theoretical calculations, ๐(NH2) and ๐(OH) modes are the dominant modes in the spectrum located below 1000 cmโ1. Note that the ๐(C โ O) modes do appear above 1000 cmโ1 as well. Unfortunately, the largest peak in the spectrum below 1000 cmโ1 corresponds to ammonia, NH3, and can be seen at 968.1 cmโ1. NH3 is a well-known thermal product of ๐ผ amino acids, which is clearly evident with the peaks growing intensity with an increase in sublimation temperature. Additional NH3 bands located in the experimental spectrum are 976.9, 986.7, 1608.3, 1608.4, 1641.6, 1641.8, and 1655.67 cmโ1, all of which have been reported in previous studies.111,112 A complete list of known NH3 bands can be found in the corresponding references, and additional impurities present in the experimental spectra are discussed at the end of this subsection. As seen in Figure 3.4F, these bands are relatively sharp, and as such, have little effect on the spectral assignment. Figure 3.4 Zoomed in sections of the 700 โ 1000 cm-1 regions of the experimental ๐ผ-serine spectrum in solid parahydrogen. The TSub = 430 K spectrum has been artificially situated above the other spectra in C, D, and E such that the larger amounts of noise do not interfere with observing other spectral features. Additionally, the TSub = 430 K spectrum was not included in A and B due to an overwhelming amount of noise. The ๐(C โ O)OH and ๐(C โ O)COOH stretching modes are located with the 950 โ1250 cmโ1 region of the experimental spectrum. Compared to other normal modes, these modes 54 are spread out over a relatively large range which is due to the intramolecular interactions present in the conformers. In general, the ๐(C โ O)OH mode is located in the 950 โ 1075 cmโ1 region, and the bands are more intense than those observed for the ๐(C โ O)COOH modes, which are located in the 1150 โ 1250 cmโ1 region. This region of the spectrum also contains bands corresponding to ๐(C4 โ N) stretching modes, along with ๐ฟ(OH)OH and ๐ฟ(OH)COOH bending modes. All of these modes, with the exception of ๐(C4 โN), are dependent on the noncovalent interactions present in the conformer. Bands assigned to the ๐(C โ O)OH mode were assigned to conformers 1, 2, 3, 4, 5, 6, 8, and 15. Conformers 3 and 8 were assigned to the same band, appearing at 1044.2 cmโ1, which was done on the basis of intensity, and the broadness of the band. That is, this band has a FWHM of 16.3 cmโ1, which is substantially broader than the approximately 2 cmโ1 observed for all other conformers. Additionally, bands associated with conformers 7 and 10 were not detected, which is due to the low intensity of these transitions. The strongest band in this region is located at 1108.7 cmโ1 and has been assigned to conformer 1 on the basis of its intensity, and its agreement with the calculated values. A weaker band on the shoulder of this peak is located at 1112.2 cmโ1, and has been assigned to conformer 4 on the basis of its predicted wavenumber and intensity. That is, it is predicted that this transition will be much weaker for conformer 4 than conformer 1, and this transition is, not only on the shoulder of the transition for conformer 1, but also about one third the intensity, which is within reasonable agreement with the calculated values. As mentioned above, the ๐(C โ O)COOH modes are located in the 1150 โ 1250 cmโ1 region of the experimental spectrum. Due to the lower intensity of this transition, only conformers 1, 4, and 5 were detected. Interestingly, the intensity of this transition is strongly dependent on the noncovalent interactions present in the conformer, indicating that the C โ O bond is weaker in these conformers. A quick inspection of the experimental spectrum, Figure 3.5, one will see that the strongest band is located at 1772.2 cmโ1, which is characteristic of the carbonyl stretching mode, ๐C=O. Comparison with the theoretical spectra indeed confirms this. A closer look shows five distinct bands located at 1755.5, 1772.7, 1777.9, 1782.3, and 1790.3 cmโ1. The band at 1755.5 cmโ1 has been assigned to both conformers 7 and 10. Based on the theoretical frequencies and Gibbโs free energies, along with the experimental intensity of the band, it is clear that this band belongs to conformers 7 and 10. That is, no other conformers absorb around 1755.5 cmโ1, and the intensity of the band is too high to belong to either one of them, however, it matches up with their combined 55 intensities. Whereas the band at 1772.7 cmโ1 has been assigned to conformers 1 and 5. Based on the shape of the band, it is clear that there are overlapping bands, however, due to the large number of bands in the region skewing the baseline, this particular overlap could not be fit. As such, the 1772.7 cmโ1 band has been assigned to conformer 1, with it being noted that the conformer 5 band is overlapping but was not resolved. Likewise, the band at 1777.9cmโ1 has been assigned to conformers 3 and 4, with conformer 4 being the larger contributor. Unlike the other bands, the bands at 1782.3 and 1790.3 cmโ1 were assigned to conformers 6 and 2. These bands are much sharper than the others mentioned, and do not contain any overlap, hence, their intensities are much lower. These assignments not only agree with previous experiments,23 but also with the predicted frequencies. Figure 3.5 Zoomed in sections of the amine (N-H), A, and carbonyl (C=O), B, stretching regions for ๐ผ-serine in solid parahydrogen. The TSub = 430 K spectrum has been artificially situated above the other spectra in A such that the larger amounts of noise do not interfere with observing other spectral features. The NH2 symmetric stretch, ๐(NH2)๐ ๐ฆ๐, and asymmetric stretch, ๐(NH2)๐๐ ๐ฆ๐, modes of vibration can be seen above in Figure 3.5. It is worth noting that this region is situated in-between two sections of hydroxyl stretching bands, however, bands corresponding to ๐(NH2) modes are easily distinguished from those corresponding to ๐(OH), as are typically finer (i.e., smaller FWHM) and are much lower in intensity. Due to the decreased intensity of these bands, it is 56 expected that they will only be observed for the conformers present in fairly large amounts. In general, ๐(NH2)๐ ๐ฆ๐ bands are slightly red shifted compared to ๐(NH2)๐๐ ๐ฆ๐ bands. As a result of the low intensity and increased noise present in this region, only bands for conformers 1, 2, 3, 4, and 6 were identified. The methyl stretching, ๐(C4H), ๐(CH2)๐ ๐ฆ๐, and ๐(CH2)๐๐ ๐ฆ๐, region of the experimental spectrum spans the 2800 โ 3000 cmโ1 range, as seen in Figure 3.6. This region is particularly difficult to assign for highly flexible hydrocarbons as a result of the large number of overlapping weak transitions. For instance, based on the predicted frequencies, there should be 37 vibrational bands, most of which are strong enough to be experimentally observed. Due to the high density of transitions, it is highly unlikely than any single band arises from one conformer, but is rather comes from multiple, making any definitive assignment in this region tremendously difficult. As a result of these difficulties, only one other group attempted to assign this portion of the experimental spectrum,21 and their results are within reasonable agreement to those shown here. Bands from every conformer were assigned in this region, however, due to the large amount of overlap, these assignments should be viewed as tentative. Figure 3.6 Zoomed in sections of the methyl stretching regions for ๐ผ-serine in solid parahydrogen. The TSub = 430 K spectrum has been artificially situated above the other spectra in A such that the larger amounts of noise do not interfere with observing other spectral features. 57 ๐ผ-serine has two distinct O-H stretching vibrational modes from the hydroxyl and carboxyl functional groups. Both the hydroxyl stretching, ๐(OH)OH, and carboxyl stretching, ๐(OH)COOH, bands have been identified for conformers 1, 2, 3, 4, 5, 6, 7, 8, 10, and 15. That is, they were identified for all present conformers, and can be seen below in Figure 3.7. Interestingly, the ๐(OH)OH band corresponds to the first normal mode for every conformer, with the exception of conformers 1, 4, and 5, for which, the ๐(OH)COOH band is the first normal mode. Unfortunately, the ๐(OH)OH band for conformer 1, based on the predicted frequencies and previous experiments in argon, occurs at approximately 3542.0 cmโ1, which is the exact located of waters ๐(OH) stretching band in parahydrogen. Additionally, the water bands located at 1631.6 cmโ1 are too intense to be accurately measured, hence, assignment of the conformer 1 band cannot be done indirectly based on whether or not the band at 3542.0 cmโ1 is more intense than one would expect. The ๐(OH)OH band is located just below that for conformer 1, and as such, there is a chance that the intensity has been artificially increased as a result. In addition, the ๐(OH)OH bands for conformers 3 and 5 were previously identified at 3516.9 and 3518.6 cmโ1 respectively, in argon. Here, the presence of a single broad band was located at 3512.6 cmโ1, and was tentatively assigned to both conformers 3 and 5 based on integration of the band, A๐H2. The conformers can be partitioned into two distinct groups based on the ๐(OH)COOH mode. Firstly, the ๐(OH)COOH mode is the first normal mode for conformers 1, 4, and 5, which is located around 3550 cmโ1. All of these conformers exhibit a free carboxyl hydrogen, with the oxygen lone pair being a hydrogen bond acceptor from either an amine or alkyl hydrogen. Secondly, the ๐(OH)COOH mode is the fourth normal mode for conformers 2, 3, 6, and 8, which is located in the 3200 cmโ1 region. The carboxyl hydrogen in these conformers sufficiently interacting with the amine lone pair such that not only a critical point is found between the two, but also the presence of a ring critical point is detected. Thereby stabilizing the stretching motion and explaining the observed red shift. Lastly, the ๐(OH)COOH mode is the second normal mode for conformers 7, 10, and 15, also being located in the 3550 cmโ1 region. Unlike the other two groups, the carboxyl hydroxyl group here have no intramolecular interactions, and can be viewed as being โfreeโ. Since the first and third groups show very similar stretching frequencies, it is obvious that the intramolecular interactions present in the conformers of the first group are fairly weak. A quick 58 inspection of the QTAIM results reveals that this is in fact, true, as no critical points were identified for said interactions. Figure 3.7 Zoomed in sections of the hydroxyl stretching bands for ๐ผ-serine in solid parahydrogen. The legend shown in Figure A (top right) is the same for all of the figures shown within. The TSub = 430 K spectrum has been artificially situated above the other spectra in B through E such that the larger amounts of noise do not interfere with observing other spectral features. Table 3.3 Experimental assignment of the mid-IR vibrational bands for ๐ผ-serine in solid parahydrogen. Assignment Matrix Isolation GVPT2 Solid Parahydrogenb Argonc MP2 Conf. Modea ๐๐๐๐ ๐๐๐๐ ๐ก๐๐๐ ๐
๐๐๐ ๐๐๐ซ ๐ฐ๐๐ซ ๐๐๐๐ ๐ฐ๐๐๐ 8 ๐(OH)OH 3667.7 0.001 0.002 0.32 - 3693.2 45 15 ๐(OH)OH 3666.8 0.006 0.006 0.71 - 3696.9 46 6 ๐(OH)OH 3666.1 0.011 0.005 1.58 3658.1 s 3689.4 50 10 ๐(OH)OH 3638.1 0.021 0.004 3.21 - 3634.9 58 2 ๐(OH)OH 3630.1 0.058 0.008 4.38 3626.2 s 3647.8 46 7 ๐(OH)OH 3614.7 0.044 0.004 6.38 3616.3 60 59 Assignment Matrix Isolation GVPT2 Solid Parahydrogenb Argonc MP2 Conf. Modea ๐๐๐๐ ๐๐๐๐ ๐ก๐๐๐ ๐
๐๐๐ ๐๐๐ซ ๐ฐ๐๐ซ ๐๐๐๐ ๐ฐ๐๐๐ 4 ๐(OH)COOH 3574.5 0.310 0.040 5.00 3566.2 vs 3587.0 80 15 ๐(OH)COOH 3573.4 0.005 0.007 0.783 - 3583.7 65 10 ๐(OH)COOH 3573.2 0.096 0.016 3.57 - 3577.5 76 5 ๐(OH)COOH 3564.0 0.164 0.050 2.17 3561.0 s 3577.5 76 1 ๐(OH)COOH 3562.0 0.304 0.050 3.95 3556.8 s 3577.5 71 1 ๐(OH)OH 3542.0 3540.4 w 3544.5 66 4 ๐(OH)OH 3527.0 0.053 0.003 11.61 3538.6 w 3546.0 58 5 ๐(OH)OH 3512.6 0.013 0.002 5.13 3516.9 w 3551.2 66 3 ๐(OH)OH 3518.6 s 3528.4 112 1 ๐(NH2)sym 3440.6 0.024 0.002 7.34 - 3450.1 11 2 ๐(NH2)sym 3434.4 0.019 0.006 2.14 - 3457.3 14 3 ๐(NH2)sym 3426.9 0.001 0.002 0.42 - 3462.5 12 4 ๐(NH2)asym 3423.8 0.008 0.004 1.71 3426.4 w 3473.3 16 2 ๐(NH2)asym 3421.6 0.037 0.004 6.73 3418.6 vw 3457.3 25 6 ๐(NH2)asym 3421.6 0.030 0.004 6.73 3406.0 m 3448.0 20 1 ๐(NH2)asym 3411.5 0.031 0.008 2.44 3403.6 w 3450.1 11 3 ๐(NH2)sym 3345.5 0.005 0.001 2.57 3317.6 vs 3373.1 1 6 ๐(NH2)sym 3328.9 0.002 0.002 0.57 3317.9 vs 3377.3 8 2 ๐(NH2)sym 3364.4 7 8 ๐(OH)COOH 3263.1 0.167 0.002 47.60 - 3235.6 237 3 ๐(OH)COOH 3230.4 0.341 0.004 50.02 - 3227.0 191 6 ๐(OH)COOH - 3202.7 241 2 ๐(OH)COOH 3187.6 0.372 0.008 30.12 - 3203.6 210 1 ๐(CH2)asym 2992.7 0.013 0.003 3.32 - 3030.2 18 4 ๐(CH2)asym - 3037.8 15 2 ๐(CH2)asym 2987.3 0.027 0.005 3.79 - 3003.0 28 3 ๐(CH2)asym - 2995.3 27 1 ๐(CH2)sym 2982.2 0.032 0.003 6.70 - 3002.9 6 2975.4 0.011 0.001 5.91 - 2 ๐(C4H) 2972.4 0.028 0.004 4.57 - 2961.9 10 3 ๐(C4H) - 2977.7 10 15 ๐(C4H) 2963.1 0.056 0.003 11.19 - 2966.7 16 7 ๐(C4H) - 2976.1 17 7 ๐(CH2)sym 2952.6 0.056 0.005 7.77 - 2968.0 28 - 2944.7 0.021 0.002 6.04 - - 2938.4 0.073 0.006 7.89 - 4 ๐(CH2)sym 2928.6 0.077 0.003 15.25 - 2924.0 20 6 ๐(CH2)sym 2919.1 0.010 0.001 6.44 - 2919.0 16 10 ๐(CH2)sym - 2901.5 22 4 ๐(C4H) 2891.9 0.006 0.002 2.5 - 2899.6 38 5 ๐(CH2)sym 2888.6 0.018 0.003 4.2 - 2887.1 35 3 ๐(CH2)sym - 2877.7 10 60 Assignment Matrix Isolation GVPT2 Solid Parahydrogenb Argonc MP2 Conf. Modea ๐๐๐๐ ๐๐๐๐ ๐ก๐๐๐ ๐
๐๐๐ ๐๐๐ซ ๐ฐ๐๐ซ ๐๐๐๐ ๐ฐ๐๐๐ 1 ๐(C4H) 2882.5 0.089 0.005 11.5 - 2885.6 33 - 2857.9 0.136 0.008 11.2 - 8 ๐(C = O) 1803.9 0.023 0.003 5.78 - 1811.0 107 2 ๐(C = O) 1790.3 0.275 0.042 4.15 1789.2 vs 1801.6 66 15 ๐(C = O) - 1795.4 140 6 ๐(C = O) 1782.2 0.096 0.011 5.64 1788.9 vs 1811.4 75 4 ๐(C = O) 1777.9 0.158 0.024 4.20 1778.4 vs 1788.5 136 3 ๐(C = O) 1778.3 vs 1786.3 8 1 ๐(C = O) 1772.7 0.519 0.129 2.56 1773.0 vs 1782.7 187 5 ๐(C = O) 1771.4 vs 1789.5 156 10 ๐(C = O) 1755.5 0.263 0.081 2.23 1754.9 w 1768.9 130 7 ๐(C = O) 1767.2 125 1 ๐ฟ(COH)OH 1408.4 0.199 0.009 14.59 1410.1 m 1420.0 22 4 ๐(CH2) 1411.3 vw 1405.0 13 1 ๐ฟ(NCH) 1399.3 0.210 0.024 5.55 1400.5 w 1405.9 11 2 ๐ฟ(COH)COOH 1400.0 s 1413.2 66 5 ๐ฟ(COH)OH 1400.5 w 1415.4 21 2 ๐(CH2) 1390.3 0.089 0.007 8.11 1385.6 w 1373.5 51 4 ๐ฟ(COH)OH 1394.9 w 1349.0 2 3 ๐ฟ(COH)COOH 1376.2 0.155 0.038 2.57 1371.4 vs 1391.6 231 6 ๐ฟ(COH)COOH 1376.3 vs 1384.9 259 2 ๐ฟ(NCH) 1345.1 0.073 0.007 6.42 1344.7 m 1364.8 108 6 ๐ฟ(C4C6H) 1341.3 0.004 0.002 1.43 1339.8 w 1346.4 14 1 ๐ฟ(COH)COOH 1331.4 0.068 0.005 8.83 1328.3 w 1334.1 7 4 ๐ฟ(COH)COOH 1308.0 0.025 0.003 4.48 1305.6 w 1317.0 11 1 ๐ฟ(C4C6H) 1277.6 0.006 0.001 4.13 1278.9 w 1273.6 5 1 ๐(CH2) 1267.6 0.066 0.003 12.27 1264.8 vw 1219.2 5 5 ๐ฟ(NH2)twis 1265.1 w 1266.9 12 3 ๐(C โ O)COOH 1203.4 0.183 0.007 17.54 1202.0 w 1205.3 4 6 ๐ฟ(COH)OH 1208.9 m 1205.3 6 4 ๐ฟ(NH2)twis 1185.7 0.012 0.002 4.03 1186.1 vw 1187.2 10 5 ๐ฟ(C4C6H) 1164.9 0.069 0.011 4.15 1168.0 m 1169.3 7 1 - 1162.8 s 1163.7 9 1 ๐ฟ(NH2)twis 1151.0 0.047 0.012 2.42 1149.6 s 1154.6 62 5 ๐(C โ O)COOH 1149.6 s 1132.1 106 10 - 1139.1 0.174 0.024 4.63 1131.3 139 2 ๐ฟ(NH2)twis 1130.2 0.319 0.019 10.24 1139.7 w 1158.9 7 7 - 1136.0 w 1124.8 36 4 ๐(C โ O)COOH 1137.4 s 1132.7 30 4 ๐(C4 โ N) 1112.3 0.062 0.016 2.50 1107.8 vs 1115.3 25 1 ๐(C โ O)COOH 1108.8 0.300 0.084 2.28 1105.5 vs 1107.7 171 5 ๐(C4 โ N) 1099.8 1105.4 s 1131.8 9 61 Assignment Matrix Isolation GVPT2 Solid Parahydrogenb Argonc MP2 Conf. Modea ๐๐๐๐ ๐๐๐๐ ๐ก๐๐๐ ๐
๐๐๐ ๐๐๐ซ ๐ฐ๐๐ซ ๐๐๐๐ ๐ฐ๐๐๐ 2 ๐(C4 โ C6) 1093.2 0.047 0.015 1.95 1095.4 m 1101.8 36 3 ๐(C โ O)OH 1072.2 0.087 0.014 3.97 - 1083.0 60 5 ๐(C โ O)OH 1073.5 vs 1072.4 93 1 ๐(C โ O)OH 1066.2 0.320 0.048 4.21 1066.2 vs 1073.5 82 15 ๐(C โ O)OH 1064.6 0.019 0.009 1.38 - 1050.3 50 4 ๐(C โ O)OH 1054.7 0.039 0.018 1.42 1056.3 s 1067.2 78 6 ๐(C โ O)OH 1058.9 0.063 0.021 1.93 1056.9 s 1061.0 55 3 ๐(C โ O)OH 1044.2 0.181 0.007 16.32 1042.4 s 1083.0 60 3 ๐(C4 โ C6) 1031.7 0.004 0.001 1.91 1037.7 s 1047.4 68 6 ๐(CH2) 995.1 0.003 0.005 0.41 1002.8 w 1061.0 55 2 ๐(C โ O)OH 975.9 0.025 0.015 1.05 975.2 m 982.8 43 6 ๐(C4 โ C6) 942.9 0.025 0.006 2.51 942.7 m 950.7 49 5 ๐(NH2) 932.7 0.040 0.003 9.75 932.6 s 909.8 77 1 ๐(NH2) 928.3 0.040 0.007 3.84 927.0 w 940.6 41 2 ๐(NH2) 915.3 0.034 0.005 4.08 914.5 w 908.4 61 1 ๐(C4 โ C6) 901.3 0.052 0.006 5.57 - 893.4 19 4 - - 899.4 vw 908.2 18 3 - 874.0 0.025 0.009 1.70 871.1 w 907.1 114 6 ๐(C4 โ C6) 872.8 m 859.9 21 4 ๐(C4 โ C6) 863.1 0.037 0.005 4.82 860.6 s 871.6 101 2 ๐(OH)COOH 847.7 0.053 0.010 3.30 846.4 w 841.6 15 3 - 832.4 0.024 0.005 3.10 832.2 s 754.1 73 6 - 824.0 0.014 0.006 1.47 824.2 w 799.7 64 1 ๐(C4 โ N) 819.0 0.334 0.041 5.20 817.5 vs 838.3 106 5 ๐(NH2) 815.9 vs 814.5 106 6 ๐(C4 โ N) 799.7 64 4 ๐(NH2) 784.5 0.046 0.009 3.12 784.8 s 808.8 54 4 ๐(OH)COOH 748.4 0.045 0.024 1.21 747.5 s 751.2 34 6 - 725.7 0.026 0.011 1.48 725.3 w 735.3 19 5 ๐(OH)COOH 722.2 w 734.5 7 1 ๐(OH)COOH 718.6 0.077 0.010 4.67 721.8 w 728.5 20 a Total energy distribution for each normal mode can be found in the Appendix. b Italicized frequencies indicate that the band cannot be definitively assigned due to overlap with strong impurity bands. c Argon assignments are taken from Najbauer et al. (reference 23). It is well known that amino acids are thermally instable and readily decompose, which is evident by the presence of a number of different impurities present in the experimental spectrum. All other studies of gas-phase neutral amino acids also contain these impurities, and it is essential to find the optimal sublimation temperature, at which, one maximizes the amount of neutral amino acid, while minimizing the amount of thermal decomposition. The thermal decomposition products 62 for ๐ผ-serine present in our experimental spectrum include ammonia, NH3, carbon dioxide, CO2, and water, H2O. Vibrational frequencies for these molecules in solid parahydrogen have been thoroughly studied,113 and as such, it is easy for one to distinguish between bands corresponding to those products, and those corresponding to ๐ผ-serine. Multiple previous investigations on ๐ผ-serine in argon and nitrogen matrices have indicated the presence of ethanolamine as a thermal product.19,20,21,23 In order to ensure that vibrational bands belonging to ethanolamine were not accidentally assigned to ๐ผ-serine, the experimental spectrum of ethanolamine in solid parahydrogen was recorded and has been tentatively assigned. From this spectrum, it is clear that no ethanolamine is present in the ๐ผ-serine spectrum, and the spectrum can be found in Appendix B. Hence, we are confident that all spectral bands assigned to ๐ผ-serine do, in fact, come from the conformers of ๐ผ-serine. Unfortunately, with experimental spectra as complex as the ones presented in this section, one can never be absolutely certain of every band assigned. That being said, as discussed above, all possible measures were taken to ensure that the assignment is as accurate as possible. It should be noted that the Brendan Moore generously recorded and provided the experimental spectrum of ethanolamine. In total ten conformers of ๐ผ-serine have been identified in the solid parahydrogen matrix, including the first experimental detection of conformer 15. Jarmelo et al. also identified conformer 9, which has not been experimentally observed in any other experiments, and based on the Boltzmann percentage, it was not expected to be observed here.20 However, based on their assignment of conformer 9, it is possible that conformer 15 was present in their experimental spectrum and accidentally assigned as conformer 9, due to the fact that conformer 15 was not discovered until 11 years after the publication of this study. Additionally, according to the relative Gibbโs free energies, conformers 1 and 2 should be in the largest abundance, however, the most intense bands have been assigned to conformers 1 and 4. This is in agreement with previous studies, and likely arises as a result of the weaker absorptions displayed by conformer 2. 3.2.3. ๐ถ-serine Conformational Dependence on Sublimation Temperature Spectral figures presented in the previous section show that ๐ผ-serine has little conformational dependence on sublimation temperature. However, as seen in Table 3.4, the relative Gibbs free energies do change a fair bit with the corresponding change in temperature. Due to the large number of overlapping peaks, it is extremely difficult to correlate the change in peak area with the 63 change in conformational distribution. That being said, the thermal decomposition of serine increases substantially in the 430 K spectrum compared to that observed in the 420 K spectrum. An investigation into this behavior is currently ongoing. Table 3.4 Tabulation of the relative Gibbs free energies at 400, 410, 420, and 430 K for the 15 lowest energy conformers of ๐ผ-serine. a Anharmonic thermochemical corrections were calculated as described in the main text with MP2. Electronic energies were calculated at the CCSD(T)-F12/cc-pVTZ-F12 level of theory. The temperature values in the right most column denote the relative Gibbโs free energy at that temperature. 3.3. Excited State Dynamics of ๐ถ-alanine 3.3.1. Computational Procedures for ๐ถ-alanine The basis of the computational methods employed here are described in detail in Section 2 of this thesis. This section details the computational methods and experimental techniques required to perform the analysis, which is discussed in the following subsection. The 13 lowest energy conformers of ๏ก-alanine, as depicted in Figure 3.8, were selected from the work of Cao et al.,114 were used for the re-assessment of the ground state conformational analysis. The nomenclature employed in this thesis is simply related to the relative electronic energies of the conformers. Geometry optimizations and anharmonic vibrational frequency calculations were performed on all 13 conformers at the DFT and MP2 levels of theory, with the Gaussian 16 suite of programs.77 In addition, geometry optimizations and harmonic frequency calculations were performed on all 13 conformers at the CCSD and CCSD(T) levels of theory, with the CFOUR 2.1 suite of programs.105,106,107,108 For the sake of consistency, all calculations were performed on the L-enantiomer of ๏ก-alanine, although, the calculations between L- and D- enantiomers would yield identical results within the non-relativistic limit. The frozen core approximation was utilized in all calculations. Interestingly, all 13 conformers are local minima on the ground state PES at the CCSD(T) level of theory, indicating that the assumption that B3LYP accurately describes the physicochemical characteristics of ๏ก-alanine is a bad assumption. a 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 400 K 0.00 1.30 5.11 3.14 5.83 7.83 6.01 7.46 8.63 6.73 6.65 8.72 10.36 14.08 6.40 410 K 0.00 1.38 5.23 3.19 5.80 7.88 6.05 7.46 8.70 6.72 6.62 8.66 10.31 14.03 6.36 420 K 0.00 1.46 5.34 3.24 5.77 7.93 6.09 7.47 8.77 6.70 6.59 8.59 10.27 13.98 6.33 430 K 0.00 1.54 5.45 3.29 5.74 7.99 6.13 7.48 8.83 6.69 6.56 8.53 10.22 13.93 6.29 64 DFT and MP2 calculations performed for ๏ก-alanine are identical to those for ๏ก-serine, which were discussed in Section 3.2. In short, the B3LYP functional with the inclusion of D3 dispersion corrections78 and BJ damping79 was employed with the def2-TZVPD basis set for all DFT calculations.80,81 Since ๏ก-alanine is chemically very similar to ๏ก-serine, it should come of no surprise that a pruned integration grid consisting of 175,974 points for hydrogen, and 250,974 points for carbon, nitrogen, and oxygen was required for the accurate calculation of thermochemical properties. Unlike ๏ก-serine, the MP2 geometry optimizations and anharmonic frequency calculations were performed with the aug-cc-pVTZ basis set,82 which, due to ๏ก-alanineโs smaller size, was computationally more affordable. All geometry optimizations at the MP2 and DFT levels of theory were converged to an RMS force value of 1 ร 10โ12, and the stationary points were verified through the use of an analytical hessian calculation. Once again, the Berny algorithm83 was applied for geometry optimizations, and DVPT2 theory was utilized for the anharmonic corrections.84,85,86,87 The ORCA 4 suite of programs89,90 was employed for the calculation of highly accurate CCSD(T)-F12 energies,88 utilizing the cc-pVTZ-F12,91 cc-pVTZ-F12/CABS,92 and cc-pVQZ/C basis sets.93 To save on computational time, the RI approximation was applied to the F12 explicitly correlated terms.94 CCSD(T)-F12 energies were calculated for geometries optimized at the DFT, MP2, and CCSD(T) levels of theory. Lastly, as previously mentioned, analytical CCSD(T) geometry optimizations and harmonic frequency calculations were performed with the cc-pVTZ basis set82 within the CFOUR 2.1 suite of programs. Vertical excitation energies for the five lowest lying singlet electronic excited states were calculated using the single-reference EOM-CCSD115,116 method, and calculations with the SAC-CI117,118 method are currently ongoing. All 13 lowest energy conformers of ๐ผ-alanine were considered, Figure 3.8. These calculations were performed using the CCSD(T) optimized ground state geometries. Additionally, further analysis with multi-reference coupled cluster (MRCC)119 and fully internally contracted multi-reference configuration interaction (MRCI)120 techniques are currently ongoing. The EOM-CCSD vertical excitation energies were calculated with the ORCA 4 suite of programs,90 whereas the SAC-CI vertical excitation energies are being calculated with the Gaussian 16 suite of programs.77 All calculations utilize the aug-cc-pVTZ basis set, which was described in the above section.82 For the sake of comparison, vertical excitation energies were also calculated with CIS(D).121,122 All CIS(D) calculations were performed with the ORCA 4 suite of programs. 65 Initial geometry optimizations for the three lowest lying singlet electronic excited states were calculated with the configuration interaction singles (CIS)123 method. The excited state geometries were verified to be stationary states through the use of analytical harmonic frequency calculations. These geometries are currently being reoptimized at the EOM-CCSD and SAC-CI level of theory with the aug-cc-pVTZ basis set. It is worth mentioning that, unlike the vertical excitation energies, the EOM-CCSD geometry optimizations are being performed with the CFOUR 2.1 suite of programs.103 Lastly, the first three singlet excited states were also characterized at the TD-DFT level of theory via an electron-hole analysis as proposed by Liu et al.,124 utilizing the hybrid exchange-correlation functional proposed by Perdew-Burke-Ernzerhof (PBE0).125 The TD-DFT calculations were performed within the ORCA 4 suite of programs, while the MultiWFN suite of programs was employed for the electron-hole analysis. All TD-DFT calculations utilized Alhrichs def2-TZVP basis set, and the grid properties were identical to those for the ground state calculations. Geometry optimizations were not performed at the TD-DFT level of theory. 66 SAla-1 SAla-2 SAla-3 SAla-4 SAla-5 SAla-6 SAla-7 SAla-8 SAla-9 SAla-10 SAla-11 SAla-12 SAla-13 Figure 3.8 Three-dimensional molecular structures of the 13 lowest energy conformers of ๐ผ-alanine. 67 3.3.2. Lowest Lying Electronic States for ๐ถ-alanine The photodissociation observed in ๐ผ-alanine is consistent with the Norrish Type 1 reactions for the photolysis of carbonyl compounds. In this reaction, the initial extension of the C-C bond adjacent to the carbonyl group, is termed as ๐ผ-cleavage. Under collision-free conditions, such as those presented in solid parahydrogen, bond rupture is the most likely outcome of photoirradiation. However, the presence of caging effects may result in the rearrangement of excited products to form stable products further complicating the spectral assignment, although this effect is minimized through the use of a solid parahydrogen matrix. Norrish type I reactions typically proceed via a spin-allowed electronic excitation to a singlet excited state, followed by an intersystem crossing (ISC) to a near-lying dissociative triplet state. However, recent investigations29 have proposed that this dissociation proceeds via the first electronically excited singlet state. An electron-hole analysis was utilized to examine the change in electron density for each of the first three singlet excited states. Within this analysis, a single electron excitation is described as an electron leaving a โholeโ going to an โelectronโ. If said excitation can be entirely described as a HOMO LUMO transition, then the โholeโ is the HOMO, and the โelectronโ is the LUMO. However, very few excitations can be completely represented by a single pair of molecular orbitals (MOs). As such, we will describe the transition as the sum of the weighted MO wavefunctions, which is given by the following expression: ๐hole(๐) =โโ(๐ค๐๐)2๐๐(๐)๐๐(๐)๐ฃ๐๐๐ก๐๐๐๐๐โ ๐๐ocalhole+โโโ๐ค๐๐๐ค๐๐๐๐(๐)๐๐(๐)๐ฃ๐๐๐ก๐๐๐๐๐โ ๐๐๐๐๐โ ๐๐๐๐๐ ๐ hole, (3.3) ๐ele(๐) = โโ(๐ค๐๐)2๐๐(๐)๐๐(๐)๐ฃ๐๐๐ก๐๐๐๐๐โ ๐๐ocalele+โโโ๐ค๐๐๐ค๐๐๐๐(๐)๐๐(๐)๐ฃ๐๐๐ก๐โ ๐๐ฃ๐๐๐ก๐๐๐๐๐โ ,๐๐๐๐๐ ๐ ele (3.4) where ๐ is the MO wavefunctions and ๐ค is the weighting coefficient. Since this analysis is based upon the electron density as opposed to the wavefunction, the โholeโ and โelectronโ do not have phase. Due to the orthonormality requirement for MOs, the sum of the square of all configuration coefficients must be one, that is, โซ๐hole(๐)๐๐ = 1, (3.5) 68 โซ๐ele(๐)๐๐ = 1. (3.6) This property allows one to check that the same amount of electron density is in the โholeโ and the โelectronโ, ensuring that an adequate number of MOs were considered to describe the excitation. From this it is possible to calculate the overlap through the use of the following function, ๐(๐) = โ๐hole(๐)๐ele(๐). (3.7) The analysis presented within this thesis utilizes TD-DFT, and as such, one must also account for de-excitations. Doing so yields the following expressions for the local terms, ๐๐๐๐๐๐hole (๐) = โโ(๐ค๐๐)2๐๐(๐)๐๐(๐)๐ฃ๐๐๐ก๐๐๐๐๐โโโ(๐ค๐โฒ๐)2๐๐(๐)๐๐(๐)๐๐๐๐๐ฃ๐๐๐ก๐, (3.8) ๐๐๐๐๐๐ele (๐) = โโ(๐ค๐๐)2๐๐(๐)๐๐(๐)๐ฃ๐๐๐ก๐๐๐๐๐โโโ(๐ค๐โฒ๐)2๐๐(๐)๐๐(๐)๐๐๐๐๐ฃ๐๐๐ก๐. (3.9) Similarly, the cross terms are given by the following expressions, ๐๐๐๐๐ ๐ hole (๐) =โโโ๐ค๐๐๐ค๐๐๐๐(๐)๐๐(๐)๐ฃ๐๐๐ก๐๐๐๐๐โ ๐๐๐๐๐โโโโ๐ค๐โฒ๐๐ค๐โฒ๐๐๐(๐)๐๐(๐)๐๐๐๐โ ๐๐๐๐๐๐ฃ๐๐๐ก๐, (3.10) ๐๐๐๐๐ ๐ ele (๐) = โโโ๐ค๐๐๐ค๐๐๐๐(๐)๐๐(๐)๐ฃ๐๐๐ก๐โ ๐๐ฃ๐๐๐ก๐๐๐๐๐โโโโ๐ค๐โฒ๐๐ค๐โฒ๐๐๐(๐)๐๐(๐)๐๐๐๐๐ฃ๐๐๐ก๐โ ๐๐ฃ๐๐๐ก๐. (3.11) This analysis was performed on the four lowest energy conformers of ๐ผ-alanine, which can be seen below in Figure 3.9. 69 SAla-1 SAla-2 SAla-3 SAla-4 S0 โ S1 S0 โ S2 S0 โ S3 Figure 3.9 Electron-hole analysis for the four lowest energy conformers of ๐ผ-alanine. Blue mesh corresponds to the โelectronโ, and red mesh corresponds to the โholeโ. Looking at the first excited state, it is obvious that electron density is primarily moved from the bond between the carboxyl and chiral carbons, to the carboxyl carbon. This particular pattern is typical for Norrish Type 1 photodissociationโs,126 in which the C-C bond breaks forming two radicals, those being the HOCO and ethylamine radicals in this case. That being said, one cannot claim that photodissociation will occur based only upon an electron-hole analysis. However, to accurately model photodissociation processes one must employ multi-reference techniques. Said investigation is currently ongoing at the MRCI and MRCC levels of theory. Interestingly, the second excited state for conformers 2 and 3 show substantial loss of electron density between this same C-C bond, indicating that these conformers may undergo photodissociation forming HOCO and ethylamine radicals via either the first or second excited states. An initial look at the fourth conformer may convince the reader that it also has two dissociative states, however, the electron density is evenly distributed between both the carboxyl and chiral carbons, with the loss of electron density coming from the middle of the bond. Since the electron density is localized on the two carbons, it is highly unlikely that this state is a dissociative state. An initial investigation of the 70 photodissociation was attempted through the use of geometry optimizations for the first three singlet excited states for all 13 low energy conformers of ๐ผ-alanine at the CIS/aug-cc-pVTZ level of theory. These optimizations indicate that the first three singlet states are bound states, in disagreement with the previous results30 obtained with TD-DFT, and further confirming the need for multi-reference techniques. That being said, photochemical processes can rarely be described by single-reference methods. The vertical excitation energies are shown below in Table 3.5. In addition to the excitation energies, the singles character is reported for the EOM-CCSD calculations. A quick comparison of the two values indicates that CIS(D) slightly underestimates the excitation energy, which is a known downside to the CIS(D) method.122 However, due to the reduced cost of CIS(D) compared to that of EOM-CCSD, the values are very good. The only experimental electronic spectrum of ๐ผ-alanine was recorded in methanol by Abu-Eittah et al. in 1988.127 This spectrum, which can be found in the reference, was recorded for the 200 to 350 nm region. It contains two peaks, one at around 200 nm, starting before the spectral recording, and the second at roughly 265 nm. Unfortunately, ๐ผ-alanine in methanol forms a substantial number of stabilizing intermolecular interactions with the methanol. As such, it is expected that the electronic excitation levels will show a large shift. However, the experimental values are in reasonable agree with the energies calculated for the S1 and S2 electronic states. Table 3.5 Vertical excitation energies for the first three singlet excited states for all 13 conformers of ๐ผ-alanine at the EOM-CCSD/aug-cc-pVTZ and CIS(D)/aug-cc-pVTZ levels of theory. EOM-CCSD CIS(D) Conf. ๐๐ (eV) ๐๐ (eV) ๐๐ (eV) ๐๐ (eV) ๐๐ (eV) ๐๐ (eV) 1 5.95 (90.42)a 6.68 (90.36) 7.38 (88.73) 5.85 6.42 7.17 2 5.75 (92.00) 6.69 (92.37) 7.44 (91.86) 5.65 6.45 7.18 3 5.82 (91.86) 6.75 (92.37) 7.45 (91.81) 5.71 6.51 7.05 4 5.71 (91.53) 6.41 (92.20) 7.19 (91.63) 5.60 6.14 6.96 5 5.68 (91.61) 6.45 (92.12) 7.20 (91.53) 5.58 6.18 6.94 6 6.03 (92.05) 6.43 (92.16) 7.18 (90.91) 5.91 6.17 6.92 7 5.97 (92.05) 6.29 (92.05) 7.27 (91.37) 5.87 6.02 6.98 8 5.83 (92.09) 6.39 (92.19) 7.33 (91.59) 5.74 6.12 7.16 9 5.89 (92.08) 6.42 (92.25) 7.32 (91.66) 5.78 6.16 7.01 10 5.85 (92.08) 6.65 (92.22) 7.14 (91.72) 5.76 6.41 6.74 71 EOM-CCSD CIS(D) Conf. ๐๐ (eV) ๐๐ (eV) ๐๐ (eV) ๐๐ (eV) ๐๐ (eV) ๐๐ (eV) 11 5.70 (91.91) 6.11 (92.25) 7.17 (92.35) 5.59 5.86 6.85 12 5.94 (92.06) 6.48 (92.00) 7.03 (91.34) 5.83 6.23 6.68 13 5.89 (92.04) 6.32 (91.94) 7.28 (91.45) 5.80 6.08 6.9 a Values in parentheses correspond to the percentage of singles character for each of the states. The NCI in the ground state ๐ผ-alanine conformers were investigated through the use of NCI and QTAIM. Projections of these analyses onto the molecular structure can be seen below in Figure 3.10 and Figure 3.11. It is evident that conformers 2 and 3 display the strongest intramolecular interactions. An initial analysis indicates that the excitation energies are not dependent on the NCIs. This is contrary to the electron-hole analysis, which showed a conformational dependence on dissociative states. 72 Figure 3.10 NCI analysis for the conformers of ๐ผ-alanine. The colour bar on the right denotes whether the interaction is attractive, repulsive, or weak, denoted by red, blue, green respectively. SAla-1 SAla-2 SAla-3 SAla-4 SAla-5 SAla-6 SAla-7 SAla-8 SAla-9 SAla-10 SAla-11 SAla-12 SAla-13 73 Figure 3.11 Projection of QTAIM critical points onto ๐ผ-serine. Nuclear critical points are purple, ring critical points are blue, and bond critical points are yellow. SAla-1 SAla-2 SAla-3 SAla-4 SAla-5 SAla-6 SAla-7 SAla-8 SAla-9 SAla-10 SAla-11 SAla-12 SAla-13 74 4. Conclusion The conformational dynamics of ๐ผ-serine were investigated through the use of ab initio calculations and MI-FTIR spectroscopy with parahydrogen matrices, in which ten conformers of ๐ผ-serine were identified confirming the superiority of the parahydrogen matrix for studying the conformational landscape of highly flexible species. That is, only nine conformers were identified within argon matrices, which is one less than was observed here. Interestingly, ethanolamine was not observed as a thermal product, despite being categorized as one in previous studies. The three lowest lying singlet excited states for the four lowest energy conformers of ๐ผ-alanine were investigated through the use of an electron-hole analysis, which indicated that the first state is likely a dissociative state for all conformers. Additionally, for conformers exhibiting strong hydrogen bonds between the carboxyl hydrogen and amine lone pair, the second excited state may also be a dissociative state. All discussed dissociative states were shown to yield the HOCO and ethylamine radicals. Additionally, vertical excitation energies were calculated for all 13 low energy conformers of ๐ผ-alanine, in which, the CIS(D) and EOM-CCSD values appear to be in good agreement and confirmed the need for experimental data on the excited states of neutral alanine. 75 Bibliography 1 M. Adrian-Scotto, S. Antonczak, J. H. Bredehรถft, S. V. Hoffmann and U. J. Meierhenrich, Chiroptical Properties of Amino Acids: A Density Functional Theory Study, Symmetry, 2010, 2, 935โ949. 2 I. Wagner and H. Musso, New Naturally Occurring Amino Acids, Angewandte Chemie International Edition in English, 1983, 22, 816โ828. 3 V. Ya. Yablokov, I. L. Smelโtsova, I. A. Zelyaev and S. V. Mitrofanova, Studies of the rates of thermal decomposition of glycine, alanine, and serine, Russ J Gen Chem, 2009, 79, 1704โ1706. 4 P. Ehrenfreund, M. P. Bernstein, J. P. Dworkin, S. A. Sandford and L. J. Allamandola, The Photostability of Amino Acids in Space, ApJ, 2001, 550, L95. 5 B. Barbier, A. Chabin, D. Chaput and A. Brack, Photochemical processing of amino acids in Earth orbit, Planetary and Space Science, 1998, 46, 391โ398. 6 S. G. Stepanian, A. Yu. Ivanov, D. A. Smyrnova and L. Adamowicz, UV-induced isomerization of ฮฒ-alanine isolated in argon matrices, Journal of Molecular Structure, 2012, 1025, 6โ19. 7 P. V. Johnson, R. Hodyss, V. F. Chernow, D. M. Lipscomb and J. D. Goguen, Ultraviolet photolysis of amino acids on the surface of icy Solar System bodies, Icarus, 2012, 221, 800โ805. 8 A. J. Barnes, Matrix isolation vibrational spectroscopy as a tool for studying conformational isomerism, Journal of Molecular Structure, 1984, 113, 161โ174. 9 C. Van Alsenoy, J. N. Scarsdale, H. L. Sellers and L. Schรคfer, Ab initio studies of structural features not easily amenable to experiment: The molecular structures of two low-energy forms of unionized serine, Chemical Physics Letters, 1981, 80, 124โ126. 10 C. Van Alsenoy, S. Kulp, K. Siam, V. J. Klimkowski, J. D. Ewbank and L. Schรคfer, Ab initio studies of structural features not easily amenable to experiment: Part 63. Conformational analysis and structural study of serine, Journal of Molecular Structure: THEOCHEM, 1988, 181, 169โ178. 11 P. Tarakeshwar and S. Manogaran, Conformational effects on vibrational frequencies of cysteine and serine: an ab initio study, Journal of Molecular Structure: THEOCHEM, 1994, 305, 205โ224. 12 M. Noguera, L. Rodrฤฑฬguez-Santiago, M. Sodupe and J. Bertran, Protonation of glycine, serine and cysteine. Conformations, proton affinities and intrinsic basicities, Journal of Molecular Structure: THEOCHEM, 2001, 537, 307โ318. 76 13 A. Lazreg, I. N. Taleb-Mokhtari, N. Yousfi and M. Sekkal-Rahal, Density Functional Theory Investigations on Vibrational Spectra, Molecular Structure, and Properties of the L-Serine, L-Cysteine, and L-Aspartic Acid Molecules, Journal of the Chinese Chemical Society, 2017, 64, 503โ521. 14 K.-D. Wang, P.-F. Ma and X. Shan, Conformation effects on the molecular orbitals of serine, Chinese Phys. B, 2011, 20, 033102. 15 A. Singh, V. Gupta, A. Misra and P. Tandon, Quantum-chemical approach to serine formation in the interstellar medium: A possible reaction pathway, A&A, 2014, 563, A55. 16 S. Gronert and R. A. J. OโHair, Ab Initio Studies of Amino Acid Conformations. 1. The Conformers of Alanine, Serine, and Cysteine, J. Am. Chem. Soc., 1995, 117, 2071โ2081. 17 R. Miao, C. Jin, G. Yang, J. Hong, C. Zhao and L. Zhu, Comprehensive Density Functional Theory Study on Serine and Related Ions in Gas Phase: Conformations, Gas Phase Basicities, and Acidities, J. Phys. Chem. A, 2005, 109, 2340โ2349. 18 K. He and W. D. Allen, Conformers of Gaseous Serine, J. Chem. Theory Comput., 2016, 12, 3571โ3582. 19 B. Lambie, R. Ramaekers and G. Maes, Conformational Behavior of Serine: An Experimental Matrix-Isolation FT-IR and Theoretical DFT(B3LYP)/6-31++G** Study, J. Phys. Chem. A, 2004, 108, 10426โ10433. 20 S. Jarmelo, L. Lapinski, M. J. Nowak, P. R. Carey and R. Fausto, Preferred Conformers and Photochemical (ฮป > 200 nm) Reactivity of Serine and 3,3-Dideutero-Serine In the Neutral Form, J. Phys. Chem. A, 2005, 109, 5689โ5707. 21 S. Jarmelo and R. Fausto, Entropy effects in conformational distribution and conformationally dependent UV-induced photolysis of serine monomer isolated in solid argon, Journal of Molecular Structure, 2006, 786, 175โ181. 22 S. Blanco, M. E. Sanz, J. C. Lรณpez and J. L. Alonso, Revealing the multiple structures of serine, PNAS, 2007, 104, 20183โ20188. 23 E. E. Najbauer, G. Bazsรณ, R. Apรณstolo, R. Fausto, M. Biczysko, V. Barone and G. Tarczay, Identification of Serine Conformers by Matrix-Isolation IR Spectroscopy Aided by Near-Infrared Laser-Induced Conformational Change, 2D Correlation Analysis, and Quantum Mechanical Anharmonic Computations, J. Phys. Chem. B, 2015, 119, 10496โ10510. 77 24 M. Cao, S. Q. Newton, J. Pranata and L. Schรคfer, Ab initio conformational analysis of alanine, Journal of Molecular Structure: THEOCHEM, 1995, 332, 251โ267. 25 A. G. Csรกszรกr, Conformers of Gaseous ฮฑ-Alanine, J. Phys. Chem., 1996, 100, 3541โ3551. 26 R. M. Balabin, Conformational equilibrium in alanine: Focal-point analysis and ab initio limit, Computational and Theoretical Chemistry, 2011, 965, 15โ21. 27 H. M. Jaeger, H. F. Schaefer, J. Demaison, A. G. Csรกszรกr and W. D. Allen, Lowest-Lying Conformers of Alanine: Pushing Theory to Ascertain Precise Energetics and Semiexperimental Re Structures, J. Chem. Theory Comput., 2010, 6, 3066โ3078. 28 R. Maul, M. Preuss, F. Ortmann, K. Hannewald and F. Bechstedt, Electronic Excitations of Glycine, Alanine, and Cysteine Conformers from First-Principles Calculations, J. Phys. Chem. A, 2007, 111, 4370โ4377. 29 M. Tia, B. Cunha de Miranda, S. Daly, F. Gaie-Levrel, G. A. Garcia, L. Nahon and I. Powis, VUV Photodynamics and Chiral Asymmetry in the Photoionization of Gas Phase Alanine Enantiomers, J. Phys. Chem. A, 2014, 118, 2765โ2779. 30 A. J. Fleisher, B. J. Bjork, T. Q. Bui, K. C. Cossel, M. Okumura and J. Ye, Mid-Infrared Time-Resolved Frequency Comb Spectroscopy of Transient Free Radicals, J. Phys. Chem. Lett., 2014, 5, 2241โ2246. 31 B. A. Tom, S. Bhasker, Y. Miyamoto, T. Momose and B. J. McCall, Producing and quantifying enriched para-H2, Review of Scientific Instruments, 2009, 80, 016108. 32 S. Tam and M. E. Fajardo, Ortho/para hydrogen converter for rapid deposition matrix isolation spectroscopy, Review of Scientific Instruments, 1999, 70, 1926โ1932. 33 F. Neese, M. Atanasov, G. Bistoni, D. Maganas and S. Ye, Chemistry and Quantum Mechanics in 2019: Give Us Insight and Numbers, J. Am. Chem. Soc., 2019, 141, 2814โ2824. 34 R. Hoffmann and J.-P. Malrieu, Simulation vs. Understanding: A Tension, in Quantum Chemistry and Beyond. Part A. Stage Setting, Angewandte Chemie International Edition, 2020, 59, 12590โ12610. 35 R. Hoffmann and J.-P. Malrieu, Simulation vs. Understanding: A Tension, in Quantum Chemistry and Beyond. Part B. The March of Simulation, for Better or Worse, Angewandte Chemie, 2020, 132, 13256โ13278. 78 36 A. D. Becke, A new mixing of HartreeโFock and local densityโfunctional theories, The Journal of Chemical Physics, 1993, 98, 1372โ1377. 37 A. D. Becke, Densityโfunctional thermochemistry. III. The role of exact exchange, The Journal of Chemical Physics, 1993, 98, 5648โ5652. 38 N. Chรฉron, D. Jacquemin and P. Fleurat-Lessard, A qualitative failure of B3LYP for textbook organic reactions, Phys. Chem. Chem. Phys., 2012, 14, 7170โ7175. 39 H. Kruse, L. Goerigk and S. Grimme, Why the Standard B3LYP/6-31G* Model Chemistry Should Not Be Used in DFT Calculations of Molecular Thermochemistry: Understanding and Correcting the Problem, J. Org. Chem., 2012, 77, 10824โ10834. 40 C. E. Check and T. M. Gilbert, Progressive Systematic Underestimation of Reaction Energies by the B3LYP Model as the Number of CโC Bonds Increases: Why Organic Chemists Should Use Multiple DFT Models for Calculations Involving Polycarbon Hydrocarbons, J. Org. Chem., 2005, 70, 9828โ9834. 41 L. Lu, H. Hu, H. Hou and B. Wang, An improved B3LYP method in the calculation of organic thermochemistry and reactivity, Computational and Theoretical Chemistry, 2013, 1015, 64โ71. 42 V. Fock, Nรคherungsmethode zur Lรถsung des quantenmechanischen Mehrkรถrperproblems, Z. Physik, 1930, 61, 126โ148. 43 D. R. Hartree, The Wave Mechanics of an Atom with a Non-Coulomb Central Field. Part I. Theory and Methods, Mathematical Proceedings of the Cambridge Philosophical Society, 1928, 24, 89โ110. 44 A. D. Becke, Perspective: Fifty years of density-functional theory in chemical physics, J. Chem. Phys., 2014, 140, 18A301. 45 Chr. Mรธller and M. S. Plesset, Note on an Approximation Treatment for Many-Electron Systems, Phys. Rev., 1934, 46, 618โ622. 46 F. Coester, Bound states of a many-particle system, Nuclear Physics, 1958, 7, 421โ424. 47 J. ฤรญลพek, On the Correlation Problem in Atomic and Molecular Systems. Calculation of Wavefunction Components in UrsellโType Expansion Using QuantumโField Theoretical Methods, J. Chem. Phys., 1966, 45, 4256โ4266. 48 A. Grรผneis, S. Hirata, Y. Ohnishi and S. Ten-no, Perspective: Explicitly correlated electronic structure theory for complex systems, J. Chem. Phys., 2017, 146, 080901. 79 49 S. Ehrlich, J. Moellmann and S. Grimme, Dispersion-Corrected Density Functional Theory for Aromatic Interactions in Complex Systems, Acc. Chem. Res., 2013, 46, 916โ926. 50 B. T. Sutcliffe and R. G. Woolley, On the quantum theory of molecules, J. Chem. Phys., 2012, 137, 22A544. 51 W. Liu, Essentials of relativistic quantum chemistry, J. Chem. Phys., 2020, 152, 180901. 52 M. Born and R. Oppenheimer, Zur Quantentheorie der Molekeln, Annalen der Physik, 1927, 389, 457โ484. 53 T. Jecko, On the mathematical treatment of the Born-Oppenheimer approximation, Journal of Mathematical Physics, 2014, 55, 053504. 54 D. A. McQuarrie, Quantum chemistry, University Science Books, Sausalito, Calif, 2nd ed., 2008. 55 F. Jensen, Introduction to computational chemistry, John Wiley & Sons, Chichester, England; Hoboken, NJ, 2nd ed., 2007. 56 J. C. Slater, The Theory of Complex Spectra, Phys. Rev., 1929, 34, 1293โ1322. 57 C. C. J. Roothaan, New Developments in Molecular Orbital Theory, Rev. Mod. Phys., 1951, 23, 69โ89. 58 G. G. Hall, The molecular orbital theory of chemical valency VIII. A method of calculating ionization potentials, Proc. R. Soc. Lond. A, 1951, 205, 541โ552. 59 K. Raghavachari and J. B. Anderson, Electron Correlation Effects in Molecules, J. Phys. Chem., 1996, 100, 12960โ12973. 60 W. Kutzelnigg, How many-body perturbation theory (MBPT) has changed quantum chemistry, International Journal of Quantum Chemistry, 2009, 109, 3858โ3884. 61 J. F. Stanton and R. J. Bartlett, The equation of motion coupledโcluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties, J. Chem. Phys., 1993, 98, 7029โ7039. 62 H. Nakatsuji and K. Hirao, Cluster expansion of the wavefunction. Symmetryโadaptedโcluster expansion, its variational determination, and extension of openโshell orbital theory, J. Chem. Phys., 1978, 68, 2053โ2065. 63 I. SHAVITT, The history and evolution of configuration interaction, Molecular Physics, 1998, 94, 3โ17. 80 64 R. S. Ruoff, T. D. Klots, T. Emilsson and H. S. Gutowsky, Relaxation of conformers and isomers in seeded supersonic jets of inert gases, The Journal of Chemical Physics, 1990, 93, 3142โ3150. 65 C. D. Carlson, N. A. Seifert, M. Heger, F. Xie, J. Thomas and Y. Xu, Conformational dynamics of 1-phenyl-2,2,2-trifluoroethanol by rotational spectroscopy and ab initio calculations, Journal of Molecular Spectroscopy, 2018, 351, 62โ67. 66 J. Thomas, I. Peรฑa, C. D. Carlson, Y. Yang, W. Jรคger and Y. Xu, Structural and dynamical features of the 2,2,2-trifluoroethanolโฏammonia complex, Physical Chemistry Chemical Physics, 2020, 22, 23019โ23027. 67 J. Laane, Ed., Frontiers and advances in molecular spectroscopy, Elsevier, Amsterdam, Netherlands, 2018. 68 J. W. Ochterski, Thermochemistry in gaussian [White paper], Gaussian Inc, Pittsburgh, PA, 2000. 69 D. G. Truhlar and A. D. Isaacson, Simple perturbation theory estimates of equilibrium constants from force fields, The Journal of Chemical Physics, 1991, 94, 357โ359. 70 J. Bloino, M. Biczysko and V. Barone, General Perturbative Approach for Spectroscopy, Thermodynamics, and Kinetics: Methodological Background and Benchmark Studies, J. Chem. Theory Comput., 2012, 8, 1015โ1036. 71 A. J. Barnes, Matrix isolation vibrational spectroscopy as a tool for studying conformational isomerism, Journal of Molecular Structure, 1984, 113, 161โ174. 72 M. Tsuge and L. Khriachtchev, Tunneling Isomerization of Small Carboxylic Acids and Their Complexes in Solid Matrixes: A Computational Insight, J. Phys. Chem. A, 2015, 119, 2628โ2635. 73 G. Bazsรณ, G. Magyarfalvi and G. Tarczay, Tunneling Lifetime of the ttc/VIp Conformer of Glycine in Low-Temperature Matrices, J. Phys. Chem. A, 2012, 116, 10539โ10547. 74 Y. Grenie, J. Lassegues and C. GarrigouโLagrange, Infrared Spectrum of MatrixโIsolated Glycine, J. Chem. Phys., 1970, 53, 2980โ2982. 75 A. Y. Ivanov, G. Sheina and Yu. P. Blagoi, FTIR spectroscopic study of the UV-induced rotamerization of glycine in the low temperature matrices (Kr, Ar, Ne), Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 1998, 55, 219โ228. 81 76 H. Rostkowska, L. Lapinski and M. J. Nowak, Photochemical Generation of Benzoazetinone by UV Excitation of Matrix-Isolated Precursors: Isatin or Isatoic Anhydride, J. Phys. Chem. A, 2020, 124, 4106โ4114. 77 Gaussian 16, Revision C.01, M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. V. Marenich, J. Bloino, B. G. Janesko, R. Gomperts, B. Mennucci, H. P. Hratchian, J. V. Ortiz, A. F. Izmaylov, J. L. Sonnenberg, D. Williams-Young, F. Ding, F. Lipparini, F. Egidi, J. Goings, B. Peng, A. Petrone, T. Henderson, D. Ranasinghe, V. G. Zakrzewski, J. Gao, N. Rega, G. Zheng, W. Liang, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, K. Throssell, J. A. Montgomery, Jr., J. E. Peralta, F. Ogliaro, M. J. Bearpark, J. J. Heyd, E. N. Brothers, K. N. Kudin, V. N. Staroverov, T. A. Keith, R. Kobayashi, J. Normand, K. Raghavachari, A. P. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, J. M. Millam, M. Klene, C. Adamo, R. Cammi, J. W. Ochterski, R. L. Martin, K. Morokuma, O. Farkas, J. B. Foresman, and D. J. Fox, Gaussian, Inc., Wallingford CT, 2016. 78 S. Grimme, S. Ehrlich and L. Goerigk, Effect of the damping function in dispersion corrected density functional theory, Journal of Computational Chemistry, 32, 1456โ1465. 79 A. D. Becke and E. R. Johnson, A density-functional model of the dispersion interaction, The Journal of Chemical Physics, 2005, 123, 154101. 80 F. Weigend and R. Ahlrichs, Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy, Phys. Chem. Chem. Phys., 2005, 7, 3297โ3305. 81 D. Rappoport and F. Furche, Property-optimized Gaussian basis sets for molecular response calculations, The Journal of Chemical Physics, 2010, 133, 134105. 82 T. H. Dunning, Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen, J. Chem. Phys., 1989, 90, 1007โ1023. 83 H. B. Schlegel, Optimization of equilibrium geometries and transition structures, Journal of Computational Chemistry, 3, 214โ218. 84 V. Barone, Anharmonic vibrational properties by a fully automated second-order perturbative approach, J. Chem. Phys., 2004, 122, 014108. 82 85 J. Bloino and V. Barone, A second-order perturbation theory route to vibrational averages and transition properties of molecules: General formulation and application to infrared and vibrational circular dichroism spectroscopies, J. Chem. Phys., 2012, 136, 124108. 86 J. Bloino, A VPT2 Route to Near-Infrared Spectroscopy: The Role of Mechanical and Electrical Anharmonicity, J. Phys. Chem. A, 2015, 119, 5269โ5287. 87 M. Page, C. Doubleday and J. W. McIver, Following steepest descent reaction paths. The use of higher energy derivatives with ab initio electronic structure methods, J. Chem. Phys., 1990, 93, 5634โ5642. 88 T. B. Adler, G. Knizia and H.-J. Werner, A simple and efficient CCSD(T)-F12 approximation, J. Chem. Phys., 2007, 127, 221106. 89 F. Neese, F. Wennmohs, U. Becker and C. Riplinger, The ORCA quantum chemistry program package, J. Chem. Phys., 2020, 152, 224108. 90 F. Neese, Software update: the ORCA program system, version 4.0, WIREs Computational Molecular Science, 2018, 8, e1327. 91 F. Weigend, A. Kรถhn and C. Hรคttig, Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculations, J. Chem. Phys., 2002, 116, 3175โ3183. 92 J. G. Hill and K. A. Peterson, Correlation consistent basis sets for explicitly correlated wavefunctions: Pseudopotential-based basis sets for the post-d main group elements GaโRn, J. Chem. Phys., 2014, 141, 094106. 93 S. Kritikou and J. G. Hill, Auxiliary Basis Sets for Density Fitting in Explicitly Correlated Calculations: The Atoms HโAr, J. Chem. Theory Comput., 2015, 11, 5269โ5276. 94 H.-J. Werner and F. R. Manby, Explicitly correlated second-order perturbation theory using density fitting and local approximations, J. Chem. Phys., 2006, 124, 054114. 95 R. F. W. Bader, A quantum theory of molecular structure and its applications, Chem. Rev., 1991, 91, 893โ928. 96 E. R. Johnson, S. Keinan, P. Mori-Sรกnchez, J. Contreras-Garcรญa, A. J. Cohen and W. Yang, Revealing Noncovalent Interactions, J. Am. Chem. Soc., 2010, 132, 6498โ6506. 97 T. Lu and F. Chen, Multiwfn: A multifunctional wavefunction analyzer, Journal of Computational Chemistry, 33, 580โ592. 83 98 E. F. Pettersen, T. D. Goddard, C. C. Huang, G. S. Couch, D. M. Greenblatt, E. C. Meng and T. E. Ferrin, UCSF ChimeraโA visualization system for exploratory research and analysis, Journal of Computational Chemistry, 25, 1605โ1612. 99 E. F. Pettersen, T. D. Goddard, C. C. Huang, E. C. Meng, G. S. Couch, T. I. Croll, J. H. Morris and T. E. Ferrin, UCSF ChimeraX: Structure visualization for researchers, educators, and developers, Protein Sci, 2021, 30, 70โ82. 100 T. D. Goddard, C. C. Huang, E. C. Meng, E. F. Pettersen, G. S. Couch, J. H. Morris and T. E. Ferrin, UCSF ChimeraX: Meeting modern challenges in visualization and analysis, Protein Sci, 2018, 27, 14โ25. 101 Michal. J. Jamrรณz, Vibrational Energy Distribution Analysis VEDA 4, Warsaw, 2004-2010. 102 M. H. Jamrรณz, Vibrational Energy Distribution Analysis (VEDA): Scopes and limitations, Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2013, 114, 220โ230. 103 D. A. Matthews, L. Cheng, M. E. Harding, F. Lipparini, S. Stopkowicz, T.-C. Jagau, P. G. Szalay, J. Gauss and J. F. Stanton, Coupled-cluster techniques for computational chemistry: The CFOUR program package, J. Chem. Phys., 2020, 152, 214108. 104 J. F. Stanton, J. Gauss, L. Cheng, M. E. Harding, D. A. Matthews and P. G. Szalay, CFOUR, Coupled-Cluster techniques for Computational Chemistry, a quantum-chemical program package, 2019. For the current version, see http://www.cfour.de 105 G. E. Scuseria, Analytic evaluation of energy gradients for the singles and doubles coupled cluster method including perturbative triple excitations: Theory and applications to FOOF and Cr2, J. Chem. Phys., 1991, 94, 442โ447. 106 M. E. Harding, T. Metzroth, J. Gauss and A. A. Auer, Parallel Calculation of CCSD and CCSD(T) Analytic First and Second Derivatives, J. Chem. Theory Comput., 2008, 4, 64โ74. 107 M. Kรกllay and J. Gauss, Analytic second derivatives for general coupled-cluster and configuration-interaction models, J. Chem. Phys., 2004, 120, 6841โ6848. 108 H. Koch, H. J. Aa. Jensen, P. Jo/rgensen, T. Helgaker, G. E. Scuseria and H. F. Schaefer, Coupled cluster energy derivatives. Analytic Hessian for the closedโshell coupled cluster singles and doubles wave function: Theory and applications, J. Chem. Phys., 1990, 92, 4924โ4940. 84 109 R. A. Shaw and J. G. Hill, A Linear-Scaling Method for Noncovalent Interactions: An Efficient Combination of Absolutely Localized Molecular Orbitals and a Local Random Phase Approximation Approach, J. Chem. Theory Comput., 2019, 15, 5352โ5369. 110 R. F. W. Bader and H. Essรฉn, The characterization of atomic interactions, J. Chem. Phys., 1984, 80, 1943โ1960. 111 M. Ruzi and D. T. Anderson, Matrix Isolation Spectroscopy and Nuclear Spin Conversion of NH3 and ND3 in Solid Parahydrogen, J. Phys. Chem. A, 2013, 117, 9712โ9724. 112 M. Ruzi and D. T. Anderson, Fourier Transform Infrared Studies of Ammonia Photochemistry in Solid Parahydrogen, J. Phys. Chem. A, 2013, 117, 13832โ13842. 113 W. R. Wonderly and D. T. Anderson, Transient HDO rovibrational satellite peaks in solid parahydrogen: Evidence of hydrogen atoms or vacancies, Low Temperature Physics, 2012, 38, 673โ678. 114 M. Cao, S. Q. Newton, J. Pranata and L. Schรคfer, Ab initio conformational analysis of alanine, Journal of Molecular Structure: THEOCHEM, 1995, 332, 251โ267. 115 J. F. Stanton and R. J. Bartlett, The equation of motion coupledโcluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties, J. Chem. Phys., 1993, 98, 7029โ7039. 116 M. Kรกllay and J. Gauss, Calculation of excited-state properties using general coupled-cluster and configuration-interaction models, J. Chem. Phys., 2004, 121, 9257โ9269. 117 H. Nakatsuji and K. Hirao, Cluster expansion of the wavefunction. Symmetryโadaptedโcluster expansion, its variational determination, and extension of openโshell orbital theory, J. Chem. Phys., 1978, 68, 2053โ2065. 118 R. Cammi, R. Fukuda, M. Ehara and H. Nakatsuji, Symmetry-adapted cluster and symmetry-adapted cluster-configuration interaction method in the polarizable continuum model: Theory of the solvent effect on the electronic excitation of molecules in solution, J. Chem. Phys., 2010, 133, 024104. 119 F. A. Evangelista, Perspective: Multireference coupled cluster theories of dynamical electron correlation, J. Chem. Phys., 2018, 149, 030901. 120 M. Saitow, Y. Kurashige and T. Yanai, Fully Internally Contracted Multireference Configuration Interaction Theory Using Density Matrix Renormalization Group: A Reduced- 85 Scaling Implementation Derived by Computer-Aided Tensor Factorization, J. Chem. Theory Comput., 2015, 11, 5120โ5131. 121 A perturbative correction to restricted open shell configuration interaction with single substitutions for excited states of radicals, Chemical Physics Letters, 1995, 246, 114โ121. 122 Martin Head-Gordon, R. J. Rico, M. Oumi and T. Lee, A doubles correction to electronic excited states from configuration interaction in the space of single substitutions, Chemical Physics Letters, 1994, 219, 21โ29. 123 J. B. Foresman, M. Head-Gordon, J. A. Pople and M. J. Frisch, Toward a systematic molecular orbital theory for excited states, J. Phys. Chem., 1992, 96, 135โ149. 124 Z. Liu, T. Lu and Q. Chen, An sp-hybridized all-carboatomic ring, cyclo[18]carbon: Electronic structure, electronic spectrum, and optical nonlinearity, Carbon, 2020, 165, 461โ467. 125 J. P. Perdew, M. Ernzerhof and K. Burke, Rationale for mixing exact exchange with density functional approximations, J. Chem. Phys., 1996, 105, 9982โ9985. 126 B. Marchetti, T. N. V. Karsili and M. N. R. Ashfold, Exploring Norrish type I and type II reactions: an ab initio mechanistic study highlighting singlet-state mediated chemistry, Phys. Chem. Chem. Phys., 2019, 21, 14418โ14428. 127 R. Abu-Eittah, A. Obaid, S. Basahl and E. Diefallah, Molecular Orbital Treatment of Some Amino Acids, BCSJ, 1988, 61, 2609โ2613. 86 Appendices A. Supplemental Material for the Experimental and Computational Methods Algorithm 1. Hartree-Fock Variational Self-Consistent Field Procedure 1 Input a guess of the molecular orbital coefficients. 2 Form the Fock matrix. 3 Diagonalize the Fock matrix. 4 Obtain the new density matrix. 5A If the new density matrix is not the same as the initial density matrix, repeat steps 2 through 4. 5B If the new density matrix is equal to the initial density matrix (within some value) then the procedure is completed. 6 Check the converged energies and orbital coefficients. 87 B. Supplemental Material for the Conformational Analysis of Serine The following section contains the raw data from the conformational analysis of L-serine presented in Section 3. This includes the full FTIR spectra taken at each sublimation temperature, the theoretically calculated infrared frequencies and their associated intensities, and the QTAIM and NCI figures for the higher energy conformers. Figure B.1 FTIR spectrum of L-serine immediately after deposition in the solid parahydrogen matrix (flow rate = 5 ccm, TSub = 400 K, TDep. = 4 K). 88 Figure B.2 FTIR spectrum of L-serine immediately after deposition in the solid parahydrogen matrix (flow rate = 5 ccm, TSub = 410 K, TDep. = 4 K). 89 Figure B.3 FTIR spectrum of L-serine immediately after deposition in the solid parahydrogen matrix (flow rate = 5 ccm, TSub = 420 K, TDep. = 4 K). 90 Figure B.4 FTIR spectrum of L-serine immediately after deposition in the solid parahydrogen matrix (flow rate = 5 ccm, TSub = 430 K, TDep. = 4 K). Table B.1 Nomenclature conversion table for the conformers of serine in select recent investigations. Here Allena Blancob Lambiec Jarmelod Najbauere 1 II Ia 1 1 1 2 I IIb 2 2 2 3 III IIc 3 3 3 4 IV Iโb 4 4 4 5 VIII III๏กa 5 6 5 6 VI IIa - 5 6 7 V III๏ขb 6 7 7 8 IX - 7 8 8 9 VII - - 9 9 10 XII III๏ขc 8 - 10 11 X Ib 9 - 11 12 - Ic 10 - 12 13 - - 11 - 13 14 - - 12 - 14 15 XI - - - - a Reference 18. b Reference 22. c Reference 19. d References 20 and 21. e Reference 23. 91 Table B.2 Definition of internal coordinates, Si, for use in the normal mode analysis and total energy distribution reported for each normal mode. Atom numbering corresponds to that shown in Figure 3.1. Bond Length Bond Angle Dihedral Angle ๐1 ๐(N1 โ H2) ๐14 ๐(C4 โ N1 โ H2) ๐26 ๐(H2 โN1 โ C4 โ H5) ๐2 ๐(N1 โ H3) ๐15 ๐(C4 โ N1 โ H3) ๐27 ๐(H2 โN1 โ C4 โ C6) ๐3 ๐(N1 โ C4) ๐16 ๐(N1 โ C4 โ H5) ๐28 ๐(H2 โN1 โ C4 โ C7) ๐4 ๐(C4 โ H5) ๐17 ๐(N1 โ C4 โ C6) ๐29 ๐(C7 โ C4 โN1 โ H3) ๐5 ๐(C4 โ C6) ๐18 ๐(N1 โ C4 โ C7) ๐30 ๐(N1 โ C4 โ C6 โ H8) ๐6 ๐(C4 โ C7) ๐19 ๐(C4 โ C6 โH8) ๐31 ๐(N1 โ C4 โ C6 โ H9) ๐7 ๐(C6 โ H8) ๐20 ๐(C4 โ C6 โH9) ๐32 ๐(N1 โ C4 โ C6 โ O10) ๐8 ๐(C6 โ H9) ๐21 ๐(C4 โ C6 โO10) ๐33 ๐(N1 โ C4 โ C7 โ O11) ๐9 ๐(C6 โ O10) ๐22 ๐(C4 โ C7 โO11) ๐34 ๐(C7 โ C4 โ C6 โ O10) ๐10 ๐(C7 โ O11) ๐23 ๐(C6 โ O10 โ H12) ๐35 ๐(C4 โ C7 โ O13 โH14) ๐11 ๐(C7 โ O13) ๐24 ๐(C4 โ C7 โO13) ๐36 ๐(O13 โ C7 โ C4 โO11) ๐12 ๐(O10 โH12) ๐25 ๐(C7 โ O13 โ H14) ๐13 ๐(O13 โH14) Table B.3 Conversion of internal coordinates, Si, to those presented by He and Allen (Ref. 18). Here He and Allen Here He and Allen Here He and Allen ๐1 ๐1 ๐14 ๐3 ๐26 ๐6 ๐2 ๐34 ๐15 ๐35 ๐27 ๐9 ๐3 ๐2 ๐16 ๐5 ๐28 ๐12 ๐4 ๐4 ๐17 ๐8 ๐29 ๐36 ๐5 ๐7 ๐18 ๐11 ๐30 ๐15 ๐6 ๐10 ๐19 ๐14 ๐31 ๐18 ๐7 ๐16 ๐20 ๐17 ๐32 ๐21 ๐8 ๐13 ๐21 ๐20 ๐33 ๐24 ๐9 ๐19 ๐22 ๐23 ๐34 ๐33 ๐10 ๐22 ๐23 ๐26 ๐35 ๐27 ๐11 ๐28 ๐24 ๐29 ๐36 ๐30 ๐12 ๐25 ๐25 ๐32 ๐13 ๐31 92 Table B.4 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 1 and 2, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. Mode Ser-1 Ser-2 MP2 B3LYP-D3BJ MP2 B3LYP-D3BJ ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐1 3577.5 71 3532.8 51 3647.8 46 3602.8 38 ๐2 3544.5 66 3505.9 59 3457.3 14 3408.1 13 ๐3 3450.1 11 3402.4 8 3364.4 7 3333.4 4 ๐4 3366.8 7 3341.6 4 3203.6 210 3183.1 231 ๐5 3030.2 18 2984.1 15 3017.8 11 2948.1 15 ๐6 3002.9 6 2935.6 12 3003.0 28 2954.6 22 ๐7 2885.6 33 2827.5 40 2961.9 10 2908.4 5 ๐8 1782.7 187 1773.8 231 1801.6 66 1796.8 152 ๐9 1595.7 23 1572.9 7 1586.0 24 1579.4 22 ๐10 1485.6 1 1471.3 2 1467.6 8 1454.9 4 ๐11 1420.0 22 1403.7 19 1413.2 66 1390.8 41 ๐12 1405.9 11 1388.8 5 1373.5 51 1342.9 32 ๐13 1351.5 4 1338.2 2 1364.8 108 1356.0 99 ๐15 1334.1 7 1317.7 8 1346.0 23 1345.6 15 ๐16 1273.6 5 1262.2 5 1278.2 5 1267.7 9 ๐17 1219.2 5 1219.2 12 1215.7 3 1200.7 14 ๐18 1163.7 9 1157.3 3 1197.2 13 1169.3 13 ๐19 1154.6 62 1126.9 56 1158.9 7 1150.1 5 ๐20 1107.7 171 1086.6 98 1101.8 36 1078.2 37 ๐21 1073.5 82 1050.1 68 1073.0 9 1051.6 17 ๐22 1011.4 4 996.6 2 982.8 43 964.5 74 ๐23 940.6 41 909.2 22 925.8 21 911.0 2 ๐24 893.4 19 870.6 52 761.4 48 873.6 89 ๐25 838.3 106 806.3 116 908.4 61 738.4 53 ๐26 728.5 20 722.9 4 841.6 15 825.8 3 ๐27 658.3 10 657.6 14 764.8 6 759.5 10 ๐28 578.0 75 547.8 32 594.8 4 586.3 3 ๐29 488.4 69 439.5 82 553.4 5 543.5 5 ๐30 550.5 28 581.3 69 511.0 3 509.2 3 ๐31 407.9 21 403.8 18 279.1 58 270.3 73 ๐31 345.2 2 334.9 1 420.9 13 422.8 7 ๐32 307.8 2 316.1 2 361.9 18 329.3 21 ๐33 210.8 6 154.7 8 298.5 15 288.5 13 ๐34 178.9 6 211.7 1 225.9 3 221.1 2 ๐35 133.5 10 124.2 10 108.6 8 102.8 7 ๐36 54.0 1 50.5 1 69.8 4 69.2 6 93 Table B.5 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 5 and 4, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. Mode Ser-3 Ser-4 MP2 B3LYP-D3BJ MP2 B3LYP-D3BJ ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐1 3528.4 112 3471.2 135 3587.0 80 3546.9 65 ๐2 3462.5 12 3412.7 9 3546.0 58 3513.1 48 ๐3 3373.1 1 3343.8 0 3473.3 16 3428.7 20 ๐4 3227.0 192 3203.6 212 3380.7 12 3352.8 7 ๐5 2995.3 27 2904.9 21 3037.8 15 2978.6 19 ๐6 2977.7 10 2919.9 10 2924.0 20 2866.5 38 ๐7 2877.7 27 2813.4 50 2899.6 38 2828.7 50 ๐8 1786.3 8 1779.1 98 1788.5 136 1775.2 128 ๐9 1567.2 7 1635.3 34 1569.9 3 1623.7 22 ๐10 1481.5 0 1471.3 1 1482.6 2 1466.3 1 ๐11 1436.0 14 1416.2 9 1417.8 13 1404.0 15 ๐12 1391.6 231 1344.9 60 1405.0 35 1382.4 16 ๐13 1358.1 93 1367.0 58 1349.9 2 1338.2 1 ๐15 1347.3 4 1340.9 53 1317.0 11 1298.3 9 ๐16 1281.1 4 1279.8 9 1275.8 1 1261.0 3 ๐17 1212.2 2 1200.5 2 1205.8 7 1200.0 22 ๐18 1205.3 4 1178.8 7 1187.2 10 1179.5 11 ๐19 1163.7 2 1161.0 8 1132.7 30 1122.4 101 ๐20 1110.4 4 1082.2 10 1115.3 25 1090.9 46 ๐21 1083.0 60 1058.8 106 1067.2 78 1042.1 85 ๐22 1047.4 68 1014.8 42 994.6 1 979.9 1 ๐23 985.7 12 979.0 16 908.2 18 884.7 14 ๐24 907.1 114 884.7 70 871.6 101 838.6 24 ๐25 754.1 73 727.5 34 808.8 54 772.3 59 ๐26 810.4 1 792.6 7 751.2 34 740.1 58 ๐27 733.5 5 733.1 5 621.1 57 614.1 52 ๐28 476.5 84 684.2 1 655.4 21 619.4 10 ๐29 663.7 6 433.6 88 472.3 35 485.7 22 ๐30 548.8 4 541.3 3 524.4 31 515.3 32 ๐31 418.9 3 416.7 2 478.7 42 449.8 118 ๐31 381.5 7 375.5 9 365.0 2 359.2 3 ๐32 328.8 21 325.2 13 296.1 4 251.0 6 ๐33 246.1 18 242.7 14 230.1 15 309.2 5 ๐34 229.9 14 227.5 14 191.4 21 138.5 7 ๐35 162.3 4 152.6 4 109.7 5 103.3 5 ๐36 62.7 2 63.4 2 69.5 1 68.4 1 94 Table B.6 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 5 and 4, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. Mode Ser-5 Ser-6 MP2 B3LYP-D3BJ MP2 B3LYP-D3BJ ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐1 3577.5 76 3538.7 60 3689.4 50 3642.7 42 ๐2 3551.2 66 3511.7 56 3448.0 20 3399.4 16 ๐3 3462.7 10 3414.6 8 3377.3 8 3352.2 2 ๐4 3382.1 5 3357.4 2 3202.7 241 3180.8 193 ๐5 3026.1 14 2982.2 13 3003.3 14 2924.3 24 ๐6 3006.8 5 2951.1 8 2977.2 8 2906.4 8 ๐7 2887.1 35 2834.9 36 2919.0 16 2871.9 19 ๐8 1789.5 156 1773.4 127 1811.4 75 1809.7 98 ๐9 1580.0 8 1600.4 8 1604.4 10 1597.9 23 ๐10 1486.4 1 1472.0 2 1490.7 3 1475.0 4 ๐11 1415.4 21 1398.0 26 1420.6 1 1404.2 1 ๐12 1377.6 9 1367.6 11 1384.9 259 1356.1 32 ๐13 1345.6 0 1332.1 2 1346.4 14 1335.6 19 ๐15 1298.9 32 1332.3 5 1310.5 3 1290.1 1 ๐16 1266.9 12 1260.1 25 1282.0 8 1278.7 13 ๐17 1247.9 1 1244.4 9 1227.6 5 1216.6 14 ๐18 1169.3 7 1160.7 8 1205.3 6 1177.0 4 ๐19 1132.1 106 1112.9 79 1165.0 12 1150.6 13 ๐20 1131.8 9 1098.2 7 1103.5 40 1078.5 37 ๐21 1072.4 93 1048.9 52 1061.0 55 1030.4 57 ๐22 1001.9 5 987.7 1 1011.1 14 993.5 12 ๐23 909.8 77 878.2 34 950.7 49 924.7 45 ๐24 897.2 15 868.6 26 920.9 54 859.8 71 ๐25 814.5 106 777.8 118 799.7 64 897.4 11 ๐26 735.3 19 730.0 16 859.9 21 750.6 60 ๐27 661.9 27 649.3 28 734.5 5 726.5 9 ๐28 557.2 41 548.7 37 649.1 7 646.8 8 ๐29 459.5 69 593.0 3 557.3 6 553.7 8 ๐30 587.7 19 432.1 74 422.8 9 417.5 6 ๐31 412.7 11 406.8 11 402.5 10 340.2 6 ๐31 345.9 2 337.0 2 288.3 37 269.0 41 ๐32 297.4 9 284.1 2 320.6 2 363.5 9 ๐33 235.7 5 235.7 8 194.5 72 153.4 61 ๐34 182.0 2 184.7 2 192.6 12 199.5 8 ๐35 134.4 5 127.7 4 135.9 9 127.3 10 ๐36 30.7 1 29.9 1 78.6 0 70.4 0 95 Table B.7 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 7 and 8, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. Mode Ser-7 Ser-8 MP2 B3LYP-D3BJ MP2 B3LYP-D3BJ ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐1 3616.3 60 3553.5 60 3693.2 45 3645.9 29 ๐2 3575.5 75 3528.2 57 3452.4 12 3407.9 9 ๐3 3458.2 11 3414.4 9 3369.6 1 3339.1 0 ๐4 3369.0 3 3342.7 4 3235.6 237 3191.7 223 ๐5 3029.5 16 2966.6 21 3008.1 14 2922.4 18 ๐6 2968.0 28 2911.8 28 2975.6 9 2911.6 1 ๐7 2976.1 17 2855.4 22 2909.7 24 2860.2 27 ๐8 1767.2 125 1751.9 21 1811.0 107 1802.8 140 ๐9 1573.9 9 1624.3 11 1600.6 25 1592.0 23 ๐10 1468.4 6 1455.8 3 1492.4 2 1478.2 3 ๐11 1408.2 22 1398.3 7 1416.8 2 1403.7 0 ๐12 1382.8 20 1366.5 6 1388.4 54 1349.2 245 ๐13 1370.4 14 1350.3 11 1359.1 75 1374.8 42 ๐15 1353.4 5 1348.5 12 1292.2 13 1285.6 9 ๐16 1286.0 1 1250.2 7 1257.5 3 1243.8 7 ๐17 1229.9 5 1219.4 9 1214.8 10 1189.8 8 ๐18 1174.9 0 1169.8 0 1191.1 5 1169.5 8 ๐19 1124.8 36 1109.7 234 1168.4 3 1156.5 5 ๐20 1117.4 41 1092.0 27 1108.1 20 1083.4 17 ๐21 1052.1 6 1025.2 10 1078.9 18 1050.8 22 ๐22 998.7 31 982.4 48 1041.7 86 1012.5 81 ๐23 916.4 9 892.7 2 949.9 38 933.0 39 ๐24 851.9 84 807.5 35 885.4 74 897.9 71 ๐25 826.6 50 781.7 110 905.5 6 842.2 10 ๐26 747.4 26 748.4 23 755.6 54 740.2 30 ๐27 699.1 13 693.9 16 727.4 6 725.4 6 ๐28 572.4 94 577.9 91 612.5 15 610.3 10 ๐29 646.8 47 644.0 45 525.1 0 521.3 1 ๐30 407.8 118 416.7 121 425.2 8 420.6 6 ๐31 387.8 5 376.3 6 388.2 11 366.1 2 ๐31 396.9 10 394.6 7 320.3 7 309.7 10 ๐32 274.6 9 242.1 22 344.2 44 214.2 10 ๐33 243.9 14 259.9 2 114.7 52 118.2 52 ๐34 218.8 1 213.3 1 189.2 12 321.8 34 ๐35 159.4 8 151.1 8 119.2 12 115.7 11 ๐36 40.0 1 38.9 1 63.2 2 57.1 1 96 Table B.8 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 9 and 10, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. Mode Ser-9 Ser-10 MP2 B3LYP-D3BJ MP2 B3LYP-D3BJ ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐1 3673.8 38 3621.4 31 3634.9 58 3584.0 46 ๐2 3449.2 20 3400.0 15 3577.5 76 3534.3 56 ๐3 3376.7 11 3352.0 3 3455.2 3 3410.3 0 ๐4 3217.6 251 3201.4 181 3371.1 2 3348.3 1 ๐5 3021.3 13 2955.2 18 3031.6 17 2969.6 19 ๐6 3001.1 17 2876.6 15 2983.1 18 2925.5 16 ๐7 2952.5 12 2889.3 12 2901.5 22 2854.1 31 ๐8 1807.0 106 1795.7 227 1768.9 130 1754.3 232 ๐9 1598.0 24 1629.1 6 1570.6 21 1640.6 12 ๐10 1478.5 3 1465.2 3 1476.6 2 1460.0 4 ๐11 1376.3 252 1384.0 183 1399.4 15 1395.0 15 ๐12 1397.2 12 1350.8 262 1377.4 3 1358.0 4 ๐13 1352.6 13 1344.9 0 1364.4 24 1343.6 8 ๐15 1333.8 8 1323.2 22 1355.2 3 1340.7 13 ๐16 1284.0 2 1280.4 3 1280.6 12 1265.9 13 ๐17 1248.4 1 1244.4 1 1226.7 5 1215.3 17 ๐18 1200.8 13 1172.8 23 1171.0 2 1156.8 4 ๐19 1132.3 14 1117.1 25 1131.3 139 1115.6 183 ๐20 1107.0 37 1085.7 21 1117.9 72 1082.2 19 ๐21 1057.5 42 1029.7 95 1080.1 23 1048.5 108 ๐22 990.6 15 976.4 9 1020.3 19 1003.5 22 ๐23 942.8 19 911.3 36 987.4 3 971.3 5 ๐24 902.5 26 868.7 61 838.6 126 833.8 13 ๐25 781.2 65 880.2 28 807.6 40 735.6 129 ๐26 883.0 49 741.9 60 738.0 11 733.3 10 ๐27 721.4 10 716.6 12 642.1 19 647.1 12 ๐28 653.5 4 647.5 3 544.1 52 527.0 78 ๐29 557.9 2 553.9 2 607.3 48 603.9 25 ๐30 424.5 3 423.1 23 375.4 79 338.0 105 ๐31 306.5 58 391.1 4 448.3 6 447.8 6 ๐31 359.6 21 333.4 3 289.7 6 286.4 4 ๐32 416.1 21 253.3 16 368.5 1 356.5 2 ๐33 184.5 39 221.5 100 168.8 28 176.8 25 ๐34 188.4 7 191.8 8 208.5 1 203.8 1 ๐35 133.2 12 123.5 13 159.7 2 143.8 3 ๐36 81.3 3 76.1 2 44.9 1 40.9 2 97 Table B.9 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 11 and 12, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. Mode Ser-11 Ser-12 MP2 B3LYP-D3BJ MP2 B3LYP-D3BJ ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐1 3674.3 34 3624.4 26 3653.1 48 3597.0 39 ๐2 3573.7 65 3531.2 45 3566.3 86 3522.8 62 ๐3 3437.6 10 3394.0 8 3442.5 5 3397.0 5 ๐4 3368.6 3 3344.1 2 3372.8 2 3352.3 1 ๐5 3034.5 11 2971.9 15 3032.8 11 2988.2 9 ๐6 3004.9 24 2957.9 27 2989.8 9 2922.0 13 ๐7 2940.7 10 2866.3 7 2891.6 30 2840.7 38 ๐8 1791.0 160 1780.6 237 1790.4 93 1785.9 124 ๐9 1604.3 15 1604.8 11 1603.0 15 1655.2 2 ๐10 1469.4 3 1459.8 4 1488.6 1 1472.2 0 ๐11 1388.8 12 1378.6 7 1406.2 18 1395.6 13 ๐12 1369.7 9 1357.1 5 1377.4 2 1369.9 0 ๐13 1357.1 2 1341.4 2 1353.9 8 1335.7 5 ๐15 1337.7 0 1319.2 3 1302.2 9 1283.5 6 ๐16 1278.3 1 1270.5 2 1251.3 6 1244.9 3 ๐17 1239.5 2 1224.5 4 1211.6 4 1202.2 11 ๐18 1171.4 14 1166.5 15 1188.4 1 1176.7 7 ๐19 1131.7 48 1116.5 35 1143.6 13 1118.0 23 ๐20 1105.9 203 1084.6 203 1101.4 214 1081.6 198 ๐21 1047.2 23 1018.9 30 1075.0 29 1038.6 40 ๐22 996.2 47 976.0 49 1017.8 5 1001.5 13 ๐23 908.7 60 878.5 6 985.6 2 959.8 4 ๐24 893.8 77 848.2 233 868.2 133 813.3 163 ๐25 831.0 10 816.6 14 792.2 10 771.6 5 ๐26 737.3 11 732.1 21 735.7 8 730.7 12 ๐27 632.1 82 622.6 67 669.6 5 656.8 6 ๐28 590.1 15 581.0 21 560.0 24 554.5 79 ๐29 500.5 19 491.4 18 512.9 21 508.2 31 ๐30 475.1 7 471.5 13 345.9 48 344.8 39 ๐31 475.4 9 306.9 7 417.3 27 410.8 26 ๐31 146.0 74 437.8 16 317.7 5 308.6 23 ๐32 173.4 7 141.8 22 290.3 6 294.7 8 ๐33 345.4 21 194.7 58 216.3 7 213.1 9 ๐34 213.8 1 210.9 6 148.0 14 134.7 13 ๐35 98.9 4 101.3 4 139.2 28 120.3 37 ๐36 35.4 0 31.8 1 45.7 10 47.1 7 98 Table B.10 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformers 13 and 14, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. Mode Ser-13 Ser-14 MP2 B3LYP-D3BJ MP2 B3LYP-D3BJ ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐1 3698.7 46 3653.4 40 3679.0 33 3625.3 26 ๐2 3576.4 68 3536.7 53 3580.7 73 3540.7 60 ๐3 3460.0 9 3413.1 7 3483.2 13 3441.3 13 ๐4 3384.4 4 3364.8 3 3392.4 2 3367.5 2 ๐5 3005.7 18 2935.9 22 3028.9 15 2962.3 19 ๐6 2968.2 12 2895.1 18 3000.9 17 2958.3 16 ๐7 2919.5 11 2874.2 9 2848.6 40 2777.6 30 ๐8 1785.1 168 1766.6 213 1794.9 45 1783.8 158 ๐9 1597.0 21 1663.7 5 1611.7 14 1593.5 24 ๐10 1481.9 5 1465.6 7 1471.2 2 1456.4 2 ๐11 1409.1 1 1401.5 0 1400.7 4 1381.4 2 ๐12 1371.6 13 1355.8 0 1380.1 20 1367.9 6 ๐13 1339.0 17 1320.0 7 1352.9 9 1340.5 3 ๐15 1309.5 20 1290.3 15 1318.8 10 1303.8 3 ๐16 1283.0 17 1275.1 26 1292.9 0 1285.9 1 ๐17 1230.7 17 1218.8 12 1230.1 8 1222.5 6 ๐18 1122.8 9 1143.0 26 1166.6 28 1157.8 8 ๐19 1133.0 133 1104.7 143 1148.5 13 1129.2 16 ๐20 1169.6 8 1125.5 5 1108.9 245 1090.0 231 ๐21 1048.2 64 1031.2 72 1062.2 23 1032.2 44 ๐22 1023.7 20 998.1 29 1012.6 44 995.8 66 ๐23 910.0 64 884.1 1 930.8 10 901.1 21 ๐24 879.6 195 819.7 125 865.8 5 851.0 4 ๐25 828.3 14 809.0 42 789.2 120 720.4 20 ๐26 748.0 25 736.2 27 713.3 54 699.6 151 ๐27 640.9 112 634.9 102 663.1 19 658.1 18 ๐28 560.5 29 556.6 35 532.7 57 528.4 43 ๐29 524.6 13 518.3 21 492.5 4 491.4 11 ๐30 482.7 4 475.8 4 426.2 9 421.5 9 ๐31 375.6 4 362.3 2 361.5 5 364.8 4 ๐31 155.5 25 149.6 25 390.5 16 357.5 21 ๐32 299.3 5 342.7 7 143.5 38 126.8 33 ๐33 276.6 34 244.1 6 246.0 63 258.5 60 ๐34 186.7 131 194.0 90 202.9 1 205.8 6 ๐35 113.7 2 58.4 23 130.7 4 122.4 4 ๐36 46.3 4 40.9 2 48.5 0 51.7 0 99 Table B.11 Anharmonic vibrational frequencies (cm-1) and intensities (km mol-1) for ๐ผ-serine conformer 15, calculated at the MP2/cc-pVTZ and B3LYP-D3BJ/def2-TZVPD levels of theory. Mode Ser-15 - MP2 B3LYP-D3BJ MP2 B3LYP-D3BJ ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐/ cm-1 Int. ๐1 3696.9 46 3654.8 39 ๐2 3583.7 65 3538.3 46 ๐3 3441.2 9 3395.1 8 ๐4 3372.7 3 3346.4 2 ๐5 3005.2 20 2929.8 28 ๐6 2966.7 16 2894.2 13 ๐7 2918.4 9 2859.5 7 ๐8 1795.4 140 1782.4 100 ๐9 1605.2 22 1608.9 6 ๐10 1482.9 1 1467.9 6 ๐11 1412.6 1 1406.7 3 ๐12 1358.2 0 1360.0 0 ๐13 1373.8 9 1335.5 4 ๐15 1296.3 11 1270.5 1 ๐16 1277.3 1 1283.1 21 ๐17 1225.4 9 1219.6 9 ๐18 1170.8 2 1163.4 3 ๐19 1141.9 6 1119.5 8 ๐20 1101.8 211 1087.4 156 ๐21 1050.3 50 1030.9 50 ๐22 1024.7 16 1003.6 42 ๐23 902.7 118 894.1 4 ๐24 919.3 1 851.0 130 ๐25 834.0 13 817.9 11 ๐26 743.1 21 732.0 10 ๐27 638.2 19 626.4 31 ๐28 563.6 53 572.8 66 ๐29 508.8 13 500.4 21 ๐30 470.9 14 466.8 14 ๐31 349.6 47 352.7 31 ๐31 368.2 15 211.8 12 ๐32 148.7 32 366.2 24 ๐33 208.9 28 123.5 22 ๐34 199.3 90 178.5 98 ๐35 113.6 1 108.3 1 ๐36 60.9 1 55.9 1 "@en ;
edm:hasType "Thesis/Dissertation"@en ;
vivo:dateIssued "2021-05"@en ;
edm:isShownAt "10.14288/1.0396979"@en ;
dcterms:language "eng"@en ;
ns0:degreeDiscipline "Chemistry"@en ;
edm:provider "Vancouver : University of British Columbia Library"@en ;
dcterms:publisher "University of British Columbia"@en ;
dcterms:rights "Attribution-NonCommercial-NoDerivatives 4.0 International"@* ;
ns0:rightsURI "http://creativecommons.org/licenses/by-nc-nd/4.0/"@* ;
ns0:scholarLevel "Graduate"@en ;
dcterms:title "Exploring the conformational landscape of amino acids in solid parahydrogen matrices through ab initio calculations"@en ;
dcterms:type "Text"@en ;
ns0:identifierURI "http://hdl.handle.net/2429/77997"@en .