Reliable structural, thermodynamic, and spectroscopic properties of organic molecules adsorbed on silicon surfaces from computational modeling: the case of glycine@Si(100)

First-Principles Calculations for Glycine Adsorption Dynamics and Surface-Enhanced Raman Spectroscopy on Diamond Surfaces

Based on first-principles calculations, the stability of three adsorption configurations of glycine on the (100) surface of diamonds was studied, leading to an investigation into the surface-enhanced Raman scattering (SERS) effect of the diamond substrate. The results showed that the carboxyl-terminated adsorption configuration (CAR) was the most stable and shortest interface distance compared to other configurations. This stability was primarily attributed to the formation of strong polar covalent bonds between the carboxyl O atoms and the surface C atoms of the (100) surface of diamonds. These results were further corroborated by first-principles molecular dynamics simulations. Within the temperature range of 300 to 500 K, the glycine molecules in the carboxyl-terminated adjacent-dimer phenyl-like (CAR) configuration exhibited only simple thermal vibrations with varying amplitudes. In contrast, the metastable ATO and carboxyl-terminated trans-dimer phenyl-like ring (CTR) configurations were observed to gradually transform into benzene-ring-like structures akin to the CAR configuration. After adsorption, the intensity of glycine's characteristic peaks increased substantially, accompanied by a blue shift phenomenon. Notably, the characteristic peaks related to the carboxyl and amino groups exhibited the highest enhancement amplitude, exceeding 200 times, with an average enhancement amplitude exceeding 50 times. The diamond substrate, with its excellent adsorption properties and strong surface Raman spectroscopy characteristics, represents a highly promising candidate in the field of biomedicine.

Describing the Disulfide Bond: From the Density Functional Theory and Back through the "Lego Brick" Approach

Selected molecular species containing the disulfide bond, RSSR, have been considered, these ranging from hydrogen disulfide, H2S2 (R = H), to diphenyl disulfide with R = C6H5. The aim of this work is two-fold: (i) to investigate different computational approaches in order to derive accurate equilibrium structures at an affordable cost, (ii) to employ the results from the first goal in order to benchmark cheaper methodologies rooted in the density functional theory. Among the strategies used for the accurate geometrical determinations, the semiexperimental approach has been exploited in combination with a reduced-dimensionality VPT2 model, without however obtaining satisfactory results. Instead, the so-called "Lego brick" approach turned out to be very effective despite the flexibility of the systems investigated. Concerning the second target of this work, the focus was mainly on the S-S bond and the structural parameters related to it. Among those tested, PBE0(-D3BJ), M06-2X(-D3) and DSD-PBEP86-D3BJ have been found to be the best-performing functionals.

Tracing the Primordial Chemical Life of Glycine: A Review from Quantum Chemical Simulations

Glycine (Gly), NH2CH2COOH, is the simplest amino acid. Although it has not been directly detected in the interstellar gas-phase medium, it has been identified in comets and meteorites, and its synthesis in these environments has been simulated in terrestrial laboratory experiments. Likewise, condensation of Gly to form peptides in scenarios resembling those present in a primordial Earth has been demonstrated experimentally. Thus, Gly is a paradigmatic system for biomolecular building blocks to investigate how they can be synthesized in astrophysical environments, transported and delivered by fragments of asteroids (meteorites, once they land on Earth) and comets (interplanetary dust particles that land on Earth) to the primitive Earth, and there react to form biopolymers as a step towards the emergence of life. Quantum chemical investigations addressing these Gly-related events have been performed, providing fundamental atomic-scale information and quantitative energetic data. However, they are spread in the literature and difficult to harmonize in a consistent way due to different computational chemistry methodologies and model systems. This review aims to collect the work done so far to characterize, at a quantum mechanical level, the chemical life of Gly, i.e., from its synthesis in the interstellar medium up to its polymerization on Earth.

Computational vibrational spectroscopy of molecule-surface interactions: what is still difficult and what can be done about it

Interactions of molecules with solid surfaces are responsible for key functionalities for a range of currently actively pursued technologies, including heterogeneous catalysis for synthesis or decomposition of molecules, sensitization, surface...

Formation of Germa-ketenimine on the Ge(100) Surface by Adsorption of tert-Butyl Isocyanide

Reactions of the (100) surfaces of Ge and Si with organic molecules have been generally understood within the concept of "dimers" formed by the 2 × 1 surface reconstruction. In this work, the adsorption of tert-butyl isocyanide on the Ge(100)-2 × 1 surface at large exposures is investigated under ultrahigh vacuum conditions. A combination of infrared spectroscopy, X-ray photoelectron spectroscopy, and temperature-programmed desorption experiments along with dispersion-corrected density functional theory calculations is used to determine the surface products. Upon adsorption of a dense monolayer of tert-butyl isocyanide, a product whose structure resembles a germa-ketenimine (N=C=Ge) with σ donation toward and π back-donation from the Ge(100) surface appears. Formation of this structure involves divalent-type surface Ge atoms that arise from cleavage of the Ge(100)-2 × 1 surface dimers. Our results reveal an unprecedented class of reactions of organic molecules at the Ge(100) surface.

Reliable vibrational wavenumbers for C=O and N-H stretchings of isolated and hydrogen-bonded nucleic acid bases

The accurate prediction of vibrational wavenumbers for functional groups involved in hydrogen-bonded bridges remains an important challenge for computational spectroscopy. For the specific case of the C=O and N-H stretching modes of nucleobases and their oligomers, the paucity of experimental reference values needs to be compensated by reliable computational data, which require the use of approaches going beyond the standard harmonic oscillator model. Test computations performed for model systems (formamide, acetamide and their cyclic homodimers) in the framework of the second order vibrational perturbation theory (VPT2) confirmed that anharmonic corrections can be safely computed by global hybrid (GHF) or double hybrid (DHF) functionals, whereas the harmonic part is particularly challenging. As a matter of fact, GHFs perform quite poorly and even DHFs, while fully satisfactory for C=O stretchings, face unexpected difficulties when dealing with N-H stretchings. On these grounds, a linear regression for N-H stretchings has been obtained and validated for the heterodimers formed by 4-aminopyrimidine with 6-methyl-4-pyrimidinone (4APM-M4PMN) and by uracil with water. In view of the good performance of this computational model, we have built a training set of B2PLYP-D3/maug-cc-pVTZ harmonic wavenumbers (including linear regression scaling for N-H) for six-different uracil dimers and a validation set including 4APM-M4PMN, one of the most stable hydrogen-bonded adenine homodimers, as well as the adenine-uracil, adenine-thymine, guanine-cytosine and adenine-4-thiouracil heterodimers. Because of the unfavourable scaling of DHF harmonic wavenumbers with the dimensions of the investigated systems, we have optimized a linear regression of B3LYP-D3/N07D harmonic wavenumbers for the training set, which has been next checked against the validation set. This relatively cheap model, which shows very good agreement with experimental data (average errors of about 10 cm(-1)), paves the route toward a reliable analysis of spectroscopic signatures for larger polynucleotides.

The Virtual Multifrequency Spectrometer: a new paradigm for spectroscopy

On going developments of hardware and software are changing computational spectroscopy from a strongly specialized research area to a general tool in the inventory of most researchers. Increased interactions between experimentally-oriented users and theoretically-oriented developers of new methods and models would result in more robust, flexible and reliable tools and studies for the systems of increasing complexity, which are of current scientific and technological interest. This is the philosophy behind this review, which presents the development of a so-called virtual multi-frequency spectrometer (VMS) including state-of-the-art approaches in a user-friendly frame. The current status of the VMS tool will be illustrated by a number of case studies with special reference to infrared and UV-vis regions of the electro-magnetic spectrum including also chiral spectroscopies. Only the basic theoretical background will be provided avoiding explicit equations as far as possible, and pointing out the most recent advancements beyond the standard rigid-rotor harmonic-oscillator model coupled to vertical electronic excitation energies.

Anharmonic Effects on Vibrational Spectra Intensities: Infrared, Raman, Vibrational Circular Dichroism, and Raman Optical Activity

The aim of this paper is 2-fold. First, we want to report the extension of our virtual multifrequency spectrometer (VMS) to anharmonic intensities for Raman optical activity (ROA) with the full inclusion of first- and second-order resonances for both frequencies and intensities in the framework of the generalized second-order vibrational perturbation theory (GVPT2) for all kinds of vibrational spectroscopies. Then, from a more general point of view, we want to present and validate the performance of VMS for the parallel analysis of different vibrational spectra for medium-sized molecules (IR, Raman, VCD, ROA) including both mechanical and electric/magnetic anharmonicity. For the well-known methyloxirane benchmark, careful selection of density functional, basis set, and resonance thresholds permitted us to reach qualitative and quantitative agreement between experimental and computed band positions and shapes. Next, the whole series of halogenated azetidinones is analyzed, showing that it is now possible to interpret different spectra in terms of mass, electronegativity, polarizability, and hindrance variation between closely related substituents, chiral spectroscopies being particular effective in this connection.

Quantum Chemistry Meets Spectroscopy for Astrochemistry: Increasing Complexity toward Prebiotic Molecules

For many years, scientists suspected that the interstellar medium was too hostile for organic species and that only a few simple molecules could be formed under such extreme conditions. However, the detection of approximately 180 molecules in interstellar or circumstellar environments in recent decades has changed this view dramatically. A rich chemistry has emerged, and relatively complex molecules such as C60 and C70 are formed. Recently, researchers have also detected complex organic and potentially prebiotic molecules, such as amino acids, in meteorites and in other space environments. Those discoveries have further stimulated the debate on the origin of the building blocks of life in the universe. Many efforts continue to focus on the physical, chemical, and astrophysical processes by which prebiotic molecules can be formed in the interstellar dust and dispersed to Earth or to other planets.Spectroscopic techniques, which are widely used to infer information about molecular structure and dynamics, play a crucial role in the investigation of planetary atmosphere and the interstellar medium. Increasingly these astrochemical investigations are assisted by quantum-mechanical calculations of structures as well as spectroscopic and thermodynamic properties, such as transition frequencies and reaction enthalpies, to guide and support observations, line assignments, and data analysis in these new and chemically complicated situations. However, it has proved challenging to extend accurate quantum-chemical computational approaches to larger systems because of the unfavorable scaling with the number of degrees of freedom (both electronic and nuclear).In this Account, we show that it is now possible to compute physicochemical properties of building blocks of biomolecules with an accuracy rivaling that of the most sophisticated experimental techniques, and we summarize specific contributions from our groups. As a test case, we present the underlying computational machinery through the investigation of oxirane. We describe how we determine the molecular structure and then how we characterize the rotational and IR spectra, the most important issues for a correct theoretical description and a proper comparison with experiment. Next, we analyze the spectroscopic properties of representative building blocks of DNA bases (uracil and pyrimidine) and of proteins (glycine and glycine dipeptide analogue).Solvation, surface chemistry (dust fraction, adsorption, desorption), and inter- and intramolecular interactions, such as self-organization and self-interaction, are important molecular processes for understanding astrochemistry. Using the specific cases of uracil dimers and glycine adsorbed on silicon grains, we also illustrate approaches in which we treat different regions, interactions, or effects at different levels of sophistication.

A new Gaussian MCTDH program: implementation and validation on the levels of the water and glycine molecules

We report the main features of a new general implementation of the Gaussian Multi-Configuration Time-Dependent Hartree model. The code allows effective computations of time-dependent phenomena, including calculation of vibronic spectra (in one or more electronic states), relative state populations, etc. Moreover, by expressing the Dirac-Frenkel variational principle in terms of an effective Hamiltonian, we are able to provide a new reliable estimate of the representation error. After validating the code on simple one-dimensional systems, we analyze the harmonic and anharmonic vibrational spectra of water and glycine showing that reliable and converged energy levels can be obtained with reasonable computing resources. The data obtained on water and glycine are compared with results of previous calculations using the vibrational second-order perturbation theory method. Additional features and perspectives are also shortly discussed.

Dispersion corrected DFT approaches for anharmonic vibrational frequency calculations: nucleobases and their dimers

Computational spectroscopy techniques have become in the last few years an effective means to analyze and assign infrared (IR) spectra of molecular systems of increasing dimensions and in different environments. However, transition from compilation of harmonic data to fully anharmonic simulations of spectra is still underway. The most promising results for large systems have been obtained, in our opinion, by perturbative vibrational approaches based on potential energy surfaces computed by hybrid (especially B3LYP) density functionals and medium size (e.g. SNSD) basis sets. In this framework, we are actively developing a comprehensive and robust computational protocol aimed at quantitative reproduction of the spectra of nucleic acid base complexes and their adsorption on solid supports (organic/inorganic). In this contribution we report the essential results of the first step devoted to isolated monomers and dimers. It is well known that in order to model the vibrational spectra of weakly bound molecular complexes dispersion interactions should be taken into proper account. In this work we have chosen two popular and inexpensive approaches to model dispersion interactions, namely the semi-empirical dispersion correction (D3) and pseudopotential based (DCP) methodologies both in conjunction with the B3LYP functional. These have been used for simulating fully anharmonic IR spectra of nucleobases and their dimers through generalized second order vibrational perturbation theory (GVPT2). We have studied, in particular, isolated adenine, hypoxanthine, uracil, thymine and cytosine, the hydrogen-bonded and stacked adenine and uracil dimers, and the stacked adenine-naphthalene heterodimer. Anharmonic frequencies are compared with standard B3LYP results and experimental findings, while the computed interaction energies and structures of complexes are compared to the best available theoretical estimates.

Fully anharmonic IR and Raman spectra of medium-size molecular systems: accuracy and interpretation

Computation of full infrared (IR) and Raman spectra (including absolute intensities and transition energies) for medium- and large-sized molecular systems beyond the harmonic approximation is one of the most interesting challenges of contemporary computational chemistry. Contrary to common beliefs, low-order perturbation theory is able to deliver results of high accuracy (actually often better than those issuing from current direct dynamics approaches) provided that anharmonic resonances are properly managed. This perspective sketches the recent developments in our research group toward the development of a robust and user-friendly virtual spectrometer rooted in second-order vibrational perturbation theory (VPT2) and usable also by non-specialists essentially as a black-box procedure. Several examples are explicitly worked out in order to illustrate the features of our computational tool together with the most important ongoing developments.

A robust and effective time-independent route to the calculation of Resonance Raman spectra of large molecules in condensed phases with the inclusion of Duschinsky, Herzberg-Teller, anharmonic, and environmental effects

We present an effective time-independent implementation to model vibrational resonance Raman (RR) spectra of medium-large molecular systems with the inclusion of Franck-Condon (FC) and Herzberg-Teller (HT) effects and a full account of the possible differences between the harmonic potential energy surfaces of the ground and resonant electronic states. Thanks to a number of algorithmic improvements and very effective parallelization, the full computations of fundamentals, overtones, and combination bands can be routinely performed for large systems possibly involving more than two electronic states. In order to improve the accuracy of the results, an effective inclusion of the leading anharmonic effects is also possible, together with environmental contributions under different solvation regimes. Reduced-dimensionality approaches can further enlarge the range of applications of this new tool. Applications to imidazole, pyrene, and chlorophyll a1 in solution are reported, as well as comparisons with available experimental data.

Anharmonic theoretical simulations of infrared spectra of halogenated organic compounds

The recent implementation of the computation of infrared (IR) intensities beyond the double-harmonic approximation [J. Bloino and V. Barone, J. Chem. Phys. 136, 124108 (2012)] paved the route to routine calculations of infrared spectra for a wide set of molecular systems. Halogenated organic compounds represent an interesting class of molecules, from both an atmospheric and computational point of view, due to the peculiar chemical features related to the halogen atoms. In this work, we simulate the IR spectra of eight halogenated molecules (CH2F2, CHBrF2, CH2DBr, CF3Br, CH2CHF, CF2CFCl, cis-CHFCHBr, cis-CHFCHI), using two common hybrid and double-hybrid density functionals in conjunction with both double- and triple-ζ quality basis sets (SNSD and cc-pVTZ) as well as employing the coupled-cluster theory with basis sets of at least triple-ζ quality. Finally, we compare our results with available experimental spectra, with the aim of checking the accuracy and the performances of the computational approaches.

Accurate structure, thermodynamic and spectroscopic parameters from CC and CC/DFT schemes: the challenge of the conformational equilibrium in glycine

The structures, relative stabilities, and infrared spectra of the six low-energy conformers of glycine have been characterized using a state-of-the-art quantum-mechanical approach allowing the bond distances, conformational enthalpies and vibrational frequencies to be determined well within the chemical accuracy. Transition state structures governing interconversion among the different energy minima have also been characterized. In detail, the gas-phase thermodynamic properties (at 15 K and 410 K) of the glycine conformers considered have been obtained with a 1 kJ mol(-1) accuracy, and it has been shown that the employment of DFT geometries usually reduces such accuracy by at most 0.1 kJ mol(-1). Regarding molecular structures, the use of two different composite schemes allowed us to further confirm the suitability of a rather cost-effective approach and provide geometrical parameters with an overall accuracy better than 0.002 Å for distances and 1 degree for angles. Thanks to a hybrid CC/DFT approach, the infrared spectra of all conformers considered and of several deuterated isotopologues have been reproduced (when experimental data were available) or predicted with an accuracy of 10 cm(-1). Finally, the joint thermodynamic and spectroscopic investigation allowed us to shed some light on the possible observation of elusive conformers. On the whole, the high accuracy of the computational results allows us to draw a fully consistent interpretation of the available experimental data and to obtain a more complete characterization of the potential energy surface of glycine.

Theoretical investigation on the structural and thermodynamic properties of FeSe at high pressure and high temperature

A theoretical investigation on structural and thermodynamic properties of 11-type iron-based superconductor FeSe at high pressure and high temperature was performed by employing the first-principles method based on the density functional theory. Some structural parameters of FeSe in both tetragonal and hexagonal phases are reported. According to the fourth-order Birch-Murnaghan equation of states, the transition pressure P(t) of FeSe from the PbO-type phase to the NiAs-type phase was determined. The calculated results are found to be in good agreement with the available experimental data. Based on the quasi-harmonic Debye model, the pressure and temperature dependence of the thermodynamic properties for hexagonal phase FeSe were investigated. Our theoretical calculations suggest that the pressure and temperature have significant effects on the heat capacity, vibrational internal energy, vibrational entropy, vibrational Helmholtz free energy, thermal expansion coefficient and Debye temperature. Even though few theoretical reports on the structural properties of FeSe are found in the current literature, to our knowledge, this is a novel theoretical investigation on the structural and thermodynamic properties of FeSe at high temperature. We hope that the theoretical results reported here can give more insight into the structural and thermodynamic properties of other iron-based superconductors at high temperature.