Meta-hybrid density functional theory calculations using M06-2X/6-31+G(d,p) and M06-2X/6-311+G(d,p) levels of theory have been performed to understand the strength of C-H(…)π interactions of two possible types for benzene-acetylene, 1,3,5-trifluorobenzene-acetylene and coronene-acetylene complexes. Our study reveals that the C-H(...)π interaction complex where acetylene located above to the center of benzene ring (classical T-shaped) is the lowest energy structure. This structure is twice more stable than the configuration characterized by H atom of benzene interacting with the π-cloud of acetylene. The binding energy of 2.91 kcal/mol calculated at the M06-2X/6-311+G(d,p) level for the lowest energy configuration (1A) is in very good agreement with the experimental binding energy of 2.7 ± 0.2 kcal/mol for benzene-acetylene complex. Interestingly, the C-H(...)π interaction of acetylene above to the center of the aromatic ring is not the lowest energy configuration for 1,3,5-trifluorobenzene-acetylene and coronene-acetylene complexes. The lowest energy configuration (2A) for the former complex possesses both C-H(...)π interaction and C-H(...)F hydrogen bond, while the lowest energy structure for the coronene-acetylene complex involves both π-π and C-H(...)π interactions. C-H stretching vibrational frequencies and the frequency shifts are reported and analyzed for all of the configurations. We observed red-shift of the vibrational frequency for the stretching mode of the C-H bond that interacts with the π-cloud. Acetylene in the lowest-energy structures of the complexes exhibits significant red-shift of the C-H stretching frequency and change in intensity of the corresponding vibrational frequency, compared to bare acetylene. We have examined the molecular electrostatic potential on the surfaces of benzene, 1,3,5-trifluorobenzene, coronene and acetylene to explain the binding strengths of various complexes studied here.

Calculations revealed that only one of the three possible tautomers of AFL is stable, both in the gas phase and water. The electronic properties of aflatoxicol are calculated as similar to aflatoxin B1 and this may be an explanation of similar carcinogenicity and toxicity of these compounds, which has been proved by experimental results.

We have evaluated the performance of two of the recently developed density functionals (M06-2X and B2PLYP-D), which are widely used, by considering three important prototype systems, including benzene-acetylene, benzene-methane, and benzene-ammonia, possessing C-H···π or N-H···π interactions. Computational results are compared with the available experimental data. Considered density functionals are from two different classes: hybrid meta density functional (M06-2X) and double hybrid density functional (B2PLYP-D). The performance of a range of basis sets (6-31G(d), 6-31+G(d), 6-31+G(d,p), 6-311G(d,p), 6-311+G(d,p), aug-cc-pVXZ (X = D, T, Q)) with the above-mentioned two density functionals was evaluated. Comparison of the results includes Pople's basis sets versus Dunning's correlation consistent basis sets with the M06-2X and B2PLYP-D functionals considered in this study. The basis set effect on geometrical parameters, dissociation energies, and selected vibrational frequency shifts was thoroughly analyzed. We have addressed whether the counterpoise corrections with geometry optimizations and vibrational frequencies are important. Our computational study reveals that calculations carried out with smaller basis sets very well reproduce the reported experimental values of dissociation energies. The present study also shows that using the very large Dunning's correlation consistent basis set worsens the results. The necessity of including counterpoise correction for binding energies depends on the system and the type of method used. In general, vibrational frequency calculations using these DFT functionals generate characteristic red shifts for the C-H···π or N-H···π interactions in the complexes.

Dimers of cytosine and its N¹-methylated counterpart were investigated in gas-phase and in various solvents including chloroform, dimethylsulfoxide, and water. The studies were performed at DFT/M06-2X/6-31+G(d,p) level of theory. Relative stabilities of tautomers of cytosine solvated explicitly by a small number of solvent molecules were evaluated. Further solvation effect calculations for homodimers were carried out with conductor-like polarizable continuum model (CPCM). H-NMR shifts and IR frequencies for optimized structures were calculated and compared with available experimental data.