Spectral Tuning and Excitation-Energy Transfer by Unique Carotenoids in Diatom Light-Harvesting Antenna.

The light-harvesting antennae of diatoms and spinach are composed of similar chromophores; however, they exhibit different absorption wavelengths. Recent advances in cryoelectron microscopy have revealed that the diatom light-harvesting antenna fucoxanthin chlorophyll a/c-binding protein (FCPII) forms a tetramer and differs from the spinach antenna in terms of the number of protomers; however, the detailed molecular mechanism remains elusive. Herein, we report the physicochemical factors contributing to the characteristic light absorption of the diatom light-harvesting antenna based on spectral calculations using an exciton model. Spectral analysis reveals the significant contribution of unique fucoxanthin molecules (fucoxanthin-S) in FCPII to the diatom-specific spectrum, and further analysis determines their essential role in excitation-energy transfer to chlorophyll. It was revealed that the specificity of these fucoxanthin-S molecules is caused by the proximity between protomers associated with the tetramerization of FCPII. The findings of this study demonstrate that diatoms employ fucoxanthin-S to harvest energy under the ocean in the absence of long-wavelength sunlight and can provide significant information about the survival strategies of photosynthetic organisms to adjust to their living environment.

Sulfur-Bridged Cationic Diazulenomethenes: Formation of Charge-Segregated Assembly with High Charge-Carrier Mobility.

Sulfur-bridged cationic diazulenomethenes were synthesized and exhibited high stability even under basic conditions due to the delocalization of positive charge over the whole π-conjugated skeleton. As a result of the effective delocalization and the absence of orthogonally oriented bulky substituents, the cationic π-conjugated skeletons formed a π-stacked array with short interfacial distances. A derivative with SbF6- as a counter anion formed a charge-segregated assembly in the crystalline state, rather than the generally favored charge-by-charge arrangement of oppositely charged species based on electrostatic interactions. Theoretical calculations suggested that the destabilization caused by electrostatic repulsion between two positively charged π-conjugated skeletons is compensated by the dispersion forces. In addition, the counter anion SbF6- played a role in regulating the molecular alignment through F⋯H-C and F-S interactions, which resulted in the charge-segregated alignment of the cationic π-skeletons. This characteristic assembled structure gave rise to a high charge-carrier mobility of 1.7 cm2 V-1 s-1 as determined using flash-photolysis time-resolved microwave conductivity.

Quasi-degenerate extension of local N-electron valence state perturbation theory with pair-natural orbital method based on localized virtual molecular orbitals.

Chemical phenomena involving near-degenerate electronic states, such as conical intersections or avoided crossing, can be properly described using quasi-degenerate perturbation theory. This study proposed a highly scalable quasi-degenerate second-order N-electron valence state perturbation theory (QD-NEVPT2) using the local pair-natural orbital (PNO) method. Our recent study showed an efficient implementation of the PNO-based state-specific NEVPT2 method using orthonormal localized virtual molecular orbitals (LVMOs) as an intermediate local basis. This study derived the state-coupling (or off-diagonal) terms to implement QD-NEVPT2 in an alternative manner to enhance efficiency based on the internally contracted basis and PNO overlap matrices between different references. To facilitate further acceleration, a local resolution-of-the-identity (RI) three-index integral generation algorithm was developed using LMOs and LVMOs. Although the NEVPT2 theory is considered to be less susceptible to the intruder-state problem (ISP), this study revealed that it can easily suffer from ISP when calculating high-lying excited states. We ameliorated this instability using the imaginary level shift technique. The PNO-QD-NEVPT2 calculations were performed on small organic molecules for the 30 lowest-lying states, as well as photoisomerization involving the conical intersection of 1,1-dimethyldibenzo[b,f] silepin with a cis-stilbene skeleton. These calculations revealed that the PNO-QD-NEVPT2 method yielded negligible errors compared to the canonical QD-NEVPT2 results. Furthermore, we tested its applicability to a large photoisomerization system using the green fluorescent protein model and the ten-state calculation of the large transition metal complex, showcasing that off-diagonal elements can be evaluated at a relatively low cost.

Nonplanar Nanographene: A Hydrocarbon Hole-Transporting Material That Competes with Triarylamines.

Hole-transporting materials (HTMs) are essential for optoelectronic devices, such as organic light-emitting diodes (OLEDs), dye-sensitized solar cells, and perovskite solar cells. Triarylamines have been employed as HTMs since they were introduced in 1987. However, heteroatoms or side chains embedded in the core skeleton of triarylamines can cause thermal and chemical stability problems. Herein, we report that hexabenzo[a,c,fg,j,l,op]tetracene (HBT), a small nonplanar nanographene, functions as a hydrocarbon HTM with hole transport properties that match those of triarylamine-based HTMs. X-ray structural analysis and theoretical calculations revealed effective multidirectional orbital interactions and transfer integrals for HBT. In-depth experimental and theoretical analyses revealed that the nonplanarity-inducing annulative π-extension can achieve not only a stable amorphous state in bulk films, but also a higher increase in the highest occupied molecular orbital level than conventional linear or cyclic π-extension. Furthermore, an in-house manufactured HBT-based OLED exhibited excellent performance, featuring superior curves for current density-voltage, external quantum efficiency-luminance, and lifetime compared to those of representative triarylamine-based OLEDs. A notable improvement in device lifetime was observed for the HBT-based OLED, highlighting the advantages of the hydrocarbon HTM. This study demonstrates the immense potential of small nonplanar nanographenes for optoelectronic device applications.

In Silico Screening and Experimental Verification of Near-Infrared-Emissive Two-Boron-Doped Polycyclic Aromatic Hydrocarbons.

Embedding two boron atoms into a polycyclic aromatic hydrocarbon (PAH) leads to the formation of a neutral analogue that is isoelectronic to the corresponding dicationic PAH skeleton, which can significantly alter its electronic structure. Based on this concept, we explore herein the identification of near-infrared (NIR)-emissive PAHs with the aid of an in silico screening method. Using perylene as the PAH scaffold, we embedded two boron atoms and fused two thiophene rings to it. Based on this design concept, all possible structures (ca. 2500 entities) were generated using a comprehensive structure generator. Time-dependent DFT calculations were conducted on all these structures, and promising candidates were extracted based on the vertical excitation energy, transition dipole moment, and atomization energy per bond. One of the extracted dithieno-diboraperylene candidates was synthesized and indeed exhibited emission at 724 nm with a quantum yield of 0.40 in toluene, demonstrating the validity of this screening method. This modification was further applied to other PAHs, and a series of thienobora-modified PAHs was synthesized.

Zwitterionic Acridinium Amidate: A Nitrogen-Centered Radical Catalyst for Photoinduced Direct Hydrogen Atom Transfer.

The development of small organic molecules that can convert light energy into chemical energy to directly promote molecular transformation is of fundamental importance in chemical science. Herein, we report a zwitterionic acridinium amidate as a catalyst for the direct functionalization of aliphatic C-H bonds. This organic zwitterion absorbs visible light to generate the corresponding amidyl radical in the form of excited-state triplet diradical with prominent reactivity for hydrogen atom transfer to facilitate C-H alkylation with a high turnover number. The experimental and theoretical investigations revealed that the noncovalent interactions between the anionic amidate nitrogen and a pertinent hydrogen-bond donor, such as hexafluoroisopropanol, are crucial for ensuring the efficient generation of catalytically active species, thereby fully eliciting the distinct reactivity of the acridinium amidate as a photoinduced direct hydrogen atom transfer catalyst.

A quantum chemical study on the anti-SARS-CoV-2 activity of TMPRSS2 inhibitors.

Nafamostat and camostat are known to inhibit the spike protein-mediated fusion of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by forming a covalent bond with the human transmembrane serine protease 2 (TMPRSS2) enzyme. Previous experiments revealed that the TMPRSS2 inhibitory activity of nafamostat surpasses that of camostat, despite their structural similarities; however, the molecular mechanism of TMPRSS2 inhibition remains elusive. Herein, we report the energy profiles of the acylation reactions of nafamostat, camostat, and a nafamostat derivative by quantum chemical calculations using a combined molecular cluster and polarizable continuum model (PCM) approach. We further discuss the physicochemical relevance of their inhibitory activity in terms of thermodynamics and kinetics. Our analysis attributes the strong inhibitory activity of nafamostat to the formation of a stable acyl intermediate and its low activation energy during acylation with TMPRSS2. The proposed approach is also promising for elucidating the molecular mechanisms of other covalent drugs.

Algorithm for analytic nuclear energy gradients of state averaged DMRG-CASSCF theory with newly derived coupled-perturbed equations.

We present an algorithm for evaluating analytic nuclear energy gradients of the state-averaged density matrix renormalization group complete-active-space self-consistent field (SA-DMRG-CASSCF) theory based on the newly derived coupled-perturbed (CP) DMRG-CASSCF equations. The Lagrangian for the conventional SA-CASSCF analytic gradient theory is extended to the SA-DMRG-CASSCF variant that can fully consider a whole set of constraints on the parameters of multi-root canonical matrix product states formed at all the DMRG block configurations. An efficient algorithm to solve the CP-DMRG-CASSCF equations for determining the multipliers was developed. The complexity of the resultant analytic gradient algorithm is overall the same as that of the unperturbed SA-DMRG-CASSCF algorithm. In addition, a reduced-scaling approach was developed to directly compute the SA reduced density matrices (SA-RDMs) and their perturbed ones without calculating separate state-specific RDMs. As part of our implementation scheme, we neglect the term associated with the constraint on the active orbitals in terms of the active-active rotation in the Lagrangian. Thus, errors from the true analytic gradients may be caused in this scheme. The proposed gradient algorithm was tested with the spin-adapted implementation by checking how accurately the computed analytic energy gradients reproduce numerical gradients of the SA-DMRG-CASSCF energies using a common number of renormalized bases. The illustrative applications show that the errors are sufficiently small when using a typical number of the renormalized bases, which is required to attain adequate accuracy in DMRG's total energies.

Local N-electron valence state perturbation theory using pair-natural orbitals based on localized virtual molecular orbitals.

Second-order N-electron valence state perturbation theory (NEVPT2) is an exactly size-consistent and intruder-state-free multi-reference theory. To accelerate the NEVPT2 computation, Guo and Neese combined it with the local pair-natural orbital (PNO) method using the projected atomic orbitals (PAOs) as the underlying local basis [Guo et al., J. Chem. Phys. 144, 094111 (2016)]. In this paper, we report the further development of the PNO-NEVPT2 method using the orthonormal and non-redundant localized virtual molecular orbitals (LVMOs) instead of PAOs. The LVMOs were previously considered to perform relatively poor compared to PAOs because the resulting orbital domains were unacceptably large. Our prior work, however, showed that this drawback can be remedied by re-forming the domain construction scheme using differential overlap integrals [Saitow et al., J. Chem. Phys. 157, 084101 (2022)]. In this work, we develop further refinements to enhance the feasibility of using LVMOs. We first developed a two-level semi-local approach for screening out so-called weak-pairs to select or truncate the pairs for PNO constructions more flexibly. As a refinement specific to the Pipek-Mezey localization for LVMOs, we introduced an iterative scheme to truncate the Givens rotations using varying thresholds. We assessed the LVMO-based PNO-NEVPT2 method through benchmark calculations for linear phenyl alkanes, which demonstrate that it performs comparably well relative to the PAO-based approach. In addition, we evaluated the Co-C bond dissociation energies for the cobalamin derivatives composed of 200 or more atoms, which confirms that the LVMO-based method can recover more than 99.85% of the canonical NEVPT2 correlation energy.

Machine-Learning Classification for the Prediction of Catalytic Activity of Organic Photosensitizers in the Nickel(II)-Salt-Induced Synthesis of Phenols.

Catalytic systems using a small amount of organic photosensitizer for the activation of an inorganic (on-demand ligand-free) nickel(II) salt represent a cost-effective method for cross-coupling reactions, while C(sp2 )-O bond formation remains less developed. Herein, we report a strategy for the synthesis of phenols with a nickel(II) salt and an organic photosensitizer, which was identified via an investigation into the catalytic activity of 60 organic photosensitizers consisting of various electron donor and acceptor moieties. To examine the effect of multiple intractable parameters on the catalytic activity of photosensitizers, machine-learning (ML) models were developed, wherein we embedded descriptors representing their physical and structural properties, which were obtained from DFT calculations and RDKit, respectively. The study clarified that integrating both DFT- and RDKit-derived descriptors in ML models balances higher "precision" and "recall" across a wide range of search space relative to using only one of the two descriptor sets.

Thermally Stable Array of Discrete C60s on a Two-Dimensional Crystalline Adlayer of Macrocycles both in Vacuo and under Ambient Pressure.

A periodic monolayer array of discrete C60s was generated on an atomically flat Au(111) surface with the aid of a template adlayer. The template was a two-dimensional (2D) array of molecular pits prepared on an Au(111) surface through 2D crystallization of shape-persistent macrocycles composed of four carbazole and four salphens/Ni-salphens with a 1 nm hollow. Scanning tunneling microscopy imaging under ultra-high vacuum revealed that the square-shaped macrocycles, with 1.5 nm sides, were arranged with a periodic spacing of approximately 4.0 nm on the Au(111) surface, where the orientation and periodicity of the macrocycles were dependent on their chemical structures. After sublimation of C60s onto the adlayer, a single C60 molecule was entrapped in each pit, and an ordered molecular array of C60s was attained with a pattern similar to that of the macrocycles. The periodic pattern of C60s on the surface was thermally stable up to approximately 200 °C, even under ambient pressure. Scanning tunneling spectroscopy suggested the existence of an electronic interaction between the C60s and the Au(111) surface that was influenced by the macrocycle template on the surface.

Impact of Hydrophobic/Hydrophilic Balance on Aggregation Pathways, Morphologies, and Excited-State Dynamics of Amphiphilic Diketopyrrolopyrrole Dyes in Aqueous Media.

We report the thermodynamic and kinetic aqueous self-assembly of a series of amide-functionalized dithienyldiketopyrrolopyrroles (TDPPs) that bear various hydrophilic oligoethylene glycol (OEG) and hydrophobic alkyl chains. Spectroscopic and microscopic studies showed that the TDPP-based amphiphiles with an octyl group form sheet-like aggregates with J-type exciton coupling. The effect of the alkyl chains on the aggregated structure and the internal molecular orientation was examined via computational studies combining MD simulations and TD-DFT calculations. Furthermore, solvent and thermal denaturation experiments provided a state diagram that indicates the formation of unexpected nanoparticles during the self-assembly into nanosheets when longer OEG side chains are introduced. A kinetic analysis revealed that the nanoparticles were obtained selectively as an on-pathway intermediate state toward the formation of thermodynamically controlled nanosheets. The metastable aggregates were used for seed-initiated supramolecular assembly, which allowed establishing control over the assembly kinetics and the aggregate size. The sheet-like aggregates prepared using the seeding method exhibited coherent vibration in the excited state, indicating a well-ordered orientation of the TDPP units. These results underline the significance of fine tuning of the hydrophobic/hydrophilic balance in the molecular design to kinetically control the assembly of amphiphilic π-conjugated molecules into two-dimensional nanostructures in aqueous media.

Structure-Function Study of a Novel Inhibitor of Cyclin-Dependent Kinase C in Arabidopsis.

The circadian clock, an internal time-keeping system with a period of about 24 h, coordinates many physiological processes with the day-night cycle. We previously demonstrated that BML-259 [N-(5-isopropyl-2-thiazolyl) phenylacetamide], a small molecule with mammal CYCLIN DEPENDENT KINASE 5 (CDK5)/CDK2 inhibition activity, lengthens Arabidopsis thaliana (Arabidopsis) circadian clock periods. BML-259 inhibits Arabidopsis CDKC kinase, which phosphorylates RNA polymerase II in the general transcriptional machinery. To accelerate our understanding of the inhibitory mechanism of BML-259 on CDKC, we performed structure-function studies of BML-259 using circadian period-lengthening activity as an estimation of CDKC inhibitor activity in vivo. The presence of a thiazole ring is essential for period-lengthening activity, whereas acetamide, isopropyl and phenyl groups can be modified without effect. BML-259 analog TT-539, a known mammal CDK5 inhibitor, did not lengthen the period nor did it inhibit Pol II phosphorylation. TT-361, an analog having a thiophenyl ring instead of a phenyl ring, possesses stronger period-lengthening activity and CDKC;2 inhibitory activity than BML-259. In silico ensemble docking calculations using Arabidopsis CDKC;2 obtained by a homology modeling indicated that the different binding conformations between these molecules and CDKC;2 explain the divergent activities of TT539 and TT361.

Phosphorylation of RNA Polymerase II by CDKC;2 Maintains the Arabidopsis Circadian Clock Period.

The circadian clock is an internal timekeeping system that governs about 24 h biological rhythms of a broad range of developmental and metabolic activities. The clocks in eukaryotes are thought to rely on lineage-specific transcriptional-translational feedback loops. However, the mechanisms underlying the basic transcriptional regulation events for clock function have not yet been fully explored. Here, through a combination of chemical biology and genetic approaches, we demonstrate that phosphorylation of RNA polymerase II by CYCLIN DEPENDENT KINASE C; 2 (CDKC;2) is required for maintaining the circadian period in Arabidopsis. Chemical screening identified BML-259, the inhibitor of mammalian CDK2/CDK5, as a compound lengthening the circadian period of Arabidopsis. Short-term BML-259 treatment resulted in decreased expression of most clock-associated genes. Development of a chemical probe followed by affinity proteomics revealed that BML-259 binds to CDKC;2. Loss-of-function mutations of cdkc;2 caused a long period phenotype. In vitro experiments demonstrated that the CDKC;2 immunocomplex phosphorylates the C-terminal domain of RNA polymerase II, and BML-259 inhibits this phosphorylation. Collectively, this study suggests that transcriptional activity maintained by CDKC;2 is required for proper period length, which is an essential feature of the circadian clock in Arabidopsis.

Substituent and Solvent Effects on the Photoisomerization of Cinnamate Derivatives: An XMS-CASPT2 Study.

Cinnamate derivatives show a variety of photo-induced reactions. Among them is trans-cis photoisomerization, which may involve the nonradiative decay (NRD) process. The extended multistate complete active space second-order perturbation (XMS-CASPT2) method was used in this study as a suitable theory for treating multireference electronic nature, which was frequently manifested in the photoisomerization process. The minimum energy paths of the trans-cis photoisomerization process of cinnamate derivatives were determined, and the activation energies were estimated using the resulting intrinsic reaction coordinate (IRC) paths. Natural orbital analysis revealed that the transition state's (TS) electronic structure is zwitterionic-like, elucidating the solvent and substituent effect on the energy barrier of photoisomerization paths. Furthermore, it was found that the charge on the pyramidalized carbon atom at the TS structure was strongly correlated with the activation energy barrier for the cinnamate derivatives. Thus, it seemingly provided a physical picture of the photoisomerization of cinnamates and was a good descriptor potentially applicable to molecular design for controlling the rate constant of the photoisomerization reaction.

A pentanuclear iron catalyst designed for water oxidation.

Although the oxidation of water is efficiently catalysed by the oxygen-evolving complex in photosystem II (refs 1 and 2), it remains one of the main bottlenecks when aiming for synthetic chemical fuel production powered by sunlight or electricity. Consequently, the development of active and stable water oxidation catalysts is crucial, with heterogeneous systems considered more suitable for practical use and their homogeneous counterparts more suitable for targeted, molecular-level design guided by mechanistic understanding. Research into the mechanism of water oxidation has resulted in a range of synthetic molecular catalysts, yet there remains much interest in systems that use abundant, inexpensive and environmentally benign metals such as iron (the most abundant transition metal in the Earth's crust and found in natural and synthetic oxidation catalysts). Water oxidation catalysts based on mononuclear iron complexes have been explored, but they often deactivate rapidly and exhibit relatively low activities. Here we report a pentanuclear iron complex that efficiently and robustly catalyses water oxidation with a turnover frequency of 1,900 per second, which is about three orders of magnitude larger than that of other iron-based catalysts. Electrochemical analysis confirms the redox flexibility of the system, characterized by six different oxidation states between Fe(II)5 and Fe(III)5; the Fe(III)5 state is active for oxidizing water. Quantum chemistry calculations indicate that the presence of adjacent active sites facilitates O-O bond formation with a reaction barrier of less than ten kilocalories per mole. Although the need for a high overpotential and the inability to operate in water-rich solutions limit the practicality of the present system, our findings clearly indicate that efficient water oxidation catalysts based on iron complexes can be created by ensuring that the system has redox flexibility and contains adjacent water-activation sites.

A local pair-natural orbital-based complete-active space perturbation theory using orthogonal localized virtual molecular orbitals.

The multireference second-order perturbation theory (CASPT2) is known to deliver a quantitative description of various complex electronic states. Despite its near-size-consistent nature, the applicability of the CASPT2 method to large, real-life systems is mostly hindered by large computational and storage costs for the two-external tensors, such as two-electron integrals, amplitudes, and residuum. To this end, Menezes and co-workers developed a reduced-scaling CASPT2 scheme by incorporating the local pair-natural orbital (PNO) representation of the many-body wave functions using non-orthonormal projected atomic orbitals (PAOs) into the CASPT theory [F. Menezes et al., J. Chem. Phys. 145, 124115 (2016)]. Alternatively, in this paper, we develop a new PNO-based CASPT2 scheme using the orthonormal localized virtual molecular orbitals (LVMOs) and assess its performance and accuracy in comparison with the conventional PAO-based counterpart. Albeit the compactness, the LVMOs were considered to perform somewhat poorly compared to PAOs in the local correlation framework because they caused enormously large orbital domains. In this work, we show that the size of LVMO domains can be rendered comparable to or even smaller than that of PAOs by the use of the differential overlap integrals for domain construction. Optimality of the MOs from the CASSCF treatment is a key to reducing the LVMO domain size for the multireference case. Due to the augmented Hessian-based localization algorithm, an additional computational cost for obtaining the LVMOs is relatively minor. We demonstrate that the LVMO-based PNO-CASPT2 method is routinely applicable to large, real-life molecules such as Menshutkin SN2 reaction in a single-walled carbon nanotube reaction field.

Stochastic evaluation of four-component relativistic second-order many-body perturbation energies: A potentially quadratic-scaling correlation method.

A second-order many-body perturbation correction to the relativistic Dirac-Hartree-Fock energy is evaluated stochastically by integrating 13-dimensional products of four-component spinors and Coulomb potentials. The integration in the real space of electron coordinates is carried out by the Monte Carlo (MC) method with the Metropolis sampling, whereas the MC integration in the imaginary-time domain is performed by the inverse-cumulative distribution function method. The computational cost to reach a given relative statistical error for spatially compact but heavy molecules is observed to be no worse than cubic and possibly quadratic with the number of electrons or basis functions. This is a vast improvement over the quintic scaling of the conventional, deterministic second-order many-body perturbation method. The algorithm is also easily and efficiently parallelized with 92% strong scalability going from 64 to 4096 processors.

Fluorescent Organic π-Radicals Stabilized with Boron: Featuring a SOMO-LUMO Electronic Transition.

We report on the fluorescence properties of a new class of emissive and stable π-radicals that contain a boron atom at a position distant from the radical center. A fully planarized derivative exhibited an intense red fluorescence with high fluorescence quantum yields (ΦF >0.67) even in polar solvents. To elucidate the origin of this phenomenon, we synthesized another boron-stabilized radical that contains a bulky aryl group on the boron atom. A comparison of these derivatives, as well as with conventional donor-π-acceptor (D-π-A)-type emissive π-radicals, unveiled several characteristic features in their photophysical properties. A theoretical analysis revealed that the SOMO-LUMO electronic transition generates an emissive D1 state. Unlike conventional D-π-A-type π-radicals, this state does not undergo significant structural relaxation. The boron-stabilized π-radicals demonstrated promising potential for organic light-emitting diodes as an emitting material.

Investigating the Nonradiative Decay Pathway in the Excited State of Silepin Derivatives: A Study with Second-Order Multireference Perturbation Wavefunction Theory.

The fluorescence quantum yield for fluorescent organic molecules is an important molecular property, and tuning it up is desired for various applications. For the computational estimation of the fluorescence quantum yield, the theoretical prediction of the nonradiative decay rate constant has become an attractive subject of study. The rate constant of thermally activated nonradiative decay is related to the activation energy in the photoreaction; thus, the accuracy and reliability of the excited-state potential energies in the quantum chemical computation are critical. In this study, we employed a second-order multireference perturbation wavefunction theory for studying the thermally activated decay via conical intersection (CI) of 1,1-dimethyldibenzo[b,f]silepin derivatives. The correlation between the computed activation energy to reach the CI geometry in the S1 state and the experimentally determined fluorescence quantum yield implied that silepins nonradiatively decay via the CI triggered by the twisting of the central C-C bond. Geometry optimization of the transition state using multireference perturbation theory drastically reduced the estimated activation energy. Our computation gave reasonable predictions of the activation free energies of photoexcited 1,1-dimethyldibenzo[b,f]silepin. The energy profiles and geometry optimizations using proper quantum chemical methods played a critical role in reliable estimation of the rate constant and fluorescence quantum yield.