Optimal Tuning of Range-Separated Hybrids for Solvated Molecules with Time-Dependent Density Functional Theory

Machine-Learning-Assisted Rational Design of Si─Rhodamine as Cathepsin-pH-Activated Probe for Accurate Fluorescence Navigation

High-performance fluorescent probes stand as indispensable tools in fluorescence-guided imaging, and are crucial for precise delineation of focal tissue while minimizing unnecessary removal of healthy tissue. Herein, machine-learning-assisted strategy to investigate the current available xanthene dyes is first proposed, and a quantitative prediction model to guide the rational synthesis of novel fluorescent molecules with the desired pH responsivity is constructed. Two novel Si─rhodamine derivatives are successfully achieved and the cathepsin/pH sequentially activated probe Si─rhodamine─cathepsin-pH (SiR─CTS-pH) is constructed. The results reveal that SiR─CTS-pH exhibits higher signal-to-noise ratio of fluorescence imaging, compared to single pH or cathepsin-activated probe. Moreover, SiR─CTS-pH shows strong differentiation abilities for tumor cells and tissues and accurately discriminates the complex hepatocellular carcinoma tissues from normal ones, indicating its significant application potential in clinical practice. Therefore, the continuous development of xanthene dyes and the rational design of superior fluorescent molecules through machine-learning-assisted model broaden the path and provide more advanced methods to researchers.

Unveiling the photophysical and excited state properties of multi-resonant OLED emitters using combined DFT and CCSD method

Multi-resonance thermally-activated delayed fluorescence (MR-TADF) is predominantly observed in organoboron heteroatom-embedded molecules, featuring enhanced performance in organic light-emitting diodes (OLEDs) with high color purity, chemical stability, and excellent photoluminescence quantum yields. However, predicting the impact of any chemical change remains a challenge. Computational methods including density functional theory (DFT) still require accurate descriptions of photophysical properties of MR-TADF emitters. To circumvent this drawback, we explored recent investigations on the CzBX (Cz = carbazole, X = O, S, or Se) molecule as a central building block. We constructed a series of MR-TADF molecules by controlling chalcogen atom embedding, employing a combined approach of DFT and coupled-cluster (CCSD) methods. Our predicted results for MR-TADF emitter molecules align with the reported experimental data in the literature. The variation in the positions of chalcogen atoms embedded within the CzBX2X framework imparts unique photophysical properties.

Machine-learning screening of luminogens with aggregation-induced emission characteristics for fluorescence imaging

Due to the excellent biocompatible physicochemical performance, luminogens with aggregation-induced emission (AIEgens) characteristics have played a significant role in biomedical fluorescence imaging recently. However, screening AIEgens for special applications takes a lot of time and efforts by using conventional chemical synthesis route. Fortunately, artificial intelligence techniques that could predict the properties of AIEgen molecules would be helpful and valuable for novel AIEgens design and synthesis. In this work, we applied machine learning (ML) techniques to screen AIEgens with expected excitation and emission wavelength for biomedical deep fluorescence imaging. First, a database of various AIEgens collected from the literature was established. Then, by extracting key features using molecular descriptors and training various state-of-the-art ML models, a multi-modal molecular descriptors strategy has been proposed to extract the structure-property relationships of AIEgens and predict molecular absorption and emission wavelength peaks. Compared to the first principles calculations, the proposed strategy provided greater accuracy at a lower computational cost. Finally, three newly predicted AIEgens with desired absorption and emission wavelength peaks were synthesized successfully and applied for cellular fluorescence imaging and deep penetration imaging. All the results were consistent successfully with our expectations, which demonstrated the above ML has a great potential for screening AIEgens with suitable wavelengths, which could boost the design and development of novel organic fluorescent materials.

Spatial Decay and Limits of Quantum Solute-Solvent Interactions

Molecular excitations in the liquid-phase environment are renormalized by the surrounding solvent molecules. Herein, we employ the GW approximation to investigate the solvation effects on the ionization energy of phenol in various solvent environments. The electronic effects differ by up to 0.4 eV among the five investigated solvents. This difference depends on both the macroscopic solvent polarizability and the spatial decay of the solvation effects. The latter is probed by separating the electronic subspace and the GW correlation self-energy into fragments. The fragment correlation energy decays with increasing intermolecular distance and vanishes at ∼9 Å, and this pattern is independent of the type of solvent environment. The 9 Å cutoff defines an effective interacting volume within which the ionization energy shift per solvent molecule is proportional to the macroscopic solvent polarizability. Finally, we propose a simple model for computing the ionization energies of molecules in an arbitrary solvent environment.

Probing aqueous ions with non-local Auger relaxation

Non-local analogues of Auger decay are increasingly recognized as important relaxation processes in the condensed phase. Here, we explore non-local autoionization, specifically Intermolecular Coulombic Decay (ICD), of a series of aqueous-phase isoelectronic cations following 1s core-level ionization. In particular, we focus on Na+, Mg2+, and Al3+ ions. We unambiguously identify the ICD contribution to the K-edge Auger spectrum. The different strength of the ion-water interactions is manifested by varying intensities of the respective signals: the ICD signal intensity is greatest for the Al3+ case, weaker for Mg2+, and absent for weakly-solvent-bound Na+. With the assistance of ab initio calculations and molecular dynamics simulations, we provide a microscopic understanding of the non-local decay processes. We assign the ICD signals to decay processes ending in two-hole states, delocalized between the central ion and neighbouring water. Importantly, these processes are shown to be highly selective with respect to the promoted water solvent ionization channels. Furthermore, using a core-hole-clock analysis, the associated ICD timescales are estimated to be around 76 fs for Mg2+ and 34 fs for Al3+. Building on these results, we argue that Auger and ICD spectroscopy represents a unique tool for the exploration of intra- and inter-molecular structure in the liquid phase, simultaneously providing both structural and electronic information.

Poisoning density functional theory with benchmark sets of difficult systems

Large benchmark sets like GMTKN55 [Goerigk et al., Phys. Chem. Chem. Phys., 2017, 19, 32184] let us analyse the performance of density functional theory over a diverse range of systems and bonding types. However, assessing over a large and diverse set can miss cases where approaches fail badly, and can give a misleading sense of security. To this end we introduce a series of 'poison' benchmark sets, P30-5, P30-10 and P30-20, comprising systems with up to 5, 10 and 20 atoms, respectively. These sets represent the most difficult-to-model systems in GMTKN55. We expect them to be useful in developing new approximations, identifying weak points in existing ones, and to aid in selecting appropriate DFAs for computational studies involving difficult physics, e.g. catalysis.

Linear fractional charge behavior in density functional theory through dielectric tuning of conductor-like polarizable continuum model

A property of exact density functional theory is linear fractional charge behavior as electrons are added or removed from a molecule. Typical density functional approximations (DFAs) exhibit delocalization error, which overstabilizes this fractional charge. Conversely, solvent corrections have been shown to erroneously destabilize this fractional charge. This work will show that an implicit solvent correction with a tuned dielectric can be used as an ad hoc correction to offset the delocalizing character of DFAs and achieve linear fractional charge behavior. While desirable, in principle, we find that this linear charge behavior degrades the vertical ionization energies reported by DFAs. Our results reveal that the localizing character of the solvent correction and the Hartree-Fock (HF) exchange offset each other. This helps explain the decreased ratios of HF exchange to DFA exchange in long-range hybrid tuning studies that use a solvent correction.

The power of trichlorosilylation: isolable trisilylated allyl anions, allyl radicals, and allenyl anions

Treatment of hexachloropropene (Cl₂C Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> C(Cl)–CCl₃) with Si₂Cl₆ and [nBu₄N]Cl (1 : 4 : 1) in CH₂Cl₂ results in a quantitative conversion to the trisilylated, dichlorinated allyl anion salt [nBu₄N][Cl₂CC(SiCl₃)–C(SiCl₃)₂] ([nBu₄N][1]). Tetrachloroallene Cl₂CCCCl₂ was identified as the first intermediate of the reaction cascade. In the solid state, [1]⁻ adopts approximate C_s symmetry with a dihedral angle between the planes running through the olefinic and carbanionic fragments of [1]⁻ of CC–Si//Si–C–Si = 78.3(1)°. One-electron oxidation of [nBu₄N][1] with SbCl₅ furnishes the distillable blue radical 1˙. The neutral propene Cl₂CC(SiCl₃)–C(SiCl₃)₂H (2) was obtained by (i) protonation of [1]⁻ with HOSO₂CF₃ (HOTf) or (ii) H-atom transfer to 1˙ from 1,4-cyclohexadiene. Quantitative transformation of all three SiCl₃ substituents in 2 to Si(OMe)₃ (2^OMe) or SiMe₃ (2^Me) substituents was achieved by using MeOH/NMe₂Et or MeMgBr in CH₂Cl₂ or THF, respectively. Upon addition of 2 equiv. of tBuLi, 2^Me underwent deprotonation with subsequent LiCl elimination, 1,2-SiMe₃ migration and Cl/Li exchange to afford the allenyl lithium compound Me₃Si(Li)CCC(SiMe₃)₂ (Li[4]), which is an efficient building block for the introduction of Me, SiMe₃, or SnMe₃ (5) groups. The trisilylated, monochlorinated allene Cl₃Si(Cl)CCC(SiCl₃)₂ (6), was obtained from [nBu₄N][1] through Cl⁻-ion abstraction with AlCl₃ and rearrangement in CH₂Cl₂ (1˙ forms as a minor side product, likely because the system AlCl₃/CH₂Cl₂ can also act as a one-electron oxidant).

Treatment of hexachloropropene (Cl₂CC(Cl)–CCl₃) with Si₂Cl₆ and [nBu₄N]Cl (1 : 4 : 1) in CH₂Cl₂ results in a quantitative conversion to the trisilylated, dichlorinated allyl anion salt [nBu₄N][Cl₂CC(SiCl₃)–C(SiCl₃)₂] ([nBu₄N][1]).

Efficient treatment of molecular excitations in the liquid phase environment via stochastic many-body theory

Accurate predictions of charge excitation energies of molecules in the disordered condensed phase are central to the chemical reactivity, stability, and optoelectronic properties of molecules and critically depend on the specific environment. Herein, we develop a stochastic GW method for calculating these charge excitation energies. The approach employs maximally localized electronic states to define the electronic subspace of a molecule and the rest of the system, both of which are randomly sampled. We test the method on three solute-solvent systems: phenol, thymine, and phenylalanine in water. The results are in excellent agreement with the previous high-level calculations and available experimental data. The stochastic calculations for supercells containing up to 1000 electrons representing the solvated systems are inexpensive and require ≤1000 central processing unit hrs. We find that the coupling with the environment accounts for ∼40% of the total correlation energy. The solvent-to-solute feedback mechanism incorporated in the molecular correlation term causes up to 0.6 eV destabilization of the quasiparticle energy. Simulated photo-emission spectra exhibit red shifts, state-degeneracy lifting, and lifetime shortening. Our method provides an efficient approach for an accurate study of excitations of large molecules in realistic condensed phase environments.

Machine Learning Enables Highly Accurate Predictions of Photophysical Properties of Organic Fluorescent Materials: Emission Wavelengths and Quantum Yields

The development of functional organic fluorescent materials calls for fast and accurate predictions of photophysical parameters for processes such as high-throughput virtual screening, while the task is challenged by the limitations of quantum mechanical calculations. We establish a database covering >4300 solvated organic fluorescent dyes with 3000 distinct compounds and develop a new machine learning approach aimed at efficient and accurate predictions of emission wavelength and photoluminescence quantum yield (PLQY). Our feature engineering has given rise to a functionalized structure descriptor (FSD) and a comprehensive general solvent descriptor (CGSD), whereby a highly black-box computational framework is realized with consistently good accuracy across different dye families, ability of describing substitution effects and solvent effects, efficiency for large-scale predictions, and workability with on-the-fly learning. Evaluations with unseen molecules suggest a remarkable mean absolute error of 0.13 for PLQY and 0.080 eV for emission energy, the latter comparable to time-dependent density functional theory (TD-DFT) calculations. An online prediction platform was constructed based on the ensemble model to make predictions in various solvents. Our statistical learning methodology will complement quantum mechanical calculations as an efficient alternative approach for the prediction of these parameters.

Performance of TDDFT Vertical Excitation Energies of Core-Substituted Naphthalene Diimides

We have evaluated the performance of various density functionals, covering generalized gradient approximation (GGA), global hybrid (GH) and range‐separated hybrid (RSH), using time dependent density functional theory (TDDFT) for computing vertical excitation energies against experimental absorption maximum (λ_max) for a set of 10 different core‐substituted naphthalene diimides (cNDI) recorded in dichloromethane. The computed excitation in case of GH PBE0 is most accurate while the trend is most systematic with RSH LCY‐BLYP compared to λ_max. We highlight the importance of including solvent effects for optimal agreement with the λ_max. Increasing the basis set size from TZ2P to QZ4P has a negligible influence on the computed excitation energies. Notably, RSH CAMY‐B3LYP gave the least error for charge‐transfer excitation. The poorest agreement with λ_max is obtained with semi‐local GGA functionals. Use of the optimally‐tuned RSH LCY‐BLYP* is not recommended because of the high computational cost and marginal improvement in results.

Synthesis of Bis-β-Diketonate Lanthanide Complexes with an Azobenzene Bridge and Studies of Their Reversible Photo/Thermal Isomerization Properties

The ligand, bis-β-diketone with an azobenzene bridge (4,4'-(4,4,4-trifluoro-1,3-butanedione)azobenzene, H 2 L), was prepared for the synthesis of a series of dinuclear lanthanide complexes with the formula [Ln 2 L 3 (DMSO) 4 ] (Ln = Eu3+, Gd3+, Tb3+, and DMSO = dimethyl sulfoxide). X-ray crystallographic analysis reveals that the three complexes are triple-stranded dinuclear structures formed by three bis-β-diketonate ligands with two lanthanide ions (Ln3+). The trans-to-cis photoisomerization rates of the azobenzene group of the three [Ln 2 L 3 (DMSO) 4 ] complexes in ethanol and acetonitrile solutions are similar to those of the pure H 2 L ligand and other azobenzene-containing mononuclear lanthanide complexes, but the trans-to-cis quantum yields (Φt→c = 10-3) are 1 order of magnitude smaller. The first-order rate constant for the cis-to-trans thermal isomerization at 50 °C of the H 2 L ligand is similar to those of azobenzene derivatives, while those for the [Ln 2 L 3 (DMSO) 4 ] complexes (k iso = 10-4 s-1) are higher than those of the mononuclear azobenzene-containing lanthanide complexes. Furthermore, as the lanthanide ionic radius becomes smaller from Eu3+ to Gd3+ to Tb3+, the thermal isomerization rate constant decreases and the half-life increases. All these results are proposed to arise from the rigidity at both ends of the azo group by coordination to the dinuclear lanthanide ions and the different isomerization mechanisms. These are the first examples of bis-β-diketonate dinuclear lanthanide complexes with an azobenzene bridge and help illustrate the mechanism of azobenzene isomerization.

Beyond Koopmans' theorem: electron binding energies in disordered materials

The topical review focuses on calculating ionization energies (IE), or electronic polarons in quasi-particle terminology, in large disordered systems, e.g. for a solute dissolved in a molecular solvent. The simplest estimate of the ionization energy is provided by one-electron energies in the Hartree-Fock theory, but the calculated quantities are not accurate. Density functional theory as many-body theory provides a principal opportunity for calculating one-electron energies including correlation and relaxation effects, i.e. the true energies of electronic polarons. We argue that such a principal possibility materializes within the concept of optimally tuned range-separated hybrid functionals (OT-RSH). We describe various schemes for optimal tuning. Importantly, the OT-RSH scheme is investigated for systems capped with dielectric continuum, providing a consistent picture on the QM/dielectric boundary. Finally, some limitations and open issues of the OT-RSH approach are addressed.

Do water's electrons care about electrolytes?

Ions have a profound effect on the geometrical structure of liquid water and an aqueous environment is known to change the electronic structure of ions. Here we combine photoelectron spectroscopy measurements from liquid microjets with molecular dynamical and quantum chemical calculations to address the reverse question, to what extent do ions affect the electronic structure of liquid water? We study aqueous solutions of sodium iodide (NaI) over a wide concentration range, from nearly pure water to 8 M solutions, recording spectra in the 5 to 60 eV binding energy range to include all water valence and the solute Na+ 2p, I- 4d, and I- 5p orbital ionization peaks. We observe that the electron binding energies of the solute ions change only slightly as a function of electrolyte concentration, less than 150 ± 60 meV over an ∼8 M range. Furthermore, the photoelectron spectrum of liquid water is surprisingly mildly affected as we transform the sample from a dilute aqueous salt solution to a viscous, crystalline-like phase. The most noticeable spectral changes are a negative binding energy shift of the water 1b2 ionizing transition (up to -370 ± 60 meV) and a narrowing of the flat-top shape water 3a1 ionization feature (up to 450 ± 90 meV). A novel computationally efficient technique is introduced to calculate liquid-state photoemission spectra using small clusters from molecular dynamics (MD) simulations embedded in dielectric continuum. This theoretical treatment captured the characteristic positions and structures of the aqueous photoemission peaks, reproducing the experimentally observed narrowing of the water 3a1 feature and weak sensitivity of the water binding energies to electrolyte concentration. The calculations allowed us to attribute the small binding energy shifts to ion-induced disruptions of intermolecular electronic interactions. Furthermore, they demonstrate the importance of considering concentration-dependent screening lengths for a correct description of the electronic structure of solvated systems. Accounting for electronic screening, the calculations highlight the minimal effect of electrolyte concentration on the 1b1 binding energy reference, in accord with the experiments. This leads us to a key finding that the isolated, lowest-binding-energy, 1b1, photoemission feature of liquid water is a robust energetic reference for aqueous liquid microjet photoemission studies.

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.