Interaction-induced electric properties and cooperative effects in model systems

From Oligo(Phenyleneethynylene) Monomers to Supramolecular Helices: The Role of Intermolecular Interactions in Aggregation

Supramolecular helices that arise from the self-assembly of small organic molecules via non-covalent interactions play an important role in the structure and properties of the corresponding materials. Here we study the supramolecular helical aggregation of oligo(phenyleneethynylene) monomers from a theoretical point of view, always guiding the studies with experimentally available data. In this way, by systematically increasing the number of monomer units, optimized n-mer geometries are obtained along with the corresponding absorption and circular dichroism spectra. For the geometry optimizations we use density functional theory together with the B3LYP-D3 functional and the 6–31G** basis set. For obtaining the spectra we resort to time-dependent density functional theory using the CAM-B3LYP functional and the 3–21G basis set. These combinations of density functional and basis set were selected after systematic convergence studies. The theoretical results are analyzed and compared to the experimentally available spectra, observing a good agreement.

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.

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.

Accurate calculation of optically induced birefringences in chiral systems using efficient polarized basis sets

Using state-of-the-art ab initio methodology, we evaluate universal molecular parameters entering the expressions for various optically induced birefringences in chiral fluids. For this, we use the single and double excitation coupled cluster (CCSD) theory together with Dunning's augmented correlation consistent polarized basis sets of increasing size. As this is the first time these parameters are evaluated for chiral molecules using the CCSD approach, we choose possibly small test systems: a model asymmetric methane molecule, and (R)-fluoro-oxirane. With this choice, the convergence of the molecular parameters with the increase of the basis set size is investigated in detail. The results are compared to those obtained with the LPol-n (n = ds, dl, fs) and the ORP basis sets, and to the corresponding Density Functional Theory (DFT) counterparts. We can conclude that the DFT medium constant values are considerably far from coupled cluster, and therefore, aware of the known excellent performance of the CCSD method in the evaluation of various dynamic responses, we do not recommend the former methodology for accurate evaluation of the present properties. Regarding basis set convergence, the performace of the LPol-ds basis set is more efficient than that of the d-aug-cc-pVDZ set, and therefore, the former basis set can be a good choice when dealing with the study of larger systems.

Interaction-induced electric (hyper)polarizability in the dihydrogen-neon pair: basis set and electron correlation effects

We investigated the interaction (hyper)polarizability of neon-dihydrogen pairs by performing high-level ab initio calculations with atom/molecule-specific, purpose-oriented Gaussian basis sets. We obtained interaction-induced electric properties at the SCF, MP2, and CCSD levels of theory. At the CCSD level, for the T-shaped configuration, around the respective potential minimum of 6.437 a₀, the interaction-induced mean first hyperpolarizability varies for 5 <  R/a₀ < 10 as[Formula: see text]Again, at the CCSD level, but for the L-shaped configuration around the respective potential minimum of 6.572 a₀, this property varies for 5 <  R/a₀ < 10 as[Formula: see text] Graphical Abstract Interaction-induced mean dipole polarizability ([Formula: see text]) for the T-shaped configuration of H₂-Ne calculated at the SCF, MP2, and CCSD levels of theory.

Partitioning of interaction-induced nonlinear optical properties of molecular complexes. I. Hydrogen-bonded systems

Interaction-induced electronic and vibrational (hyper)polarizabilities were decomposed into different interaction types (electrostatic, exchange, induction and dispersion).

On the physical origins of interaction-induced vibrational (hyper)polarizabilities

This paper presents the results of a pioneering exploration of the physical origins of vibrational contributions to the interaction-induced electric properties of molecular complexes. In order to analyze the excess nuclear relaxation (hyper)polarizabilities, a new scheme was proposed which relies on the computationally efficient Bishop-Hasan-Kirtman method for determining the nuclear relaxation contributions to electric properties. The extension presented herein is general and can be used with any interaction-energy partitioning method. As an example, in this study we employed the variational-perturbational interaction-energy decomposition scheme (at the MP2/aug-cc-pVQZ level) and the extended transition state method by employing three exchange-correlation functionals (BLYP, LC-BLYP, and LC-BLYP-dDsC) to study the excess properties of the HCN dimer. It was observed that the first-order electrostatic contribution to the excess nuclear relaxation polarizability cancels with the negative exchange repulsion term out to a large extent, resulting in a positive value of Δα(nr) due to the contributions from the delocalization and the dispersion terms. In the case of the excess nuclear relaxation first hyperpolarizability, the pattern of interaction contributions is very similar to that for Δα(nr), both in terms of their sign as well as relative magnitude. Finally, our results show that the LC-BLYP and LC-BLYP-dDsC functionals, which yield smaller values of the orbital relaxation term than BLYP, are more successful in predicting excess properties.

On the performance of long-range-corrected density functional theory and reduced-size polarized LPol-n basis sets in computations of electric dipole (hyper)polarizabilities of π-conjugated molecules

Static longitudinal electric dipole (hyper)polarizabilities are calculated for six medium-sized π-conjugated organic molecules using recently developed LPol-n basis set family to assess their performance. Dunning's correlation-consistent basis sets of triple-ζ quality combined with MP2 method and supported by CCSD(T)/aug-cc-pVDZ results are used to obtain the reference values of analyzed properties. The same reference is used to analyze (hyper)polarizabilities predicted by selected exchange-correlation functionals, particularly those asymptotically corrected.

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.