Synthesis, characterization and NLO properties of 1,4-phenylenediamine-based Schiff bases: a combined theoretical and experimental approach

In the current study, three novel 1,4-phenylenediamine-based chromophores (3a-3c) were synthesized and characterized and then their nonlinear optical (NLO) characteristics were explored theoretically. The characterization was done by spectroscopic analysis, i.e. FT-IR, UV-Visible, and NMR spectroscopy, and elemental analysis. Notably, these chromophores exhibited UV-Visible absorption within the range of 378.635-384.757 nm in acetonitrile solvent. Additionally, the FMO findings for 3a-3c revealed the narrowest band gap (4.129 eV) for 3c. The GRPs for these chromophores were derived from HOMO-LUMO energy values, which showed correspondence with FMO results by depicting a minimum hardness (2.065 eV) for 3c. Among these compounds, 3c displayed the highest nonlinear behavior with maximum μtot, βtot and γtot values of 4.79 D, 8.00 × 10-30 and 8.13 × 10-34 a.u., respectively. Our findings disclosed that the synthesized 1,4-phenylenediamine chromophores may be considered promising candidates for nonlinear optical materials, showing potential applications in the realm of optoelectronic devices.

Theoretical study on porphyrin arch-tapes of carbonyl-inserted seven-membered rings with high nonlinear optical properties

Porphyrin tapes have attracted extensive attention because their fully conjugated π-networks act as nonlinear optical (NLO) materials. A family of Ni(II) and Zn(II) porphyrin arch-tapes that are connected by varying bridge (B) ligands (meso-meso β-β doubly linked dimer 1, meso-meso β-β β-β triply linked dimer 3, methylene-inserted dimer 2 and trimer 5, carbonyl-inserted dimer 4, trimer 6, and Zn(II) trimer 7) have been synthesized by a density functional theory (DFT) method. The results show that carbonyl-inserted arch-tapes significantly enhance second hyperpolarizability (γ), indicating that the remarkably contorted structure incorporated seven-membered ring(s) directly affect their NLO properties of our focus. Moreover, the electronic absorption spectra calculated for all studied complexes with time-dependent DFT theory (TDDFT) predict that carbonyl-inserted complex 4 contributes to a red-shift of the Q-band (160 nm) for the meso-meso β-β doubly linked complex 1. The third-order NLO responses and the electron transition properties strongly depend on the nature of the bridge (B) ligand, which means that an active involvement of the carbonyl group presents an advantage for its application in NLO materials.

A quantum chemical investigation of the second hyperpolarizability of p-nitroaniline

Recent measurements of the third harmonic scattering responses of molecules have given a new impetus for computing molecular second hyperpolarizabilities (γ) and for deducing structure-property relationships. This paper has employed a variety of wavefunction and density functional theory methods to evaluate the second hyperpolarizability of the p-nitroaniline prototypical push-pull π-conjugated molecule, addressing also numerical aspects, such as the selection of an integration grid and the impact of the order of differentiation vs the achievable accuracy by using the Romberg quadrature. The reliability of the different methods has been assessed by comparison to reference Coupled-Cluster Singles and Doubles with perturbative treatment of the Triples results. On the one hand, among wavefunction methods, the MP2 scheme offers the best accuracy/cost ratio for computing the static γ. On the other hand, using density functional theory, γ remains a challenging property to compute because all conventional, global hybrid or range-separated hybrid, exchange-correlation functionals underestimate static γ values by at least 15%. Even tuning the range-separating parameter to minimize the delocalization errors does not enable to improve the γ values. Nevertheless, the original double-hybrid B2-PLYP functional, which benefits from 27% of PT2 correlation and 53% Hartree-Fock exchange, provides accurate estimates of static γ values. Unfortunately, the best performing exchange-correlation functionals for γ are not necessarily reliable for the first hyperpolarizability, β, and vice versa. In fact, the β of p-nitroaniline (pNA) could be predicted, with a good accuracy, with several hybrid exchange-correlation functionals (including by tuning the range-separating parameter), but these systematically underestimate γ. As for γ, the MP2 wavefunction method remains the best compromise to evaluate the first hyperpolarizability of pNA at low computational cost.

Accurate prediction of global-density-dependent range-separation parameters based on machine learning

In this work, we develop an accurate and efficient XGBoost machine learning model for predicting the global-density-dependent range-separation parameter, ωGDD, for long-range corrected functional (LRC)-ωPBE. This ωGDDML model has been built using a wide range of systems (11 466 complexes, ten different elements, and up to 139 heavy atoms) with fingerprints for the local atomic environment and histograms of distances for the long-range atomic correlation for mapping the quantum mechanical range-separation values. The promising performance on the testing set with 7046 complexes shows a mean absolute error of 0.001 117 a0-1 and only five systems (0.07%) with an absolute error larger than 0.01 a0-1, which indicates the good transferability of our ωGDDML model. In addition, the only required input to obtain ωGDDML is the Cartesian coordinates without electronic structure calculations, thereby enabling rapid predictions. LRC-ωPBE(ωGDDML) is used to predict polarizabilities for a series of oligomers, where polarizabilities are sensitive to the asymptotic density decay and are crucial in a variety of applications, including the calculations of dispersion corrections and refractive index, and surpasses the performance of all other popular density functionals except for the non-tuned LRC-ωPBE. Finally, LRC-ωPBE (ωGDDML) combined with (extended) symmetry-adapted perturbation theory is used in calculating noncovalent interactions to further show that the traditional ab initio system-specific tuning procedure can be bypassed. The present study not only provides an accurate and efficient way to determine the range-separation parameter for LRC-ωPBE but also shows the synergistic benefits of fusing the power of physically inspired density functional LRC-ωPBE and the data-driven ωGDDML model.

Ultrafast photoinduced dynamics of a donor-([Formula: see text])bridge-acceptor based merocyanine dye

Merocyanine dyes are of great interest amongst researchers due to their nonlinear optical (NLO) properties and solvatochromism. Molecular structure of these dyes constitutes conjugated pathway between the donor and acceptor substituents, with lowest energy transition of [Formula: see text]-[Formula: see text]* character. To rationalize the design of these dyes and deduce structure-property relationship, it is eminent to unravel the excited state dynamics in these complex molecular structures in different solvents. Here we have studied excited state dynamics of a merocyanine dye known as HB194, which has shown commendable efficiency in small molecule based bulk heterojuction solar cells. We have employed femtosecond transient absorption in combination with the quantum chemistry calculations to unravel the solvent dependent charge transfer dynamics of HB194. The excited state decays of the HB194 in different solvents show multi-exponential components. The analysis of the time-resolved data reveals that the polar solvents induce conformationally relaxed intramolecular charge transfer state. In non-polar solvent cyclohexane, only solvent-stabilized ICT state is observed. Additionally, we observe an anomalously red-shifted emission in ethylene glycol centred at [Formula: see text] 750 nm. Our computational calculations suggest the presence of molecular dimers resulting into observed red-shifted emission band. Our work therefore underscores the importance of gathering molecular-level insight into the system-bath interactions for designing next generation merocynanine-based solvatochromic dyes.

Charge-Transfer Excitation within a Hybrid-(G)KS Framework through Cartesian Grid DFT

Organic molecules that exhibit charge-transfer (CT) excited states are known to play an important role in processes linked to electron transfer properties and molecular conductance. In this article, we present a simple technique based on "Becke's excitation theorem" that offers an accurate picture of these electronic states. It expresses the correlated energy splitting between triplet and its corresponding singlet states by a two-electron integral, which is numerically evaluated by our recently developed strategy on Cartesian grid. We first examine the consistency of our adopted numerical strategy to evaluate the integral with the original prescribed technique. Then we assess the method on weakly bound CT complexes with three different functionals (BLYP, B3LYP, and LC-BLYP). The accuracy on asymptotic limit of CT excitation is also explored. Finally in order to illustrate the strength and feasibility, it is further extended to a few "challenging" molecules. The method, when employed with hybrid B3LYP functional, turns out to be quite accurate to describe CT excitation energy.

Physical mechanisms of photoinduced charge transfer in neutral and charged donor-acceptor systems

In this paper, we provide visualization methods to reveal the physical mechanisms of photoinduced charge transfer in neutral and charged donor-acceptor systems. These visualization methods use the charge density difference and transition density matrix, which can promote deeper understanding of photoinduced charge transfer in donor-acceptor systems.

Density Functional Theory Study of Substitution Effects on the Second-Order Nonlinear Optical Properties of Lindquist-Type Organo-Imido Polyoxometalates

Density functional theory and time-dependent density functional theory have been enacted to investigate the effects of donor and acceptor on the first hyperpolarizability of Lindquist-type organo-imido polyoxometalates (POMs). These calculations employ a range-separated hybrid exchange-correlation functional (ωB97X-D), account for solvent effects using the implicit polarizable continuum model, and analyze the first hyperpolarizabilities by using the two-state approximation. They highlight the beneficial role of strong donors as well as of π-conjugated spacers (CH=CH rather than C≡C) on the first hyperpolarizabilities. Analysis based on the unit sphere representation confirms the one-dimensional push-pull π-conjugated character of the POMs substituted by donor groups and the corresponding value of the depolarization ratios close to 5. Furthermore, the use of the two-state approximation is demonstrated to be suitable for explaining the origin of the variations of the first hyperpolarizabilities as a function of the characteristics of a unique low-energy charge-transfer excited state and to attribute most of the first hyperpolarizability changes to the difference of dipole moment between the ground and that charge-transfer excited state.

Influence of a Substituted Methyl on the Photoresponsive Third-Order Nonlinear-Optical Properties Based on Azobenzene Metal Complexes

For studying the effect of a substituted group on the photoresponsive third-order nonlinear-optical (NLO) properties, photosensitive azobenzene derivative H₂L1 was first selected to construct metal complexes {[Zn₂(L1)₂(H₂O)₃]·2DMA)}_n (1) and {[Cd(L1)(4,4'-bpy)H₂O]·H₂O}_n (2). Then H₂L2 with a substituted methyl on the azobenzene ring was used to construct complexes {[Zn(L2)(4,4'-bpy)(H₂O)]}_n (3) and {[Cd(L2)(4,4'-bpy)(H₂O)]}_n (4). When the azobenzene moiety of the complexes is trans, the NLO behaviors of the complexes are the same. However, after the azobenzene moiety is excited by ultraviolet (UV) light to change from trans to cis, the substituted methyl increases the repulsion between two azobenzene rings in 3 and 4, thereby affecting their NLO behaviors. Therefore, the nonlinearity of the two types of complexes is different after UV irradiation. Density functional theory calculations support this result. The substituted methyl has a significant influence on the nonlinear absorption behaviors of 3 and 4. This work not only reports the examples of photoresponsive NLO materials based on metal complexes but also provides a new idea to deeply explore NLO properties.

Thermally activated delayed fluorescence: A critical assessment of environmental effects on the singlet-triplet energy gap

The effective design of dyes optimized for thermally activated delayed fluorescence (TADF) requires the precise control of two tiny energies: the singlet-triplet gap, which has to be maintained within thermal energy, and the strength of spin-orbit coupling. A subtle interplay among low-energy excited states having dominant charge-transfer and local character then governs TADF efficiency, making models for environmental effects both crucial and challenging. The main message of this paper is a warning to the community of chemists, physicists, and material scientists working in the field: the adiabatic approximation implicitly imposed to the treatment of fast environmental degrees of freedom in quantum-classical and continuum solvation models leads to uncontrolled results. Several approximation schemes were proposed to mitigate the issue, but we underline that the adiabatic approximation to fast solvation is inadequate and cannot be improved; rather, it must be abandoned in favor of an antiadiabatic approach.

NLOphoric imidazole-fused fluorescent anthraquinone dyes

The Z-scan and DFT techniques were explored to investigate the non-linear optical properties of anthraquinone fused imidazole-based D-π-A dyes. With the help of UV-visible spectral analysis, pH study, CV analysis, HOMO-LUMO interaction, MEP plots, and quinoidal character the influence of intramolecular charge transfer characteristics and H-bonding process on electronic and photophysical studies of anthraquinone derivatives were understood. The dyes 2-(4-(diethylamino)-2-hydroxyphenyl)-3H-anthra[1,2-d]imidazole-6,11-dione (AQ1) and 2-(2-hydroxynaphthalen-1-yl)-1H-anthra[1,2-d]imidazole-6,11-dione (AQ2) displayed a single emission with pronounced Stokes shift and thermal stability (upto 290 °C). The dye AQ1 exhibited strong charge transfer character which is explained by the ICT process leading to high nonlinear susceptibility χ⁽³⁾ in AQ1 3.43 × 10^-13 e.s.u relative to the dye AQ2 which has only ESIPT core. But, the dye AQ2 6.27 J cm^-2 showed better optical limiting value. The NLO properties of AQ1 and AQ2 were computed by DFT functionals based on Hartree Fork (HF) percentage exchange. The dye AQ1 exhibits noticeable NLO properties. The global hybrid functionals with HF composition beyond 50% (BHHLYP, M06-HF) and the long rang corrected functionals (CAM-B3LYP, ωB97, and ωB97X) demonstrated comparable NLO properties relative to B3LYP and PBE0. It was observed that the combined effect of ICT and the H-bonding cores enhance the NLO properties. The experimental findings (Z-scan) were successfully correlated with theoretical (DFT) results.

Computational study on the optical and NLO properties of donor–acceptor interaction molecules containing a triazolobenzothiadiazole or benzothiadiazole central acceptor core

The nonlinear optical properties of the studied compounds were studied with the aid of DFT calculations.

The inspiration and challenge for through-space charge transfer architecture: from thermally activated delayed fluorescence to non-linear optical properties

Organic molecules consisting of electron donor (D) and electron acceptor (A) subunits linked by π-conjugated bridges are promising building blocks for thermally activated delayed fluorescence (TADF) and non-linear optics (NLO) materials due to their intramolecular charge transfer (CT) processes in response to external stimuli. According to the electron interaction pattern, the CT process in D-π-A architectures can be divided into two categories, through-bond/-space charge transfer (TB/TSCT). To date, research into the TADF properties of TSCT characteristic molecules has since seen significant growth. In fact, TSCT characteristic materials show great advantages in such NLO responses. In this perspective, we first briefly introduced the basic principles of NLO and TADF effects. Successively, we discuss the influence of TBCT and TSCT patterns on NLO and TADF properties, especially for TSCT characteristic. In the final part, we address the diversity and potential advantages of TSCT characteristic molecules as high-performance NLO materials. With these, it is expected that the greater structural flexibility of spatial conjugation can bring more functionality to NLO materials in the future.

Understanding TADF: a joint experimental and theoretical study of DMAC-TRZ

TADF offers a promising way to harvest triplets in OLED for improved efficiency. To concurrently optimize the dye inside the matrix, a thorough experimental and theoretical study is presented of a the TADF dye addressing environmental effects.

Polymorphism Dependent 9-Phosphoanthracene Derivative Exhibiting Thermally Activated Delayed Fluorescence: A Computational Investigation

Polymorphs of anthracene derivatives exhibit diverse photophysical properties that can help to develop efficient organic-based photovoltaic devices. 10-Anthryl-9-phosphoanthracene (10-APA) shows different photophysical behaviors for the solid state due to its variety in crystalline arrangement. Herein, we investigate the ground and excited-state properties of the monomer and two different polymorphs of 10-APA from first-principles. Calculations reveal that strong spin-orbit coupling (SOC) between first excited singlet state (S₁) and triplet manifolds at their S₁-optimized geometries enabling the reverse intersystem crossing (RISC). The electron-vibration coupling (Huang-Rhys factor) in the excited state is the most relevant factor here. For both ISC and RISC, a similarity in Huang-Rhys factors for the molecular vibration along the π···π stacking at low-frequency region makes the rates effective. On the other side, the nonvanishing vibronic relaxation modes provide a relatively slower RISC rate in the red crystal. However, for the red crystal, small reorganization energy (λ) and large Huang-Rhys factor toward S₁ → S₀ conversion reduce nonradiative decay, leading to a prompt fluorescence. As the feasibility of S₁ ↔ T₁ conversion increases in the yellow dimer, it allows a delay in fluorescence emission, leading to thermally activated delayed fluorescence (TADF).

Dipole moments of conjugated donor-acceptor substituted systems: calculations vs. experiments

We find that quantum mechanical calculations using B3LYP/aug-cc-pVTZ model chemistry involving anharmonic correction on simple conjugated organic compounds without rotating moieties provide the dipole moment values and molecular geometries with high accuracy. In the presence of one or two conjugated electron donating or accepting substituents capable of hindered rotation, the calculated dipole moments reproduce the experimental results equally well only in the cases when the experiments were done at the temperatures at which rotation of substituents remains hindered. In order to reproduce the experimental dipole moments determined at higher temperatures, a model assuming free (unhindered) rotation should be applied. In these cases, the contribution of each rotamer is equal and using anharmonic correction is not necessary. The APFD functional produces similar results and the M062X functional yields larger deviations from the experimental data. The other methods, like HF and MP2, are the least accurate with the basis sets usually employed for interpreting the experimental data.

B3LYP functional can reproduce the experimental dipole moments of conjugated organic molecules with high accuracy.

Singlet fission relevant energetics from optimally tuned range-separated hybrids

As a promising idea to design high-efficiency organic photovoltaics, singlet fission (SF) mechanism, i.e., generating two triplet excitons out of a single photon absorption, has recently come into the spotlight. Even though much effort has been devoted to this arena, accurately accounting for the SF process from the theoretical perspective has proven to be challenging. Herein, the SF energetics have thoroughly been investigated with the help of optimally tuned range-separated hybrid functionals (OT-RSHs) in both gas and solvent phases. Taking a series of experimentally known SF chromophores as working models, we have proposed and validated several variants of OT-RSH approximations for the reliable prediction of the energy levels which match the crucial criteria for the SF process, namely, the negative singlet-triplet and triplet-triplet energy gaps. We scrutinize the role of the OT-RSH ingredients, i.e., the underlying density functional approximations, short- and long-range exact-like exchange, as well as the range-separation parameter, for our purpose. The newly designed OT-RSHs outperform the standard RSHs and other related schemes such as screened-exchange approximations as well as other density functionals from different rungs for describing the SF energetics. More importantly, it is unveiled that although the OT-RSH coupled with the polarizable continuum model, OT-RSH-PCM, as well as the screened versions, OT-SRSHs, which account for the screening effect by the electron correlation through the scalar dielectric constant have some advantages over gas-phase computations using OT-RSHs, the energetics criteria of the SF process may not necessarily be satisfied. This in turn corroborates the idea of performing both the optimal tuning procedure and subsequent computations of the SF relevant energetics using OT-RSHs as a more reliable and affordable framework, at least for the present purpose. The applicability of the proposed models is also put into broader perspective, where they are used for the computational design of several chromophores as promising candidates prone to utilization in the SF-based materials. Hopefully, our recommended OT-RSHs can function as efficient models for both the theoretical modeling of SF chromophores and confirming the experimental observations in the field.

Unveiling the role of short-range exact-like exchange in the optimally tuned range-separated hybrids for fluorescence lifetime modeling

We propose and validate several variants of the optimally tuned range-separated hybrid functionals (OT-RSHs) including different density functional approximations for predicting the fluorescence lifetimes of different categories of fluorophores within the time-dependent density functional theory (TD-DFT) framework using both the polarizable continuum and state-specific solvation models. Our main idea originates from performing the optimal tuning in the presence of a contribution of the exact-like exchange at the short-range part, which, in turn, leads to the small values of the range-separation parameter, and computing the fluorescence lifetimes using the models including no or small portions of the short-range exact-like exchange. Particular attention is also paid to the influence of the geometries of emitters on fluorescence lifetime computations. It is shown that our developed OT-RSHs along with the polarizable continuum model can be considered as the promising candidates within the TD-DFT framework for the prediction of fluorescence lifetimes for various fluorophores. We find that the proposed models not only outperform their standard counterparts but also provide reliable data better than or comparable to the conventional hybrid functionals with both the fixed and interelectronic distance-dependent exact-like exchanges. Furthermore, it is also revealed that when the excited state geometries come into play, more accurate descriptions of the fluorescence lifetimes can be achieved. Hopefully, our findings can give impetus for future developments of OT-RSHs for computational modeling of other characteristics in fluorescence spectroscopy as well as for verification of the related experimental observations.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Electronic structure and second-order nonlinear optical properties of lemniscular [16]cycloparaphenylene compounds

Chiral organic compounds are excellent second-order nonlinear optical (NLO) materials due to their inherent non-symmetric electronic structures combined with the advantages of organic compounds. At present, density functional theory (DFT) has become a powerful tool for predicting the properties of novel materials. In this paper, based on chiral lemniscular [16]cycloparaphenylene, three novel compounds are designed by introduction of donor/acceptor units and their combinations. The geometrical/electronic structure, electronic absorption, and the second-order NLO properties of these compounds have been systematically investigated by DFT/TDDFT theory. The simulated UV-Vis/CD spectra of compound 1 are in good agreement with the experimental ones, enabling us to assign their electronic transition characteristics and absolute configuration with high confidence. The investigations show that energy gaps, absorption wavelength and second-order NLO response may be effectively tuned by the introduction of the donor or acceptor units or their combinations. For instance, the second-order NLO value of compound 4 is about 207 times as large as the average second-order polarizability of the organic molecule urea. Thus, the studied compounds are expected to be potential large second-order NLO materials.

Performance of DFT functionals for calculating the second-order nonlinear optical properties of dipolar merocyanines

Evolution of the static HRS hyperpolarizability of a tricyanopropylidene-based merocyanine dye with the length of the polyenic bridge, as calculated using various ab initio and DFT approximations.

Large-scale comparison of 3d and 4d transition metal complexes illuminates the reduced effect of exchange on second-row spin-state energetics

The origin of distinct 3d vs. 4d transition metal complex sensitivity to exchange is explored over a large data set.

A new tuned range-separated density functional for the accurate calculation of second hyperpolarizabilities

Among the nine functionals benchmarked, the most accurate γ are obtained by Tα-LC-BLYP, reducing about half the errors of LC-BLYP.

Toward photophysical characteristics of triplet-triplet annihilation photon upconversion: a promising protocol from the perspective of optimally tuned range-separated hybrids

The photon upconversion (UC) process assisted by the triplet-triplet annihilation (TTA) mechanism has recently come into the spotlight. Given the rich collection of efforts in this area, theoretical explorations regarding TTA-UC are relatively limited and have proven to be challenging for its control in devices. In this contribution, the photophysical properties crucial for TTA-UC, such as triplet excited state energies and triplet-triplet energy transfer gaps of the essential ingredients involved in the process, namely sensitizers, annihilators and their pairs, have theoretically been investigated using optimally tuned range-separated hybrid functionals (OT-RSHs) and their screened exchange counterparts, OT-SRSHs. Taking a series of experimentally proven-to-work sensitizer/annihilator pairs as working models, we have constructed and validated several variants of OT-RSHs using both full time-dependent and Tamm-Dancoff formalisms for a reliable description of the TTA-UC photophysics. Given the bimolecular biphotonic nature of the TTA-UC process under study, particular attention is paid to the influence of the factors like the underlying density functional approximations and the tunable parameters such as short- and long-range exact-like exchanges as well as the range-separation parameter for both the sensitizers and annihilators separately. Dissecting all the aspects and relying on the appropriate choices from the tested models, we propose an OT-RSH with the correct asymptotic behavior as a cost-effective yet useful tool for this purpose. Not only against the standard RSHs but also in comparison to the conventional hybrids, the newly developed OT-RSH yields a more reliable description for the TTA-UC energetics in the gas phase and dielectric medium. Accountability of the proposed model has further been confirmed for several theoretically designed sensitizer/annihilator pairs prone to be used in the TTA-UC process. Summing up, in light of this study additional pieces of convincing evidence on the quality of OT-(S)RSHs for computational modeling and experimental verifications of the photophysics of the photon UC based on TTA and other possible technologies are showcased.

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.

Theoretical study on the electronic structure and second-order nonlinear optical properties of benzannulated or selenophene-annulated expanded helicenes

Currently, discovering new materials with superior second-order nonlinear optical (NLO) performance has become a very hot research topic in the fields of chemistry and materials science. Now, density functional theory (DFT) has become a powerful tool to predict the performance of novel materials. In this paper, based on benzannulated and selenophene-annulated expanded helicenes, twenty-six helicenes are designed by introduction donor/acceptor moieties and their combinations at different substituent positions. The geometrical/electronic structures, electronic transition, and second-order NLO properties of these helicenes are full investigated by DFT/TDDFT theory. The investigations show that these helicenes have large first hyperpolarizability values (β HRS). For instance, the β HRS value (29.95 × 10-30 esu) of helicene H24 is about 7 times larger than that of the highly π-delocalized phenyliminomethyl ferrocene complex. In addition, the introduction of acceptor NO2 unit at R7 and R8 positions for helicenes H1 and H15 can obtain the largest β HRS value, which is attributed to the enhancement of electron acceptor ability. In view of large NLO response and intrinsic asymmetric structures, the studied helicenes have the possibility to be excellent second-order NLO materials.

Structural diradical character

A reliable first-principles description of singlet diradical character is essential for predicting nonlinear optical and magnetic properties of molecules. As diradical and closed-shell electronic structures differ in their distribution of single, double, triple, and aromatic bonds, modeling electronic diradical character requires accurate bond-length patterns, in addition to accurate absolute bond lengths. We therefore introduce structural diradical character, which we suggest as an additional measure for comparing first-principles calculations with experimental data. We employ this measure to identify suitable exchange-correlation functionals for predicting the bond length patterns and electronic diradical character of a biscobaltocene with the potential for photoswitchable nonlinear optical activity. Out of four popular approximate exchange-correlation functionals with different exact-exchange admixtures (BP86, TPSS, B3LYP, TPSSh), the two hybrid functionals TPSSh and B3LYP perform best for diradical bond length patterns, with TPSSh being best for the organometallic validation systems and B3LYP for the organic ones (for which the D3 dispersion correction was included). Still, none of the functionals is suitable for correctly describing relative bond lengths across the range of molecules studied, so that none can be recommended for predictive studies of (potential) diradicals without reservation. © 2018 Wiley Periodicals, Inc.

Photo-tunable linear and nonlinear optical response of cyclophanediene-dihydropyrene photoswitches

Cyclophanediene (CPD)-dihydropyrene (DHP) is a negative T-type photochrome pair having a thermodynamically stable colored form, i.e., DHP. Interconversion between cyclophanediene and dihydropyrene is associated with significant changes in dipole moment, absorption wavelength and polarizability, which can impart substantial linear and nonlinear optical response. In this study, phototunable linear and nonlinear optical response of cyclophanediene-dihdyropyrene photoswitches is described. Cyclophanedienes and dihydropyrenes are functionalized at the internal position for maximum changes in volume and polarizability. The UV-Vis spectra are calculated at ɷB97XD, which was validated through a benchmark approach. An excellent correlation is observed between theoretical and experimental absorption spectra. Several CPD-DHP pairs have been recognized for clean interconversion in UV-Vis light without formation of a photostationary state. Nonlinear optical response of dihydropyrenes is remarkably higher than that of cyclophanedienes. In general, the calculated hyperpolarizability values of dihydropyrenes are about two to three orders of magnitude higher than those for cyclophanedienes. The trends in calculated hyperpolarizabilities are rationalized through two level method. The high nonlinear optical response of dihydropyrenes stems from low excitation energies. The remarkable difference in hyperpolarizabilities of these isomeric forms paves path for the design of phototunable nonlinear optical materials.

Transport and Optical Gaps in Amorphous Organic Molecular Materials

The standard procedure to identify the hole- or electron-acceptor character of amorphous organic materials used in OLEDs is to look at the values of a pair of basic parameters, namely, the ionization potential (IP) and the electron affinity (EA). Recently, using published experimental data, the present authors showed that only IP matters, i.e., materials with IP > 5.7 (<5.7) showing electron (hole) acceptor character. Only three materials fail to obey this rule. This work reports ab initio calculations of IP and EA of those materials plus two materials that behave according to that rule, following a route which describes the organic material by means of a single molecule embedded in a polarizable continuum medium (PCM) characterized by a dielectric constant ε . PCM allows to approximately describe the extended character of the system. This "compound" system was treated within density functional theory (DFT) using several combinations of the functional/basis set. In the preset work ε was derived by assuming Koopmans' theorem to hold. Optimal ε values are in the range 4.4⁻5.0, close to what is expected for this material family. It was assumed that the optical gap corresponds to the excited state with a large oscillator strength among those with the lowest energies, calculated with time-dependent DFT. Calculated exciton energies were in the range 0.76⁻1.06 eV, and optical gaps varied from 3.37 up to 4.50 eV. The results are compared with experimental data.

Dielectric Screening Meets Optimally Tuned Density Functionals

A short overview of recent attempts at merging two independently developed methods is presented. These are the optimal tuning of a range-separated hybrid (OT-RSH) functional, developed to provide an accurate first-principles description of the electronic structure and optical properties of gas-phase molecules, and the polarizable continuum model (PCM), developed to provide an approximate but computationally tractable description of a solvent in terms of an effective dielectric medium. After a brief overview of the OT-RSH approach, its combination with the PCM as a potentially accurate yet low-cost approach to the study of molecular assemblies and solids, particularly in the context of photocatalysis and photovoltaics, is discussed. First, solvated molecules are considered, with an emphasis on the challenge of balancing eigenvalue and total energy trends. Then, it is shown that the same merging of methods can also be used to study the electronic and optical properties of molecular solids, with a similar discussion of the pros and cons. Tuning of the effective scalar dielectric constant as one recent approach that mitigates some of the difficulties in merging the two approaches is considered.