Tuned Range-Separated Time-Dependent Density Functional Theory Applied to Optical Rotation

Simplified Tuning of Long-Range Corrected Time-Dependent Density Functional Theory

Range-separated hybrid functionals have dramatically improved the description of charge-transfer excitations in time-dependent density functional theory (TD-DFT), especially when the range-separation parameter is adjusted in order to satisfy the ionization energy (IE) criterion, εHOMO = -IE. However, this "optimal tuning" procedure is molecule-specific, inconvenient, expensive for large systems, and problematic in extended or periodic systems. Here, we consider an alternative procedure known as global density-dependent (GDD) tuning, which sets the range-separation parameter in an automated way based on properties of the exchange hole. In small molecules, we find that long-range corrected functionals with either IE or GDD tuning afford remarkably similar TD-DFT excitation energies, for both valence and charge-transfer excitations. However, GDD tuning is more efficient and is well-behaved even for large systems. It provides a black-box solution to the optimal-tuning problem that can replace IE tuning for many applications of TD-DFT.

Conformational and Solvent Effects on the Photoinduced Electron Transfer Dynamics of a Zinc Phthalocyanine-Benzoperylenetriimide Conjugate: A Nonadiabatic Dynamics Simulation

Herein, we employed a combination of static electronic structure calculations and nonadiabatic dynamics simulations at linear-response time dependent density functional theory (LR-TDDFT) level with the optimally tuned range-separated hybrid (OT-RSH) functional to explore the ultrafast photoinduced dynamics of a zinc phthalocyanine-benzoperylenetriimide (ZnPc-BPTI) conjugate. Due to the flexibility of the linker, we identified two major conformations: the stacked conformation (ZnPc-BPTI-1) and the extended conformation (ZnPc-BPTI-2). Since the charge transfer states are much lower than the lowest local excitation in ZnPc-BPTI-1, which is contrary to ZnPc-BPTI-2, the ultrafast electron transfer (~3.6 ps) is only observed in the nonadiabatic simulations of ZnPc-BPTI-1 upon local excitation around the absorption maximum of ZnPc. However, when considering the solvent effects in benzonitrile: the lowest S1 states are both charge transfer states from ZnPc to BPTI for different conformers. Subsequent nonadiabatic dynamics simulations indicate that both conformers experience ultrafast electron transfer in benzonitrile with two time constants of 90 [100] fs and 1.40 [1.43] ps. Our present work not only agrees well with previous experimental study, but also points out the important role of conformational changes and solvent effects in regulating the photodynamics of organic donor-acceptor conjugates.

Impact of doping with organic dopants and mixed doping with alkali metals and organic dopants on the absorption, electronic, optoelectronic, thermodynamic and nonlinear optical properties of dibenzo[b,def]chrysene in gaseous media: DFT and TD-DFT studies

CONTEXT

In this study, we evaluate the geometrical, absorption, optoelectronic, electronic, nonlinear optical (NLO) and thermodynamic properties of dibenzo[b,def]chrysene molecule derivatives by means of DFT and TD-DFT simulations. In view of the aim of producing new high-performance materials for non-linear optics (NLO) by doping test, two types of doping were used. We obtained six derivatives by doping with organic dopants (Nitro, amide and ticyanoethenyl) and mixed alkali metal (potassium) and organic dopants. Doping with organic dopants produced molecules A, B and C, respectively when substituting one hydrogen with nitro (NO2), amide (CONH2) and tricyanoethenyl (C5N3) groups, while mixed doping involved considering A, B and C and then substituting two hydrogens with two potassiums to obtain compounds D, E and F respectively. The negative values of the various interaction energies calculated for all the doped molecules show that they are all stable, but also that molecules C and F are the most stable in the case of both dopings. The gap energies calculated at the B3LYP level of theory are all below 3 eV, which means that all the molecules obtained are semiconductors. Better still, compounds C and F, with gap energies of 1.852 eV and 1.204 eV, respectively, corresponding to decreases of 35.67% and 58.18% in gap energy compared with the pristine molecule, are more reactive than the other doped molecules. Mixed doping is therefore a highly effective way of narrowing the energy gap and boosting the semiconducting character and reactivity of organic materials. Optoelectronic properties have also been improved, with refractive index values higher than those of the reference material, glass. This shows that our compounds could be used under very high electric field conditions of the order of 4.164 × 109 V.m-1 for C and 7.410 × 109 V.m-1 for F the highest values at the B3LYP level of theory. The maximum first-order hyperpolarizability values for both types of doping are obtained at the CAM-B3LYP level of theory by C:β mol = 92.088 × 10-30esu and by F:β mol = 129.449 × 10-30esu, and second-order values are also given by these same compounds. These values are higher than the reference value, which is urea, making our compounds potential candidates for high-performance NLO applications. In dynamic mode and at a frequency of 1064 nm, at the CAM-B3LYP level of theory, the highest dynamic hyperpolarizability coefficients were obtained by C and F. Hyper-Rayleigh scattering β HRS , coefficients of the electro-optical Pockel effect (EOPE), EFISHG, third-order NLO-response degree four-wave mixingγ DFWM , quadratic nonlinear refractive index n2 were also calculated. The maximum values of n2 are obtained by C (6.13 × 10-20 m2/W) and F (6.60 × 10-20 m2/W), these values are 2.24 times higher than that of fused silica which is the reference for degenerate four-wave mixing so our molecules could also have applications in optoelectronics as wavelength converters, optical pulse modulators and optical switches.

METHODS

Using the DFT method, we were able to determine the optimized and stable electronic structures of doped dibenzo[b,def]chrysene derivatives in the gas phase. We limited ourselves to using the proven B3LYP and CAMB3LYP levels of theory for calculating electronic properties, and non-linear optics with the 6-311G +  + (d,p) basis set, which is a large basis set frequently used for these types of compound. Gaussian 09 software was used to run our calculations, and Gauss View 6.0.16 was used to visualize the output files. TD-DFT was also used to determine absorption properties at the B3LYP level of theory, using the same basis set.

Understanding Photophysical Properties of Molecules Relevant in Organic Semiconductor Laser Diodes from Electron Localization Function-Tuned and Solvent-Tuned Range-Separated Functionals

Organic semiconductor laser diodes (OSLDs) are prevalent in optoelectronics because of their sustainable energy applications. Organic molecules used in such diodes are usually large; hence, their studies are computationally challenging with high-end benchmark methods. Computational methods with reliable accuracy and efficiency are always indispensable. In the present work, we have applied our computationally inexpensive, nonempirically tuned [electron localization function (ELF*) and solvent (Sol*)] range-separated (RS) functionals to study five molecules used in OSLDs. The emission energies in three different environments [toluene, CBP (4,4'-bis(n-carbazolyl)-1,1'-biphenyl) film, and gas] have been computed with the tuned functionals and compared with the experimental emission energies. ELF* and Sol* functionals can accurately reproduce emission energies in toluene and CBP film environments. On the other hand, both ELF* and IP-tuned functionals with excited-state geometry (IP*) perform better in the gas phase. In addition, a comparative study is performed between time-dependent density functional theory and the Tamm-Dancoff approximation. Along with the emission energy, oscillator strength values have also been reported. Different IP-tuned RS parameters were obtained with the ground- and excited-state geometries. Interestingly, it has been observed that the optimally tuned RS parameter with excited-state geometry (IP*) performs better compared to that with ground-state geometries (IP). Fractional occupation calculations show that the tuned functionals exhibit less localization and delocalization error. The study envisages that ELF* and Sol* functionals can be used to design future candidates for OSLDs.

Solvent effects on the photoinduced charge separation dynamics of directly linked zinc phthalocyanine-perylenediimide dyads: a nonadiabatic dynamics simulation with an optimally tuned screened range-separated hybrid functional

Herein, we have employed a combination of the optimally tuned screened range-separated hybrid (OT-SRSH) functional, the polarizable continuum model (PCM), and nonadiabatic dynamics (NAMD) simulations to investigate the photoinduced dynamics of directly linked donor-acceptor dyads formed using zinc phthalocyanine (ZnPc) and perylenediimide (PDI), in which ZnPc is the donor while PDI is the acceptor. Our simulations aim to analyze the behavior of these dyads upon local excitation of the ZnPc moiety in the gas phase and in benzonitrile. Our findings indicate that the presence of a solvent can significantly influence the excited state dynamics of ZnPc-PDI dyads. Specifically, the polar solvent benzonitrile effectively lowers the vertical excitation energies of the charge transfer (CT) state from ZnPc to PDI. As a result, the energetic order of the locally excited (LE) states of ZnPc and the CT states is reversed compared to the gas phase. Consequently, the photoinduced electron transfer (PET) dynamics from ZnPc to PDI, which is absent in the gas phase, takes place in benzonitrile with a time constant of 10.4 ps. Importantly, our present work not only qualitatively agrees with experimental results but also provides in-depth insights into the underlying mechanisms responsible for the photoinduced dynamics of ZnPc-PDI. Moreover, this study emphasizes the importance of appropriately considering solvent effects in NAMD simulation of organic donor-acceptor systems, taking into account the distinct excited state dynamics observed in the gas phase and benzonitrile. Furthermore, the combination of the OT-SRSH functional, the PCM solvent model, and nonadiabatic dynamics simulations shows promise as a strategy for investigating the complex excited state dynamics of organic donor-acceptor systems in solvents. These findings will be valuable for the future design of novel organic donor-acceptor structures with improved performance.

On the choice of coordinate origin in length gauge optical rotation calculations

In this work, we explore the issue of origin dependence in optical rotation (OR) calculations in the length dipole gauge (LG) using standard approximate methods belonging to density functional theory (DFT) and coupled cluster (CC) theory. We use the origin-invariant LG approach, LG(OI), that we recently proposed as reference for the calculations, and we study whether a proper choice of coordinate origin and molecular orientation can be made such that diagonal elements of the LG-OR tensor match those of the LG(OI) tensor. Using a numerical search algorithm, we show that multiple spatial orientations can be found where the LG and LG(OI) results match. However, a simple analytical procedure provides a spatial orientation where the origin of the coordinate system is close to the center of mass of the molecule. At the same time, we also show that putting the origin at the center of mass is not an ideal choice for every molecule (relative errors in the OR up to 70% can be obtained in out test set). Finally, we show that the choice of coordinate origin based on the analytical procedure is transferable across different methods and it is superior to putting the origin in the center of mass or center of nuclear charge. This is important because the LG(OI) approach is trivial to implement for DFT, but not necessarily for nonvariational methods in the CC family. Therefore, one can determine an optimal coordinate origin at DFT level and use it for standard LG-CC response calculations.

X-ray and Optical Circular Dichroism as Local and Global Ultrafast Chiral Probes of [12]Helicene Racemization

Chirality is a fundamental molecular property that plays a crucial role in biophysics and drug design. Optical circular dichroism (OCD) is a well-established chiral spectroscopic probe in the UV-visible regime. Chirality is most commonly associated with a localized chiral center. However, some compounds such as helicenes (Figure 1) are chiral due to their screwlike global structure. In these highly conjugated systems, some electric and magnetic allowed transitions are distributed across the entire molecule, and OCD thus probes the global molecular chirality. Recent advances in X-ray sources, in particular the control of their polarization and spatial profiles, have enabled X-ray circular dichroism (XCD), which, in contrast to OCD, can exploit the localized and element-specific nature of X-ray electronic transitions. XCD therefore is more sensitive to local structures, and the chirality probed with it can be referred to as local. During the racemization of helicene, between opposite helical structures, the screw handedness can flip locally, making the molecule globally achiral while retaining a local handedness. Here, we use the racemization mechanism of [12]helicene as a model to demonstrate the capabilities of OCD and XCD as time-dependent probes for global and local chiralities, respectively. Our simulations demonstrate that XCD provides an excellent spectroscopic probe for the time-dependent local chirality of molecules.

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.

Accurate Non-Adiabatic Couplings from Optimally-Tuned Range-Separated Hybrid Functionals

Precise theoretical calculations of non-adiabatic couplings, which describe the interaction between two Born-Oppenheimer surfaces, are important for the modeling of radiationless decay mechanisms in photochemical processes. Here, we demonstrate that accurate non-adiabatic couplings can be calculated in the framework of linear-response time-dependent density functional theory by using non-empirical, optimally-tuned range-separated hybrid (OT-RSH) functionals. We focus on molecular radicals, in which ultrafast non-radiative decay plays a crucial role, to find that the OT-RSH functional compares well to wave-function based reference data and competes with the accuracy of semi-empirical CAM-B3LYP calculations. Our findings show that the OT-RSH approach provides very accurate non-adiabatic couplings and, therefore, provides a computationally efficient alternative to wave-function based techniques.

Excited-State Properties of Some Thermally Activated Delayed Fluorescence Emitters: Quest for an Accurate and Reliable Computational Method

Thermally activated delayed fluorescence (TADF) finds application in organic light-emitting diodes. The molecules exhibiting TADF are characterized by small singlet-triplet energy gaps that help reverse intersystem crossing. Recently, ionization potential (IP)-tuned range-separated (RS) density functionals have been well accepted for studying excited-state properties. In the present work, two efficient descriptor-based tuning schemes [electron localization function (ELF) and Sol] of RS density functionals have been used to accurately reproduce the excited-state properties of TADF emitters by performing a single self-consistent field calculation. The lowest singlet vertical excitation energies (E_VA(S₁)) and the vertical singlet-triplet energy gaps (ΔE_VST) are computed with ELF-, Sol-, and IP-tuned RS functionals (LC-BLYP, ωB97, ωB97X, and ωB97XD). Encouraging mean absolute deviations from the experimental values with ELF*-, Sol*-, and IP-tuned functionals are observed. Consistent performance of the non-empirical tuned functionals is noted in different solvent dielectrics. In addition to these, fractional occupation calculations have shown that our tuned functionals almost satisfy the energy linearity curve. Thus, ELF*- and Sol*-tuned functionals are promising and reliable alternatives in computing the excited-state properties. Considering the small experimental singlet-triplet gap, we recommend ELF* to calculate E_VA(S₁) and Sol* to calculate ΔE_VST.

Finding reliable methodology for optical rotation and correct predictions of (s)-methyloxirane and (1R,5R)-β-pinene

Theoretical calculations of optical rotation (OR), although important to predict absolute configurations (ACs) and corroborate experiments, require efficient methodology able to reproduce enantiomer specificity and real OR values. Also, troublesome molecules are recurring in the literature, such as (S)-methyloxirane and (1R,5R)-β-pinene. This study evaluates DFT functionals B3LYP, CAM-B3LYP, ωB97X-D, M06-2X, and PBE0 considering basis sets aug-cc-pVDZ, aug-cc-pVTZ, 6-311++G(2d,p), and 6-311++G(3df,2p) in OR prediction of 42 rigid organic molecules assessing cases with wrong enantiomeric determination comparing to available experimental data at wavelengths 355, 589, and 633 nm. Functionals CAM-B3LYP and ωB97X-D with aug-cc-pVTZ are indicated here to reproduce experimental values more accurately considering fewer number of wrong AC predictions, normalized RMSD values below 0.70, and a good approximation to experimental values in hierarchical cluster analysis. Methyloxirane AC was reproduced in CAM-B3LYP and PBE0, with [ α ] 355 = 6 . 94 $$ {\left[\alpha \right]}_{355}=6.94 $$ for CAM-B3LYP/aug-cc-pVTZ close to experimental value [ α ] 355 = 7 . 49 ± 0 . 30 $$ {\left[\alpha \right]}_{355}=7.49\pm 0.30 $$ . Good results were found for AC of β-pinene in M06-2X, CAM-B3LYP, and ωB97X-D while the latter in 6-311++G(3df,2p) obtained OR values of [ α ] 589 = 3 . 44 $$ {\left[\alpha \right]}_{589}=3.44 $$ and [ α ] 689 = 4 . 20 $$ {\left[\alpha \right]}_{689}=4.20 $$ close to experimental values [ α ] 589 = 2 . 8 $$ {\left[\alpha \right]}_{589}=2.8 $$ and [ α ] 689 = 4 . 66 ± 0 . 60 $$ {\left[\alpha \right]}_{689}=4.66\pm 0.60 $$ . The two molecules aforementioned are, for the first time, reported to give valid theoretical OR values in such simple methodologies. OR calculations were all performed after geometry optimization at the same level of theory, and analysis of different functional combinations for each step in β-pinene showed it can interfere with AC prediction even in rigid molecules.

Comparative Study of Vibrational Raman Optical Activity with Different Time-Dependent Density Functional Approximations: The VROA36 Database

A new database, VROA36, is introduced to investigate the performance of computational approaches for vibrational Raman optical activity (VROA) calculations. The database is composed of 36 molecules with known experimental VROA spectra. It includes 93 conformers. Normal modes calculated with B3LYP-D3(BJ)/def2-TZVP are used to compute the VROA spectra with four functionals, B3LYP-D3(BJ), ωB97X-D, M11, and optimally tuned LC-PBE, as well as several basis sets. SimROA indices and frequency scaling factors are used to compare calculated spectra with each other and with experimental data. The four functionals perform equally well independently of the basis set and usually achieve good agreement with the experimental data. For molecules in near- or at-resonance conditions, the inclusion of a complex (damped) linear response approach is important to obtain physically meaningful VROA intensities. The use of any of the tested functional approximations with the def2-SVPD Gaussian-type basis set, or a basis of similar flexibility, can be recommended for efficient and reliable theoretical VROA studies.

Basis Set Dependence of Optical Rotation Calculations with Different Choices of Gauge

In this work, the basis set dependence of optical rotation (OR) calculations is examined for various choices of gauge/level of theory. The OR is calculated for a set of 50 molecules using B3LYP and CAM-B3LYP and 17 molecules using coupled cluster with single and double excitations (CCSD). The calculations employ the correlation-consistent basis sets, aug-cc-pVζZ with ζ = D, T, Q. An inverse-power extrapolation formula is then utilized to obtain OR values at the complete basis set (CBS) limit. We investigate the basis set convergence for these methods and three choices of gauge: length gauge (with gauge-including atomic orbitals, LG(GIAOs), for DFT), the origin-invariant length gauge [LG(OI)], and the modified velocity gauge (MVG). The results show that all methods converge smoothly to the CBS limit and that the LG(OI) approach has a slightly faster convergence rate than the other choices of gauge. While the DFT methods reach gauge invariance at the CBS limit, CCSD does not. The significant difference between the MVG and LG(OI) results at the CBS limit, 26%, indicates that CCSD is not quite at convergence in the description of electron correlation for this property. On the other hand, gauge invariance at the CBS limit for DFT does not lead to the same OR values for the two density functionals, which is also due to electron correlation incompleteness. A limited comparison to gas-phase experimental OR values for the DFT methods shows that CAM-B3LYP seems more accurate than B3LYP. Overall, this study shows that the LG(OI) approach with the aug-cc-pVTZ basis set for DFT, and with the CBS(DT) extrapolation for CCSD, provides a good cost/accuracy balance.

Automated Quantum Chemistry-Based Calculation of Optical Rotation for Large Flexible Molecules

The calculation of optical rotation (OR, [α]_D) for nonrigid molecules was limited to small systems due to the challenging problem of generating reliable conformer ensembles, calculating accurate Boltzmann populations and the extreme sensitivity of the OR to the molecules' three-dimensional structure. Herein, we describe and release the crenso workflow for the automated computation of conformer ensembles in solution and corresponding [α]_D values for flexible molecules. A comprehensive set of 28 organic drug molecules (28-144 atoms) with experimentally determined values is used in our assessment. In all cases, the correct OR sign is obtained with an overall mean relative deviation of 72% (mean absolute deviation of 82 °[dm(g/cm³)]^-1 for experimental values in the range -160 to 287 °[dm(g/cm³)]^-1). We show that routine [α]_D computations for very flexible, biologically active molecules are both feasible and reproducible in about a day of computation time on a standard workstation computer. Furthermore, we observed that the effect of energetically higher-lying structures in the ensemble on the OR is often averaged out and that in 23 out of 28 cases, the correct OR sign is obtained by just considering only the lowest free energy conformer. In four example cases, we show that the approach can also describe the OR of pairs of flexible diastereomers properly. In summary, even very sensitive, multifactorial physicochemical properties appear reliably predictable with minimal user input from efficiently automated quantum chemical methods.

Valence-shell photoelectron circular dichroism of ruthenium(III)-tris-(acetylacetonato) gas-phase enantiomers

Chiral transition-metal complexes are of interest in many fields ranging from asymmetric catalysis and molecular materials science to optoelectronic applications or fundamental physics including parity violation effects. We present here a combined theoretical and experimental investigation of gas-phase valence-shell photoelectron circular dichroism (PECD) on the challenging open-shell ruthenium(III)-tris-(acetylacetonato) complex, Ru(acac)₃. Enantiomerically pure Δ- or Λ-Ru(acac)₃, characterized by electronic circular dichroism (ECD), were vaporized and adiabatically expanded to produce a supersonic beam and photoionized by circularly-polarized VUV light from the DESIRS beamline at Synchrotron SOLEIL. Photoelectron spectroscopy (PES) and PECD experiments were conducted using a double imaging electron/ion coincidence spectrometer, and compared to density functional theory (DFT) and time-dependent DFT (TDDFT) calculations. The open-shell character of Ru(acac)₃, which is not taken into account in our DFT approach, is expected to give rise to a wide multiplet structure, which is not resolved in our PES signals but whose presence might be inferred from the additional striking features observed in the PECD curves. Nevertheless, the DFT-based assignment of the electronic bands leads to the characterisation of the ionized orbitals. In line with other recent works, the results confirm that PECD persists independently on the localization and/or on the achiral or chiral nature of the initial orbital, but is rather a probe of the molecular potential as a whole. Overall, the measured PECD signals on Ru(acac)₃, a system exhibiting D₃ propeller-type chirality, are of similar magnitude compared to those on asymmetric-carbon-based chiral organic molecules which constitute the vast majority of species investigated so far, thus suggesting that PECD is a universal mechanism, inherent to any type of chirality.

A Density Functional Theory Study of Optical Rotation in Some Aziridine and Oxirane Derivatives

We present time-dependent density functional theory (TDDFT) calculations of the electronic optical rotation (ORP) for seven oxirane and two aziridine derivatives in the gas phase and in solution and compare the results with the available experimental values. For seven of the studied molecules it is the first time that their optical rotation was studied theoretically and we have therefore investigated the influence of several settings in the TDDFT calculations on the results. This includes the choice of the one-electron basis set, the exchange-correlation functional or the particular polarizable continuum model (PCM). We can confirm that polarized quadruple zeta basis sets augmented with diffuse functions are necessary for converged results and find that the aug-pc-3 basis set is a viable alternative to the frequently employed aug-cc-pVQZ basis set. Based on our study, we cannot recommend the generalized gradient functional KT3 for calculations of the ORP in these compounds, whereas the hybrid functional PBE0 gives results quite similar to the long-range correct CAM-B3LYP functional. Finally, we observe large differences in the solvent effects predicted by the integral equation formalism of PCM and the SMD variant of PCM. For the majority of solute/solvent combinations in this study, we find that the SMD model in combination with the PBE0 functional and the aug-pc-3 basis set gives the best agreement with the experimental values.

Full optical rotation tensor at coupled cluster with single and double excitations level in the modified velocity gauge

This work presents the first simulations of the full optical rotation (OR) tensor at coupled cluster with single and double excitations (CCSD) level in the modified velocity gauge (MVG) formalism. The CCSD-MVG OR tensor is origin independent, and each tensor element can in principle be related directly to experimental measurements on oriented systems. We compare the CCSD results with those from two density functionals, B3LYP and CAM-B3LYP, on a test set of 22 chiral molecules. The results show that the functionals consistently overestimate the CCSD results for the individual tensor components and for the trace (which is related to the isotropic OR), by 10%-20% with CAM-B3LYP and 20%-30% with B3LYP. The data show that the contribution of the electric dipole-magnetic dipole polarizability tensor to the OR tensor is on average twice as large as that of the electric dipole-electric quadrupole polarizability tensor. The difficult case of (1S,4S)-(-)-norbornenone also reveals that the evaluation of the former polarizability tensor is more sensitive than the latter. We attribute the better agreement of CAM-B3LYP with CCSD to the ability of this functional to better reproduce electron delocalization compared with B3LYP, consistent with previous reports on isotropic OR. The CCSD-MVG approach allows the computation of reference data of the full OR tensor, which may be used to test more computationally efficient approximate methods that can be employed to study realistic models of optically active materials.

Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional Theory

The interplay of electronic excitations and structural changes in molecules impacts nonradiative decay and charge transfer in the excited state, thus influencing excited-state lifetimes and photocatalytic reaction rates in optoelectronic and energy devices. To capture such effects requires computational methods providing an accurate description of excited-state potential energy surfaces and geometries. We suggest time-dependent density functional theory using optimally tuned range-separated hybrid (OT-RSH) functionals as an accurate approach to obtain excited-state molecular geometries. We show that OT-RSH provides accurate molecular geometries in excited-state potential energy surfaces that are complex and involve an interplay of local and charge-transfer excitations, for which conventional semilocal and hybrid functionals fail. At the same time, the nonempirical OT-RSH approach maintains the high accuracy of parametrized functionals (e.g., B3LYP) for predicting excited-state geometries of small organic molecules showing valence excited states.

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).

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.

Approximate versus Exact Embedding for Chiroptical Properties: Reconsidering Failures in Potential and Response

We investigate the suitability of subsystem time-dependent density-functional theory (sTDDFT) for describing chiroptical properties with a focus on optical rotation parameters. Our starting point is a new implementation of the recently proposed projection-based, coupled frozen-density embedding (FDEc) framework. We adapt the generalized, non-Hermitian formulation of TDDFT and derive corresponding expressions for regular and damped response properties from subsystem TDDFT. We verify that our implementation of this "exact" formulation allows to reproduce supermolecular results of electronic circular dichroism (ECD) spectra, of optical rotatory dispersion, and of polarizabilities. We present a systematic test of the main approximations typically introduced in practical frozen-density embedding (FDE) calculations of response properties: (i) the use of approximate nonadditive kinetic-energy (NAKE) functionals, which can be avoided through projection techniques, (ii) the use of monomer (subsystem) basis sets rather than supersystem basis sets, and (iii) the neglect of intersubsystem response coupling within the so-called uncoupled FDE (or FDEu) approximation. While approximation (i) is known to generally lead to large errors for covalently bound subsystems, we present cases in which either the basis set or the coupling step are similarly or even (much) more important. In particular, we explicitly demonstrate by comparison to a fully coupled calculation that missing intersubsystem response couplings are responsible for the failure of FDE reported in a previous study [ J. Chem. Theory Comput. 2015, 11, 5305-5315]. We show that good agreement with reference results can be obtained in this case even with standard NAKE approximations for the FDE potentials and efficient monomer basis sets, making calculations for larger systems well accessible.

The effect of N-heterocyclic carbene units on the absorption spectra of Fe(ii) complexes: a challenge for theory

The absorption spectra of five Fe(ii) homoleptic and heteroleptic complexes containing strong sigma-donating N-heterocyclic carbene (NHC) and polypyridyl ligands have been theoretically characterized using a tuned range-separation functional.

Piecewise linearity, freedom from self-interaction, and a Coulomb asymptotic potential: three related yet inequivalent properties of the exact density functional

Three properties of the exact energy functional of DFT are important in general and for spectroscopy in particular, but are not necessarily obeyed by approximate functionals. We explain what they are, why they are important, and how they are related yet inequivalent.

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.

Phosphahelicenes: From Chiroptical and Photophysical Properties to OLED Applications

Herein, the experimental physicochemical and chiroptical properties of a series of phosphahelicenes are reported, focusing on their UV/Vis absorption, luminescence, electronic circular dichroism, optical rotations, and circularly polarized luminescence. Furthermore, detailed analysis of absorption and ECD spectra performed with the help of quantum-chemical calculations allowed us to highlight general features of these helicenic phosphines. Finally, due to well-suited electrochemical properties and thermal stability, the systems were successfully used as emitters in organic light-emitting diodes.

Delocalization Errors in Density Functional Theory Are Essentially Quadratic in Fractional Occupation Number

Approximate functionals used in practical density functional theory (DFT) deviate from the piecewise linear behavior of the exact functional for fractional charges. This deviation causes excess charge delocalization, which leads to incorrect densities, molecular properties, barrier heights, band gaps, and excitation energies. We present a simple delocalization function for characterizing this error and find it to be almost perfectly linear vs the fractional electron number for systems spanning in size from the H atom to the C₁₂H₁₄ polyene. This causes the delocalization energy error to be a quadratic polynomial in the fractional electron number, which permits us to assess the comparative performance of 47 popular and recent functionals through the curvature. The quadratic form further suggests that information about a single fractional charge is sufficient to eliminate the principal source of delocalization error. Generalizing traditional two-point information like ionization potentials or electron affinities to account for a third, fractional charge-based data point could therefore permit fitting/tuning of functionals with lower delocalization error.

Synthesis and Combined Experimental and Theoretical Characterization of Dihydro-tetraaza-acenes

We present a combined experimental and theoretical study of electronic and optical properties of dihydro-tetraaza-acenes (DHTAn). Using solvent-free condensation, we are able to synthesize not only DHTA5 but also the longer DHTA6 and DHTA7 molecules. We then investigate their gas-phase electronic structures by means of ab initio density functional calculations employing an optimally tuned range-separated hybrid functional. By comparing with the parent linear oligoacenes (nA) and based on computed ionization potentials and electron affinities, we predict DHTAn molecules to be more stable than acenes of the same length, where we expect DHTAn molecules to be persistent at least up to n = 7 rings. We further exploit the analogy with nA by analyzing the entire intramolecular π-band structure of the DHTAn molecules. This clearly reveals that the additional two electrons donated by the dihydropyrazine group are delocalized over the entire molecule and contribute to its π-electron system. As a consequence, the symmetry of the frontier orbitals of DHTAn differs from that of the parent nA molecule. This also affects the UV-vis absorption spectra which have been measured for DHTA5, 6, and 7 dissolved in dimethyl sulfoxide and analyzed by means of excited state calculations within a time-dependent density functional theory framework.

TDDFT calculations and photoacoustic spectroscopy experiments used to identify phenolic acid functional biomolecules in Brazilian tropical fruits in natura

Time-dependent density functional theory (TDDFT) calculations of electronic transitions have been widely used to determine molecular structures. The excitation wavelengths and oscillator strengths obtained with the hybrid exchange-correlation functional B3LYP in conjunction with the ADZP basis set are employed to simulate the UV-Vis spectra of eight phenolic acids. Experimental and theoretical UV-Vis spectra reported previously in the literature are compared with our results. The fast, sensitive and non-destructive technique of photoacoustic spectroscopy (PAS) is used to determine the UV-Vis spectra of four Brazilian tropical fresh fruits in natura. Then, the PAS along with the TDDFT results are for the first time used to investigate and identify the presence of phenolic acids in the fruits studied in this work. This theoretical method with this experimental technique show to be a powerful and cheap tool to detect the existence of phenolic acids in fruits, vegetables, cereals, and grains. Comparison with high performance liquid chromatography results, when available, is also carried out.

Using non-empirically tuned range-separated functionals with simulated emission bands to model fluorescence lifetimes

Fluorescence lifetimes were evaluated using TD-DFT under different approximations for the emitting molecule and various exchange-correlation functionals, such as B3LYP, BMK, CAM-B3LYP, LC-BLYP, M06, M06-2X, M11, PBE0, ωB97, ωB97X, LC-BLYP*, and ωB97X* where the range-separation parameters in the last two functionals were tuned in a non-empirical fashion. Changes in the optimised molecular geometries between the ground and electronically excited states were found to affect the quality of the calculated lifetimes significantly, while the inclusion of vibronic features led to further improvements over the assumption of a vertical electronic transition. The LC-BLYP* functional was found to return the most accurate fluorescence lifetimes with unsigned errors that are mostly within 1.5 ns of experimental values.

Calculation of linear and nonlinear optical properties of azobenzene derivatives with Kohn–Sham and coupled-cluster methods

A non-empirically tuned generalized Kohn–Sham functional allows the prediction of accurate low-energy excitation energies and linear polarizabilities. Second hyperpolarizabilities are not improved when compared to coupled-cluster benchmark data.

Synthesis and Chiroptical Properties of Hexa-, Octa-, and Deca-azaborahelicenes: Influence of Helicene Size and of the Number of Boron Atoms

Four members of a new class of cycloborylated hexa-, octa-, and deca-helicenes (1 a-d) have been prepared in enantiopure form, along with two cycloplatinated deca-helicenes (1 d', 1 d¹ ), further extending the family of cycloplatinated hexa- and octa-helicenes reported previously. The azabora[n]helicenes display intense electronic circular dichroism and large optical rotations; the dependence of the optical activity on the size of the helix (n=6, 8, 10) and the number of boron atoms (1 or 2) has been examined in detail both experimentally and theoretically. The photophysical properties (nonpolarized and circularly polarized luminescence) of these new fluorescent organic helicenes have been measured and compared with the corresponding organometallic phosphorescent cycloplatinated derivatives (1 a¹ -d¹ ).

Optical Rotation Calculations for a Set of Pyrrole Compounds

Optical rotation of 14 molecules containing the pyrrole group is calculated by employing both time-dependent density functional theory (TDDFT) with the CAM-B3LYP functional and the second-order approximate coupled-cluster singles and doubles (CC2) method. All optical rotations have been provided using the aug-cc-pVDZ basis set at λ = 589 nm. The two methods predict similar results for both sign and magnitude for the optical rotation of all molecules. The obtained signs are consistent with experiments as well, although several conformers for four molecules needed to be studied to reproduce the experimental sign. We have also calculated excitation energies and rotatory strengths for the six lowest lying electronic transitions for several conformers of the two smallest molecules and found that each rotatory strength has various contributions for each conformer which can cause different optical rotations for different conformers of a molecule. Our results illustrate that both methods are able to reproduce the experimental optical rotations, and that the CAM-B3LYP functional, the least computationally expensive method used here, is an applicable and reliable method to predict the optical rotation for these molecules in line with previous studies.