Near-Infrared Laser Induced Conformational Change of Alanine in Low-Temperature Matrixes and the Tunneling Lifetime of Its Conformer VI

NIR and UV induced transformations of indazole-3-carboxylic acid isolated in low temperature matrices

The molecular structure and NIR and UV induced phototransformations of indazole-3-carboxylic acid were studied in low temperature argon and nitrogen matrices by FTIR spectroscopy and B3LYP/6-311++G(2d,2p) calculations. Eleven minima of IC were located on the S₀ potential energy surface. The three most stable structures were detected experimentally in both matrices after deposition. Upon NIR irradiation a dominant process was rotamerization within the carboxylic group leading to changes in the population of the trans 1HIC1 and cis 1HIC2 forms. In turn, at UV irradiation at λ = 260 nm two new tautomers (2HIC2 and 2HIC3) were generated indicating that the hydrogen atom transfer in pyrazole ring took place. Based on the obtained kinetic curves, differences in the phototransformation rates in different matrices were indicated.

Infrared Spectra and Phototransformations of meta-Fluorophenol Isolated in Argon and Nitrogen Matrices

Monomers of meta-fluorophenol (mFP) were trapped from the gas phase into cryogenic argon and nitrogen matrices. The estimated relative energies of the two conformers are very close, and in the gas phase they have nearly equal populations. Due to the similarity of their structures (they only differ in the orientation of the OH group), the two conformers have also similar predicted vibrational signatures, which makes the vibrational characterization of the individual rotamers challenging. In the present work, it has been established that in an argon matrix only the most stable trans conformer of mFP exists (the OH group pointing away from the fluorine atom). On the other hand, the IR spectrum of mFP in a nitrogen matrix testifies to the simultaneous presence in this matrix of both the trans conformer and of the higher-energy cis conformer (the OH group pointing toward the fluorine atom), which is stabilized by interaction with the matrix gas host. We found that the exposition of the cryogenic N₂ matrix to the Globar source of the infrared spectrometer affects the conformational populations. By collecting experimental spectra, either in the full mid-infrared range or only in the range below 2200 cm^-1, we were able to reliably distinguish two sets of experimental bands originating from individual conformers. A comparison of the two sets of experimental bands with computed infrared spectra of the conformers allowed, for the first time, the unequivocal vibrational identification of each of them. The joint implementation of computational vibrational spectroscopy and matrix-isolation infrared spectroscopy proved to be a very accurate method of structural analysis. Some mechanistic insights into conformational isomerism (the quantum tunneling of hydrogen atom and vibrationally-induced conformational transformations) have been addressed. Finally, we also subjected matrix-isolated mFP to irradiations with UV light, and the phototransformations observed in these experiments are also described.

IR-induced and tunneling reactions in cryogenic matrices: the (incomplete) story of a successful endeavor

In this article, IR-induced and tunneling-driven reactions observed in cryogenic matrices are described in a historical perspective, the entangling of the two types of processes being highlighted. The story of this still ongoing fascinating scientific endeavor is presented here following closely our own involvement in the field for more than 30 years, and thus focuses mostly on our work. It is, because of this reason, also an incomplete story. Nevertheless, it considers a large range of examples, from very selective IR-induced conformational isomerizations to IR-induced bond-breaking/bond-forming reactions and successful observations of rare heavy atom tunneling processes. As a whole, this article provides a rather general overview of the major progress achieved in the field.

Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case

The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)^MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm^-1, has been obtained by combining CCSD(T)/(CBS+CV)^MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.

Conformational Space, IR-Induced, and UV-Induced Chemistry of Carvacrol Isolated in a Low-Temperature Argon Matrix

In this work, monomers of carvacrol (5-isopropyl-2-methylphenol), a natural monoterpene exhibiting wide range bioactivity, were trapped in a cryogenic argon matrix and characterized by infrared spectroscopy, while quantum chemical calculations at the B3LYP and MP2 levels were employed to characterize the conformational landscape of the isolated molecule. Four conformers have been localized on the potential energy surface, and the factors accounting for their relative stability were analyzed. The two most stable conformers of carvacrol, differing in the relative orientation of the isopropyl group and both having the OH group pointing away from the vicinal methyl fragment, were identified in the cryomatrix for the first time. The individual spectral signatures of the two conformers were distinguished based on the change in their relative abundance induced by exposing the matrix to broadband infrared light. Matrix-isolated carvacrol was also irradiated with broadband UV light (λ > 200 nm), which resulted in the cleavage of the OH group. Recombination of the released H atom at the ortho- or para-position of the ring resulted in the formation of alkyl-substituted cyclohexadienones. These were found to undergo subsequent valence and open-ring isomerizations, leading, respectively, to the formation of a Dewar isomer and open-chain conjugated ketenes. Decarbonylation of the photoproducts was also observed for longer irradiation times. A mechanistic analysis of the observed photochemical transformations is presented.

Bond-Breaking/Bond-Forming Reactions by Vibrational Excitation: Infrared-Induced Bidirectional Tautomerization of Matrix-Isolated Thiotropolone

Infrared vibrational excitation is a promising approach for gaining exceptional control of chemical reactions, in ways that cannot be attained via thermal or electronic excitation. Here, we report an unprecedented example of a bond-breaking/bond-forming reaction by vibrational excitation under matrix isolation conditions. Thiotropolone monomers were isolated in cryogenic argon matrices and characterized by infrared spectroscopy and vibrational computations (harmonic and anharmonic). Narrowband near-infrared irradiations tuned at frequencies of first CH stretching overtone (5940 cm^-1) or combination modes (5980 cm^-1) of the OH tautomer, the sole form of the compound that exists in the as-deposited matrices, led to its conversion into the SH tautomer. The tautomerization in the reverse direction was achieved by vibrational excitation of the SH tautomer with irradiation at 5947 or 5994 cm^-1, corresponding to the frequencies of its CH stretching combination and overtone modes. This pioneer demonstration of bidirectional hydroxyl ↔ thiol tautomerization controlled by vibrational excitation creates prospects for new advances in vibrationally induced chemistry.

Conformational isomerizations triggered by vibrational excitation of second stretching overtones

Vibrational excitation using frequency-tunable IR laser light has been developed as a powerful tool for selective manipulation of molecular conformations. In this methodology, vibrational excitation has been typically applied to the first stretching overtones (∼80 kJ mol^-1) but also to the fundamental modes (∼40 kJ mol^-1). Here, we demonstrate that selective conformational isomerizations are also achieved using excitation to second stretching overtones (∼120 kJ mol^-1). The extremely weak absorptions of the second stretching overtones of molecules isolated in low-temperature matrices were measured for the first time; here using three prototype molecules: hydroxyacetone (HA), glycolic acid (GAc) and glycolamide (GAm). Benchmarking of computed anharmonic IR spectra showed that the B3LYP/SNSD method provides the best agreement with experimental frequencies of the ν(OH), 2ν(OH) and 3ν(OH) modes for the studied molecules in argon matrices. Selective irradiation at the 3ν(OH) frequencies (9850-10 500 cm^-1) of HA, GAc and GAm monomers in argon matrices at 15 K successfully triggers their conformational isomerization. These results open the door to extend control over conformations separated by higher barriers and to induce other transformations not energetically accessible by excitation to the fundamental or first stretching overtone modes.

UV-promoted radical formation, and near-IR-induced and spontaneous conformational isomerization in monomeric 9-methylguanine isolated in low-temperature Ar matrices

Three low-energy isomers of 9-methylguanine, the amino-oxo (AO) form and two amino-hydroxy (AH1 and AH2) conformers, were trapped from the gas phase into low-temperature argon matrices. The AH1 and AH2 isomers, differing in the orientation of the OH group, were found to transform into each other upon excitation with near-IR light. The population of the AO form of the compound was not changed upon any near-IR irradiation of the matrix samples. Using monochromatic near-IR light, generated by a frequency-tunable laser source, it was possible to selectively induce the AH1 → AH2 or AH2 → AH1 conversion. Photoreversibility of this conformational transformation was then demonstrated. Exposure of matrix-isolated monomers of 9-methylguanine to broadband near-IR light also led to conformational conversions within the amino-hydroxy tautomeric form; the final stage of this process was always the same photostationary state independent of the initial ratio of AH1 and AH2 populations. Spontaneous conformational conversion, transforming the higher-energy AH2 form into the lower-energy AH1 isomer, was observed for matrix-isolated monomers of 9-methylguanine kept in the dark. The mechanism of this process must rely on quantum tunneling of the light hydrogen atom. Irradiation of matrix-isolated 9-methylguanine with UV laser light at λ = 288 or 285 nm led to a substantial consumption of the two AH forms, while the amount of AO isomer remained unchanged. On the other hand, a decrease in the population of the AO isomer occurred upon excitations at shorter wavelengths, λ = 280 or 275 nm. The spectral changes observed after UV-irradiation suggest the generation (and stabilization in the matrix) of a radical species, resulting from the photocleavage of the O-H or N1-H bonds, in the AH or AO isomer, respectively.

Shape and core-excited resonances in electron scattering from alanine

We present detailed ab initio scattering calculations using the R-matrix method for electron collisions with the most stable conformer of α-alanine. The shape resonances that we identify are in good agreement with earlier calculations and experiments. Core-excited and mixed-character resonances are identified and characterized computationally for the first time. Dissociative electron attachment results are discussed in relation to the resonances identified.

S–H rotamerizationviatunneling in a thiol form of thioacetamide

Rotamerization of the S–H groupviahydrogen tunneling is reported for the first time.

Nunes C, Tarczay G, Fausto R. Selective conformational control by excitation of NH imino vibrational antennas

We provide experimental evidence for the occurrence of selective and reversible conformational control over the SH group by vibrational excitation of remote NH groups. Using an imino group that acts as a molecular antenna has no precedents.

Intricate Conformational Tunneling in Carbonic Acid Monomethyl Ester

Disentangling internal and external effects is a key requirement for understanding conformational tunneling processes. Here we report the s- trans/ s- cis tunneling rotamerization of carbonic acid monomethyl ester (1) under matrix isolation conditions and make comparisons to its parent carbonic acid (3). The observed tunneling rate of 1 is temperature-independent in the 3-20 K range and accelerates when using argon instead of neon as the matrix material. The methyl group increases the effective half life (τ_eff) of the energetically disfavored s- trans-conformer from 3-5 h for 3 to 11-13 h for 1. Methyl group deuteration slows the rotamerization further (τ_eff ≈ 35 h). CCSD(T)/cc-pVQZ//MP2/aug-cc-pVTZ computations of the tunneling probability suggest that the rate should be almost unaffected by methyl substitution or its deuteration. Thus the observed relative rates are puzzling, and they disagree with previous explanations involving fast vibrational relaxation after the tunneling event facilitated by the alkyl rotor.

Nunes C, Fausto R, Reva I. Conformational control over an aldehyde fragment by selective vibrational excitation of interchangeable remote antennas

We apply vibrational antennas (OH or NH₂ group) to achieve unprecedented conformational control over the heavy aldehyde fragment in 2-formyl-2H-azirine, using selective vibrational excitations of the OH or NH₂ stretching overtones and combination modes.

Hydrogen-atom tunneling through a very high barrier; spontaneous thiol → thione conversion in thiourea isolated in low-temperature Ar, Ne, H2and D2matrices

Hydrogen-atom tunneling in thiourea involves breaking of the S–H chemical bond and the formation of a new N–H bond.

Computational Chemistry: The Fate of Current Methods and Future Challenges

"Where do we go from here?" is the underlying question regarding the future (perhaps foreseeable) developments in computational chemistry. Although this young discipline has already permeated practically all of chemistry, it is likely to become even more powerful with the rapid development of computational hard- and software.

Experimental Evidence of Long-Range Intramolecular Vibrational Energy Redistribution through Eight Covalent Bonds: NIR Irradiation Induced Conformational Transformation of E-Glutaconic Acid

Long-range intramolecular vibrational energy redistribution (IVR) driven conformational changes were investigated in a matrix-isolated open-chain, asymmetrical dicarboxylic acid, E-glutaconic acid. Although the analysis was challenging due to the presence of multiple backbone conformers and short lifetimes of the prepared higher energy cis conformers, it was shown that the selective excitation of the O-H stretching overtone of one of the carboxylic groups can induce the conformational change (trans to cis) of the other carboxylic group, located at the other end of the E-glutaconic acid molecule. This is a direct proof that the IVR process can act through eight covalent bonds in a flexible molecule before the excess energy completely dissipates into the matrix. The lifetime of the prepared higher energy conformers (averaged over the different backbones) was measured to be 12 s.

Thermally and vibrationally induced conformational isomerizations, infrared spectra, and photochemistry of gallic acid in low-temperature matrices

Near-infrared (near-IR) narrowband selective vibrational excitation and annealing of gallic acid (3,4,5-trihydroxybenzoic acid) isolated in cryogenic matrices were used to induce interconversions between its most stable conformers. The isomerizations were probed by infrared spectroscopy. An extensive set of quantum chemical calculations, carried out at the DFT(B3LYP)/6-311++G(d,p) level of approximation, was used to undertake a detailed analysis of the ground state potential energy surface of the molecule. This investigation of the molecule conformational space allowed extracting mechanistic insights into the observed annealing- or near-IR-induced isomerization processes. The infrared spectra of the two most stable conformers of gallic acid in N2, Xe, and Ar matrices were fully assigned. Finally, the UV-induced photochemistry of the matrix isolated compound was investigated.

Evidence of Coupling between the Motions of Water and Peptides

Studies of protein dynamics at low temperatures are generally performed on hydrated powders and not in biologically realistic solutions of water because of water crystallization. However, here we avoid the problem of crystallization by reducing the size of the biomolecules. We have studied oligomers of the amino acid l-lysine, fully dissolved in water, and our dielectric relaxation data show that the glass transition-related dynamics of the oligomers is determined by the water dynamics, in a way similar to that previously observed for solvated proteins. This implies that the crucial role of water for protein dynamics can be extended to other types of macromolecular systems, where water is also able to determine their conformational fluctuations. Using the energy landscape picture of macromolecules, the thermodynamic criterion for such solvent-slaved macromolecular motions may be that the macromolecules need the entropy contribution from the solvent to overcome the enthalpy barriers between different conformational substates.

Conformational changes in matrix-isolated 6-methoxyindole: Effects of the thermal and infrared light excitations

Conformational changes induced thermally or upon infrared excitation of matrix-isolated 6-methoxyindole were investigated. Narrowband near-infrared excitation of the first overtone of the N-H stretching vibration of each one of the two identified conformers is found to induce a selective large-scale conversion of the pumped conformer into the other one. This easily controllable bidirectional process consists in the intramolecular reorientation of the methoxy group and allowed a full assignment of the infrared spectra of the two conformers. Matrices with different conformational compositions prepared by narrow-band irradiations were subsequently used to investigate the effects of both thermal and broadband infrared excitations on the conformational mixtures. Particular attention is given to the influence of the matrix medium (Ar vs. Xe) and conformational effects of exposition of the sample to the spectrometer light source during the measurements.

Tunnelling in carbonic acid

The cis,trans-conformer of carbonic acid (H₂CO₃), generated by near-infrared radiation, undergoes an unreported quantum mechanical tunnelling rotamerization with half-lives in cryogenic matrices of 4–20 h, depending on temperature and host material.

Conformers, infrared spectrum, UV-induced photochemistry, and near-IR-induced generation of two rare conformers of matrix-isolated phenylglycine

The conformational space of α-phenylglycine (PG) have been investigated theoretically at both the DFT/B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) levels of approximation. Seventeen different minima were found on the investigated potential energy surfaces, which are characterized by different dominant intramolecular interactions: type I conformers are stabilized by hydrogen bonds of the type N-H···O=C, type II by a strong O-H···N hydrogen bond, type III by weak N-H···O-H hydrogen bonds, and type IV by a C=O···H-C contact. The calculations indicate also that entropic effects are relevant in determining the equilibrium populations of the conformers of PG in the gas phase, in particular in the case of conformers of type II, where the strong intramolecular O-H···N hydrogen bond considerably diminishes entropy by reducing the conformational mobility of the molecule. In consonance with the relative energies of the conformers and barriers for conformational interconversion, only 3 conformers of PG were observed for the compound isolated in cryogenic Ar, Xe, and N2 matrices: the conformational ground state (ICa), and forms ICc and IITa. All other significantly populated conformers existing in the gas phase prior to deposition convert either to conformer ICa or to conformer ICc during matrix deposition. The experimental observation of ICc had never been achieved hitherto. Narrowband near-IR irradiation of the first overtone of νOH vibrational mode of ICa and ICc in nitrogen matrices (at 6910 and 6930 cm(-1), respectively) led to selective generation of two additional conformers of high-energy, ITc and ITa, respectively, which were also observed experimentally for the first time. In addition, these experiments also provided the key information for the detailed vibrational characterization of the 3 conformers initially present in the matrices. On the other hand, UV irradiation (λ = 255 nm) of PG isolated in a xenon matrix revealed that PG undergoes facile photofragmentation through two photochemical pathways that are favored for different initial conformations of the reactant: (a) decarboxylation, leading to CO2 plus benzylamine (the dominant photofragmentation channel in PG cis-COOH conformers ICa and ICc) and (b) decarbonylation, with generation of CO plus benzonitrile, H2O and H2 (prevalent in the case of the trans-COOH conformer, IITa).

CC/DFT Route toward Accurate Structures and Spectroscopic Features for Observed and Elusive Conformers of Flexible Molecules: Pyruvic Acid as a Case Study

The structures and relative stabilities as well as the rotational and vibrational spectra of the three low-energy conformers of pyruvic acid (PA) have been characterized using a state-of-the-art quantum-mechanical approach designed for flexible molecules. By making use of the available experimental rotational constants for several isotopologues of the most stable PA conformer, Tc-PA, the semiexperimental equilibrium structure has been derived. The latter provides a reference for the pure theoretical determination of the equilibrium geometries for all conformers, thus confirming for these structures an accuracy of 0.001 Å and 0.1 deg for bond lengths and angles, respectively. Highly accurate relative energies of all conformers (Tc-, Tt-, and Ct-PA) and of the transition states connecting them are provided along with the thermodynamic properties at low and high temperatures, thus leading to conformational enthalpies accurate to 1 kJ mol(-1). Concerning microwave spectroscopy, rotational constants accurate to about 20 MHz are provided for the Tt- and Ct-PA conformers, together with the computed centrifugal-distortion constants and dipole moments required to simulate their rotational spectra. For Ct-PA, vibrational frequencies in the mid-infrared region accurate to 10 cm(-1) are reported along with theoretical estimates for the transitions in the near-infrared range, and the corresponding infrared spectrum including fundamental transitions, overtones, and combination bands has been simulated. In addition to the new data described above, theoretical results for the Tc- and Tt-PA conformers are compared with all available experimental data to further confirm the accuracy of the hybrid coupled-cluster/density functional theory (CC/DFT) protocol applied in the present study. Finally, we discuss in detail the accuracy of computational models fully based on double-hybrid DFT functionals (mainly at the B2PLYP/aug-cc-pVTZ level) that avoid the use of very expensive CC calculations.

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

The conformers of α-serine were investigated by matrix-isolation IR spectroscopy combined with NIR laser irradiation. This method, aided by 2D correlation analysis, enabled unambiguously grouping the spectral lines to individual conformers. On the basis of comparison of at least nine experimentally observed vibrational transitions of each conformer with empirically scaled (SQM) and anharmonic (GVPT2) computed IR spectra, six conformers were identified. In addition, the presence of at least one more conformer in Ar matrix was proved, and a short-lived conformer with a half-life of (3.7 ± 0.5) × 10(3) s in N2 matrix was generated by NIR irradiation. The analysis of the NIR laser-induced conversions revealed that the excitation of the stretching overtone of both the side chain and the carboxylic OH groups can effectively promote conformational changes, but remarkably different paths were observed for the two kinds of excitations.

Understanding the conformational behaviour of Ac-Ala-NHMe in different media. A joint NMR and DFT study

The conformational behaviour of Ac-Ala-NHMe was studied in the gas-phase and in solution by theoretical calculations (B3LYP-D3/aug-cc-pVDZ level) and experimental (1)H NMR. The conformational preferences of this compound were shown to result from a complex interplay between the strengths of possible intramolecular hydrogen bonds, steric interactions, hyperconjugation, entropy effects and the overall dipole moments. The Ac-Ala-N(Me)2 derivative was studied in addition, to design a system akin to Ac-Ala-NHMe, but with disrupted intramolecular hydrogen bonds involving the -NHMe group, mimicking the effect of polar protic solvents.

Conformational switching in pyruvic acid isolated in Ar and N₂ matrixes: spectroscopic analysis, anharmonic simulation, and tunneling

Monomers of pyruvic acid (PA) isolated in cryogenic argon and nitrogen matrixes were characterized by mid- and near-infrared spectroscopy. Interpretation of the experiments was aided by fully anharmonic calculations of the fundamental modes, overtones, and combinations up to two quanta, including their infrared intensities. The initially dominating PA conformer (Tc) has a cis CCOH arrangement and is stabilized by a strong intramolecular H-bond. Selective near-infrared excitation of Tc at the first OH overtone (6630 cm(-1) in Ar, 6643 cm(-1) in N2) induced a large scale conformational conversion to the higher-energy conformer (Tt) with trans CCOH arrangement. Tt was then converted back to Tc by selective NIR irradiation at the first Tt OH overtone (6940 cm(-1) in Ar, 6894 cm(-1) in N2). In N2 matrix, the Tt form was stabilized due to interaction between the OH group and the matrix molecules. This stabilization manifested itself in the absence of Tt → Tc relaxation and in a considerable change of the vibrational Tt signature upon going from argon to nitrogen matrix. In argon, the Tt form spontaneously decayed back to Tc in the dark (characteristic lifetime +16 h). In the presence of broad-band near-infrared light, the Tt → Tc relaxation speed considerably increased. The decay mechanisms are discussed.

Near-infrared laser-induced generation of three rare conformers of glycolic acid

Structural transformations were induced in conformers of glycolic acid by selective excitation with monochromatic tunable near-infrared laser light. For the compound isolated in Ar matrixes, near-IR excitation led to generation of two higher-energy conformers (GAC; AAT) differing from the most stable SSC form by 180° rotation around the C-C bond. A detailed investigation of this transformation revealed that one conformer (GAC) is produced directly from the near-IR-excited most stable conformer. The other higher-energy conformer (AAT) was effectively generated only upon excitation of the primary photoproduct (GAC) with another near-IR photon. Once these higher-energy conformers of glycolic acid were generated in an Ar matrix, they could be subsequently transformed into one another upon selective near-IR excitations. Interestingly, no repopulation of the initial most stable SSC conformer occurred upon near-IR excitation of the higher-energy forms of the compound isolated in solid Ar. A dramatically different picture of near-IR-induced conformational transformations was observed for glycolic acid isolated in N2 matrixes. In this case, upon near-IR excitation, the most stable SSC form converted solely into a new conformer (SST), where the acid OH group is rotated by 180°. This conformational transformation was found to be photoreversible. Moreover, SST conformer, photoproduced in the N2 matrix, spontaneously converted to the most stable SSC form of glycolic acid, when the matrix was kept at cryogenic temperature and in the dark.

Performance of M06, M06-2X, and M06-HF density functionals for conformationally flexible anionic clusters: M06 functionals perform better than B3LYP for a model system with dispersion and ionic hydrogen-bonding interactions

We present a comparative assessment of the performance of the M06 suite of density functionals (M06, M06-2X, and M06-HF) against an MP2 benchmark for calculating the relative energies and geometric structures of the Cl(-)·arginine and Br(-)·arginine halide ion-amino acid clusters. Additional results are presented for the popular B3LYP density functional. The Cl(-)·arginine and Br(-)·arginine complexes are important prototypes for the phenomenon of anion-induced zwitterion formation. Results are presented for the canonical (noncharge separated) and zwitterionic (charge separated) tautomers of the clusters, as well as the numerous conformational isomers of the clusters. We find that all of the M06 functions perform well in terms of predicting the general trends in the conformer relative energies and identifying the global minimum conformer. This is in contrast to the B3LYP functional, which performed significantly less well for the canonical tautomers of the clusters where dispersion interactions contribute more significantly to the conformer energetics. We find that the M06 functional gave the lowest mean unsigned error for the relative energies of the canonical conformers (2.10 and 2.36 kJ/mol for Br(-)·arginine and Cl(-)·arginine), while M06-2X gave the lowest mean unsigned error for the zwitterionic conformers (0.85 and 1.23 kJ/mol for Br(-)·arginine and Cl(-)·arginine), thus providing insight into the types of physical systems where each of these functionals should perform best.