The mechanism behind the photochromism and photomagnetism of type II biindenylidenediones: multiconfigurational, perturbative and density functional theory studies.

Biindenylidenediones (BIDs) are a family of compounds that have been studied for a relatively short time. The crystals of these compounds are yellowish, and become purplish when they are irradiated and return back to their original color slowly in the dark or quickly when they are heated up. BIDs can be classified into different subfamilies depending on the nature of their substituents. BID-II crystals show a thermally dependent electron paramagnetic resonance (EPR) signal that is a characteristic of chemical species with unpaired electrons. These properties make BIDs very attractive for industrial applications but the mechanisms responsible for their photochromism and photomagnetism are still under debate. In this article, a computational study focused on the BID-II subfamily is presented. A variety of multiconfigurational methods (CASSCF, CASPT2 and IDDCI) have been used to study exhaustively the topography of the potential energy surfaces of the lower electronic states of a single BID molecule. Methods based on density functional theory (DFT) were then used to model the most important structures in a periodic crystal system. Our results suggest that δ-hydrogen abstraction could explain the observed experimental phenomena. After the initial excitation to the 1ππ* state, non-symmetric nπ* minima are populated, which are adiabatically connected to the photoproduct zone through a barrier along the reaction coordinate. Based on our set of results, we propose that an epoxide constitutes the most stable and accessible photoproduct preceded by the population of a triplet biradical of πOπ* nature which has only small geometrical differences in comparison with the reactant. The spin-orbit coupling indicates that the EPR signal arises due to the population of a low energy triplet through a thermally accessible intersystem crossing in the photoproduct zone.

Theoretical studies on the energy structures and optical properties of copper cysteamine - a novel sensitizer.

Copper cysteamine (Cu-Cy) is a new type of photosensitizer, which can be activated not only by ultraviolet light, but also by X-rays, microwaves and ultrasound to generate reactive oxygen species for treating cancer and infection diseases. Moreover, copper cysteamine has a strong luminescence, which can be used for both therapeutics and imaging. In addition, it can also be used for solid state lighting, radiation detection and sensing. However, its electronic structures, and particularly its excited states, are not yet clear. Here, we present a computational study aiming to determine the nature of the excited states involved in the photophysical processes that lead to the luminescence of this compound. This study has been conducted using density functional theory (DFT), using both hybrid functionals and time-dependent DFT. It is found that both absorption and emission involve the replacement of an electron among the 3d and 4s orbitals of one or the other of the two types of Cu atoms found in the system. Our computed results compared well with the experimental absorption and emission results. These results are very helpful for the understanding of the experimental observations.

Photosensitization Versus Photocyclization: Competitive Reactions of Phenylphenalenone in Its Role as Phytoanticipins in Plant Defense Strategies.

Phenalenone (PN) derivatives are involved in plant defense strategies, producing molecular singlet oxygen in a photosensitization process. Many experimental and theoretical studies determined that PN can performe this process with a quantum yield close to 1. However, it has been observed that the efficiency of some of its derivatives is much lower. This is the case of 9-phenylphenalenone (9-PhPN). To elucidate the factors that determine the different photochemistry of PN and its derivate 9-PhPN, we developed a complete active space self-consistent field/multi-configurational second-order perturbation theory study where several deactivation paths through the lowest excited states were explored. We found that the characteristics of the low-lying excited states are similar for both PN and 9-PhPN in the areas near the geometry of excitation. Consequently, the first processes that take place immediately after absorption are possible in both systems, including the population of the triplet state responsible for oxygen sensitization. However, 9-PhPN can also undergo cyclization by a bond formation between the carbonyl oxygen and a carbon atom of the phenyl substituent. This process competes favorably with population of triplet states and is responsible for the decrease of the quantum yield of oxygen sensitization in 9-PhPN relative to PN.

Phenylazopyridine as Switch in Photochemical Reactions. A Detailed Computational Description of the Mechanism of Its Photoisomerization.

Azo compounds are organic photochromic systems that have the possibility of switching between cis and trans isomers under irradiation. The different photochemical properties of these isomers make azo compounds into good light-triggered switches, and their significantly different geometries make them very interesting as components in molecular engines or mechanical switches. For instance, azo ligands are used in coordination complexes to trigger photoresponsive properties. The light-induced trans-to-cis isomerization of phenylazopyridine (PAPy) plays a fundamental role in the room-temperature switchable spin crossover of Ni-porphyrin derivatives. In this work, we present a computational study developed at the SA-CASSCF/CASPT2 level (State Averaged Complete Active Space Self Consistent Field/CAS second order Perturbation Theory) to elucidate the mechanism, up to now unknown, of the cis-trans photoisomerization of 3-PAPy. We have analyzed the possible reaction pathways along its lowest excited states, generated by excitation of one or two electrons from the lone pairs of the N atoms of the azo group (nazoπ*² and nazo²π*² states), from a π delocalized molecular orbital (ππ* state), or from the lone pair of the N atom of the pyridine moiety (npyπ* state). Our results show that the mechanism proceeds mainly along the rotation coordinate in both the nazoπ* and ππ* excited states, although the nazo²π*² state can also be populated temporarily, while the npyπ* does not intervene in the reaction. For rotationally constrained systems, accessible paths to reach the cis minimum along planar geometries have also been located, again on the nazoπ* and ππ* potential energy surfaces, while the nazo²π*² and npyπ* states are not involved in the reaction. The relative energies of the different paths differ from those found for azobenzene in a previous work, so our results predict some differences between the reactivities of both compounds.

Excited-State Decay in the Photoisomerisation of Azobenzene: A New Balance between Mechanisms.

The mechanism of the photoisomerisation of azobenzene has been studied by means of multiconfigurational ab initio calculations. Our results show that it is necessary to account for the dynamic electron correlation in the location of the critical points (CASPT2 optimizations) to obtain a correct description of the topography of the potential energy surfaces of the low energy singlet excited states. By using this methodology, we have found that the state populated by the initial excitation is the S2 (ππ*) state, which decays very efficiently to the S1 (nπ*) state at a pedal-like non-rotated geometry. In the S1 state, relaxation leads to a rotated geometry where the system decays to the ground state, in which further relaxation can lead to either the trans or cis geometries. However, the S1 /S0 conical intersection seam also extends to planar geometries, so this reaction path is also accessible for rotation-constrained systems. Our results explain the experimental observations satisfactorily.

Intramolecular charge transfer in aminobenzonitriles and tetrafluoro counterparts: fluorescence explained by competition between low lying excited states and radiationless deactivation. Part II: influence of substitution on luminescence patterns.

In this paper, we study the mechanisms of charge transfer, luminescence and radiationless decay of three derivatives of 4-aminobenzonitrile (ABN): dimethyl-ABN (DMABN) and the tetrafluorinated derivatives, ABN-4F and DMABN-4F. Our CASSCF/CASPT2 computations explain the different luminescence patterns observed in these three compounds and in comparison with the parent system, ABN, in spite of their similar architecture. We have found that the modifications made by the different substitutions in ABN tune the relative energies of the locally excited (LE) and charge transfer (CT) excited states due to electronic and structural factors. In all cases, the only potentially emitting species of CT character is the twisted-ICT. The increasing stabilization of this later species in the series formed by ABN-4F, DMABN and DMABN-4F explains the increasing intensity of the anomalous emission band in these compounds. Nevertheless, other factors like probability of emission vs. nonradiative decay must have also been taken into account. In fact fluoro-substitution increases the accessibility to conical intersections of the excited states with the ground state, opening an internal conversion channel that decreases the fluorescence quantum yield in the fluorinated derivatives. Our results also show that the involvement of the π-σ* state in the CT process is only possible in ABN-4F, but even in this case it is not probable.

Intramolecular charge transfer in aminobenzonitriles and tetrafluoro counterparts: fluorescence explained by competition between low-lying excited states and radiationless deactivation. Part I: A mechanistic overview of the parent system ABN.

Recent theoretical and experimental studies on the Intramolecular Charge Transfer (ICT) reaction of some members of the aminobezonitrile family (ABN) suggest the involvement of a (π-σ*) excited state (called ICT(CN) in this work) in the ICT process and the existence of a partially twisted ICT species that could be responsible for the anomalous fluorescence observed. These suggestions made us to revise our previous study on the photophysics of ABN and dimethyl-ABN (DMABN), based on the analysis of the potential energy surfaces of the low-lying excited states by means of ab initio calculations, using the CASSCF/CASPT2 protocol. We have first focused our attention to ABN. We have found that the (π-σ*) excited state can be in fact an intermediary state in the path to populate the ICT bright state, although its involvement in the process is not very probable. Our results suggest that the ICT most stable species is the twisted ICT(TICT) and that the partially twisted ICT minimum found in previous studies could be an artefact of the computational method. We have also found that radiationless deactivation is a competitive reaction that must be taken into account to explain the fluorescence patterns of these systems. To confirm our theories, we have also studied other systems with a similar architecture but with a very different luminescence behaviour: dimethyl-ABN, and the 2,3,4,5-tetrafluoro derivatives of ABN and DMABN (ABN-4F and DMABN-4F). The extension of the work and the different approaches in the study of the parent system and of the derivatives make the division of the work in two parts advisable. Part I collects the characterization of the minima and reaction paths connecting the critical points of the potential energy surfaces of the states involved in the ICT reaction of ABN. We have obtained, for the first time, the pathways of radiationless deactivation for this compound. We have also computed transition energies from the excited minima, to interpret the excited state absorption (ESA) spectra obtained experimentally. This information helps in the elucidation of the mechanism of ICT. In Part II we show an analoguous study for DMABN and ABN-4F and DMABN-4F and analyse and compare these results and those of Part I to explain different luminescence behaviours of the four systems studied.

Intramolecular Charge Transfer in 4-Aminobenzonitrile Does Not Need the Twist and May Not Need the Bend.

A study combining accurate quantum chemistry and full-dimensional quantum dynamics is presented to confirm the existence of an ultrafast radiationless decay channel from the charge-transfer state to the locally excited state in 4-aminobenzonitrile. This intramolecular charge-transfer pathway proceeds through a newly found planar conical intersection, and it is shown to be more efficient in the presence of acetonitrile than in the gas phase. Our results are consistent with recent experimental observations.

Insight into the mechanisms of luminescence of aminobenzonitrile and dimethylaminobenzonitrile in polar solvents. An ab initio study.

The dual fluorescence of 4-(dimethylamino)benzonitrile (DMABN) has been intensively studied in the last decades, but surprisingly there is not any detailed theoretical study of its photochemistry in polar solvents. In this work, we rationalize the different luminescent behavior of 4-aminobenzonitrile (ABN) and DMABN in acetonitrile by a computational study developed at the CASSCF/CASPT2 level and using the Polarized Continuum Model to reproduce the solvent environment. We present here the critical geometries, energies, and connections between the potential energy surfaces of the low-lying excited states: the locally excited state (LE) and several intramolecular charge transfer states (ICT). The computational results show that the topology of the potential energy surfaces (PES) does not change substantially when the effect of a polar solvent is included, in comparison with the gas phase. For DMABN, though, polar solvents stabilize preferentially the ICT states in such a way that the different interplay with the LE state induces strong qualitative changes in the photochemistry of this compound. Specifically, the planar ICT (PICT) species located on the S2 surface in the gas phase is, in acetonitrile, located on the S1 surface. that is, at the geometry of the PICT minimum, the LE state is higher in energy than the ICT. Now LE and PICT minima are practically degenerate and, given that both correspond to first excited state species, emission can take place from both of them. However, the twisted ICT (TICT) species is still the most favored thermodynamically so it is expected that this species would be preferentially populated. On the other hand, for ABN the equilibrium lies in favor of LE, as the TICT species was found at a much higher energy with a low reaction barrier toward LE. This explains why dual fluorescence cannot be observed in ABN, even in polar solvents.

Computational study of the mechanism of the photochemical and thermal ring-opening/closure reactions and solvent dependence in spirooxazines.

The spirooxazine/merocyanine couple constitutes a photochromic system that can change from the colorless spirooxazine to the intensely colored merocyanine by thermal or photochemical activation by a reaction that opens the spiro ring of the oxazine. The mechanisms of the ring-opening/closure reactions that interconnect these two isomers have been elucidated by means of a computational study. First, we have used the CASSCF/CASPT2 method to determine in detail these mechanisms in the gas phase for a small model that contains the photoactive part of the whole system. We have found that the state of spirooxazine excited by the initial absorption changes gradually to a lower excited state that is involved in a conical intersection that connects it with the ground state of merocyanine. The same conical intersection is involved in the backward photochemical reaction. Second, using a larger model that includes all the heteroatoms of the system and using the DFT (B3LYP) method, we have studied the influence of a solvent environment on the thermal equilibrium between the open and the closed species. It has been observed experimentally that the thermal equilibrium between these forms is practically unaltered by polar aprotic solvents, but it can be displaced toward the colored form in mixtures of polar protic and aprotic solvents, even if the first one is found in a very small proportion. To reproduce the experimental environments, we have taken into account the long-range effect of the polar aprotic solvent considering it a polarizable continuum, as done in the PCM method, and the short-range effect of the protic solvent including some explicit water molecules in the cluster studied at the atomic level. The results obtained are in good agreement with the experimental observations and explain the reason for this peculiar behavior.

Mechanism of the photochemical process of singlet oxygen production by phenalenone.

Phenalenone (PN) is a very efficient singlet oxygen sensitiser in a wide range of solvents. This work uses ab initio quantum chemical calculations (CASSCF/CASPT2 protocol) to study the mechanism for populating the triplet state of PN responsible for this reaction, the (3)(π-π*) state. To describe in detail this reaction path, the singlet and triplet low-lying excited states of PN have been studied, the critical points of the potential energy surfaces corresponding to these states located and the vertical and adiabatic energies calculated. Our results show that, after the initial population of the S(2) excited state of (π-π*) character, the system undergoes an internal conversion to the (1)(n-π*) state. After populating the dark S(1) state, the system relaxes to the (1)(n-π*) minimum, but rapidly populates the triplet manifold through a very efficient intersystem crossing to the (3)(π-π*) state. Although the population of the minimum of this triplet state is strongly favoured, a conical intersection with the (3)(n-π*) surface opens an internal conversion channel to this state, a path accessible only at high temperatures. Radiationless deactivation processes are ruled out on the basis of the high-energy barriers found for the crossings between the excited states and the ground state. Our computational results satisfactorily explain the experimental findings and are in very good agreement with the experimental data available. In the case of the frequency of fluorescence, this is the first time that these data have been theoretically predicted in good agreement with the experimental results.

The peculiar spectral properties of amino-substituted uracils: a combined theoretical and experimental study.

A detailed experimental and computational study of the absorption and fluorescence spectra of 5-aminouracil (5 AU) and 6-aminouracil (6 AU) in aqueous solution is reported. The lowest energy band of the steady-state absorption spectra of 5 AU is considerably red-shifted, noticeably less intense, and broader than its counterpart in uracil (U). On the contrary, the 6 AU lowest energy absorption peak is close in energy to that of U, but it is much narrower and the transition is much more intense. The emission properties of 5 AU, 6 AU, and U are also very different. Both amino-substituted compounds exhibit indeed a much larger Stokes shift as compared to U, and the emission band of 5 AU is much narrower than that of 6 AU. Those features are fully rationalized with the help of PCM/TD-PBE0 calculations in aqueous solution and MS-CASPT2/CASSCF calculations in the gas phase. A stable minimum on the potential energy surface of the lowest energy bright state is found for 5 AU, both in the gas phase and in aqueous solution. For 6 AU a barrierless path leads to the conical intersection with the ground electronic state, but a nonplanar plateau region is predicted in aqueous solution, which is responsible for the very large Stokes shift. Some general considerations on the excited-state dynamics of uracil derivatives are also reported.

On the mechanism of the photoinduced magnetism in copper octacyanomolybdates.

Ab initio calculations show that a possible mechanism for the photomagnetism in copper octacyanomolybdate compounds consists of the initial excitation of the diamagnetic Cu(II)-Mo(IV-CS) pair to a Cu(II)-Mo(IV-T) state, whose geometry relaxation stabilizes the magnetic doublet and quartet states.

The decay from the dark npi* excited state in uracil: an integrated CASPT2/CASSCF and PCM/TD-DFT study in the gas phase and in water.

By integrating the results of MS-CASPT2/CASSCF and TD-PBE0 calculations, we propose a mechanism for the decay of the excited dark state in pyrimidine, fully consistent with all the available experimental results. An effective conical intersection (CI-npi) exists between the spectroscopic pi/pi* excited state (Spi) and a dark n/pi* state (Sn), and a fraction of the population decays to the minimum of Sn (Sn-min). The conical intersection between Sn and the ground-state is not involved in the decay mechanism, because of its high energy gap with respect to Sn-min. On the other hand, especially in hydrogen bonding solvents, the energy gap between Sn-min and CI-npi is rather small. After thermalization in Sn-min, the system can thus recross CI-npi and then quickly proceed on the Spi barrierless path toward the conical intersection with the ground state.

Unified model for the ultrafast decay of pyrimidine nucleobases.

Ultrafast decay processes detected after absorption of UV radiation in gas-phase pyrimidine nucleobases uracil, thymine, and cytosine are ascribed to the barrierless character of the pathway along the low-lying 1(pipi*) hypersurface connecting the Franck-Condon region with an out-of-plane distorted ethene-like conical intersection with the ground state. Longer lifetime decays and low quantum yield emission are on the other hand related to the presence of a 1(pipi*) state planar minimum on the S1 surface and the barriers to access other conical intersections. A unified model for the three systems is established on the basis of accurate multiconfigurational CASPT2 calculations, whereas the effect of the different levels of theory on the results is carefully analyzed.

Theoretical investigation of luminescence behavior as a function of alkyl chain size in 4-aminobenzonitrile alicyclic derivatives.

There is still controversy about the structure of the intramolecular charge transfer (ICT) emitting species in pi-electron donor-acceptor systems that show dual fluorescence. Although the twisted ICT model is quite generally accepted, the planar ICT model is not ruled out because firm experimental evidence supports it. Among these it is the fact that some rigidized systems such as bicyclic 4-aminobenzonitrile derivatives exhibit dual fluorescence. We present here an ab initio CASSCF/CASPT2 study of a series of these compounds with the alicyclic chain ranging from 5 to 7 carbon atoms and compare their ICT mechanism with the more flexible 4-aminobenzonitrile (ABN) and 4-(dimethylamino)benzonitrile (DMABN). We present the energetics, geometries, and valence bond structures of the critical points of the potential-energy surfaces of the ground, local excited (LE), and ICT states. Our results show that the photophysical differences of the studied systems may be rationalized by two factors: the position of the ICT and LE potential-energy surfaces at the first stages of the ICT reaction and the relative energies of the excited-state minima. Computational evidence is presented that a twisted ICT structure can be adopted in some molecules such as NXC6 and NXC7 and that the anomalous band of the fluorescence spectra of these systems is emitted from a twisted ICT species.

Efficient photochemical merocyanine-to-spiropyran ring closure mechanism through an extended conical intersection seam. A model CASSCF/CASPT2 study.

A mechanism of the thermal and photochemical bleaching of merocyanine to spiropyran is proposed on the basis of CASSCF/CASPT2 calculations on the 6-(2-propenyliden)cyclohexadienone model system. Our results suggest that this photochemical transformation takes place in two steps. First, the initially pumped 1(pi-pi) S2 undergoes radiationless decay to 1(n-pi) S1 via an extended S2/S1 conical intersection seam that runs approximately parallel to the trans-to-cis isomerization coordinate, a few kilocalories per mole higher in energy. Thus, S2 --> S1 internal conversion is possible at all values of the S2 trans-to-cis reaction coordinate. Second, on the S1 potential energy surface, there is a barrierless ring closure reaction path from the S1 cis minimum that leads to a peaked S1/S0 conical intersection where the deactivation to the ground state takes place. The inertia of the moving nuclei then drives the system toward the ground-state minimum of the 2H-chromene product. Thus, the extended seam topology of the S2/S1 conical intersection and the coordinate of the branching space of the S1/S0 conical intersection are essential to explain the efficiency and high speed of this reaction.

Intramolecular charge transfer in 4-aminobenzonitriles does not necessarily need the twist.

In electron donor/acceptor species such as 4-(dimethylamino)benzonitrile (DMABN), the excitation to the S(2) state is followed by internal conversion to the locally excited (LE) state. Dual fluorescence then becomes possible from both the LE and the twisted intramolecular charge-transfer (TICT) states. A detailed mechanism for the ICT of DMABN and 4-aminobenzonitrile (ABN) is presented in this work. The two emitting S(1) species are adiabatically linked along the amino torsion reaction coordinate. However, the S(2)/S(1) CT-LE radiationless decay occurs via an extended conical intersection "seam" that runs almost parallel to this torsional coordinate. At the lowest energy point on this conical intersection seam, the amino group is untwisted; however, the seam is accessible for a large range of torsional angles. Thus, the S(1) LE-TICT equilibration and dual fluorescence will be controlled by (a) the S(1) torsional reaction path and (b) the position along the amino group twist coordinate where the S(2)/S(1) CT-LE radiationless decay occurs. For DMABN, population of LE and TICT can occur because the two species have similar stabilities. However, in ABN, the equilibrium lies in favor of LE, as a TICT state was found at much higher energy with a low reaction barrier toward LE. This explains why dual fluorescence cannot be observed in ABN. The S(1)-->S(0) deactivation channel accessible from the LE state was also studied.

Thermal and photochemical rearrangement of bicyclo[3.1.0]hex-3-en-2-one to the ketonic tautomer of phenol. Computational evidence for the formation of a diradical rather than a zwitterionic intermediate.

The ground state (S(0)) and lowest-energy triplet state (T(1)) potential energy surfaces (PESs) concerning the thermal and photochemical rearrangement of bicyclo[3.1.0]hex-3-en-2-one (8) to the ketonic tautomer of phenol (11) have been extensively explored using ab initio CASSCF and CASPT2 calculations with several basis sets. State T(1) is predicted to be a triplet pipi lying 66.5 kcal/mol above the energy of the S(0) state. On the S(0) PES, the rearrangement of 8 to 11 is predicted to occur via a two-step mechanism where the internal cyclopropane C-C bond is broken first through a high energy transition structure (TS1-S(0)()), leading to a singlet intermediate (10-S(0)()) lying 25.0 kcal/mol above the ground state of 8. Subsequently, this intermediate undergoes a 1,2-hydrogen shift to yield 11 by surmounting an energy barrier of only 2.7 kcal/mol at 0 K. The rate-determining step of the global rearrangement is the opening of the three-membered ring in 8, which involves an energy barrier of 41.2 kcal/mol at 0 K. This high energy barrier is consistent with the fact that the thermal rearrangement of umbellulone to thymol is carried out by heating at 280 degrees C. Regarding the photochemical rearangement, our results suggest that the most efficient route from the T(1) state of 8 to ground state 11 is the essentially barrierless cleavage of the internal cyclopropane C-C bond followed by radiationless decay to the S(0) state PES via intersystem crossing (ISC) at a crossing point (S(0)()/T(1)()-1) located at almost the same geometry as TS1-S(0)(), leading to the formation of 10-S(0)() and the subsequent low-barrier 1,2-hydrogen shift. The computed small spin-orbit coupling between the T(1) and S(0) PESs at S(0)()/T(1)()-1 (1.2 cm(-)(1)) suggests that the ISC between these PESs is the rate-determining step of the photochemical rearrangement 8 --> 11. Finally, computational evidence indicates that singlet intermediate 10-S(0)() should not be drawn as a zwitterion, but rather as a diradical having a polarized C=O bond.