Intramolecular Charge Transfer with 1-tert-Butyl-6-cyano-1,2,3,4-tetrahydroquinoline (NTC6) and Other Aminobenzonitriles. A Comparison of Experimental Vapor Phase Spectra and Crystal Structures with Calculations

A mountaineering strategy to excited states: Accurate vertical transition energies and benchmarks for substituted benzenes

In an effort to expand the existing QUEST database of accurate vertical transition energies [Véril et al. WIREs Comput. Mol. Sci. 2021, 11, e1517], we have modeled more than 100 electronic excited states of different natures (local, charge-transfer, Rydberg, singlet, and triplet) in a dozen of mono- and di-substituted benzenes, including aniline, benzonitrile, chlorobenzene, fluorobenzene, nitrobenzene, among others. To establish theoretical best estimates for these vertical excitation energies, we have employed advanced coupled-cluster methods including iterative triples (CC3 and CCSDT) and, when technically possible, iterative quadruples (CC4). These high-level computational approaches provide a robust foundation for benchmarking a series of popular wave function methods. The evaluated methods all include contributions from double excitations (ADC(2), CC2, CCSD, CIS(D), EOM-MP2, STEOM-CCSD), along with schemes that also incorporate perturbative or iterative triples (ADC(3), CCSDR(3), CCSD(T)(a)  ⋆ , and CCSDT-3). This systematic exploration not only broadens the scope of the QUEST database but also facilitates a rigorous assessment of different theoretical approaches in the framework of a homologous chemical series, offering valuable insights into the accuracy and reliability of these methods in such cases. We found that both ADC(2.5) and CCSDT-3 can provide very consistent estimates, whereas among less expensive methods SCS-CC2 is likely the most effective approach. Importantly, we show that some lower order methods may offer reasonable trends in the homologous series while providing quite large average errors, and vice versa. Consequently, benchmarking the accuracy of a model based solely on absolute transition energies may not be meaningful for applications involving a series of similar compounds.

Linearized Pair-Density Functional Theory

Multiconfiguration pair-density functional theory (MC-PDFT) is a post-SCF multireference method that has been successful at computing ground- and excited-state energies. However, MC-PDFT is a single-state method in which the final MC-PDFT energies do not come from diagonalization of a model-space Hamiltonian matrix, and this can lead to inaccurate topologies of potential energy surfaces near locally avoided crossings and conical intersections. Therefore, in order to perform physically correct ab initio molecular dynamics with electronically excited states or to treat Jahn-Teller instabilities, it is necessary to develop a PDFT method that recovers the correct topology throughout the entire nuclear configuration space. Here we construct an effective Hamiltonian operator, called the linearized PDFT (L-PDFT) Hamiltonian, by expanding the MC-PDFT energy expression to first order in a Taylor series of the wave function density. Diagonalization of the L-PDFT Hamiltonian gives the correct potential energy surface topology near conical intersections and locally avoided crossings for a variety of challenging cases including phenol, methylamine, and the spiro cation. Furthermore, L-PDFT outperforms MC-PDFT and previous multistate PDFT methods for predicting vertical excitations from a variety of representative organic chromophores.

Assessing the Role of the Kohn-Sham Density in the Calculation of the Low-Lying Bethe-Salpeter Excitation Energies

We adopt the GW many-body perturbation theory in conjunction with the Bethe-Salpeter equation (BSE) to compute 57 excitation energies of a set of 37 molecules. By using the PBEh global hybrid functional and a self-consistent scheme on the eigenvalues in GW, we show a strong dependence of the BSE energy on the starting Kohn-Sham (KS) density functional. This arises from both the quasiparticle energies and the spatial localization of the frozen KS orbitals employed to compute the BSE. In order to address the arbitrariness in the mean field choice, we adopt an orbital-tuning scheme where the amount of Fock exchange, α, is tuned to impose that the KS HOMO matches the GW quasiparticle eigenvalue, thus fulfilling the ionization potential theorem in DFT. The performance of the proposed scheme yields excellent results and it is similar to M06-2X and PBEh with α = 75%, consistent with tuned values of α ranging between 60% and 80%.

Ultrafast Carrier Dynamics of Non-fullerene Acceptors with Different Planarity: Impact of Steric Hindrance

Most high-performance non-fullerene acceptors are of the acceptor-donor-acceptor (A-D-A)-type structure. Under photoexcitation, the intramolecular charge transfer effect on the A-D-A framework results in a large dipole moment change, facilitating the efficient generation of charge carriers. Achieving more efficient intramolecular charge transfer by adjusting the molecular structure is one of the current research ideas. Recently, we found that the power conversion efficiency can be improved from 4.41 to 13.13% by tuning the planarity of the non-fused ring electron acceptor backbone through steric hindrance of lateral substituents. We found that the planar backbone can effectively improve the intramolecular charge transfer, which has a great influence on the power conversion efficiency of the device. Our results demonstrate that charge transfer dynamics can be controlled by optimizing steric hindrance, which plays a crucial role in the photovoltaic performance of organic solar cells.

Simulating the Nonadiabatic Relaxation Dynamics of 4-(N,N-Dimethylamino)benzonitrile (DMABN) in Polar Solution

The compound 4-(N,N-dimethylamino)benzonitrile (DMABN) represents the archetypal system for dual fluorescence, a rare photophysical phenomenon in which a given fluorophore shows two distinct emission bands. Despite extensive studies, the underlying mechanism remains the subject of debate. In the present contribution, we address this issue by simulating the excited-state relaxation process of DMABN as it occurs in polar solution. The potential energy surfaces for the system are constructed with the use of the additive quantum mechanics/molecular mechanics (QM/MM) method, and the coupled dynamics of the electronic wave function and the nuclei is propagated with the semiclassical fewest switches surface hopping method. The DMABN molecule, which comprises the QM subsystem, is treated with the use of the second-order algebraic diagrammatic construction (ADC(2)) method with the imposition of spin-opposite scaling (SOS). It is verified that this level of theory achieves a realistic description of the excited-state potential energy surfaces of DMABN. The simulation results qualitatively reproduce the main features of the experimentally observed fluorescence spectrum, thus allowing the unambiguous assignment of the two fluorescence bands: the normal band is due to the near-planar locally excited (LE) structure of DMABN, while the so-called "anomalous" second band arises from the twisted intramolecular charge transfer (TICT) structure. The transformation of the LE structure into the TICT structure takes place directly via intramolecular rotation, and is not mediated by another excited-state structure. In particular, the oft-discussed rehybridized intramolecular charge transfer (RICT) structure, which is characterized by a bent nitrile group, does not play a role in the relaxation process.

Discrete and continuum modeling of solvent effects in a twisted intramolecular charge transfer system: The 4-N,N-dimethylaminobenzonitrile (DMABN) molecule

The 4-N,N-dimethylaminobenzonitrile (DMABN) molecule is a prototypical system displaying twisted intramolecular (TICT) charge transfer effects. The ground and the first four electronic excited states (S₁-S₄) in gas phase and upon solvation were studied. Charge transfer values as function of the torsion angle between the donor group (dimethylamine) and the acceptor moiety (benzonitrile) were explicitly computed. Potential energy curves were also obtained. The algebraic diagrammatic construction method at the second-order [ADC(2)] ab initio wave function was employed. Three solvents of increased polarities (benzene, DMSO and water) were investigated using discrete (average solvent electrostatic configuration - ASEC) and continuum (conductor-like screening model - COSMO) models. The results for the S₃ and S₄ excited states and the S₁-S₄ charge transfer curves were not previously available in the literature. Electronic gas phase and solvent vertical spectra are in good agreement with previous theoretical and experimental results. In the twisted (90°) geometry the optical oscillator strengths have negligible values even for the S₂ bright state. Potential energy curves show two distinct pairs of curves intersecting at decreasing angles or not crossing in the more polar solvents. Charge transfer and electric dipole values allowed the rationalization of these results. The former effects are mostly independent of the solvent model and polarity. Although COSMO and ASEC solvent models mostly lead to similar results, there is an important difference: some crossings of the excitation energy curves appear only in the ASEC solvation model, which has important implications to the photochemistry of DMABN.

D-A type sensor array for differentiation and identification of white wine varieties based on specific solvent effect activated by CT-LE transition

In this work, we synthesize a series of compounds with electron donor (D) and acceptor (A) units. They show general solvent effect in aprotic solvents, suggesting a charge transfer (CT) process. While in protic solvents including water, ethanol and methanol, the spectra exert no polarity-dependence but a remarkable hypochromatic shift together with the fading of CT band. Dynamic analysis implies that intermolecular hydrogen bond will be formed between carboxylic acid and protic solvent, boosting another deactivation pathway that jumps off a bigger energy gap, in other words, favoring the locally excited (LE) state emission. The CT-LE transition involves variations in both absorption and emission spectra, and further poses competition with other mechanisms including activated/restricted intramolecular rotation (IR/RIR). Inspired by the cross-reactivity, we turn our attention to the development of sensor array, in order to identify white wine varieties. The differential spectral responses are recorded, generating multiple factors including absorption wavelength (λ_ab), emission wavelength (λ_em), absorbance (Abs.) and emission intensity (Int.). These factors are processed with principal component analysis (PCA), creating a three-dimensional fingerprint data base for white wines. The data points in the coordinate system are clustered into 10 different groups, demonstrating a clear differentiation of all the white wines. More importantly, as our final test for whether the sensor array can identify the counterfeits, an adulterated liquor sample, which is provided by police officers, is fingerprinted on the three-dimensional diagram. Its canonical factors fall into an area distinct from the adulterated wine, indicating a clear identification.

Absence of Intramolecular Charge Transfer with 4-Fluoro-N,N-dimethylaniline (DMA4F), Contrary to an Experimental Report Supported by Computations

With 4-fluoro-N,N-dimethylaniline (DMA4F), only a single fluorescence from a locally excited (LE) state is observed, irrespective of solvent polarity, temperature, and excitation wavelength. The relatively small excited state dipole moment μ_e = 7.3 D confirms the identification as LE. The single exponential fluorescence decays in the nonpolar n-hexane (2.04 ns) and in the strongly polar acetonitrile (5.73 ns) are a further support. Similar results are obtained with 4-chloro-N,N-dimethylaniline (DMA4Cl), having a chlorobenzene subgroup, a somewhat better electron acceptor than the fluorobenzene moiety in DMA4F. The absence of intramolecular charge transfer (ICT) with DMA4F is in accord with its large energy gap ΔE(S₁,S₂) of 8300 cm^-1 in n-hexane between the two lowest singlet excited states, which is even larger than that (6300 cm^-1) of N,N-dimethylaniline (DMA), for which an LE → ICT reaction likewise does not occur. The results with DMA4F are in contradiction with a publication by Fujiwara et al. ( Chem. Phys. Lett. 2013 , 586 , 70 ), in which the appearance of dual LE + ICT emission is reported for DMA4F in n-hexane and MeCN at room temperature. The ICT/LE fluorescence quantum yield ratio Φ'(ICT)/Φ(LE) reached a maximum value of ∼2, in n-hexane and surprisingly also in MeCN, as the excitation wavelength approaches the red-edge of the absorption spectrum. These, in our opinion, erroneous observations were supported by time-dependent density functional theory (TDDFT) calculations, which compute a perpendicularly twisted lowest ICT state (TICT) state. This is a further example of the general tendency of computations to find a TICT conformation for the lowest excited singlet state of electron donor/acceptor molecules such as p-substituted anilines.

Multiconfiguration Pair-Density Functional Theory Is as Accurate as CASPT2 for Electronic Excitation

A correct description of electronically excited states is critical to the interpretation of visible-ultraviolet spectra, photochemical reactions, and excited-state charge-transfer processes in chemical systems. We have recently proposed a theory called multiconfiguration pair-density functional theory (MC-PDFT), which is based on a combination of multiconfiguration wave function theory and a new kind of density functional called an on-top density functional. Here, we show that MC-PDFT with a first-generation on-top density functional performs as well as CASPT2 for an organic chemistry database including valence, Rydberg, and charge-transfer excitations. The results are very encouraging for practical applications.

Tailoring of spectral response and intramolecular charge transfer in β-enaminones through band gap tuning: synthesis, spectroscopy and quantum chemical studies

In this paper, we investigate the synthetic tailoring of the spectral response and intramolecular charge transfer (ICT) of β-enaminones through bandgap modulation.

Two-State Intramolecular Charge Transfer (ICT) with 3,5-Dimethyl-4-(dimethylamino)benzonitrile (MMD) and Its Meta-Isomer mMMD. Ground State Amino Twist Not Essential for ICT

From X-ray structure analysis, amino twist angles of 90.0° for 2,4-dimethyl-3-(dimethylamino)benzonitrile (mMMD), 82.7° for 4-(di-tert-butylamino)benzonitrile (DTABN), and 88.7° for 6-cyanobenzoquinuclidine (CBQ) are determined, all considerably larger than the 57.4° of 3,5-dimethyl-4-(dimethylamino)benzonitrile (MMD). This large twist leads to lengthening of the amino-phenyl bond, 143.5 pm (mMMD), 144.1 pm (DTABN), 144.6 pm (CBQ), and 141.4 pm (MMD), as compared with 136.5 pm for the planar 4-(dimethylamino)benzonitrile (DMABN). As a consequence, the electronic coupling between the amino and phenyl subgroups in mMMD, DTABN, CBQ, and MMD is much weaker than in DMABN, as seen from the strongly reduced molar absorption coefficients. The fluorescence spectrum of MMD in n-hexane at 25 °C consists of two emissions, from a locally excited (LE) and an intramolecular charge transfer (ICT) state, with a fluorescence quantum yield ratio Φ'(ICT)/Φ(LE) of 12.8. In MeCN, a single ICT emission is found. With mMMD in n-hexane, in contrast, only LE fluorescence is observed, whereas the spectrum in MeCN originates from the ICT state. These differences are also seen from the half-widths of the overall fluorescence bands, which in n-hexane are larger for MMD than for mMMD, decreasing with solvent polarity for MMD and increasing for mMMD, reflecting the disappearance of LE and the onset of ICT in the overall spectra, respectively. From solvatochromic measurements the dipole moments μe(ICT) of MMD (16 D) and mMMD (15 D) are obtained. Femtosecond excited state absorption (ESA) spectra at 22 °C, together with the dual (LE + ICT) fluorescence, reveal that MMD in n-hexane undergoes a reversible LE ⇄ ICT reaction, with LE as the precursor, with a forward rate constant ka = 5.6 × 10(12) s(-1) and a back-reaction kd ∼ 0.05 × 10(12) s(-1). With MMD in the strongly polar solvent MeCN, ICT is faster: ka = 10 × 10(12) s(-1). In the case of mMMD in n-hexane, the ESA spectra show that ICT does not take place, contrary to MeCN, in which ka = 2.5 × 10(12) s(-1). The ICT reactions with MMD and mMMD are much faster than that of the parent compound DMABN in MeCN, with ka = 0.24 × 10(12) s(-1). Because of the very short ICT reaction times of 180 fs (MMD, n-hexane), 100 fs (MMD, MeCN), and 400 fs (mMMD, MeCN), it is clear that the picosecond fluorescence decays of these systems appear to be single exponential, due to the insufficient time resolution of 3 ps. It is concluded that the faster LE → ICT reaction of MMD as compared with DMABN (ka = 0.24 × 10(12) s(-1) in MeCN) is caused by a smaller energy gap ΔE(S1,S2) between the lowest singlet excited states and not by the large amino twist angle. Similarly, the larger ΔE(S1,S2) of mMMD as compared with MMD is held responsible for its smaller ICT efficiency (no reaction in n-hexane).

Intramolecular Charge-Transfer Excited-State Processes in 4-(N,N-Dimethylamino)benzonitrile: The Role of Twisting and the πσ* State

The structural processes leading to dual fluorescence of 4-(dimethylamino)benzonitrile in the gas phase and in acetonitrile solvent were investigated using a combination of multireference configuration interaction (MRCI) and the second-order algebraic diagrammatic construction (ADC(2)) methods. Solvent effects were included on the basis of the conductor-like screening model. The MRCI method was used for computing the nonadiabatic interaction between the two lowest excited ππ* states (S₂(L_a, CT) and S₁(L_b, LE)) and the corresponding minimum on the crossing seam (MXS) whereas the ADC(2) calculations were dedicated to assessing the role of the πσ* state. The MXS structure was found to have a twisting angle of ∼50°. The branching space does not contain the twisting motion of the dimethylamino group and thus is not directly involved in the deactivation process from S₂ to S₁. Polar solvent effects are not found to have a significant influence on this situation. Applying C_s symmetry restrictions, the ADC(2) calculations show that CCN bending leads to a strong stabilization and to significant charge transfer (CT). Nevertheless, this structure is not a minimum but converts to the local excitation (LE) structure on releasing the symmetry constraint. These findings suggest that the main role in the dynamics is played by the nonadiabatic interaction of the LE and CT states and that the main source for the dual fluorescence is the twisted internal charge-transfer state in addition to the LE state.

The driving forces for twisted or planar intramolecular charge transfer

The driving forces for twisting or planarization of excited D–A (donor–acceptor)-type chromophore have been investigated.

Triple emission from p-dimethylaminobenzonitrile-cucurbit[8]uril triggers the elusive excimer emission

The intriguing dual-emission behavior of p- dimethylaminobenzonitrile (DMABN) and the identity of the associated excited states is, arguably, the most extensively investigated and also controversially discussed molecule- specific phenomenon of modern photochemistry. We have now found a new, third fluorescence band when DMABN is encapsulated within the water-soluble molecular container cucurbit[8]uril (CB8). It is centered between the previously observed emissions and assigned to the elusive excimer emission from DMABN through 1:2 CB8:DMABN complex formation. Heating of the CB8⋅(DMABN)2 complex from 0 to 100 °C results in the dissociation of the ternary complex and restoration of the dual-emission properties of the monomer. Alternatively, monomer emission can be obtained by selecting cucurbit[7]uril (CB7), a host homologue that is too small to accommodate two DMABN molecules, or by introducing ethyl instead of methyl groups at the amino terminus of the aminobenzonitrile guest.

Red-shifted fluorescent aminated derivatives of a conformationally locked GFP chromophore

A novel class of fluorescent dyes based on conformationally locked GFP chromophore is reported. These dyes are characterized by red-shifted spectra, high fluorescence quantum yields and pH-independence in physiological pH range. The intra- and intermolecular mechanisms of radiationless deactivation of ABDI-BF2 fluorophore by selective structural locking of various conformational degrees of freedom were studied. A unique combination of solvatochromic and lipophilic properties together with "infinite" photostability (due to a dynamic exchange between free and bound dye) makes some of the novel dyes promising bioinspired tools for labeling cellular membranes, lipid drops and other organelles.

Intramolecular charge transfer with crystal violet lactone in acetonitrile as a function of temperature: reaction is not solvent-controlled

Intramolecular charge transfer (ICT) with crystal violet lactone (CVL) in the excited singlet state takes place in solvents more polar than n-hexane, such as ethyl acetate, tetrahydrofuran, and acetonitrile (MeCN). In these solvents, the fluorescence spectrum of CVL consists of two emission bands, from a locally excited (LE) and an ICT state. The dominant deactivation channel of the lowest excited singlet state is internal conversion, as the quantum yields of fluorescence (0.007) and intersystem crossing (0.015) in MeCN at 25 °C are very small. CVL is a weakly coupled electron donor/acceptor (D/A) molecule, similar to an exciplex (1)(A(-)D(+)). A solvatochromic treatment of the LE and ICT emission maxima results in the dipole moments μe(LE) = 17 D and μe(ICT) = 33 D, much larger than those previously reported. This discrepancy is attributed to different Onsager radii and spectral fluorimeter calibration. The LE and ICT fluorescence decays of CVL in MeCN are double exponential. As determined by global analysis, the LE and ICT decays at 25 °C have the times τ2 = 9.2 ps and τ1 = 1180 ps, with an amplitude ratio of 35.3 for LE. From these parameters, the rate constants ka = 106 × 10(9) s(-1) and kd = 3.0 × 10(9) s(-1) of the forward and backward reaction in the LE ⇄ ICT equilibrium are calculated, resulting in a free enthalpy difference ΔG of -8.9 kJ/mol. The amplitude ratio of the ICT fluorescence decay equals -1.0, which signifies that the ICT state is not prepared by light absorption in the S0 ground state, but originates exclusively from the directly excited LE precursor. From the temperature dependence of the fluorescence decays of CVL in MeCN (-45 to 75 °C), activation energies E(a) = 3.9 kJ/mol (LE → ICT) and E(d) = 23.6 kJ/mol (ICT → LE) are obtained, giving an enthalpy difference ΔH (= E(a) - E(d)) of -19.7 kJ/mol, and an entropy difference ΔS = -35.5 J mol(-1) K(-1). These data show that the ICT reaction of CVL in MeCN is not barrierless. The ICT reaction time of 9.2 ps is much longer than the mean solvent relaxation time of MeCN (0.26 ps), indicating, in contrast with earlier reports in the literature, that the reaction is not solvent controlled. This conclusion is supported by the observation of double exponential LE and ICT fluorescence with the same decay times.

Dual Fluorescence of Fluorazene in Solution: A Computational Study

The fluorazene molecule presents dual fluorescence in polar solvents. Its absorption and emission properties in gas phase and in acetonitrile solution have been studied theoretically using the complete active space second-order perturbation//complete active space self-consistent field quantum methodology and average solvent electrostatic potential from molecular dynamics for the solvent effects. In gas phase, two optimized excited-state geometries were obtained, one of them corresponds to a local excitation (LE), and the other is an intramolecular charge transfer (ICT) and lies higher in energy. In acetonitrile solution, a second ICT structure where the molecule remains planar is found, and the energy differences are reduced. Fluorescence energies from LE and the planar ICT have a good agreement with the experimental bands, but emission from the bent ICT has too low an energy.