The gas-phase structures of cationized glycine (Gly), including complexes with Li + , Na + , K + , Rb + , and Cs + , are examined using infrared multiple-photon dissociation (IRMPD) spectroscopy utilizing light generated by a free electron laser, in conjunction with ab initio calculations. To identify the structures present in the experimental studies, measured IRMPD spectra are compared to spectra calculated at B3LYP/6-311+G(d,p) for the Li + , Na + , and K + complexes and at B3LYP/def2TZVP for the Rb + and Cs + complexes. Single-point energy calculations were carried out at the B3LYP, B3P86, and MP2(full) levels using the 6-311+G(2d,2p) basis set for Li + , Na + , K + and the def2TZVPP basis set for Rb + and Cs + . The Li + and Na + complexes are identified as metal cation coordination to the amino nitrogen and carbonyl oxygen, [N,CO]-tt, although Na + (Gly) may have contributions from additional structures. The heavier metal cations coordinate to either the carbonyl oxygen, [CO]-cc, or the carbonyl oxygen and hydroxy oxygen, [CO,OH]-cc, with the former apparently preferred for Rb + and Cs + and the latter for K + . These two structures reside in a double-well potential and different levels of theory predict very different relative stabilities. Some experimental evidence is provided that MP2(full) theory provides the most accurate relative energies.
more »
« less
Vibrational scaling factors for transition metal carbonyls
Vibrational frequencies for a selected set of transition metal carbonyl complexes are computed with various forms of density functional theory (B3LYP, BP86, M06, and M06-L), employing several different basis sets. The computed frequencies for the carbonyl stretches are compared to the experimental values obtained from gas phase infrared spectra of isolated neutrals and ions. Recommended carbonyl-stretch scaling factors which are developed vary significantly for different functionals, but there is little variation with basis set. Scaled frequencies compared to experimental spectra for cobalt and tantalum carbonyl cations reveal additional variations in multiplet patterns and relative band intensities for different functionals.
more »
« less
- Award ID(s):
- 1505095
- PAR ID:
- 10016700
- Date Published:
- Journal Name:
- Chemical physics letters
- Volume:
- 640
- ISSN:
- 0009-2614
- Page Range / eLocation ID:
- 175-179
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Jackson, George; Head-Gordon, Martin; Helgaker, Trygve; Liu, Wenjian; Osterwalder, Adreas (Ed.)We computed vibrational circular dichroism (VCD) and Raman optical activity (ROA) spectra for a test set of six chiral compounds using two standard density-functionals and an array of basis sets. We analysed the performance of property-oriented basis sets using a quadruple-zeta basis as a reference against four key metrics. We find little qualitative difference between the spectra produced by the larger basis sets (ORP, LPolX, aug-cc-pVTZ, and aug-cc-pVQZ), though their quantitative metrics exhibit wide variations. The smaller basis sets (rDPS, augD-3-21G, augT3-3-21G, Sadlej-pVTZ, and aug-cc-pVDZ) performed better for VCD rotatory strengths than for the corresponding ROA circular intensity differences (CIDs). However, this trend diminishes as the basis-set size is increased, lending validity to the conclusion that more robust property-oriented basis sets are required for ROA spectral generation than that of VCD. We observed improved performance in the mid-infrared region compared to the high-frequency regime, as well as overestimation of VCD rotatory strengths in the latter region as compared to the reference. We conclude that the ORP and LPol-ds basis sets are the most efficient and effective choices of basis set for the prediction of VCD and ROA spectra, as they provide both highly accurate results at reduced computational expense.more » « less
-
Abstract We have performed a series of highly accurate calculations between CO2and the 20 naturally occurring amino acids for the investigation of the attractive noncovalent interactions. Different nucleophilic groups present in the amino acid structures were considered (α‐NH2, COOH, side groups), and the stronger binding sites were identified. A database of accurate reference interactions energies was compiled as computed by explicitly‐correlated coupled‐cluster singles‐and‐doubles, together with perturbative triples extrapolated to the complete‐basis‐set limit. The CCSD(F12)(T)/CBS reference values were used for comparing a variety of popular density functionals with different basis sets. Our results show that most density functionals with the triple‐zeta basis set def2‐TZVPP align with the CCSD(F12)(T)/CBS reference values, but errors range from 0.1 kcal/mol up to 1.0 kcal/mol.more » « less
-
Complexes of 18-crown-6 ether (18C6) with four protonated amino acids (AAs) are examined using infrared multiple photon dissociation (IRMPD) action spectroscopy utilizing light generated by the infrared free electron laser at the Centre Laser Infrarouge d’Orsay (CLIO). The AAs examined in this work include glycine (Gly) and the three basic AAs: histidine (His), lysine (Lys), and arginine (Arg). To identify the (AA)H + (18C6) conformations present in the experimental studies, the measured IRMPD spectra are compared to spectra calculated at the B3LYP/6-311+G(d,p) level of theory. Relative energies of various conformers and isomers are provided by single point energy calculations carried out at the B3LYP, B3P86, M06, and MP2(full) levels using the 6-311+G(2p,2d) basis set. The comparisons between the IRMPD and theoretical IR spectra indicate that 18C6 binds to Gly and His via the protonated backbone amino group, whereas protonated Lys prefers binding via the protonated side-chain amino group. Results for Arg are less definitive with strong evidence for binding to the protonated guanidino side chain (the calculated ground conformer at most levels of theory), but contributions from backbone binding to a zwitterionic structure are likely.more » « less
-
Our ability to understand and simulate the reactions catalyzed by iron depends strongly on our ability to predict the relative energetics of spin states. In this work, we studied the electronic structures of Fe 2+ ion, gaseous FeO and 14 iron complexes using Kohn–Sham density functional theory with particular focus on determining the ground spin state of these species as well as the magnitudes of relevant spin-state energy splittings. The 14 iron complexes investigated in this work have hexacoordinate geometries of which seven are Fe( ii ), five are Fe( iii ) and two are Fe( iv ) complexes. These are calculated using 20 exchange–correlation functionals. In particular, we use a local spin density approximation (LSDA) – GVWN5, four generalized gradient approximations (GGAs) – BLYP, PBE, OPBE and OLYP, two non-separable gradient approximations (NGAs) – GAM and N12, two meta-GGAs – M06-L and M11-L, a meta-NGA – MN15-L, five hybrid GGAs – B3LYP, B3LYP*, PBE0, B97-3 and SOGGA11-X, four hybrid meta-GGAs – M06, PW6B95, MPW1B95 and M08-SO and a hybrid meta-NGA – MN15. The density functional results are compared to reference data, which include experimental results as well as the results of diffusion Monte Carlo (DMC) calculations and ligand field theory estimates from the literature. For the Fe 2+ ion, all functionals except M11-L correctly predict the ground spin state to be quintet. However, quantitatively, most of the functionals are not close to the experimentally determined spin-state splitting energies. For FeO all functionals predict quintet to be the ground spin state. For the 14 iron complexes, the hybrid functionals B3LYP, MPW1B95 and MN15 correctly predict the ground spin state of 13 out of 14 complexes and PW6B95 gets all the 14 complexes right. The local functionals, OPBE, OLYP and M06-L, predict the correct ground spin state for 12 out of 14 complexes. Two of the tested functionals are not recommended to be used for this type of study, in particular M08-SO and M11-L, because M08-SO systematically overstabilizes the high spin state, and M11-L systematically overstabilizes the low spin state.more » « less
An official website of the United States government

