Many-body interactions and the electric response of hydrogen-bonded molecular chains

DFT and spatial confinement: a benchmark study on the structural and electrical properties of hydrogen bonded complexes

An extended set of 37 exchange correlation functionals, representing different DFT approximations, has been evaluated on a difficult playground represented by the dipole moment (μ_z), polarizability (α_zz), first hyperpolarizability (β_zzz), and the corresponding interaction-induced electrical properties (Δμ_z, Δα_zz, Δβ_zzz) of spatially confined hydrogen bonded (HB) dimers. A two-dimensional harmonic oscillator potential was used to exert the effect of spatial restriction. The performance of DFT methods in predicting hydrogen bond lengths in the studied molecular complexes upon confinement has also been examined. The data determined using a high-level CCSD(T) method serve as a reference. The conducted analyses allow us to conclude that methods rooted in DFT constitute a precise tool for the calculation of μ_z and α_zz as well as Δμ_z and Δα_zz, as most of the tested functionals provide results affected by rather small relative errors. On the other hand, an accurate description of the nonlinear optical response of the studied HB systems remains a great challenge for most of the analyzed DFT functionals, both in vacuum and in the presence of an analytical confining potential. Some of the tested DFT methods are found to be prone to catastrophic failure in the prediction of β_zzz as well as Δβ_zzz. The obtained results indicate that there is no great chasm in performance between functionals belonging to different DFT approximations or functionals including different amount of Hartree-Fock exchange when the values of dipole moment and first hyperpolarizability as well as the corresponding interaction-induced electrical properties are considered. However, a higher fraction of Hartree-Fock exchange improves the quality of predictions of α_zz and Δα_zz. Additionally, it has been shown that only three functionals from the examined set, namely B2PLYP, B3LYP and ωB97X-D, provide highly accurate structural parameters for the investigated systems. Of significant importance is the conclusion that the ωB97X-D functional, representing a modern and highly parametrized range-separated hybrid, demonstrates the most coherent behavior, showing rather small deviations from the reference data in the case of μ_z, α_zz, Δμ_z and Δα_zz as well as the structural parameters of the studied HB dimers. Moreover, our results indicate that the presence of spatial confinement has a rather small effect on the performance of DFT methods.

Hydrogen bonding inside and outside carbon nanotubes: HF dimer as a case study

In this theoretical work we analyze the noncovalent interactions of molecular complexes formed between the hydrogen bonded HF dimer and single-walled carbon nanotubes (SWCNTs) of different diameters. In particular, the interaction energies of: (i) spatially confined hydrogen fluoride molecules and (ii) HF dimer and the exterior or interior of SWCNTs are investigated. The computations are carried out in a supermolecular manner using the M06-2X exchange-correlation functional. In order to establish the influence of mutual orientation of the hydrogen fluoride dimer and molecular carbon cages on the analyzed energetic parameters energy scans are performed. Furthermore, changes in the charge distribution of the investigated endo- and exohedral complexes are studied employing the Natural Bond Orbital analysis. Among others, the position of the HF dimer with respect to the carbon cages proves to have a significant influence on the analyzed quantities. The results of our study also indicate that the HF dimer interacts stronger with the interior rather than the exterior of SWCNTs. Moreover, a substantial enhancement of the basis set superposition error is disclosed.

On the potential application of DFT methods in predicting the interaction-induced electric properties of molecular complexes. Molecular H-bonded chains as a case of study

A detailed analysis of the selected DFT functionals for the calculations of interaction-induced dipole moment, polarizability and first-order hyperpolarizability has been carried out. The hydrogen-bonded model chains consisting of HF, H(2)CO and H(3)N molecules have been chosen as a case study. The calculations of the components of the static electric properties using the diffuse Dunning's basis set (aug-cc-pVDZ) have been performed employing different types of density functionals (B3LYP, LC-BLYP, PBE0, M06-2X and CAM-B3LYP). Obtained results have been compared with those gained at the CCSD(T) level of theory. The counterpoise correction scheme, namely site-site function counterpoise, has been applied in order to eliminate basis set superposition error. The performed tests allow to conclude that the DFT functionals can provide a useful tool for prediction of the interaction-induced electric properties, however a caution has to be urged to their decomposition to the two- and many-body terms.