Abstract
Several studies of the time varying photon absorption effects, which occur during the photoinitiation process involving in photopolymer materials, have been presented. Three primary mechanisms have been identified: (i) the dye absorption, (ii) recovery, and (iii) bleaching. Based on an analysis of these mechanisms, the production of primary radicals can be physically described and modelled. In free radical photopolymerization systems, the excited dye molecules induce the production of the primary radicals, , which are key in determining how many monomers are polymerized. This, in turn, is closely related to the refractive index modulation formed during holographic recording. In this paper, to avoid the complexities involved in estimating the rate constant of intersystem crossing, , in going from the excited singlet state dye to the excited triplet state dye, we introduce two rates, and these are the proposed rate constants of photon absorption in going from the ground state to the singlet and triplet states, respectively. Using the resulting model, four kinds of Xanthene dyes: Erythrosin B; Eosin Y; Phloxine B, Rose Bengal, and one Thiazine dye: Methylene Blue, are experimentally characterised for use in an AA/PVA photopolymer.
1. Introduction
Photopolymer materials and the photochemical kinetics associated with them [1β3] are being actively studied for practical applications [4β8] such as 3D hybrid optoelectronic circuits, holographic storage [9β13], photoembossing (including the manufacturing of refractive and diffractive optical elements), metrology, 3D displays, and the self-trapping of light [14].
During the photopolymerisation process, the initial step, which involves the absorption of light by a dye, eventually results in the generation of radicals, and this initiation process plays a critical role in grating formation. The theoretical model of the photochemical processes has been developed to include the effects of photosensitizer recovery and bleaching [15β18], in order to permit the accurate prediction of the resulting material behaviour.
In a previous model of this material, in order to simplify the modelling of intersystem crossing from the excited singlet state dye to the excited triplet state dye, a bulk parameter was introduced [18]. It was assumed that during exposure a constant fraction of the ground state molecules present was always available to be converted into the singlet state, while () of them are available to be converted into the triplet excited state. In [19], a more physical representation of intersystem crossing, occurring at a rate constant, , was introduced, but this rate is difficult to estimate by the method described in [19] in that the monomer will be quickly crystallised when triethanolamine is not included in the layer. In order to overcome the limitations of these models, in this paper, we simplify our approach by introducing two new rates, and , into the model, which are the rate constant of photon absorption from ground state to singlet state and triplet state, respectively. Then using this simplified model, the behaviours of five different photosensitizers, in an acrylamide polyvinyl alcohol (AA/PVA) photopolymer material, are studied. Detailed experimental and theoretical analyses are presented.
The paper is organized as follows. In Section 2, we introduce the photochemical reactions that take place and the NPDD model, including descriptions of the processes of initiation, propagation, termination, and inhibition. In Section 3, the material preparation is presented. In Section 4, the photoabsorptive behaviours of four green-sensitive types photosensitizer: Erythrosine B (EB); Eosin Y (EY); Phloxine B (PB), Rose Bengal (RB), of the same molarity in equivalent PVA/AA layers, and one red-sensitive photosensitizer, Methylene Blue (MB), of a higher molarity (due to the weaker absorptivity at the wavelength of 633βnm), also in a PVA/AA layer, are experimentally examined. Based on these results, the recovery and bleaching behaviours of these dyes are also estimated. In Section 5, physical parameter values extracted by fitting the measured results using our model are presented and analysed. Finally, in Section 6, a brief conclusion is given.
2. Theoretical Analysis of Photopolymer Material
In free radical photopolymerization systems, the excited dye molecules induce the production of the primary radicals, the concentration of which is key in determining how much monomers are polymerized. This, in turn, is closely related to the refractive index modulation formed during holographic recording.
In this paper, nonlocal photopolymerization driven diffusion (NPDD) model is developed. This model includes the effects of (a) the kinetics of the major photochemical processes; (b) the temporally and spatially varying photon absorption; (c) the nonlocal material response.
2.1. Photochemical Reactions
In the case of free radical photopolymerization systems, the kinetic model describing what takes place involves four main processes [2, 16, 20β23], (I) initiation, (II) propagation, (III) termination, and (IV) Inhibition. We note that each may involve several physicochemical reactions. We highlight the major chemical reactions in each process [3, 20, 24β27]. A summary of the discussion in this section is given in the flow chart in Figure 1. The excitation of the Xanthenes dyes is similar and is due to the breaking up of the C=O double bond when exposed to the green beam [28]. Since the structure of Methylene Blue is different from the Xanthene dye, the excitation is also different, that is, due to the breaking up of the C=N bond when exposed to the red beam [29]. However the reactions between Methylene Blue and the electron donor are similar to the case of Xanthenes dyes. The detailed chemical scheme of the photopolymerisation is described in [20, 28, 29]. The reactions highlighted within the dashed box will be studied experimentally in Section 4. In particular the introduction of and should be noted when comparing this model to previous models.
(I) Initiation
During illumination, the reaction between the photosensitiser and the electron donor (coinitiator) leads to the production of initiator radicals, , which can react with the monomers to produce chain initiators, .where represents the photosensitizer molecule, indicates the energy absorbed from a photon, , ( and are the quantum efficiencies of the reactions in which the ground state dye is converted into the singlet and triplet dye states, is the dye molar absorptivity, is the thickness of the layer, and Einsteinβs/cm3βs is the incident intensity [27]) are the rate constants of photon absorption in ground state, is the recovery rate of the dye from singlet excited state to ground state, including the processes of (1) radiationless energy transfer to another molecule such as triethanolamine which acts as an electron donor (ED), the emission of a photon by fluorescence [16, 18]. The recovery of triplet state dye is ignored for simplicity and also due to the even lower recovery rate in this process [18, 19], is the electron donor, and is the rate constant of electron donation through which becomes a free radical, , represents a radicalized dye, which abstracted a hydrogen ion from the electron donor, is the rate constant of the photobleaching process, that is, the rate formation of dihydro dye, , is an intermediate form of the electron donor, is the chain initiation kinetic constant, and represents a monomer molecule. We note that (1a) and (1b) are new differing from the models proposed in [18, 19].
(II) Propagation
The chain initiator, will attach itself to another monomer molecule, , by addition to the C=C bond yielding a growing polymer radical with an active tip. Through propagation polymer chain growth then follows [20].
where is the rate constant of propagation and and are the growing macroradical chains of length and () repeat monomeric units ().
(III) Termination
Termination can occur in three ways. Two of these, disproportionation and combination, involve two growing macroradicals interacting, that is, the bimolecular termination mechanism.where and are the rate constants of combination and disproportionation termination respectively. , , and represent terminated chains which have no radical tip, that is, a dead polymer. In this analysis, and will be treated as a single lumped rate constant, (β), as the mode of termination does not effect the polymerization kinetics [25].
A third possible termination mechanism involves primary radical termination [19, 26].
where is the rate constant of primary radical termination. In this step, a growing macroradical chain reacts with a primary radical (initiator radical) leading once again to the production of inactive or dead polymer chains [25].
(IV) Inhibition
Inhibitors are chemicals which react with the initiating and propagating radical species by rapidly removing or scavenging these radicals. Polymerization is completely halted until they are all consumed [24]. Several possible inhibitor reaction mechanisms are listed in the following:where is the inhibitor species, for example, oxygen, is the excited photosensitiser, is the concentration of singlet oxygen [21, 24, 27], and , and are the rate constants of inhibition of the primary radicals [15β17, 21, 30], the photosensitiser and the macroradicals respectively. Inhibition leads to a dead band at the start of exposure, that is, stopping grating formation during the initial exposure. The effects of inhibitors are especially significant when lower exposure energies (dose) are used, for example, when large areas must be exposed (low intensities) or short duration exposures must be used [27].
2.2. Inclusion in the NPDD Model
Based on the analysis present above and summarised in Figure 1, a set of coupled differential equations representing the spatial and temporal evolutions of the material concentrations associated with (1a)β(5c) can be derived using the same notation and following the same methodology as in references [24β26]. where , , and are the concentrations of free-monomer, polymer, and macroradical. , , and are the diffusion rates of monomer, polymer, and inhibitor respectively. Equation (6) arises due to the new proposed reactions in (1a) and (1b). Furthermore, we note that a rate equation governing the total bleached photosensitizer concentration can be obtained: where denotes the total (leuco and dihydro) bleached photosensitizer concentration [19]. (see Figure 1).
In (10), is the nonlocal material spatial response function given by [24]: in which is the constant nonlocal response parameter [20, 31]. The nonlocal spatial response function represents the effect of initiation at a location on the amount of monomer polymerized at location [31].
2.3. Summary
In this section, we have shown how the new dye model can be included in the NPDD model, in which the processes of initiation, propagation, termination and inhibition are included. A method to overcome the limitations of the use of the parameter in [18] and the intersystem crossing rate constant in the previous model [19] is presented. We introduce and into the model equations (1a) and (1b), which are the rate constant of photon absorption from ground state to singlet state and triplet state respectively. This leads to the coupled differential equation (6). This new model simplifies the modelling and data fitting of the experimental data in Section 4.
3. Material Preparation
Photopolymer materials are made sensitive to a particular wavelength using photosensitizing dye. The photosensitisers examined in this paper are that is, Xanthene dyes {Erythrosin B (EB), Eosin Y (EY), Phloxine B (PB) and Rose Bengal (RB)}. These allow holographic recording to be carried out using a 532βnm Solid-State Crystal Laser. We have also sensitized our AA/PVA material in the red (HeNe 633βnm) using Methylene Blue (MB), which is a Thiazine dye.
Our PVA/AA material was made using the components listed in Table 1. The material was prepared as follows:(a)10βg of PVA was added to 100βcm3 of deionised water and dissolved using a heater/stirrer. This solution is then allowed to cool and then 70βcm3 of this solution was transferred into a beaker.(b)8βcm3 of Triethanolamine () was added to the PVA solution and stirred thoroughly.(c)2.4βg of Acrlyamide and 0.8βg of Bis-acrylamide were added to the PVA solution under a fume cupboard and stirred until completely dissolved.(d)16βcm3 of βM Xanthene dye (EB/EY/FB/RB) or 16βcm3 of βM Methylene Blue (MB) was then added to the beaker. As indicated in Section 1, a higher concentration of MB is used due to its weaker absorptivity at the wavelength of 633βnm. This step and subsequent steps were carried out under a safety light, as the material is now sensitive to green/red light. (e)The solution is then made up to 100βcm3 in a volumetric flask with de-ionised water.(f)The solution is then stored in the dark ready for plate preparation.
To prepare dry material layers for holographic exposure the solution prepared above is used as follows:(a)The microscope glass slice on which the material is to be deposited (75βmm Γ 25βmm) is cleaned thoroughly using de-ionised water and Acetone. Once cleaned the plates are placed on a level surface so that the photopolymer layers would adhere to the glass evenly, producing a layer of uniform thickness.(b)1.3β1.4βmLs of the photopolymer solution is then deposited evenly over the area of the glass plate using a syringe.(c)Using this method the typical material thickness is microns. Different thicknesses can be obtained by depositing (drop casting) different quantities of material. The thickness and uniformity of these layers can be measured using a micrometer screw gauge. (d)The plates are then left in the dark for approximately 24 hours until dry. Drying times are dependent on the thickness of the material and the relative humidity.
4. Behaviors of the Five Photosensitizers in AA/PVA
4.1. Transmission Spectra
The dyes examined absorb light and act as photoinitiator, of the polymerization process. Figure 2 shows the transmittance of the four Xanthenes dyes, that is. EB, EY, PB, and RB over a range of visible wavelength.
The normalized transmittance can be expressed as: where is the transmission fraction which corrects for the boundary and scatter losses, is the molar absorption of the dye, is the time varying dye concentration, and is the thickness of the material layer. Taking , (14) reduces to: where is the transmittance value for one particular wavelength, , is presented in Figure 2, for an initial dye concentration βmol/cm3. We define to be the maximum transmittance value. Using this maximum value the molar absorption, , for a wavelength of βnm, that is, the exposing wavelength in all the following experiments, can be estimated.
We also note that the transmittance at a wavelength of 633βnm, for all the dyes is over 59%, taking into account the increased scatter as a result of the crystallisation of the monomer with the time for this case and the thicker thickness of the layer, this means that 633βnm laser light can be used as the probe beam during recording. The resulting parameter values estimated are listed in Tables 2 and 3.
4.2. Transmission
In this subsection, we begin by examining the transmission process of the material. In all cases intensities of 10βmW/cm2 (for EB, EY, PB, and RB), and 4.03βmW/cm2 (for MB) are used to illuminate the layers. First, we need to examine the effects of adding different dyes to the standard AA/PVA material using simple experiments. In all cases the set-ups involve: uniform plane waves of wavelength 532βnm (for EB, EY, PB, and RB) or of 633βnm (for MB), pass normally through the material and the transmitted intensity is measured. The area of illumination is 0.25βcm2. Normalised transmittance curve for the same dye concentrations for EB, EY, PB, and RB, and a higher dye concentration for MB, are given in Figure 3. Using the model proposed in Section 2 dye parameter values are estimated by fitting these curves. Table 4 lists the values of molar absorptivity, , and the quantum efficiency of the reactions in which the ground state dye is converted into the singlet and triplet dye states with the quantum efficiencies and respectively, as estimate by fitting the transmittance curves. We note that the transmittance curves all decrease once exposure ends. This takes place due to the recovery process of the dye once the illuminating light is switched off [18]. βs, βs, βs, βs, and βs are the times when illumination ends for: (a) EB, (b) EY, (c) PB, (d) RB, and (e) MB, respectively. It is worth noting that the values of molar absorptivity, , estimated here agree with those presented in Section 4.1. The value in Section 4.1 is calculated directly by measuring the spectrum and in Section 4.2 the numbers are estimated by fitting the experimental data. However given to the material and experimental environment difference, that is, temperature and humidity, the value could be slightly different. However we can see the agreement is still surprisingly good, and demonstrates that the absorptivity value extracted using our model is reasonable.
4.3. Recovery
In this section we estimate the rate, , at which the photosensitizer recovers back to its initial ground state after it has been excited. In order to determine a value for , a set of experiments is performed, which enables the photosensitizer concentration to be measured at any time, , after a given exposure time, . To do so we use the results presented in Figure 4.
Using the same setup describe in Section 4.2, the recovery process is examined experimentally. The material is illuminated for , after the laser is switched off for time . After [19], the material is again illuminated and the transmittance immediately measured. This measurements repeat several times for different values, so that we obtain the normalised transmission post-exposure, , which is related to the dye concentration post-exposure. As a result, by fitting the experimental data the recovery rate, , can be estimated.
Based on the best fits achieved to the experimental results, using the model described in Section 2, Table 5 shows the estimated values of the recovery rate, . The mean squared error (MSE) values achieved for the fitting procedure are also presented in the table.
4.4. Bleaching
In this section the estimation of the rate constant of photobleaching of the photosensitizer, , the rate constant of electron donation (with which becomes a free radical, ), , and the rate constants of inhibition of the photosensitiser, , are discussed.
The setup for this experiment is the same as that used in Sections 4.2 and 4.3. The material is illuminated for , after the illumination is switched off for . The material is left unilluminated for a long time, in order to allow the dye to fully recover. In practise, we assume that the dye is fully recovered after hours. After hours the illumination is once again switched on and the transmittance measured. This procedure is repeated several times with different values for all the dyes.
The amount of photosensitizer concentration bleached by an exposure of duration, , can be estimated by taking the differences between the initial concentration and , where is the initial photosensitizer concentration and is the final value. A typical set of experimental results are shown in Figure 5, and the resulting parameter values estimated are listed in Table 6.
4.5. Summary
In this section, the photoabsorptive behaviours of five photosensitizers, in almost equivalent PVA/AA layers, are examined using the model derived in Section 2. For each dye, the recovery and bleaching parameters are estimated by fitting reproducible experimental transmittance data.
It can be noted from the quality of the fitting results that the RB is not as well modelled as the other dyes, in that at a certain time , the sum of the concentration of the recovered dye and the bleached dye is less than the original dye concentration value . Furthermore the quality of the fit is poorer as indicated by the fact that the MSE value is relatively higher. One possibility is that this could be a result of our neglect of the recovery of the triplet excited state dye, back into either the ground state, Dye,or the singlet state, . However, following some analyses we note that even if such recovery mechanisms are included into the model the quality of the fit does not improve significantly. There is a connection between the recovery and the bleaching process. Using the same experiment data, as is presented in Figures 6 and 7 it is found that when we get the best fitting for the transmission and recovery of RB, the theoretical predicted values are higher than the experimental values for bleaching, especially for the longest values of time. Table 7 lists the parameter values extracted by applying this procedure. In conclusion, based on these observations, we expect that there are other reactions taking place in the RB system, which are of less significance for the other dyes studied.
5. Simulation Model Predictions
Based on the theoretical analysis and experimental work presented in Sections 2 to 4, the temporal evolution of the primary radical concentration, (or ), can be predicted. It is found that the behaviours of the five dyes are significantly different.
From Section 2, we know that the primary radical concentration, , is produced by the reaction between and , at the rate constant of electron donation, , see (1d). Simultaneously, is also being removed as it initiates photopolymer chains, see (1f), as it reacts with , see (4), and is scavenged by the dissolved inhibitor, see (5a). Applying the full model of the photoinitiation process, using the parameter values estimated above, the temporal evolution of the concentration of the primary radicals, , for (a) EB, (b) EY, (c) PB, (d) RB, and (e) MB, can be predicted using (8).
The result is presented in Figure 8. As can be seen in the initial stages of the exposure, the generation of is significantly inhibited by the original dissolved oxygen, see (5a). Following the resulting initial inhibition period, the generation of progressively increases towards a maximum value. However, after this maximum the concentration value decreases gradually. This is due to the photoinitiation process described by (1f), and the reaction with during the termination process described in (4).
From Figure 8 we can see that for the four Xanthene dyes examined, EB reaches the peak concentration value in a very short time and the peak concentration is also highest, which means that the inhibitor has less effect for the EB case. However the inhibition effect can be more accurately measured using lower exposure intensity, and this is a topic for future work. In addition, for the corresponding MB case, the general behaviour of the concentration is similar to that of the Xanthene dyes.
For PB and MB, since fewer molecules are generated, fewer are available to react with the electron donor, that is, (PB) and (MB) are smaller, thus less is produced. Therefore: (1) the inhibitor scavenging effect takes less time to end; at the same corresponding time in Figure 8, the amount of dihydro dye formed is reduced, that is, there is less bleaching of the PB and MB dyes for the same exposure time, see Figure 5.
Furthermore, comparing the results for these five types of photosensitizer, it is clear that the generation of primary radicals can be affected by the different reactivities between the molecules and the various excited dye molecules. This is important in the photoinitiation process, especially if one wishes to optimize the interaction between the electron donor and the photosensitizer. Therefore, and in agreement with previous results in the literature, in order to optimize a free-radical polymerization process [1, 16, 22, 25, 26, 32, 33], it is clear that a detailed model of dye kinetics, that is, photon absorption, photosensitizer recovery and bleaching is necessary. Once again we note that the behaviour of RB dye is different from that of the others, and once again this may be due to the presence of some other process.
6. Conclusion and Analysis
One of the aims of this paper has been the development of a new dye model which is both simple enough to use easily but also of practical value and sufficiently flexible to be used for different photosensitisers. In this paper, key dye parameter values for five photosensitizers, (i.e., molar absorptivity, , quantum efficiency of the reaction, and , recovery rate, , and bleaching rate, ), have been extracted by fitting the experiment data using our new simplified dye model. This is achieved by measuring the transmission of the material layer, for certain exposure intensities. The recovery behaviours of the dyes are thus experimentally compared under identical conditions. Following this, the bleaching processes [19] were examined using the parameter values already estimated. In this way a full theoretical treatment of the photoinitiation processes was presented. The results are very useful, that is, PB should be used to record gratings with low energetic exposure, in that the inhibition has a less effect on it.
While significant progress has been made in the work presented here, much remains to be done. It appears that some other reactions are taking placing during and after-exposure in the RB system case, indicating that a more accurate model is needed. The resulting accurate description of the photoinitiation processes must then be incorporated into the full nonlocal photopolymerization driven diffusion (NPDD) model, allowing a more complete modelling of free radical photopolymerization to take place. Such a model can then be applied to the study of photopolymer use for such application as holographic data storage and 3D hybrid optoelectronic circuit fabrication.
A new dye model, which includes the effect of photosensitizer recovery and bleaching has been developed and compared to experimental results. The photoabsorptive behaviours of five types of photosensitizer: Erythrosine B (EB); Eosin Y (EY); Phloxine B (PB); and Rose Bengal (RB), of the same molarity in equivalent PVA/AA layers, and Methylene Blue (MB) of higher molarity in PVA/AA layers, are experimentally examined. We modify the previous dye model by introducing and into the model, which are the rate constant of photon absorption from ground state dye to singlet state and triplet state respectively. A description of how to incorporate this dye model into a full NPDD material model, in order to simulate photopolymer material behaviour, has been presented and discussed.
Future work must include characterising the spatial frequency response of the AA/PVA material. There is a large difference between the molecular weights between five dyes. It is reasonable to assume that, barring an unusual molecular structure, the smaller the molecule, the quicker it will diffuse through a given polymer matrix [34]. It is expected that the faster the rate of diffusion, the more total polymerisation that can take place in the exposed areas. This should lead to a larger saturation value of refractive index modulation and a stronger material response. In summary the effect of the five different polarities, shapes and size of the Xanthene and Thiazine dyes on the spatial frequency response of the AA/PVA material must be quantitatively examined [35, 36]. Since the existence and diffusion of oxygen is non-ignorable in low intensity exposure cases, the inhibition process will be studied by using low exposure intensity in the future.
Acknowledgments
The authors would like to acknowledge the funding support of the EU ERUSMAS Mundus fund. They would also like to acknowledge the support of Enterprise Ireland and Science Foundation Ireland under the national development fund and the Irish Research Council for Science, Engineering and Technology (IRCSET) through the EMPOWER postgraduate and postdoctoral fellowship schemes and the Ministerio de EconomΓa y Competitividad of Spain under Project FIS2011-29803-C02-01 as well.