Toward spectroscopic studies of biologically relevant systems: Vibrational spectrum of adenine as a test case for performances of long-range/dispersion corrected density functionals

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers

Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.

Comparative assessment of structure-property relationships of new Cu(II) complex in selected density functionals

In order to evaluate the structure-property relationships of Cu(II) complex by using DFT methods, the structure of the newly synthesized Cu(II) complex, [Cu(6-Brpic)2(bpy)], was investigated by XRD, FTIR, UV-Vis, and fluorescence spectroscopic methods. In addition, Hirshfeld surface and NBO analyses were fulfilled to identify possible interactions in the intermolecular and coordination environment. The five different DFT methods (HCTH, M06L, TPSSTPSS, B3LYP, and CAM-B3LYP levels), having four different functionalities (the GGA, meta-GGA, hybrid-GGA, and range-separated hybrid), were carried out so as to investigate the structure-property relationship, considering the geometric parameters (bond lengths and angles), vibrational frequencies, electronic absorption wavelengths, electronic transitions, and linear and nonlinear optical parameters. The R2 for structural and vibrational parameters, as well as MPD%, MAD, an optimal scaling factor (λ) and overall root mean square (RMS) deviation, were considered only at vibration frequencies. While it was determined that M06-L and TPSSTPSS levels gave the best results for the bond lengths and angles of the Cu(II) complex, the best results for vibrational frequencies were obtained in the HCTH method along with these methods. In NLO parameters, the static and dynamic first-order hyperpolarizability (<β(0;0,0)> and β(-ω;ω,0)/<β(-2ω;ω,ω)>) values, the largest values were obtained in the HCTH method (38.817 × 10-30 and 437.86 × 10-30/201.55 × 10-30 esu), whereas the smallest values were found to be in the CAM-B3LYP/TPSSTPSS levels (6.118 × 10-30 esu, 8.270 × 10-30/11.730 × 10-30 esu). By regarding the static γ (<γ(0;0,0,0)>) and dynamic (<γ(-ω;ω,0,0)> parameters, the largest values were calculated in the M06L (232.101 × 10-36) and HCTH (1711.52 × 10-36) methods and the smallest values were obtained in the CAM-B3LYP (43.281 × 10-36 and 60.844 × 10-36) method. In fact, it is obviously seen that the β and γ values obtained by the aforementioned DFT levels are many times higher than that of the standard molecule of urea. These results indicate that the Cu(II) complex may be used as a potential NLO material to evolve optoelectronic devices.

Simulated NIR spectra as sensitive markers of the structure and interactions in nucleobases

Near-infrared (near-IR; NIR) spectroscopy is continuously advancing in biophysical and biochemical fields of investigation. For instance, recent progresses in NIR hyperspectral imaging of biological systems may be noted. However, interpretation of NIR bands for biological samples is difficult and creates a considerable barrier in exploring the full potential of NIR spectroscopy in bioscience. For this reason, we carried out a systematic study of NIR spectra of adenine, cytosine, guanine, and thymine in polycrystalline state. Interpretation of NIR spectra of these nucleobases was supported by anharmonic vibrational analysis using Deperturbed Vibrational Second-Order Perturbation Theory (DVPT2). A number of molecular models of nucleobases was applied to study the effect of the inter-molecular interactions on the NIR spectra. The accuracy of simulated NIR spectra appears to depend on the intra-layer interactions; in contrast, the inter-layer interactions are less influential. The best results were achieved by combining the simulated spectra of monomers and dimers. It is of particular note that in-plane deformation bands are far more populated than out-of-plane ones and the importance of ring modes is relatively small. This trend is in contrast to that observed in mid-IR region. As shown, the local, short-range chemical neighborhood of nucleobase molecules influence their NIR spectra more considerably. This suggests that NIR spectra are more sensitive probe of the nucleobase pairing than mid-IR ones. The obtained results allow, for the first time, to construct a frequency correlation table for NIR spectra of purines and pyrimidines.

Vibrational investigation of nucleobases by means of divide and conquer semiclassical dynamics

In this work, we report a computational study of the vibrational features of four different nucleobases employing the divide-and-conquer semiclassical initial value representation molecular dynamics method. Calculations are performed on uracil, cytosine, thymine, and adenine. Results show that the overall accuracy with respect to experiments is within 20 wavenumbers, regardless of the dimensionality of the nucleobase. Vibrational estimates are accurate even in the complex case of cytosine, where two relevant conformers are taken into account. These results are promising in the perspective of future studies on more complex systems, such as nucleotides or nucleobase pairs.

Can Density Functional Theory Be Trusted for High-Order Electric Properties? The Case of Hydrogen-Bonded Complexes

This work reports on an extensive assessment of the performance of a wide palette of density functional approximations in predicting the (high-order) electric properties of hydrogen-bonded complexes. To this end, we compute the electronic and vibrational contributions to the electric polarizability and the first and second hyperpolarizabilities, using the CCSD(T)/aug-cc-pVTZ level of theory as reference. For all the studied properties, the average absolute errors below 20% can only be obtained using the CAM-B3LYP functional, while LC-BLYP and MN15 are shown to be only slightly less accurate (average absolute errors not exceeding 30%). Among Minnesota density functionals, i.e., M06, M06-2X, and MN15, we only recommend the latter one, which quite accurately predicts the electronic and vibrational (hyper)polarizabilities. We also analyze the optimal tuning of the range-separation parameter μ for the LC-BLYP functional, finding that this approach does not bring any systematic improvement in the predictions of electronic and vibrational (hyper)polarizabilities and the accuracy of computed properties is largely system-dependent. Finally, we report huge errors in predicting the vibrational second hyperpolarizability by ωB97X, M06, and M06-2X functionals. Based on the explicit evaluation of anharmonic terms contributing to the second hyperpolarizability, this failure is traced down to a poor determination of third- and fourth-order energy derivatives with respect to normal modes. These results reveal serious flaws of some density functional approximations and suggest caution in selecting the appropriate functional to calculate not only electronic and vibrational (hyper)polarizabilities but also other molecular properties that contain vibrational anharmonic contributions.

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.

Vibrational and quantum chemical investigation of cyclization of thiosemicarbazide group in 1-benzoyl-4-phenyl-3-thiosemicarbazide

1-Benzoyl-4-phenyl-3-thiosemicarbazide (H3bpt) was treated with acid - base in one sequence and base - acid in other sequence, both of which lead to ring formation of thiosemicarbazide group, giving N-phenyl-5-phenyl-1,3,4-thiadiazol-2-amine (Hppta) in the first case and 4,5-diphenyl-2,4-dihydro-1,2,4-triazole-3-thione (Hdptt) in the second case. The primary (H3bpt) as well as the resulting compounds (Hppta & Hdptt) has been characterized by elemental analyses, NMR, FTIR and Raman spectroscopic techniques. The quantum chemical calculations of the compounds are performed using DFT/B3LYP/6311G(d,p) method for geometry optimizations and also for prediction of the molecular properties. The cyclization is confirmed by disappearance of many bands belonging to the open chain subgroups of H3bpt such as; NH stretching, NH bending, CN stretching, NH puckering, CO stretching etc. The ring formation of 1-benzoyl-4-phenyl-3-thiosemicarbazide (H3bpt) has been further confirmed by the appearance of many bands belonging to the closed ring of thiosemicarbazide in the resulting compounds Hppta and Hdptt.

Quantum chemical and spectroscopic investigations of 3-methyladenine

FTIR, FT-Raman and UV-Vis spectra of 3-methyladenine have been recorded and investigated using quantum chemical calculations. The molecular geometry and vibrational spectra of 3-methyladenine in the ground state are computed by using HF and DFT methods with 6-311G(d,p) basis set. VSCF, CC-VSCF methods based on 2MR-QFF and PT2 (Barone method) have been utilized for computing anharmonic vibrational frequencies. These methods yield results that are in remarkable agreement with the experimental data. The magnitudes of coupling between pair of modes have been also computed. Vibrational modes are assigned with the help of visual inspection of atomic displacements. The electronic spectra, simulated at TD-B3LYP/6-311++G(d,p) level of theory, are compared to the experiment. The global quantities: electronic chemical potential, electrophilicity index, chemical hardness and softness based on HOMO and LUMO energy eigenvalues are also computed at B3LYP/6-311++G(d,p) level of theory.

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.

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.

Investigation of coordination properties of isolated adenine to copper metal: a systematic spectroscopic and DFT study

The coordination properties of copper with adenine have been studied by the analyzing the changes in Fourier Transform Infra-red (FTIR) and Raman spectra of adenine and adenine-copper complex. The geometry of adenine and adenine copper complex were optimized and theoretical Infra-red and Raman spectra of the optimized structures were calculated using Density Functional Theory (DFT). During synthesis of adenine-copper complex specific procedure was adopted to attach the Cu atom with particular N-atom of adenine (N9). The results of Raman and DFT confirmed the attachment. The Raman bands at 625, 330 and 230 cm(-1) of adenine-copper complex contain significant contribution of the vibrational motions of Cu metal coordinated to N9 and Cl atoms. The DFT calculations give additional vibrational modes containing the Cu, N9 and N9* atoms, which are not observed in FTIR and Raman spectra. The Raman, IR and DFT study confirm that Cu metal has good binding affinity to the isolated adenine base.

Theoretical studies on the dimerization of substituted paraphenylenediamine radical cations

Organic radical cations form dicationic dimers in solution, observed experimentally as diamagnetic species in temperature-dependent EPR and low temperature UV/Vis spectroscopy. Dimerization of paraphenylenediamine, N,N-dimethyl-paraphenylenediamine and 2,3,5,6-tetramethyl-paraphenylenediamine radical cation in ethanol/diethylether mixture was investigated theoretically according to geometry, energetics and UV/Vis spectroscopy. Density Functional Theory including dispersion correction describes stable dimers after geometry optimization with conductor-like screening model of solvation and inclusion of the counter-ion. Energy corrections were done on double-hybrid Density Functional Theory with perturbative second-order correlation (B2PLYP-D) including basis set superposition error (BSSE), and multireference Møller-Plesset second-order perturbation theory method (MRMP2) based on complete active space method (CASSCF(2,2)) single point calculation, respectively. All three dication π-dimers exhibit long multicenter π-bonds around 2.9±0.1Å with strongly interacting orbitals. Substitution with methyl groups does not influence the dimerization process substantially. Dispersion interaction and electrostatic attraction from counter-ion play an important role to stabilize the dication dimers in solution. Dispersion-corrected double hybrid functional B2PLYP-D and CASSCF(2,2) can describe the interaction energetics properly. Vertical excitations were computed with Tamm-Dancoff approximation for time-dependent Density Functional Theory (TDA-DFT) at the B3LYP level with the cc-pVTZ basis set including ethanol solvent molecules explicitly. A strong interaction of the counter-ion and the solvent ethanol with the monomeric species is observed, whereas in the dimers the strong interaction of both radical cation species is the dominating factor for the additional peak in UV/Vis spectra.

Environmental effects in computational spectroscopy: accuracy and interpretation

Spectroscopic techniques are valuable tools for understanding the structure and dynamics of complex systems, such as biomolecules or nanomaterials. Most of the current research is devoted to the development of new experimental techniques for improving the intrinsic resolution of different spectra. However, the subtle interplay of several different effects acting at different length and time scales still makes the interpretation and analysis of such spectra a very difficult task. In this respect, computational spectroscopy is becoming a needful and versatile tool for the assignment and interpretation of experimental spectra. It is in fact possible nowadays to model with relatively high accuracy the physical-chemical properties of complex molecules in different environments, and to link spectroscopic evidence directly to the structural and dynamical properties of optically or magnetically active solvated probes. In this Review, significant steps toward the simulation of entire spectra in condensed phases are presented together with some basic aspects of computational spectroscopy, which highlight how intramolecular and intermolecular degrees of freedom influence several spectroscopic parameters.

Harmonic and Anharmonic Vibrational Frequency Calculations with the Double-Hybrid B2PLYP Method: Analytic Second Derivatives and Benchmark Studies

This work aims to provide reliable benchmark data on the accuracy of harmonic and anharmonic vibrational frequencies computed with the B2PLYP double-hybrid density functional method. The exchange-correlation contributions required for the B2PLYP analytical second derivatives are presented here, which allow for the effective calculation of harmonic frequency as well as cubic and semidiagonal quartic force fields. The latter, in turn, are necessary to compute the anharmonic vibrational frequencies with the perturbative approach (VPT2). The quality of harmonic vibrational frequencies computed in conjunction with basis sets of double- to quadruple-ζ quality has been checked against reference data from the F38 benchmark set. Then, for an additional set of small closed- and open-shell systems, both harmonic frequencies and anharmonic contributions computed at the B2PLYP/N07D and the B2PLYP/aug-cc-pVTZ levels have been compared to their CCSD(T) counterparts. Moreover, for selected medium-size molecules (furan, pyrrole, thiophene, uracil, anisole, phenol, and pyridine), anharmonic frequencies have been compared to well established experimental results. Such benchmark studies have shown that the B2PLYP/N07D model provides good quality harmonic frequencies and describes correctly anharmonic contributions, the latter being of similar accuracy to their B3LYP/N07D counterparts, but obtained at significantly larger computational cost. Additionally, increased accuracy can be obtained by adopting hybrid models where the B2PLYP/N07D anharmonic contributions are combined with harmonic frequencies computed with more accurate quantum mechanical (QM) approaches or by B2PLYP with larger basis sets. This work confirmed also that most of the recently developed density functionals are significantly less suited for vibrational computations, while the B2PLYP method can be recommended for spectroscopic studies where a good accuracy of vibrational properties is required.