Hyperpolarizabilities of Push-Pull Chromophores in Solution: Interplay between Electronic and Vibrational Contributions.

Contemporary design of new organic non-linear optical (NLO) materials relies to a large extent on the understanding of molecular and electronic structure-property relationships revealed during the years by available computational approaches. The progress in theory-hand-in-hand with experiment-has enabled us to identify and analyze various physical aspects affecting the NLO responses, such as the environmental effects, molecular vibrations, frequency dispersion, and system dynamics. Although it is nowadays possible to reliably address these effects separately, the studies analyzing their mutual interplay are still very limited. Here, we employ density functional theory (DFT) methods in combination with an implicit solvent model to examine the solvent effects on the electronic and harmonic as well as anharmonic vibrational contributions to the static first hyperpolarizability of a series of push-pull α,ω-diphenylpolyene oligomers, which were experimentally shown to exhibit notable second-order NLO responses. We demonstrate that the magnitudes of both vibrational and electronic contributions being comparable in the gas phase significantly increase in solvents, and the enhancement can be, in some cases, as large as three- or even four-fold. The electrical and mechanical anharmonic contributions are not negligible but cancel each other out to a large extent. The computed dynamic solute NLO properties of the studied systems are shown to be in a fair agreement with those derived from experimentally measured electric-field-induced second-harmonic generation (EFISHG) signals. Our results substantiate the necessity to consider concomitantly both solvation and vibrational effects in modeling static NLO properties of solvated systems.

HeH+ under Spatial Confinement.

In the present study, the influence of spatial confinement on the bond length as well as dipole moment, polarizability and (hyper)polarizabilities of HeH+ ion was analyzed. The effect of spatial confinement was modelled by cylindrically symmetric harmonic oscillator potential, that can be used to mimic high pressure conditions. Based on the conducted research it was found that the spatial confinement significantly affects the investigated properties. Increasing the confinement strength leads to a substantial decrease of their values. This work may be of particular interest for astrochemistry as HeH+ is believed to be the first compound to form in the Universe.

Much of a Muchness: On the Origins of Two- and Three-Photon Absorption Activity of Dipolar Y-Shaped Chromophores.

The molecular origin of two- (2PA) and three-photon absorption (3PA) activity in three experimentally studied chromophores, prototypical dipolar systems, is investigated. To that end, a generalized few-state model (GFSM) formula is derived for the 3PA transition strength for nonhermitian theories and employed at the coupled-cluster level of theory. Using various computational techniques such as molecular dynamics, linear and quadratic response theories, and GFSM, an in-depth analysis of various optical channels involved in 2PA and 3PA processes is presented. It is found that the four-state model involving the second and third excited singlet states as intermediates is the smallest model among all considered few-state approximations that produces 2PA and 3PA transition strengths (for S0 → S1 transition) close to the reference results. By analyzing various optical channels appearing in these models and involved in studied multiphoton processes, we found that the 2PA and 3PA activities in all the three chromophores are dominated and hence controlled by the dipole moment of the final excited state. The similar origins of the 2PA and the 3PA in these prototypical dipolar chromophores suggest transferability of structure-property relations from the 2PA to the 3PA domain.

A Look at the Spatial Confining Effect on the Molecular Electrostatic Potential (MEP)-A Case Study of the HF and BrCN Molecules.

In this theoretical study, we report on the molecular electrostatic potential (MEP) of titled molecules confined by repulsive potentials of cylindrical symmetry mimicking a topology. Our calculations show that the spatial restriction significantly changes the picture of the MEP of molecules in a quantitative and qualitative sense. In particular, the drastic changes in the MEP as a function of the strength of spatial confinement are observed for the BrCN molecule. This preliminary study is the first step in the investigation of the behavior of the MEP of molecular systems under orbital compression.

An effective potential for Frenkel excitons.

Excitation energy transfer (EET) is a ubiquitous process in life and materials sciences. Here, a new and computationally efficient method of evaluating the electronic EET couplings between interacting chromophores is introduced that is valid in a wide range of intermolecular distances. The proposed approach is based on the effective elimination of electron repulsion integrals from the excitonic Hamiltonian matrix elements via the density-fitting approach and distributed multipole approximation. The excitonic Hamiltonian represented in a basis including charge transfer (CT) states is re-cast in terms of the effective one-electron potential functions (EOPs) and adapted into the effective fragment parameter (EFP) framework. Calculations for model systems indicate that the speedup of at least three orders of magnitude, as compared to the state-of-the-art methods, can be achieved while maintaining the accuracy of the EET couplings even at short intermolecular distances.

Ab initio effective one-electron potential operators: Applications for charge-transfer energy in effective fragment potentials.

The concept of effective one-electron potentials (EOPs) has proven to be extremely useful in efficient description of electronic structure of chemical systems, especially extended molecular aggregates such as interacting molecules in condensed phases. Here, a general method for EOP-based elimination of electron repulsion integrals is presented, that is tuned toward the fragment-based calculation methodologies such as the second generation of the effective fragment potentials (EFP2) method. Two general types of the EOP operator matrix elements are distinguished and treated either via the distributed multipole expansion or the extended density fitting (DF) schemes developed in this work. The EOP technique is then applied to reduce the high computational costs of the effective fragment charge-transfer (CT) terms being the bottleneck of EFP2 potentials. The alternative EOP-based CT energy model is proposed, derived within the framework of intermolecular perturbation theory with Hartree-Fock noninteracting reference wavefunctions, compatible with the original EFP2 formulation. It is found that the computational cost of the EFP2 total interaction energy calculation can be reduced by up to 38 times when using the EOP-based formulation of CT energy, as compared to the original EFP2 scheme, without compromising the accuracy for a wide range of weakly interacting neutral and ionic molecular fragments. The proposed model can thus be used routinely within the EFP2 framework.

Infrared Spectra of Hydrogen-Bonded Molecular Complexes Under Spatial Confinement.

Infrared (IR) spectroscopy is commonly used in chemical laboratories to study the geometrical structure of molecules and molecular complexes. The analysis of experimental IR spectra can nowadays be reliably supported by the results of quantum-chemical computations as vibrational frequencies and corresponding vibrational transition intensities are routinely calculated using harmonic approximation by virtually all quantum chemistry packages. In the present study we combine the methodology of computing vibrational spectra using high-level electron correlation treatments with an analytical potential-based approach to take into account spatial confinement effects. Using this approach, we perform a pioneering analysis of the impact of the spatial confinement caused by a cylindrical harmonic oscillator potential on the harmonic vibrational transition intensities and frequencies of two hydrogen-bonded complexes: HCN…HCN and HCN…HNC. The emphasis is put on the largest-intensity bands, which correspond to the stretching vibrations. The obtained results demonstrate that embedding the molecular complexes in an external confining potential causes significant changes of transition intensities and vibrational frequencies.

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.

Assessment of DFT for endohedral complexes' dipole moment: PNO-LCCSD-F12 as a reference method.

We present a systematic evaluation of the performance of a wide range of exchange-correlation functionals and related dispersion correction schemes for the computation of dipole moments of endohedral complexes, formed through the encapsulation of an AB molecule (AB = LiF, HCl) inside carbon nanotubes (CNTs) of different diameter. The consistency and accuracy of (i) generalized gradient approximation, (ii) meta GGA, (iii) global hybrid, and (iv) range-separated hybrid density functionals are assessed. In total, 37 density functionals are tested. The results obtained using the highly accurate pair natural orbitals based explicitly correlated local coupled cluster singles doubles (PNO-LCCSD-F12) method of Werner and co-workers [Schwilk et al., J. Chem. Theory Comput., 2017, 13, 3650; Ma et al., J. Chem. Theory Comput., 2017, 13, 4871] with the aug-cc-pVTZ basis set serve as a reference. The static electric dipole moment is computed via the finite field response or, when possible, as the expectation value of the dipole operator. Among others, it is shown that functionals belonging to the class of range-separated hybrids, provide results closest to the coupled cluster reference data. In particular, the ωB97X as well as the M11 functional may be considered as a promising choice for computing electric properties of noncovalent endohedral complexes. On the other hand, the worst performance was found for the functionals which do not include the Hartree-Fock exchange. The analysis of both the coupled cluster and the DFT results indicates a strong coupling of dispersion and polarization that may also explain why lower level DFT methods, as well as Hartree-Fock and MP2, cannot yield dipole moments beyond a qualitative agreement with the higher order reference data. Interestingly, the much smaller and less systematically constructed basis sets of Pople of moderate size provide results of accuracy at least comparable with the extended Dunning's aug-cc-pVTZ basis set.

Benchmarking the Performance of Exchange-Correlation Functionals for Predicting Two-Photon Absorption Strengths.

The present work investigates the performance of exchange-correlation functionals in the prediction of two-photon absorption (2PA) strengths. For this purpose, we considered six common functionals used for studying 2PA processes and tested these on six organoboron chelates. The set consisted of two semilocal (PBE and BLYP), two hybrid (B3LYP and PBE0), and two range-separated (LC-BLYP and CAM-B3LYP) functionals. The RI-CC2 method was chosen as a reference level and was found to give results consistent with the experimental data that are available for three of the molecules considered. Of the six exchange-correlation functionals studied, only the range-separated functionals predict an ordering of the 2PA strengths that is consistent with experimental and RI-CC2 results. Even though the range-separated functionals predict correct relative trends, the absolute values for the 2PA strengths are underestimated by a factor of 2-6 for the molecules considered. An in-depth analysis, on the basis of the derived generalized few-state model expression for the 2PA strength for a coupled-cluster wave function, reveals that the problem with these functionals can be linked to underestimated excited-state dipole moments and, to a lesser extent, overestimated excitation energies. The semilocal and hybrid functionals exhibit less predictable errors and a variation in the 2PA strengths in disagreement with the reference results. The semilocal and hybrid functionals show smaller average errors than the range-separated functionals, but our analysis reveals that this is due to fortuitous error cancellation between excitation energies and the transition dipole moments. Our results constitute a warning against using currently available exchange-correlation functionals in the prediction of 2PA strengths and highlight the need for functionals that correctly describe the electron density of excited electronic states.

Spectral and physicochemical properties of difluoroboranyls containing N,N-dimethylamino group studied by solvatochromic methods.

The solvatochromism of the dyes was analyzed based on the four-parameter scale including: polarizability (SP), dipolarity (SdP), acidity (SA) and basicity (SB) parameters by method proposed by Catalán. The change of solvent to more polar caused the red shift of absorption and fluorescence band position. The frequency shifts manifest the change in the dipole moment upon excitation. The ground-state dipole moment of the difluoroboranyls was estimated based on changes in molecular polarization with temperature. Moreover, the Stokes shifts were used to calculate the excited state dipole moments of the dyes. For the calculation, the ground-state dipole moments and Onsager cavity radius were also determined theoretically using density functional theory (DFT). The experimentally determined excited-state dipole moments for the compounds are higher than the corresponding ground-state values. The increase in the dipole moment is described in terms of the nature of the excited state.

On the methodology of the determination of charge concentration dependent mobility from organic field-effect transistor characteristics.

We developed a new methodology for determining charge concentration dependent mobility from organic field-effect transistor (OFET) characteristics, applicable for semiconducting polymers with structural and energy disorder. We show that basic formulae recommended by the "IEEE Standard for Test Methods for the Characterization of Organic Transistors and Materials" for the determination of the field-effect mobility as obtained from the slope ISD1/2vs. VSG (in the saturation regime) or from the transconductance dISD/dVSG (in the linear regime) are not suitable for materials with concentration dependent charge carrier mobility. We propose alternative expressions, which can be directly analytically derived from the drift-diffusion equation with the mobility explicitly dependent on the charge concentration. This methodology for mobility determination was used for analysis of the experimental data obtained for a poly(3-hexylthiophene)-based OFET with the bottom gate-bottom SD electrode configuration.

Nonempirical Simulations of Inhomogeneous Broadening of Electronic Transitions in Solution: Predicting Band Shapes in One- and Two-Photon Absorption Spectra of Chalcones.

We have examined several approaches relying on the Polarizable Embedding (PE) scheme to predict optical band shapes for two chalcone molecules in methanol solution. The PE-TDDFT and PERI-CC2 methods were combined with molecular dynamics simulations, where the solute geometry was kept either as rigid, flexible or partly-flexible (restrained) body. The first approach, termed RBMD-PE-TDDFT, was employed to estimate the inhomogeneous broadening for subsequent convolution with the vibrationally-resolved spectra of the molecule in solution determined quantum-mechanically (QM). As demonstrated, the RBMD-PE-TDDFT/QM-PCM approach delivers accurate band widths, also reproducing their correct asymmetric shapes. Further refinement can be obtained by the estimation of the inhomogeneous broadening using the RBMD-PERI-CC2 method. On the other hand, the remaining two approaches (FBMD-PE-TDDFT and ResBMD-PE-TDDFT), which lack quantum-mechanical treatment of molecular vibrations, lead to underestimated band widths. In this study, we also proposed a simple strategy regarding the rapid selection of the exchange-correlation functional for the simulations of vibrationally-resolved one- and two-photon absorption spectra based on two easy-to-compute metrics.

Quantifying the Performances of DFT for Predicting Vibrationally Resolved Optical Spectra: Asymmetric Fluoroborate Dyes as Working Examples.

This article aims at a quantitative assessment of the performances of a panel of exchange-correlation functionals, including semilocal (BLYP and PBE), global hybrids (B3LYP, PBE0, M06, BHandHLYP, M06-2X, and M06-HF), and range-separated hybrids (CAM-B3LYP, LC-ωPBE, LC-BLYP, ωB97X, and ωB97X-D), in predicting the vibrationally resolved absorption spectra of BF2-carrying compounds. To this end, for 19 difluoroborates as examples, we use, as a metric, the vibrational reorganization energy (λvib) that can be determined based on the computationally efficient linear coupling model (a.k.a. vertical gradient method). The reference values of λvib were determined by employing the CC2 method combined with the cc-pVTZ basis set for a representative subset of molecules. To validate the performances of CC2, comparisons with experimental data have been carried out as well. This study shows that the vibrational reorganization energy, involving Huang-Rhys factors and normal-mode frequencies, can indeed be used to quantify the reliability of functionals in the calculations of the vibrational fine structure of absorption bands, i.e., an accurate prediction of the vibrational reorganization energy leads to absorption band shapes better fitting the selected reference. The CAM-B3LYP, M06-2X, ωB97X-D, ωB97X, and BHandHLYP functionals all deliver vibrational reorganization energies with absolute relative errors smaller than 20% compared to CC2, whereas 10% accuracy can be achieved with the first three functionals. Indeed, the set of examined exchange-correlation functionals can be divided into three groups: (i) BLYP, B3LYP, PBE, PBE0, and M06 yield inaccurate band shapes (λvib,TDDFT < λvib,CC2), (ii) BHandHLYP, CAM-B3LYP, M06-2X, ωB97X, and ωB97X-D provide accurate band shapes (λvib,TDDFT ≈ λvib,CC2), and (iii) LC-ωPBE, LC-BLYP, and M06-HF deliver rather poor band topologies (λvib,TDDFT > λvib,CC2). This study also demonstrates that λvib can be reliably estimated using the CC2 model and the relatively small cc-pVDZ basis set. Therefore, the linear coupling model combined with the CC2/cc-pVDZ level of theory can be used as a very efficient approach to determine λvib values that can be used to select the most adequate functional for more accurate vibronic calculations, e.g., including more refined models and environmental effects.

Vibrational nonlinear optical properties of spatially confined weakly bound complexes.

This study focuses on the theoretical description of the influence of spatial confinement on the electronic and vibrational contributions to (hyper)polarizabilities of two dimeric hydrogen bonded systems, namely HCNHCN and HCNHNC. A two-dimensional analytical potential is employed to render the confining environment (e.g. carbon nanotube). Based on the results of the state-of-the-art calculations, performed at the CCSD(T)/aug-cc-pVTZ level of theory, we established that: (i) the influence of spatial confinement increases with increasing order of the electrical properties, (ii) the effect of spatial confinement is much larger in the case of the electronic than vibrational contribution (this holds for each order of the electrical properties) and (iii) the decrease in the static nuclear relaxation first hyperpolarizability upon the increase of confinement strength is mainly due to changes in the harmonic term, however, in the case of nuclear relaxation second hyperpolarizability the anharmonic terms contribute more to the drop of this property.

On the calculations of the static electronic dipole (hyper)polarizability for the free and spatially confined H.

The quality of the static electronic dipole (hyper)polarizabilities (α and γ) calculated using the coupled cluster with singles and doubles (CCSD) method and various basis sets for the isolated hydrogen anion (H-) was verified against highly accurate results obtained within the variation-perturbation scheme with the optimized explicitly correlated Gaussian functions. Moreover, the influence of spatial confinement on α and γ of H- was analyzed. The effect of confinement (orbital compression) was modeled by the spherically symmetric harmonic oscillator potential. Our results demonstrate that in order to obtain the value of α with a relative error less than 10% and correct order of magnitude of γ for the unconfined H- at the CCSD level of theory, at least triply augmented correlation consistent basis set (t-aug-cc-pVXZ, X = D,T,Q,5,6) has to be used. Furthermore, it was found that the spatial confinement significantly affects α and γ of H-.

About the nature of halogen bond interaction under the spatial confinement.

Nowadays, much attention is put toward the description of noncovalent complexes exposed to the high pressure or embedded in confining environments. Such conditions may strongly modify the physical and chemical properties of molecular systems. This study focuses on the theoretical description of the confinement induced changes in geometry and energetic parameters of the halogen bonded FCl⋯CNF complex. A model analytical potential is applied to render the effect of orbital compression. In order to analyze the nature of halogen bond interaction, in the presence of spatial confinement, the supermolecular approach together with the symmetry-adapted perturbation theory is used. Furthermore, a thorough analysis of topological parameters, characterizing the halogen bond upon orbital compression, is performed within the quantum theory of atoms in molecules. The calculations are carried out using the ωB97x and CCSD(T) methods in connection with the aug-cc-pVTZ basis set. Among others, the obtained results indicate that the spatial confinement not only modifies the nature of halogen bond interaction but also induces the appearance of a completely new form of the studied FCl⋯CNF system.

Modelling of the charge carrier mobility in disordered linear polymer materials.

We introduced a molecular-scale description of disordered on-chain charge carrier states into a theoretical model of the charge carrier transport in polymer semiconductors. The presented model combines the quantum mechanical approach with a semi-classical solution of the inter-chain charge hopping. Our model takes into account the significant local anisotropy of the charge carrier mobility present in linear conjugated polymers. Contrary to the models based on the effective medium approximation, our approach allowed avoiding artefacts in the calculated concentration dependence of the mobility originated in its problematic configurational averaging. Monte Carlo numerical calculations show that, depending on the degree of the energetic and structural disorder, the charge carrier mobility increases significantly with increasing charge concentration due to trap filling. At high charge carrier concentrations, the effect of the energetic disorder disappears and the mobility decreases slightly due to the lower density of unoccupied states available for the hopping transport. It could explain the experimentally observed mobility degradation in organic field-effect transistors at high gate voltage.

Two-photon absorption of the spatially confined LiH molecule.

In the present contribution we study the influence of spatial restriction on the two-photon dipole transitions between the X1Σ+ and A1Σ+ states of lithium hydride. The bond-length dependence of the two-photon absorption strength is also analyzed for the first time in the literature. The highly accurate multiconfiguration self-consistent field (MCSCF) method and response theory are used to characterize the electronic structure of the studied molecule. In order to render the effect of orbital compression we apply a two-dimensional harmonic oscillator potential, mimicking the topology of cylindrical confining environments (e.g. carbon nanotubes, quantum wires). Among others, the obtained results provide evidence that at large internuclear distances the TPA response of lithium hydride may be significantly enhanced and this effect is much more pronounced upon embedding of the LiH molecule in an external confining potential. To understand the origin of the observed variation in the two-photon absorption response a two-level approximation is employed.

Two-photon absorption of BF2-carrying compounds: insights from theory and experiment.

This communication presents a structure-property study of a few novel pyridine-based difluoroborate compounds with a N-BF2-O core, which exhibit outstanding fluorescence properties. To exploit their potential for two-photon bioimaging, relationships between the two-photon action cross section and systematic structural modifications have been investigated and unravelled.