Fundamental vibrational frequencies of small polyatomic molecules from density-functional calculations and vibrational perturbation theory

Unusual Electronic Structures of an Electron Transfer Series of [Cr(μ-η1 : η1 -N2 )Cr]0/1+/2

In dinitrogen (N2 ) fixation chemistry, bimetallic end-on bridging N2 complexes M(μ-η1  : η1 -N2 )M can split N2 into terminal nitrides and hence attract great attention. To date, only 4d and 5d transition complexes, but none of 3d counterparts, could realize such a transformation. Likewise, complexes {[Cp*Cr(dmpe)]2 (μ-N2 )}0/1+/2+ (1-3) are incapable to cleave N2 , in contrast to their Mo congeners. Remarkably, cross this series the N-N bond length of the N2 ligand and the N-N stretching frequency exhibit unprecedented nonmonotonic variations, and complexes 1 and 2 in both solid and solution states display rare thermally activated ligand-mediated two-center spin transitions, distinct from discrete dinuclear spin crossovers. In-depth analyses using wave function based ab initio calculations reveal that the Cr-N2 -Cr bonding in complexes 1-3 is distinguished by strong multireference character and cannot be described by solely one electron configuration or Lewis structure, and that all intriguing spectroscopic observations originate in their sophisticate multireference electronic structures. More critical is that such multireference bonding of complexes 1-3 is at least a key factor that contributes to their kinetic inertness toward N2 splitting. The mechanistic understanding is then used to rationalize the disparate reactivity of related 3d M(μ-η1  : η1 -N2 )M complexes compared to their 4d and 5d analogs.

Isomerism effects in relaxation dynamics of Au24(SR)16thiolate-protected gold nanoclusters

Understanding the excited state behavior of isomeric structures of thiolate-protected gold nanoclusters is still a challenging task. In this paper, based on grand unified model and ring model for describing thiolate-protected gold nanoclusters, we have predicted four isomers of Au₂₄(SR)₁₆nanoclusters. Density functional theory calculations show that the total energy of one of the predicted isomers is 0.1 eV lower in energy than previously crystallized isomer. The nonradiative relaxation dynamics simulations of Au₂₄(SH)₁₆isomers are performed to reveal the effects of structural isomerism on relaxation process of the lowest energy states, in which that most of the low-excited states consist of core states. In addition, crystallized isomer possesses the shorter e-h recombination time, whereas the most stable isomer has the longer recombination time, which may be attributed to the synergistic effect of nonadiabatic coupling and decoherence time. Our results could provide practical guidance to predict new gold nanoclusters for future experimental synthesis, and stimulate the exploration of atomic structures of same sized gold nanoclusters for photovoltaic and optoelectronic devices.

Comparative Analysis of Uncoupled Mode Approximations for Molecular Thermochemistry and Kinetics

The accurate prediction of thermochemistry and kinetic parameters is an important task for reaction modeling. Unfortunately, the commonly used harmonic oscillator model is often not accurate enough due to the absence of anharmonic effects. In this work, we improve the representation of an anharmonic potential energy surface (PES) using uncoupled mode (UM) approximations, which model the full-dimensional PES as a sum of one-dimensional potentials of each mode. We systematically analyze different PES sampling schemes and coordinate systems for constructing the one-dimensional potentials, and benchmark the performance of UM methods on data sets of molecular thermochemistry and kinetic properties. The results show that the accuracy of the UM approach strongly depends on how the one-dimensional potentials are defined. If one-dimensional potentials are constructed by sampling along normal mode directions (UM-N) or along the directions that minimize intermode coupling (E- and E'-optimized), the accuracies of the predicted properties are not significantly improved compared to the harmonic oscillator model. However, significant improvements can be achieved by sampling the torsional modes separately from the vibrational modes (UM-T and UM-VT). In this work, three types of coordinate systems are examined, including redundant internal coordinates (RIC), hybrid internal coordinates (HIC), and translation-rotation-internal coordinates (TRIC). The HIC and TRIC coordinate systems can outperform RIC since transition state species may contain large-amplitude interfragmentary motions that regular internal coordinates can not describe adequately. Among all the methods we examined, the activation energies and pre-exponential factors calculated using UM-VT with either TRIC or HIC best agree with the reference values. Since UM-VT requires only a number of additional single point energy calculations for each independent mode, the scaling of computational costs of UM-VT is the same as that of the standard harmonic oscillator model, making UM-VT an appealing way of calculating the thermochemistry and kinetic properties for large-size systems.

Hydrogen Atom Transfer Thermodynamics of Homologous Co(III)- and Mn(III)-Superoxo Complexes: The Effect of the Metal Spin State

Systematic investigations on H atom transfer (HAT) thermodynamics of metal O₂ adducts is of fundamental importance for the design of transition metal catalysts for substrate oxidation and/or oxygenation directly using O₂. Such work should help elucidate underlying electronic-structure features that govern the OO-H bond dissociation free energies (BDFEs) of metal-hydroperoxo species, which can be used to quantitatively appraise the HAT activity of the corresponding metal-superoxo complexes. Herein, the BDFEs of two homologous Co^III- and Mn^III-hydroperoxo complexes, 3-Co and 3-Mn, were calculated to be 79.3 and 81.5 kcal/mol, respectively, employing the Bordwell relationship based on experimentally determined pK _a values and redox potentials of the one-electron-oxidized forms, 4-Co and 4-Mn. To further verify these values, we tested the HAT capability of their superoxo congeners, 2-Co and 2-Mn, toward three different substrates possessing varying O-H BDFEs. Specifically, both metal-superoxo species are capable of activating the O-H bond of 4-oxo-TEMPOH with an O-H BDFE of 68.9 kcal/mol, only 2-Mn is able to abstract a H atom from 2,4-di-tert-butylphenol with an O-H BDFE of 80.9 kcal/mol, and neither of them can react with 3,5-dimethylphenol with an O-H BDFE of 85.6 kcal/mol. Further computational investigations suggested that it is the high spin state of the Mn^III center in 3-Mn that renders its OO-H BDFE higher than that of 3-Co, which features a low-spin Co^III center. The present work underscores the role of the metal spin state being as crucial as the oxidation state in modulating BDFEs.

Interfacing single-atom catalysis with continuous-flow organic electrosynthesis

The global warming crisis has sparked a series of environmentally cautious trends in chemistry, allowing us to rethink the way we conduct our synthesis, and to incorporate more earth-abundant materials in our catalyst design. "Single-atom catalysis" has recently appeared on the catalytic spectrum, and has truly merged the benefits that homogeneous and heterogeneous analogues have to offer. Further still, the possibility to activate these catalysts by means of a suitable electric potential could pave the way for a true integration of diverse synthetic methodologies and renewable electricity. Despite their esteemed benefits, single-atom electrocatalysts are still limited to the energy sector (hydrogen evolution reaction, oxygen reduction, etc.) and numerous examples in the literature still invoke the use of precious metals (Pd, Pt, Ir, etc.). Additionally, batch electroreactors are employed, which limit the intensification of such processes. It is of paramount importance that the field continues to grow in a more sustainable direction, seeking new ventures into the space of organic electrosynthesis and flow electroreactor technologies. In this piece, we discuss some of the progress being made with earth abundant homogeneous and heterogeneous electrocatalysts and flow electrochemistry, within the context of organic electrosynthesis, and highlight the prospects of alternatively utilizing single-atom catalysts for such applications.

Modeling Gas-Phase Unimolecular Dissociation for Bond Dissociation Energies: Comparison of Statistical Rate Models within RRKM Theory

The Rice-Ramsperger-Kassel-Marcus (RRKM) theory provides a simple yet powerful rate theory for calculating microcanonical rate constants. In particular, it has found widespread use in combination with gas-phase kinetic experiments of unimolecular dissociations to extract experimental bond dissociation energies (BDEs). We have previously found several discrepancies between the computed BDE values and the respective experimental ones, obtained with our empirical rate model, named L-CID. To investigate the reliability of our rate model, we conducted a theoretical analysis and comparison of the performance of conventional rate models and L-CID within the RRKM framework. Using the previously published microcanonical rate data as well as reaction cross-section data, we show that the BDE values obtained with the L-CID model agree with the ones from the other rate models within the expected uncertainty bounds. Based on this agreement, we discuss the possible rationalization of the good performance of the L-CID model.

Synthesis, crystal structure, Hirshfeld surface analysis and energy framework calculations of trans-3,7,9,9-tetra-methyl-10-(prop-2-yn-1-yl)-1,2,3,4,4a,9,9a,10-octa-hydro-acridine

The title heterocyclic compound, C20H27N, has been prepared in good yield (72%) via a BiCl3-catalyzed cationic Povarov reaction between N-propargyl-4-methyl-aniline and (±)-citronellal. The X-ray single-crystal study indicates that the structure consists of mol-ecules connected by C-H⋯π contacts to produce chains, which pack in a sandwich-herringbone fashion along the b-axis direction. Hirshfeld surface analysis indicates that H⋯H inter-actions dominate by contributing 79.1% to the total surface. Energy frameworks and DFT calculations indicate a major contribution of dispersive forces to the total inter-action energy.

Quantitative Comparison of Experimental and Computed IR-Spectra Extracted from Ab Initio Molecular Dynamics

Experimentally measured infrared spectra are often compared to their computed equivalents. However, the accordance is typically characterized by visual inspection, which is prone to subjective judgment. The primary challenge for a similarity-based analysis is that the artifacts introduced by each approach are very different and, therefore, may require preprocessing steps to determine and correct impeding irregularities. To allow for automated objective assessment, we propose a practical and comprehensive workflow involving scaling factors, a novel baseline correction scheme, and peak smoothing. The resulting spectra can then easily be compared quantitatively using similarity measures, for which we found the Pearson correlation coefficient to be the most suitable. The proposed procedure is then applied to compare the agreement of the experimental infrared spectra from the NIST Chemistry Web book with the calculated spectra using standard harmonic frequency analysis and spectra extracted from ab initio molecular dynamics simulations at different levels of theory. We conclude that the direct, quantitative comparison of calculated and measured IR spectra might become a novel, sophisticated approach to benchmark quantum-chemical methods. In the present benchmark, simulated spectra based on ab initio molecular dynamics show in general better agreement with the experiment than static calculations.

Multiresolution continuous wavelet transform for studying coupled solute-solvent vibrations via ab initio molecular dynamics

Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.

Toward a First-Principles Evaluation of Transport Mechanisms in Molecular Wires

Understanding charge transport through molecular wires is important for nanoscale electronics and biochemistry. Our goal is to establish a simple first-principles protocol for predicting the charge transport mechanism in such wires, in particular the crossover from coherent tunneling for short wires to incoherent hopping for longer wires. This protocol is based on a combination of density functional theory with a polarizable continuum model introduced by Kaupp et al. for mixed-valence molecules, which we had previously found to work well for length-dependent charge delocalization in such systems. We combine this protocol with a new charge delocalization measure tailored for molecular wires, and we show that it can predict the tunneling-to-hopping transition length with a maximum error of one subunit in five sets of molecular wires studied experimentally in molecular junctions at room temperature. This suggests that the protocol is also well suited for estimating the extent of hopping sites as relevant, for example, for the intermediate tunneling-hopping regime in DNA.

Computing Bulk Phase Resonance Raman Spectra from ab Initio Molecular Dynamics and Real-Time TDDFT

We present our novel approach for computing resonance Raman (RR) spectra of periodic bulk phase systems from ab initio molecular dynamics, including solvent influence and some anharmonic effects. Based on real-time time-dependent density functional theory, we obtain the RR spectra for all laser wavelengths in one pass. We compute the RR spectrum of uracil in aqueous solution, which is in good agreement with experiment. This is the first simulation of a bulk phase RR spectrum.

Uncertainty quantification for quantum chemical models of complex reaction networks

For the quantitative understanding of complex chemical reaction mechanisms, it is, in general, necessary to accurately determine the corresponding free energy surface and to solve the resulting continuous-time reaction rate equations for a continuous state space. For a general (complex) reaction network, it is computationally hard to fulfill these two requirements. However, it is possible to approximately address these challenges in a physically consistent way. On the one hand, it may be sufficient to consider approximate free energies if a reliable uncertainty measure can be provided. On the other hand, a highly resolved time evolution may not be necessary to still determine quantitative fluxes in a reaction network if one is interested in specific time scales. In this paper, we present discrete-time kinetic simulations in discrete state space taking free energy uncertainties into account. The method builds upon thermo-chemical data obtained from electronic structure calculations in a condensed-phase model. Our kinetic approach supports the analysis of general reaction networks spanning multiple time scales, which is here demonstrated for the example of the formose reaction. An important application of our approach is the detection of regions in a reaction network which require further investigation, given the uncertainties introduced by both approximate electronic structure methods and kinetic models. Such cases can then be studied in greater detail with more sophisticated first-principles calculations and kinetic simulations.

The parameter uncertainty inflation fallacy

Statistical estimation of the prediction uncertainty of physical models is typically hindered by the inadequacy of these models due to various approximations they are built upon. The prediction errors caused by model inadequacy can be handled either by correcting the model's results or by adapting the model's parameter uncertainty to generate prediction uncertainties representative, in a way to be defined, of model inadequacy errors. The main advantage of the latter approach (thereafter called PUI, for Parameter Uncertainty Inflation) is its transferability to the prediction of other quantities of interest based on the same parameters. A critical review of implementations of PUI in several areas of computational chemistry shows that it is biased, in the sense that it does not produce prediction uncertainty bands conforming to model inadequacy errors.

A 4 K FT-ICR cell for infrared ion spectroscopy

We present the design of the newly constructed cryogenic Fourier-transform ion cyclotron resonance (FT-ICR) ion trap for infrared ion spectroscopy. Trapped ions are collisionally cooled by the pulsed introduction of buffer gas into the cell. Using different buffer gases and cell temperatures, we record action spectra of weakly bound neutral gas-analyte complexes with an IR laser source. We show for the first time that ion-He complexes can be observed in an ICR cell at temperatures around 4 K. We compare the experimental vibrational spectra of Ag(PPh₃)₂⁺ obtained by tagging with different neutral gases: He, Ne, Ar, H₂, and N₂ to computed vibrational spectra. Furthermore, the conditions necessary for the formation of neutral tags within an ICR ion trap are studied.

Calculation of Ligand Dissociation Energies in Large Transition-Metal Complexes

The accurate calculation of ligand dissociation (or equivalently, ligand binding) energies is crucial for computational coordination chemistry. Despite its importance, obtaining accurate ab initio reference data is difficult, and density-functional methods of uncertain reliability are chosen for feasibility reasons. Here, we consider advanced coupled-cluster and multiconfigurational approaches to reinvestigate our WCCR10 set of 10 gas-phase ligand dissociation energies [ J. Chem. Theory Comput. 2014, 10, 3092]. We assess the potential multiconfigurational character of all molecules involved in these reactions with a multireference diagnostic [ Mol. Phys. 2017, 115, 2110] in order to determine where single-reference coupled-cluster approaches can be applied. For some reactions of the WCCR10 set, large deviations of density-functional results including semiclassical dispersion corrections from experimental reference data had been observed. This puzzling observation deserves special attention here, and we tackle the issue (i) by comparing to ab initio data that comprise dispersion effects on a rigorous first-principles footing and (ii) by a comparison of density-functional approaches that model dispersion interactions in various ways. For two reactions, species exhibiting nonnegligible static electron correlation were identified. These two reactions represent hard problems for electronic structure methods and also for multireference perturbation theories. However, most of the ligand dissociation reactions in WCCR10 do not exhibit static electron correlation effects, and hence, we may choose standard single-reference coupled-cluster approaches to compare with density-functional methods. For WCCR10, the Minnesota M06-L functional yielded the smallest mean absolute deviation of 13.2 kJ mol^-1 out of all density functionals considered (PBE, BP86, BLYP, TPSS, M06-L, PBE0, B3LYP, TPSSh, and M06-2X) without additional dispersion corrections in comparison to the coupled-cluster results, and the PBE0-D3 functional produced the overall smallest mean absolute deviation of 4.3 kJ mol^-1. The agreement of density-functional results with coupled-cluster data increases significantly upon inclusion of any type of dispersion correction. It is important to emphasize that different density-functional schemes available for this purpose perform equally well. The coupled-cluster dissociation energies, however, deviate from experimental results on average by 30.3 kJ mol^-1. Possible reasons for these deviations are discussed.

Complete analytic anharmonic hyper-Raman scattering spectra

We present the first computational treatment of the complete second-order vibrational perturbation theory applied to hyper-Raman scattering spectroscopy. The required molecular properties are calculated in a fully analytic manner using a recently developed program [Ringholm, Jonsson and Ruud, J. Comp. Chem., 2014, 35, 622] that utilizes recursive routines. For some of the properties, these calculations are the first analytic calculations of their kind at their respective levels of theory. We apply this approach to the calculation of the hyper-Raman spectra of methane, ethane and ethylene and compare these to available experimental data. We show that the anharmonic corrections have a larger effect on the vibrational frequencies than on the spectral intensities, but that the inclusion of combination and overtone bands in the anharmonic treatment can improve the agreement with the experimental data, although the quality of available experimental data limits a detailed comparison.

Anharmonic vibrational analysis of s-trans and s-cis conformers of acryloyl fluoride using numerical-analytic Van Vleck operator perturbation theory

A new gas-phase infrared (IR) spectrum of acryloyl fluoride (ACRF, CH₂CHCFO) with a resolution of 0.1cm^-1 in the range 4000-450cm^-1 was measured. Theoretical ab initio molecular structures, full quartic potential energy surfaces (PES), and cubic surfaces of dipole moments and polarizability tensor components (electro-optical properties, EOP) of the s-trans and s-cis conformers of the ACRF were calculated by the second-order Møller-Plesset electronic perturbation theory with a correlation consistent Dunning triple-ζ basis set. The numerical-analytic implementation of the second-order operator canonical Van Vleck perturbation theory was employed for predicting anharmonic IR and Raman scattering (RS) spectra of ACRF. To improve the anharmonic predictions, harmonic frequencies were replaced by their counterparts evaluated with the higher-level CCSD(T)/cc-pVTZ model, to form a "hybrid" PES. The original operator representation of the Hamiltonian is analytically reduced to a quasi-diagonal form, integrated in the harmonic oscillator basis and diagonalized to account for strong resonance couplings. Double canonical transformations of EOP expansions enabled prediction of integral intensities of both fundamental and multi-quanta transitions in IR/RS spectra. Enhanced band shape analysis reinforced the assignments. A thorough interpretation of the new IR experimental spectra and existing matrix-isolation literature data for the mixture of two conformers of ACRF was accomplished, and a number of assignments clarified.

Actinide embedded nearly planar gold superatoms: structural properties and applications in surface-enhanced Raman scattering (SERS)

Actinide embedded in a gold ring and applications in surface-enhanced Raman scattering (SERS).

Femtosecond infrared spectroscopy reveals the primary events of the ferrioxalate actinometer

Ultrafast vibrational spectroscopy unravels the primary processes of Hatchard and Parker's ferrioxalate actinometer – a widely used tool in the photochemistry community.

Probing Intermolecular Vibrational Symmetry Breaking in Self-Assembled Monolayers with Ultrahigh Vacuum Tip-Enhanced Raman Spectroscopy

Ultrahigh vacuum tip-enhanced Raman spectroscopy (UHV-TERS) combines the atomic-scale imaging capability of scanning probe microscopy with the single-molecule chemical sensitivity and structural specificity of surface-enhanced Raman spectroscopy. Here, we use these techniques in combination with theory to reveal insights into the influence of intermolecular interactions on the vibrational spectra of a N-N'-bis(2,6-diisopropylphenyl)-perylene-3,4:9,10-bis(dicarboximide) (PDI) self-assembled monolayer adsorbed on single-crystal Ag substrates at room temperature. In particular, we have revealed the lifting of a vibrational degeneracy of a mode of PDI on Ag(111) and Ag(100) surfaces, with the most strongly perturbed mode being that associated with the largest vibrational amplitude on the periphery of the molecule. This work demonstrates that UHV-TERS enables direct measurement of molecule-molecule interaction at nanoscale. We anticipate that this information will advance the fundamental understanding of the most important effect of intermolecular interactions on the vibrational modes of surface-bound molecules.

SERRS and absorption spectra of pyridine on Au m Ag n (m + n = 6) bimetallic nanoclusters: substrate composition and applied electric field effects

Surface-enhanced Raman scattering (SERS) and absorption spectra of the pyridine molecule adsorbed on Au _m Ag _n (m + n = 6) bimetallic clusters are theoretically investigated by time-dependent density functional theory. The contributions of static chemical enhancement to the ground-state system are analyzed, and the static Raman intensity of Py-Au _m Ag _n complexes are enhanced by an order of 10. A method of visualization on charge transfer is used to distinguish the contributions of charge-transfer enhancement and electromagnetic enhancement. The intensity of surface-enhanced resonance Raman scattering (SERRS) spectroscopy of Py-Au _m Ag _n is strongly enhanced by an order of 10³-10⁵, compared to the static Raman intensity of pyridine. The influence of the static external electric field on SERS is investigated by calculating the optical properties of the Py-Au₃Ag₃ complex. The intensity of SERRS spectra and normal Raman spectra can be significantly enhanced by the positive electric fields, and the intensities of specific Raman vibrational modes could be selectively enhanced or weakened by tuning the direction and strength of the static electric field applied on Py-Au₃Ag₃.

Efficient and Accurate Born-Oppenheimer Molecular Dynamics for Large Molecular Systems

An efficient scheme for the calculation of Born-Oppenheimer molecular dynamics (BOMD) simulations is introduced. It combines the corrected small basis set Hartree-Fock (HF-3c) method by Sure and Grimme [J. Comput. Chem. 2013, 43, 1672], extended Lagrangian BOMD (XL-BOMD) by Niklasson et al. [J. Chem. Phys. 2009, 130, 214109], and the calculation of the two electron integrals on graphics processing units (GPUs) [J. Chem. Phys. 2013, 138, 134114; J. Chem. Theory Comput. 2015, 11, 918]. To explore the parallel performance of our strong scaling implementation of the method, we present timings and extract, as its validation and first illustrative application, high-quality vibrational spectra from simulated trajectories of β-carotene, paclitaxel, and liquid water (up to 500 atoms). We conclude that the presented BOMD scheme may be used as a cost-efficient and reliable tool for computing vibrational spectra and thermodynamics of large molecular systems including explicit solvent molecules containing 500 atoms and more. Simulating 50 ps of maitotoxin (nearly 500 atoms) employing time steps of 0.5 fs requires ∼3 weeks on 12 CPUs (Intel Xeon E5 2620 v3) with 24 GPUs (AMD FirePro 3D W8100).

Computing Bulk Phase Raman Optical Activity Spectra from ab initio Molecular Dynamics Simulations

We present our novel methodology for computing Raman optical activity (ROA) spectra of liquid systems from ab initio molecular dynamics (AIMD) simulations. The method is built upon the recent developments to obtain magnetic dipole moments from AIMD and to integrate molecular properties by using radical Voronoi tessellation. These techniques are used to calculate optical activity tensors for large and complex periodic bulk phase systems. Only AIMD simulations are required as input, and no time-consuming perturbation theory is involved. The approach relies only on the total electron density in each time step and can readily be combined with a wide range of electronic structure methods. To the best of our knowledge, these are the first computed ROA spectra for a periodic bulk phase system. As an example, the experimental ROA spectrum of liquid (R)-propylene oxide is reproduced very well.

Reliable Estimation of Prediction Uncertainty for Physicochemical Property Models

One of the major challenges in computational science is to determine the uncertainty of a virtual measurement, that is the prediction of an observable based on calculations. As highly accurate first-principles calculations are in general unfeasible for most physical systems, one usually resorts to parameteric property models of observables, which require calibration by incorporating reference data. The resulting predictions and their uncertainties are sensitive to systematic errors such as inconsistent reference data, parametric model assumptions, or inadequate computational methods. Here, we discuss the calibration of property models in the light of bootstrapping, a sampling method that can be employed for identifying systematic errors and for reliable estimation of the prediction uncertainty. We apply bootstrapping to assess a linear property model linking the ⁵⁷Fe Mössbauer isomer shift to the contact electron density at the iron nucleus for a diverse set of 44 molecular iron compounds. The contact electron density is calculated with 12 density functionals across Jacob's ladder (PWLDA, BP86, BLYP, PW91, PBE, M06-L, TPSS, B3LYP, B3PW91, PBE0, M06, TPSSh). We provide systematic-error diagnostics and reliable, locally resolved uncertainties for isomer-shift predictions. Pure and hybrid density functionals yield average prediction uncertainties of 0.06-0.08 mm s^-1 and 0.04-0.05 mm s^-1, respectively, the latter being close to the average experimental uncertainty of 0.02 mm s^-1. Furthermore, we show that both model parameters and prediction uncertainty depend significantly on the composition and number of reference data points. Accordingly, we suggest that rankings of density functionals based on performance measures (e.g., the squared coefficient of correlation, r², or the root-mean-square error, RMSE) should not be inferred from a single data set. This study presents the first statistically rigorous calibration analysis for theoretical Mössbauer spectroscopy, which is of general applicability for physicochemical property models and not restricted to isomer-shift predictions. We provide the statistically meaningful reference data set MIS39 and a new calibration of the isomer shift based on the PBE0 functional.

Computational study of Th(4+) and Np(4+) hydration and hydrolysis of Th(4+) from first principles

The aqueous solvation of Th and Np in the IV oxidation state was examined using cluster models generated by Monte Carlo simulations and density functional theory embedded within the COSMO continuum model to approximate the effect of bulk water. Our results suggest that the coordination number (CN) for both Th(IV) and NP(IV) should be 9, in accordance to some experimental and theoretical results from the literature. The structural values for average oxygen-metal distances are within 0.01 Å compared to experimental data, and also within the experimental error. The calculated ΔG _Sol⁰ are in very good agreement with experimental reported values, with deviations at CN = 9 lower than 1% for both Th(IV) and Np(IV). The hydrolysis constants are also in very good agreement with experimental values. Finally, this [corrected] methodology has the advantage of using a GGA functional (BP86) that not only makes the calculations more affordable computationally than hybrid functional or ab initio molecular dynamics simulations (Car-Parrinello) calculations, but also opens the perspective to use resolution of identity (RI) calculations for more extended systems.

On the benefits of localized modes in anharmonic vibrational calculations for small molecules

Anharmonic vibrational calculations can already be computationally demanding for relatively small molecules. The main bottlenecks lie in the construction of the potential energy surface and in the size of the excitation space in the vibrational configuration interaction (VCI) calculations. To address these challenges, we use localized-mode coordinates to construct potential energy surfaces and perform vibrational self-consistent field and L-VCI calculations [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365 (2014)] for all vibrational modes of two prototypical test cases, the ethene and furan molecules. We find that the mutual coupling between modes is reduced when switching from normal-mode coordinates to localized-mode coordinates. When using such localized-mode coordinates, we observe a faster convergence of the n-mode expansion of the potential energy surface. This makes it possible to neglect higher-order contributions in the n-mode expansion of the potential energy surface or to approximate higher-order contributions in hybrid potential energy surfaces, which reduced the computational effort for the construction of the anharmonic potential energy surface significantly. Moreover, we find that when using localized-mode coordinates, the convergence with respect to the VCI excitation space proceeds more smoothly and that the error at low orders is reduced significantly. This makes it possible to devise low-cost models for obtaining a first approximation of anharmonic corrections. This demonstrates that the use of localized-mode coordinates can be beneficial already in anharmonic vibrational calculations of small molecules and provides a possible avenue for enabling such accurate calculations also for larger molecules.

Theoretic Study on Dispersion Mechanism of Boron Nitride Nanotubes by Polynucleotides

Due to the unique electrical and mechanical properties of boron nitride nanotubes (BNNT), BNNT has been a promising material for many potential applications, especially in biomedical field. Understanding the dispersion of BNNT in aqueous solution by biomolecules is essential for its use in biomedical applications. In this study, BNNT wrapped by polynucleotides in aqueous solution was investigated by molecular dynamics (MD) simulations. Our results demonstrated that the BNNT wrapped by polynucleotides could greatly hinder the aggregation of BNNTs and improve the dispersion of BNNTs in aqueous solution. Dispersion of BNNTs with the assistance of polynucleotides is greatly affected by the wrapping manner of polynucleotides on BNNT, which mainly depends on two factors: the type of polynucleotides and the radius of BNNT. The interaction between polynucleotides and BNNT(9, 9) is larger than that between polynucleotides and BNNT(5, 5), which leads to the fact that dispersion of BNNT(9, 9) is better than that of BNNT(5, 5) with the assistance of polynucleotides in aqueous solution. Our study revealed the molecular-level dispersion mechanism of BNNT with the assistance of polynucleotides in aqueous solution. It shades a light on the understanding of dispersion of single wall nanotubes by biomolecules.

Bisboronic Acids for Selective, Physiologically Relevant Direct Glucose Sensing with Surface-Enhanced Raman Spectroscopy

This paper demonstrates the direct sensing of glucose at physiologically relevant concentrations with surface-enhanced Raman spectroscopy (SERS) on gold film-over-nanosphere (AuFON) substrates functionalized with bisboronic acid receptors. The combination of selectivity in the bisboronic acid receptor and spectral resolution in the SERS data allow the sensors to resolve glucose in high backgrounds of fructose and, in combination with multivariate statistical analysis, detect glucose accurately in the 1-10 mM range. Computational modeling supports assignments of the normal modes and vibrational frequencies for the monoboronic acid base of our bisboronic acids, glucose and fructose. These results are promising for the use of bisboronic acids as receptors in SERS-based in vivo glucose monitoring sensors.

Exploring the relationship between vibrational mode locality and coupling using constrained optimization

The use of alternate coordinate systems as a means to improve the efficiency and accuracy of anharmonic vibrational structure analysis has seen renewed interest in recent years. While normal modes (which diagonalize the mass-weighted Hessian matrix) are a typical choice, the delocalized nature of this basis makes it less optimal when anharmonicity is in play. When a set of modes is not designed to treat anharmonicity, anharmonic effects will contribute to inter-mode coupling in an uncontrolled fashion. These effects can be mitigated by introducing locality, but this comes at its own cost of potentially large second-order coupling terms. Herein, a method is described which partially localizes vibrations to connect the fully delocalized and fully localized limits. This allows a balance between the treatment of harmonic and anharmonic coupling, which minimizes the error that arises from neglected coupling terms. Partially localized modes are investigated for a range of model systems including a tetramer of hydrogen fluoride, water dimer, ethene, diphenylethane, and stilbene. Generally, partial localization reaches ∼75% of maximal locality while introducing less than ∼30% of the harmonic coupling of the fully localized system. Furthermore, partial localization produces mode pairs that are spatially separated and thus weakly coupled to one another. It is likely that this property can be exploited in the creation of model Hamiltonians that omit the coupling parameters of the distant (and therefore uncoupled) pairs.

A polarizable QM/MM approach to the molecular dynamics of amide groups solvated in water

The infrared (IR) spectra of polypeptides are dominated by the so-called amide bands. Because they originate from the strongly polar and polarizable amide groups (AGs) making up the backbone, their spectral positions sensitively depend on the local electric fields. Aiming at accurate computations of these IR spectra by molecular dynamics (MD) simulations, which derive atomic forces from a hybrid quantum and molecular mechanics (QM/MM) Hamiltonian, here we consider the effects of solvation in bulk liquid water on the amide bands of the AG model compound N-methyl-acetamide (NMA). As QM approach to NMA we choose grid-based density functional theory (DFT). For the surrounding MM water, we develop, largely based on computations, a polarizable molecular mechanics (PMM) model potential called GP6P, which features six Gaussian electrostatic sources (one induced dipole, five static partial charge distributions) and, therefore, avoids spurious distortions of the DFT electron density in hybrid DFT/PMM simulations. Bulk liquid GP6P is shown to have favorable properties at the thermodynamic conditions of the parameterization and beyond. Lennard-Jones (LJ) parameters of the DFT fragment NMA are optimized by comparing radial distribution functions in the surrounding GP6P liquid with reference data obtained from a "first-principles" DFT-MD simulation. Finally, IR spectra of NMA in GP6P water are calculated from extended DFT/PMM-MD trajectories, in which the NMA is treated by three different DFT functionals (BP, BLYP, B3LYP). Method-specific frequency scaling factors are derived from DFT-MD simulations of isolated NMA. The DFT/PMM-MD simulations with GP6P and with the optimized LJ parameters then excellently predict the effects of aqueous solvation and deuteration observed in the IR spectra of NMA. As a result, the methods required to accurately compute such spectra by DFT/PMM-MD also for larger peptides in aqueous solution are now at hand.

Hydrogenase Enzymes and Their Synthetic Models: The Role of Metal Hydrides

Hydrogenase enzymes efficiently process H2 and protons at organometallic FeFe, NiFe, or Fe active sites. Synthetic modeling of the many H2ase states has provided insight into H2ase structure and mechanism, as well as afforded catalysts for the H2 energy vector. Particularly important are hydride-bearing states, with synthetic hydride analogues now known for each hydrogenase class. These hydrides are typically prepared by protonation of low-valent cores. Examples of FeFe and NiFe hydrides derived from H2 have also been prepared. Such chemistry is more developed than mimicry of the redox-inactive monoFe enzyme, although functional models of the latter are now emerging. Advances in physical and theoretical characterization of H2ase enzymes and synthetic models have proven key to the study of hydrides in particular, and will guide modeling efforts toward more robust and active species optimized for practical applications.

Nanoscale Chemical Imaging of a Dynamic Molecular Phase Boundary with Ultrahigh Vacuum Tip-Enhanced Raman Spectroscopy

Nanoscale chemical imaging of a dynamic molecular phase boundary has broad implications for a range of problems in catalysis, surface science, and molecular electronics. While scanning probe microscopy (SPM) is commonly used to study molecular phase boundaries, its information content can be severely compromised by surface diffusion, irregular packing, or three-dimensional adsorbate geometry. Here, we demonstrate the simultaneous chemical and structural analysis of N-N'-bis(2,6-diisopropylphenyl)-1,7-(4'-t-butylphenoxy)perylene-3,4:9,10-bis(dicarboximide) (PPDI) molecules by UHV tip-enhanced Raman spectroscopy. Both condensed and diffusing domains of PPDI coexist on Ag(100) at room temperature. Through comparison with time-dependent density functional theory simulations, we unravel the orientation of PPDI molecules at the dynamic molecular domain boundary with unprecedented ∼4 nm spatial resolution.

Optimizing Vibrational Coordinates To Modulate Intermode Coupling

The choice of coordinate system strongly affects the convergence properties of vibrational structure computations. Two methods for efficient generation of improved vibrational coordinates are presented and justified by analysis of a model anharmonic two-mode Hessian and numerical computations on polyatomic molecules. To produce optimal coordinates, metrics which quantify off-diagonal couplings over a grid of Hessian matrices are minimized through unitary rotations of the vibrational basis. The first proposed metric minimizes the total squared off-diagonal coupling, and the second minimizes the total squared change in off-diagonal coupling. In this procedure certain anharmonic modes tend to localize, for example X-H stretches. The proposed methods do not rely on prior fitting of the potential energy, vibrational structure computations, or localization metrics, so they are unique from previous vibrational coordinate generation algorithms and are generally applicable to polyatomic molecules. Fitting the potential to the approximate n-mode representation in the optimized bases for all-trans polyenes shows that off-diagonal anharmonic couplings are substantially reduced by the new choices of coordinate system. Convergence of vibrational energies is examined in detail for ethylene, and it is shown that coupling-optimized modes converge in vibrational configuration interaction computations to within 1 cm(-1) using only 3-mode couplings, where normal modes require 4-mode couplings for convergence. Comparison of the vibrational configuration interaction convergence with respect to excitation level for the two proposed metrics shows that minimization of the total off-diagonal coupling is most effective for low-cost vibrational structure computations.

Nonresonant chemical mechanism in surface-enhanced Raman scattering of pyridine on M@Au12 clusters

By employing density functional theory (DFT), this study presents a detailed analysis of nonresonant surface-enhanced Raman scattering (SERS) of pyridine on M@Au12 (M = V(-), Nb(-), Ta(-), Cr, Mo, W, Mn(+), Tc(+), and Re(+))-the stable 13-atom neutral and charged gold buckyball clusters. Changing the core atom in M@Au12 enabled us to modulate the direct chemical interactions between pyridine and the metal cluster. The results of our calculations indicate that the ground-state chemical enhancement does not increase as the binding interaction strengthens or the transfer charge increases between pyridine and the cluster. Instead, the magnitude of the chemical enhancement is governed, to a large extent, by the charged properties of the metal clusters. Pyridine on M@Au12 anion clusters exhibits strong chemical enhancement of a factor of about 10(2), but the equivalent increase for pyridine adsorbed on M@Au12 neutral and cation clusters is no more than 10. Polarizability and deformation density analyses clearly show that compared with the neutral and cation clusters, the anion clusters have more delocalized electrons and occupy higher energy levels in the pyridine-metal complex. Accordingly, they produce larger polarizability, leading to a stronger nonresonant enhancement effect.