In order to use photopolymers in the recording of holographic memories, high physical thickness is required. This generates many problems associated with the attenuation of light in the recording due to Beer’s law. One of the more significant disadvantages is the fact that there are differences between the physical thickness of the material and the optical thickness of the holograms recorded. The optical thickness characterizes the angular selectivity of the holograms and determines the separation between two consecutive holograms in angular peristrophic multiplexing. In this work we propose a new method to record many holograms multiplexed with similar diffraction efficiency values taking into account the different effective optical thickness of each hologram.
© 2007 Optical Society of America
Recently, great advances have been made in the use of photopolymers as data storage media [1-3]. Photopolymer materials exhibit good optical properties (low losses, low scatter, high refractive index modulations and high density information storage) and are relatively inexpensive. Due to these reasons, some companies have presented the first holographic disks based on photopolymers [4-5]. One of the requirements to obtain competitive holographic memories is that the material must have a thickness of 500 µm or more . This gives rise to differences between the physical thickness of the material and the effective optical thickness of the holograms recorded. The effective optical thickness determines the Bragg selectivity and the separation between two consecutive holograms in some multiplexing techniques . Some of these aspects have been studied in hologram formation in these thick layers [7-8]. On the other hand, schedules have been proposed based on experimental and theoretical studies. For example J.T. Sheridan and co-workers proposed an interesting schedule based on a two dimensional non-local diffusion model to obtain holograms with the same diffraction efficiency using these materials . In these studies the authors assume that the diffraction efficiency only depends on the refractive index modulation of the grating because the thickness of the gratings is considered constant; in other words the physical and effective optical thickness have the same value. This assumption is valid for thin layers, but for thick photopolymer layers with high light absorption we need new methods to study the multiplexing schedule.
In our recent experimental studies  we observed some problems when multiplexing many holograms using iterative algorithms . These problems are due to the variation in the effective optical thickness of the gratings recorded during the multiplexing process. In Fig. 1 the experimental effective optical thickness of 60 gratings as a function of the hologram number is shown for peristrophic multiplexing. We fitted the angular response of each grating to obtain the effective optical thickness. The classical algorithm useful for non-overlapped gratings is the one proposed by Kogelnik . The first 10 gratings have very low effective optical thickness and it is difficult to read each grating separately and obtain the hologram thickness. The average hologram diffraction efficiency is 2%; afterwards, the transmission of the layer at the end of the multiplexing process is around 50%. In Fig. 1 we can see two important effects. Firstly we can see the low effective optical thickness of the first 10 holograms stored. Based on this result, more angular separation between these holograms is required and the capacity of the material decreases. The second interesting effect is observed in the analysis of the last 10 holograms. We can see a maximum effective optical thickness of around 470 µm, and then the thickness of the following holograms is lower. This can be explained if all the monomer is consumed on the illuminated polymer side, and the following gratings are recorded deeper in the material.
Therefore we present in this paper a model to optimize the scheduling, taking into account the different optical thickness of the gratings stored in the thick layers. We simulated the exposure time to record more than 200 gratings for different chemical compositions and characteristics of the dried layer. Nevertheless, for our simulations we use values of parameters characteristic of the thick layers based on PVA/acrylamide developed by our research group (such as the absorption, inhibition process, etc.). We studied the differences between the effective optical thicknesses of each grating and the influence of different parameters on the optimal times to obtain gratings with the same diffraction efficiencies.
2. Theoretical model
In previous studies we analyzed the importance of making some assumptions in order to present a simplified model which describes the photopolymer behavior. The first important approximation is that we only consider two harmonics in the concentration of monomer and polymer. We have demonstrated in previous papers that in gratings recorded in the thick layers of our material the higher harmonics are more than 20 times lower than the two main harmonics, which occurs because the monomer diffusion is 10-10 cm2 s-1 or higher . Another important assumption is that we are working with the local case, we disregard the non-local polymerization introduced by Sheridan and co-workers in the polymerization process . Then the equation that governs the gratings’ formation can be written as:
Where [M] is the monomer concentration, [P] is the polymer concentration, D is the monomer diffusivity and γ indicates the relationship between intensity and polymerization rate (kR). Assuming Beer’s law, which explains the attenuation of the light intensity with depth, the incident light inside the material can be written as:
Where α is the coefficient of the light attenuation. The initial value of α [α (t=0)=α0] can be obtained if the transmittance and the physical thickness of the layer before the recording process are known and depends on the dye concentration, the drying process and the recording wavelength. I0 takes the value of 5 mW cm-2 in this work.
We assume that the polymerization rate increases very quickly due to the Trommsdorff effect  and then decreases because it is limited by the viscosity of the material when the polymer is generated. It is given by:
Where φκ is the attenuation coefficient of the polymerization (we take φκ=5 in this work). Therefore the monomer diffusion inside the material decreases when the polymer is generated as follows:
Where φD is the attenuation coefficient of the monomer diffusivity (we take φD=5 in this work, slow diffusivity attenuation).
One important novelty introduced in the study of thick layers is the cut-off intensity presented by some materials for very low intensities. In PVA/AA materials used in our research laboratory the lower cut-off intensity is around 50 µW/cm2 . The importance of the inhibition period for the first hologram, due to the existence of oxygen and impurities inside the layer, should be borne in mind and we took into account this effect for each sub-layer. In our material we measured the energy density necessary to start the grating formation, which is around 6 mJ/cm2.
Another important approximation made in our model is that the monomer diffusion along the Z axis can be disregarded. The validity of this approximation is shown in reference  where a 3-dimensional model is solved by the finite differences method. Our simple approach to solve these equations in order to obtain qualitative information about the behaviour of the material in depth is as follows: The photopolymer film will be divided into G different sub-films each of thickness dg (Fig. 2); the total thickness of the photopolymer, d, is the sum of the thickness of the different sub-films.
where G is the number of sub-films used.
The exposure intensity in each sub-grating can be calculated using eq. (3) and the initial monomer concentration of each is M0, the same as that of the global grating. Knowing the initial conditions for each layer, the 2-dimensional model can be used to examine the 3-D problem .
In order to consider a particular photopolymer film useful for holographic memory applications, its thickness must be greater than 500 µm. In our research group the thickness of typical photopolymer layers is between 500 µm and 1000 µm. For this reason the minimum value of G is chosen to be 100. In this work the value taken for the thickness is 800 µm and G=800.
In this paper the model used to solve the diffusion equations in each grating will be a first harmonic diffusion model, because the presence of the higher harmonics can often be disregarded in photopolymers [12, 16].
In reference  the time between two consecutive gratings necessary to obtain a uniform distribution of monomer in the material again is studied. This time is around 10 seconds; in the schedule presented in this work, we assume that the monomer concentration is homogeneous before recording a new grating.
In order to record many holograms with the same diffraction efficiency (DE), it is important to analyze the dependence of the diffraction efficiency on the monomer and polymer concentrations. Classically this dependence can be written as follows :
We assume that the explicit equation of Kogelnik’s theory can be used to obtain the relation between the DE and n1 (average refractive index modulation in the layer). Nevertheless, there are two important facts to take into account in the parameters deff and n1. Firstly, it is important to understand that the thickness used is the effective optical thickness of the gratings (deff) . Secondly, the refractive index modulation is attenuated in depth, so an average value should be used . On the other hand, the refractive index modulation depends on the distribution of the monomer and polymer, but if at the end of the schedule we consume all the monomer or the residual monomer is homogenously distributed. We can write:
where ndark is the refractive index of the layer before exposure (ndark=1.478), nb is the binder refractive index (nb=1.474) np is the polymer refractive index and depends on the final length of the polymer chains (np=1.512 when BMA is used as crosslinker), and [P]1 is the first harmonic of the polymer concentration.
These values were obtained using the Lorentz-Lorenz equation. The calculations are based on refractometer measurements using water solutions.
Therefore if we want to multiplex many holograms with the same diffraction efficiency we need the product n1 deff to be constant. This condition has been introduced in our model when DE is fixed. On the other hand, the effective optical thickness deff determines the width of the central lobe in the gratings recorded and it is the parameter that governs the separation needed between two consecutive holograms so that they may be read separately. This was a problem in the first experimental results obtained by us in previous studies  where the first holograms have a very small effective thickness and are overlapped. An accurate prediction of the effective optical thickness and performance over time necessary to obtain equal diffraction efficiencies in each grating are the main goals of this work.
3. Results and discussions
The importance of the separation between two consecutive holographic gratings or data pages for them to be reproduced separately is fundamental in order to maximize the capacity of the information stored in photopolymers. In the first approximation of the problem of multiplexing many holograms in the same volume of light absorbent material, it is easy to understand that the effective optical thickness of the holograms will be greater initially and the limit of the effective optical thickness is the physical thickness of the layer. This first part of the work consists of the simulation of this type of materials and discussion of the importance of the bleaching process and absorption coefficient variation. Secondly, we present the analysis of some materials where the effective optical thickness never achieves values close to the physical thickness of the layer. This occurs in materials where the polymerization velocity is faster than the dye consumption rate.
3.1 High velocities of the bleaching process
The initial value of the light attenuation coefficient in depth of each photopolymeric material can be directly determined measuring the initial transmittance of the layer for the recording wavelength. In this section we use the measured valued α0=0.008 µm-1 for 800 µm thick layers recorded with a wavelength of 514 nm, with yellowish eosin as the dye . The variation in the attenuation of the light inside the material has a critical importance in the evolution of the effective optical thickness of the multiplexed holograms. It may be studied using only one incident beam and fitting the variation in the layer transmittance as a function of the time (for each incident intensity). This evolution has different velocities for different drying process times, drying temperatures, chemical compositions, etc. Therefore, very similar materials can present different behaviors.
The absorption of the layer changes as a function of the time when the dye is consumed; in a first approximation it can be written as:
where β is a constant that determines the influence on the intensity as the dye is consumed and depends on many factors (chemical composition of the material, temperature, humidity, etc.). The absorption decay depends on the dye and the intensity used, and Kα defines the rate at which the dye is consumed. The normal evolution of the transmittance of the layer for a power intensity of 5 mW/cm2 and β=1 is plotted in Fig. 3, for two different values of Kα, 0.001 cm2mW-1s-1 and 0.0005 cm2mW-1s-1
Bleaching velocity simulations of the 250 gratings recorded with a diffraction efficiency of 0.4% for these two cases were carried out. In this simulation, the total sum of the diffraction efficiency at each Bragg’s angle is around 100%. This type of photopolymer is similar to the material analyzed in references [19-20], but with lower light absorption. For both light absorption variations we can observe that there is a similar behavior as regards the time required to record the first 100 holograms (Fig. 4). After that, we can see an exponential increase in the time needed to achieve the diffraction efficiency of 0.4%, because the materials are near saturation. This effect can be seen when the photopolymer is near saturation, as can be seen in the experimental papers [9-10]. The common parameters used in the two simulations are kR=0.019 cm mW-1/2s-1, D0=10-10s cm-2 and γ=0.5.
The last results are simply the reproduction of the well-known experimental results [9-10]; the important part of this work is the prediction of the effective optical thickness to minimize the separation between the holograms and to maximize the material storage capacity. For the two materials simulated in this section, the effective optical thickness is the same as the physical thickness of the layer after 67 holograms are recorded in the first case and 124 in the second case (Fig. 5). This effect can be observed because the dye consumption is faster than the hologram formation, and can be demonstrated fitting the angular responses of each hologram.
3.2 High velocities of grating formation
Let us analyze in this section the importance of the velocity of the grating formation in the time scheme and in the effective optical thickness of the multiplexed gratings. The polymerization rate and the monomer diffusion inside the material depend on many factors, and the same compositions of material can present important differences if the drying conditions of the layer are different. On the other hand, when the photopolymer composition presents an excessive dye concentration or elevated photopolymerization rates (this case occurs when crosslinker monomers are used in the material composition ), two interesting effects may be seen. The time required to record holograms with 0.4% diffraction efficiency is lower (see Fig. 6, where kR=0.026cm2mW-1s-1, D0=3×10-10scm-2 and γ=1) as could be expected. Nevertheless, the most interesting effect can be seen in Fig. 7, where the effective optical thickness of the gratings recorded never achieves the value of the physical thickness of the grating. This effect can be seen in previous experimental papers  where the maximum effective optical thickness achieved after the recording of 60 gratings with an average diffraction efficiency of 2.4 % was 550 µm.
It is possible to make a deeper analysis of the last effect observed. In order to better appreciate the limit of the effective optical thickness, we simulated 600 grating with 0.4% diffraction efficiency (Fig. 8). But now the value of the initial light absorption considered is higher α0=0.012 µm-1.. This value can be obtained when a diode-pumped frequency-doubled Nd:YAG  laser at 532 nm is employed to store the multiplexed gratings in one layer as analyzed in Fig. 1. As we can see in Fig. 8, at the end of the multiplexing process the layer is near saturation and the time needed to multiplex one grating with 0.4% diffraction efficiency grows exponentially.
An important fact can be seen in Fig. 9, the evolution of effective optical thickness as a function of the multiplexed holograms presents a maximum around 600 µm after 350 holograms are recorded and then the effective optical thickness decays. This effect has been observed before in PVAAA based polymer and is due to the consumption of 90% of the monomer in the first 190 µm of the layer.
In order to better understand Fig. 9, let us analyze the monomer distribution simulation after 300 holograms and 600 holograms are recorded. In Fig. 10 we can see the gradient in depth of the monomer and diffusion concentration. After 300 holograms all the monomer in the first 160 µm has been consumed. At this moment, the effective optical thickness of the layers recorded is around 540 µm. Lastly, it is important to note that after 600 gratings are recorded only in the last 400 µm of the layer is there some free monomer.
Finally, let us explain better the effect analyzed and predicted in this work using the diagram presented in Fig. 11. Initially, for a small hologram number, the hologram’s effective optical thickness (the length of the zone were the hologram exists) is very low due to the short penetration depth of the light inside the material and only the monomer that is near the incident recording wave is consumed. When the hologram number increases, the monomer in this zone has been consumed and the following holograms are recorded deeper in the layer but the hologram’s effective optical thickness decreases.
The work has a great importance for the design of the multiplexing scheme in holographic memories based on photopolymers. The method presented can predict the times required to obtain many holograms multiplexed with equal diffraction efficiency and maximum storage capacity. Different types of materials have been studied in order to analyze the different dependence of the effective optical thickness on the number of holograms multiplexed. The analysis of the evolution of the effective optical thickness as a function of the exposure time can be used to characterize the hologram formation and to calculate some parameters such as the polymerization rate of the hologram.
This work was supported by the “Ministerio de Educación y ciencia”(Spain) under projects FIS2005-05881-C02-01 and FIS2005-05881-C02-02, and by Generalitat Valenciana under projects ACOMP/2007/020 and ACOMP/2007/027. We thank the referees for their interesting and very helpful comments, in particular regarding Fig. 11.
References and links
2. Y. Tomita, K. Furushima, K. Ochi, K. Ishizu, A. Tanaka, M. Ozawa, M. Hidaka, and K. Chikama, “Organic nanoparticle (hyperbranched polymer)-dispersed photopolymers for volume holographic storage,” Appl. Phys. Lett. 88, 071103 (2006). [CrossRef]
3. J. T. Sheridan, F. T. O’Neill, and J. V. Kelly, “Holographic data storage: optimized scheduling using the non local polymerization-driven diffusion model,” J. Opt. Soc. Am. B 20, 1443–1451 (2004). [CrossRef]
4. D. A. Waldman, C. J. Butler, and D. H Raguin, “CROP holographic storage media for optical data storage at grater than 100 bits/µm2,” Proc. SPIE 5216, 10–25 (2003). [CrossRef]
5. W. L. Wilson, K. R. Curtis, K. Anderson, M. C. Tackitt, A. J. Hill, M. Pane, C. Stanhope, T. Earhart, W. Loechel, C. Bergman, K. Wolfgang, C. Shuman, G. Hertrich, K. Parris, K. Malang, B. Riley, and M. Ayer, “Realization of high performance holographic data storage: The inPhase Technologies demonstration platform,” Proc. SPIE , 5216, 178–191 (2003). [CrossRef]
6. H. J. Coufal and D. PsaltisG. T. Sincerbox, eds., Holographic Data Storage, Springer Series in Optical Sciences, (Springer-Verlag, Berlin, 2000).
7. S. Gallego, M. Ortuño, C. Neipp, A. Márquez, A. Beléndez, I. Pascual, J. V. Kelly, and J. T. Sheridan, “Physical and effective optical thickness of holographic diffraction gratings recorded in photopolymers,” Opt. Express 13, 1939–1950 (2005). [CrossRef]
8. S. Gallego, M. Ortuño, C. Neipp, A. Márquez, A. Beléndez, I. Pascual, J. V. Kelly, and J. T. Sheridan, “3 dimensional analysis of holographic photopolymers based memories,” Opt. Express 13, 3543–3554 (2005). [CrossRef]
9. E. Fernández, C. García, M. Ortuño, S. Gallego, A. Belendez, and I. Pascual, “Optimization of a thick PVA/acrylamide photopolymer for data storage using a combination of angular and peristrophic holographic multiplexing,” Appl. Opt. 45, 7661–7666 (2006). [CrossRef]
11. H. Kogelnik, “Coupled wave theory for thick hologram gratings,” Bell Syst. Tech. J. 48, 2909–2947 (1969).
12. S. Gallego, C. Neipp, M. Ortuño, A. Benléndez, E. Fernández, and I. Pascual, “Analysis of monomer diffusion in depth in photopolymer materials,” Opt. Commun. 274, 43 (2007). [CrossRef]
13. J. T. Sheridan and J. R. Lawrence, “Nonlocal-response diffusion model of holographic recording in photopolymer,” J. Opt. Soc. Am. A 17, 1008–1014 (2000). [CrossRef]
14. J. V. Kelly, F. T. O’ Neill, C. Neipp, S. Gallego, M. Ortuno, and J. T. Sheridan, “Holographic photopolymer materials: non-local polymerisation driven diffusion under non-ideal kinetic conditions” J. Opt. Soc. of Am. B 22, 407–406 (2005). [CrossRef]
15. S. Gallego, M. Ortuño, C. Neipp, A. Márquez, A. Beléndez, E. Fernández, and I. Pascual “3-dimensional characterization of thick grating formation in PVA/AA based photopolymer,” Opt. Express 14, 5121–5128, (2006). [CrossRef]
16. C. Neipp, A. Beléndez, S. Gallego, M. Ortuño, I. Pascual, and J. T. Sheridan, “Angular responses of the first and second diffracted orders in transmission diffraction grating recorded on photopolymer material,” Opt. Express 11, 1835–1843 (2003). [CrossRef]
17. C. Neipp, J. T. Sheridan, S. Gallego, M. Ortuño, I. Pascual, and A. Beléndez “Effect of a depth attenuated refractive index profile in the angular responses of the efficiency of higher orders in volume gratings recorded in a PVA/Acrylamide photopolymer,” Opt. Comun. 233, 311–322 (2004). [CrossRef]
18. M. Ortuño, S. Gallego, C. García, C. Neipp, A. Beléndez, and I. Pascual, “Optimization of a 1 mm thick PVA/acrylamide recording material to obtain holographic memories: method of preparation and holographic properties,” Appl. Phys. B 76, 851–857 (2003). [CrossRef]
19. M. Ortuño, S. Gallego, C. García, C. Neipp, and I. Pascual, “Holographic characteristics of a 1 mm thick photopolymer to be used in holographic memories,” Appl. Opt 42, 7008–7012 (2003). [CrossRef]
20. C. Neipp, S. Gallego, M. Ortuño, A. Márquez, A. Beléndez, and I. Pascual, Characterization of a PVA/acrylamide photopolymer. Influence of a cross-linking monomer in the final characteristics of the hologram” Optics Commun. 224, 27–34 (2003). [CrossRef]