Useful Spectrokinetic Methods for the Investigation of Photochromic and Thermo-Photochromic Spiropyrans

Maafi, Mounir

doi:10.3390/molecules13092260

Open AccessReview

Useful Spectrokinetic Methods for the Investigation of Photochromic and Thermo-Photochromic Spiropyrans

by

Mounir Maafi

Leicester School of Pharmacy, De Montfort University, The Gateway, Leicester, LE1 9BH, UK

Molecules 2008, 13(9), 2260-2302; https://doi.org/10.3390/molecules13092260

Submission received: 14 August 2008 / Revised: 19 September 2008 / Accepted: 19 September 2008 / Published: 25 September 2008

(This article belongs to the Special Issue Spiro Compounds)

Download

Browse Figures

Versions Notes

Abstract

:

This review reports on the main results of a set of kinetic elucidation methods developed by our team over the last few years. Formalisms, procedures and examples to solve all possible AB photochromic and thermophotochromic kinetics are presented. Also, discussions of the operating conditions, the continuous irradiation experiment, the spectrokinetic methods testing with numerical integration methods, and the identifiability/distinguishability problems, are included.

Keywords:

Spectrokinetic methods; kinetic elucidation; spiropyrans; photochromic and thermochromic compounds

1. Introduction

Chromism often describes a reversible colour change of a material that can be induced photochemically (photochromism) or thermally (thermochromism) [1,2,3,4]. It translates the transformation of chemical species between two (or more) forms which have different absorption spectra (Figure 1). In the case where only two species can be monitored during the photo- or thermochromism, the reaction is named an AB system (vide infra section 2, section 4 and section 5). When both the latter processes responsible for the colour change take place, the AB system is called a thermo-photochromic system.

The initial species A is usually the thermodynamically stable and colourless (or weakly coloured) form whereas the product species B is the coloured form. Although these processes can be observed both in solution and solid state, this review considers only the investigations carried out in solutions.

Figure 1. Typical absorption spectra for AB photochromic or thermochromic systems.

Organic photochromic, thermochromic and thermo-photochromic compounds of the AB type, like spiropyrans, have known an increasing interest over the last few decades because of their effective and/or potential applications in a number of technologically important fields [1,2,3,4,5,6,7]. Such a position boosted not only the synthetic methods, that has allowed tens of new derivatives to be prepared every year, but also reinforced the need for the development of reliable and easy-to-implement kinetic methods that provide a better insight in the thermophotochromic and photochromic kinetic properties of such compounds [1,2,3,4,5,6,7].

Since the 1960s, continuous research efforts have been conceded for the elaboration of elucidation methods. Dozens of approaches have been proposed thus far for the various AB reaction cases [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44]. These methods met with varying success. Most, if not all, were originally constructed on simplifying hypotheses, did not solve the identifiability/distinguishability problems (vide infra § 8) and/or required extensive experimental or calculation means. Comprehensive analyses and discussions of these methods and their limitations can be found in several publications [2,3,7,18,45,46,47].

From this perspective, we can identify at least two important interests in developing new and reliable kinetic elucidation methods for AB systems (which should present a better performance and fewer limitations than those already proposed). On one hand, this concern addresses a real need for such efficient methods that do not necessitate coercive experimental conditions and/or extensive data treatments while capable of generating reliable and accurate kinetic and spectroscopic data, i.e., these methods deliver the true kinetic solution. On the other hand, such reliable kinetic results can be put to advantage in the design of new materials endowed with the specific features that are required for particular technological and/or scientific applications.

In order to build up new reliable kinetic investigation tools and useful kinetic elucidation methods, it is imperative that a few general but important requirements are met. For example, the approaches must be based on tight mathematical formalisms, use only observables and avoid as much as possible to rely on assumptions and simplifying hypotheses. They also must clearly identify whether pure kinetic data are sufficient to achieve a complete elucidation (otherwise, the type and the number of supplementary information required for a total elucidation of the kinetics at hand, must be specified). Furthermore, it is necessary that such methods address, if not solve, the identifiability/ distinguishability problems due to the direct incidence of these problems on the reliability of the method itself. Finally, they must be tested against data obtained by an independent technique, such as numerical integration methods (§ 7), before considering their application to experimental data.

Such a strategy has been adopted for the development of the elucidation methods that are reviewed in this paper, which encompass all possible reaction schemes that can govern organic photochromic and thermphotochromic AB systems.

2. Physical properties, reactivity and kinetic study of spiropyrans

Background

The chemical structure of spiropyrans involves a benzopyran and an indole moiety joined at a spiro centre (C_spiro). The two possible enantiomers are separated by an energy barrier of approximately 86 kJ mol^-1 [48]. Depending on the molecular structure, the reactivity of spiropyran derivatives is believed to be governed by reversible thermal and/or photochemical mechanisms (Scheme 1). Hence, either or both A and B species can be thermally and photochemically reactive and the reaction sequence might include up to two photochemical and two thermal reaction steps.

Scheme 1. Closed and open forms of spiropyrans.

The ring opening reaction of spiropyrans (SP) yields a planar merocyanine form (MC) which is characterized by an extended conjugated system (Scheme 1). Both quiniodal and zwitterionic forms can be postulated for MC (Scheme 2). The stability of these resonance forms is sensitive to solvent polarity, with the former likely to be more stable in non-polar media.

The increased conjugation of both resonance forms shifts the electronic absorption of the trans conformer, MC, towards the visible range of the electromagnetic spectrum (Figure 1). Therefore, MC isomers are deeply coloured compared to colourless or light yellow closed-form isomers (SP). This chromism is also accompanied by a significant increase of the ground state dipole moment of the species from 4D for SP to 18D for MC [49].

Both thermal and photochemical transformations of SP into MC are responsible of the breakage of the C_spiro–O bond. The activation energy of the thermal conversion ranges between 80 and 130 kJ mol^-1 [46].

Scheme 2. Quiniodal and zwitterionic forms of MC.

Time resolved studies of the forward SP → MC photoreaction, have demonstrated that the bond cleavage originated from the first excited singlet state of SP (¹SP*) within less than 200 fs of the excitation pulse [50]. The reaction stages subsequent to the initial bond cleavage have been extensively studied. They ultimately lead to eight possible “cis” and “trans” isomers of MC. The four cis conformations are much less stable than their trans counterparts (Scheme 3). Also, according to computational calculations, TTC is the most stable isomer for SP molecules bearing NO₂ groups [51,52]. Nonetheless, the overall thermal isomerization pathway, the major resonance forms and/or the most stable isomer remain subject to debate.

Scheme 3. Some intermediates of SP and MC interconversions.

From a photophysical viewpoint, the fate of the excited state, ¹SP^*, produced after light excitation, is dominated by either a singlet or a mixed singlet-triplet pathways depending on whether or not the SP molecule possesses a nitro substituent. The successive transformations leading to the ground state MC (¹SP^* → ¹MC₀), which include cis and trans isomer intermediates (X), have various time scales ranging between tens of picoseconds to a few milliseconds [2,3,50,51,52] (Scheme 4).

The photobleaching process leading to ring closure (i.e. the decay of excited MC, ¹MC^* → SP) has also been investigated by time resolved techniques [53,54]. The phototransformation of MC is believed to proceed via similar types of decay routes (i.e. which can involve a mixed singlet-triplet and/or a singlet pathway depending on whether a nitro substituent is present in the MC molecule). Here as well, the longest lifetime recorded for the intermediates does not exceed a few ms.

The fast photochemical processes observed for both photocolouring and photobleaching of these photochromes has been interpreted by the occurrence of a conical intersection between the ground and excited state potential energy surfaces [55,56].

The thermal recyclisation (MC → SP), on the other hand, obeys first order kinetics in solution. Lifetimes of thermocolouring and thermobleaching reactions were found to be in the sec/min ranges.

Scheme 4. A simplified representation of SP → MC phototransformation.

Considerations for continuous irradiation experiments

The progress of spiropyran photreactions can also be monitored while the reaction medium is continuously subjected to irradiation. In continuous irradiation experiments, conversely to time resolved studies, the excited state is gradually populated and only long living species can be observed.

The results obtained from time resolved studies clearly showed that amongst the species manifold involved in the photoprocesses of spiropyrans and their merocyanine isomers, only SP and MC (Scheme 1) have lifetimes of over a second. Accordingly, only these two species (and not any of the other intermediates) are expected to be observed and monitored during continuous irradiation experiments (the time resolution of modern diode array spectrophotometers is around 0.5 s). Hence, and even though the real photochemical and thermal mechanisms of such materials are complex, the observed kinetics for these materials can be considered to be similar to that involving bimolecular AB systems. It is however worth noting that in such a situation, the obtained photochemical efficiencies will represent overall rather than unimolecular quantum yields.

There are also a few considerations to be pointed out regarding thermal reactions. Indeed, since both forward and reverse thermal reactions (SP→MC and MC→SP) proceed in sec/min time spans, it is unrealistic to expect separating either of the thermally reactive species (MC and SP) by a physical method. As a matter of fact, the electronic absorption spectra of such transient species are not experimentally accessible and only their kinetic data can be used for elucidation purposes. Therefore, the determination of both species spectra must be considered, for the cases where thermal reactions are involved. This is one of the targets that ought to be met by the kinetic elucidation method (vide infra Table 2 and Table 3).

Finally, it should be kept in mind that the irradiation of the initial species with UV light may also initiate a photoreaction of the formed photoproduct (since SP and MC usually share absorption bands in the UV-range). The latter photoreactions’ efficiencies might be irradiation dependent as well. Both these issues cannot be ignored in the absence of tangible information.

3. AB systems with similar kinetic behaviour to spiropyrans

The mechanistic considerations, the kinetic behaviour and the spectroscopic features reviewed above with relation to spiropyran kinetics are also common to a number of families of photochromic and thermo-photochromic materials (Figure 2).

Figure 2. A selection of photochromic and/or thermochromic AB systems possessing a kinetic behaviour similar to that of spiropyrans.

Hence, the kinetic elucidation methods reviewed below (§ 9) are equally applicable to these materials when their solutions are subjected to continuous irradiation. It is also worth noting that most of these kinetic elucidation methods are useful to investigate non-chromic AB systems as well (Figure 3).

Figure 3. Some non-chromic AB systems.

4. Typical kinetic mechanisms

As stipulated earlier in Section 2, the mechanism governing each photochromic and/or thermophotochromic AB system, can encompass up to four reaction steps, of which two are pure photochemical and two are pure thermal processes. Therefore, a total of eleven sequences can be postulated for these reactions (Table 1). Most of these reaction schemes have effectively been postulated for either spiropyrans or other AB systems of similar kinetic behaviour (Figure 2). However, a few less frequent schemes (S₄, S₈ and S₁₀), are included here to complete the series of kinetic elucidation methods for all possible cases of AB reactions.

Photo and thermochromic interconversions might lead to undesirable by products and ultimately to depletion of the active chromic species. Resistance to such degradation (also called fatigue resistance) varies between photochromes of the same series and between families of photochromes. From a mechanistic point of view, this means that at least one reaction step is added to the sequences given in Table 1 (e.g. B→ C with C being the degradation product of species B). Mechanisms including fatigue are not treated in the present review.

Table 1. Possible reaction schemes for AB systems.

**Table 1.** Possible reaction schemes for AB systems.
Reaction	Label	References to Spiropyran compounds
A unimolecular photoreaction; AB(1ϕ)	S₁	61,76,77
A unimolecular thermal reaction; AB(1k)	S₂	61
A photochemical reaction reversed by a thermal reaction; AB(1ϕ,1k)	S₃	77, 78
A thermal reaction reversed by a photochemical reaction; AB(1k,1ϕ)	S₄	61,78,79,80
Pure thermal opposed reactions; AB(2k)	S₅	41,81,82,83
Pure photochemical opposed reactions; AB(2ϕ)	S₆	79,84
Photochemical opposed reactions coupled with a reverse thermal reaction. AB(2ϕ,1k)	S₇	32,80, 85,86
Photochemical opposed reactions coupled with a forward thermal reaction. AB(2ϕ,1k)	S₈	32,80,85,86
Thermal opposed reactions coupled with a forward photochemical reaction. AB(1ϕ,2k)	S₉	82,87
Thermal opposed reactions coupled with a reverse photochemical reaction. AB(1ϕ,2k)	S₁₀	–
Thermal and photochemical opposed reactions. AB(2k,2ϕ)	S₁₁	87,88,89

5. Observables and Experimental Conditions

Photochemical Reactors

The photochemical reactor is a quartz cell where the sample solution is continuously and uniformly stirred and maintained at a constant temperature. The dimensions of the cell may vary (Scheme 4) but usually spectrophotometric cuvettes are used (with a total volume of 1 × 1 × 4 cm³). Ususally, the reactor serves also as a sample holder for spectrophotometric measurements during reaction progress.

Irradiation Conditions

Kinetic investigations using continuous irradiation require that the irradiation (or excitation) beam is monochromatic, of low photon flux and delivered by a highly powered lamp (usually, lasers are not used for this type of experiments). The beam is directed on the sample at an angle to the probing light of the spectrophotometer to avoid interference at the spectrophotometer (or diode array) detector. For this type of experiments, the directions of the probing and irradiation beams are most often perpendicular, i.e. the irradiation is carried out either laterally or top-down ((a) and (b) respectively in Scheme 4). Polychromatic beams are not suitable for continuous irradiation experiments because both the incident light intensity and the photochemical quantum yield values (I₀ and ϕ, respectively) may vary with wavelength.

Scheme 4. The photoreactor and the possible paths of irradiation and probing lights.

The optical path length of both the probing light (l_probe) and lateral irradiation beam (l_irr) correspond to the dimension of the photoreactor (1 cm for the classical spectrophotometric cuvettes). However, for top irradiations (Scheme 4b), l_irr will be equal to the length of the liquid solution inside the photochemical reactor that is subjected to the irradiation (e.g. 0 cm < l_irr < 4 cm, for a cuvette). In all cases, it is important to identify the (experimental) values of l_probe and l_irr as they are included in the formalisms of the kinetic elucidation methods (see next section).

The sample can be subjected to two types of irradiations, namely isosbestic and non-isosbetic irradiations. The former corresponds to the situation where the wavelength of the excitation beam (λ_irr) is identical to that of an isosbestic point (or isosbestic region, see next section) whereas the latter means that λ_irr is different from that of an isosbestic point. The elucidation methods developed hereafter may require either or both types of irradiations.

Observables

Not all fundamental kinetic features and species’ attributes can be measured. The experimentally accessible reaction attributes are called observables. The remaining features are the unknowns. In a typical continuous irradiation experiment, the initial concentration (C₀), the intensity of the irradiation beam (I₀), the medium temperature, the optical path lengths l_probe and l_irr are all measurable quantities which can be defined before the start of the experiment. The accessible kinetic data collected on a given photochromic and/or thermochromic reactive system are represented by plots of the measurable variation of the absorbance (M) with reaction time (t). It is important to underline that spectrophotometric techniques of analysis allow only one absorbance (labelled M), monitored by the probing light at the condition of observation (i.e. where the optical path length is l_probe), to be recorded. The absorbance of the medium in the condition of irradiation (along the optical path length l_irr) is usually not directly accessible (Scheme 4), and hence it is not considered to be an observable.

A particular recorded M(t) curve is specific to an irradiation (λ_irr) and an observation (λ_obs) wavelength as well as to l_probe, l_irr, I₀ and the medium temperature (T). Such curves are called kinetic traces. They are obtained either during the irradiation or the thermal relaxation of the reactive medium. The photochemical traces relate to the progress of both photochemical and thermal reactions of a given sequence while thermal curves are exclusively due to thermal processes (these traces are recorded in the dark, i.e. I₀= 0). Also, it is worth noting that traces for individual species are most often accessible if the reaction medium is observed at a wavelength (λ_obs) where only one species absorbs (most likely in the visible region of the spectrum because the spectra of reactant and product often overlap throughout the UV region).

A number of other observables may be obtained by alternative (spectroscopic and/or physico-chemical) methods. However, the methods presented hereafter have been built on the hypothesis that kinetic data and the observables listed above are the only output of the experiment. This is because alternative techniques and physical methods are time consuming, expensive and more importantly they are difficult to set up for a majority of sequences (Table 1) due to the transient character of one or both species involved in the reaction. Reference to the use of alternative techniques will be made where appropriate in the following sections.

As alluded to above, in order to perform isosbestic irradiations it is necessary to observe an isosbestic point on the absorption spectra of the reactive medium during its evolution with irradiation time. Thus far, the presence of an isosbestic point was confirmed for all known photochromes and thermo-photochromes. Incidentally, this feature is often used to indicate the smooth progress of the reaction and the absence of by-products. Alternatively, similar irradiations can be carried out using wavelengths situated in isosbestic regions where the spectra of A and B are the same (where the species’ molar extinction coefficients are equal, ε_A = ε_B, over a spectral range).

According to the properties of this type of AB reactions, the photostationary state (PSS) and/or the state of thermal equilibrium (STE) are reached within minutes. They correspond to plateau regions on the kinetic traces recorded at long reaction times. Some of the kinetic elucidation methods reported here require data relative to PSS and STE, and hence these states are considered to be observables.

The unknowns

The fundamental features that are not directly accessible from the continuous irradiation experiment (the unknowns) are of two types, kinetic and spectroscopic quantities. The unknown kinetic reaction attributes are the quantum yields of the photochemical steps (ϕ_AB or ϕ_BA) and the rate constants of the thermal reactions (k_AB and k_BA). The inaccessible spectroscopic data are the molar extinction coefficients of species A and B (ε_A and ε_B respectively). Quantum yields and molar extinction coefficients are wavelength dependent while rate constants are temperature dependent. Of course, the total number of unknowns varies with the sequence as indicated in Table 2. In this respect, it is relevant to consider that ε_A is unknown for all thermally reactive initial species (A) because the progress of such reactions is immediate following the preparation of the stock solutions. (Reliable spectrophotometric measurements cannot be carried out promptly enough due to relatively fast rate constants, k_AB [2,3]). In addition to the above unknowns, the species equilibrium concentrations (at the PSS, C_i(PSS), and at the STE, C_i(STE); with i = A or B) are not measurable in the type of experimental conditions previously set and therefore they must be accurately defined by the methods as well.

Table 2. Sets of unknowns specific to each kinetic case of Table 1.

**Table 2.** Sets of unknowns specific to each kinetic case of Table 1.
Reaction scheme	Possible unknowns for the sequence
S₁	ϕ_AB and ε_B
S₂	k_AB, ε_A and ε_B
S₃	ϕ_AB, k_BA and ε_B
S₄	k_AB, ϕ_BA, ε_A and ε_B
S₅	k_AB, k_BA, ε_A and ε_B
S₆	ϕ_AB, ϕ_BA and ε_B
S₇	ϕ_AB, ϕ_BA, k_BA and ε_B
S₈	ϕ_AB, ϕ_BA, k_AB and ε_B
S₉	k_AB, k_BA, ϕ_AB, ε_A and ε_B
S₁₀	k_AB, k_BA, ϕ_BA, ε_A and ε_B
S₁₁	k_AB, k_BA, ϕ_AB, ϕ_BA, ε_A and ε_B

The kinetic solution sought by each elucidation method is represented by the true set of absolute values determined by the method for the unknowns corresponding to the kinetic case considered (see Section 8 and Section 9)

6. Fundamental Kinetic Laws

The mathematical formalism required for the build up of the elucidation methods is considered in this section. The labelling employed here (such as a₄, a₅, a₁₉, M, m₀….) is similar to that used in the original papers [90,91,92,93,94,95]. A glossary of the symbols and labels adopted in this review is included in Section 13.

In the following, irr, isos and obs stand respectively, for the irradiation, isosbestic and observation wavelengths (where irr ≠ isos, and, obs may or may not be equal to irr). The absorbance (M) and the initial reaction velocity (m₀) are referred to an excitation and an observation wavelength (λ_exc/λ_obs), i.e. either irr or isos for the excitation of the photochemical reactions, and either obs, irr or an arbitrary λ, for the monitoring wavelength. The photokinetic factor (F) which is always defined by the irradiation features ( Molecules 13 02260 i073

,

, l_irr...etc) can be time dependent F^irr(t) or constant when measured at t = 0, F^irr(t) = Molecules 13 02260 i075

, or corresponding to either the PSS, F^irr(PSS) = Molecules 13 02260 i076

, or the STE, F^irr(STE) = Molecules 13 02260 i077

. The absorbance Molecules 13 02260 i078

monitored in the observation conditions (l_probe), is labelled Molecules 13 02260 i079

throughout.

The Fundamental Kinetic Equations for AB Systems

The continuous irradiation experiments must be carried out in a photoreactor where a homogeneous sample is continuously stirred and irradiated at a constant temperature. The fundamental differential equation describing the variation of species concentrations with reaction time (Eq.1) is set for the case where the wavelength of the monochromatic irradiation beam (λ_irr) is supposed to correspond to a spectral region where species A and B absorb different amounts of the incident light, and the concentration of the excited state remains negligible. C_A(t) and C_B(t) are the concentrations of species A and B, respectively. Eq.1 corresponds to the most inclusive reaction mechanism S₁₁ (in Table 1). Hence, individual differential equations for the remaining S₁-S₁₀ sub-sequences can easily be inferred from Eq. 1.

(1)

The F^irr is the variable photokinetic factor, which is expressed as a function of the total absorbance at the irradiation wavelength (l_irr),

(2)

and which formula is

(3a)

Constant values of the photokinetic factor can be calculated for various situations including for the initial absorbance recorded at the start of a photochemical reaction ( Molecules 13 02260 i075

), the absorbance reached at PSS ( Molecules 13 02260 i076

) or STE ( Molecules 13 02260 i077

), and the isosbestic absorbance (F^isos). These photokinetic factors are given by the following equation where only observables are used (with, in Eq. 3b, init = 0, STE or PSS),

(3b, c)

The fundamental equation (Eq. 1) can be rearranged to express the rate of change of the total observed absorbance

(4)

where the sample is irradiated at λ_irr and observed at λ_obs (with λ_irr ≠ λ_obs), as

(5)

From Eq. 5 we can derive the general expression of the initial velocity ( Molecules 13 02260 i080

) of the reaction, as

(6)

where, depending on the reaction considered, the argument “init” can be equal to 0, PSS or STE.

The mass balance equation gives the relationship between the species concentrations at time t, the initial time (t = 0), STE and PSS and the total concentration C₀, as

C₀ = C_A(t) + C_B(t) = C_A(0) + C_B(0) = C_A(STE) + C_B(STE) = C_A(PSS) + C_B(PSS)

(7)

It is also possible to extract the individual analytical expressions of the species concentrations at PSS and STE (at these states, the absolute value of the concentration is considered invariant with reaction time, hence, dC_A/dt = 0).

(8.a)

(8.b)

and hence C_B(PSS) and C_B(STE) can be worked out from Eqs. 7 and 8. By using the latter equations in the equation of Molecules 13 02260 i081

, Eq. 4, we can write Molecules 13 02260 i082

formula, as

(9)

From a mathematical viewpoint, Eqs. 1 and 5 cannot be solved through a closed-form integration. The difficulty is due to the time dependence of the photokinetic factor (F as given in Eq. 3a). The only known exception is that of the unimolecular photoreaction AB(1ϕ) whose differential equation can be analytically integrated [93]. The corresponding kinetic model is reviewed at the beginning of Section 9.

However, the differential equations (Eqs. 1 and 5) are readily solved for thermal reactions and also for photoreactions if F remains constant. This condition is met whenever the wavelength of the monochromatic irradiation beam is equal to that of an isosbestic point (§ 6). The general equations relative to these cases are presented in the next section.

Basic Kinetic Laws for AB Kinetics When F is Constant

The kinetic laws, representing the temporal variation of the species concentrations (C_A and C_B), for all AB(nk,nϕ) systems studied under isosbestic irradiations (with Molecules 13 02260 i083

its intensity), whether they involve pure photochemical, pure thermal or a combination of photochemical and thermal reaction steps, have been defined [90] as:

(10a)

(10b)

where init = 0, STE or PSS, as previously set. The a_i coefficients are functions of the kinetic parameters [90]; they are defined for the photochemical reaction (where both the thermal and the photochemical processes operate, and indexed by the letter P, a_iP) and the thermal reaction (indexed T, a_iT, where only the thermal processes are active, i.e. after the irradiation is cut off and Molecules 13 02260 i083

= 0). Hence, for an AB(2k,2ϕ) system as given by Sequence S₁₁ (Table 2), we obtain :

(11a,b)

a_4T = k_AB a_5T = −k_AB

(12a,b)

Coefficient γ_isos in Eqs. 11 depends on the attributes of the isosbestic point ( Molecules 13 02260 i084

,

), the irradiation condition (l_irr= l_isos) and F^isos, the photokinetic factor for the isosbestic irradiation (Eq. 3b).

(13)

The overall observed rate constants of the photochemical and the thermal reactions are related to the exponential factors in Eqs. 10, as follows [90]

(14a,b)

The cumulative absorbances (M) of the species A and B at an arbitrary observation wavelength considered at init = 0, PSS and STE, are respectively defined by

(15)

Similarly, the cumulative initial velocities of the various thermal (Eq. 16a) and photochemical (Eq. 16b) reactions are given by

(16a)

(16b)

The equilibrium concentrations of the species (that might be used in Eqs. 10, 15 and 16) are defined on the basis of the respective kinetic coefficients a_i (Eqs. 11 and 12) and the total initial concentration.

(17a,b)

(17c.d)

On the basis of Eqs. 11-17, we can rearrange Eqs. 10 so that a general formula (Eq. 18) is set for the time dependence of the cumulative absorbances (i.e., where the two species are simultaneously monitored at the observation wavelength; with symbols Ω = P or T, and θ = isos for photoreactions).

(18)

Eqs. 10-18 can be used for any sequence in Table 2 provided that the missing parameters and coefficients in each kinetic reaction are given the zero value. They also represent the model equations employed to fitting simulated and experimental traces.

7. Data Simulation, Fitting of the Kinetic Traces and Testing the Methods

Data simulation

Numerical integration methods (NIM) are important tools to calculate the values of integrals, which are specially crucial for non-solvable differential equations (e.g. Eqs. 1 and 5). Hence, kinetic traces for any system amongst those listed in Table 2 can readily be constructed by feeding the numerical integration of Eq. 1 with chosen realistic and physically meaningful values of the involved parameters. Such calculations may include cases of constant and variable photokinetic factors as well as photoreactions for which quantum yields are wavelength dependent. The simulated data presented hereafter were obtained with a fifth-order Runge-Kutta numerical integration method within Mathcad premium software (version 2001i).

Fitting of the Kinetic Traces

The data treatment for such kinetics starts with fitting the traces with the appropriate model equations (Eqs. 18). Each simulated and/or experimental trace is fitted and the corresponding fitting coefficients for the model-equation are extracted. These are Molecules 13 02260 i086

(0 or PSS), Molecules 13 02260 i087

and

for thermal reactions, and, Molecules 13 02260 i089

(0 or STE), Molecules 13 02260 i090

and

for photoreactions. The kinetic curves may be increasing or decreasing. The fitting coefficients are used further to carry out the elucidation method corresponding to the considered kinetic system.

A number of available software programmes can readily perform such non-linear curve fittings. In this paper a Levenberg-Marquardt iterative programme within the Origin 6.0 software package were used for the determination of the best fit curves and the extraction of the fitting coefficients.

Testing the Methods

Prior to applying the methods to experimental data, it is necessary to confirm their validity. One way to achieve this confirmation would be to test the methods with data generated by a completely independent mathematical approach, i.e. the numerical integration methods (see data simulation section above). It is a reasonable strategy because the mathematically analytical elucidation methods should stand against any independently generated data like those provided by NIMs. The advantage here is the perfect knowledge of all reaction attributes. Indeed, the method is accepted as valid if it is capable to retrieve those values of the attributes originally used for NIM.

Furthermore, as NIMs afford the possibility of simulating a great number of reaction cases and experimental conditions, they can represent an efficient tool for testing the elucidation methods in a large set of possible situations. This has advantage over the limited number of data afforded by experiments on a particular reaction case. Nonetheless, after the methods have been positively tested against NIM data, they should be applied, in a subsequent stage, to particular experimental results.

The testing procedure is as follows: let set A be the set of parameter values (ϕ, k and ε) that has been used to feed the numerical integration of Eq. 1 and led to the generation of photochemical and thermal kinetic traces (e.g. Tr_P and Tr_T, respectively). Then, in order to achieve the conditions of a realistic experiment, Tr_P and Tr_T together with C₀, I₀, l_probe and l_irr are considered as observables of known values whereas set A will represent the experimentally inaccessible set of unknowns that the method must define (i.e. set A, including the values of ϕ, k and ε, will remain unknown during the rest of the treatment). If data are obtained by an isosbestic irradiation, the traces Tr_P and Tr_T are fitted to the appropriate model equations and their respective fitting coefficients determined. In the case of non-isosbestic irradiation, the equilibrium absorbances and/or the initial velocities are determined. The adequate spectrokinetic method is then applied and an elucidation solution is found. The latter elucidation solution, corresponding to a set of values of the unknown parameters, is called set B.

Since the spectrokinetic methods have been built from mathematically analytical formalisms, sets A and B should incorporate the same value for each parameter. A high degree of precision should characterize these values (theoretically they should be exactly identical). Otherwise, any discrepancy between the values of sets A and B would strongly indicate that the spectrokinetic method at hand is flawed and therefore cannot be considered as an elucidation method for real kinetic cases.

It is worth noting that this strategy differs from that using the experimental data and traces as a control set while the parameter values feeding the RK integration are changed until experimental and calculated data show a good agreement. The difference lies in the fact that the “good agreement” found in such a way may only represent a convergence to a local minimum or to one of several identical minima. Therefore, this classical approach does not necessarily lead to a unique kinetic solution, but more importantly, it does not generate the true solution as an exclusive output (several solutions are found if several identical minima are possible). This issue is developed further in Section 8.

8. Distinguishability and Identifiability Problems

The purpose of a kinetic elucidation method is to determine with certainty the reaction model and the values of the unknown kinetic and spectroscopic parameters defining the dynamics. This means that the only successful elucidation methods are those able to properly address and lift any ambiguity relating to the uniqueness of both the specific attributes of a kinetic problem and the mechanism governing its reaction [96,97,98,99,100].

Typically, ambiguity occurs when different sets of parameters’ values and/or reaction models lead to an equal goodness of fit between the experimental data and the chosen kinetic model. The root cause for such a situation is the excess number of unknown parameters compared to the number of the linearly independent equations describing the kinetics. The problem is even more acute if the latter equations are non-linear (even though they are linearly independent). As a consequence, the experimental data collected on a reaction dynamics, and which is used for the purpose of kinetic elucidation, can be interpreted in different ways, i.e. leading to a number of kinetic solutions sharing a similar probability of occurrence. The total number of such solutions can exceed the number of equations or even be infinite. This set of kinetic solutions represents what is called a degenerate solution [98,100].

The degeneracy of the solution can be seen through what is generally known as the distinguishability and identifiability problems.

In the cases where more than one reaction scheme reproduces the experimental data with equal goodness of fit, the system is said to be non-distinguishable. For example, kinetics relative to schemes S₃ and S₇ which involve respectively two and three reaction steps, are non-distinguishable from one another if not solved analytically. Indeed, both are represented by the same model Eqs.18, i.e. their experimental traces can both be fitted by the same Eqs. 18 and in both cases the corresponding fitting coefficients are extracted (§ 7). Therefore, at such a simple analysis level (i.e. if the kinetic analysis is only limited to observing a good fit of the traces with Eqs. 18), the distinguishability problem will remain undoubtedly unsolved. The requirement for more powerful elucidation methods is then required.

The identifiability problem is raised with respect to a given reaction scheme for which more than one set of values for the unknown parameters, enables the kinetic model to fit the experimental data with high accuracy. The kinetics of S₁₁ is an obvious example for identifiability when only kinetic data is used (§ 9).

Therefore, a successful elucidation method is not the one that achieves a unique solution (i.e. any option amongst the set of the degenerate solution), but is the one capable of finding the kinetic option that stands for the true solution. Hence, solving the identifiability/distinguishability problems is a central issue for such kinetic studies.

In this respect, the elaboration of elucidation methods that allow singling out the true solution is fundamental for studying the kinetics of AB systems.

For most spiropyrans kinetics and similar AB systems (Table 2), the question is whether fully efficient kinetic methods are achievable? The fulfilment of this goal will require that the elucidation procedures are built upon assumption free strategies, and employing tight formalisms. The methods presented in the next section, have been developed with these requirements in mind.

The testing strategy based on NIM simulations (§ 7) represents a good way of tackling identifiability and distinguidshability problems. When the true solution of a kinetic system is unequivocally met with such tests, it means that the identifiability problem is solved and the reliability of the method employed is proven.

9. The Spectrokinetic Elucidation Methods

The spectrokinetic methods presented below have been obtained through tight mathematical analyses, based on the equations of Section 6. This set of equations is sufficient to reach analytical kinetic solutions for a majority of the reaction cases revealed in Table 1 and which involve, for each reaction, the maximum number of possible unknowns (Table 2). For the purpose of clarity, only the main results will be reported here. The detailed formalisms of the elucidation methods can be found in the original papers [92,93,94,95]. The variability recorded on the number of both unknowns and accessible kinetic traces for the various reactions, makes that each reaction sequence requires a particular elucidation method, as has been the case in the literature. So far, a unique elucidation method capable of dealing with all kinetic cases and solving their identifiability/ distinguishability problems, is yet to be invented.

AB(1ϕ) reactions, S₁

For the unimolecular photoreaction, the progress of the reaction is induced by the irradiation of the sample with a light that selectively hits species A (and not B), i.e. for this reaction case a non-isosbestic irradiation is necessary.

As pointed out above (§ 6), the fundamental differential equation Eq. 5 has been solved for S₁-type reactions [93]. The closed-form integration of Eq. 5 is achieved with a variable separation method after proceeding with a change of variables. The irradiation and observation wavelengths (λ_irr = λ_obs) must be equal. The resulting kinetic “Log-exp” model is [93]:

(19a)

where, Log is the base 10 logarithm and the coefficient Molecules 13 02260 i092

represents the wavelength-dependent overall photoreaction rate, defined as

(19b)

with Ln the neperian logarithm. The parameter Molecules 13 02260 i092

is the fitting coefficient of Eq. 19a. The photoreaction efficiency, Molecules 13 02260 i093

, is the unknown and all the remaining quantities in Eqs. 19 are experimentally accessible.

A good agreement characterises the results of the model when tested against simulated RK data. It fits well the simulated kinetic traces for the disappearance of the reactant (A) while the same set of RK data is poorly reproduced by a monoexponential model (Figure 4).

Figure 4. Comparison between fitting simulated RK data to (a) monoexponential and (b) analytical model Eq.19a [reprinted from Ref. 93 with permission of the Journal of Photochemistry and Photobiology: A Chemistry].

Also, Eq. 19a has been successfully applied to the kinetic transformation of the open form of a photochromic diarylethene derivative into its closed form using visible light (Scheme 5). Both isomers are thermally stable. Photoirradiations at two wavelengths, 437 and 517 nm, were performed separately on the open isomer. The kinetic traces for both processes were fitted with Eq. 19a. The good agreement between the experimental data and the model is shown in Figure 5 [93].

The determination of the fitting parameters ( Molecules 13 02260 i094

and

) has allowed the determination of the absolute value of the quantum yield (i.e. the unknown system) at each irradiation wavelength, Molecules 13 02260 i096

= 0.28 and Molecules 13 02260 i097

=0.24.

Scheme 5. Unimolecular phototransformation of the closed into the open form of a diarylethene.

Figure 5. Photokinetic traces in hexane solution (1.82 10^-5 M (ο) and 3.28 10^-5 M (□); 15°C) at two irradiation wavelengths in the visible. Experimental data (circles and squares) are readily fitted by the theoretical model, Eq.19a (lines) [reprinted from Ref. 93 with permission of the Journal of Photochemistry and Photobiology: A Chemistry].

The model is believed to be useful to the classical example of the photoreduction of potassium ferrioxalate, whose kinetics is considered to occur via a single photochemical step (Fe³⁺―hv→ Fe²⁺). So far, the experimental data of the ferrioxalate actinometery were usually analysed by Runge-Kutta or Simpson’s numerical integration methods [2,3].

AB(1k) Systems, S₂

This is a classical first order thermal reaction which is fully described by a monoexponential model Eq. 20 (extracted from Eq. 18). It is a typical AB reaction that involves a single thermally active species (either A or B). The reaction corresponds to the thermal relaxation of one isomer in the dark, which ends when the thermally stable species had been fully produced in the reaction medium.

(20a)

Both increasing and decreasing traces are readily fitted by the model equation (Eqn. 20a). Since the exponential coefficient a_19T, which is either equal to (– k_AB) for A→B transformations or (– k_BA) for B→A reactions (Eq. 14b), is easily obtained from the curve fitting, the determination of the reaction rate constant is straightforward.

Figure 6. Examples of unimolecular thermal reaction traces.

(a) The reaction starts from an initial concentration of the reactive species, and (b) it starts from a state where A and B are already in a photoequilibrium.

In realistic situations however, the recording of the early stages of the reaction is difficult to achieve in solution if the initial thermochromic material (A) is the species monitored (because most often short lifetimes characterise such thermochromic interconversions).

In the cases of sequences S₃, S₇ and S₈ (Table 2), the thermal relaxation of the species starts from the PSS (Figure 6b) and hence a similar formula to Eq. 20a is used to fit the data but where the initial absorbance is replaced by that observed at the PSS (Eq. 20b and Figure 6b).

(20b)

AB(1ϕ,1k) Systems, S₃

This reaction involves a phototransformation reversed by a thermal relaxation (which are characterized by a quantum yield ϕ_AB and a rate constant k_BA, respectively). The reaction progress, under irradiation, leads to a PSS from which the thermal relaxation takes place in the dark (Figure 7). When an isosbestic irradiation is used, the model equation for the photochemical reaction is [92]

(21a)

The fitting of the photochemical and thermal traces to Eqs. 21a and 20b, respectively, yield the following fitting parameters: Molecules 13 02260 i098

,

,

, a_19T and Molecules 13 02260 i101

.

Scheme 6. Reaction of spiropyran NO₂-BIPS.

Figure 7. Kinetic traces of the photocolouring and the thermal fading reactions of NO₂-BIPS (Scheme 6) in ethyl acetate (plots represent Absorbance vs. Time (in s)). The experimental data (squares) are fitted by model Eqs. 21a and 20b (solid lines). T = 14˚C, C₀ = 5.83 10^-5 M, λ_isos = 339 nm, λ_obs = 580 nm [reprinted from Ref. 92 with permission of the International Journal of Chemical Kinetics].

The elucidation of this kinetics can be achieved by using only pure kinetic data. The true solution (consisting of the absolute values for the three unknowns, Table 2) as well as the full spectrum of transient species B can be worked out as follows.

The value of the thermal reaction rate constant is obtained from the fitting coefficient a_19T (which is equal to (– k_BA) as described for data treatment corresponding to S₂ kinetics), and the absolute value of the quantum yield is obtained as

(21b)

Once these two parameters ( Molecules 13 02260 i102

and k_BA) are known, the calculation of species concentrations at PSS is achieved on the basis of Eqs. 17c and 17d. This allows the determination of the exact value of the molar extinction coefficient of the transient species ( Molecules 13 02260 i103

at the observation wavelength) from the equation and value of Molecules 13 02260 i100

(Eq. 15, with init = PSS). Accordingly, the kinetic elucidation is completed.

The reconstruction of the spectrum of species B is achieved by applying either terms of Eqs. 21c to different observation wavelengths throughout the absorption region (Figure 8). (λ is an arbitrary wavelength that takes any value where A and B absorb).

(21c)

Figure 8. Absorption spectra of species A (plain line) and B (dashed line) for Scheme 6. The spectrum of B has been reconstructed using data from the elucidation method and Eq.21c. (ε in L.mol^-1.cm^-1 and λ in nm) [reprinted from Ref. 92 with permission of the International Journal of Chemical Kinetics].

AB(1k,1ϕ) Systems, S₄

The elucidation of this kinetic case is much simpler than that of S₃ because species B can be obtained in a stable form and hence its spectra can be recorded. When species B is the coloured form [78], species A can be regenerated by applying a non-isosbestic irradiation to B and the efficiency of the photoreaction at a given irradiation wavelength ( Molecules 13 02260 i104

) is determined as per S₁. Once the phototransformation of B is complete (i.e. Molecules 13 02260 i105

) and while keeping the photoirradiation light on, the full spectrum of species A can be recorded. The thermal rate constant k_AB is subsequently obtained by fitting the thermal relaxation trace of species A (recorded in the dark) to the model Eq. 20a (the same treatment used for S₂ applies). The photochemical efficiencies at wavelengths where both species absorb can be derived from Eq. 6 (where Molecules 13 02260 i104

= k_BA = 0).

In the case where species B is the colourless form, the thermal transformation of species A into B will generate a similar situation to that described for S₃ (i.e. the elucidation method developed for the kinetic case S₃ applies except that the labelling of the equations must take into account that now the initial species is B).

The elucidation of this kinetic case is simpler than the preceding one (S₃) because the full spectra of both species are easily accessible.

AB(2k) Systems, S₅

The systems involving pure and opposed thermal reactions cannot be solved on the basis of kinetic data alone if the early stages of the starting reaction are not accessible (S₂ kinetics and Table 2). Alternative physical data can be considered. In this instance, the determination of the equilibrium concentration of one or both species can be obtained, for example, by NMR. This allows the calculation of the equilibrium constant, K_e = k_AB/k_BA (as equal to C_B(STE)/C_A(STE)). However, the determination of the individual rate constants is impossible if no complete kinetic traces (at least one corresponding to either A , B or A and B) are available. This issue is discussed further with relation to the kinetics of thermo-photochromes (S₁₁).

Figure 9. Examples of conmplete kinetic traces for AB(2k). (a) Observation of individual traces of species A and B. (b) Observation of the temporal evolution of the total absorbance of the same medium as in (a).

Reversible Photochemical AB(2ϕ) and AB(2ϕ,1k) Systems, S₆, S₇ and S₈

The two sequences S₆ and S₇ are presented together in this section because they are solved by the same kinetic elucidation method [94]. We consider here the typical spectral features of photochromes where the spectra of the species overlap throughout the UV region while B has an individual absorption in the visible range.

There are three kinetic situations relating to the properties of the photoreactions that can be reviewed separately. The first kinetics concerns the sequence involving opposed pure photochemical reaction steps (S₆) where each step is exclusively sensitive to either UV (e.g. Molecules 13 02260 i106

) or visible light (e.g. Molecules 13 02260 i104

, Scheme 7). The second case is that of sequences S₆ and S₇ whose photoreactions are simultaneously initiated by UV light and are independent of the irradiation wavelength. In the third kinetics we will consider the preceding case but with forward photoreaction efficiencies dependent on irradiation. The much complex kinetics where both direct and reverse quantum yields are irradiation dependent is also briefly discussed.

A simple test has been developed to check the variability of forward quantum yields with irradiation. It is given by a factor α (Eq. 22), that uses only accessible experimental data [94].

(22)

where, the value of Molecules 13 02260 i108

is obtained by differentiation of the non-isosbestic kinetic trace (irr/irr) and Molecules 13 02260 i109

is a fitting parameter of Eq. 21a. The photokinetic factors, F^isos and Molecules 13 02260 i075

, are calculated using Eqs. 3c and 3b, respectively.

For example, the quantum yields at an isosbestic ( Molecules 13 02260 i102

) and a non-isosbestic ( Molecules 13 02260 i093

) UV-irradiation wavelength are identical if α is equal to unity. Otherwise, for α ≠ 1, they are necessarily different.

Pure photochemical opposed reactions AB(2ϕ) which are responsive to different light ranges

For this reaction, one of the species is sensitive to UV irradiation (e.g. species A) and its coloured photoproduct (e.g. species B) is only photoreactive if subjected to visible light (Scheme 7). Each reaction ends by the depletion of the irradiated species and no PSS is observed.

Scheme 7. AB(2ϕ) reactions involving photoreactions activated by different lights.

Since both species spectra can be recorded, the extraction of the photoreactions quantum yields can be performed by using Eq. 21a for the isosbestic UV-irradiation and Eq. 19a for the non isosbestic irradiations in the visible. It is worth noting that Eq. 19a can be rearranged to be applicable when the reaction medium is subjected to visible light (where only B absorbs) and monitored at a wavelength where both species absorb (i.e. in the UV range).

(23)

Even though such a mechanism has been proposed for photochromes [84], it is nonetheless reasonable to question the assumption employed here. Thus, unless tangible evidence is provided, the UV irradiation cannot be supposed to initiate selectively the photoreaction of species A. It is then necessary to take into account the reactivity of both species under UV-irradiation. These cases will be considered in the following sections.

AB(2ϕ) and AB(2ϕ,1k) systems where both photochemical reactions are responsive to the same excitation beam and their quantum yields are independent of the irradiation wavelength (e.g. UV).

The pioneering work of Fischer [8] on the kinetic elucidation of AB(2ϕ) systems (S₆) has demonstrated that the kinetic solution for this case can only be achieved provided that the ratio of the forward to reverse quantum yields, at any two non-isosbestic irradiations, is constant. This constant ratio hypothesis is a reasonable assumption according to Kasha’s rule. However, the mathematical formalism adopted for this method does not hold if an additional thermal reaction step is incorporated in the reaction scheme (S₇). The new approach presented hereafter overcomes this problem and elucidates equally S₆ and S₇ kinetics.

The working assumption for this kinetic case supposes the invariability of the quantum yield values with the irradiation wavelength (α = 1, Eq. 22).

(24a)

The unknowns of each kinetic system are given in Table 2. The elucidation method relies on an isosbestic and a non-isosbestic irradiations in the UV-range. The corresponding kinetic traces, obtained by observing the system at a unique wavelength λ_irr (Figure 10), are used for the treatments.

Figure 10. Examples of simulated (a) AB(1k,2ϕ) kinetic traces using RK integration for the photochemical reaction performed at two excitation conditions (here at exc/irr wavelength, Molecules 13 02260 i110

) and (b) AB(2 ϕ) photoreaction ( Molecules 13 02260 i111

at the observation wavelength) [reprinted from Ref. 94 with permission of Photochemical and Photobiological Sciences].

Figure 10. Examples of simulated (a) AB(1k,2ϕ) kinetic traces using RK integration for the photochemical reaction performed at two excitation conditions (here at exc/irr wavelength, Molecules 13 02260 i110

) and (b) AB(2 ϕ) photoreaction ( Molecules 13 02260 i111

at the observation wavelength) [reprinted from Ref. 94 with permission of Photochemical and Photobiological Sciences].

The elucidation procedure for this kinetic cases starts with the determination of the thermal rate constant (k_BA). It is a fitting parameter (–a_19T) of Eq. 20b to the thermal relaxation trace. The mathematical formalism yields the analytical expression of the extinction coefficient of species B at the irradiation wavelength. The formula of Molecules 13 02260 i105

is exclusively defined by experimentally available coefficients.

(24b)

The previous information ( Molecules 13 02260 i105

, Eq. 24b) allows the determination of the absolute value of the forward quantum yield,

(24c)

which, in turn, serves the determination of the reverse quantum yield.

(24d)

The reconstruction of the whole spectrum of species B ( Molecules 13 02260 i112

) is realised through the calculation of PSS concentrations of the species as given by Eqs. 18c and 18d. These values can then be considered in the expression of the observed total absorbance (Eq.4 at PSS, i.e. Molecules 13 02260 i113

or

) to extract the extinction coefficient values ( Molecules 13 02260 i112

).

(24e)

Table 3. Initial parameters^a experimentally accessible data and fitting results for Fig 10a,b.

**Table 3.** Initial parameters^a experimentally accessible data and fitting results for Fig 10a,b.
Initial parameters
Figure #	C_o(M)	(M^-1 cm^-1)	(M^-1 cm^-1)	(Einst.s^-1 dm^-1)	(Einstein.s^-1 dm^-1)	k_BA (s^-1)
10a	1.9 10^-5	22558	9536	2 10^-6	9.1 10^-7	0.04
10b	8.5 10^-6	18617	58690	6 10^-6	2.0 10^-6	0
Fitting results for the simulated kinetic traces and all experimentally accessible data
Figure #	(0) ^b	^c	^c	γ_isos ^c	(PSS)
10a	0.1812	0.0234	-0.0889	0.0660	0.2527	1.75
10b	0.4988	-0.0673	-0.2651	0.2150	0.1285	1.99

^a Throughout l_irr = l_probe = 1 cm. ^b: Molecules 13 02260 i115

. ^c: These parameters are given in s^-1 [reprinted from Ref. 94 with permission of Photochemical and Photobiological Sciences].

Table 4. Spectroscopic and kinetic values of the unknown parameters calculated by the spectrokinetic method (Eqs. 24) using data of Table 1.

**Table 4.** Spectroscopic and kinetic values of the unknown parameters calculated by the spectrokinetic method (Eqs. 24) using data of Table 1.
Figure #		ϕ_AB	ϕ_BA	C_A(PSS)
Figure #		ϕ_AB	ϕ_BA	λ_isos ^a	λ_irr ^a
1	55000	0.41	0.33	1.32 10^-5	1.74 10^-5
2	10359	0.76	0.47	3.25 10^-6	8.36 10^-7

^a: wavelength of the irradiation beam responsible for yielding C_A(PSS) in solution [reprinted from Ref. 94 with permission of Photochemical and Photobiological Sciences].

The above procedure leads to the true kinetic solution for the cases targeted, for which the independence of the quantum yields from irradiation is known with certainty. If however, for a lack of convincing information, there is doubt, as it happens for a majority of new photochromes and other AB kinetic systems, then the application of the above method may lead to parameter values which have no consistent physical meaning (e.g. a quantum yield value higher than unity or negative). For these cases, the working hypothesis is invalidated; meaning that at least one of the quantum yields is irradiation dependent. Therefore, it is necessary to consider the variability of the quantum yields with irradiation wavelengths.

AB(2ϕ) and AB(2ϕ,1k) systems where both photochemical reactions are responsive to the same excitation beam and their forward quantum yields are dependent on the irradiation wavelength.

The same experimental conditions and a set of two kinetic traces, as described in the preceding spectrokinetic method, are used for the elucidation of this dynamics [94]. For the present case we consider that ϕ_AB is irradiation dependent while ϕ_BA is constant with irradiation wavelength (Eq. 25a), i.e. α ≠ 1 (Eq. 22).

(25a)

The formalism of this kinetic elucidation case yields a second order algebraic equation for the determination of Molecules 13 02260 i105

as the first unknown parameter.

(25b)

where the constants c_i take the formulae below.

(25c)

(25d)

(25e)

with

(25f)

(25g)

(25h)

The solution of Eq. 25b will depend on the value of the quantity, Molecules 13 02260 i116

; there are three possibilities available.

A negative value for Δ means that in the actual conditions there is no possible solution to the problem considered under the working assumption (see next section). However, for positive Δ and/or when Δ is equal to zero, the roots of Eq. 25b are given by,

(26a)

Once the value of the extinction coefficient is obtained, the successive use of the following expressions allows the determination of the remaining unknown parameters and the completion of the kinetic elucidation of this case.

(26b)

(26c)

(26d)

The set of Eqs. 26 represents the kinetic solution of the systems under consideration. For those situations where a solution exists, the complete spectrum of species B is recovered in a similar manner to that proposed in the preceding method (employing Eq. 24e).

The application of this elucidation method to the kinetic data of a diarylethene derivative, DAE, (Scheme 8) irradiated in the UV range at λ_isos = 317 nm and λ_irr = 345 nm (Figure 11), yielded the kinetic solution given in Table 5.

Figure 11. Experimental traces of DAE (5.6 10^-5 M) in hexane at 15 °C. The kinetic traces are labelled by the wavelengths at which the irradiation and the observation were carried out (λ_exc / λ_obs). Two of the curves (317/317 and 317/345, in circles) have been performed using the isosbestic irradiation and the third resulted from a non-isosbestic irradiation (345/345, in squares). The line linking the points of trace 317/345 represents the best fit to the monoexponential model, Eq.21a. The line shown on the trace 345/345 has been calculated by a RK integration using the data supplied by the elucidation method (Table 5) [reprinted from Ref. 94 with permission of Photochemical and Photobiological Sciences].

It is worth noting that the application of the methods developed for the kinetics of S₃, and, S₆ and S₇ with ϕ ≠ f(λ_exc), do not yield consistent results indicating that the system under consideration does not obey AB(1ϕ,1k) or AB(2ϕ,1k) kinetics with irradiation independent quantum yields. Furthermore, the use of the second root obtained from Eq. 26a ( Molecules 13 02260 i117

= 2094 M^-1 cm^-1, since Δ > 0, Molecules 13 02260 i118

is given in Table 5 as Molecules 13 02260 i119

) in Eqs. 26b-26d generates a negative value for ϕ_BA and values over unity for Molecules 13 02260 i102

and

.

Scheme 8. Open and closed forms of DAE.

Table 5. Calculated and experimental spectroscopic and kinetic values for the unknown parameters of the system.

**Table 5.** Calculated and experimental spectroscopic and kinetic values for the unknown parameters of the system.
Parameters	^a			ϕ_BA	^b	^b	^b
Results of the Elucidation Method	8887.8	0.23	0.419	0.095	1.6 10^-5	3.4 10^-5

^a: units: mol^-1·L.cm^-1. ^b: The concentration of the initial species at the PSS (in mol·L^-1); where the irradiation wavelength used for the experiment is also indicated.

Such results prove that the identifiability/distinguishability problem can be readily solved for this kinetics. In general terms it is crucial to review the different possibilities that can emerge from the signs and values of the roots (Δ ≥ 0) in relation to the identifiability and distinguishability issues.

In the case where the spectrokinetic method yields a unique solution, i.e. a single positive root that generates plausible results for the remaining parameters, the latter will represent the true solution for the studied kinetics since there are no identifiability or distinguishability problems. However, if both extinction coefficient values are positive ( Molecules 13 02260 i124

and

> 0), and if both lead to plausible sets ofparameter values then the identifiability and perhaps also (if ϕ_BA is zero for one solution) the distinguishability problems arise because of the occurrence of two options. In this case, the complete solution would require further information to discern the true solution.

AB(2ϕ) and AB(2ϕ,1k) systems where both photochemical reactions are responsive to the same excitation beam and both their forward and reverse quantum yields are dependent on the irradiation wavelength.

The most inclusive kinetic case for AB(2ϕ,1k) systems is the one where both quantum yields are irradiation dependent, i.e. Molecules 13 02260 i106

≠

and

≠

. For this case, where Molecules 13 02260 i093

= f(λ_exc) and Molecules 13 02260 i104

= g(λ_exc), an analytical (assumption-free) method based on the use of pure kinetic data is yet to be developed. The difficulty in finding such a useful kinetic elucidation method stems from the fact that the number of unknowns is higher than that of the linearly independent equations that describe the reactive system. Progress on a better physical description of the kinetics by the integration of the fundamental differential equation of such kinetics would certainly represent a hope to complete this kinetic elucidation story.

Nonetheless, we ought to consider the situation where physical information is available for the system. Indeed, for AB(2ϕ) kinetics it is possible to determine the PSS concentrations of the species (e.g. by HPLC) since reactant and product are both thermally stable. The kinetic elucidation in this case is rather straightforward and requires a single trace obtained either under an isosbestic or a non-isosbestic irradiation.

The spectrum of species B is determined, in a first elucidation step, by using Eq. 24e, its PSS concentration, the spectrum of the initial species and that of the medium recorded at PSS. The forward quantum yield value can subsequently be determined in various ways, for instance, from the equation of the initial velocity; either using Eq. 27a or 27b for an isosbetic and a non-isobestic irradiation, respectively. The numerical value of the latter is obtained by differentiation of the corresponding kinetic trace as indicated above (Eq. 6).

(27a)

(27b)

Finally, the value of the reverse quantum yield can, for example, be extracted from the expression of C_A(PSS) (Eqs. 8 and 17c, respectively for non-isosbestic and isosbestic irradiations).

This procedure may be repeated with data corresponding to a set of irradiations to determine the absolute values of both quantum yields throughout the UV-range. It can be combined with the application of Eq. 19a to achieve a thorough study, by determining the irradiation-dependence of species B quantum yield in the visible range.

AB(2ϕ,1k) systems, S₈

The spectrum of A is supposed to be inaccessible to the experimentalist, in this case. Hence, once it has completely transformed into B the kinetics can be treated as for that of S₇ and the corresponding elucidation method applies.

Thermophotochromic AB (1ϕ,2k) systems, S₉ and S₁₀

The kinetics of S₈ and S₉ are characterised by opposing thermal reaction steps and a single photochemical step responsible for the transformation of either A into B (S₈) or B into A (S₉). The total number and the identity of the five unknown parameters is the same for both kinetics (Table 2) except that each reaction sequence involves a different quantum yield (ϕ_BAreplaces ϕ_AB for S₉). Since the lifetimes of the thermal reactions are relatively short for photochromes, recording individual spectra of initial species or the early stages of the kinetic traces by spectrophotometry is not possible. Therefore, for such kinetics the photochemical reaction is realised after the thermal equilibrium (STE) has been reached. In these circumstances, the two reaction sequences do not differ on principle in such a way that the kinetic elucidation method developed for S₈ [95] hereafter can similarly be used for the treatment of S₉ provided that the appropriate modifications due to the new reaction sequence are taken into account in the labelling of the equations.

The full elucidation of these kinetics is possible using pure kinetic data only and the formalism shows that an optimal number of one photochemical and one thermal traces is sufficient to identify all five unknowns including the spectrum of the initial species.

The kinetic behaviour of AB(1ϕ,2k) systems involves in theory three traces (Figure 12). The first one, is due to the thermal transformation of the initial species until the reaction reaches the state of thermal equilibrium. This stage is labelled Thermochromism 1 (or shortly Th₁). As previously stated, this first trace in generally not experimentally accessible and its specific equations are not considered in the elucidation formalism. The second kinetic trace starting from STE, corresponds to the photochemical transformation due to an isosbestic irradiation (labelled P for photochromism). When light is switched off, the system relaxes thermally from PSS to STE during a last stage called thermochromism 2 (Th₂).

Figure 12. Typical kinetic traces depicting the three successive phases (Th1 (ο), P (▪) and Th2 (+)) of a thermophoto-reactive system (with M_tot the total absorbance of the medium at the observation wavelength; irr. on and irr. off represent respectively the start and end of the isosbestic irradiation) [reprinted from Ref. 95 with permission of the International Journal of Chemical Kinetics].

The kinetic laws obeyed by the useful thermal and the photochemical reactions ( Molecules 13 02260 i126

and

respectively), are:

(28a)

(28b)

The fitting results relative to the overall reactions’ rates are used to determine the absolute value of the quantum yield ( Molecules 13 02260 i106

) by using Eq. 21b. The knowledge of this value facilitates the calculation of the value of the molar extinction coefficient of species B at a given observation wavelength (λ).

(28c)

This expression proves that it is possible to determine the spectrum of the coloured species at any observation wavelength prior to knowledge of the spectrum of the starting material and the thermal rate constants. It also shows that a previous conclusion relative to these kinetic cases relating to the necessity of knowing C_B(STE) in order to determine the value of Molecules 13 02260 i112

[2,3], is not valid.

In order to determine k_AB, the absorbance of species B at the thermal equilibrium ( Molecules 13 02260 i128

) is used. (Note that in the visible region, Molecules 13 02260 i129

= 0).

(28d)

And from there, the value of the reverse thermal rate constant is easily worked out as

k_BA = −a_19T − k_AB

(28e)

Since all three kinetic parameter values have been defined, the four equilibrium concentrations are obtained by the application of Eqs. 17.

Finally, the kinetic elucidation is completed by deriving the extinction coefficients values at different wavelengths for species A; using the spectrum of B, the concentrations at an equilibrium state (either STE or PSS) and the full spectrum reached at that equilibrium state, as follows:

(28f)

Thermophotochromic AB (2ϕ, 2k) systems, S₁₁

The kinetics represented by a thermal and a photochemical equilibrium has not been solved by a mathematically analytical method. The solution was not achieved thus far irrespective of the data used (i.e., using only kinetic or a mixted set of kinetic and physical data).

The reason that hampers devising a mathematically analytical elucidation remains, here as well, an excessive number of unknowns compared to that of the available equations. The serious identifiability problem, which is inherent to this kinetic elucidation, is unfortunately unsolvable.

However, the thermal rate constants can be obtained if kinetic and physical information are combined. For example, let the concentration of the species at STE as predicted for AB(2k) systems (S₅), be obtained by a physical method (e.g. NMR, realised at the same temperature and the same solvent as that of the kinetic study). The kinetic data required are represented by a set of one thermal and one photochemical trace (realised under an isosbestic irradiation). In these conditions, the system will formally behave similarly to Figure 12.

The kinetic traces are fitted with Eqs. 28a and 28b, and their respective fitting coefficients defined. According to Eqs. 17a and 17b, the ratio of STE species concentrations is equal to the ratio of the rate constants, as:

(29a)

Combining Eq.29a and the fitting coefficient a_19T (Eq.14b) will allow the extraction of the absolute values for the individual rate constants of the thermal reactions (k_AB and k_BA). Also, the value of the sum of quantum yields is obtained as follows:

(29b)

For this case, unfortunately, the remaining individual values of quantum yields and extinction coefficients for both species cannot be obtained analytically.

10. Determination of reaction quantum yields at any irradiation wavelength

The knowledge of the possible variability of the quantum yield with irradiation wavelengths is important for both application and theoretical purposes. The methods proposed above for the elucidation of AB kinetics provide the absolute values for the quantum yields at specific wavelengths (generally either at λ_isos or at both λ_isos and λ_irr). This is true for the solved kinetics involving photoreactions (S₁, S₃, S₄, S₆-S₁₀).

Since the kinetic laws for non-isosbestic irradiations are not known with the exception of S₁ (Eq. 19), the determination of the values of the quantum yields at different irradiations must be carried out after the kinetic elucidation is completed in the conditions set in Section 9. In these circumstances, the spectra of both species as well as the thermal rate constants are all available. This is important because if the reaction, subjected to the new non-isosbestic irradiation wavelength (labelled irr_x), is performed with all experimental conditions similar to those that served for the kinetic elucidation (of Section 9) then the photochemical quantum yields values are the only unknowns of the new kinetics. Accordingly, the absolute value of the forward quantum yield ( Molecules 13 02260 i130

) can be derived by several ways; one of which being the use of the initial velocity Eq. 6, as:

(30)

For reaction sequences involving a second photochemical step (S₆–S₈ and S₁₀), the value of the reverse quantum yield ( Molecules 13 02260 i131

) can be worked out from the equation of the total absorbance at the PSS (Eq. 9). Hence, individual dependence of the efficiencies ( Molecules 13 02260 i130

and

) on the irradiation wavelength can easily be studied.

11. Applicability of the kinetic elucidation methods to non-chromic systems

Those methods described in Section 9 that do not require in their elucidation procedure data relative to a chromic species (e.g. species B absorbing alone in the visible) are applicable to the elucidation of the kinetics of non-chromic systems.

12. Conclusions

The paper reviews all possibilities of kinetic elucidations relating to AB systems. Simple and easy-to-implement procedures are provided. The elucidation methods presented in this paper tackle and in most cases solve the complex identifiability/distinguishability problem. Their performance and usefulness for photochromes and non-chromic materials have also been shown through a number of examples.

13. Glossary

Labelling

AB (nk,m)	Bimolecular kinetic system
A	Species A (initial reactant)
B	Species B (thermo- or photoproduct)
AB	Transformation of species A into B
BA	Transformation of species B into A
λ	An arbitrary wavelength
obs	Observation wavelength (_obs)
irr	Irradiation wavelength (_irr)
isos	Isosbestic-point wavelength (_isos)
0	Index or argument corresponding to time zero of the reaction
PSS	Index or argument relative to the reaction at the Photo-Stationary State
STE	Index or argument relative to the reaction at the State of thermal Equilibrium
init	Index or argument equal to 0, STE or PSS
P	Photochemical process
T	Purely thermal process
Ω	Index standing for either P or T letters
f(λ_exc)	A function dependent on the excitation wavelength
g(λ_exc)	A function dependent on the excitation wavelength

Irradiation and observation conditions

λ_exc	An arbitrary wavelength used to irradiate the sample; it includes _isos which is the wavelength of an isosbestic point and _irr (≠ _isos) an irradiation wavelength where (and both ≠ 0).
λ_obs	An arbitrary observation wavelength (which might be equal to _irr)
	Incident light intensity of the excitation beam at λ_ιρρ( $I_{0}^{irr}$ ) or _isos ( $I_{0}^{isos}$ )
l_probe	Optical path length of the spectrophotometer (probing) light inside the sample
l_irr	Optical path length of the irradiation light inside the sample

Concentrations

C_i(t)	Concentration of species i = A or B at time t
C_i(PSS)	Concentration of species i = A or B at time PSS
C_i(STE)	Concentration of species i = A or B at time STE
C₀	Total concentration of the species in the reactive medium

Kinetic parameters

	Generic quantum yield of the transformation j when the reaction medium is subjected to an excitation beam whose wavelength is λ_εξχ (= λ_ισοσ or λ_irr).
	Quantum yield of the transformation of A into B at the excitation wavelength λ_exc
	Quantum yield of the transformation of B into A at the excitation wavelength λ_exc
k_BA	First-order rate constant of the thermal transformation of B into A or A into B
F^irr(t)	Photokinetic factor at the irradiation wavelength and time t
F^isos	(or F^isos(t) = F^isos) Constant photokinetic factor at λ_isos
	(or F^irr(0)) Constant photokinetic factor at λ_irr and t = 0
	(or F^irr(PSS)) Constant photokinetic factor at λ_irr at PSS
	(or F^irr(STE)) Constant photokinetic factor at λ_irr at STE
	Initial velocity of reaction Ω when irradiation has been carried out using λ_exc (= isos or irr) and observed at λ_obs (= irr)
a_19Ω	Overall first-order reaction rate for the process Ω (λ_exc = λ_isos)
a_4Ω	First-order direct reaction rate for the process Ω (λ_exc = λ_isos)
a_5Ω	First-order reverse reaction rate for the process Ω (λ_exc = λ_isos)
γ_isos	Constant factor that multiplies quantum yields when λ_exc = λ_isos

Spectroscopic parameters

	Extinction coefficient of species i (A or B) at λ (= isos, irr or obs)
	The total absorbance of the medium measured at various excitation/observation wavelengths’ combinations (e.g. , and respectively at isos/isos, isos/irr and irr/irr)
	Particular constant values of the total absorbance measured at init = 0, PSS, STE

Acronyms

HPLC	High Performance Chromatography
MC	Merocyanine molecule
NMR	Nuclear Magnetic Resonance
RK	Runge-Kutta numerical integration method
SP	Spiropyran molecule
TTC	Trans-Trans-Cis isomer of spiropyran (the same meaning for T and C applies for TTT, CTT and CTC)
UV	Ultraviolet spectrophotometry or range
V_is	Spectrophotometry or spectral range relating to the visible range of the electromagnetic spectrum

Acknowledgements

Financial assistance by the Leicester School of Pharmacy at DeMontfort University is gratefully acknowledged.

References

Bouas-Laurent, H.; Dürr, H. Organic photochromism (IUPAC Technical Report). Pure Appl. Chem. 2001, 73, 639–665. [Google Scholar]
Brown, G.H. Photochromism: molecules and systems; Dürr, H., Bouas-Laurent, H., Eds.; Elsevier: Amsterdam, 2003. [Google Scholar]
Bertelson, R.C. Topics in Applied Chemistry "Organic Photochromic and Thermochromic Compounds"; Crano, J.C., Guglielmetti, R.J., Eds.; Plenum Press: New York, USA, 1999; Volumes 1 and 2. [Google Scholar]
Feringa, B.L. (Ed.) Molecular Switches; Wiley-VCH Verlag: Berlin, Germany, 2001.
Irie, M. Diaryethenes for memories and switches. Chem. Rev. 2000, 100, 1685–1716. [Google Scholar] [CrossRef]
Minkin, V.I. Photo-, thermo-, solvo-, and electrochromic spiroheterocyclic compounds. Chem. Rev. C 2004, 104, 2751–2776. [Google Scholar] [CrossRef]
Kuhn, H.J.; Braslavsky, S.E. Chemical actinometry (IUPAC Technical Report). Pure Appl. Chem. 2004, 76, 2105–2146. [Google Scholar]
Fischer, E. The calculation of photostationary state in systems A B when only A is known. J. Phys. Chem. 1967, 71, 3704–3706. [Google Scholar] [CrossRef]
Flanney, J.B., Jr. The photo- and thermochromic transients from substituted 1 ', 3', 3'-trimethyl-indolinobenzospiropyrans. J. Am. Chem. Soci. 1968, 90, 5660–5671. [Google Scholar] [CrossRef]
Blanc, J.; Ross, D.L. A procedure for determining the absorption spectra of mixed photochromic isomers not requiring their separation. J. Phys. Chem. 1968, 72, 2817–2824. [Google Scholar] [CrossRef]
Tomlinson, W.J.; Chandross, E.A.; Fork, R.L.; Pryde, C.A.; Lamola, A.A. Reversible Photodimerization: a New Type of Photochromism. Appl. Optics 1972, 11, 533–548. [Google Scholar] [CrossRef]
Wyman, G. M.; Zarnegar, B.M. Excited state chemistry of indigo dyes. Fluoescence vs. cis-trans isomerization. J. Phys. Chem. 1973, 77, 831–837. [Google Scholar] [CrossRef]
Wyman, G.M. A new method for calculating the absorption spectra of the components of photochromic systems A B. Mol. Photochem. 1974, 6, 81–90. [Google Scholar]
Appriou, P.; Gugglielmetti, R.; Garnier, F. Étude du processus photochimique impliqué dans la réaction d’ouverture du cycle benzopyrannique. J. Photochem. 1978, 8, 145–165. [Google Scholar]
Chu, N.Y.C. Photochromism of spiroindolinonaphthoxazine. 1. Photopysical properties. Can. J. Chem. 1983, 61, 300–305. [Google Scholar] [CrossRef]
Brauer, H.-D.; Schmidt, R. A new reusable chemical actinometer for UV irradiation in the 248-334 nm range. Photochem. Photobiol. 1983, 37, 587–591. [Google Scholar] [CrossRef]
Brauer, H.-D.; Schmidt, R.; Gauglitz, G.; Hubig, S. Chemical actinometry in the visible (475–610 nm) by meso-diphenylhelianthrene. Photochem. Photobiol. 1983, 37, 595–598. [Google Scholar] [CrossRef]
Rau, H.; Greiner, G. Determination of quantum yields of the X/Y isomerization system from spectroscopic data. Europ. Photochem. Assoc. Newsl. 1991, 41, 40–55. [Google Scholar]
Hanazawa, M.; Sumiya, R.; Horikawa, Y.; Irie, M. Thermally irreversible photochromic systems. Reversible photocyclization of 1,2-bis (2-methylbenzo[b]thiophen-3-yl)perfluorocyclocoalkene derivatives. J. Chem. Soc. Chem. Commun. 1992, 206–207. [Google Scholar]
Sheah, J.R.; Diamond, J.J. Measurement of quantum yields for photochemical systems interconverted by light and heat. J. Phys. Chem. 1992, 96, 4938–4944. [Google Scholar] [CrossRef]
Gauglitz, G.; Scheerer, E. Method for the determination of absorption coefficients, reaction rate constants and thermodynamic data in the system A B. J. Photochem. Photobiol. A: Chem. 1993, 71, 205–212. [Google Scholar] [CrossRef]
Favaro, G.; Malatesta, V.; Mazzucato, U.; Ottavi, G.; Romani, A. Thermally reversible photoconversion of spiroindoline-naphtho-oxazines to photomerocyanines: a photochemical and kinetic study. J. Photochem. Photobiol. A: Chem. 1995, 87, 235–241. [Google Scholar] [CrossRef]
Gauglitz, G.; Stoβel, P.; Meier, H.; Rau, H. Photokinetic examination of (Z,E,E)-4,4′-distyrylazobenzene. J. Photochem. Photobiol. A: Chem. 1995, 85, 207–211. [Google Scholar] [CrossRef]
Gade, R.; Porada, Th. Determination of quantum yields and product spectra by studying reversible photoisomerization in solution. J. Photochem. Photobiol. A: Chem. 1997, 109, 27–34. [Google Scholar] [CrossRef]
Pimienta, V.; Froute, C.; Deniel, M.H.; Lavabre, D.; Guglielmetti, R.; Micheau, J.C. Kinetic modelling of the photochromism and photodegradation of a spiro[indoline-naphthoxazine]. J. Photochem. Photobiol. A: Chem. 1999, 122, 199–204. [Google Scholar] [CrossRef]
Yamashita, T. Quantitative Kinetic Analysis of Photochromic Reactions. J. Photopolym. Sci. Technol. 1999, 12, 257–262. [Google Scholar] [CrossRef]
Strashnikova, N.; Papper, V.; Parkhomyuk, P.; Likhtenshtein, G.I.; Ratner, V.; Marks, R. Local medium effects in the photochemical behavior of substituted stilbenes immobilized on quartz surfaces. J. Photochem. Photobiol. A: Chem. 1999, 122, 133–142. [Google Scholar] [CrossRef]
Matsuda, K.; Irie, M. Photochromism of diarylethenes having nitronyl nitroxides. Tetrahedron Lett. 2000, 41, 2577–2580. [Google Scholar] [CrossRef]
Häupl, T.; Zimmermann, T.; Hermann, R.; Brede, O. On the Photochromism of Spiro[cyclohexadiene-dihydroacridines]. Photochem. Photobiol. 2000, 71, 294–299. [Google Scholar] [CrossRef]
Brookman., J.; Chacon, J.N.; Sinclair, R.S. Some photophysical studies of cis- and trans-urocanic acid. Photochem. Photobiol. Sci. 2002, 1, 327–332. [Google Scholar] [CrossRef]
Tonne, J.; Prinzbach, H.; Michl, J. Double actinometry for rigorous evaluation of quantum yields of clean photoreversible photochemical reactions. Photochem. Photobiol. Sci. 2002, 1, 105–110. [Google Scholar] [CrossRef]
Atabekyan, L.S.; Zakharova, G.V.; Ogienko, V.N.; Chibisov, A.K. Photochromic transformations of spiro compounds. 2. The kinetics of photocoloration under continuous irradiation. High Energ. Chem. 2002, 36, 322–325. [Google Scholar] [CrossRef]
Sheepwash, M.A.L.; Ward, T.R.; Wang, Y.; Bandyopadhyay, S.; Mitchell, R.H.; Bohne, C. Mechanistic studies on the photochromism of [e]-annelated dimethyldihydropyrenes. Photochem. Photobiol. Sci. 2003, 2, 104–112. [Google Scholar] [CrossRef]
Sun, F.; Zhang, F.; Zhao, F.; Zhou, X.; Pu, S. Preparation and photochromic properties of diarylethene nanoparticles. Chem. Phys. Lett. 2003, 380, 206–212. [Google Scholar] [CrossRef]
Belfield, K.D.; Bondar, M.V.; Liu, Y.; Przhonska, O.V. A new blue light-emitting oligofluorene glass: synthesis, characterization and photophysical properties. J. Phys. Org. Chem. 2003, 16, 69–78. [Google Scholar] [CrossRef]
Airinei, A.; Barboiu, V.; Rusu, E.; Timpu, D. Photolysis of 4-azidodiphenyl in solution studied by UV absorption spectroscopy. J. Photochem. Photobiol. A: Chem. 2004, 162, 579–584. [Google Scholar] [CrossRef]
Zelentsov, S.V.; Simdyanov, I.V. Simulation of exposure of light-sensitive layers: Processes involving reversible reactions. High Energ. Chem. 2005, 39, 179–183. [Google Scholar] [CrossRef]
Borisenko, V.; Woolley, G.A. Reversibility of conformational switching in light-sensitive peptides. J. Photochem. Photobiol. A: Chem. 2005, 173, 21–28. [Google Scholar] [CrossRef]
Jeong, Y.-C.; Yang, S.I.; Kim, E.; Ahn, K.-K. Development of highly fluorescent photochromic material with high fatigue resistance. Tetrahedron 2006, 62, 5855–5861. [Google Scholar] [CrossRef]
Pu, S.; Yang, T.; Li, G.; Xu, J.; Chen, B. Syntheses and properties of new photochromic diarylethene derivatives having a pyrazole unit. Tetrahedron Lett. 2006, 47, 3167–3171. [Google Scholar]
Attia, M.S.; Khalil, M.M.H.; Abdel-Mottaleb, M.S.A.; Lukyanova, M.B.; Alekseenko, Y.A.; Lukyanov, B. Effect of complexation with lanthanide metal ions on the photochromism of (1,3,3-trimethyl-5’-hydroxy-6’-formyl-indoline-spiro2,2_-[2H]chromene) in different media. Int. J. Photoenergy 2006, 1–9. [Google Scholar]
Corredor, C.C.; Belfield, K.D.; Bondar, M.V.; Przhonska, O.V.; Hernandez, F.E.; Kachkovsky, O.D. One- and two-photon photochemical stability of linear and branched fluorene derivatives. J. Photochem. Photobiol. A: Chem. 2006, 184, 177–183. [Google Scholar] [CrossRef]
Bossi, M.L.; Rodriguez, J.B.; Aramendia, P.F. Photoinduced changes of absorption and circular dichroism in a chiral nematic phase containing a photochromic fulgide. J. Photochem. Photobiol. A: Chem. 2006, 179, 35–41. [Google Scholar] [CrossRef]
Andrey, Kh.V.; Menshykau, D. Kinetics of photochemical reactions in optically dense media with reagent diffusion. J. Photochem. Photobiol. A: Chem. (in press).
Mauser, H.; Gauglitz, G. Photokinetics: Theoretical Fundamentals and Applications. In Comprehensive Chemical Kinetics; Compton, R.G., Hancock, G., Eds.; Elsevier: Amsterdam, 1998; Volume 36. [Google Scholar]
Görner, H. Photochromism of nitrospiropyrans: effects of structure, solvent and temperature. Phys. Chem. 2001, 3, 416–423. [Google Scholar]
Atabekyan, L.S.; Zakharova, G.V.; Chibisov, A.K. Photochromic transformations of spiro compounds: 1. A new approach to determining the photocoloration quantum yield. High Energ. Chem. 2001, 35, 349–388. [Google Scholar] [CrossRef]
Kießwetter, R.; Pustet, N.; Brandl, F.; Mannschreck, A. 1′,3′,3′-Trimethyl-6-nitrospiro[2H-1 benzopyran-2,2′-indoline]: its thermal enantiomerization and the equilibration with its merocyanine. Tatrahedron Asymmetr. 1999, 10, 4677–4687. [Google Scholar] [CrossRef]
Bletz, M.; Pfeifer-Fukumura, U.; Kolb, U.; Baumann, W. Ground- and first-excited-singlet-state electric dipole moments of some photochromic spirobenzopyrans in their spiropyran and merocyanine forms. J. Phys. Chem. 2002, 106, 2232–2236. [Google Scholar]
Zhang, J.S.; King, J.; Harris, C. Ultrafast studies of photochromic spiropyrans in solution. J. Am. Chem. Soc. 1992, 114, 10921–10927. [Google Scholar] [CrossRef]
Sheng, G.; Leszczynski, J.; Garcia, A.A.; Rosario, R.; Gust, D; Springer, J. Comprehensive Theoretical Study of the Conversion Reactions of Spiropyrans: Substituent and Solvent Effects. J. Phys. Chem. B 2004, 108, 16233–16243. [Google Scholar]
Cottone, G.; Noto, R.; La Manna, G. Theoretical study of spiropyran–merocyanine thermal isomerization. Chem. Phys. Lett. 2004, 388, 218–222. [Google Scholar] [CrossRef]
Hobley, J.; Pfeifer-Fukumura, U.; Bletz, M.; Asahi, T.; Masuhara, H.; Fukumura, H. Ultrafast Photo-Dynamics of a Reversible Photochromic Spiropyran. J. Phys. Chem. 2002, 106, 2265–2270. [Google Scholar]
Wohl, C. J.; Kuciauskas, D. Excited-state dynamics of spiropyran-derived merocyanine isomers. J. Phys. Chem. B 2005, 109, 22186–22191. [Google Scholar] [CrossRef]
Celani, P.; Bernardi, F.; Olivucci, M.; Robb, M.A. Conical Intersection Mechanism for Photochemical Ring Opening in Benzospiropyran Compounds. J. Am. Chem. Soc. 1997, 119, 10815–10820. [Google Scholar] [CrossRef]
Hunt, P.A.; Robb, M.A. Systematic Control of Photochemistry: The Dynamics of photoisomerization of a model cyanine dye. J. Am. Chem. Soc. 2005, 127, 5720–5726. [Google Scholar] [CrossRef]
Sheepwash, M.A.L.; Mitchell, R.H.; Bohne, C. Mechanistic insights into the photochromism of trans-10,10c-methyl-10b,10c-dihydropyrene derivatives. J. Am. Chem. Soc. 2002, 124, 4693–4700. [Google Scholar] [CrossRef]
Xie, N.; Chen, Y.; Yao, B.; Lie, M. Photochromic diarylethene for reversible holographic recording. Mat. Sci. Eng. B 2007, 138, 210–213. [Google Scholar] [CrossRef]
Biedermann, P.U.; Stezowski, J.J.; Agranat, I. Polymorphism versus thermochromism: interrelation of color and conformation in overcrowded bistricyclic aromatic enes. Chem.: A Europ. J. 2006, 12, 3345–3354. [Google Scholar] [CrossRef]
Kaneda, K.; Arai, T. Mechanistic approach to the cyclization reaction of a 2’-hydroxychalcone analogue with light and solvent. Org. Biomol. Chem. 2003, 1, 2041–2043. [Google Scholar] [CrossRef]
Gomes, R.; Parola, J.; Laia, C.A.T.; Pina, F. Efficient photochromism from the network of chemical reactions of 7,4’-dihydroxyflavylium in CTAB micelles. Photochem. Photobiol. Sci. 2007, 6, 1003–1009. [Google Scholar] [CrossRef]
Ballardini, R.; Ceroni, P.; Credi, A.; Gandolfi, M.T.; Maestri, M.; Semararo, M.; Venturi, M.; Balzani, V. Molecular photochemionics. Adv. Funct. Mat. 2007, 17, 740–750. [Google Scholar] [CrossRef]
Wang, P.Y.; Wu, C.J. Photochromic behaviour of some phenoxyanthraquinone dyes in solution and on polyester substrate. Dyes Pigm. 1997, 35, 279–288. [Google Scholar] [CrossRef]
Fu, Z.-S.; Sun, B.-B.; Chen, J.; Yuan, L. Preparation and photochromism of carboxymethyl chitin derivatives containing spirooxazine moiety. Dyes Pigm. 2008, 76, 515–518. [Google Scholar] [CrossRef]
Vicaro, J.; Katsonis, N.; Ramon, B.S.; Bastiaansen, C.W.M.; Broer, D.J.; Feringa, B.L. Nanomotor rotates microscale objects. Nature 2006, 440, 163. [Google Scholar]
Zhao, L.; Hou, Q.; Sui, D.; Wang, Y.; Jiang, S. Multistate/multifunctional switches based on photochromic Schiff base. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2007, 67, 1120–1125. [Google Scholar] [CrossRef]
Bayrakçeken, F.; Sevinç, P.C. Photophysical properties and photochemical rearrangement of single crystal salicylideneaniline in solution. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2007, 66, 184–187. [Google Scholar] [CrossRef]
Tomasuloa, M.; Sortino, S.; Raymo, F.M. A new family of photochromic compounds based on the photoinduced opening and thermal closing of [1,3]oxazine rings. J. Photochem. Photobiol. A: Chem. (in press).
Asiri, A.M.; Ismaiel, N.H. Novel photochromic system derived from tetracyanoquinodimethane and pyrazoles. Pigm. Resin Technol. 2006, 35, 147–150. [Google Scholar] [CrossRef]
Tatewaki, H.; Baden, N.; Momotake, A.; Arai, T.; Terazima, M. Dynamics of Water-Soluble Stilbene Dendrimers upon Photoisomerization. J. Phys. Chem. B 2004, 108, 12783–12789. [Google Scholar] [CrossRef]
Knoll, H. Effects of unidirectional and mutual interactions between microstructures and azo dyes as kinetic probe molecules on cis → trans isomerization rate constants in aqueous P85 and F88 triblock copolymer solutions. Int. J. Chem. Kinet. 2008, 40, 59–65. [Google Scholar] [CrossRef]
Maafi, M.; Aaron, J.J.; Lion, C. Direct valence isomerization of newly synthesized norbornadiene aromatic derivatives. A kinetic and photophysical study. Chem. Lett. 1994, 1865–1868. [Google Scholar] [CrossRef]
Maafi, M.; Aaron, J.J.; Lion, C. Synthesis, photophysical properties and valence photoisomerization of new fluorescent aromatic norbornadienes. New J. Chem. 1996, 20, 559–570. [Google Scholar]
Maafi, M.; Aaron, J.J.; Lion, C. Electrosynthesis and photoisomerization of a new fluorescent aromatic polymer bearing a pendant substituted norbornadiene group. Synthet. Metal 1996, 167–169. [Google Scholar] [CrossRef]
D'Aniello, C.; Musto, P.; Venditto, V.; Guerra, G. Photoisomerization patterns based on molecular complex phases of syndiotactic polystyrene. J. Mat. Chem. 2007, 17, 531–535. [Google Scholar] [CrossRef]
Wang, S.; Choi, M.-S.; Kim, S.-H. Multiple switching photochromic poly(N-isopropylacrylamide) with spironaphthoxazine hydrogel. Dyes Pigm. 2008, 78, 8–14. [Google Scholar] [CrossRef]
Sakai, H.; Ebana, H.; Sakai, K.; Tsuchiya, K.; Ohkubo, T.; Abe, M. Photo-isomerization of spiropyran-modified cationic surfactants. J. Colloid Interface Sci. 2007, 316, 1027–1030. [Google Scholar] [CrossRef]
Keum, S.-R.; Roh, S.-J.; Ahn, S.-M.; Lim, S.-S.; Kim, S.-H.; Koh, K. Solvatochromic behavior of non-activated indolinobenzospiropyran-6-carboxylates in aqueous binary solvent mixtures: part 2. Dyes Pigm. 2007, 74, 343–347. [Google Scholar] [CrossRef]
Li, L.; Yu, M.X.; Li, F.Y.; Yi, T.; Huang, C.H. Inhibit logic gate based on spiropyran sensitized semiconductor electrode. Colloids Surf. A Physicochem. Eng. Asp. 2007, 304, 49–53. [Google Scholar] [CrossRef]
Suh, H.-J.; Lim, W.-T.; Cui, J.-Z.; Leeb, H.-S.; Kim, G.-H.; Heo, N.-H.; Kim, S.-H. The crystalline-state photochromism, thermochromism and X-ray structural characterization of a new spiroxazine. Dyes Pigm. 2002, 57, 149–159. [Google Scholar]
Keum, S.-R.; Choi, Y.-K.; Kim, S.-H.; Yoon, C.-M. Symmetric and unsymmetric indolinobenzospiropyran dimers: synthesis and characterization. Dyes Pigm. 1999, 41, 41–47. [Google Scholar] [CrossRef]
Chernyshev, A.V.; Voloshin, N.A.; Raskita, I.M.; Metelista, A.V.; Minkin, V.I. Photo- and ionochromism of 5’-(4,5-diphenyl-1,3-oxazol-2-yl) substituted spiro[indoline-naphthopyrans]. J. Photochem. Photobiol. A: Chem. 2006, 184, 189–297. [Google Scholar]
Chibisov, A.K.; Görner, H. Complexes of spiropyran-derived merocyanines with metal ions: relaxation kinetics, photochemistry and solvent effects. Chem. Phys. 1998, 237, 425–442. [Google Scholar] [CrossRef]
Moniruzzaman, M.; Fernando, G.F.; Bellamy, A.J. Synthesis and structural investigation of 1’,3’,3’-trimethyl-6-hydroxy-spiro(2H-1benzopyran-2,2’-indoline),1’,3’,3’-trimethyl-6-metha cryloxy-spiro(2H-1-benzopyran-2,2’-indoline) and a copolymer with methyl methacrylate by 1D and 2D NMR spectroscopy. Europ. Polym. J. 2006, 42, 1455–1466. [Google Scholar] [CrossRef]
Li, X.; Li, J.; Wang, Y.; Matsuura, T.; Meng, J. Synthesis of functionalized spiropyran and spirooxazine derivatives and their photochromic properties. J. Photochem. Photobiol. A: Chem. 2004, 161, 201–213. [Google Scholar] [CrossRef]
Toman, P.; Bartkowiak, W.; Nešpůrek, S.; Sworakowski, J.; Zaleśny, R. Quantum-chemical insight into the design of molecular optoelectrical switch. Chem. Phys. 2005, 316, 267–278. [Google Scholar] [CrossRef]
Godsi, O.; Peskin, U.; Kapon, M.; Natan, E.; Eichen, Y. Site effects in controlling the chemical reactivity in crystals: solid-state photochromism of N-(n-propyl)nitrospiropyrane. Chem. Commun. 2001, 2132–2133. [Google Scholar]
Kawanishi, Y.; Seki, K.; Tamaki, T.; Sakuragi, M.; Suzuki, Y. Tuning reverse ring closure in the photochromic and thermochromic transformation of 1’,3’,3’-trimethyl-6-nitrospiro[2H-1-benzopyran-2,2’-indole] analogues by ionic moieties. J. Photochem. Photobiol. A: Chem. 1997, 109, 237–242. [Google Scholar] [CrossRef]
Hadjoudis, E.; Mavridis, I.M. Photochromism and thermochromism of Shiff bases in the solid state: structural aspects. Chem. Soc. Rev. 2004, 33, 279–588. [Google Scholar]
Maafi, M.; Brown, R.G. General analytical solutions for the kinetics of AB (k,ϕ) and ABC(k,ϕ) systems. Int. J. Chem. Kinet. 2005, 37, 162–174. [Google Scholar] [CrossRef]
Maafi, M. Comment on “photokinetic methods: A mathematical analysis of the rate equations in photochromic systems” by G.Ottavi, F.Ortica, and G. Favaro (Int J Chem Kinet, 1999, 31, 303-313). Int. J. Chem. Kinet. 2008, 40, 524–525. [Google Scholar] [CrossRef]
Maafi, M.; Brown, R.G. An analytical solution for the kinetics of AB (1k,1ϕ) systems and its application to a spirobenzopyran. Int. J. Chem. Kinet. 2007, 39, 539–545. [Google Scholar] [CrossRef]
Maafi, M.; Brown, R.G. The kinetic model for AB(1ϕ) systems. A closed-form integration of the differential equation with a variable photokinetic factor. J. Photochem. Photobiol. A: Chem. 2007, 187, 319–324. [Google Scholar] [CrossRef]
Maafi, M.; Brown, R.G. Kinetic analysis and elucidation options for AB (1k, 2ϕ) systems. New spectrokinetic methods for photochromes. Photochem. Photobiol. Sci. (in press).
Maafi, M.; Brown, R.G. A mathematical analysis and elucidation of AB (2k, 1ϕ) thermophoto-chromic kinetics. Analytical solutions using pure kinetic data and considering as unknown the rate constants, quantum yield and both species’ spectra. Int. J. Chem. Kinet. (in press).
Vajda, S.; Rabitz, H. Identifiability and distinguishability of first-order reaction systems. J. Phys. Chem. 1988, 92, 701–707. [Google Scholar] [CrossRef]
Schwaab, M.; Monteiro, J.L.; Pinto, J.C. Sequential experimental design for model discrimination taking into account the posterior covariance matrix of differences between model predictions. Chem. Eng. Sci. 2008, 63, 2408–2419. [Google Scholar] [CrossRef]
Maafi, M.; Brown, R.G. Analysis of diarylnaphthopyran kinetics. Degeneracy of the kinetic solution. Int. J. Chem. Kinet. 2005, 37, 717–727. [Google Scholar] [CrossRef]
Maafi, M.; Brown, R.G. Analysis of diarylnaphthopyran kinetics. Degeneracy of the kinetic solution. Int. J. Chem. Kinet. 2005, 37, 717–727. [Google Scholar] [CrossRef]
Maafi, M.; Brown, R.G. A fundamental system of equations to describe ABC (3k, 6ϕ) dynamics. A contribution to the elucidation of [2H]-naphthopyran kinetics. Int. J. Chem. Kinet. 2008, 40, 269–281. [Google Scholar]

Sample Availability: Samples are not available.

© 2008 by the authors. Licensee Molecular Diversity Preservation International, Basel, Switzerland. This article is an open-access article distributed under the terms and conditions of the Creative Commons Attribution license ( http://creativecommons.org/licenses/by/3.0/).

Share and Cite

MDPI and ACS Style

Maafi, M. Useful Spectrokinetic Methods for the Investigation of Photochromic and Thermo-Photochromic Spiropyrans. Molecules 2008, 13, 2260-2302. https://doi.org/10.3390/molecules13092260

AMA Style

Maafi M. Useful Spectrokinetic Methods for the Investigation of Photochromic and Thermo-Photochromic Spiropyrans. Molecules. 2008; 13(9):2260-2302. https://doi.org/10.3390/molecules13092260

Chicago/Turabian Style

Maafi, Mounir. 2008. "Useful Spectrokinetic Methods for the Investigation of Photochromic and Thermo-Photochromic Spiropyrans" Molecules 13, no. 9: 2260-2302. https://doi.org/10.3390/molecules13092260

APA Style

Maafi, M. (2008). Useful Spectrokinetic Methods for the Investigation of Photochromic and Thermo-Photochromic Spiropyrans. Molecules, 13(9), 2260-2302. https://doi.org/10.3390/molecules13092260

Article Menu

Useful Spectrokinetic Methods for the Investigation of Photochromic and Thermo-Photochromic Spiropyrans

Abstract

1. Introduction

2. Physical properties, reactivity and kinetic study of spiropyrans

Background

Considerations for continuous irradiation experiments

3. AB systems with similar kinetic behaviour to spiropyrans

4. Typical kinetic mechanisms

5. Observables and Experimental Conditions

Photochemical Reactors

Irradiation Conditions

Observables

The unknowns

6. Fundamental Kinetic Laws

The Fundamental Kinetic Equations for AB Systems

Basic Kinetic Laws for AB Kinetics When F is Constant

7. Data Simulation, Fitting of the Kinetic Traces and Testing the Methods

Data simulation

Fitting of the Kinetic Traces

Testing the Methods

8. Distinguishability and Identifiability Problems

9. The Spectrokinetic Elucidation Methods

AB(1ϕ) reactions, S1

AB(1k) Systems, S2

AB(1ϕ,1k) Systems, S3

AB(1k,1ϕ) Systems, S4

AB(2k) Systems, S5

Reversible Photochemical AB(2ϕ) and AB(2ϕ,1k) Systems, S6, S7 and S8

Pure photochemical opposed reactions AB(2ϕ) which are responsive to different light ranges

AB(2ϕ) and AB(2ϕ,1k) systems where both photochemical reactions are responsive to the same excitation beam and their quantum yields are independent of the irradiation wavelength (e.g. UV).

AB(2ϕ) and AB(2ϕ,1k) systems where both photochemical reactions are responsive to the same excitation beam and their forward quantum yields are dependent on the irradiation wavelength.

AB(2ϕ) and AB(2ϕ,1k) systems where both photochemical reactions are responsive to the same excitation beam and both their forward and reverse quantum yields are dependent on the irradiation wavelength.

AB(2ϕ,1k) systems, S8

Thermophotochromic AB (1ϕ,2k) systems, S9 and S10

Thermophotochromic AB (2ϕ, 2k) systems, S11

10. Determination of reaction quantum yields at any irradiation wavelength

11. Applicability of the kinetic elucidation methods to non-chromic systems

12. Conclusions

13. Glossary

Labelling

Irradiation and observation conditions

Concentrations

Kinetic parameters

Spectroscopic parameters

Acronyms

Acknowledgements

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI

AB(1ϕ) reactions, S₁

AB(1k) Systems, S₂

AB(1ϕ,1k) Systems, S₃

AB(1k,1ϕ) Systems, S₄

AB(2k) Systems, S₅

Reversible Photochemical AB(2ϕ) and AB(2ϕ,1k) Systems, S₆, S₇ and S₈

AB(2ϕ,1k) systems, S₈

Thermophotochromic AB (1ϕ,2k) systems, S₉ and S₁₀

Thermophotochromic AB (2ϕ, 2k) systems, S₁₁