Photochemical formation of the elusive Dewar isomers of aromatic systems: why are substituted azaborines different?

Photochemical reactions enabling efficient transformation of aromatic systems into energetic but stable non-aromatic isomers have a long history in organic chemistry. One recently discovered reaction in this realm is that where derivatives of 1,2-azaborine, a compound isoelectronic with benzene in which two adjacent C atoms are replaced by B and N atoms, form the non-hexagon Dewar isomer. Here, we report quantum-chemical calculations that explain both why 1,2-azaborine is intrinsically more reactive toward Dewar formation than benzene, and how suitable substitutions at the B and N atoms are able to increase the corresponding quantum yield. We find that Dewar formation from 1,2-azaborine is favored by a pronounced driving force that benzene lacks, and that a large improvement in quantum yield arises when the reaction of substituted 1,2-azaborines proceeds without involvement of an intermediary ground-state species. Overall, we report new insights into making photochemical use of the Dewar isomers of aromatic compounds.


Introduction
][4][5] Around the same time and later, these findings triggered the search for valence-bond isomers of other aromatic and heteroaromatic compounds, leading to the detection of Dewar isomers of pyridines, 7-10 naphthalenes 11 and anthracenes 12 as photolysis products of their corresponding parent compounds.
9][20] As for the possibility to pursue photochemical synthesis of Dewar species, quantum- This journal is © the Owner Societies 2024 chemical modelling has also shown that the photoisomerization reactions through which they are formed proceed in singlet excited states, 18,20 as originally suggested for Dewar benzene. 3hus, if such synthesis is attempted, addition of triplet sensitizers or removal of triplet quenchers is not required.However, any potential application of Dewar species that involves their photochemical formation is challenged by the fact that these processes generally exhibit rather low quantum yields.For example, for Dewar benzene and Dewar pyridine, the quantum yields are as low as 0.006 2 and 0.07, 7 respectively.Hence, for applications, other arenes or heteroarenes should likely be considered.
1,2-Dihydro-1,2-azaborine, from now on azaborine, is a wellknown compound 21,22 isoelectronic with benzene in which two adjacent C atoms are replaced by B and N atoms (see Fig. 1b).4][25] Then, in 2009, Dixon and Liu and their co-workers presented the first synthetic route giving access to the pristine compound, which turned out to be stable (consistent with substantial aromatic character) and showed no appreciable degradation when subjected to high-temperature conditions for several days. 26hese researchers also compared the electronic properties of azaborine to those of benzene and borazine (another benzeneisoelectronic compound in which all six C atoms are replaced by alternating B and N atoms) through measurements of UV-vis absorption spectra in pentane.Relative to these systems, the absorption of azaborine is strongly red-shifted, with the most intense bands peaking at 269 (azaborine), 208 (benzene) and 203 nm (borazine), respectively. 26hortly thereafter, in 2012, the first photolysis experiments of azaborine were reported. 27Using low-temperature Ne/Ar/Xe matrices and a low-pressure Hg lamp (254 nm), conversion of azaborine into a photoproduct was observed within 105 minutes. 27y comparing the IR spectrum of the photoproduct with calculated IR spectra of potential structures, the photoproduct was identified as the Dewar isomer of azaborine featuring a single bond between the C 3 and C 6 atoms (those adjacent to the B and N atoms, see Fig. 1b). 27A few years later, Liu and Bettinger and their co-workers presented two azaborine derivatives, bearing a silyl group at the N atom and a mesityl group or Cl atom at the B atom, capable of achieving high quantum yields (up to 0.46) for Dewar formation. 28Combined with the finding that both of the corresponding Dewar isomers show sizable barriers for the thermal ring-opening reaction back to the parent isomer, these results point to the possibility of using azaborines for future applications in solar-energy storage. 280][31][32] However, despite the existence of several promising MOST systems (based on, e.g., azobenzene, norbornadiene and dihydroazulene photoswitches), 30 it has proven very difficult for any single one of them to simultaneously meet all photophysical and thermochemical requirements for viable applications in solar-energy storage. 30,31,33,34Besides showing a high quantum yield for the photoisomerization that produces the energystoring isomer and a large barrier for the thermal backreaction (so as to enable long storage times), other desirable characteristics include a high energy density, UV-vis absorption in the violet-blue regime, and minimal spectral overlap between the parent and energy-storing isomers. 31,32iven the fundamental importance of an efficient photoisomerization step for the utility of a MOST system, optimizing azaborines for such applications will likely require a better understanding of their photochemistry than currently available.Based on quantum-chemical calculations, it has been found that the photoisomerization of azaborine in the first singlet excited state (S 1 ) that produces the C 3 -C 6 bond is more favorable than photoisomerizations forming a N-C 4 or B-C 5 bond. 35urthermore, it has been shown that azaborine features an easily accessible conical intersection (CI) between the S 1 and ground S 0 states, which is connected to both the parent species and a prefulvene-like structure containing a B-C 6 bond. 36The latter structure is quite unstable, and can evolve through small energy barriers to either the parent species or the C 3 -C 6 Dewar isomer. 36ecently, these results were corroborated by means of nonadiabatic molecular dynamics (NAMD) simulations, from which it was concluded that the prefulvene-like structure is indeed a key mechanistic element for the photochemistry of azaborine, mediating both the regeneration of the parent species and the formation of the Dewar isomer. 37However, the quantum yield for Dewar formation extractable from these simulations appears much too low (only a few percent 37 ) compared to the experimental observations. 27espite the many insights provided by the aforementioned studies, a more thorough appreciation of the unique photochemistry of azaborines and their potential for future applications will require resolving a number of fundamental issues, which are yet to be explored.Firstly, it is crucial to understand exactly how and why the photochemistry of azaborine is different from that of archetypal benzene with respect to Dewar formation.Secondly, it is equally important to explain the origin Fig. 2 Molecules studied in this work (from left to right and from top to bottom): benzene, azaborine, 1-silyl-2-chlorobenzene and 1-silyl-2chloroazaborine.The atom-numbering schemes displayed are used throughout the paper.1-Silyl-2-chloroazaborine is very similar to one of the substituted azaborines studied by Liu and Bettinger, 28 except that it harbors the simplest possible silyl group (SiH 3 ) instead of a tertbutyldimethylsilyl group for computational expedience.We do not expect the neglected alkyl groups to have any significant impact on the formation of the Dewar isomer.
of the favorable (for Dewar azaborine formation) substituent effects observed by Liu and Bettinger in 2018. 28As a first step towards filling these gaps, we here present a detailed comparison of the photochemistry of benzene, azaborine and their 1-silyl-2chloro-substituted counterparts (see Fig. 2) in the S 1 state using both static quantum-chemical calculations and semi-classical NAMD simulations 38,39 based on time-dependent density functional theory (TD-DFT). 40Moreover, the TD-DFT modelling is complemented by calculations based on complete active space self-consistent field (CASSCF) theory. 41Through these endeavors, we are able to identify key electronic and structural features that make azaborine more prone than benzene to form the Dewar isomer, and allow substituted azaborines to attain much higher quantum yields for Dewar formation than the parent species.

UV-vis absorption spectra
In order to assess the accuracy of the adopted DFT-based computational methodology relative to both previous modelling studies and experimental data, our first task was to calculate vertical transition energies for the five lowest singlet excited states of each of the four molecules in Fig. 2. As can be seen from Table 1, all transitions lie in the UV-C regime (o280 nm, 44.43 eV), with the S 0 -S 1 transition for the two azaborines being red-shifted by ca.0.35 eV relative to the S 0 -S 1 transition for their benzene counterparts.Furthermore, from the corresponding oscillator strengths, it is clear that whereas the S 0 -S 1 transition is strongly allowed for the C s -symmetric azaborines, it is forbidden for the benzenes (D 6h and C s ).Notably, the S 0 -S 1 energy of 5.18 eV obtained for the plain azaborine agrees very well with previous high-level, ab initio calculations, which have yielded values of 5.15 (using the multireference configuration interaction method 36 ) and 5.16 eV (using the equation-of-motion coupled cluster singles and doubles method 37 ) for this transition.
Next, in order to also evaluate the computational methodology relative to experimental data, 26 we calculated the UV-vis absorption spectra of the four molecules by describing their nuclear motion around the respective S 0 minimum through Wigner sampling from harmonic vibrational frequencies. 42,43The resulting spectra, shown in Fig. 3, reveal that the breaking of planar symmetry allows some forbidden transitions to take place.For benzene and azaborine, the absorption maxima occur at 208 and 246 nm, respectively, which are reasonably close to the positions of 208 and 269 nm recorded experimentally in a pentane solvent 26 (notwithstanding the fact that we have here chosen to not include solvent effects in the modelling).For 1-silyl-2-chloroazaborine, the absorption maximum is further red-shifted to 262 nm, which reflects the observation in Table 1 that the silyl group and Cl atom exert a red-shifting effect on each of the five lowest excited states.Notably, owing to the molecular motion, 1-silyl-2-chloroazaborine is even predicted to exhibit some residual absorption in the UV-B regime (280-315 nm).

NAMD simulations
As noted in the Introduction, a central goal of this work is to gain an understanding of why the S 1 state of azaborines is more reactive towards Dewar formation than the S 1 state of benzene.Thus, we are concerned with identifying key features of this state, and its interplay with the ground S 0 state, that promote Dewar formation.For this reason, all modelling performed to address these issues focused exclusively on the S 1 and S 0 states, so as to ensure a balanced comparison of the four systems under study.Although such an approach does not account for the proposal that Dewar benzene might be formed from the second singlet excited state (S 2 ), 3 this incongruence is of no consequence given the very low quantum yield (0.006) with which Dewar benzene is produced. 2That this approach is sound is also attested by the many previous studies devoted Table 1 Calculated vertical transition energies (DE and l, in eV and nm, respectively) and oscillator strengths (f) for the molecules studied in this work  This journal is © the Owner Societies 2024 to investigating the photochemical reactivity of aromatic compounds in the S 1 state. 44aving tested the performance of the DFT-based computational methodology, we then proceeded to model the photodynamics of the four systems by means of NAMD simulations.For each system, the simulations were run for ca. 1 ps with 100 different initial nuclear geometries and velocities starting from the vertically excited S 1 Franck-Condon (FC) points of these geometries.At any given time during the simulations, the classical populations P 0 and P 1 of the S 0 and S 1 states were obtained (following the usual convention) as the proportions of trajectories evolving in the respective state at that particular time.From these populations, the S 1 lifetime -denoted t 1/2was estimated as the time it takes before the two states become equally populated.
For each system, Fig. 4 plots the time evolution of the classical populations of the S 0 and S 1 states predicted by the NAMD simulations.Interestingly, from these results, it can be inferred that the photophysical behavior of the azaborines is rather different from that of their benzene counterparts.As for the latter, only 9 (benzene) and 13% (1-silyl-2-chlorobenzene) of the trajectories have decayed back to the S 0 state after 1 ps of simulation.In sharp contrast, the former show full (100%, azaborine) or almost full (96%, 1-silyl-2-chloroazaborine) S 0 recovery within the same time frame.Furthermore, from the observation of a rather substantial difference in the S 1 lifetimes of azaborine (111 fs) and 1-silyl-2-chloroazaborine (305 fs), the photophysical behavior of the plain azaborine appears more controllable through functionalization than what is the case for benzene.
Next, we analyzed the distributions of species formed upon internal conversion (IC) to the S 0 state.The corresponding results are shown as histograms in Fig. 5, in which each molecule is represented by a specific color.After 1 ps, all 9 trajectories (= 9% of the trajectories) for benzene (red) and all 13 trajectories for 1-silyl-2-chlorobenzene (yellow) that did decay to the S 0 state have reformed the parent benzene and 1-silyl-2-chlorobenzene species, respectively, without producing Dewar or other photoproducts.Loosely, the complete absence of Dewar photoproducts in the simulations of these compounds is consistent with the minuscule quantum yield assigned to Dewar benzene in benzene photolysis. 2 Contrarily, among the 96 trajectories for 1-silyl-2-chloroazaborine that decayed to the S 0 state, 41 of them formed the corresponding Dewar photoproduct.Notably, this result corroborates the finding by Liu and Bettinger that 1,2-disubstituted azaborines can achieve rather high quantum yields for Dewar formation (e.g., 0.46 for an azaborine bearing a tert-butyldimethylsilyl group at the N atom and a mesityl group at the B atom). 28mong the remaining 96 À 41 = 55 trajectories, 52 reformed the parent 1-silyl-2-chloroazaborine species, whereas only 3 formed another photoproduct (a dihydroborole derivative) than the Dewar one.
Turning to the unsubstituted azaborine in Fig. 5, among the 100 trajectories (= 100% of the trajectories) that underwent comparatively faster (than in 1-silyl-2-chloroazaborine) S 0 recovery, 19 of them formed the Dewar photoproduct.Thereby, the associated quantum yield can be estimated to be about twofold lower than that for 1-silyl-2-chloroazaborine.Moreover, for the unsubstituted azaborine, the proportion of alternative (to Dewar) photoproducts produced is higher (8 vs. 3% for 1silyl-2-chloroazaborine), with four such photoproducts detected.Taken together, the NAMD results in Fig. 5 support two main conclusions that are consistent with experimental  observations. 27,28First, whereas photoexcitation to the S 1 state does not promote formation of Dewar benzene, photoisomerization of azaborine into Dewar azaborine is possible in this state and occurs with non-negligible quantum yield. 27Second, this quantum yield can be improved by introducing suitable substituents at the B and N atoms. 28However, the NAMD simulations provide no rationale for these findings.Therefore, in the next section, we will present electronic and structural factors, identified from S 0 and S 1 potential energy surfaces (PESs) calculated with static DFT methods, that do explain these results.

Potential energy surfaces
Fig. 6 summarizes the regions of the S 0 and S 1 PESs of benzene, azaborine, 1-silyl-2-chlorobenzene and 1-silyl-2-chloroazaborine relevant for the photoisomerization reactions into their Dewar forms.Given that Dewar formation is indeed the central focus of this work, Fig. 6 does not include any of the alternative reaction pathways that are possible for these molecules following IC to the S 0 state, as evidenced by the appearance of other (than Dewar) photoproducts in the NAMD simulations.Nonetheless, for the case of the unsubstituted azaborine, we will return to a discussion of such pathways below.Besides Fig. 6, results from complementary CASSCF-based calculations (which corroborate the results in Fig. 6) are presented in Section S1 of the ESI.† Starting with benzene and 1-silyl-2-chlorobenzene, their PESs are qualitatively very similar (Fig. 6, left column).Following photoexcitation to the respective S 1 FC point, in both systems subsequent relaxation populates a planar S 1 minimum slightly below the FC point (0.06 eV below for benzene and 0.09 eV below for 1-silyl-2-chlorobenzene).Furthermore, in both systems, reaching the S 1 /S 0 CI that acts as a funnel for the formation of the S 0 Dewar isomer from the S 1 minimum is an energetically uphill process.In fact, for benzene/1-silyl-2-chlorobenzene, this funnel lies 0.17/0.34eV above the S 1 minimum and 0.11/0.25 eV above the S 1 FC point.For this reason, photoexcitation to the S 1 state does not seem to provide sufficient energy for accessing these structures, which explains the observation from the NAMD simulations that the benzenes are unlikely to form the S 0 Dewar isomer from this state (see Fig. 5).The actual geometries of the CIs (including those for azaborine and 1-silyl-2-chloroazaborine) are shown in Section S2 of the ESI.† While the reaction coordinates leading to the CIs are multidimensional, which is generally the case in photochemistry 45 and explains why no particular reaction coordinate is specified in Fig. 6, Section S2 of the ESI † also presents some data for assessing relevant structural changes along the respective coordinate.
Continuing with the azaborines, their PESs (Fig. 6, right column) show a major difference relative to the benzenes in that the S 1 /S 0 CI now lies well below the S 1 FC point (1.03 eV below for azaborine and 0.72 eV below for 1-silyl-2-chloroazaborine), and thus is much more easily accessible.This difference provides a rationale why the likelihood of Dewar formation in the NAMD simulations is so much higher for the azaborines than for the benzenes (see Fig. 5).As for the unsubstituted azaborine, the S 1 /S 0 CI is reached in a barrierless fashion, which is consistent with its short S 1 lifetime of 111 fs (see Fig. 4).This CI is connected to a prefulvene-like S 0 minimum (at 3.69 eV) in one direction, and to the S 0 parent species in the other.The prefulvene-like minimum, which was also found in previous studies, 36,37 can evolve through a small energy barrier to either the S 0 Dewar isomer (barrier of 0.06 eV) or the S 0 parent species (0.12 eV).Accordingly, this minimum is a sort of pivotal point steering the competition between these two pathways.Taken together, these results identify a route through which formation of the S 0 Dewar isomer from the S 1 state of an unsubstituted azaborine is an energetically favorable process.Altogether, the S 0 dynamics following IC at the S 1 /S 0 CI has at least three different components (formation of the Dewar isomer through the prefulvenelike minimum, formation of the parent species through the prefulvene-like minimum, and direct formation of the parent species), which are further discussed in terms of calculated NAMD trajectories in Section S3 of the ESI.† As for the 1-silyl-2-chloroazaborine, its substituents have two effects on the calculated PESs in Fig. 6 relative to those of the unsubstituted azaborine.The first effect is the appearance of a minimum on the S 1 PES, 0.33 eV below the FC point.Featuring a slight ring puckering at the N atom, this minimum is not equivalent to the aforementioned planar minima on the S 1 PESs of benzene and 1-silyl-2-chlorobenzene.Furthermore, being separated from the S 1 /S 0 CI by an energy barrier of just 0.02 eV, this minimum is very shallow and does not appear to have a negative influence on Dewar formation, as judged by the

PCCP Paper
fact that the number of NAMD trajectories forming the Dewar isomer are about twice as many for 1-silyl-2-chloroazaborine than for azaborine (see Fig. 5).However, albeit tiny, this barrier may be responsible for the longer S 1 lifetime (305 vs. 111 fs) shown by 1-silyl-2-chloroazaborine. The second effect attributable to the substituents of 1-silyl-2-chloroazaborine is the disappearance of the prefulvene-like minimum on the S 0 PES (in fact, all attempts to locate such a minimum produced the Dewar isomer).In the absence of such a minimum, the S 1 /S 0 CI is instead connected directly to the Dewar isomer in one reaction channel, and to the parent species in the other.This direct route to the Dewar isomer, combined with the fact that the absence of a prefulvene-like minimum also eliminates an alternative reaction channel back to the parent species, offers a simple explanation why 1-silyl-2-chloroazaborine shows a much higher quantum yield for Dewar formation than the unsubstituted azaborine (see Fig. 5).
Given that our work focuses on photoconversion of benzene, azaborine, 1-silyl-2-chlorobenzene and 1-silyl-2-chloroazaborine into their Dewar isomers, the PESs in Fig. 6 do not include any of the alternative reaction pathways that are possible for these molecules following IC to the S 0 state.However, for the case of the unsubstituted azaborine, Fig. 7 illustrates the regions of the S 0 PES that, starting from the prefulvene-like minimum (labeled I1), describe such pathways for the formation of alternative (to Dewar) photoproducts detected by the NAMD simulations.Here, for reference, the pathway leading to the Dewar isomer (labeled P1) and the pathway leading back to the parent azaborine (labeled S 0 FC), which were shown already in Fig. 6, are displayed in black color, whereas the alternative pathways are displayed in red, green and blue colors, respectively.Notably, the I1 -P1 reaction is indeed the kinetically preferred pathway, having an energy barrier that is smaller (0.06 eV) than that for the I1 -S 0 FC (0.12 eV) process and those for all other competing pathways.Furthermore, the high (0.95 eV) barrier for the P1 -I1 back-reaction ensures that the Dewar isomer is stable.As a complement to Fig. 7, the optimized geometries of the potential-energy minima along the different pathways are shown in Section S4 of the ESI.† Besides the I1 -P1 and I1 -S 0 FC reactions in Fig. 7, a third pathway is that producing the P2 5,6-E-azaborine, which is displayed in red color in Fig. 7 and has a barrier of 0.32 eV.However, even if formed, this structure appears to be rather unstable, reverting back to I1 via a very small barrier of 0.05 eV.A fourth pathway is the formation of the P3 dihydroborole, which is shown in green color in Fig. 7 and proceeds in two steps.In the first, a N-C 5 bond is formed to produce the benzvalene-like I2 intermediate through a barrier of 0.43 eV.This intermediate is then cleaved to yield P3 in a nearly barrierless (0.01 eV) process.The fifth and final pathway considered is the four-step formation of the P4 dihydroborete represented in blue color in Fig. 7.In the first step, the C 5 -C 6 bond of I1 is cleaved to form the singletbiradical I3 intermediate via a very large barrier of 2.84 eV.Rotation of the B-C 3 bond and cleavage of the B-N bond subsequently produce, through a barrier of 0.15 eV, the linear singlet-biradical I4 intermediate, which evolves to the closedshell I5 intermediate by means of hydrogen migration (barrier of 0.44 eV).Finally, the P4 dihydroborete is furnished by an electrocyclic ring-closure reaction (barrier of 0.13 eV).However, because of the considerable barrier for the first step, the likelihood that the P4 species is formed should be very small, which is also what the NAMD simulations in Fig. 5 of the main text predict (1%).Altogether, then, the results in Fig. 7 suggest that the two predominant reaction channels originating from the I1 prefulvene-like minimum are those producing the Dewar isomer (I1 -P1) or the parent azaborine (I1 -S 0 FC), which, again, is consistent with the NAMD results in Fig. 5 (indeed, 92% of the NAMD trajectories form either P1 or S 0 FC).Finally, we also sought to provide a further explanation for the difference between the azaborines and the benzenes in their reactivity towards Dewar formation rooted in the principle of least motion.First proposed by Rice and Teller in 1938, this principle asserts that ''there shall be least change in atomic position and least change in electronic configuration during an elementary reaction''. 460][51][52][53][54] Here, as we have explored one and the same reaction, the changes in atomic positions along the pathways mapped out by the calculated PESs are very similar.For this reason, we instead focused on changes in electron density.Specifically, using densities calculated with TD-DFT, we evaluated the changes in S 1 natural bond orbital (NBO) charges 55 at the atomic positions along the S 1 PESs of the four systems, as they evolve from the FC point to the CI.Using the corresponding NBO charges at the FC point as reference, these results are presented in Fig. 8, in which a red color indicates an increase in electron density (i.e., a decrease in atomic charge) and a green color indicates a decrease in electron density (i.e., an increase in atomic charge).
Interestingly, it can be seen that the changes in NBO charges are very different in the two groups of molecules.For the azaborines, the changes are largely monotonous, as exemplified by the decrease in electron density at the C 4 and C 6 positions (these positions go from being black/dark-green-colored to being bright-green-colored as the CI is approached), and the increase in electron density at the C 5 position (this position goes from being black/dark-red-colored to being bright-red-colored as the CI is approached).In striking contrast, for the benzenes, the changes are much more variable and show no uniform trends, with some positions exhibiting a decrease/increase in electron density at one stage of the reaction and an increase/decrease at the next.Qualitatively, this means that the limited accessibility to the CI mediating Dewar formation experienced by the benzenes (as compared to the azaborines, see Fig. 6) can indeed be understood in terms of the principle of least motion.

Conclusions
In summary, using a combination of static quantum-chemical calculations and semi-classical NAMD simulations, we have investigated the S 1 photochemistry of benzene, azaborine and their 1-silyl-2-chloro-substituted counterparts, with two particular goals.The first is to explain the experimental observation that azaborine, as opposed to benzene, is able to form the Dewar isomer with a non-negligible quantum yield. 27Notably, this observation is reproduced by the NAMD simulations in that 19 (azaborine) vs. 0% (benzene) of the calculated trajectories furnish the Dewar species.The second goal is to, similarly, rationalize the experimental observation that the quantum yield can be increased through suitable substitutions at the B and N atoms of azaborine, 28 which the NAMD simulations also account for by showing that the fraction of trajectories forming the Dewar isomer is appreciably larger for 1-silyl-2-chloroazaborine (41%) than for azaborine (19%).
As for the first goal, the higher propensity of azaborine to form the Dewar isomer can be explained by the finding that the associated S 1 /S 0 CI occurs well below (by 1.03 eV) the S 1 FC point of the parent species, which provides a considerable driving force for reaching the former structure upon photoexcitation.For benzene, on the other hand, the corresponding S 1 /S 0 CI lies above (by 0.11 eV) the S 1 FC point.As for the second goal, the improvement in quantum yield achieved by 1,2disubstituted azaborines can be rationalized by the absence of the prefulvene-like minimum that is present on the S 0 PES of the unsubstituted azaborine.Rather than being linked to such a minimum, the S 1 /S 0 CI is instead connected directly to the Dewar isomer, which is clearly advantageous for the formation of this species.Furthermore, the absence of this minimum also eliminates an undesirable reaction channel back to the parent 1,2-disubstituted azaborine.All in all, these results help identify strategies for tailoring the efficiency of a well known photochemical reaction for possible applications in solarenergy storage.

Computational details
S 0 and S 1 PESs were mapped by performing DFT and TD-DFT calculations, respectively, with the CAM-B3LYP range-separated Fig. 8 Changes in S 1 NBO charges along the S 1 PESs in Fig. 6 (below each structure is given its electronic energy in eV relative to the corresponding S 0 global minimum).A red color indicates an increase in electron density (i.e., a decrease in atomic charge) and a green color indicates a decrease in electron density (i.e., an increase in atomic charge) at that position relative to the electron density at the S 1 FC point.As the S 1 PES of azaborine does not feature any minimum or TS (see Fig. 6), structures obtained by replacing the silyl group and Cl atom at the minimum and the TS on the S 1 PES of 1-silyl-2-chloroazaborine with two H atoms (and reoptimizing just the positions of these H atoms) were considered for the corresponding calculations.

PCCP Paper
Open hybrid density functional. 56While potential-energy minima, TSs and minimum-energy S 1 /S 0 CIs were calculated with the cc-pVDZ basis set, 57 more accurate electronic energies of the corresponding structures were obtained through single-point calculations with the larger cc-pVTZ basis set. 57By calculating harmonic vibrational frequencies at the CAM-B3LYP/cc-pVDZ level of theory, all potential-energy minima and TSs were characterized as stationary points having only real frequencies or one imaginary frequency along the relevant normal mode, respectively.Moreover, in order to further ensure that the located TSs do indeed connect the associated reactant and product species, intrinsic reaction coordinate calculations 58 were also carried out at this level of theory.As a means to verify the TD-DFT results, single-point calculations were additionally performed with the multistate complete active space second-order perturbation theory (MS-CASPT2) 59,60 method.These calculations also employed the cc-pVTZ basis set, and were based on CASSCF 41 reference wave functions optimized in a state-average fashion, utilizing equal weights for the S 0 and S 1 states.For each of the four molecules studied, the active space for the CASSCF calculations was chosen so as to include the full p-system (i.e., 6 electrons in 6 orbitals).Vertical transition energies (at S 0 molecular geometries) and their oscillator strengths were calculated at the TD-CAM-B3LYP/ cc-pVTZ level of theory.In order to model UV-vis absorption spectra, for each molecule, 100 different molecular geometries around the respective S 0 minimum were generated by means of Wigner sampling.Accordingly, as thoroughly described by Barbatti and Sen, 43 the nuclear phase-space is assumed to follow a Wigner distribution 42 of S 0 harmonic vibrational frequencies, from which normal modes can be sampled and converted to Cartesian coordinates.For each of the 100 sampled sets of nuclear geometries, vertical transition energies and oscillator strengths were calculated for the five lowest singlet excited states.Based on these results, the UV-vis absorption spectrum for each individual geometry was then derived by convoluting the transitions with Gaussian functions having a half-width of 0.05 eV.Finally, by summing the resulting spectra for all geometries, the spectra shown in Fig. 3 were obtained.
NAMD simulations were carried out at the TD-CAM-B3LYP/ cc-pVDZ level of theory (CAM-B3LYP/cc-pVDZ in the S 0 state) using Tully's fewest-switches surface hopping algorithm 61 as implemented in TURBOMOLE. 62,63For each system, the simulations were started from the S 1 FC points of the 100 different nuclear geometries obtained through the aforementioned Wigner sampling.For each of these geometries, initial nuclear velocities were generated randomly from a Maxwell-Boltzmann distribution at 298 K.For the classical propagation of the nuclei in Tully's method, the leapfrog Verlet algorithm 64 was used to integrate Newton's equations of motion with a fixed time step of 40 a.u.(ca. 1 fs).Propagating the nuclear trajectories through 1024 steps, the total simulation time was ca. 1 ps.For the quantum mechanical description of the electrons in Tully's method, full technical details are provided in the original literature [61][62][63] or in recent applications of the method. 65,66e CASSCF and MS-CASPT2 calculations were performed with the OpenMolcas 19.11 suites of programs. 67The NAMD simulations were carried out with the TURBOMOLE 7.5 suites of programs. 63CIs were optimized with a penalty-function approach 68 as implemented in the SHARC 2.0 package 69 interfaced with the ORCA 4.2.1 70 and Gaussian 16 71 packages.All other calculations were done with Gaussian 16.

Fig. 3
Fig. 3 Calculated UV-vis absorption spectra of the molecules studied in this work.Intensities are given relative to that of the main band of azaborine (whose intensity is set to 1).

Fig. 5
Fig.5Photoproduct distributions after 1 ps of NAMD simulations.The photoproduct structures represented by the different bins are given below the bins (the leftmost bin is for the parent species, the second leftmost bin is for the Dewar photoproduct, and the other bins are for alternative photoproducts).For each system, the percentage of trajectories that remain in the S 1 state is obtained by subtracting the sum of all photoproduct percentages from 100%.

Fig. 6
Fig. 6 Calculated S 0 (green) and S 1 (blue) PESs of benzene (top left), azaborine (top right), 1-silyl-2-chlorobenzene (bottom left) and 1-silyl-2chloroazaborine (bottom right), with dashed vertical lines indicating photoexcitation to the respective S 1 FC point and local maxima corresponding to transition structures (TSs).S 1 /S 0 CIs are indicated with hourglass symbols.For each system, electronic energies of the relevant structures (stationary points, the S 1 FC point, and the S 1 /S 0 CI) are given in eV relative to the corresponding S 0 global minimum.

Fig. 7
Fig. 7 Calculated S 0 PES of azaborine comparing different reaction pathways starting from the I1 prefulvene-like minimum formed upon IC at the S 1 /S 0 CI.Electronic energies of stationary points (minima and TSs) are given in eV relative to the S 0 global minimum (labeled S 0 FC).