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.

The influence of the spatial confinement on the electronic and vibrational contributions to longitudinal electric-dipole properties of model linear molecules including HCN, HCCH and CO ₂ is analyzed.

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 $\alpha$,$\omega$-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.